MATLAB--绘制高分辨率地形图

m_map的基本参数选取

m_proj 设置投影

m_grid设置网格线和边框

m_ghhs 添加海岸线、国界线和河流，可通过help m_ghhs查看

一般的画图命令在m_map里都有所涵盖，不过要改成m_plot等，主要函数包括：

m_plot(LONG,LAT,...line properties)

m_line(LONG,LAT,...line properties)      % 线条

m_text(LONG,LAT,‘string’)                % 文本 

m_quiver(LONG,LAT,U,V,S)                % 矢量箭头

m_patch(LONG,LAT,..patch properties)    % 色块

m_annotation(‘line’,LON,LAT)            % 注释

m_contour(LONG,LAT,VALUES)         % 等值线

m_contourf(LONG,LAT,VALUES)         % 等值线填充

m_image(LON,LAT,DATA) % 影像

m_pcolor(LON,LAT,DATA) % 色块

[IM,X,Y]=m_image(LON,LAT,DATA); % 将经纬度坐标转为XY坐标

m_shadedrelief(X,Y,IM,'coords','map’) % 地形渲染

m_etopo2(OPTION)              % 画地形图

m_ruler([.5 .8],.9,'tickdir','out','ticklen',[.007 .007]);    % 添加距离比例尺

m_northarrow(-123-4.5/60,49+19.5/60,1/60,'type',4,'aspect',1.5);  % 添加指北针

Code：

clear;clc;

station_city_lat = [30.2861,30.17];

station_city_lon = [109.471,109.16];

station_name = {'ESRS','ES'};

%中国省界线文件（含九段线）

China_bound    = shaperead('S:\5.模式与编程学习\重要的画图数据\国家基础地理数据\bou1_4m\bou1_4l.shp');

China_province = shaperead('S:\5.模式与编程学习\重要的画图数据\国家基础地理数据\bou2_4m\bou2_4l.shp');

China_diqu    = shaperead('S:\5.模式与编程学习\重要的画图数据\国家基础地理数据\bou3_4m\diquJie_polyline.shp');

bou1_4lx=[China_bound(:).X];bou1_4ly=[China_bound(:).Y];

bou2_4lx=[China_province(:).X];bou2_4ly=[China_province(:).Y];

bou3_4lx=[China_diqu(:).X];bou3_4ly=[China_diqu(:).Y];

fontsize = 16;

figure('position',[100 100 800 600])

m_proj('miller','lon',[72,136],'lat',[18,54]);hold on;  %设置投影方式和区域

load('S:\5.模式与编程学习\MATLAB\colorbar\OceanLakeLandSnow.mat');

colormap(OceanLakeLandSnow);caxis([0 4000]);

m_etopo2('shadedrelief','gradient',200);hold on;  %绘制地形图，200m间隔

ch = colorbar('southoutside');

set(get(ch,'title'),'string','ALT (m)','position',[220 -40],'FontSize',14);

axis off;

m_gshhs('fr1','color','b')    %添加河流

% m_gshhs('fb2','color','k'); %添加国界线

% m_gshhs('hc','color','k');  %添加海岸线

m_plot(bou1_4lx,bou1_4ly,'k','linewidth',1);hold on; %绘制中国国界

m_plot(bou2_4lx,bou2_4ly,'k','linewidth',1);hold on; %绘制中国省界

% m_plot(bou3_4lx,bou3_4ly,'k','linewidth',1);hold on; %绘制中国地市界

% 添加站点位置

m_plot(station_city_lon(1),station_city_lat(1),'og','color',[0.1 0.1 0.4],'LineWidth',4,'MarkerSize',2);hold on;

m_plot(station_city_lon(2),station_city_lat(2),'or','color',[0.1 0.1 0.4],'LineWidth',4,'MarkerSize',2);hold on;

% 添加研究区域位置

m_line([106 112],[28 28],'color',[0.1 0.1 0.9],'linewidth',2);hold on;

m_line([106 112],[32 32],'color',[0.1 0.1 0.9],'linewidth',2);hold on;

m_line([106 106],[28 32],'color',[0.1 0.1 0.9],'linewidth',2);hold on;

m_line([112 112],[28 32],'color',[0.1 0.1 0.9],'linewidth',2);hold on;

set(gca,'FontSize',fontsize,'FontName','Times New Roman');

m_grid('linestyle','none'); %设置投影的网格与轴