爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 17381|回复: 3

[讨论] 使用m_streamline出现问题

[复制链接]

新浪微博达人勋

发表于 2020-7-20 10:19:32 | 显示全部楼层 |阅读模式
20金钱
这里的um,vm都是lon*lat的格式,后边的那个10是看二爷给的示例里写的。

                               
登录/注册后可看大图
然后就开始报错。函数里的

                               
登录/注册后可看大图
这一句计算出来的X,Y有nan。


                               
登录/注册后可看大图


有没有成功使用过这个函数的人教一下,该如何正确调用。感谢

test.mat

314.81 KB, 下载次数: 0, 下载积分: 金钱 -5

最佳答案

查看完整内容

Update 2020.05.06 bug:经查明,本函数在面临以下情况时,会出现问题: 1.涉及到180~360°的经度范围; 2.涉及到大于85°的纬度; 当数据面临以上情况时,matlab会报错,本人后期将进行相应修正,目前建议:1.经度范围设为-180~180°范围内,如果是180-360范围,请减去180之后进行相应设置; 2.纬度范围建议设置在-85~85°范围内,由于圆柱投影在南北两级变形较大,无法正确绘制出该区域图像。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2020-7-20 10:19:33 | 显示全部楼层
Update 2020.05.06
bug:经查明,本函数在面临以下情况时,会出现问题:
1.涉及到180~360°的经度范围;
2.涉及到大于85°的纬度;
当数据面临以上情况时,matlab会报错,本人后期将进行相应修正,目前建议:1.经度范围设为-180~180°范围内,如果是180-360范围,请减去180之后进行相应设置;
2.纬度范围建议设置在-85~85°范围内,由于圆柱投影在南北两级变形较大,无法正确绘制出该区域图像。
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-7-20 10:24:26 | 显示全部楼层
这是放不上截图的吗。  
这是调用时候的代码:
[LAT,LON]=meshgrid(lats_u,lons_u);
figure;
m_proj('miller','lons',[20 120],'lats',[-60 -20]);
m_streamline(LON,LAT,um,vm,10); hold on
m_coast('patch',[.7 .7 .7]);
m_grid('fontsize',15,'ytick',[-60:10:-20]);

报错信息:
错误使用 griddedInterpolant
输入点的坐标必须为有限值;不允许使用 Inf 和 NaN。

出错 interp1 (line 161)
        F = griddedInterpolant(X,V,method);

出错 stream2 (line 62)
    sxi=interp1(xx(:),1:szu(2),sx(k));

出错 streamslice>nicestreams (line 309)
            vertsf = stream2(x,y, u, v,xstart,ystart,streamoptions);

出错 streamslice (line 134)
        [vertsout, avertsout] = nicestreams(x,y,u,v,density,arrows);

出错 m_streamline (line 79)
h=streamslice(X,Y,mU,mV,varargin{:});


返回值存在nan使得插值有问题的一句
[X,Y]=m_ll2xy(long,lat,'clip','point');
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-8-17 09:39:24 | 显示全部楼层
使用时不能包含nan值,只考虑海洋内的就好了
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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