m_map绘制地形图时，虽然自带有1°的地形图以及从NOAA下载的1分的地形图（详见：Matlab下地形图绘图包m_map安装与使用），但有时需要对地形图分辨率的要求更高，便无法满足。
此时，需要导入本地地形数据文件，例如多波束测深。
这里以从https://www.ncei.noaa.gov/maps/grid-extract/下载的15秒地形数据为例，展示如何导入本地地形数据并绘图。
关键函数为m_shadedrelief。
1. 1分的分辨率地图
m_proj('mercator','long',[118 124],'lat',[20 24]);
caxis([-6000 6000])
colormap(slanCM('terrain'))
hc=colorbar;
set(get(hc,'title'),'string','Elevation(m)')
set(hc,'tickdir','out')
m_etopo2('shadedrelief','lightangle',45,'gradient',5);
m_gshhs('ic','color','k')
m_grid('box','on','tickdir','out','gridlines','no')
set(gcf,'position',[10 10 1200 800])
figname='test1'
print('-dpng','-r1000',[figname,'.png']) % 导出png图片
结果如下：
2. 15秒的分辨率地图
load('scs_elevation.mat')
elev=flipud(elev);
lonlim=[extent(1:2)]
latlim=[extent(3:4)]
[nlat,nlon]=size(elev);
Lon=[linspace(lonlim(1),lonlim(2),nlon)];
Lat=[linspace(latlim(1),latlim(2),nlat)]';
m_proj('mercator','long',[118 124],'lat',[20 24]);
caxis([-6000 6000])
colormap(slanCM('terrain'))
hc=colorbar;
set(get(hc,'title'),'string','Elevation(m)')
set(hc,'tickdir','out')
m_shadedrelief(Lon,Lat,elev,'lightangle',45,'gradient',5)
m_gshhs('ic','color','k')
m_grid('box','on','tickdir','out','gridlines','no')
set(gcf,'position',[10 10 1200 800])
figname='test2'
print('-dpng','-r600',[figname,'.png']) % 导出png图片
- 可见15秒分辨率地图比1分地图能显示更多细节。
注：绘图色标参考自：https://blog.csdn.net/slandarer/article/details/127719784
