芝麻开门分享 http://blog.sciencenet.cn/u/opensesame 周锋,研究员

博文

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

已有 9798 次阅读 2009-1-15 13:28 |个人分类:软件技巧|系统分类:科研笔记| 中国, 台湾, 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!!!!!!!!!!!!!!!

http://blog.sciencenet.cn/blog-113657-210235.html

上一篇:[转载]人类活动加剧近海水体缺氧
下一篇:[转发]卫星海洋环境动力学国家重点实验室基金申请

1 王铮

发表评论 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2021-2-28 10:07

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部