数据来自Etopo1全球地形和水深数据。
其分为两个版本,Ice Surface和Bedrock,两个版本基本一致。不同之处在于在处理南极洲和Greenland地形时,前者给出的是加上冰盖层之后的高程,后者给出的是岩床的高程。
在每个子版本下又有多种数据格式, netCDF,binary, xyz,tiff。
结合M_Map画南海水深地形图。
clear;
clf;
latlim = [-4,26]; lonlim = [99,126]; [Z, refvec] = etopo('etopo1_bed_c_f4.flt', 1, latlim, lonlim); Z(Z>0)=NaN; lon=linspace(lonlim(1),lonlim(2),size(Z,2)); lat=linspace(latlim(1),latlim(2),size(Z,1)); [lon,lat]=meshgrid(lon,lat); m_proj('miller','lon',lonlim,'lat',latlim) m_gshhs_i('patch',[0.7 0.7 0.7]); %c'box','fancy'); title('南海地形图','fontsize',15,'fontweight','bold') hold on;
[c,h]=m_contour(lon,lat,Z,[0,-100,-200,-500,-1000,-2000,-5000,-10000],'linewidth',0.5); colormap(summer) colorbar;
等值线图为:
用颜色填充:
[c,h]=m_contourf(lon,lat,Z,[0,-100,-200,-500,-1000,-2000,-5000,-10000]);
colormap(summer)