MeteoInfoLab脚本示例：轨迹节点数

MeteoInfo

在TrajStat中进行轨迹聚类要求每条轨迹都有同样的节点数（某些轨迹可能因为气象数据等原因没有计算全），但是TrajStat中没有现成的节点数检查功能，暂时也懒得增加。可以通过写一个MeteoInfoLab小脚本来进行轨迹节点数检查。

下面的脚本打印出每条轨迹的节点数：

1. # Open trjactory shape file
   
2. print 'Open trajectory shape file ...'
   
3. trajfn = 'D:/Temp/traj/Sample/JiuquanSpr.shp'
   
4. trajLayer = shaperead(trajfn)
   
5. 
   
6. # Chek endpoint number of each trajectory
   
7. print 'Chek endpoint number of each trajectory...'
   
8. idx = 0
   
9. for tline in trajLayer.shapes():
   
10.     t = trajLayer.cellvalue('Date', idx)
    
11.     h = trajLayer.cellvalue('Hour', idx)   
    
12.     t.replace(hour=h)
    
13.     print t.strftime('%Y-%m-%d %H') + ': %i' % tline.getPointNum()
    
14.     idx += 1
    
15. 
    
16. print 'Finish...'

# Open trjactory shape file<br /> print 'Open trajectory shape file ...'<br /> trajfn = 'D:/Temp/traj/Sample/JiuquanSpr.shp'<br /> trajLayer = shaperead(trajfn)<br /> <br /> # Chek endpoint number of each trajectory<br /> print 'Chek endpoint number of each trajectory...'<br /> idx = 0<br /> for tline in trajLayer.shapes():<br />     t = trajLayer.cellvalue('Date', idx)<br />     h = trajLayer.cellvalue('Hour', idx)    <br />     t.replace(hour=h)<br />     print t.strftime('%Y-%m-%d %H') + ': %i' % tline.getPointNum()<br />     idx += 1<br /> <br /> print 'Finish...'

也可以设定一个“正确”的节点数，检查每条轨迹的节点数是否正确，打印出有问题的轨迹：

1. # Open trjactory shape file
   
2. print 'Open trajectory shape file ...'
   
3. trajfn = 'D:/Temp/traj/Sample/JiuquanSpr.shp'
   
4. trajLayer = shaperead(trajfn)
   
5. 
   
6. # Chek endpoint number of each trajectory
   
7. print 'Chek endpoint number of each trajectory...'
   
8. # Set a correct number and print the trajectories with more or less endpoints
   
9. epn = 73
   
10. idx = 0
    
11. n = 0
    
12. for tline in trajLayer.shapes():
    
13.     if tline.getPointNum() != 73:
    
14.         t = trajLayer.cellvalue('Date', idx)
    
15.         h = trajLayer.cellvalue('Hour', idx)   
    
16.         t.replace(hour=h)
    
17.         print t.strftime('%Y-%m-%d %H') + ': %i' % tline.getPointNum()
    
18.         n += 1
    
19.     idx += 1
    
20. 
    
21. if n == 0:
    
22.     print 'All trajectories have same endpoint number!'

# Open trjactory shape file<br /> print 'Open trajectory shape file ...'<br /> trajfn = 'D:/Temp/traj/Sample/JiuquanSpr.shp'<br /> trajLayer = shaperead(trajfn)<br /> <br /> # Chek endpoint number of each trajectory<br /> print 'Chek endpoint number of each trajectory...'<br /> # Set a correct number and print the trajectories with more or less endpoints<br /> epn = 73<br /> idx = 0<br /> n = 0<br /> for tline in trajLayer.shapes():<br />     if tline.getPointNum() != 73:<br />         t = trajLayer.cellvalue('Date', idx)<br />         h = trajLayer.cellvalue('Hour', idx)    <br />         t.replace(hour=h)<br />         print t.strftime('%Y-%m-%d %H') + ': %i' % tline.getPointNum()<br />         n += 1<br />     idx += 1<br /> <br /> if n == 0:<br />     print 'All trajectories have same endpoint number!'

wuxy17

王老师，上午跟您讲聚类分析出现问题，您说可以试试查找节点的方式，我按照您的代码试了一下，出了一点点问题，不好意思又麻烦您了，我没有学习过metoinfo的脚本语言，所以不太懂，请您指教一下

wuxy17

王老师，刚加入了QQ群，代码错误问题已解决，但是聚类分析时尽管轨迹集合中各轨迹均为73节点，还是出现问题了

nayaxin

wuxy17 发表于 2017-11-13 18:02
王老师，刚加入了QQ群，代码错误问题已解决，但是聚类分析时尽管轨迹集合中各轨迹均为73节点，还是出现问题 ...

您好，我也遇到了聚类分析时JAVA报错的问题 请问您是如何解决的呢

txxrk

请问，按上述方法知道节点数出错后，怎么处理数据？

leileilei

wuxy17 发表于 2017-11-13 18:02
王老师，刚加入了QQ群，代码错误问题已解决，但是聚类分析时尽管轨迹集合中各轨迹均为73节点，还是出现问题 ...

您好，我在calculate的时候也出现了类似的问题，请问您怎么解决的呢？

古月睢风

wuxy17 发表于 2017-11-13 18:02
王老师，刚加入了QQ群，代码错误问题已解决，但是聚类分析时尽管轨迹集合中各轨迹均为73节点，还是出现问题 ...

您好，可以分享一下qq群吗，我在聚类的时候显示轨迹数小于最大簇数，不知道该怎么解决，想求教

shoegaze

txxrk 发表于 2019-8-1 17:16
请问，按上述方法知道节点数出错后，怎么处理数据？

同问，知道节点数出错后，该如何处理数据呢？？？题主知道如何解决了吗？？？