matlab画地图的工具：worldmap和m_map

*用MATLAB画世界地图并加散点

clear;clc;
worldmap([-55 -30],[-110 -60]);
h1=axesm('MapProjection','eqdcylin','maplatlimit',[-58 -30],'maplonlimit',[-110 -60],'frame','on','parallellabel','on','meridianlabel','on','flinewidth',1,'plabellocation',5,'mlabellocation',10);
setm(h1,'fedgecolor',[.5 .6 .6], 'fontname','Arial','fontsize',7);
grid off;
h=geoshow('landareas.shp');
x=

y=

scatterm(y,x,'filled','k')

加等值线

clear;clc;
A=xlsread('2010作图2','5','N2:P26');
[x,y]=meshgrid(-106.5:-81.5,-46.5:-35.5);
z=abs(griddata(A(:,1),A(:,2),A(:,3),x,y,'v4'));
subplot(2,2,1)%figure(1)
worldmap([-55 -30],[-110 -60]);
h1=axesm('MapProjection','eqdcylin','maplatlimit',[-58 -30],'maplonlimit',[-110 -60],'frame','on','parallellabel','on','meridianlabel','on','flinewidth',1,'plabellocation',5,'mlabellocation',10);
setm(h1,'fedgecolor',[.5 .6 .6], 'fontname','Arial','fontsize',7);
grid off;
h=geoshow('landareas.shp');
scatterm(A(:,2),A(:,1),'filled','k')
[c,h] =contourm(y,x,z,[0.6,0.8,0.9],'k-');  %等值线
hold on;[cc,hh]=contourm(y,x,z,[0.1,0.2,0.3,0.4,0.5,0.7],'k-');
clabel(c,h)
title('5月')

1、用matlab画世界地图

[matlab自带的例子]

ax = worldmap('World');
setm(ax, 'Origin', [0 180 0])
land = shaperead('landareas', 'UseGeoCoords', true);
geoshow(ax, land, 'FaceColor', [0.5 0.7 0.5])
lakes = shaperead('worldlakes', 'UseGeoCoords', true);
geoshow(lakes, 'FaceColor', 'blue')
rivers = shaperead('worldrivers', 'UseGeoCoords', true);
geoshow(rivers, 'Color', 'blue')
cities = shaperead('worldcities', 'UseGeoCoords', true);
geoshow(cities, 'Marker', '.', 'Color', 'red')

http://s7/bmiddle/4cf8aad3471811f31c1c6&690 

2、matlab的m_map工具箱及添加行政边界底图

from http://blog.163.com/ahualian@126/blog/static/2350491720093238409154/

对习惯使用matlab的人来说，m_map是一个很好的绘制地图的免费工具箱。可以选择的投影种类近20种，包括常用的Lambert、Mercator、UTM等。可以测量距离（ m_lldist， m_xydist ），绘制等值线（m_contour），等值线填充图（m_contourf），矢量图（m_quiver），栅格图（m_pcolor）等，并与相应的matlab函数语法类似，很容易使用。

      m_map通过m_coast提供1/4 degree分辨率的全球海岸线，通过下载GSHHS可以得到更高分辨率的海岸线数据。

     除海岸线外，还可以利用已有的GIS行政边界资料，在地图上添加行政边界底图。具体作法为：下载.shp格式的行政边界文件（ 国家基础地理信息系统的下载服务），将下载的.shp文件通过mapinfo转换成.dxf文件，利用已编译好的fortran程序，读取经纬度信息，输出.dat文件。导入matlab空间，可以直接m_plot，也可以调用m_plotbndry（）。

通过以上方法可以实现.shp在matlab绘图中的应用。

m_map工具箱的下载和介绍

fortran程序及.mat文件下载

3、[原创ZHOU Feng]在matlab中利用worldmap画中国区域图时加上台湾和钓鱼岛

from http://www.sciencenet.cn/m/Print.aspx?id=210235

Matlab是我们常用的一个画图和计算、仿真工具，在我们海洋科研中，经常在画图时需要加上底图（譬如海岸线，国界、省界线等等）。常用的一个工具包是M_map。但这里我要讲的是利用matlab自带的一个画图工具包话底图，这个工具包就是worldmap.
http://www.sciencenet.cn/upload/blog/images/2009/1/2009115132755162.jpg

    worldmap的一般用法约为：
   >> figure; worldmap('china'); polcmap;
    或者
    >> figure; worldmap china; polcmap;
    如果要加入颜色的画，一般可以这样：
    >> figure
        worldmap('china','patch')
        scaleruler

    这时候问题就出来了。因为鬼子偷偷的讲台湾和大陆用两种颜色表示；甚至，如果你放大图片的画，会发现钓鱼岛也是不同颜色。这个是我们不能忍受的！！！！！
    于是我就费了一点时间，琢磨这个画图，然后修改。下面是这个脚本程序，在matlab中执行就可，这时候台湾和钓鱼岛和大陆就是一个颜色了^_^。


    用兴趣的朋友可以把这个用法举一反三（(ZHOU Feng)zhoufeng@sio.org.cn。matlab版本是6.5）。
不过老实说，我不太用matlab自带的这个画图包，用M_map比较多一点。试验一下，感觉还可以用用。
画上述图的代码如下：

% 把台湾和大陆合成一个文件保存起来，这样画图用patch就是一种颜色（Zhou Feng, 2008-06-30, SOED, Hangzhou）。
%
% by ZHOU Feng
% zhoufeng@sio.org.cn
% SOED, 2nd Institute of Oceanography
% 2008-06-30

s1 = worldhi('china');
s2 = worldhi('taiwan');
disp(s1);

% add Taiwan together (ZHOU Feng)zhoufeng@sio.org.cn
s = s1;
s.lat = [s1.lat;NaN; s2.lat];
s.long= [s1.long; NaN;s2.long];

% add the Diaoyu Island (钓鱼岛)
% 钓鱼岛群岛由钓鱼岛、黄尾岛、赤尾岛、南小岛、北小岛、大南小岛、大北小岛和飞濑岛等岛屿组成，总面积约7平方公里。
% 地理位置：东经123°-124°34′北纬25°40′-26°。
%
% ---这里的分辨率只有两块 --
s3 = worldhi('japan');  %(ZHOU Feng)zhoufeng@sio.org.cn
x = s3.long;
y = s3.lat;
idx = find(x>123.0 & x<124.5);
idy = find(y> 25.5 & y< 26.0);
m = length(idy);
id=[];
for i=1:m
    tmp=find(idx == idy(i));
    if isempty(tmp)
    else
        id = [id; idx(tmp)];    %(ZHOU Feng)zhoufeng@sio.org.cn
    end
end
%%longd = x(id);
%%latd  = y(id);
%
% find nan
dtmp = find(diff(id)>1);
if isempty(dtmp)
    disp('no change')
elseif length(dtmp)==1
    idnew = [id(1:dtmp); id(dtmp+1)-1;id(dtmp+1:end)];
else
    for j=1:length(dtmp)
        idnew = [id(1:dtmp(j)+j-1); id(dtmp(j)+1)+j-1; id(dtmp(j)+j+1:end)];        % 未试验，可能有误
    end
end
longd = x(idnew);
latd  = y(idnew);

s.lat = [s.lat; NaN; latd];
s.long= [s.long; NaN; longd];

worldmap china
h = displaym(s);
polcmap

OK!!!!!!!!!!!!!!!