摘要

GPS速度场反映区域的地壳变形特征,而目前主流的GIS平台缺少GPS速度场的制图表达,不利于和现有的矢量和栅格数据进行叠加分析。利用ArcGIS为Python提供的Arcpy/Arcgisscripting模块实现速度矢量图形的输出。具体处理过程需要编写代码对选定站点速度进行欧拉矢量计算,实现框架旋转变换,以突显差异运动特征;还包括指定置信度的误差椭圆计算。输出shape文件包括GPS点文件,其属性表为输入速度数据,以及矢量箭头线文件和误差椭圆线文件。最后以中国大陆及周边GPS连续观测站的计算结果,旋转变换到欧亚框架下,实现GPS速度场在ArcGIS中的制图表达。