大气成分统计分析利器-气团轨迹统计分析软件-TrajStat_1.2.2.2下载【9月23日更新】 - 气象人自己的资料站
&&|&&&&|&&&&|&&&&|&&&&|&&&&|&&&&|&&&&|&&&&|&& |&&
您的位置: &>&&>&
大气成分统计分析利器-气团轨迹统计分析软件-TrajStat_1.2.2.2下载【9月23日更新】
软件大小：10.56 MB
软件类型：
软件语言：简体中文
软件授权：共享软件
软件提交：mofangbao
更新时间： 13:18:00
软件作者：
软件等级：
官方网址：
程序演示：
运行环境：/all
下载统计：总下载：0，本月下载：0，本周下载：0，今日下载：0
大气成分统计分析利器-气团轨迹统计分析软件最新版-TrajStat1.2.2R下载
该软件的主要功能有：
计算（使用NOAA HYSPLIT模式中气团轨迹计算模块）、显示、查询批量的气团轨迹，进行轨迹聚类分析，可以将观测数据加入轨迹文件中，从而分析观测站点大气成分路径及来源（PSCF和CWT方法）。
下面这篇文章包含了软件的介绍：Wang, Y.Q., Zhang, X.Y. and Draxler, R., 2009. TrajStat: GIS-based software that uses various trajectory statistical analysis methods to identify potential sources from long-term air pollution measurement data. Environmental Modelling & Software, 24: 938-939.
该软件为英文界面，帮助文件也为英文
按字母检索
Powered by EmpireDown 2.5工具类服务
编辑部专用服务
作者专用服务
对聚类及聚类评价若干问题的研究
聚类是数据挖掘中的一项重要任务，在很多应用中起着重要的作用。现在已经有很多成熟的聚类方法，每种方法的适用范围不尽相同。对于同样的数据，不同聚类方法得到的结果也不尽相同。为了在具体的应用中选择最恰当的聚类方法，人们从不同的角度提出了一些聚类评价方法。然而不同的评价方法对于同一结果的评价并不完全一致，每种评价方法都有自己的适用范围。超出了方法本身的适用范围将会导致方法的失效。如何保证聚类及聚类评价的有效性成为聚类应用的一个关键问题。本文从聚类的基本概念出发对聚类及聚类评价中的相关问题进行了系统的研究，主要包括以下几方面内容：　　
一、提出了聚类相似性的一般形式。根据聚类的概念，聚类是一个依据给定相似性对数据集合进行划分的过程。相似性是聚类的关键所在。距离、密度等许多种模型都可以用来表示相似性，但是这些模型所刻画的并不是相似性的本质。为了准确反映聚类的本质，本文首先对相似性的形式进行讨论；而后，从相似性的直观含义出发，只保留那些直观的、明确的性质，得出相似性的一般形式。　　
二、构造了聚类假设空间。假设空间是机器学习的重要理论基础。为了进一步研究聚类相关问题，本文在聚类相似性一般形式基础上构造了聚类的假设空间。并利用聚类假设空间分析了导致聚类及聚类评价方法失效的主要原因。　　
三、建立了聚类的模态逻辑表示方法。聚类相似性、数据集合、簇指派是聚类的三个重要组成部分。本文在聚类相似性一般形式基础上将聚类对应于模态逻辑中Kripke语义结构，进而通过模态逻辑公式对聚类的结果进行表示，为深入分析聚类及聚类相关问题建立了理论基础。　　
四、提出了具有通用性的聚类评价方法。本文采用模态逻辑公式对聚类结果中各个数据所反映的信息进行刻画，根据数据反映的不同信息，本文提出了一种基于模态逻辑的聚类代表点方法。在此基础上，本文进一步提出了基于代表点的聚类评价方法，这种方法不受相似性计算方式的限制，具有较强的通用性。此外，代表点本身除了可对聚类结果进行定量的评判外，还可以给出有关聚类结果的定性分析。　　
五、提出了增量聚类的风险理论及评价方法。聚类可以看作是归纳学习，而归纳有“失真”的风险。本文提出并证明了归纳应该就近进行的原则。聚类的模态逻辑表示方法不但有助于常规聚类的评价，还可以用于对增量聚类的评价。不同增量数据中所反映的不同逻辑公式表明了该数据进行归纳的风险。通过降低归纳风险可以得到较好的增量聚类结果；而通过计算增量聚类结果的归纳风险则可以对增量聚类的结果进行评价。　　
最后，本文通过具体的应用实例进一步验证了以上内容。其中合并小样本的应用中分别验证了聚类假设空间的可行性及在此基础上构造的聚类评价方法；航班延误定级应用验证了基于代表点的聚类评价方法的可行性，以及该方法所得到的定量与定性双重结果在实际工程中的应用价值。
学科专业：
授予学位：
学位授予单位：
导师姓名：
学位年度：
在线出版日期：
本文读者也读过
相关检索词
万方数据知识服务平台--国家科技支撑计划资助项目（编号：2006BAH03B01）(C)北京万方数据股份有限公司
万方数据电子出版社后使用快捷导航没有帐号？
只需一步, 快速开始
只需一步，快速开始
查看: 8989|回复: 100
【求助】HYSPLIT在线版数据格式与TrajStat的转换
在线时间 小时
登录后查看更多精彩内容~
才可以下载或查看，没有帐号？
各位前辈好！我想实现的功能是利用HYSPLIT模拟一个点过去50年之内的若干天（可能每隔五年模拟三个月，每个月模拟四天）的气团运动轨迹。由于PC版需要单独下载数据，而我想实现的这种功能下载数据量可能过大，于是选择了在线版。
然而，在线版可以生成的文件类型如下图，而我想要实现的功能是把模拟生成的若干条轨迹导入TrajStat进行聚类分析，画出一条具有统计学意义的轨迹。
问题1.jpg (96.1 KB, 下载次数: 0)
10:59 上传
几种格式我都尝试了一下，只有转换成GIS数据包，其中的SHP文件可以在meteoinfo中打开，如下图二。（还没有安装TrajStat，但是考虑到两款软件都是王老师团队开发的，所以猜测HYSPLIT数据兼容性相似）但是如果是SHP文件的话，似乎就没有办法进一步统计分析了。而在线版生成的彩色图样（图三）却无法在meteoinfo中显示。
问题2.jpg (32.01 KB, 下载次数: 1)
10:59 上传
问题四.jpg (40.09 KB, 下载次数: 0)
11:03 上传
故不知想要实现TraiStat的统计功能，是否一定要用PC版？如何让轨迹显示出彩色并可以分别编辑呢？
另外我尝试安装了一下HYSPLIT的PC版，按照官网提示下载了如下相关软件，其中第三款软件无法正常安装
最后HYSPLIT安装目录下显示的是这样子：
问题五.jpg (19.57 KB, 下载次数: 0)
11:03 上传
综合的结果就是HYSPLIT PC版安装没有成功。。。{:soso_e154:}
综上总共两个问题，还请前辈们不吝赐教！万分感激！！~~~
已有1人关注本帖
密码修改失败请联系qq：
在线时间 小时
额。。。图的位置好像有点错了。。。最后一张图是官网提示的安装HYSPLIT应该安装的系列软件，其中第三款gs905安装失败了。。。
密码修改失败请联系qq：
在线时间 小时
HYSPLIT与TrajStat是做什么的软件呢？以前没有见过呀？
密码修改失败请联系qq：
在线时间 小时
kongfeng0824 发表于
HYSPLIT与TrajStat是做什么的软件呢？以前没有见过呀？
第一个可以进行气团轨迹模拟，第二个可以进行气团轨迹分析~~
密码修改失败请联系qq：
在线时间 小时
lvcha1874 发表于
第一个可以进行气团轨迹模拟，第二个可以进行气团轨迹分析~~
只做轨迹的话不用安装PC的HYSPLIT，TrajStat里有计算轨迹的模块（来自HYSPLIT的PC版）。
MeteoInfo（TrajStat也是一样）对轨迹的显示设置是很丰富的：
Image00365.png (150.34 KB, 下载次数: 0)
14:00 上传
密码修改失败请联系qq：
在线时间 小时
网上计算的轨迹结果数据也可以在MeteoInfo中打开，注意是Trajectory endpoints file。在MeteoInfo的Meteo Data对话框打开数据按钮下HYSPLIT Data -& Trajectory菜单。
Image00366.png (279.84 KB, 下载次数: 1)
14:19 上传
密码修改失败请联系qq：
在线时间 小时
MeteoInfo 发表于
网上计算的轨迹结果数据也可以在MeteoInfo中打开，注意是Trajectory endpoints file。在MeteoInfo的Meteo D ...
果然可以！！！谢谢版主老师！！！
密码修改失败请联系qq：
在线时间 小时
MeteoInfo 发表于
网上计算的轨迹结果数据也可以在MeteoInfo中打开，注意是Trajectory endpoints file。在MeteoInfo的Meteo D ...
版主老师您好，下午我尝试用TrajStat做分析的时候遇到了如下几个问题：
1, 之前我用HYSPLIT对一个站点，一次大气运动的气团轨迹模拟了三个高度，而我想在这三个高度上分别针对时间的变化拟合出三条具有统计学意义的轨迹。不知要是这样的话是否需要三个高度层分别模拟，各自生成一条轨迹呢？
2，我理解的想要生成一条统计学意义的轨迹，应该使用TrajStat里面的Cluster Analysis功能。我注意到user guide里面说要先把traj轨迹转换成.lcd格式。请问在转换时，是不是应该对所有参与cluster analysis计算的轨迹要打包转换，而不能单独转换呢？参与Cluster Analysis的是否也应该是包含所有参与轨迹的打包好的.lcd文件？
3，我截取了HYSPLIT模拟轨迹生成的一份.txt文件如下图，这个文件包括了一个站点固定时间内三个高度层的模拟数据。
我尝试用如下两种方式进行.lcd文件的转换，一种是一起打开六天的轨迹，生成一个traj层；另一种是只打开一天的轨迹，生成一个traj层，但是它们都遇到了如下问题：
这个错误不知是否说明应该给层添加CSV日期信息，但是我还是编辑了一个CSV日期文件，试图添加到一次打开六天的文件生成的层内，但是错误显示跟上面一致。
date.jpg (7.88 KB, 下载次数: 0)
18:44 上传
请问我如果想要实现1当中所列的功能，哪种才是比较好的解决办法呢？
真诚的感谢版主！！！
密码修改失败请联系qq：
在线时间 小时
MeteoInfo 发表于
网上计算的轨迹结果数据也可以在MeteoInfo中打开，注意是Trajectory endpoints file。在MeteoInfo的Meteo D ...
似乎图片的显示出现问题了。。补充下问题3的几张图{:soso_e136:}：
数据1.jpg (30 KB, 下载次数: 0)
18:52 上传
站点三个高度层HYSPLIT模拟txt图
问题1.jpg (13.96 KB, 下载次数: 0)
18:52 上传
两种方式创建层
问题2.jpg (59.83 KB, 下载次数: 0)
18:52 上传
密码修改失败请联系qq：
在线时间 小时
lvcha1874 发表于
似乎图片的显示出现问题了。。补充下问题3的几张图：
站点三个高度层HYSPLIT模拟txt图
从HYSPLIT计算的轨迹数据需要在TrajStat软件中先生成TGS文件（用'Trajectory -& 'Convert To TGS File'菜单），然后再由TGS文件生成shape文件（用'Trajectory -& Convert To Shape File'菜单）。
密码修改失败请联系qq：后使用快捷导航没有帐号？
只需一步, 快速开始
只需一步，快速开始
查看: 2139|回复: 18
求助，如何利用NECP FNL资料、MetoInfo里的插件TrajStat，绘制气流轨迹，聚类分析？
在线时间 小时
登录后查看更多精彩内容~
才可以下载或查看，没有帐号？
已经下载了NECP FNL资料，grib2格式的，安装好了MeteoInfo Java 版和TrajStat插件，已经看了TrajStat的帮助文档，没找着该怎么生成第一步的气流轨迹，困在第一步，恳请知道的老师、同学帮忙解答一下，不胜感激啊~刚接触，确实啥也不懂，看了半天也没明白{:soso_e109:}
密码修改失败请联系qq：
在线时间 小时
HYSPLIT模式需要专门格式的气象数据（ARL格式），grib2格式的气象数据需要先转换。建议你直接下载ARL格式的气象数据来做轨迹，详见此贴：
如果一定要用其它格式气象数据，也可以通过写MeteoInfo脚本程序进行格式转换。
密码修改失败请联系qq：
在线时间 小时
老师，格式转换麻烦吗？刚才试了一下下载ARL格式的，想分析一年的，下载速度好像太慢了
密码修改失败请联系qq：
在线时间 小时
HYSPLIT模式需要专门格式的气象数据（ARL格式），grib2格式的气象数据需要先转换。建议你直接下载ARL格式的 ...
老师，转换格式麻烦吗? 刚才试了一下ARL格式的，我想分析全年的，下载速度好像太慢了点{:soso_e105:}，
密码修改失败请联系qq：
在线时间 小时
老师，转换格式麻烦吗? 刚才试了一下ARL格式的，我想分析全年的，下载速度好像太慢了点，
之前写过的一个脚本程序，根据你自己的需要修改：
#--------------------------------------------------------& && &&&
# Author: Yaqiang Wang& && && && && && && && && && && && && && &
# Date: & && && && && && && && && && && && && && && &&&
# Purpose: Convert GRIB data to ARL data&&
# Note: Sample& && && && && && && && && && && && && && && && && &
#-----------------------------------------------------------
from org.meteoinfo.data.meteodata import MeteoDataInfo
from org.meteoinfo.data.meteodata.arl import ARLDataInfo
from org.meteoinfo.data.meteodata.arl import DataLabel
import os
#---- Set directories
dataDir = &D:/Temp&
#---- Set output data file
outFile = os.path.join(dataDir, 'arl/test2.arl')
#---- Read a GRIB data file
mydata = MeteoDataInfo()
infile = os.path.join(dataDir, 'grib/.pgbh06.gdas.00105.grb2')
print infile
mydata.openNetCDFData(infile)
print 'GRIB file has been opened...'
#---- Set output ARL data info
arlDI = ARLDataInfo()
#---- Set variable and level list
gvars = ['Pressure_surface','Temperature_height_above_ground',\
&&'u-component_of_wind_height_above_ground','v-component_of_wind_height_above_ground',\
&&'Geopotential_height_isobaric','Temperature_isobaric',\
&&'u-component_of_wind_isobaric','v-component_of_wind_isobaric','Vertical_velocity_pressure_isobaric',\
&&'Relative_humidity_isobaric']
avars = ['PRSS','T02M','U10M','V10M','HGTS','TEMP','UWND','VWND','WWND','RELH']
levels = [0,,925,900,875,850,825,800,775,750,700,\
&&650,600,550,500,450,400,350,300,250,225,200,175,150,\
&&125,100,70,50,30,20,10,7,5,3,2,1]
for l in levels:
&&arlDI.levels.add(l)
&&if l == 0:
& & arlDI.LevelVarList.add(['PRSS','T02M','U10M','V10M'])
&&else:
& & arlDI.LevelVarList.add(['HGTS','TEMP','UWND','VWND','WWND','RELH'])
#---- Write ARL data file
dataInfo = mydata.getDataInfo()
arlDI.createDataFile(outFile)
arlDI.X = dataInfo.getXDimension().getValues()
arlDI.Y = dataInfo.getYDimension().getValues()
variables = dataInfo.getVariables()
tNum = dataInfo.getTimeNum()
for t in range(0, tNum):
&&mydata.setTimeIndex(t)
&&atime = dataInfo.getTimes().get(t)
&&aDH = arlDI.getDataHead(mydata.getProjectionInfo(), 'FNL1', 2)
&&arlDI.writeIndexRecord(atime, aDH)
&&lidx = 0
&&for l in arlDI.levels:
& & print l
& & for v in arlDI.LevelVarList[lidx]:
& && &vName = gvars[avars.index(v)]
& && &print vName
& && &if lidx == 0:
& && &&&mydata.setLevelIndex(lidx)
& && &else:
& && &&&variable = dataInfo.getVariable(vName)
& && &&&nlidx = variable.getZDimension().getDimValue().indexOf(l*100.0)
& && &&&mydata.setLevelIndex(nlidx)
& && &gData = mydata.getGridData(vName)
& && &if v == 'PRSS' or v == 'WWND':
& && &&&gData = gData.divide(100)& && &
& && &aDL = DataLabel(atime)
& && &aDL.setLevel(lidx)
& && &aDL.setVarName(v)
& && &aDL.setGrid(99)
& && &aDL.setForecast(0)
& && &arlDI.writeGridData(aDL, gData)
& & lidx += 1
arlDI.closeDataFile()
print 'Finished!'
密码修改失败请联系qq：
在线时间 小时
之前写过的一个脚本程序，根据你自己的需要修改：
好的，谢谢老师！
密码修改失败请联系qq：
在线时间 小时
之前写过的一个脚本程序，根据你自己的需要修改：
老师，看不明白，
密码修改失败请联系qq：
在线时间 小时
老师，看不明白，
如果你一点Python的基础也没有，看不明白很正常。
密码修改失败请联系qq：
在线时间 小时
如果你一点Python的基础也没有，看不明白很正常。
真没接触过Python，先下载数据看看了
密码修改失败请联系qq：
在线时间 小时
好高端，帮顶
密码修改失败请联系qq：