爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 2593|回复: 4

[程序设计] 用griddata将沿轨风场数据网格化

[复制链接]
发表于 2023-4-11 17:07:56 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

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

x
用linear方法进行时出现这种图,用natural方法就是正常的图,这是为什么,我记得我当时用linear方法可以画出正常的图,但是不知道改动了那里就出现这种情况。
QQ图片20230411170725.png
密码修改失败请联系微信:mofangbao
发表于 2023-4-11 18:49:10 | 显示全部楼层
信息量太少了 无法判断出错原因
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-4-11 19:13:18 | 显示全部楼层
本帖最后由 echo111 于 2023-4-11 19:22 编辑
wjy_ecnu 发表于 2023-4-11 18:49
信息量太少了 无法判断出错原因

speed_2b = h5read('D:\风场数据\HY-2B\H2B_OPER_SCA_L2B_OR_20220401T022839_20220401T041316_17223_dps_250_21_owv.h5','/wind_speed_selection');
dir_2b = h5read('D:\风场数据\HY-2B\H2B_OPER_SCA_L2B_OR_20220401T022839_20220401T041316_17223_dps_250_21_owv.h5','/wind_dir_selection');
speed_2b(speed_2b<0)=nan;
dir_2b(dir_2b<0)=nan;
u_2b = speed_2b.*sin((dir_2b-180)*pi/180);
v_2b = speed_2b.*cos((dir_2b-180)*pi/180);
lat_2b=double(lat_2b);
lon_2b=double(lon_2b);
[xi,yi]=meshgrid(-90:0.25:90,0:0.25:360);
ucha1=griddata(lat_2b,lon_2b,u_2b,xi,yi,'linear');
数据是h5类型的沿轨风场数据,经纬度风场数据都是二维的大矩阵,矩阵某一位置对应的的经纬度和风场数据是一个散点,把他们进行网格化
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2023-4-12 09:28:10 | 显示全部楼层
双线性插值 如果NaN过多 插值结果也是NaN  你这里dir_2b小于0的设成NaN  可能是原因
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-4-14 20:14:18 | 显示全部楼层
我这个图中间不应该有值的地方被插入了值,请问这个是为什么呢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表