基于PixelGrid的全球地理信息资源建设应用
- 格式:pdf
- 大小:1.93 MB
- 文档页数:4
第32卷第1期2018年1月北京测绘Beijing Surveying and MappingVol. 32 No. 1January 2018
引文格式:韩晓霞,王保前,许彪,等.基于PbcelGdd的全球地理信息资源建设应用[J].北京测绘,2018(1):44-47.
DOI:10. 19580/i. cnki 1007-3000. 2018. 01. 011
基于PixelGrid的全球地理信息资源建设应用
韩晓霞1王保前2许彪1孙钰珊1赵栋梁3(1.中国测绘科学研究院,北京100830; 2.北京四维空间数码科技有限公司,北京100083;
3.黑龙江第二测绘工程院,黑龙江哈尔滨150025)
[摘要]针对以国产卫星影像为主要数据源,为适应全球地理信息资源建设对海量遥感数据处理规 模化、产品生产业务化和数据生产高效智能化的迫切需求,PixelGrid系统实现了从大规模影像联合区域网平
差到DSM、DEM、DOM的快速自动化处理,为全球地理信息资源建设提供技术支撑。
[关键词]高分辨率遥感影像数据一体化测图系统(PixelGrid);全球测图;大规模影像联合区域网平
差;云区域自动提取
[中图分类号]P208 P237 [文献标识码]A [文章编号]1007 —3000(2018)01 — 0044 — 4
〇引言
以“全球地理信息资源建设工程”为背景,并
以国产卫星影像为主要数据源,为适应全球地理
信息资源建设对海量遥感数据处理规模化、产品
生产业务化和数据生产高效智能化的迫切需求,
以解决海量卫星影像数据无控、精确、快速、一体
化智能处理为目标,基于超多核计算机(云服务
器)、高性能GPU/CPU集群计算机系统、高速存
储传输网络的高性能集群处理平台,PixelGrd系
统实现了从大规模影像联合区域网平差到 DEM、DOM的快速自动化处理,为全球地理信息
资源建设提供技术支撑,并已在多个测绘地理信
局的相关业务中得到了生产应用。
1关键技术及流程
1.1关键技术
1.1.1海量连接点及参考控制点自动提取
PxelGrd系统提出基于投票的匹配策略,以
公共的地理信息数据(如SRTM和谷歌影像)作
为控制参考,自动获取卫星影像间的海量连接
点,同时与谷歌影像自动匹配以获取大量的参考
控制点。投票法匹配[1]是一种整体匹配策略,用卫星
影像和底图影像上分别提取的边缘特征可能存
在的对应关系在参数空间投票直接进行影像间
快速匹配,其思想来源于数字图像处理中的
Hough变换,也可称为广义Hough变换影像匹
配。同时,为了加快匹配速度,PixelGrid系统采
用由粗到精的金字塔影像匹配策略,通常情况
下,19的金字塔影像即可兼顾效率与可靠性。
1.1.2大规模影像联合区域网平差
PixelGrid系统将有理函数模型(Rational
FunctionModel,RFM)[]和基于像方的仿射变换
补偿模型作为联合区域网平差的基础方程,通过结
合连接点、谷歌影像参考控制点、ICESat测高数据
及其他辅助控制信息,利用光束法区域网平差整体
联合求解影像的像方仿射变换定向参数及加密点
物方坐标,实现大规模影像的高精度快速定位。
为提高联合区域网平差结果的可靠性和稳
定性,PxelGrd系统采用专有的粗差探测算法,
用以自动定位并剔除像点观测和控制点观测中
可能存在的大、中、小等多种类型的粗差,为后续
处理提供高可靠性的定向参数。同时,系统利用
验后方差分量估计原理,通过计算区域网平差系
统的精度矩阵和可靠性矩阵,对网形的优化设计
[收稿日期]2017 -12—17[基金项目]测绘地理信息公益性行业科研专项(201511078)。[作者筒介]韩晓霞(1984 —),女,河北沧州人,硕士,工程师,研究方向为地图综合。 E-mail: Ixql957@sohu.
com第32卷第1期韩晓霞,王保前,许彪,等:基于PixelGrid的全球地理信息资源建设应用
及质量控制等提供技术支撑[3]。
1.1.3自动影像匹配生成DSM
PixelGrid系统基于核线几何约束条件下的
多视影像匹配策略,融合多种匹配基元(特征点、
格网点、边缘等),将影像的纹理信息、影像与影
像之间的几何信息(核线约束条件等)、物体与物
体之间的几何约束信息(物体表面为分段光滑的
假设等)运用到匹配算法中,采用由粗到精的金
字塔影像匹配、匹配参数自适应调节、概率松弛
匹配等策略,用以生成数目巨大的匹配特征,实
现基于多基线、多重匹配特征的匹配技术自动获
取 DSM[4]。
该技术的基本流程为:首先,对原始影像进
行预处理并创建金字塔影像;其次,采用基于几
何约束条件的互相关匹配GCCC (Geometrically
ConstrainedCross-Correlation)核心算法分别匹
配各级金字塔影像的匹配基元(包括特征点、格
网点、边缘关键点等),匹配的相似性度量测度基
于多个立体像对的相关系数平均和SNCC(Sum
of NormalizedCross-Correlation),可有效解决匹
配多解和遮蔽等原因引起的误匹配问题;最后,
最小二乘法精细匹配原始影像上的各个匹配基
元:单点匹配采用改进的基于几何约束条件的多
片最小二乘匹配方法、边缘匹配采用物方空间下
的最小二乘 B-Spline Snakes(LSB-Snakes)方法。
1.1.4云区域自动提取
卫星影像中的云区域对影像密集匹配易产
生成片的匹配粗差或者空洞。PixelGrid系统首
先通过DSM与参考SRTM的配准和比较提取
出疑似云区域,再通过影像分析最终识别和确定
云区域,同时统计影像中的云量。该技术主要包
含以下三个步骤:
(1) 基于3D最小二乘表面匹配泛型数学模
型自动匹配DSM与参考SRTM,并校准参考
SRTM[8];
(2) 比较DSM与校准后的参考DEM之间的
差值提取云区域种子点,并反投影到像方得到种
子点的像方坐标;
(3) 根据影像前景分析方法对种子点进行分
类,去除低亮度纹理缺乏区域的影响;然后种子点在影像上做区域增长;最后对提取的区域进行
分类,区分出云和雪两种区域,去除其中的雪区
域,得到最终的云区域。
1.1.5 DEM自动提取
DEM自动提取主要包括DSM自动滤波和
漏洞自动修补两部分,PixelGrid系统提出一种以
SRTM的坡度信息作为辅助,局部渐进优化的方
式自动滤波DSM以获取DEM。由于受云覆盖、
误匹配、纹理缺乏等因素的影响,DSM中可能会
出现粗差或者漏洞问题,系统通过结合SRTM,
自动修补DEM中由于云遮挡及误匹配等其他原
因造成的漏洞区域,兼顾精度的同时保证DEM
的完整性,减少后续人工编辑的工作量,提高生
产效率。
1.1.6 DOM制作及智能镶嵌
PxelGrid系统基于高速局域网的分布式多
核多进程和显卡加速(CPU + GPU)的并行处理
技术,实现卫星影像的准实时正射纠正及融合。
在影像色彩一致性处理方面选用基于Walls滤
波器的整体匀色解决方案,色调统计参数基于全
测区所有影像,以最大程度上保证镶嵌结果
DOM具有整体的色彩一致性。
系统利用最优生成树方法智能检测镶嵌线,
根据影像重叠区域亮度与梯度信息计算差分影
像,并将其视为带权无向图,采取Prim算法求其
最优生成树,在生成树上从叶子结点到根结点逆
序遍历直接生成镶嵌线,回避了常规的迭代搜索
过程,显著减少了计算量,极大地提高搜索镶嵌
线的效率。另外,为保证镶嵌线附近色彩过渡自
然,采用模版羽化方法对镶嵌线附近的色差进行
过渡处理以获得最终的无缝镶嵌DOM。
1.2技术流程
PxelGrid系统在处理海量卫星遥感数据方
面主要包含卫星及参考数据准备、连接点及参考控
制点自动提取、基于多源控制信息的大规模联合区
域网平差、自动提取DSM、DEM、DOM数据制作及
镶嵌等关键步骤,其整体技术流程如图1所示。
2工程应用
2.1资源三号测区
测区区域:南亚某地区,测区面积约88
万平46北京测绘第32卷第1期
方公里。
测区数据:020张资源三号三线阵立体影
像,实际作业2993张影像;共1078个立体模
型,实际作业1067个模型。删除的卫星数据
包括水域面积占80%以上或质量极差影响处
理的数据,如图2所示。
图1总体技术流程图
图2三线阵测区范围覆盖图
平面参考控制资料为GoogleEarth影像,采
用自动匹配方法共提取464万个控制点。高程
参考控制资料为30米格网的全球SRTM数据。区域网平差过程中自动剔除粗差点,最终剩余
350万个谷歌控制点,平差过程统计如表1所示,
平差后控制点平面中误差为1. 5米,高程中误差
为2.5米。
测区内共6个实测检查点,其中5个检查点
分布于2个模型6张影像上,1个检查点分布于1
个模型3张影像上。各个检查点的精度结果如
表2所示,其中,检查点的平面误差为5.65米,高
程误差为2.55米。
2.2多源数据测区
测区区域:亚洲某地区,测区面积近150万
平方公里。
测区数据:多源数据包括1766景资源三号
下视影像、2598景高分一号PMS1/2影像。该测
区使用的平面参考资料为全球15米分辨率
ETM+正射影像,高程参考资料为全球30米格
网间距ASTER DEM数据。
全自动连接点及控制点提取后共提取2515
万个连接点,840万个参考点。自动剔除粗差及
区域网平差后单位权中误差为0. 35像素,正射
纠正后影像间接边中误差为0. 41像素,最大接
边误差为2.1像素,如图3所示。
图3自动连接点提取效果图
表1统计数据步骤控制点数总连接点数检查点数 迭代次数最终中误差/像素删点阈值/像素自由网平差093115740 70.270无区域网平差464059993115746 120.7493
表2精度检查结果单位:米点号方向误差Y方向误差平面误差高程误差1-4.083.865.621.732-0.416.306.31-2.6032. 696. 206.74-1.484一 1.90-0.071.903.215-1.515.355.563.2162. 760.952.92一
0.19