爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 12318|回复: 18

[讨论] 只画中国境内数据(大家讨论下)

[复制链接]
发表于 2013-10-28 12:58:32 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
本帖最后由 小寶 于 2013-10-28 13:11 编辑

我以前是用IDL编写AOT数据,不过发现实现中国境内数据很难,而且一些坐标设置很麻烦,比如字体等,现在用MATLAB编写,查看了论坛中的几个帖子,因为才看MATLAB,所以有些函数不太明白。
我这里还是遇到麻烦所以请大家帮忙给我看看。不多说上图和代码
~RLBH7PZ(WO(E`ZIJ4ITCA7.jpg
  1. clc;clear
  2. s='F:\studyData\MODIS\Terra\L3MonthyAerosol\2004\032\MOD08_M3.A2004032.051.2010284172747.hdf';
  3. l333=hdfread(s,'Optical_Depth_Land_And_Ocean_Mean_Mean');
  4. [alon,alat]=meshgrid(72:1:137,16:1:55);
  5. l1=double(l333(90-55:90-16,72+180:137+180));  
  6. in = importdata('matlabchinain.mat');
  7. map_path = shaperead('bou1_4p.shp') ;
  8. map_X = [map_path(:).X];
  9. map_Y = [map_path(:).Y];  
  10. % in(:,:)=inpolygon(alon,alat,map_X,map_Y);
  11. china=l1.*flipud(in);
  12. imagesc(china,[0 1000]);
复制代码

到这里介绍下数据,l333是一个180*360的数组,in = importdata('matlabchinain.mat'); 就是下面% in(:,:)=inpolygon(alon,alat,map_X,map_Y);注释掉的这个in因为很耗时所以直接保存in。

in也是数组,在中国境内是1其他地方是0



imagesc显示出来了图像,效果我是比较满意的,但是还有问题请教大家:
1  这时候怎么加入中国的国界我试了几次还是不行。地图和上面的一样,南海暂时不需要或者大家有什么好的显示图的方式也告诉我下,帮帮忙
2  如何把等于0的数据显示成白色











密码修改失败请联系微信:mofangbao
发表于 2013-10-28 22:21:52 | 显示全部楼层
用MeteoInfo Java版尝试了一下MOD08_M3数据,HDF4-EOS格式,打开没问题。

Image00046.png
Image00047.png

密码修改失败请联系微信:mofangbao
 楼主| 发表于 2013-10-29 09:10:52 | 显示全部楼层
密码修改失败请联系微信:mofangbao
发表于 2013-10-29 09:13:00 | 显示全部楼层
小寶 发表于 2013-10-29 09:10
能把你那两个文件shp发给我么谢谢

你下载MeteoInfo软件里面就有。
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2013-10-29 09:19:57 | 显示全部楼层
MeteoInfo 发表于 2013-10-29 09:13
你下载MeteoInfo软件里面就有。

哦正在下试试这个软件呵呵,貌似很强大呢
密码修改失败请联系微信:mofangbao
发表于 2013-10-29 09:25:57 | 显示全部楼层
小寶 发表于 2013-10-29 09:19
哦正在下试试这个软件呵呵,貌似很强大呢

需要注意的是HDF4-EOS格式只在MeteoInfo Java版支持,C#版不支持。
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2013-10-29 09:28:27 | 显示全部楼层
MeteoInfo 发表于 2013-10-29 09:25
需要注意的是HDF4-EOS格式只在MeteoInfo Java版支持,C#版不支持。

貌似下错了,java版在? 我下的貌似是c的我说没找到呢。。。
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2013-10-29 09:34:57 | 显示全部楼层
MeteoInfo 发表于 2013-10-29 09:25
需要注意的是HDF4-EOS格式只在MeteoInfo Java版支持,C#版不支持。

哦找到了 在你们官网下呢
密码修改失败请联系微信:mofangbao
发表于 2013-10-29 09:42:35 | 显示全部楼层
(1)这个需要你先把中国地图画出来(网上应该有方法),然后用hold on命令;接着再用imagesc就可以解决了。
(2)这个只需要调整色标。用colcormap editor把小值对应的区域调整为白色即可
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2013-10-29 10:57:14 | 显示全部楼层
MeteoInfo 发表于 2013-10-28 22:21
用MeteoInfo Java版尝试了一下MOD08_M3数据,HDF4-EOS格式,打开没问题。

有这个图相关设置的帮助么? 我弄半天没弄出来只显示中国的。。。。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表