《测绘学报》
构建与学术的桥梁 拉近与权威的距离
复制链接,关注《测绘学报》抖音!
【测绘学报的个人主页】长按复制此条消息,长按复制打开抖音查看TA的更多作品##7NsBSynuc88##[抖音口令]
本文内容来源于《测绘学报》2020年第10期,审图号GS(2020)5551号。
高分辨率光学遥感卫星反射点源像点提取方法
徐伟伟1,2,张黎明1,2,李鑫1,2,杨宝云1,2
1.中国科学院安徽光学精密机械研究所, 安徽 合肥 230031;2.中国科学院通用光学定标与表征技术重点实验室, 安徽 合肥 230031
基金项目:国家自然科学基金(41601388)
摘要:像点提取精度是高分辨光学遥感卫星几何定标精度的关键因素之一。本文提出一种基于轻小型、自动化反射点源的像点提取方法,以参数化高斯模型拟合获取点源影像的像点坐标。在轨试验结果初步表明,反射点源法所提取的像点坐标的共线误差小于0.001像素,比例常量的相对标准偏差优于5‰,与质心法等常用方法提取的点源像点坐标相比较,精度优于0.05像素。反射点源不仅能够实现高精度的像点提取,而且还能够实现高分辨率光学遥感卫星的高精度几何定标与像质评价,对提高我国航天遥感立体测绘精度具有重要意义。
关键词:几何定标 像点提取 反射点源 共线检验 比例检验
![d890e7260838df56ad272e4c88da6ab7.png](https://i-blog.csdnimg.cn/blog_migrate/659a4e9d58a551c7a43766b449301cc1.jpeg)
![88601d5b12c7118b8fb7ce9e223cd559.png](https://i-blog.csdnimg.cn/blog_migrate/4e687c6f5e4edde55cb062d972cbd2a4.jpeg)
引文格式:徐伟伟, 张黎明, 李鑫, 等. 高分辨率光学遥感卫星反射点源像点提取方法. 测绘学报,2020,49(10):1295-1302. DOI: 10.11947/j.AGCS.2020.20190246.
阅读全文:http://xb.sinomaps.com/article/2020/1001-1595/20201006.htm
全文概述
高分辨率光学遥感卫星以一种非常精细的方式观测地表,所获取的遥感影像不仅可以更加清楚地表达地物目标的空间结构与表层纹理特征,而且发展到有效解译分析地球系统多圈层的物理参数与演变规律。高分辨率遥感影像在自然资源调查、生态环境保护、城市规划建设、军事侦察及测绘制图等国民经济建设与国防安全各领域有着非常广阔的应用前景及重要意义[1]。然而,高分辨率遥感影像应用的广度与深度有赖于定标技术的支持,需要定期或不定期开展高分辨率光学遥感卫星在轨几何定标,确定相机在轨运行状态的内外方位元素、安置参数、镜头畸变参数等[2],以实现高精度定位与几何质量。
高分辨率光学遥感卫星在轨几何定标的核心任务之一就是通过检测来求得相机高精度的内外方位元素。这涉及控制点影像的像点坐标提取,由于像点的位置信息包含了相机内部的几何畸变等,不可能通过改变像点位置关系的平差来提高内方位元素的检校精度,只能通过高精度的像点提取或统计平差来实现[3-4]。目前采用自然地物,如道路交叉口、房屋拐角等,作为地面控制点的测量精度已达厘米甚至毫米量级,星上位姿测量的精度也在此量级,然而自然地物控制点对应影像的像点坐标提取主要是人工选点或像点测量算法得到。人工选点依靠主观判断,误差可达0.5~1像素,商业测量软件虽然采用数字图像处理算法,但是一般的中误差也在0.3像素左右[5-8],无法与地面坐标测量精度相匹配。通常是采用多次测量统计平差来提高其精度,但存在检校链路长,平台是否稳定的问题,同时要求控制点数量多且密集,相应地增加地面测量的任务量。从地面控制点-卫星-相机三位一体的立体测绘来说,控制点影像的像点提取精度极大地限制了几何定标结果的精度与准确性。无人机遥感制作高精度DEM或DOM,通过图像匹配提供丰富的大量的地面控制点信息,但由于DOM与卫星遥感影像获取的时间同步、比例尺、配准及退化等差异,需要附加几何约束来实现图像匹配,例如手动选取少量的图像控制点,而图像控制点通常也是选自地物目标(或者地面控制点),故均匀分布的地面控制点是高精度几何定标的基本条件[9-11],另外地形地貌的高精度DEM或DOM,实时更新的代价也比较高。因此提出以轻量化、小型化、自动化的反射点源作为高分辨率光学遥感卫星在轨几何检校的地面控制点,以期提高控制点遥感影像的像点提取精度,进而提高高分辨率光学遥感卫星的在轨几何定标精度。
1 基本原理
光学遥感卫星空间分辨率的不断提高,可以在地面设置空间尺度足够“小”,辐射亮度足够“亮”的反射点源来对星载光学相机进行在轨几何定标,以检校其内外方位元素及位姿等参数。反射点源能够将近似平行入射的太阳光发散一定的角度,便于卫星观测成像,如图 1所示,根据参数设计与太阳光束0.53°发散角,由几何光学追迹与数值计算分析,光学遥感卫星所接收的反射光斑仅为点源目标镜面上尺度为厘米量级的一小区域,相对米量级分辨率的卫星相机,反射点源可作为点激冲目标。
![36a871d908e316bf52fee5b9f01ac8d6.png](https://i-blog.csdnimg.cn/blog_migrate/2f0867af1109587f1202c29e0cbc7fe5.jpeg)
图 1 反射点源法成像关系Fig. 1 Imaging relationship of the reflected point source method
图选项
根据光学原理及傅里叶光学可知,对于线性位移不变成像系统,光学遥感卫星的成像关系可表示为
(1)
式中,g(x,y)是遥感影像;f(x,y)是地物目标场景;*是卷积运算符;h(x,y)是系统点扩散函数;b是背景值。
当目标场景输入为点激冲目标(如恒星、反射点源)时[3, 12-13],由卷积定理可得光学卫星相机的成像关系可表示为
(2)
此时,遥感器成像系统的输出表现为其自身的特性即点扩散特性。
光学遥感卫星成像系统主要由前置光学子系统、电子学子系统、焦平面CCD探测器等部分,根据系统组成模型分析与相关报道IKONOS、Eros-B等高分辨率相机的恒星观测以及SPOT、Quickbird等高分相机对地面点目标成像时,遥感成像系统的点扩散特性均采用高斯模型来表示[12-15],故高分辨率光学遥感卫星的点扩散函数可近似采用高斯模型来表达
![8e87666f5f43322eb4b9473269450ddc.png](https://i-blog.csdnimg.cn/blog_migrate/9377fc7c9ac2eb5211c08aeb396dd887.jpeg)
(3)
式中,k是系数因子;(x0,y0)是中心位置即像点坐标;σ与ζ是标准偏差。
光学遥感卫星成像系统对地面反射点源目标的遥感影像可表示为
![ea766dc9e15a87925701d058b09b4459.png](https://i-blog.csdnimg.cn/blog_migrate/4c37c957abb80ed366f8c8cc02e48c57.jpeg)
(4)
由遥感成像系统的物像关系可知,光学卫星相机对反射点源的响应值是表征成像系统特性的点扩散值与背景值之和。据此,结合在轨检测的地面反射点源目标设置及其遥感影像,对于5×5像素的反射点源影像计数值,根据式(4)可列25个方程,其中6个未知数待解,通过最小二乘法进行二维高斯拟合,得反射点源遥感影像的像点坐标(x0,y0)
![7d88fe2234e3b6f328735e2684a7bfa7.png](https://i-blog.csdnimg.cn/blog_migrate/ac4bc0f3c66fd9daf81397b040e571e8.jpeg)
(5)
式中,g(xi, yj)是像素坐标为(xi, yj)的反射点源影像计数值;(x0,y0)是反射点源遥感影像的像点坐标。
2 试验与分析
2.1 数据处理
利用研制的反射点源目标于2015年5月在陕西西安对资源三号卫星开展了在轨成像质量评价试验。如图 2所示,分别沿星载遥感器的垂轨与顺轨方向集中布设了非整数像素间隔的4×4反射点源阵列(按照从左到右,从上到下的顺序标记为P1、P2、…、P15、P16),相邻反射点源间距10.25像素,相间反射点源间距20.5像素。专业测绘人员采用RTK-GPS结合地理信息动态采集系统的差分数据对反射点源目标中心的几何位置进行高精度测量,得到各反射点源几何位置的坐标数据,据此计算得各点源中心的实际间距,见表 1。
![e73016d78652ac97883497ee0b5770d2.png](https://i-blog.csdnimg.cn/blog_migrate/dd988f5bd8df79e2656feb15b0d4d4de.jpeg)
图 2 在轨试验Fig. 2 On-orbit experiment
图选项
表 1 地面反射点源中心间距Tab. 1 Center distances of the reflected point source on ground