利用MeteoInfo中TrajStat插件做条件轨迹聚类分析

davidyin

    本人一直想做降雨期间水汽的来源轨迹，但是由于时间关系一直做做停停，一直没有取得大的进展。刚开始是用的Hysplit，但是在聚类分析上遇到点问题，后面就停下了。有一天在论坛看到王老师开发的TrajStat模块，试用了一下，发现简单还实用，就试着做了分析，但是聚类分析却无法做出降水期间的聚类分析，比如6月份做轨迹后，是对每天的轨迹合起来做聚类分析，无法就选择出的降水期间的数据进行聚类分析，在跟王老师交流后，终于解决。虽然目前仍无法融会贯通，但是想将自己的一点点经验分享，希望对遇到同样问题的朋友有所帮助。

    软件：Java + MeteoInfo_Java （Java版）+ TrajStat_Plugin ;  TrajStat桌面版 （TrajStat Desktop Application）。（下载地址：http://www.meteothinker.com/Downloads.html）
   
    气象数据来源：GDAS 1-deg 2005年以来的气象数据（这个我是用Hysplit下载的，不知道是否还有其他方式下载）


    操作步骤：参照王老师编写： TrajStat_Help_V1.2

    MeteoInfo中安装TrajStat插件，具体参照 http://bbs.06climate.com/forum.p ... mp;page=1#pid509819

   具体操作如下（操作步骤本贴以桌面版TrajStat为例介绍）：
（1）计算轨迹： 选择TrajStat桌面版菜单栏中Trajectory→Calculation, MeteoInfo中菜单栏TrajStat→Calculate Trajectories，这些具体步骤在Help文件中都有详细说明，照着做就行。但是在这里要注意，选择数据的时候，如要计算6月份的轨迹，一定要把上个月最后几天的数据加进去，保持计算的轨迹长度一致，这个论坛里有帖子解疑，有兴趣可以自行搜索。各种参数设置好后，依次点击Calculate→Convert→Combine，一个月的数据就已经计算好，并已经合并成一个文件，格式为xx.tgs。
（2）将tgs文件转化为shape文件，"Trajectory→convert to shape file"。这一步就已经完成轨迹的分析。
（3）将已经准备好的气象数据，如6月降水数据，保存格式为csv，加入到轨迹文件中。“Trajectory→Add data to Trajectory”, 将6月份降水数据加入到轨迹中，这个可以在view→Attribute Data中可以看到。
（4）这一步在MeteoInfo里面完成。将添加了降水数据的Trajectory文件(.shp格式)添加到MeteoInfo中，选择菜单栏中“选择→通过属性数据选择”，弹出"Select by Attribute"菜单栏，Layer选择加入的数据，Method选择"Create a new selection"，Fields选择加入的气象数据，如precipitation，选择并双击,precipitation将出现在Select query string栏中，然后选择“＞”，并输入0，即选择降水大于0的轨迹，之后点击"Select"，降水大于0的轨迹将在轨迹图中高亮显示。
（5）这一步也在MeteoInfo中完成。将降水大于0的轨迹输出。选择菜单栏中“工具→输出地图数据”，弹出对话框“Output Map Data”，Map Layer选择加入的文件，Output format选择"shape file"，之后点击"Output"，输入新文件名字后保存。
（6）聚类分析，这一步及后面又回到TrajStat桌面版。先将刚才的降水大于0的轨迹文件加入到图层中。选择菜单栏“Trajectory→Convert to line data”，将shape文件转化为lcd文件。这一步完成之后，再选择菜单栏"Trajectory→Cluster Analysis"，弹出对话框"Cluster Calculation"，Input file选择刚才输出的lcd文件，Output file自己命名一个文件名，之后点击"Calculate Cluster"，后面的“Cluster Number”自己选择合适的轨迹条数,如3，之后点击"Calculate Mean Traj"，聚类分析即完成。结果在图层Cluster中显示。
（7）关于聚类分析结果的显示。选择分析的聚类分析结果，点击工具栏“set label”(黄色箭头LB)，出现对话框“Label Set”，Field中选择自己的输出内容，如ratio,Font是设置字体，Shadow color可以选择是否字体加背景。设置好后，点击Add，整个聚类分析就做好了。

    这个是我做的一个结果图：
   

结果图

davidyin

大雪纷飞 发表于 2020-7-9 21:07
您好 请问您的气象数据是如何获取的呢？所获取的气象数据应该是gbl之类的吗？

ftp://arlftp.arlhq.noaa.gov/pub/archives/gdas1/

Mc.Fish

你的后向追踪使用HYSPLIT还是用TrajStat？

宁静的夏天lf

楼主在帖子中说到：（5）这一步也在MeteoInfo中完成。将降水大于0的轨迹输出。选择菜单栏中“工具→输出地图数据”，弹出对话框“Output Map Data”，Map Layer选择加入的文件，Output format选择"shape file"，之后点击"Output"，输入新文件名字后保存。（6）聚类分析，这一步及后面又回到TrajStat桌面版。先将刚才的降水大于0的轨迹文件加入到图层中。选择菜单栏“Trajectory→Convert to line data”，将shape文件转化为lcd文件。我想问的是：对轨迹进行选择后，满足要求的轨迹都高亮加粗显示，然后我根据（5）的操作步骤确实输出了一个.shp的轨迹文件，但是在插件版的TrajStat中怎么加载这个.shp的轨迹文件呢？我现在无法加载这个.shp文件，因此无法创建新聚类分析的图层，接下来没办法做了，烦请解答，谢谢。你的帖子好详细啊，但就是这一步我还没看懂

MeteoInfo

非常详细的经验分享贴，赞一个。

anfly_282

今天注册，路过，问好！

ldepn

不错！{:eb502:}{:eb502:}

菲儿喜欢我

学习了~~~~~~

sharrow

谢谢分享，又弄明白了点问题

蒂小婉

楼主 trajstat里面的三个都要下载吗

davidyin

理论上讲MeteoInfo+Trajstat插件就能解决问题，但是好像现在这个版本的Trajstat插件中有几个功能没有。具体可以咨询一下王老师。

paulwei

学习中