第４１卷第５期２０１８年５月
测绘与空间地理信息
ＧＥＯＭＡＴＩＣＳ＆ＳＰＡＴＩＡＬＩＮＦＯＲＭＡＴＩＯＮＴＥＣＨＮＯＬＯＧＹ
Ｖｏｌ.４１ꎬＮｏ.５Ｍａｙꎬ２０１８
收稿日期:２０１７－１０－１１
作者简介：张㊀胜(１９９２－)ꎬ男ꎬ河南信阳人ꎬ助理工程师ꎬ硕士ꎬ２０１７年毕业于中国地质大学测绘工程专业ꎬ主要从事工程测量㊁
ＧＮＳＳ应用及数据处理工作ꎮ
ＴｒｉｍｂｌｅＲＴＸ后处理解算精度分析
张㊀胜１ꎬ杜士澍１ꎬ屈津年１ꎬ常㊀远２ꎬ翟警卫１ꎬ李常春１
(１.武警黄金第四支队ꎬ辽宁辽阳１１１０００ꎻ２.吉林省交通规划设计院ꎬ吉林长春１３００２１)
摘
要：在实际项目中实验了ＴｒｉｍｂｌｅＲＴＸ后处理在中国内陆解算精度ꎬ并且探索了基于ＲＴＸ－ＰＰ获取ＣＧＣＳ２０００坐标的可行性ꎮ实验结果显示:１)延长数据观测时长(实验数据１２.１５ｈ)ꎬＲＴＸ－ＰＰ的定位可优于１ｃｍ精度ꎻ２)ＲＴＸ－ＰＰ直接提供的ＣＧＣＳ２０００解算结果精度在ｄｍ级别ꎬ而利用更高精度的速度场模型解算的ＣＧＣＳ２０００坐标精度可以在ｃｍ级ꎮ
关键词:ＴｒｉｍｂｌｅＲＴＸꎻ定位精度ꎻＣＧＣＳ２０００ꎻ坐标转换
中图分类号:Ｐ２２６＋
.３㊀㊀㊀文献标识码:Ａ㊀㊀㊀文章编号:１６７２－５８６７(２０１８)０５－０１７０－０３
ＡｃｃｕｒａｃｙＡｎａｌｙｓｉｓｏｆｔｈｅＴｒｉｍｂｌｅＣｅｎｔｅｒＰｏｉｎｔＲＴＸＰｏｓｔ－ｐｒｏｃｅｓｓｉｎｇ
ＺＨＡＮＧＳｈｅｎｇ１ꎬＤＵＳｈｉｓｈｕ１ꎬＱＵＪｉｎｎｉａｎ１ꎬＣＨＡＮＧＹｕａｎ２ꎬＺＨＡＩＪｉｎｇｗｅｉ１ꎬＬＩＣｈａｎｇｃｈｕｎ１
(１.ＦｏｕｒｔｈＤｅｔａｃｈｍｅｎｔｏｆＧｏｌｄＨｅａｄｑｕａｒｔｅｒｓꎬＣｈｉｎｅｓｅＡｒｍｅｄＰｏｌｉｃｅＦｏｒｃｅꎬＬｉａｏｙａｎｇ１１１０００ꎬＣｈｉｎａꎻ
２.ＪｉｌｉｎＰｒｏｖｉｎｃｉａｌＣｏｍｍｕｎｉｃａｔｉｏｎＰｌａｎｎｉｎｇａｎｄＤｅｓｉｇｎＩｎｓｔｉｔｕｔｅꎬＣｈａｎｇｃｈｕｎ１３００２１ꎬＣｈｉｎａ)
Ａｂｓｔｒａｃｔ:ＩｎｐｒａｃｔｉｃａｌｐｒｏｊｅｃｔｓꎬｔｈｅａｃｃｕｒａｃｙｏｆＴｒｉｍｂｌｅＲＴＸｐｏｓｔ－ｐｒｏｃｅｓｓｉｎｇｉｎｉｎｌａｎｄＣｈｉｎａｉｓｉｎｖｅｓｔｉｇａｔｅｄꎬａｎｄｔｈｅｆｅａｓｉｂｉｌｉｔｙｏｆｏｂ￣ｔａｉｎｉｎｇＣＧＣＳ２０００ｃｏｏｒｄｉｎａｔｅｓｂａｓｅｄｏｎＲＴＸ－ＰＰｉｓｅｘｐｌｏｒｅｄ.Ｔｈｅｅｘｐｅｒｉｍｅｎｔａｌｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔ:１)ｔｏｅｘｔｅｎｄｔｈｅｔｉｍｅｏｆｏｂｓｅｒｖａ￣ｔｉｏｎｄａｔａ(ｅｘｐｅｒｉｍｅｎｔａｌｏｂｓｅｒｖａｔｉｏｎｄａｔａｌｅｎｇｔｈ１２.１５ｈ)ꎬｔｈｅｌｏｃａｌｉｚａｔｉｏｎｏｆＲＴＸ－ＰＰｃａｎｒｅａｃｈｂｅｔｔｅｒｔｈａｎ１ｃｍａｃｃｕｒａｃｙꎻ２)ＣＧＣＳ２０００ＲＴＸ－ＰＰｐｒｏｖｉｄｅｓｔｈｅａｃｃｕｒａｃｙｏｆｓｅｔｔｌｅｍｅｎｔｒｅｓｕｌｔｓｄｉｒｅｃｔｌｙａｔｔｈｅＤＭｌｅｖｅｌꎬａｎｄｔｈｅｕｓｅｏｆＣＧＣＳ２０００ｃｏｏｒｄｉｎａｔｅｐｒｅｃｉ￣
ｓｉｏｎｆｉｅｌｄｍｏｄｅｌｔｏｓｏｌｖｅｔｈｅｈｉｇｈｅｒｐｒｅｃｉｓｉｏｎｏｆｔｈｅｓｐｅｅｄｃａｎｂｅａｔｔｈｅＣＭｌｅｖｅｌ.Ｋｅｙｗｏｒｄｓ:ＴｒｉｍｂｌｅＲＴＸꎻｐｏｓｉｔｉｏｎｉｎｇａｃｃｕｒａｃｙꎻＣＧＣＳ２０００ꎻｃｏｏｒｄｉｎａｔｅｃｏｎｖｅｒｓｉｏｎ
０㊀引㊀言
早在２０１１年ꎬＴｒｉｍｂｌｅ公司推出了一种新型的实时查分扩展技术  ＣｅｎｔｅｒＰｏｉｎｔＲＴＸ(Ｒｅａｌ－ＴｉｍｅｅＸｔｅｎｄ￣ｅｄ)[１
－２]
ꎬ实时的ＧＰＳ/ＧＬＯＮＡＳＳ卫星精确轨道和钟差数
据ꎬ这种实时定位服务不依赖基准站就能提供厘米级定位精度ꎮ相比其他在线提供ＰＰＰ(精密单点定位)[１ꎬ３]服务的系统ꎬＴｒｉｍｂｌｅＲＴＸ支持日本的ＱＺＳＳ信号ꎬ以此提升了载波相位模糊度求解及其收敛能力ꎮ
随后ꎬＴｒｉｍｂｌｅ提供了新的后处理服务ＲＴＸ－ＰＰ[１]ꎮ
这个应用服务允许用户上传ＧＮＳＳ观测数据到Ｔｒｉｍｂｌｅ公司ＲＴＸ后处理服务中心ꎬ用户可以选择不同的坐标系和
构造板块来获取最终的定位成果ꎮ观测数据达到２４ｈ后ꎬＲＴＸ－ＰＰ解算结构水平和垂直方向精度能优于１ｃｍꎮ因为解算收敛时间远小于１ｈ时ꎬ因此仅仅利用１ｈ观测的数据水平方向精度通常优于２ｃｍ
[４]
ꎮ
Ｔｒｉｍｂｌｅ的ＲＴＸ技术无疑大大地简化了常规ＧＰＳ静
态作业模式ꎬ方便且快速地获得高精度定位结果ꎬ特别是在没有足够的测绘基础设施的地区ꎮ然而这种模式在中国大陆地区的可靠性少有文章提及ꎮ本文利用实测数据与国家控制点坐标进行对比分析ＲＴＸ解算精度ꎻ另一方面ꎬ现在我国测绘成果要求统一采用ＣＧＣＳ２０００坐标系ꎬ
虽然ＲＴＸ－ＰＰ服务提供ＣＧＣＳ２０００坐标ꎬ但是实际工程项目中是否可靠也是本文研究的目的ꎮ
本文利用实际项目中的观测数据验证ＲＴＸ－ＰＰ结果的可靠性ꎬＡ００８位于宁波市ꎬ是国家Ｂ级ＧＰＳ控制点ꎮ因为国家ＧＰＳ控制点只有ＣＧＣＳ２０００坐标ꎬ因此ꎬ首先利用ＧＡＭＩＴ/ＧＬＯＢＫ软件求解３个控制点的ＩＴＲＦ２００８框架下坐标ꎬ验证ＲＴＸ－ＰＰ绝对定位的准确性ꎻ其次利用ＲＴＸ－ＰＰ直接求取ＣＧＣＳ２０００坐标ꎬ与国家控制点真实的
ＣＧＣＳ２０００坐标对比ꎻ观测数据采用ＴｒｉｍｂｌｅＲ１０仪器ꎬ数据质量检测结果见表１ꎮ
表１㊀Ａ００８控制点观测数据质量
Ｔａｂ.１㊀Ａ００８ｃｏｎｔｒｏｌｐｏｉｎｔｏｂｓｅｒｖａｔｉｏｎｄａｔａｑｕａｌｉｔｙ
文件名观测时长ｈｏｕｒ(ｓ)完好性观测值比例ｍｐ１(Ｌ１载波多路径效应)ｍｐ２(Ｌ２载波多路径效应)ｏ/ｓｌｐｓ(观测值和周跳比)Ａ００８１４５０.１ｏ１２.１５９１％０.３５０.３９２０３３３
１㊀ＲＴＸ－ＰＰ定位结果精度
ＲＴＸ－ＰＰ后处理利用由Ｒｅａｌ－ｔｉｍｅＣｅｎｔｅｒＰｏｉｎｔＲＴＸ系统导出的精确轨道和时钟ꎬ并通过Ｌ１和Ｌ２频率的 电离层自由 组合消除了ＲＴＸ定位中的一阶电离层效应ꎬ同时对流层通过一个模型来处理ꎬ这个模型中附加了垂直湿延迟和在南北和东西方向的两个梯度ꎮ处理中估计ＧＰＳ和ＧＬＯＮＡＳＳ的独立接收机时钟误差ꎬ而ＱＺＳＳ卫星被认为具有与ＧＰＳ相同的接收机时钟误差ꎮ此外ꎬＲＴＸ－ＰＰ还做了如下解算更正[４]:１)卫星和接收机天线校正ꎻ２)固体潮汐ꎻ３)波浪潮ꎻ４)海潮加载效应ꎻ５)相对论修正ꎻ６)相位叠加ꎻ７)代码和载波相位偏差ꎮ
用户将ＧＮＳＳ观测数据上传到ＲＴＸ－ＰＰ在线服务中心(ｈｔｔｐ://ｗｗｗ.ｔｒｉｍｂｌｅｒｔｘ.ｃｏｍ/Ｈｏｍｅ.ａｓｐｘ)ꎬ通过选择不同的坐标系和构造板块来进行转换ꎬ最后通过发送解算结果至用户ꎬ整个过程只需要几分钟ꎮ
需要注意的是ꎬ上从结果上传的数据文件必须包含双频伪和载波相位观测(Ｌ１和Ｌ２)ꎬ数据有效历元
需保证在６０ｍｉｎ ２４ｈ以内ꎻＲＴＸ－ＰＰ支持的天线类型(ｈｔｔｐ://ｗｗｗ.ｔｒｉｍｂｌｅｒｔｘ.ｃｏｍ/ＳｕｐｐｏｒｔｅｄＤｅｖｉｃｅｓ.ａｓｐｘ)ꎬ并且天线高参考点为天线底座底部ＢＡＭ(Ｂｏｔｔｏｍｏｆａｎｔｅｎｎａ㊀㊀ｍｏｕｎｔ)ꎮ对于ＴｒｉｍｂｌｅＲ１０接收机ꎬ天线ＴＲＭＲ１０的天线底座底部ＢＡＭ与天线相位参考点ＡＲＰ为同一位置[５]ꎬ如图１所示
ꎮ
图１㊀ＴｒｉｍｂｌｅＲ１０接收机天线相位参考点示意图Ｆｉｇ.１㊀ＴｒｉｍｂｌｅＲ１０ｒｅｃｅｉｖｅｒａｎｔｅｎｎａｐｈａｓｅ
㊀㊀㊀㊀ｒｅｆｅｒｅｎｃｅｐｏｉｎｔａｗａｒｅｎｅｓｓｍａｐ
㊀㊀表２所示为ＲＴＸ－ＰＰ解算的坐标框架ＩＴＲＦ２００８㊁历元２０１７.３９３结果与搜集到的控制点坐标对比ꎮ从结果可以看出三分量较差分别为０.０１１４ｍ㊁０.０１０６ｍ㊁０.０００４ｍꎮ综合实际测试结果分析ꎬ延长数据观测时长ꎬＲＴＸ－ＰＰ的定位可优于１ｃｍ精度ꎮ
表２㊀ＲＴＸ－ＰＰ解算结果与控制点坐标对比
Ｔａｂ.２㊀ＲＴＸ－ＰＰｓｏｌｕｔｉｏｎｒｅｓｕｌｔｓａｎｄｃｏｎｔｒｏｌｐｏｉｎｔｃｏｏｒｄｉｎａｔｅｓｃｏｍｐａｒｉｓｏｎ
解类型Ｘ(ｍ)Ｙ(ｍ)Ｚ(ｍ)Ｍｘ(ｍｍ)Ｍｙ(ｍｍ)Ｍｚ(ｍｍ)控制点坐标－２９０９８８６.５５９６４７０４３８５.２９７４３１６４６２９.９２９６４.９６.５４.５ＲＴＸ－ＰＰ－２９０９８８６.５７１４７０４３８５.３０８３１６４６２９.９３０６８５
２㊀利用ＲＴＸ－ＰＰ获取ＣＧＣＳ２０００坐标
ＲＴＸ－ＰＰ在线服务已经提供ＣＧＣＳ２０００坐标成果ꎬ但是这个成果是基于整个欧亚板块欧拉矢量[７]求得ꎬ精度必然不会很高ꎬ本文提出基于ＲＴＸ－ＰＰ获取的ＩＴＲＦ２００８坐标进行框架转换ꎬ实现从ＩＴＲＦ２００８坐标转换到ＣＧＣＳ２０００坐标ꎬ如图２所示ꎬ具体步骤如下[８－１２]:１)框架变换:２０１７.３９３历元下ꎬＩＴＲＦ２００８框架转换到ＩＴＲＦ９７框架ꎻ
２)历元变换:ＩＴＲＦ９７框架中ꎬ２０１７.３９３历元坐标转换到２０００.０历元下的坐标ꎮ
从ＩＴＲＦ(ｈｔｔｐ://ｉｔｒｆ.ｅｎｓｇ.ｉｇｎ.ｆｒ/ｔｒａｎｓ＿ｐａｒａ.ｐｈｐ)可以获得ＩＴＲＦ２００８框架与其他ＩＴＲＦ框架的转换参数ꎬ见表
３ꎮ图２㊀ＩＴＲＦ２００８与ＣＧＣＳ２０００坐标转换流程图
Ｆｉｇ.２㊀ＩＴＲＦ２００８ａｎｄＣＧＣＳ２０００ｃｏｏｒｄｉｎａｔｅ
㊀㊀㊀㊀ｃｏｎｖｅｒｓｉｏｎｆｌｏｗｃｈａｒｔ１７１
第５期张㊀胜等:ＴｒｉｍｂｌｅＲＴＸ后处理解算精度分析
表３㊀ＩＴＲＦ２００８转换到ＩＴＲＦ９７的参数Ｔａｂ.３㊀ＩＴＲＦ２００８ｐａｒａｍｅｔｅｒｓｃｏｎｖｅｒｔｅｄｔｏＩＴＲＦ９７
参数单位Ｔｘ
ｍｍ
Ｔｙ
ｍｍ
Ｔｚ
ｍｍ
Ｄ
ｐｐｂ
Ｒｘ
ｍａｓ
Ｒｘ
ｍａｓ
Ｒｘ
ｍａｓ
ＥＰＯＣＨ
变化率单位Ｔ
ｘ
ｍｍ/ｙＴ
ｙ
ｍｍ/ｙＴ
ｚ
ｍｍ/ｙＤ
ｐｐｂ/ｙＲ
ｘ
ｍａｓ/ｙＲ
ｘ
ｍａｓ/ｙＲ
ｘ
ｍａｓ/ｙＩＴＲＦ９７４.８２.６
－３３.２
２.９２０.０００.０００.０６２０００.０
Ｒａｔｅｓ
０.１
－０.５－３.２
０.０９
０.００
０.００
０.０２
㊀㊀注:１ｐｐｂ＝１０－９ꎻ１ｍａｓ＝０.００１ᵡ＝４.８４８１３ˑ１０－９ｒａｄ
ＸＳＹＳＺＳéëêêêùûúúú＝ＸＹＺéëêêêùûúúú＋ＴｘＴｙＴｚéëêêêùûúúú＋Ｄ－ＲｚＲｙＲｚＤ－Ｒｘ－Ｒｙ
Ｒｘ
ＤéëêêêùûúúúＸＹＺéëêêêù
û
úúú(１)
式(１)中ꎬＸ㊁Ｙ㊁Ｚ是ＩＴＲＦ２００８框架下的坐标值ꎬＸＳ㊁ＸＹ㊁ＸＺ是ＩＴＲＦ９７框架下的坐标值ꎮ对于不同历元ꎬ还涉及历元转换ꎬ公式如下:
Ｐｔ()＝ＰＥＰＯＣＨ()＋Ｐ
ˑ(ｔ－ＥＰＯＣＨ)(２)
式(２)中ꎬＰ表示历元(ＥＰＯＣＨ)２０００.０时刻转换参
数ꎬ也就是上表中的数值ꎮＰ
表示Ｐ的速率ꎮ
历元转换涉及站点速度场ꎬ本文采用魏子卿提出的
格网平均值法[８]ꎬ参考中国大陆Ⅱ级活动块体划分[１３
－４]
ꎬ
常远
Ａ００８点位于鲁东－黄海块体与华南块体交接地带ꎬ其地壳运动背景场整体是在往东偏南运动ꎬ年速率在ｃｍ级别
(Ｖｘ＝－０.０２９９ꎬＶｙ＝－０.０１０１ꎬＶｚ＝－０.０１４８ꎬ单位:ｍ/ｙｒ)ꎮ
Ｘｔ＝２０００.０＝Ｘｔ＋(２０００.０－ｔ)Ｖ
(３)
式(３)称为历元转换ꎬ式中ꎬＸｔ＝２０００.０对应ＩＴＲＦ９７ꎬ
２０００.０历元下坐标向量ꎬ也就是ＣＧＣＳ２０００坐标ꎻＸｔ是ＩＴＲＦ９７框架下ｔ历元(本文ｔ＝２０１７.３９３)的坐标ꎻＶ表示站点速度向量ꎮ
综合表４结果ꎬＲＴＸ－ＰＰ直接提供的ＣＧＣＳ２０００解算结果精度在ｄｍ级别ꎬ而利用更高精度的速度场模型解算的ＣＧＣＳ２０００坐标精度可以在ｃｍ级ꎮ
表４㊀ＲＴＸ－ＰＰ获取ＣＧＣＳ２０００坐标结果对比
Ｔａｂ.４㊀ＲＴＸ－ＰＰｔｏｏｂｔａｉｎＣＧＣＳ２０００ｃｏｏｒｄｉｎａｔｅｒｅｓｕｌｔｓｃｏｎｔｒａｓｔ
类型Ｘ(ｍ)较差Ｙ(ｍ)较差Ｚ(ｍ)较差Ｂ级控制点－２９０９８８６.０４５００
４７０４３８５.５１０１０
３１６４６３０.１２３６０
ＲＴＸ－ＰＰ－２９０９８８６.２０００.１５５４７０４３８５.４３２－０.０７８１３１６４６３０.０９２－０.０３１６框架转换
－２９０９８８６.０５６
０.０１１４７０４３８５.４９９
－０.０１１１
３１６４６３０.１６３
０.０３９４３㊀结束语
ＲＴＸ－ＰＰ为美国天宝公司提供的后处理服务ꎬ能够快
速处理静态双频多系统观测数据(ＧＰＳ㊁北斗和伽利略)获得高精度定位结果ꎮ本文采用实测ＧＮＳＳ数据验算ＲＴＸ－ＰＰ在中国内陆地区的精度ꎬ同时利用坐标系框架转换获得ＣＧＣＳ２０００坐标ꎬ以ＲＴＸ－ＰＰ解算的坐标框架ＩＴＲＦ２００８㊁历元２０１７.３９３结果与搜集到的控制点坐标对比ꎬ从结果上可以看出三分量较差分别为０.０１１４ｍ㊁
０.０１０６ｍ㊁０.０００４ｍꎮ综合实际测试结果分析ꎬ延长数据观测时长ꎬＲＴＸ－ＰＰ的定位可优于１ｃｍ精度ꎮ另一方面ꎬ基于ＲＴＸ－ＰＰ获取的ＩＴＲＦ２００８坐标进行框架转换ꎬ实现从ＩＴＲＦ０８坐标转换到ＣＧＣＳ２０００坐标结果精度在ｄｍ级别ꎬ采用更高精度的速度场模型解算ꎬ可以实现ｃｍ级ＣＧＣＳ２０００坐标ꎮ
参考文献:
[１]㊀ＣｈｅｎＸꎬＡｌｌｉｓｏｎＴꎬＣａｏＷꎬｅｔａｌ.ＴｒｉｍｂｌｅＲＴＸꎬａｎｉｎｎｏ￣
ｖａｔｉｖｅｎｅｗａｐｐｒｏａｃｈｆｏｒｎｅｔｗｏｒｋＲＴＫ[Ｊ].ＰｒｏｃｅｅｄｉｎｇｓｏｆＩＯＮＧＮＳＳꎬ２０１１(３):２２１４－２２１９.
[２]㊀徐以厅ꎬ宋济宇.ＴｒｉｍｂｌｅＲＴＸ技术综述[Ｊ].测绘通报ꎬ
２０１４(２):１３７－１３８.[３]㊀徐永斌.在线ＰＰＰ软件的精度比较分析[Ｊ].测绘与空
间地理信息ꎬ２０１１ꎬ３４(４):１６７－１６９.[４]㊀ＴｒｉｍｂｌｅＮａｖｉｇａｔｉｏｎＬｉｍｉｔｅｄ[ＥＢ/ＯＬ].(２０１４－０８－１６)
[２０１７－０９－１０].ｈｔｔｐ://ｗｗｗ.ｔｒｉｍｂｌｅｒｔｘ.ｃｏｍ/ＵｐｌｏａｄＦｏｒｍ.ａｓｐｘ.
[５]㊀ＩＧＳ[ＥＢ/ＯＬ].(１９９４－０１－２３)[２０１７－０９－１０].ｈｔｔｐｓ://ｉｇｓｃｂ.ｊｐｌ.ｎａｓａ.ｇｏｖ/ｉｇｓｃｂ/ｓｔａｔｉｏｎ/ｇｅｎｅｒａｌ/ａｎｔｅｎｎａ.ｇｒａ.[６]㊀ＨｅｒｒｉｎｇＴＡꎬＤｏｎｇＤꎬＫｉｎｇＲＷ.Ｓｕｂ‐ｍｉｌｌｉａｒｃｓｅｃｏｎｄｄｅ￣
ｔｅｒｍｉｎａｔｉｏｎｏｆｐｏｌｅｐｏｓｉｔｉｏｎｕｓｉｎｇＧｌｏｂａｌＰｏｓｉｔｉｏｎ
ｉｎｇＳｙｓｔｅｍｄａｔａ[Ｊ].ＧｅｏｐｈｙｓｉｃａｌＲｅｓｅａｒｃｈＬｅｔｔｅｒｓꎬ２０１３ꎬ１８(１０):１８９３－１８９６.
[７]㊀吴富梅ꎬ刘光明ꎬ魏子卿.利用局域欧拉矢量法建立ＣＧＣＳ２０００速度场模型[Ｊ].武汉大学学报：信息科学版ꎬ２０１２ꎬ３７(４):４３２－４３５.
[８]㊀魏子卿ꎬ刘光明ꎬ吴富梅.２０００中国大地坐标系：中国大陆速度场[Ｊ].测绘学报ꎬ２０１１ꎬ４０(４):４０３－４１０.[９]㊀高乐.ＣＧＣＳ２０００板块运动模型的建立[Ｄ].济南：山东
科技大学ꎬ２０１１.
(下转第１７５页)
２
７１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀
测绘与空间地理信息㊀
㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０１８年
①以分幅数据为作业单元ꎬ根据地表覆盖类型㊁点云
扫描日期的差异等ꎬ适当划分交互处理单元ꎬ处理单元范围线应予以保存ꎻ②理解分类工具的计算原理ꎬ能够熟练㊁灵活㊁综合地利用相关分类工具ꎻ
③测试调整相关参数设置ꎬ确保达到最优分类效果ꎬ
减少后期人工分类工作量ꎬ由于设备㊁地形㊁季节㊁气候等多方面原因ꎬ无法确定固定参数ꎬ只能根据分类后的情况确定具体参数ꎻ
④在利用数字正射影像辅助作业时ꎬ应分析其航摄
日期与点云数据间的差异性ꎮ
５)植被降高程处理高于地面的植被降高程处理ꎬ处理前后效果对比如
图２所示
ꎮ
图２㊀效果对比Ｆｉｇ.２㊀Ｅｆｆｅｃｔｏｆｃｏｎｔｒａｓｔ
下面以ＴｅｒｒａＳｏｌｉｄ为例介绍降植被高程的方法[３]:
①量测植被点高出地面的距离ꎬ圈出要降的植被范
围ꎬ按照植被高程分块处理(本次选取一部分为例ꎬ对比处理效果)ꎮ②将该块内的植被点分回到缺省层ꎬ地面点确保留在地面点层ꎮ
③针对所选区域ꎬ使用 Ｃｌａｓｓｉｆｙｂｙｈｅｉｇｈｔｆｒｏｍ
ｇｒｏｕｎｄ 工具从缺省层中选出距离地面０.３ １.０ｍ的点ꎬ放到相应的植被扩展层中ꎮ
④使用 ｆｉｎｄ 工具ꎬ出植被扩展层中的点进行编
辑ꎬ根据量测的植被高出地面的高度ꎬ对这些点进行将高程处理ꎬ在 Ｅｄｉｔｓｅｖｅｒａｌｐｏｉｎｔｓ 工具的 Ｄｚ 中输入高程的改变值ꎬ即改变所选中点的高程值ꎮ５㊀ＤＥＭ处理
１)数据构ＴＩＮ
①点云地面点构ＴＩＮ时ꎬ点云数据应外扩３００ｍꎬ在
外扩３００ｍ时周边点要素仍不能满足要求时要继续外扩ꎬ直至能满足三角网和真实地形要求ꎮ
②采用不规则三角网(ＴＩＮ)内插ＤＥＭꎬ要保证生成
ＴＩＮ时各种参数正确ꎻ确保ＴＩＮ覆盖整个图幅范围ꎬ并向图廓外适当延伸ꎻ保证三角形构网合理ꎮ如发现异常问题ꎬ需对相关数据检查修改后重新构网ꎮ
２)内插ＤＥＭ格网间距５ｍˑ５ｍꎮ
３)ＤＥＭ裁切生成ＤＥＭ时应按照上文规定的数据范围进行裁切处理ꎮ
４)ＤＥＭ接边要求
相邻数字高程模型应接边ꎮ相邻ＤＥＭ模型重叠带
内同名格网点的高程差异若大于２倍ＤＥＭ高程中误差ꎬ则改正后接边ꎮ接边后同名格网点的高程值应保持一致ꎮ
６㊀结束语
我们使用本文的技术方法ꎬ完成了试验区ＤＥＭ数据的制作ꎬ在数据处理过程中ꎬ对一些特殊问题进行了归纳总结ꎬ形成了指导生产的操作手册和问题处理方法ꎬ这为以后批量进行基础地理信息更新积累了技术方法㊁奠定了技术基础ꎮ在以后的生产过程中ꎬ还需进一步优化生产流程㊁提高生产效率以满足批量生产的要求ꎮ
参考文献:
[１]㊀贾玉明.ＬＩＤＡＲ机载激光雷达数据制作ＤＥＭ原理及应
用[Ｊ].广东输电与变电技术ꎬ２００７(５):２８－３０.
[２]㊀李明ꎬ孙彬淳ꎬ张俊ꎬ等.基于机载ＬｉＤＡＲ数据的鄱阳湖
地区高精度ＤＥＭ生产研究[Ｊ].测绘与空间地理信息ꎬ
２０１６ꎬ３９(３):１６４－１６７ꎬ１７２.
[３]㊀苏春梅ꎬ曹殿才ꎬ段凯ꎬ等.基于机载ＬｉＤＡＲ数据制作高
精度ＤＥＭ产品研究[Ｊ].测绘与空间地理信息ꎬ２０１７ꎬ４０(２):７２－７４ꎬ７８.
[编辑：刘莉鑫]
(上接第１７２页)
[１０]㊀苗龙.ＩＴＲＦ２００８与ＣＧＣＳ２０００坐标系的转换[Ｊ].地理空间信息ꎬ２０１１ꎬ９(６):１４４－１４５.
[１１]㊀蒋志浩ꎬ张鹏ꎬ秘金钟ꎬ等.基于ＣＧＣＳ２０００的中国地壳水平运动速度场模型研究[Ｊ].测绘学报ꎬ２００９ꎬ３８(６):
４７１－４７６.
[１２]㊀吴富梅ꎬ刘光明ꎬ魏子卿.利用局域欧拉矢量法建立
ＣＧＣＳ２０００速度场模型[Ｊ].武汉大学学报：信息科学版ꎬ２０１２ꎬ３７(４):４３２－４３５.
[１３]㊀柴洪洲ꎬ崔岳ꎬ明锋.最小二乘配置方法确定中国大陆
主要块体运动模型[Ｊ].测绘学报ꎬ２００９ꎬ３８(１):６１－６５.
[１４]㊀李强ꎬ游新兆ꎬ杨少敏ꎬ等.中国大陆构造变形高精度大
密度ＧＰＳ监测  现今速度场[Ｊ].中国科学：地球科学ꎬ２０１２ꎬ４２(５):６２９－６３２.
[编辑：任亚茹]
５
７１第５期
孙㊀凯等：应用机载激光雷达数据制作ＤＥＭ的方法探讨