-
GF-7 卫星激光测高系统配备有两台激光足印相机(Laser Footprint Camera,LFC)、一台激光光轴监视相机(Laser Optical Axis Surveillance Camera,LOASC)。如图1所示,两台足印相机分别服务于两个激光器,在同步曝光模式下,足印相机同时对激光发射光斑和地物进行成像;光轴监视相机对光路内经多次扩束、反射得到的两个激光光斑成像,根据两者的位置关系实时监测仪器框架、温度等因素的变化。
如图2所示,从左至右依次为激光足印影像1&2(Laser Footprint Image, LFI)、激光监视相机影像(Laser Surveillance-camera Image, LSI),大小均为550×550 pixel。其中,红色框选区域为对应影像上波束1的激光光斑,蓝色框选区域为对应影像上波束2的激光光斑,足印影像外围的“小光斑”是曝光时激光在面阵相机上的虚像。传统质心提取方法多是针对点光源、背景单一的影像,而足印影像同时对光斑、背景地物成像,能拍摄到云、砂石和冰川等具有高反射率的地物,这类地物的强度值与激光光斑强度值相叠加,增强了光斑的饱和,直接使用传统算法难以有效提取光斑质心。对于监视相机影像而言,由于背景单一,传统算法也能在质心提取方面取得较好的效果。
文中使用自然资源部国土卫星遥感应用中心提供的激光测高卫星标准测高产品SLA03作为实验数据,将GF-7在轨实测光斑质心坐标作为中心,向四周扩展20 pixel形成40×40的足印影像、监视影像切片。将影像左上角作为原点,定义像平面坐标系,对于足印影像而言,波束1光斑切片左上角像平面坐标为(102,242)、波束2光斑切片左上角像平面坐标为(198,143);对于监视影像而言,波束1光斑切片左上角像平面坐标为(40,291),波束2光斑切片左上角像平面坐标为(208,138)。
-
灰度重心法(Gray Centroid Method, GCM)是点光源质心提取经典算法,根据目标区域内灰度分布赋予不同的权重,计算光斑质心坐标,即使发射光斑形状并非严格的圆或椭圆型也能取得较好的效果,且计算复杂度较低,但是直接使用该方法无法避免足印影像中高反射率地物的影响。计算方式如公式(1)、(2),其中(X,Y)是该光斑的质心坐标,I(i,j)代表影像上第i行第j列的灰度值,M、N代表影像行列的最大值、最小值。
$$X = \dfrac{{\displaystyle\sum\limits_{i = 1}^M {\displaystyle\sum\limits_{j = 1}^N {i\times I(i,j)} } }}{{\displaystyle\sum\limits_{i = 1}^M {\displaystyle\sum\limits_{j = 1}^N {I(i,j)} } }} $$ (1) $$Y = \dfrac{{\displaystyle\sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {j\times I(i,j)} } }}{{\displaystyle\sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {I(i,j)} } }} $$ (2) -
对于足印影像而言,光斑以及周围区域会受到地表地物的影响,高反射率地物会影响传统算法结果。文中提出一种阈值约束的椭圆拟合算法(Threshold constrained for extracting centroid of Ellipse Fitting Method, TEFM),如图3所示,分为以下几步。
(1)去除背景地物影响。提取光斑所在位置的足印影像切片,每个像元可视为光斑与地物相互叠加的效果。经过统计大量数据,发现背景区域灰度分布直方图往往在2 000~3000范围内聚集。切片影像去除地物影响,整体灰度值幅值减去2 000,如图4(a)所示;
图 4 阈值约束的椭圆拟合质心提取算法
Figure 4. Threshold constraint algorithm for extracting centroid of ellipse fitting
(2)基于阈值法初步提取光斑轮廓。利用公式(3)计算阈值T、确定光斑的初提取结果,作为第一次约束,将小于阈值的像元赋值为0,图4(b)为处理结果,存在一部分云区域被误识别:
$$ T = {1 \mathord{\left/ {\vphantom {1 {{{\rm e}^2}}}} \right. } {{\rm e^2}}} \times \max \left(\sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {I(i,j) - 2\;000} } \right) $$ (3) (3)形态学处理去除噪声。对初提取结果进行形态学腐蚀处理,作为第二次约束,去除细小噪声、孔隙、地物细节纹理特征的影响,避免对后续步骤造成影响,很大程度上保留、贴近光斑形状。如公式(4)所示,腐蚀处理为提取每个像素邻域D1内局部最小值的过程。图4(c)为处理结果,如图所示,腐蚀处理去除了红色框选区域的像元。
$$ {I_{dst}}(i,j) = \min \sum\limits_{i,j \in {D_1}} {{I_{src}}(i,j)} $$ (4) (4)椭圆拟合约束光斑形状。将腐蚀结果代入公式(5)、利用最小二乘进行椭圆拟合[16],求解目标函数f的最小值确定各系数。将椭圆轮廓作为第三次约束,将边界外的像元赋值为0,并保留光斑特征参数,如图4(d)所示;
$$ f(A,B,C,D,E) = \sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {(A{i^2} + Bij + C{j^2} + Di + Ej + F} } {)^2} $$ (5) (5)利用灰度重心法实现光斑质心提取。将步骤(4)中的结果作为输入,利用公式(1)、(2)计算质心坐标;
(6)基于光斑特征参数去除粗差。足印影像上的地物信息很有可能会掩盖光斑信息,甚至无法提取光斑轮廓形状,此情况会在业务化处理中带来大量的误差,因此根据偏心率、长短半轴等光斑特征去辨别是否成功提取光斑信息,借此去除误差。由于出光角度、仪器内部框架等因素,理想光斑偏心率往往在0.2~0.8之间。经过实际检校,激光光斑直径在17~20 m(波束1 19 m、波束2 21 m)之间、对应5~7 pixel,但是为了捕捉更多的光斑轮廓信息,将这一阈值扩大到理想发散光斑边缘,取值为20。
-
激光定位精度高度依赖于指向角精度,而指向角变化难以直接通过激光产品直接得到。通过激光光斑质心坐标变化可以分析指向角相对精度,间接监测指向角稳定度,为激光产品提供关于指向角变化相关的质量控制参数。质心坐标变化与激光指向变化之间的关系,由以下准则确定:足印相机视场角为0.19°,激光脉冲经过4倍扩束到达地面形成足印光斑,对应每像元分辨率为1.5 μrad,所以光斑质心坐标每偏移一个像素,指向角变化约0.31″。对于光轴监视相机影像而言,主要作用是基于两波束光斑质心位置的相对关系监测光学仪器内部框架是否变化。之后,会围绕足印影像、监视影像之间的几何关系,进一步探索激光两波束之间相对关系,建立光轴几何关系变化模型,实现激光光轴夹角的变化监测。
为了验证文中方法精度的可靠性,如图5所示分别采用夜间升轨、背景为复杂地物的足印影像作为实验数据,红、蓝、黄色圆点分别对应质心标定位置、文中方法质心提取结果、灰度重心法质心提取结果。如图5(a)所示,夜间条件、避免复杂背景的情况下,文中方法质心提取精度在0.05 pixel以内,灰度重心法质心提取精度在0.08 pixel左右;如图5(b)所示,在复杂地物影响下,文中质心提取精度在0.08 pixel以内,优于灰度重心法结果约2.5 pixel。灰度重心法的误差主要来源于光斑轮廓附近的曝光色散、背景地物的影响,而文中方法通过多次约束去除了这些干扰因素。在背景地物更为复杂的条件下,灰度重心法与文中算法的差距还会更大。
为了检验文中所提出的TEFM算法精度,随机抽选时间相近、地物类型丰富的足印影像共1 600张进行精度验证(见表1)。将两波束在轨标定位置作为真实值,分别计算坐标误差均值(Mean)、坐标误差极值(Range)、坐标分布均方根误差(RMSE),统计结果表明:传统算法在复杂背景下基本不具备提取光斑质心能力,文中采用的TEFM算法在稳定性、整体提取精度均保持较高水平,优于传统算法一个量级,满足复杂背景条件下亚像素级的质心提取需求。
表 1 足印影像质心提取算法精度评价(单位: pixel)
Table 1. Accuracy evaluation of centroid extraction algorithm for footprint image(Unit: pixel)
Method Mean Range RMSE X Y X Y X Y GCM 2.565 2.522 4.941 5.185 2.512 2.218 TEFM 0.062 0.071 0.092 0.088 0.052 0.081 -
对2020年3月15日至30日共计61轨足印相机影像数据进行质心提取,获取足印影像8.6万张,其中有27%的影像由于被云遮盖、过曝等原因丢失光斑轮廓信息,无法进行质心提取。实验结果如图6所示,横坐标表示X方向、纵坐标表示Y方向。图6(a)为波束1光斑质心坐标统计结果,坐标均值为(120.83,263.71),X方向在0.4 pixel左右来回震荡,Y方向在0.5 pixel左右来回震荡,平面位置变化在1 pixel以内;图6(b)为波束2光斑质心坐标统计结果,坐标均值为(216.32,162.51),X方向在0.5 pixel左右来回震荡,Y方向在0.5 pixel左右来回震荡,平面位置变化在1 pixel以内。可以得出初步结论:(1)相比于传统算法,文中所改进的算法能有效提取复杂地物背景下的光斑质心;(2)在每月内质心坐标变化比较稳定,在X、Y方向均小于1 pixel的幅度震荡。
图 6 波束1(a) & 2(b)激光足印影像光斑质心坐标统计
Figure 6. Spot centroid coordinate statistics of LFI of laser 1 (a) and laser 2 (b)
建立长时间序列激光指向监测系统,长时间序列质心提取的难点在于剔除粗差的同时保留质心坐标本身的变化趋势。尽管文中采用了多个条件约束光斑轮廓,精确定位光斑质心位置,但背景地物对光斑的影响仍不能完全去除,带来粗差,使光斑质心坐标出现跳变,破坏本有的坐标变化趋势。出现粗差的情形主要有以下两种:(1)高反射率地物完全遮盖了光斑轮廓;(2)光斑本身幅值较低、比较黯淡,光斑特征被地物特征淹没。经过大量统计,发现偏心率、半轴长度等光斑特征参数能有效识别粗差,提升整体识别精度。为了分析足印光斑质心变化及稳定性,从宏观、微观两个角度展开分析。
如表2所示,为了分析足印影像光斑质心在轨运行以来长周期内的稳定性,统计了2020年3月至2021年4月每月光斑质心坐标均值。如图7(a)、(b)所示,对于波束1而言,质心坐标在X方向上整体呈现减小趋势,变化约0.4 pixel,在Y方向上整体呈现减小趋势,变化约1 pixel,平面位置变化约1.1 pixel,对应指向角变化约0.341″;如图7(c)、(d)所示,对于波束2而言,质心坐标在X
方向上呈现先减后增大的趋势,变化约0.4 pixel,在Y方向上整体呈现减小趋势、变化约1.5 pixel,平面位置变化约1.4 pixel,对应指向角变化约0.434″。实际上,不论是波束1还是波束2,每月光斑质心变化不是朝着某一方向一直移动的,每月变化幅度在一个非常小的范围内。 表 2 激光足印影像光斑质心坐标每月均值统计结果(单位:pixel)
Table 2. Statistical results of monthly mean of spot centroid coordinates of LFI (Unit: pixel)
Time LFI Laser 1 Laser 2 X Y X Y 2020 Mar. 120.83 263.71 216.32 162.51 Apr. 120.50 264.51 216.88 162.63 May 120.35 263.88 216.85 162.56 Jun. 120.37 263.79 216.87 162.36 Jul. 120.52 263.90 216.93 162.10 Aug. 120.63 263.28 216.81 161.65 Sep. 120.51 262.26 216.72 161.53 Oct. 120.11 261.85 216.35 161.60 Nov. 120.33 262.02 216.42 161.64 Dec. 120.37 262.30 216.31 161.82 2021 Jan. 120.49 262.55 216.00 161.65 Feb. 120.36 262.55 216.81 161.17 Mar. 120.45 262.67 216.62 161.68 Apr. 120.43 262.63 216.63 161.07 -
此节使用2.2节中对应的激光光轴监视相机数据,波束1、2质心坐标统计结果如图8所示,横坐标表示X方向、纵坐标表示Y方向。相比于足印影像,监视影像质心提取算法较少受外界因素影响,仅使用传统算法就能获取精度较高的结果。左侧为波束1光斑质心坐标统计结果,坐标均值为(57.18, 313.10),X方向在0.9 pixel左右来回震荡,Y方向在0.5 pixel左右来回震荡,平面位置变化在1 pixel以内;右侧为波束2光斑质心坐标统计结果,坐标均值为(228.01,156.67),X方向在1 pixel左右来回震荡,Y方向在0.7 pixel左右来回震荡,平面位置变化在1.2 pixel以内。相比于激光足印影像,即使在每月内监视影像质心也会发生较大的偏移,为了保证均值在数据中具有代表性,之后将每月数据中比例较小的异常点排除,不计入统计。
图 8 波束1(a) & 2(b)激光监视影像光斑质心坐标统计
Figure 8. Centroid coordinate statistics of LSI of laser 1 (a) and laser 2 (b)
如表3所示,为了分析激光监视影像光斑质心在轨运行以来长周期内的稳定性,统计了2020年3月至2021年4月每月光斑质心坐标均值、影像上两光斑间距离变化,距离变化以3月波束1、2光斑质心之间的距离作为基准。对于波束1而言,长周期内X方向减小了0.7 pixel,Y方向减小了0.2 pixel,平面位置变化约0.7 pixel;对于波束2而言,长周期内X方向增大了0.8 pixel,Y方向减小了0.4 pixel,平面位置变化约0.9 pixel。如图9所示,两光斑间距先增大再减小,可能呈现一定的周期性变化。两光斑的距离变化表示仪器内部框架发生了一定变化,可能受到温度等内应力因素影响。
表 3 激光监视影像光斑质心坐标每月均值统计结果(单位:pixel)
Table 3. Statistical results of monthly mean of spot centroid coordinates of LSI (Unit: pixel)
Time LSI Distance between
two spots in the LSILaser 1 Laser 2 X Y X Y 2020 Mar. 57.18 313.10 228.01 156.67 0 Apr. 56.93 313.35 228.63 156.31 1.05 May 56.51 313.39 228.84 156.14 1.66 Jun. 56.41 313.35 228.79 156.09 1.71 Jul. 56.52 313.32 228.81 156.29 1.48 Aug. 56.08 313.68 228.77 156.18 2.09 Sep. 55.99 313.75 228.68 156.22 2.11 Oct. 55.77 313.63 228.82 156.62 2.03 Nov. 56.08 314.031 228.71 156.22 2.26 Dec. 55.93 313.74 228.67 156.50 1.95 2021 Jan. 56.22 313.62 228.56 156.62 1.49 Feb. 56.50 313.51 228.79 156.49 1.47 Mar. 56.47 313.62 228.59 156.56 1.38 Apr. 56.48 312.99 228.82 156.25 1.33
Analysis on the change of GF-7 satellite laser altimeter spot centroid position
-
摘要: 高分七号卫星搭载了我国首个具备全波形记录能力的对地观测激光测高仪,可以大范围地获取高精度三维坐标,其定位精度高度依赖于激光指向角的测量精度。针对数据特点,提出了一种阈值约束的椭圆拟合光斑质心提取算法,并建立了长周期指向角稳定性监测与分析系统。首先,用阈值法确定激光光斑轮廓的边缘;其次,通过腐蚀操作消除孔隙、噪声的影响;然后,通过椭圆拟合进一步约束激光光斑的形状,保留光斑的特征参数;最后,通过灰度重心法提取质心坐标。实验结果表明:光斑的质心位置在1.4 pixel内变化,对应的指向角度每月在0.434″内变化,较为稳定。文中相关算法、结论对我国后续激光测高卫星研制、指向稳度监测具有一定的借鉴意义。Abstract: GF-7 satellite is equipped with the first laser earth observation altimeter with full waveform recording ability in China, which can obtain a wide range of high-precision three-dimensional coordinates, and its positioning accuracy is highly dependent on the measurement accuracy of laser pointing angle. According to the characteristics of data, a centroid extraction algorithm of ellipse fitting spot with threshold constraint was proposed, and a long-period pointing angle stability monitoring system was established. Firstly, the threshold method was used to determine the edge of the laser spot contour. Secondly, the influence of porosity and noise was eliminated through corrosion operation. Then, the shape of the laser spot was further constrained by ellipse fitting, and the characteristic parameters of the laser spot were preserved. Finally, the centroid coordinates were extracted by using the gray gravity center method. The experimental results show that the centroid position of the spot changes within 1.4 pixel, and the corresponding pointing angle changes within 0.434″ per month, which is relatively stable. The related algorithms and conclusions in this paper have certain reference significance for the development of laser altimetry satellites and the monitoring of pointing stability in China.
-
表 1 足印影像质心提取算法精度评价(单位: pixel)
Table 1. Accuracy evaluation of centroid extraction algorithm for footprint image(Unit: pixel)
Method Mean Range RMSE X Y X Y X Y GCM 2.565 2.522 4.941 5.185 2.512 2.218 TEFM 0.062 0.071 0.092 0.088 0.052 0.081 表 2 激光足印影像光斑质心坐标每月均值统计结果(单位:pixel)
Table 2. Statistical results of monthly mean of spot centroid coordinates of LFI (Unit: pixel)
Time LFI Laser 1 Laser 2 X Y X Y 2020 Mar. 120.83 263.71 216.32 162.51 Apr. 120.50 264.51 216.88 162.63 May 120.35 263.88 216.85 162.56 Jun. 120.37 263.79 216.87 162.36 Jul. 120.52 263.90 216.93 162.10 Aug. 120.63 263.28 216.81 161.65 Sep. 120.51 262.26 216.72 161.53 Oct. 120.11 261.85 216.35 161.60 Nov. 120.33 262.02 216.42 161.64 Dec. 120.37 262.30 216.31 161.82 2021 Jan. 120.49 262.55 216.00 161.65 Feb. 120.36 262.55 216.81 161.17 Mar. 120.45 262.67 216.62 161.68 Apr. 120.43 262.63 216.63 161.07 表 3 激光监视影像光斑质心坐标每月均值统计结果(单位:pixel)
Table 3. Statistical results of monthly mean of spot centroid coordinates of LSI (Unit: pixel)
Time LSI Distance between
two spots in the LSILaser 1 Laser 2 X Y X Y 2020 Mar. 57.18 313.10 228.01 156.67 0 Apr. 56.93 313.35 228.63 156.31 1.05 May 56.51 313.39 228.84 156.14 1.66 Jun. 56.41 313.35 228.79 156.09 1.71 Jul. 56.52 313.32 228.81 156.29 1.48 Aug. 56.08 313.68 228.77 156.18 2.09 Sep. 55.99 313.75 228.68 156.22 2.11 Oct. 55.77 313.63 228.82 156.62 2.03 Nov. 56.08 314.031 228.71 156.22 2.26 Dec. 55.93 313.74 228.67 156.50 1.95 2021 Jan. 56.22 313.62 228.56 156.62 1.49 Feb. 56.50 313.51 228.79 156.49 1.47 Mar. 56.47 313.62 228.59 156.56 1.38 Apr. 56.48 312.99 228.82 156.25 1.33 -
[1] Shu Rong, Kong Wei. Current status and development of space-based active opto-electronics remote sensing technology [J]. Aerospace Shanghai, 2019, 36(3): 1-14. (in Chinese) [2] Abdalati W, Zwally H J, Bindschadler R. The ICESat-2 laser altimetry mission [J]. Proceedings of the IEEE, 2010, 98(5): 735-751. doi: 10.1109/JPROC.2009.2034765 [3] Abshire J B, Sun X L, Riris H, et al. Geoscience laser altimeter system (GLAS) on the ICESat mission: pre-launch and on-orbit measurement performance[C]// IEEE International Geoscience & Remote Sensing Symposium. IEEE, 2003. [4] Li Guoyuan, Yao Jiaqi, Zhao Yiming, et al. Progress and prospect of atmospheric scattering correction for laser altimetry satellite [J]. Infrared and Laser Engineering, 2020, 49(11): 20200234. (in Chinese) doi: 10.3788/IRLA20200234 [5] Li Guoyuan, Tang Xinming. Analysis and validation of ZY-3 02 satellite laser altimetry data [J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(12): 1939-1949. (in Chinese) [6] Tang X M, Xie J F, Liu R, et al. Overview of the GF-7 laser altimeter system mission [J]. Earth and Space Science, 2019, 7(2): 1-14. [7] Zuo Zhiqiang, Tang Xinming, Li Guoyuan, et al. Adaptive Gaussian filtering of the full waveform of GF-7 satellite laser altimeter [J]. Infrared and Laser Engineering, 2020, 49(11): 20200251. (in Chinese) doi: 10.3788/IRLA20200251 [8] Guo Jinquan, Li Guoyuan, Zuo Zhiqiang, et al. Full waveform data quality and characteristic analysis of GF-7 satellite laser altimeter [J]. Infrared and Laser Engineering, 2020, 49(S2): 20200387. (in Chinese) [9] Zhu Hong, Xie Junfeng, Sun Guangtong, et al. Centroid location accuracy analysis of spaceborne laser spot [J]. Remote Sensing Information, 2021, 36(2): 1-5. [10] Fan Chunbo, Li Jiancheng, Wang Dan, et al. ICESAT/GLAS laser footprint geolocation and error analysis [J]. Journal of Geodesy and Geodynamics, 2007, 27(1): 104-106. (in Chinese) [11] Bae S, Webb C, Schutz B, et al. GLAS PAD calibration using laser reference sensor data[C]//AIAA/AAS Astrodynamics Specialist Conference & Exhibit, 2004. [12] Sirota J M, Bae S, Millar P, et al. The transmitter pointing determination in the geoscience laser altimeter system [J]. Geophysical Research Letters, 2005, 32(22): L22S11. doi: https://doi.org/10.1029/2005GL024005 [13] Waerbeke L V, Mellier V, Erben T, et al. Detection of correlated galaxy ellipticities on CFHT data: First evidence for gravitational lensing by large-scale structures [J]. Astronomy & Astrophysics, 2000, 358(1): 30-44. [14] Quine B M, Tarasyuk V, Mebrahtu H, et al. Determining star-image location: A new sub-pixel interpolation technique to process image centroids [J]. Computer Physics Communi-cations, 2007, 177(9): 700-706. doi: 10.1016/j.cpc.2007.06.007 [15] Wang Haiyong, Zhou Wenrui, Lin Haoyu, et al. Parameter estimation of Gaussian gray diffusion model of static image spot [J]. Acta Optica Sinica, 2012, 32(3): 0323004. (in Chinese) doi: 10.3788/AOS201232.0323004 [16] Yuan Xiaoqi, Li Guoyuan, Tang Xinming, et al. Centroid automatic extraction of spaceborne laser spot image [J]. Acta Geodaetica et Cartographica Sinica, 2018, 47(2): 135-141. (in Chinese) [17] Yang Xiongdan, Li Guoyuan, Wang Peixian, et al. Laser pointing changes detection method for space-borne laser spot image [J]. Acta Geodaetica et Cartographica Sinica, 2020, 49(12): 86-94. (in Chinese) [18] Deng Yongtao, Li Xu, Yan Fanjiang, et al. Attitude deter-mination method and error analysis of space-borne laser range finder [J]. Spacecraft Recovery & Remote Sensing, 2015, 36(3): 86-91. (in Chinese) doi: 10.3969/j.issn.1009-8518.2015.03.012 [19] Hilbert C, Schmullius C. Influence of surface topography on ICESat/GLAS forest height estimation and waveform shape [J]. Remote Sensing, 2012, 4(8): 2210-2235. doi: 10.3390/rs4082210