1.前言
? ? ? ?m_tbase是matlab的地圖工具包m_map的一個(gè)內(nèi)置腳本者填,其提供了ETOPO的全球5'X5'的分辨率高程數(shù)據(jù)澜术,根據(jù)命令[ELEV,LONG,LAT]=m_tbase([LONG_MIN LONG_MAX LAT_MIN LAT_MAX])可以獲得指定[LONG_MIN LONG_MAX LAT_MIN LAT_MAX]范圍內(nèi)的高程數(shù)據(jù)庭再。本文通過m_tbase獲得了全球范圍的高程數(shù)據(jù)膳汪,并將海洋部分掩膜掉茄袖,只可視化陸地的高程數(shù)據(jù)渺贤。通過本文的操作也可以實(shí)現(xiàn)陸地點(diǎn)與海洋點(diǎn)的簡單區(qū)分雏胃。
2.全球高程數(shù)據(jù)
%% 繪制全球高程水深
[ELEV,LONG,Lat]=m_tbase([-180 180 -90 90]);%提取全球高程數(shù)據(jù)m_proj('robinson') ;m_pcolor(LONG,Lat,ELEV);
m_grid('linest','none','yticklabels','','xticklabels','')set(gcf,'color','w')set(gcf,'GraphicsSmoothing','off');%關(guān)閉圖像平滑
3.繪制陸地地形高程圖
繪制陸地地形高程圖志鞍,主要是在第2步的基礎(chǔ)上瞭亮,將高程小于等于0的部分掩膜掉即可。
%% 陸地地形高程Land=ELEV;Land(Land<=0)=nan; %掩膜掉高程數(shù)據(jù)<=0的即可figurem_proj('robinson')m_pcolor(LONG,Lat,Land);m_grid('linest','none','yticklabels','','xticklabels','')colorbarcolormap('summer')set(gcf,'color','w')set(gcf,'GraphicsSmoothing','off');%關(guān)閉圖像平滑
4.繪制全球水深圖
繪制陸地地形高程圖固棚,主要是在第2步的基礎(chǔ)上统翩,將高程大于0的部分掩膜掉即可仙蚜。
%% 全球水深圖Sea=ELEV;Sea(Sea>0)=nan; %掩膜掉高程數(shù)據(jù)>0的即可figurem_proj('robinson')m_pcolor(LONG,Lat,Sea);m_grid('linest','none','yticklabels','','xticklabels','')colorbarcolormap('winter')set(gcf,'color','w')set(gcf,'GraphicsSmoothing','off');%關(guān)閉圖像平滑