M_Map畫南海水深地形圖

數據來自Etopo1全球地形和水深數據。

其分為兩個版本，Ice Surface和Bedrock，兩個版本基本一致。不同之處在於在處理南極洲和Greenland地形時，前者給出的是加上冰蓋層之后的高程，后者給出的是岩床的高程。

在每個子版本下又有多種數據格式， netCDF，binary， xyz，tiff。

結合M_Map畫南海水深地形圖。

clear;
close all;
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)