首先下载得到全球叶绿素含量的数据,格式是csv格式,不含经纬度信息。注意到文件共有1800*3600,可知全球的分辨率是0.1°×0.1°。
在matlab中打开
S = csvread('MY1DMM_CHLORA_2009-05-01_rgb_3600x1800.CSV');
% imagesc(S)
lon = 0.1:0.1:360;
lat = -89.9:0.1:90;
[lon,lat] = meshgrid(lon,lat);
O.lon = lon;
O.lat = lat;
O.rg = flipud(S);
O.rg(O.rg==99999)=0;
O.rg(O.rg>5)=0;
O.rg(O.rg==0)=NaN;
O.rg(O.rg>5)=NaN;
% rg_plot(O)
m_proj('miller','lon',[280,315],'lat',[0,45]); 需要先配置 网页链接
% m_coast('patch',[.7 .7 .7],'edgecolor','none');
m_pcolor(lon,lat,O.rg);
m_grid('tickdir','out','linewi',2);
结果如下:
欢迎交流学习!1565903211