怎么样用matlab在地图中添加地理数据

houxiaogang

大家好，我有一个如附件所示的矩阵，数据表示某气体浓度，其中左上角起始经纬度已知（112.5,41.5），我怎么做出如附件（样例图）所示的图，请高手指教，以下代码是我自己写的，成图结果如附件（结果图）所示。主要问题是：1、我自己做的数据没有样例图好看，也就是分布不均匀；2、我做的colobar也不好看。以下是我编写的部分代码：
    figure('Position',[1 1 scrsz(3) scrsz(4)]);
    contourf(x,y,data,300);
     shading flat;
    subplot(1,1,1);
    sh1=shaperead('G:\pro\国界线.shp','UseGeoCoords',true);%显示矢量数据
    load coast;
    hold on;
    plot(long,lat,'-k','LineWidth',1);
    hold off;
    % colorbar('location','YTickLabel',...
    %     {'Freezing','Cold','Cool','Neutral',...
    %      'Warm','Hot','Burning','Nuclear'});
    %  colorbar('YTickLabel',...
    %     {'Freezing','Cold','Cool','Neutral',...
    %      'Warm','Hot','Burning','Nuclear'});
    %colorbar('West','Position',[0.88,0.1,0.02,0.8]);
     caxis([0,1]);
     colorbar('West','Position',[0.88,0.1,0.02,0.8],'YTickLabel',...
           {'0.1',' ','0.26',' ','0.42',' ','0.58',' ','0.74',' ','0.9'});
    set(gca,'GridLineStyle', '-','LineWidth',2,'FontSize',12,'Ylim',[33 41],'Xlim',[113 122],'Position',[0.05,0.09,0.8,0.8]...
         ,'XTick',[113:122],'XTicklabel',{'113E','114E','115E','116E','117E','118E','119E','120E','121E','122E'}...
         ,'YTick',[33:41],'YTicklabel',{'33N','34N','35N','36N','37N','38N','39N','40N','41N'});

二爷名声在外

我有个建议就是楼主换GrADS吧，分分钟完爆Matlab……

MeteoInfo

houxiaogang 发表于 2015-6-17 22:37
非常感谢二楼的帮助，是不是这种图只能用MeteoInfoLab做？

GrADS也能做。其实这里面关键的技术点有两个：格点加密插值和有大量无效值格点数据等值线分析。格点加密插值这个应该什么软件都能做，但是有大量无效值格点数据等值线分析可能不是什么软件都能实现的。我开发的wContour库可以做这个，相关算法的论文也在Computers & Geosciences上发表（http://www.meteothinker.com/Documentation.html），MeteoInfo中的等值线分析功能是基于wContour库的。MatLab我没怎么用过，更多时候是查查MatLab的函数用法供MeteoInfoLab开发参考。不过这么牛的软件做不好这种格点数据的等值线分析有些不能让人相信，也许有MatLab大牛知道怎么正确处理。

沙颖凯

MOD08_D3好象是需要投影变换的，然后海岸线精度和地图投影没做好。

MeteoInfo

貌似MatLab对于有大量无效数据的格点追踪等值线的算法有问题。MeteoInfoLab的结果参考此帖：MeteoInfoLab脚本示例 - AOD
http://bbs.06climate.com/forum.p ... 140&fromuid=106
(出处: 气象家园)

houxiaogang

非常感谢二楼的帮助，是不是这种图只能用MeteoInfoLab做？

houxiaogang

MeteoInfo 发表于 2015-6-18 09:05
GrADS也能做。其实这里面关键的技术点有两个：格点加密插值和有大量无效值格点数据等值线分析。格点加密 ...

谢谢王老师，今天看了你的空间，才知道你是MeteoInfoLab专家。我今天安装了MeteoInfoLab，借鉴你上面给出的例子，运行有错，不能解决，希望得到你的指点！

houxiaogang

MeteoInfo 发表于 2015-6-18 09:05
GrADS也能做。其实这里面关键的技术点有两个：格点加密插值和有大量无效值格点数据等值线分析。格点加密 ...

谢谢王老师，今天看了你的空间，才知道你是MeteoInfoLab专家。我今天安装了MeteoInfoLab，借鉴你上面给出的例子，运行有错，不能解决，希望得到你的指点！

houxiaogang

二爷名声在外 发表于 2015-6-17 16:14
我有个建议就是楼主换GrADS吧，分分钟完爆Matlab……

大侠，能具体点吗？

MeteoInfo

houxiaogang 发表于 2015-6-18 23:36
谢谢王老师，今天看了你的空间，才知道你是MeteoInfoLab专家。我今天安装了MeteoInfoLab，借鉴你上面给出 ...

下载最新的MeteoInfo_Java_1.2.5R1_Files.zip文件，可以在我的百度网盘上下载：http://yun.baidu.com/share/link?shareid=669776748&uk=51062435