- 积分
- 247
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-11-1
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
[lon lat]=meshgrid([73:0.1:104],[25:0.1:40]);
% s=griddata(s(:,2),s(:,1),lon,lat,'v4');
tibet=shaperead('tibetan.shp');
R=makerefmat('RasterSize',size(s),'Lonlim',[73 104],'Latlim',[25 40]);
MASK=vec2mtx(tibet.Y,tibet.X,s,R,'filled');
MASK(find(MASK>1))=nan;
MASK(find(MASK==1))=0;
%contourf(s);
contourf(lon,lat,s+MASK,30);
shading flat;
colorbar;
hold on;
plot(tibet.X,tibet.Y,'-k','linewidth',3);
hold off;
画出来后总是出现X的列数与s的列数不同。怎么修改?急求
|
|