m_map 实例
1、
clear all
m_proj(\'ortho\',\'lat\', 48,\'long\',-123\');%投影方式，范围
m_coast(\'patch\',\'r\');%红色填充
m_grid(\'linest\',\'-\',\'xticklabels\',[],\'yticklabels\',[]);%标注为空
patch(.55*[-1 1 1 -1],.25*[-1 -1 1 1]-.55,\'w\');%四个点，白色填充
text(0,-.55,\'M\_Map\',\'fontsize\',25,\'color\',\'b\',...
\'vertical\',\'middle\',\'horizontal\',\'center\');%填充内容，属性要求，垂直水平居中，斜杠转义符
set(gcf,\'units\',\'inches\',\'position\',[2 2 3 3]);%设置图片大小
set(gcf,\'paperposition\',[3 16 6 6]);%[left bottom width height].
2.
clear all
m_proj(\'lambert\',\'long\',[-160 -40],\'lat\',[30 80]);
m_coast(\'patch\',[1 .85 .7]);%填充海岸线
m_elev(\'contourf\',[500:500:6000]);%海深等高线
m_grid(\'box\',\'fancy\',\'tickdir\',\'in\');
colormap(flipud(copper));%翻转三原色的顺序排列，重新画图
3、
m_proj(\'stereographic\',\'lat\',90,\'long\',30,\'radius\',25);
m_elev(\'contour\',[-3500:1000:-500],\'edgecolor\',\'b\');
m_grid(\'xtick\',12,\'tickdir\',\'out\',\'ytick\',[70 80],\'linest\',\'-\');
m_coast(\'patch\',[.7 .7 .7],\'edgecolor\',\'r\');
4、
4.1默认的加载形式
clear all
Slongs=[-100 0;-75 25;-5 45; 25 145;45 100;145 295;100 290];
Slats= [ 8 80;-80 8; 8 80;-80 8; 8 80;-80 0; 0 80];
for l=1:7,
m_proj(\'sinusoidal\',\'long\',Slongs(l,:),\'lat\',Slats(l,:));%分区块的载入地图
m_grid;
m_coast;
end;
xlabel(\'Interrupted Sinusoidal Projection of World Oceans\');
% In order to see all the maps we must undo the axis limits set by m_grid calls:
set(gca,\'xlimmode\',\'auto\',\'ylimmode\',\'auto\');%自动加载所有数据
4.2
clear all
subplot(211);
Slongs=[-100 0;-75 25;-5 45; 25 145;45 100;145 295;100 290];
Slats= [ 8 80;-80 8; 8 80;-80 8; 8 80;-80 0; 0 80];
for l=1:7,
m_proj(\'sinusoidal\',\'long\',Slongs(l,:),\'lat\',Slats(l,:));
m_grid(\'fontsize\',6,\'xticklabels\',[],\'xtick\',[-180:30:360],...%字体尺寸，隐藏坐标，坐标等距离排列
\'ytick\',[-80:20:80],\'yticklabels\',[],\'linest\',\'-\',\'color\',[.9 .9 .9]);%线条形式实线，颜色
m_coast(\'patch\',\'g\');%填充绿色
end;
xlabel(\'Interrupted Sinusoidal Projection of World Oceans\');
% In order to see all the maps we must undo the axis limits set by m_grid calls:
set(gca,\'xlimmode\',\'auto\',\'ylimmode\',\'auto\');
subplot(212);
Slongs=[-100 43;-75 20; 20 145;43 100;145 295;100 295];
Slats= [ 0 90;-90 0;-90 0; 0 90;-90 0; 0 90];%比上面范围重叠更多
for l=1:6,
m_proj(\'mollweide\',\'long\',Slongs(l,:),\'lat\',Slats(l,:));
m_grid(\'fontsize\',6,\'xticklabels\',[],\'xtick\',[-180:30:360],...
\'ytick\',[-80:20:80],\'yticklabels\',[],\'linest\',\'-\',\'color\',\'k\');
m_coast(\'patch\',[.6 .6 .6]);
end;
xlabel(\'Interrupted Mollweide Projection of World Oceans\');
set(gca,\'xlimmode\',\'auto\',\'ylimmode\',\'auto\');%自动加载所有数据范围
5、
clear all
% Nice looking data
[lon,lat]=meshgrid([-136:2:-114],[36:2:54]);
u=sin(lat/6);%？
v=sin(lon/6);
m_proj(\'oblique\',\'lat\',[56 30],\'lon\',[-132 -120],\'aspect\',.8);%注意纬度顺序，高纬到低纬，\'aspect\'长宽比
subplot(121);
m_coast(\'patch\',[.9 .9 .9],\'edgecolor\',\'none\');
m_grid(\'tickdir\',\'out\',\'yaxislocation\',\'right\',...%边缘经纬线在外，纬度在右，经度在顶
\'xaxislocation\',\'top\',\'xlabeldir\',\'end\',\'ticklen\',.05);%边缘经纬度宽度
hold on;%图层合并
m_quiver(lon,lat,u,v);%经度，纬度，向东，向北的速度分量，转化为矢量箭头
xlabel(\'Simulated surface winds\');
subplot(122);
m_coast(\'patch\',[.9 .9 .9],\'edgecolor\',\'none\');
m_grid(\'tickdir\',\'in\',\'yticklabels\',[],...
\'xticklabels\',[],\'linestyle\',\'none\',\'ticklen\',.02);%经纬线去除
hold on;
[cs,h]=m_contour(lon,lat,sqrt(u.*u+v.*v));%画出轮廓
clabel(cs,h,\'fontsize\',10,\'Color\',\'red\',\'FontWeight\',\'bold\');%标注轮廓，字体，字体大小，颜色
xlabel(\'Simulated something else\');
6、
clear all
% Plot a circular orbit
lon=[-180:180];
lat=atan(tan(60*pi/180)*cos((lon-30)*pi/180))*180/pi;
m_proj(\'miller\',\'lat\',80);%查询，纬度（-80 80），经度（-180 180）
m_coast(\'color\',[0 .6 0]);
m_line(lon,lat,\'linewi\',2,\'color\',\'r\');%线宽，2；颜色
m_grid(\'linestyle\',\'none\',\'box\',\'fancy\',\'tickdir\',\'in\');%没有网格，边框相间，
%m_line(lon,lat,\'linewi\',2,\'color\',\'r\',\'linestyle\',\':\');控制线条格式，点画线还是直线
7、
clear all
m_proj(\'lambert\',\'lon\',[-10 20],\'lat\',[33 48]);%范围设置与投影有关系
m_tbase(\'contourf\');%画5-min数据库轮廓
m_grid(\'linestyle\',\'none\',\'tickdir\',\'out\',\'linewidth\',3);
没有安装terrainbase,所以轮廓很简陋，这个坑用到的时候再填
8、不同分辨率情况下画的图
clear all
% Example showing the default coastline and all of the different resolutions
% of GSHHS coastlines as we zoom in on a section of Prince Edward Island.
clf
axes(\'position\',[0.35 0.6 .37 .37]);%[left bottom width height]
m_proj(\'albers equal-area\',\'lat\',[40 60],\'long\',[-90 -50],\'rect\',\'on\');%方形，扇形
m_coast(\'patch\',[0 1 0]);%填充
m_grid(\'linest\',\'none\',\'linewidth\',2,\'tickdir\',\'out\',\'xaxisloc\',\'top\',\'yaxisloc\',\'right\');
m_text(-69,41,\'Standard coastline\',\'color\',\'r\',\'fontweight\',\'bold\');%经纬度，内容，格式
axes(\'position\',[.09 .5 .37 .37]);
m_proj(\'albers equal-area\',\'lat\',[40 54],\'long\',[-80 -55],\'rect\',\'on\');
m_gshhs_c(\'patch\',[.2 .8 .2]);
m_grid(\'linest\',\'none\',\'linewidth\',2,\'tickdir\',\'out\',\'xaxisloc\',\'top\');
m_text(-80,52.5,\'GSHHS\_C (crude)\',\'color\',\'m\',\'fontweight\',\'bold\',\'fontsize\',14);
axes(\'position\',[.13 .2 .37 .37]);
m_proj(\'albers equal-area\',\'lat\',[43 48],\'long\',[-67 -59],\'rect\',\'on\');
m_gshhs_l(\'patch\',[.4 .6 .4]);
m_grid(\'linest\',\'none\',\'linewidth\',2,\'tickdir\',\'out\');
m_text(-66.5,43.5,\'GSHHS\_L (low)\',\'color\',\'m\',\'fontweight\',\'bold\',\'fontsize\',14);
axes(\'position\',[.35 .05 .37 .37]);
m_proj(\'albers equal-area\',\'lat\',[45.8 47.2],\'long\',[-64.5 -62],\'rect\',\'on\');
m_gshhs_i(\'patch\',[.5 .6 .5]);
m_grid(\'linest\',\'none\',\'linewidth\',2,\'tickdir\',\'out\',\'yaxisloc\',\'right\');
m_text(-64.4,45.9,\'GSHHS\_I (intermediate)\',\'color\',\'m\',\'fontweight\',\'bold\',\'fontsize\',14);
axes(\'position\',[.55 .23 .37 .37]);
m_proj(\'albers equal-area\',\'lat\',[46.375 46.6],\'long\',[-64.2 -63.7],\'rect\',\'on\');
m_gshhs_h(\'patch\',[.6 .6 .6]);
m_grid(\'linest\',\'none\',\'linewidth\',2,\'tickdir\',\'out\',\'xaxisloc\',\'top\',\'yaxisloc\',\'right\');
m_text(-64.18,46.58,\'GSHHS\_H (high)\',\'color\',\'m\',\'fontweight\',\'bold\',\'fontsize\',14);
运行出现这样问题
网站https://www.ngdc.noaa.gov/mgg/shorelines/data/gshhs/oldversions/version1.10/下载
下载version1.10中 gshhs_1.10.zip 文件，将zip文件解压，把gshhs_*.b 文件复制到 matlab\R2009\toolbox\matlab\m_map\@private 文件夹内。打开m_gshhs_c.m，补全路径。
FILNAME=\'D:\MATLAB\toolbox\m_map\private/gshhs_c.b\'
同理将所有的m_gshhs_*.m路径都补充完整。
9、
clear all
m_proj(\'UTM\',\'long\',[-72 -68],\'lat\',[40 44]);
m_gshhs_i(\'color\',\'k\');
m_grid(\'box\',\'fancy\',\'tickdir\',\'in\');
% fake up a trackline
lons=[-71:.1:-67];
lats=60*cos((lons+115)*pi/180);%轨迹经纬度
dates=datenum(1997,10,23,15,1:41,zeros(1,41));%轨迹时间,1分钟到41分钟
m_track(lons,lats,dates,\'ticks\',0,\'times\',4,\'dates\',8,...%刻度间隔，时间间隔4分钟，日期标签间隔8分钟
\'clip\',\'off\',\'color\',\'r\',\'orient\',\'upright\');%颜色标签方向
10、范围环
clear all m_proj(\'hammer\',\'clong\',170);%缩放范围 m_grid(\'xtick\',[],\'ytick\',[],\'linestyle\',\'-\'); m_coast(\'patch\',\'g\'); m_line(100.5,13.5,\'marker\',\'square\',\'color\',\'r\');%增加了一个标记点 m_range_ring(100.5,13.5,[1000:1000:15000],\'color\',\'b\',\'linewi\',2);%范围，颜色，线宽 xlabel(\'1000km range rings from Bangkok\');%曼谷中心1000km的范围环
11、
clear all
bndry_lon=[-128.8 -128.8 -128.3 -128 -126.8 -126.6 -128.8];
bndry_lat=[ 49 50.33 50.33 50 49.5 49 49];
clf;
m_proj(\'lambert\',\'long\',[-130 -121.5],\'lat\',[47 51.5],\'rectbox\',\'on\');
m_gshhs_i(\'color\',\'k\'); % Coastline...
m_gshhs_i(\'speckle\',\'color\',\'r\'); % with speckle added,斑点添加
m_line(bndry_lon,bndry_lat,\'linewi\',2,\'color\',\'k\'); % Area outline ...
m_hatch(bndry_lon,bndry_lat,\'single\',30,5,\'color\',\'k\'); % ...with hatching added.斜线，30度，5个单位相距
m_grid(\'linewi\',2,\'linest\',\'none\',\'tickdir\',\'out\');
title(\'Speckled Boundaries for nice B&W presentation (best in postscript format)\',\'fontsize\',14);
m_text(-128,48,5,{\'Pacific\',\'Ocean\'},\'fontsize\',18);
12
clear all
m_proj(\'miller\',\'lat\',[-75 75]);
m_coast(\'patch\',[.7 1 .7],\'edgecolor\',\'none\');
m_grid(\'box\',\'fancy\',\'linestyle\',\'none\',\'backcolor\',[.9 .9 0.9]);
cities={\'Cairo\',\'Washington\',\'Buenos Aires\'};
lons=[ 30+2/60 -77-2/60 -58-22/60];
lats=[ 31+21/60 38+53/60 -34-45/60];
for k=1:3,
[range,ln,lt]=m_lldist([-123-6/60 lons(k)],[49+13/60 lats(k)],40);
%是两个位置之间分40个点，range返回距离，ln经度矩阵，lt纬度矩阵
m_line(ln,lt,\'color\',\'r\',\'linewi\',2);
m_text(ln(end),lt(end),sprintf(\'%s - %d km\',cities{k},round(range)));
end;
title(\'Great Circle Routes\',\'fontsize\',14,\'fontweight\',\'bold\');
13、 全球/地区温度图
（1）读取数据
clear all setup_nctoolbox %调用工具包 tic %计时 %% nc=ncgeodataset(\'tmpsfc.gdas.199401.grb2\'); %读文件 tem_1=nc.variables %浏览数据类型 %% a1=nc.geovariable(tem_1(1));%取得数据类型为Temperature_surface的数据 b1=a1.data(1,:,:); %第一个时间点温度数据 c1=squeeze(b1)-273.16;%删除单一维度，换为摄氏温度 %% a2=nc.geovariable(tem_1(2));%取得数据类型为lat的数据，纬度 b2=a2.data(:,1)%提取数据 %% a3=nc.geovariable(tem_1(3));%经度 b3=a3.data(:,1)% %% a4=nc.geovariable(tem_1(4));%取得数据类型为time的数据 b4=a4.data(:,1)% %% save A b2 b3 c1 toc
读取的是NCEP/CFSR数据，1994年1月的温度数据。
该数据一共四项。温度，纬度，经度，时间。
时间数据
经度数据
纬度数据
可以看出该数据集的数据组织形式，经度纬度构成世界地图，记录了744个时间点的温度数据。约等于1小时采集一次数据。
温度数据
保存为mat格式数据，为画图做准备。
（2）画图
clear all load A [Plg,Plt]=meshgrid(b3\',b2\');%形成网格 m_proj(\'hammer-aitoff\',\'clongitude\',-150);%投影模式 m_pcolor(Plg,Plt,c1); shading flat; colormap(\'jet\');%颜色选择 hold on; m_pcolor(Plg-360,Plt,c1); shading flat; %着色模式 colormap(\'jet\'); m_coast(); m_grid(\'xaxis\',\'middle\'); % h=colorbar(\'h\'); % set(get(h,\'title\'),\'string\',\'1991年1月全球温度\'); c=colorbar(\'southoutside\',\'fontsize\',12) c.Label.String = \'1994年1月全球温度\'; c.Label.FontSize = 15;
（3）找出最大、最小温度的经纬度
clear all load C Max_col=max(c1);%列最大值 Max_row=max(c1,[],2)%行最大值 Max=max(max(c1)); [x1,y1]=find(c1==max(max(c1)));%x 行，y 列 T_1=Plt(x1,y1)%纬度 T_2=Plg(x1,y1)%经度 Min_col=min(c1);%列最大值 Min_row=min(c1,[],2)%行最大值 Min=min(min(c1)); [x,y]=find(c1==min(min(c1)));%x 行，y 列 T_x=Plt(x,y)%纬度 T_y=Plg(x,y)%经度
可以看出最热52度，在澳大利亚那块（142.8123E，23.2610S）；最冷-62度，在北极圈那块（89.9999E，66.3486N）。
（4）中国（地区）温度图
clear all load A LATLIMS=[3 54]; LONLIMS=[72 134];%选定边界范围 [Plg,Plt]=meshgrid(b3\',b2\');%形成网格 m_proj(\'lambert\',\'lon\',LONLIMS,\'lat\',LATLIMS);%投影模式 m_pcolor(Plg,Plt,c1); shading flat; colormap(\'jet\');%颜色选择 hold on; m_pcolor(Plg-360,Plt,c1); shading flat; %着色模式 colormap(\'jet\'); m_coast(); m_grid(\'box\',\'fancy\',\'tickdir\',\'in\'); % h=colorbar(\'h\'); % set(get(h,\'title\'),\'string\',\'1991年1月全球温度\'); c=colorbar(\'southoutside\',\'fontsize\',12) c.Label.String = \'1994年1月中国温度\'; c.Label.FontSize = 15;
该方法是读取全球数据，只展示部分
（5）
改进的区域方法，读取该区域的数据数据，
clear all setup_nctoolbox %调用工具包 tic %计时 %% min_lat=115; max_lat=279; min_lon=231; max_lon=430; %区域经纬度范围，在数据中的位置 %% nc=ncgeodataset(\'tmp2m.gdas.199401.grb2\'); %读文件 tem_1=nc.variables %浏览数据类型 %% N1=nc.size(tem_1(1));%读取数据大小,可以看出数据的组织形式 a1=nc.geovariable(tem_1(1));%取得数据类型为Temperature_surface的数据 b1=a1.data(1,1,min_lat:max_lat,min_lon:max_lon); %第一个时间点温度数据 c1=squeeze(b1)-273.16;%删除单一维度，换为摄氏温度 %% N2=nc.size(tem_1(2)); a2=nc.geovariable(tem_1(2));%取得数据类型为lat的数据，纬度 b2=a2.data(min_lat:max_lat,1);%提取数据 %% N3=nc.size(tem_1(3)) a3=nc.geovariable(tem_1(3));%经度 b3=a3.data(min_lon:max_lon,1);% %% N4=nc.size(tem_1(4)) a4=nc.geovariable(tem_1(4));%取得数据类型为time的数据 b4=a4.data(:,1);% %% N5=nc.size(tem_1(5));%读取数据大小 a5=nc.geovariable(tem_1(5));%取得数据类型为time的数据 b5=a5.data(:,1);% %% save tem b2 b3 c1 toc
clear all load tem LATLIMS=[3 54]; LONLIMS=[72 134];%选定边界范围 [Plg,Plt]=meshgrid(b3\',b2\');%形成网格 m_proj(\'lambert\',\'lon\',LONLIMS,\'lat\',LATLIMS);%投影模式 m_pcolor(Plg,Plt,c1); shading flat; colormap(\'jet\');%颜色选择 hold on; m_pcolor(Plg-360,Plt,c1); shading flat; %着色模式 colormap(\'jet\'); m_coast(); m_grid(\'box\',\'fancy\',\'tickdir\',\'in\'); % h=colorbar(\'h\'); % set(get(h,\'title\'),\'string\',\'1991年1月全球温度\'); c=colorbar(\'southoutside\',\'fontsize\',12) c.Label.String = \'1994年1月中国温度\'; c.Label.FontSize = 15;
八、圆形区域的画图
1、
clear all LATLIMS=[14 22]; LONLIMS=[108 118];%南海边界范围 m_proj(\'miller\',\'lon\',LONLIMS,\'lat\',LATLIMS);%投影模式 m_coast(\'patch\',[0.1 0.7 0.4]);%绿色填充 m_grid(\'box\',\'fancy\',\'tickdir\',\'in\');%没有网格，边框相间，%m_line(lon,lat,\'linewi\',2,\'color\',\'r\',\'linestyle\',\':\');控制线条格式，点画线还是直线 lon=112:1:116; lat=16:1:20; m_line(lon,lat,\'linewi\',2,\'color\',\'r\');%线宽，2；颜色 [X,Y]=m_ll2xy(117,21); line(X,Y,\'marker\',\'.\',\'markersize\',24\',\'color\',\'r\')%画点
2、距离
clear all %distance用法 arclen=distance([37,0],[35,0])%返回两点间的相对球心的角度，以度为单位 d=distance([37,0],[35,0],6371)% [纬度，经度] [纬度，经度] [半径] D=(arclen/180)*pi*6371 %m_map中函数 dist=m_lldist([0 0],[35 37])%[经度 经度] [纬度 纬度]
3、
clear all
LATLIMS=[14 22];
LONLIMS=[108 118];%南海边界范围
m_proj(\'miller\',\'lon\',LONLIMS,\'lat\',LATLIMS);%投影模式
m_coast(\'patch\',[0.1 0.7 0.4]);%绿色填充
m_grid(\'box\',\'fancy\',\'tickdir\',\'in\');%没有网格，边框相间，%m_line(lon,lat,\'linewi\',2,\'color\',\'r\',\'linestyle\',\':\');控制线条格式，点画线还是直线
load EDH_south_sea_2008
load coordi_south_sea_2008
m_range_ring(114,18,[1e2:1e2:3e2],\'linewi\',2,\'color\',\'b\');%红色300km范围圆圈
% 矩形点阵
range_lat=4:24;%21N和15N对应的位置下标
range_lon=20:40;%111.5E和116.5E对应的下标
for i=1:length(range_lon)
for j=1:length(range_lat)
[X,Y]=m_ll2xy(lon_south_sea(range_lon(i)),lat_south_sea(range_lat(j)));%化为x,y坐标
line(X,Y,\'marker\',\'.\',\'markersize\',10,\'color\',\'k\')%画点
hold on
end
end
%离散圆
[X0,Y0]=m_ll2xy(114,18);%化为x,y坐标
line(X0,Y0,\'marker\',\'.\',\'markersize\',15,\'color\',\'r\');%画圆心
DIST=m_lldist([114 114],[18 19]);%纬度加1度，增加的距离
R=300;%300km
[X1,Y1]=m_ll2xy(114,18+R/DIST);%找到300km的一个点
r=sqrt((X0-X1)^2+(Y0-Y1)^2);%地图距离到图上距离转换
theta=0:pi/180:2*pi;%360度，360个点。
x=X0+r*cos(theta); %（X0,Y0）圆心
y=Y0+r*sin(theta);
plot(x,y,\'.\',\'color\',\'r\')
4、pcolor
clear all n =18; r = (0:n)\'/n; theta = pi*(-n:n)/n; X = r*cos(theta); Y = r*sin(theta); C = r*cos(2*theta); pcolor(X,Y,C) axis equal tight colorbar figure load PCOLOR %南海坐标和波导高度数据 colormap(\'jet\'); shading flat;%平滑方式 gca=pcolor(Plg,Plt,EDH_south_sea) set(gca, \'LineStyle\',\'none\');%去除网格 axis equal tight %按比例展示 colorbar %颜色条
clear all n =6; r = (0:n)\'/n;%0到6，半径上均分的数 theta = pi*(-n:n)/n;%将整个圆分成了13分。 X = r*cos(theta); Y = r*sin(theta); C = r*cos(2*theta); pcolor(X,Y,C) axis equal tight colorbar
