matlab绘制春季总降水量图

NIU凯

%%数据分析
clear;clc;
ncdisp('D:\data\cru_ts4.05.1901.2020.pre.dat.nc','/','full');
ncinfo('D:\data\cru_ts4.05.1901.2020.pre.dat.nc');
time=ncread('D:\data\cru_ts4.05.1901.2020.pre.dat.nc','time');
lat=double(ncread('D:\data\cru_ts4.05.1901.2020.pre.dat.nc','lat'));
lon=double(ncread('D:\data\cru_ts4.05.1901.2020.pre.dat.nc','lon'));
pre0=ncread('D:\data\cru_ts4.05.1901.2020.pre.dat.nc','pre');
pre0count=length(pre0);
timecount=length(time);
pre0(isnan(pre0))=0;%讲NaN转化为0

ii=2020-1901+1;
% pre0_3=zeros(m1,n1,ii);pre0_4=zeros(m1,n1,ii);pre0_5=zeros(m1,n1,ii);
for k=1:1:ii
    sum=0;
    pre0_3=sum+pre0(:,:,(k-1)*12+3);
    pre0_4=sum+pre0(:,:,(k-1)*12+4);
    pre0_5=sum+pre0(:,:,(k-1)*12+4);
end
sum_pre0=pre0_3+pre0_4+pre0_5;
sum_pre0=imrotate(sum_pre0,90);

%%开始绘图
m_proj('Mercator','lon',[100 130],'lat',[10 50]);
m_contourf(lat,lon,sum_pre0,'linestle','none');
colormap(flipud(m_colmap('jet')));
m_coast('patch',[.7,.7,.7],'linewidth',2);
m_grid('box','fancy');



这个代码有什么问题呀？
出来的图没有颜色

伽蓝鸟

你是不是没有hold on，于是只有最后的海岸线？不过我不懂m_map工具箱。。

474620961

"m_contourf(lat,lon,sum_pre0,'linestle','none')" 改成"m_contourf(lon,lat,sum_pre0','linestle','none')",不然在你投影的区域没有数值

NIU凯

474620961 发表于 2021-11-26 09:38
"m_contourf(lat,lon,sum_pre0,'linestle','none')" 改成"m_contourf(lon,lat,sum_pre0','linestle','none' ...

好的好的，我试试，多谢