MeteoInfo脚本示例：移动轨迹

MeteoInfo

江志红有一篇文章分析了2007年淮河流域强降水过程的水汽输送特征，利用了气团轨迹聚类分析技术。该文中轨迹起点不同，水平方向15个格点（2.5*2.5度），垂直方向3个高度。用了HYSPLIT自带的聚类分析工具进行聚类。其实不管HYSPLIT还是TrajStat中的聚类分析方法都要求所有参与聚类的轨迹都是同一起点（高度没有关系，目前的聚类算法中都没有考虑高度）。那么文中的聚类分析必然需要首先将不同起点的轨迹平移至同一起点（应该是研究区的中点），文章并没有相关的说明。不过从文章中图8b、10b来看5类平均轨迹应该是同一起点（图也不是很清楚）。

这种平移必定对之后的聚类分析有所影响，但从解决问题的角度来说也是可行的，可以对一个相对较大的区域的气团轨迹特征进行分析。因此写了一个脚本来实现气团轨迹的平移（假定气团轨迹已经转为shape文件，并加载为MeteoInfo的一个图层）：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2021-9-14
   
4. # Purpose: Move trajectory to a new start point
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. import os
   
8. from org.meteoinfo.geo.layer import VectorLayer
   
9. from org.meteoinfo.geometry.shape import ShapeTypes, PointZ, PolylineZShape
   
10. from org.meteoinfo.geometry.legend import LegendType
    
11. 
    
12. # Get trajectory layer
    
13. print 'Start...'
    
14. trajLayer = miapp.getMapView().getLayer('JiuquanSpr.shp')
    
15. 
    
16. #Create result trajectory layer
    
17. print 'Create result trajectry layer...'
    
18. layer = VectorLayer(trajLayer.getShapeType())
    
19. layer.setLayerName('Moved_trajectories')
    
20. for field in trajLayer.getFields():
    
21.     layer.editAddField(field.getColumnName(), field.getDataType())
    
22. 
    
23. #Move trajectories to new start position
    
24. print 'Move trajectories...'
    
25. newX = 105.0
    
26. newY = 40.5
    
27. i = 0
    
28. for shape in trajLayer.getShapes():
    
29.     point = shape.getPoints().get(0)
    
30.     xd = newX - point.X
    
31.     yd = newY - point.Y
    
32.     npoints = []
    
33.     for point in shape.getPoints():
    
34.         npoint = PointZ()
    
35.         npoint.X = point.X + xd
    
36.         npoint.Y = point.Y + yd
    
37.         npoint.M = point.M
    
38.         npoint.Z = point.Z
    
39.         npoints.append(npoint)
    
40.     polyline = PolylineZShape()
    
41.     polyline.setPoints(npoints)
    
42.     sn = layer.getShapeNum()
    
43.     layer.editInsertShape(polyline, sn)
    
44.     for j in range(0, layer.getFieldNumber()):
    
45.         layer.editCellValue(j, sn, trajLayer.getCellValue(j, i))
    
46.     i += 1
    
47. layer.updateLegendScheme(LegendType.SINGLE_SYMBOL, 'TrajID')
    
48. 
    
49. #Add the moved layer
    
50. print 'Add the layer...'
    
51. miapp.getMapDocument().getActiveMapFrame().addLayer(layer)
    
52. 
    
53. #Save the moved layer
    
54. print 'Save the moved layer...'
    
55. outDir = 'D:/Temp/test'
    
56. fn = os.path.join(outDir, 'traj_moved.shp')
    
57. #layer.saveFile(fn)
    
58. 
    
59. print 'Finish !'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2021-9-14<br /> # Purpose: Move trajectory to a new start point<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> import os<br /> from org.meteoinfo.geo.layer import VectorLayer<br /> from org.meteoinfo.geometry.shape import ShapeTypes, PointZ, PolylineZShape<br /> from org.meteoinfo.geometry.legend import LegendType<br /> <br /> # Get trajectory layer<br /> print 'Start...'<br /> trajLayer = miapp.getMapView().getLayer('JiuquanSpr.shp')<br /> <br /> #Create result trajectory layer<br /> print 'Create result trajectry layer...'<br /> layer = VectorLayer(trajLayer.getShapeType())<br /> layer.setLayerName('Moved_trajectories')<br /> for field in trajLayer.getFields():<br />     layer.editAddField(field.getColumnName(), field.getDataType())<br /> <br /> #Move trajectories to new start position<br /> print 'Move trajectories...'<br /> newX = 105.0<br /> newY = 40.5<br /> i = 0<br /> for shape in trajLayer.getShapes():<br />     point = shape.getPoints().get(0)<br />     xd = newX - point.X<br />     yd = newY - point.Y<br />     npoints = []<br />     for point in shape.getPoints():<br />         npoint = PointZ()<br />         npoint.X = point.X + xd<br />         npoint.Y = point.Y + yd<br />         npoint.M = point.M<br />         npoint.Z = point.Z<br />         npoints.append(npoint)<br />     polyline = PolylineZShape()<br />     polyline.setPoints(npoints)<br />     sn = layer.getShapeNum()<br />     layer.editInsertShape(polyline, sn)<br />     for j in range(0, layer.getFieldNumber()):<br />         layer.editCellValue(j, sn, trajLayer.getCellValue(j, i))<br />     i += 1<br /> layer.updateLegendScheme(LegendType.SINGLE_SYMBOL, 'TrajID')<br /> <br /> #Add the moved layer<br /> print 'Add the layer...'<br /> miapp.getMapDocument().getActiveMapFrame().addLayer(layer)<br /> <br /> #Save the moved layer<br /> print 'Save the moved layer...'<br /> outDir = 'D:/Temp/test'<br /> fn = os.path.join(outDir, 'traj_moved.shp')<br /> #layer.saveFile(fn)<br /> <br /> print 'Finish !'

tan8935079

王老师气团轨迹怎么转成shape文件呢？还有后边的加载为Meteoinfo的一个图层也看不大懂哦。新人请见谅啊

MeteoInfo

tan8935079 发表于 2015-1-5 18:16
王老师气团轨迹怎么转成shape文件呢？还有后边的加载为Meteoinfo的一个图层也看不大懂哦。新人请见谅啊{:eb ...

代码里把
#layer.saveFile(fn)
前面的#去掉即可保存图层为一个shape文件。

rceclx

王老师分析的好到位！

MeteoInfo

tan8935079 发表于 2015-1-5 18:16
王老师气团轨迹怎么转成shape文件呢？还有后边的加载为Meteoinfo的一个图层也看不大懂哦。新人请见谅啊{:eb ...

TrajStat插件里有将气团轨迹转为shape文件的功能

游子

同一起点？我怎么觉得是同一终点呢，轨迹是向前推的吧？

MeteoInfo

游子 发表于 2015-1-6 08:44
同一起点？我怎么觉得是同一终点呢，轨迹是向前推的吧？

轨迹可以向前推也可以向后推，即所谓的前向轨迹和后向轨迹。

即便是向前推，也可以理解为向前推的起点，很难理解吗？

yuejun0026

感谢王老师的解决

付亚男

王老师，我新下载的TrajStat对不同起点的轨迹实现了聚类，是否是新版本增加的功能？

MeteoInfo

付亚男 发表于 2015-9-2 21:22
王老师，我新下载的TrajStat对不同起点的轨迹实现了聚类，是否是新版本增加的功能？

没有，只要轨迹节点数相同，总能算出聚类结果，但是结果是否合理就不得而知了。