M_map画流线，m_streanline函数调试

落日飞车

在写课程小组作业时需要用到m_streamline画流线，但是频繁报错：
（1）[X,Y]=m_ll2xy(long,lat,'clip','point'); 之后出现大量nan值，需要把经纬度限制在 lon:-180-180 lat:-85-85就不会出现nan了；
（2）h=streamslice(X,Y,mU,mV,varargin{:}); 这里varargin用中括号会报错，改成varargin（:）就不会报错了，但我调用的时候会报错，说什么sqrt（dendity）未定义，我干脆就去掉了这一参数，这样他默认值取1；
（3）h=streamslice(X,Y,mU,mV）XY应该时一维递增的坐标序列，需要取第一列和第一行。报错说维度不匹配，把UV转置一下就好了。
我的调用内容：
m_proj('miller','lon',[-120,60],'lat',[0,80]);
[cs,h]=m_contourf(lon,lat,H2(:,:,3),50,'LineStyle','none');
hold on;
h1 =  m_streamline(lon,lat,U(:,:,１),V(:,:,１));
for i = 1:length(h1)
　　h1(i).Color = 'k';
end

出图和修改后的代码放在这里供有需要的童鞋使用～

m_streamline.m (3.66 KB, 下载次数: 22)

2022-5-1 11:54 上传

点击文件名下载附件

落日飞车

啊噢 m_grid没有改 ，请大家忽略细节~

落日飞车

诶 家人们 我发现要用这个函数，经度是0-360的话好麻烦。遂弃

落日飞车

后来试了试 可以用的，按照上面说的画出来图之后，再添加一个地图底图覆盖错误底图即可。

碧落

感谢！帮到了2年后的我！用的是ERA5的再分析资料，家园里找了二爷两个版本，注意了经纬度范围避免nan的错误后，“griddedInterpolant插值要求每个维度至少有两个样本点。”这个报错怎么都搞不定了。贴主是直接改了二爷的m_streanline函数脚本中streamslice函数的参数，直接用贴主的函数替代掉toolbox/m_map下的原版本，确实就画出图了！