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')

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

对习惯使用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绘图中的应用。

3、[原创ZHOU

Feng]在matlab中利用worldmap画中国区域图时加上台湾和钓鱼岛

Matlab是我们常用的一个画图和计算、仿真工具，在我们海洋科研中，经常在画图时需要加上底图(譬如海岸线，国界、省界线等等)。常用的一个工具包是M_map。但这里我要讲的是利用matlab自带的一个画图工具包话底图，这个工具包就是worldmap.

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!!!!!!!!!!!!!!!