基于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