请教matlab画NCEP FNL再分析资料等值线平滑问题

jiajian0512 · 发表于 2014-10-25 13:25:25

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 jiajian0512 于 2014-10-25 14:44 编辑

<p>% 850 hPa温度平流算法
clear all;
clc;
grib_struct=read_grib('g:\031200',-1); % 读取NCEP再分析资料
read_grib('g:\031200','inv'); %生成变量目录表文档
Hight_500hPa=grib_struct(154).fltarray;</p><p>U_850hPa=grib_struct(157).fltarray;
V_850hPa=grib_struct(158).fltarray;
for i=1:181;
for j=1:360;
lon(i,j)=181-j;
lat(i,j)=91-i;
T(i,j)=Hight_500hPa((i-1)*360+j)/10.0;
u(i,j)=U_850hPa((i-1)*360+j);
v(i,j)=V_850hPa((i-1)*360+j);
end
end
R=6.3781e6;
[dx, ~]=gradient(lon);
dx=dx.*(pi./180).*R.*cos(lat*pi./180;
[~, dy]=gradient(lat);
dy=dy.*(pi./180).*R;
[dTx dTy]=gradient(T);
advh=-(u.*dTx./dx+v.*dTy./dy);</p>

复制代码

画图用m_map 工具箱，格点数据181X360，出来的图等值线有明显尖点，请问如何能画出平滑的等值线？能具体的给个函数之类吗？前面试过再把网格按经纬度画细，用到meshgrid griddata，也不知道这样和grads画出的图是否一致？希望得到大家帮助，不胜感激！

jiajian0512 · 发表于 2014-10-26 00:12:45

matlab有什么好办法？

kongfeng0824 · 发表于 2014-10-25 21:19:38

建议算完之后导出成geotiff，用arcgis画图

jiajian0512 · 发表于 2014-10-25 23:57:45

arcgis好上手吗请问

故乡月(环) · 发表于 2015-5-7 16:03:51

楼主厉害！小弟的向你多请教！

jiajian0512 · 发表于 2015-5-7 16:56:43

故乡月(环) 发表于 2015-5-7 16:03
楼主厉害！小弟的向你多请教！

故乡月(环) · 发表于 2015-5-12 22:08:54

jiajian0512 发表于 2015-5-7 16:56

楼主后面怎么解决的，能否指点一下，我也在用fnl资料画图！

zqlbv · 发表于 2017-5-23 08:40:56

谢谢楼主分享

哈橙子1989 · 发表于 2017-7-30 12:44:07

谢谢分享

liuyg123 · 发表于 2017-7-31 16:11:44

谢谢分享

		自动登录	找回密码
密码			立即注册

[讨论] 请教matlab画NCEP FNL再分析资料等值线平滑问题


12 金钱	回复本帖可获得 1 金钱奖励! 每人限 1 次