用Trajstat聚类轨迹之后的shapefile文件

zhangyanlin

想读取一下Trajstat聚类分析之后几种平均大气轨迹的坐标、高度等信息，但是现在无法打开shapefile文件。是否有办法将Trajstat聚类轨迹之后的shapefile文件比如xxx.shp 导出成TGS或者CSV.TXT 类型文件？

MeteoInfo

参考此贴：MeteoInfoLab脚本示例：获取气团轨迹每个节点的气象数据
http://bbs.06climate.com/forum.p ... 463&fromuid=106
(出处: 气象家园)

打开shape文件，遍历每个图元（即每条聚类轨迹），获取每条轨迹上每个点的经度、纬度、高度信息，写成文本文件可以自己加些代码完成，具体可以在网上找Python写文本文件的例子学习。

zhangyanlin

谢谢王老师回复。如果只是实现读取轨迹的经纬度和高度信息，如何实现？王老师能帮忙写个脚本吗？

MeteoInfo

zhangyanlin 发表于 2015-10-27 09:02
谢谢王老师回复。如果只是实现读取轨迹的经纬度和高度信息，如何实现？王老师能帮忙写个脚本吗？

你自己先尝试着做吧，我太忙了。

zhangyanlin

王老师运行之后出现以下错误
NameError: name 'addfile_hytraj' is not defined

MeteoInfo

zhangyanlin 发表于 2015-10-27 13:21
王老师运行之后出现以下错误
NameError: name 'addfile_hytraj' is not defined

你是怎么运行的？

zhangyanlin

按照您给的脚本，将工作文件夹路径进行了替代。
出现了以下两条错误

NameError: name 'os' is not defined；增加i“mport os” 之后 出现 NameError: name 'addfile_hytraj' is not defined

#-----------------------------------------------------
# Author: Yaqiang Wang
# Date: 2015-9-30
# Purpose: Get meteorological data along trajectory
# Note: Sample
#-----------------------------------------------------

# Set working directory
trajDir = 'D:/Temp/HYSPLIT'
meteoDir = 'D:/Temp/arl'

# Open trjactory data file
print 'Open trajectory data file ...'
trajfn = os.path.join(trajDir, 'traj_20090731')
print 'Trajectory file: ' + trajfn
trajf = addfile_hytraj(trajfn)
# Create trajectory layer
trajLayer = trajf.trajlayer()

MeteoInfo

zhangyanlin 发表于 2015-10-27 13:51
按照您给的脚本，将工作文件夹路径进行了替代。
出现了以下两条错误

要在MeteoInfoLab里运行。

zhangyanlin

我的轨迹名称为1304.shp.
运行之后报错：

Trajectory file: C:/HYSPLIT data\1304
File not exist: C:/HYSPLIT data\1304

MeteoInfo

.shp也需要