遥感影像的时间序列(论文推荐徐伟伟)
《测绘学报》
构建与学术的桥梁 拉近与权威的距离
复制链接,关注《测绘学报》抖音!
【测绘学报的个人主页】长按复制此条消息,长按复制打开抖音查看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像素。反射点源不仅能够实现高精度的像点提取,而且还能够实现高分辨率光学遥感卫星的高精度几何定标与像质评价,对提高我国航天遥感立体测绘精度具有重要意义。
关键词:几何定标 像点提取 反射点源 共线检验 比例检验
引文格式:徐伟伟, 张黎明, 李鑫, 等. 高分辨率光学遥感卫星反射点源像点提取方法. 测绘学报,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°发散角,由几何光学追迹与数值计算分析,光学遥感卫星所接收的反射光斑仅为点源目标镜面上尺度为厘米量级的一小区域,相对米量级分辨率的卫星相机,反射点源可作为点激冲目标。
图 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],故高分辨率光学遥感卫星的点扩散函数可近似采用高斯模型来表达
(3)
式中,k是系数因子;(x0,y0)是中心位置即像点坐标;σ与ζ是标准偏差。
光学遥感卫星成像系统对地面反射点源目标的遥感影像可表示为
(4)
由遥感成像系统的物像关系可知,光学卫星相机对反射点源的响应值是表征成像系统特性的点扩散值与背景值之和。据此,结合在轨检测的地面反射点源目标设置及其遥感影像,对于5×5像素的反射点源影像计数值,根据式(4)可列25个方程,其中6个未知数待解,通过最小二乘法进行二维高斯拟合,得反射点源遥感影像的像点坐标(x0,y0)
(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。
图 2 在轨试验Fig. 2 On-orbit experiment
图选项
表 1 地面反射点源中心间距Tab. 1 Center distances of the reflected point source on ground
方向 | 线段 | 距离/m | 线段 | 距离/m | 线段 | 距离/m | 线段 | 距离/m | 线段 | 距离/m | 线段 | 距离/m |
垂轨 | P1P2 | 35.875 0 | P2P3 | 35.863 5 | P3P4 | 35.877 9 | P1P3 | 71.738 4 | P2P4 | 71.741 3 | P1P4 | 107.616 3 |
P5P6 | 35.871 0 | P6P7 | 35.867 6 | P7P8 | 35.870 2 | P5P7 | 71.738 6 | P6P8 | 71.737 6 | P5P8 | 107.608 7 | |
P9P10 | 35.851 9 | P10P11 | 35.913 4 | P11P12 | 35.876 1 | P9P11 | 71.765 3 | P10P12 | 71.789 3 | P9P12 | 107.641 1 | |
P13P14 | 35.861 0 | P14P15 | 35.874 9 | P15P16 | 35.869 7 | P13P15 | 71.735 9 | P14P16 | 71.744 6 | P13P16 | 107.605 5 | |
沿轨 | P1P5 | 35.878 2 | P5P9 | 35.874 5 | P9P13 | 36.262 9 | P1P9 | 71.751 7 | P5P13 | 72.137 2 | P1P13 | 108.013 5 |
P2P6 | 35.843 2 | P6P10 | 35.905 1 | P10P14 | 36.252 8 | P2P10 | 71.747 9 | P6P14 | 72.157 7 | P2P14 | 108.000 2 | |
P3P7 | 35.854 9 | P7P11 | 35.896 3 | P11P15 | 36.237 6 | P3P11 | 71.750 8 | P7P15 | 72.133 7 | P3P15 | 107.988 4 | |
P4P8 | 35.909 1 | P8P12 | 35.860 3 | P12P16 | 36.238 8 | P4P12 | 71.769 3 | P8P16 | 72.099 0 | P4P16 | 108.008 1 |
表选项
作为地面控制点的轻小型自动化反射点源主要由反射镜、太阳观察器、电控经纬仪及辅助设备等组成,可将入射至反射镜的太阳光反射至卫星入瞳以形成点激冲辐射,太阳反射波段的镜面反射率优于80%,空间指向精度优于0.1°。优化设计的反射镜能够反射适量的太阳光通量,以与卫星相机的动态范围相匹配并处于其动态范围的高端且不饱和;通过太阳观察器对太阳运行轨迹的观测实现反射点源的角度定标;进而结合光学遥感卫星的过顶参数,控制经纬仪自动实现太阳-反射点源-卫星相机三者间的光路自动配准,以将入射至点源目标的太阳光反射至卫星相机入瞳,准确获取反射点源遥感影像。
根据反射点源法像点提取原理,由式(5)所述,采用参数化二维高斯模型对于5×5像素的每个反射点源影像数据进行最小二乘法的二维曲面拟合,如图 3所示,以实现对反射点源目标遥感影像的像点位置坐标提取,即同时获取沿轨方向与垂轨方向的反射点源像点坐标,像点提取结果见表 2,X0与Y0即为16个反射点源影像的像点坐标。
图 3 像点提取Fig. 3 Pixel extraction
图选项
表 2 点源影像像点坐标Tab. 2 Pixel coordinate of point source images
坐标/像素 | P1 | P2 | P3 | P4 | P5 | P6 | P7 | P8 |
X0 | 633.788 4 | 644.889 7 | 656.135 4 | 667.288 8 | 633.670 5 | 644.736 7 | 655.958 0 | 667.137 5 |
Y0 | 554.142 8 | 554.335 1 | 554.491 9 | 554.699 1 | 564.265 8 | 564.450 0 | 564.711 7 | 564.909 8 |
坐标/像素 | P9 | P10 | P11 | P12 | P13 | P14 | P15 | P16 |
X0 | 633.451 8 | 644.581 9 | 655.745 2 | 666.954 3 | 633.505 2 | 644.502 7 | 655.620 0 | 666.710 9 |
Y0 | 574.443 3 | 574.610 9 | 574.885 8 | 575.051 3 | 584.903 4 | 584.885 2 | 585.131 3 | 585.298 2 |
表选项
反射点源在地面的布设设计为各目标中心横向或纵向连线为一条直线,据此根据闭环检验或共线检验原则,即卫星遥感影像上垂轨或沿轨方向相邻反射点源中心点距离之和应与相间反射点源中心点距离相等,见表 3。共线检验结果表明,无论在垂轨方向或沿轨方向,点源像点共线误差均小于0.001像素,所测卫星相机探元范围内,探测元均保持了良好的刚性与稳定性。
表 3 共线检验结果Tab. 3 Results of collineation checkout
方向 | 闭环关系 | 结果 | 闭环关系 | 结果 | 闭环关系 | 结果 | 闭环关系 | 结果 |
垂轨 | P1P2 P2P3-P1P3 | 3.18×10-5 | P1P2 P2P4-P1P4 | 4.26×10-6 | P1P3 P3P4-P1P4 | 3.25×10-5 | P2P3 P3P4-P2P4 | 6.01×10-5 |
P5P6 P6P7-P5P7 | 1.24×10-4 | P5P6 P6P8-P5P8 | 5.57×10-5 | P5P7 P7P8-P5P8 | 1.95×10-4 | P6P7 P7P8-P6P8 | 8.78×10-4 | |
P9P10 P10P11-P9P11 | 2.55×10-4 | P9P10 P10P12-P9P12 | 7.95×10-5 | P9P11 P11P12-P9P12 | 9.64×10-5 | P10P11 P11P12-P10P12 | 2.72×10-4 | |
P13P14 P14P15-P13P15 | 1.56×10-4 | P13P14 P14P16-P13P16 | 1.50×10-4 | P13P15 P15P16-P13P16 | 8.31×10-4 | P14P15 P15P16-P14P16 | 1.39×10-4 | |
沿轨 | P1P5 P5P9-P1P9 | 2.46×10-4 | P1P5 P5P13-P1P13 | 4.49×10-5 | P1P9 P9P13-P1P13 | 1.62×10-4 | P5P9 P9P13-P5P13 | 1.82×10-4 |
P2P6 P6P10-P2P10 | 2.99×10-8 | P2P6 P6P14-P2P14 | 4.57×10-5 | P2P10 P10P14-P2P14 | 1.90×10-4 | P6P10 P10P14-P6P14 | 1.45×10-4 | |
P3P7 P7P11-P3P11 | 3.22×10-5 | P3P7 P7P15-P3P15 | 2.21×10-6 | P3P11 P11P15-P3P15 | 1.63×10-4 | P7P11 P11P15-P7P15 | 1.93×10-4 | |
P4P8 P8P12-P4P12 | 2.68×10-5 | P4P8 P8P16-P4P16 | 1.27×10-4 | P4P12 P12P16-P4P16 | 1.82×10-4 | P8P12 P12P16-P8P16 | 8.24×10-5 |
表选项
光学遥感卫星同一台相机对地面的两个反射点源目标成像时,根据成像几何关系可知,点源目标的地面间距与其相应卫星遥感影像的像点间距之比为常量。据此由各反射点源影像提取的像点坐标可得其像点间距,进而与表 1中相应地面间距求得卫星相机的比例常量,结果见表 4,垂轨与顺轨方向的比例常量的相对标准偏差均优于5‰,所测卫星相机探测元保持了良好刚性与稳定性,比例检验结果与共线检验具有较好的一致性。
表 4 比例检验结果Tab. 4 Results of scale checkout
方向 | 线段 | 比例 | 线段 | 比例 | 线段 | 比例 | 线段 | 比例 | 线段 | 比例 | 线段 | 比例 |
垂轨 | P1P2 | 3.231 1 | P1P3 | 3.209 8 | P1P4 | 3.216 2 | P2P3 | 3.188 8 | P2P4 | 3.202 4 | P3P4 | 3.216 2 |
P5P6 | 3.241 0 | P5P7 | 3.218 1 | P5P8 | 3.214 8 | P6P7 | 3.195 5 | P6P8 | 3.201 8 | P7P8 | 3.208 1 | |
P9P10 | 3.220 8 | P9P11 | 3.218 5 | P9P12 | 3.212 4 | P10P11 | 3.216 2 | P10P12 | 3.208 2 | P11P12 | 3.200 3 | |
P13P14 | 3.260 8 | P13P15 | 3.243 6 | P13P16 | 3.240 3 | P14P15 | 3.226 2 | P14P16 | 3.230 0 | P15P16 | 3.233 8 | |
标准偏差 | 0.017 0 | 平均值 | 3.219 0 | 相对标准偏差 | 5‰ | |||||||
方向 | 线段 | 比例 | 线段 | 比例 | 线段 | 比例 | 线段 | 比例 | 线段 | 比例 | 线段 | 比例 |
沿轨 | P1P5 | 3.544 0 | P1P9 | 3.534 0 | P1P13 | 3.511 3 | P5P9 | 3.524 1 | P5P13 | 3.495 3 | P9P13 | 3.466 7 |
P2P6 | 3.543 2 | P2P10 | 3.538 2 | P2P14 | 3.534 9 | P6P10 | 3.533 2 | P6P14 | 3.530 8 | P10P14 | 3.528 4 | |
P3P7 | 3.507 8 | P3P11 | 3.517 6 | P3P15 | 3.524 0 | P7P11 | 3.527 4 | P7P15 | 3.532 1 | P11P15 | 3.536 7 | |
P4P8 | 3.516 4 | P4P12 | 3.525 9 | P4P16 | 3.529 2 | P8P12 | 3.535 4 | P8P16 | 3.535 5 | P12P16 | 3.535 6 | |
标准偏差 | 0.016 9 | 平均值 | 3.525 3 | 相对标准偏差 | 4.8‰ |
表选项
2.2 分析讨论
基于反射点源的高分辨率光学遥感卫星像点提取方法,采用参数化高斯模型表征卫星相机成像系统的点扩散特性,进而通过最小二乘法曲线拟合来提取像点坐标。反射点源法在轨几何定标精度的提高,有赖于作为地面控制点的遥感影像像点提取精度,需要考虑作为地面控制点的点源结构形式,将像点坐标与常用的非线性求解法、质心法或梯度法等多种方法提取的点源像点坐标相比较,以验证该点源及其像点提取方法的可行性与有效性。
作为地面控制点的轻小型、自动化的反射点源既有相对空间分辨率小得多的尺寸,又有相对地物背景高得多的亮度,便于光学遥感卫星对其观测成像且具有较高的信噪比,在自然地物目标中不易寻找。反射点源便于定标场的固定布设或移动加密,且具有高精度的地面测量精度。非线性求解法是根据光学遥感卫星系统点扩散函数可分离的假设[14],结合获取的反射点源遥感影像,采用非线性方程优化求解得系统一维离散点扩散函数即线扩散函数[3],进而高斯拟合得点源影像像点坐标Ⅰ,见表 5,Ⅱ是质心法点源像点提取结果,Ⅲ是梯度法像点提取结果,通过与反射点源法像点提取结果比较分析点源像点的提取精度。各种方法的点源像点提取结果具有较好的一致性,反射点源法像点提取结果与这些方法的像点提取结果差异,见表 6,最大误差优于0.05像素,验证了基于反射点源的高分辨率光学遥感卫星像点提取方法的可行性与有效性,并且点源像点提取精度优于自然地物目标(最优误差0.3像素)与方形十字对顶角标志(误差优于0.12像素)[6, 16],因此,轻小型自动化的反射点源作为地面控制点能够实现高分辨率光学遥感卫星的高精度几何定标。
表 5 像点提取结果Tab. 5 Results of pixel extraction
坐标/像素 | P1 | P2 | P3 | P4 | ||||||||
Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | |
X0 | 633.798 5 | 633.790 7 | 633.785 1 | 644.894 0 | 644.874 5 | 644.839 9 | 656.135 1 | 656.117 8 | 656.082 9 | 667.290 4 | 667.289 4 | 667.289 2 |
Y0 | 554.149 2 | 554.144 8 | 554.142 5 | 554.342 2 | 554.337 1 | 554.334 1 | 554.490 2 | 554.485 5 | 554.474 5 | 554.700 8 | 554.692 0 | 554.676 1 |
坐标/像素 | P5 | P6 | P7 | P8 | ||||||||
Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | |
X0 | 633.677 5 | 633.672 4 | 633.669 1 | 644.746 0 | 644.742 5 | 644.744 7 | 655.960 1 | 655.959 2 | 655.959 5 | 667.136 9 | 667.139 0 | 667.142 5 |
Y0 | 564.268 3 | 564.265 5 | 564.262 4 | 564.446 4 | 564.456 1 | 564.472 0 | 564.711 3 | 564.710 3 | 564.707 9 | 564.909 1 | 564.902 7 | 564.889 3 |
坐标/像素 | P9 | P10 | P11 | P12 | ||||||||
Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | |
X0 | 633.449 4 | 633.459 7 | 633.478 0 | 644.584 3 | 644.593 5 | 644.614 3 | 655.753 1 | 655.751 3 | 655.755 5 | 666.957 9 | 666.951 8 | 666.943 2 |
Y0 | 574.443 2 | 574.452 6 | 574.471 4 | 574.610 4 | 574.611 7 | 574.613 9 | 574.883 6 | 574.885 1 | 574.886 0 | 575.052 4 | 575.049 9 | 575.046 0 |
坐标/像素 | P13 | P14 | P15 | P16 | ||||||||
Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | |
X0 | 633.506 8 | 633.499 6 | 633.502 9 | 644.504 4 | 644.505 3 | 644.508 8 | 655.626 1 | 655.629 4 | 655.642 2 | 666.722 1 | 666.716 2 | 666.715 5 |
Y0 | 584.909 3 | 584.909 4 | 584.901 6 | 584.885 1 | 584.891 3 | 584.903 8 | 585.133 2 | 585.130 6 | 585.127 4 | 585.305 0 | 585.301 9 | 585.302 4 |
表选项
表 6 像点坐标误差Tab. 6 Errors of pixel coordinate
误差/像素 | P1 | P2 | P3 | P4 | ||||||||
Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | |
ΔX0 | 0.010 1 | 0.002 3 | -0.003 3 | 0.004 3 | -0.015 2 | -0.049 8 | -0.000 3 | -0.017 6 | -0.049 9 | 0.001 6 | 0.000 6 | 0.000 3 |
ΔY0 | 0.006 4 | 0.002 0 | -0.000 3 | 0.007 1 | 0.002 0 | -0.001 0 | -0.001 7 | -0.006 4 | -0.017 4 | 0.001 7 | -0.007 1 | -0.023 0 |
误差/像素 | P5 | P6 | P7 | P8 | ||||||||
Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | |
ΔX0 | 0.007 0 | 0.001 9 | -0.001 4 | 0.009 3 | 0.005 8 | 0.008 0 | 0.002 1 | 0.001 2 | 0.001 5 | -0.000 6 | 0.001 5 | 0.005 0 |
ΔY0 | 0.002 4 | -0.000 3 | -0.003 4 | -0.003 6 | 0.006 1 | 0.022 0 | -0.000 4 | -0.001 4 | -0.003 8 | -0.000 7 | -0.007 1 | -0.020 5 |
误差/像素 | P9 | P10 | P11 | P12 | ||||||||
Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | |
ΔX0 | -0.002 4 | 0.007 9 | 0.026 2 | 0.002 4 | 0.011 6 | 0.032 4 | 0.007 9 | 0.006 1 | 0.010 3 | 0.003 6 | -0.002 5 | -0.011 1 |
ΔY0 | -0.000 1 | 0.009 3 | 0.028 1 | -0.000 5 | 0.000 8 | 0.003 0 | -0.002 2 | -0.000 7 | 0.000 1 | 0.001 1 | -0.001 4 | -0.005 3 |
误差/像素 | P13 | P14 | P15 | P16 | ||||||||
Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | Ⅰ | Ⅱ | Ⅲ | |
ΔX0 | 0.001 6 | -0.005 6 | -0.002 3 | 0.001 7 | 0.002 6 | 0.006 1 | 0.006 1 | 0.009 4 | 0.022 2 | 0.011 2 | 0.005 3 | 0.004 6 |
ΔY0 | 0.005 9 | 0.006 0 | -0.001 8 | -0.000 1 | 0.006 1 | 0.018 6 | 0.001 9 | -0.000 7 | -0.003 9 | 0.006 8 | 0.003 7 | 0.004 2 |
表选项
反射点源影像像点坐标检验过程中,采用两个点源的地面距离与其像素距离的比值即比例检验来验证像点位置坐标的准确性,见表 4。该比例常量即为光学遥感卫星相机的地面像元分辨率,横向或纵向分辨率检测的相对标准偏差均优于5‰,具有较好的一致性。相对于光学遥感卫星相机空间分辨率在轨检测的周期性三线靶标或大面积辐射状靶标[17-18],反射点源不仅具有小型化、轻量化、自动化的特点,便于机动运输及展开或者固定布设,而且还能够定量地获取光学遥感卫星相机的地面像元分辨率。基于反射点源的在轨MTF检测方法与其遥感影像的像点提取,见表 2,将所有点源影像数据进行位置配准以获取亚像素插值的遥感器成像系统点扩散轮廓,进而以参数化高斯模型拟合得系统点扩散函数,如图 4所示,对其进行离散傅里叶变换取模并归一化处理得系统MTF。相较于其他在轨MTF检测方法[19-22],如刃边法,反射点源法是根据MTF的物理定义进行光学遥感卫星在轨MTF检测的直接方法。因此,轻小型自动化的反射点源作为地面参照目标能够实现高分辨率光学遥感卫星的高精度像质评价[23-24]。
图 4 点源法像质评价Fig. 4 Image quality evaluation by point source method
图选项
3 结论
高分辨率光学遥感卫星像点提取方法以轻量化、小型化、自动化的反射点源作为在轨几何定标的地面控制点,通过共线检验与比例检验来验证点源像点坐标的准确性,并将像点坐标与质心法或梯度法等常规像点提取结果相比较,来评估反射点源法像点提取精度。在轨试验结果初步表明,反射点源法提取的像点坐标的共线误差小于0.001像素,比例常量的相对标准偏差优于5‰,所测卫星相机探测元保持了良好刚性与稳定性。像点坐标与质心法等常用方法提取的点源像点坐标相比较,反射点源法像点提取精度优于0.05像素,相对于自然地物目标与方形十字对顶角标志,具有较高的像点提取精度,同时还能够进行光学遥感卫星空间分辨率与MTF的在轨检测。光学遥感卫星反射点源法几何定标的应用研究正在开展,目前从反射点源控制点-卫星星上位姿-相机点源影像像点提取等三位一体的高精度立体测绘来说,点源影像像点坐标提取精度的提高,将有可能实现与地面及星上测量相当的高精度,进而将有利于提高光学遥感卫星在轨几何定标精度。轻量化、小型化的反射点源便于网络化布设应用,这将有利于实现光学遥感卫星的自动化几何定标,提高定标效率。因此,轻小型、自动化的反射点源目标能够实现高分辨率光学遥感卫星高精度的自动化几何定标与像质评价,对提高我国立体测绘的定量化水平具有重要意义。
作者简介
第一作者简介:徐伟伟(1983-), 男, 博士, 副研究员, 研究方向为高分辨率光学遥感卫星辐射定标与像质评价等。E-mail: weilxu@aiofm.ac.cn
论文推荐 | 曾添, 隋立芬, 阮仁桂, 等:卫星精密定轨的三频观测量IF组合法
综述 | 郭树人, 蔡洪亮, 孟轶男, 等:北斗三号导航定位技术体制与服务性能
论文推荐 | 周乐韬, 黄丁发, 袁林果, 等:基于状态和残差的北斗基准站观测数据表达与信息分级
回顾经典 | 王虎, 党亚民, 侯阳飞, 等:PPP网解UPD模糊度固定的无基站差分大型CORS站整网快速精密解算
权威 | 专业 | 学术 | 前沿
微信、抖音小视频投稿邮箱 | song_qi_fan@163.com
欢迎加入《测绘学报》作者QQ群: 751717395
进群请备注:姓名 单位 稿件编号
,免责声明:本文仅代表文章作者的个人观点,与本站无关。其原创性、真实性以及文中陈述文字和内容未经本站证实,对本文以及其中全部或者部分内容文字的真实性、完整性和原创性本站不作任何保证或承诺,请读者仅作参考,并自行核实相关内容。文章投诉邮箱:anhduc.ph@yahoo.com