M_Map画南海水深地形图

Posted 2021-02-15 keegentang

数据来自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<l<i<h
m_grid(‘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)