登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 小寶 于 2013-10-28 13:11 编辑
我以前是用IDL编写AOT数据,不过发现实现中国境内数据很难,而且一些坐标设置很麻烦,比如字体等,现在用MATLAB编写,查看了论坛中的几个帖子,因为才看MATLAB,所以有些函数不太明白。
我这里还是遇到麻烦所以请大家帮忙给我看看。不多说上图和代码
- clc;clear
- s='F:\studyData\MODIS\Terra\L3MonthyAerosol\2004\032\MOD08_M3.A2004032.051.2010284172747.hdf';
- l333=hdfread(s,'Optical_Depth_Land_And_Ocean_Mean_Mean');
- [alon,alat]=meshgrid(72:1:137,16:1:55);
- l1=double(l333(90-55:90-16,72+180:137+180));
- in = importdata('matlabchinain.mat');
- map_path = shaperead('bou1_4p.shp') ;
- map_X = [map_path(:).X];
- map_Y = [map_path(:).Y];
- % in(:,:)=inpolygon(alon,alat,map_X,map_Y);
- china=l1.*flipud(in);
- 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的数据显示成白色
|