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