Compensation Model of GF-7 Panchromatic and Multispectral Image Registration Error
-
摘要: 高分七号是中国第一颗面向1∶10 000立体测图的遥感卫星，其后视全色影像与多光谱影像通过融合可以得到最高0.65 m分辨率的多光谱融合影像。但在实际生产过程中，发现融合影像存在局部晕边现象。针对这一问题，分析了高分七号卫星影像的误差来源，并提出了全色-多光谱配准误差补偿模型，包含用于补偿线性误差的线性项、用于补偿子线阵电荷耦合器件安置误差和镜头畸变的分段四次多项式以及用于补偿周期性误差的三角级数项。根据该模型设计了求解模式和复用模式两种融合方案。利用中国广东地区和青海地区的高分七号影像数据进行实验验证，结果表明所提模型可以使融合配准精度由最大0.6多光谱像素提升至约0.07多光谱像素，平面误差小于0.25倍多光谱像素的区域占比由不足50%提升到约98%，完全满足融合影像生产要求。Abstract:Objectives Gaofen-7 (GF-7) is the first 1∶10 000 stereo mapping satellite of China that can produce high-resolution fusion products. The fusion products are fused with the 0.65 m-panchromatic images and 2.6 m-multispectral images. However, in some fusion products, halos and shifts are existed.Methods To deal with this problem, we analyze the cross-sensor error property of push-broom optical satellite and deduce a compensation model that induces a linear term to compensate for the internal and external orientation errors, a quartic piecewise polynomial term to rectify the lens-distortion and the sub-CCD(charge coupled device) placement error, and a trigonometric series term to correct for the errors caused by platform jitter.Results Several scenes of GF-7 (6 in Guangdong and 3 in Qinghai) were used to demonstrate the compensation model. The results show that the registration RMSE(root mean square error) decreased from 0.6 multispectral-pixels to about 0.07 multispectral-pixels by solving and applying the compensation model, and proportion of areas which have less-than 0.25 multispectral-pixels registration error increased to about 98%, which meets the requirement of producing high-quality fusion products.Conclusions The quantitative results not only reveal the feasibility to improve the accuracy of GF-7 fusion products used by first-level data and achieve the theoretical limit accuracy, but also demonstrate the stability of the piecewise polynomial term allowing for its repeated application to further enhance compensation efficiency.
-
Keywords:
- GF-7 /
- image registration /
- image fusion /
- distortion calibration /
- error compensation /
- multispectral image
-
作为最为广泛使用的电离层经验模型之一，国际参考电离层模型（international reference ionosphere, IRI）对电离层研究和应用具有十分重要的意义[1]。IRI模型利用了垂测仪、非相干散射雷达、卫星观测及探空火箭数据，可以输出任意指定时间和经纬度的60~2 000 km的电子密度和总电子含量（total electron content, TEC）等多种电离层参数信息[2]。自1972年以来，随着电离层观测数据的积累、观测技术的革新和建模方法的改进，IRI模型已发布了多个版本，目前最新版本为IRI-2016 [3-4]。
峰值密度NmF2和峰值高度hmF2是表征电子密度廓线形态的重要参量。对于NmF2的计算，IRI-2016模型中包含IRI_CCIR和IRI_URSI两种模型选项。对于hmF2的计算，在旧版本模型使用的IRI_CCIR hmF2模型的基础上，IRI-2016模型中另外新增了IRI_AMTB和IRI_Shubin两种模型选项。hmF2的IRI_CCIR模型是基于1954—1958年全球分布的约150个垂测仪观测资料，利用hmF2与参数M(3 000)F2的强相关关系由M(3 000)F2计算得到。其中，M(3 000)F2为3 000 km高度处最大可用频率和F2层临界频率（foF2）的比值。这种间接计算hmF2的方法有诸多缺陷。IRI-2016模型中新增的两种计算hmF2的模型选项直接对hmF2进行建模，不依赖M(3 000)F2。其中，IRI_AMTB模型是根据全球分布的26个垂测仪（其中包括中国区域的武汉站、安阳站和三亚站）在1998—2006年的hmF2观测值，纯粹利用地基垂测仪观测值进行建模[5]；IRI_Shubin模型的建立则使用了挑战性小卫星有效载荷（Challenging Minisatellite Payload，CHAMP）、重力恢复与气候实验（Gravity Recovery and Climate Experiment，GRACE）和气象、电离层和气候卫星联合观测系统(Constellation Observing System for Meteorology, Ionosphere and Climate，COSMIC)3个掩星任务的观测资料以及分布全球的62个垂测仪于1987—2012年的观测数据[6]。
COSMIC由美国国防部和中国台湾地区联合研制并于2006-04-15发射升空，在掩星观测任务中具有里程碑的意义[7]。COSMIC由6颗低轨卫星组成，轨道高度在800 km左右，轨道倾角为72°，并配备能够接收GPS信号的接收机。掩星事件发生时，GPS信号穿过电离层时受电子密度影响而发生折射，产生附加相位延迟。在信号直线传播和球对称假设等前提下，由传播路径上的总电子含量使用Abel积分变换方法可得到电子密度廓线。COSMIC掩星反演结果已经得到了多种独立观测手段的验证[8-10]。
中国区域（70°E~140°E，15°N~55°N）跨越中低纬地区，电离层变化幅度较大。文献[11]利用中国区域的4个垂测仪台站对IRI-2016模型在中国区域的性能进行了评估，但只能对垂测仪台站周围的电离层特征进行验证分析。COSMIC掩星观测具有全球覆盖、高精度、垂直分辨率高等特点，可用于对IRI-2016模型在整个中国区域的表现进行分析。本文的研究内容包括使用2008—2014年的COSMIC掩星数据对地磁平静条件下IRI-2016模型在中国区域的性能进行验证，并对IRI-2016模型在中国区域使用时不同NmF2和hmF2模型选项的输出结果进行统计分析，以确定IRI-2016模型在中国区域输出值的最优选项。
1 数据与模型
1.1 COSMIC掩星数据
COSMIC掩星二级产品“ionPrf”数据可从COSMIC数据存档与分析中心（COSMIC Data Analysis and Archive Center, CDAAC）的网站http://cdaac-www.cosmic.ucar.edu下载得到。该数据包含了掩星事件发生的时间、位置、NmF2、hmF2和电子密度廓线等信息。在分析地磁平静条件下IRI模型的表现时，需要先剔除3 h ap指数大于15的时间段的COSMIC掩星数据[8]。另一方面，尽管CDAAC在反演电子密度廓线时采取了一定的质量检核措施，但通过其检核的产品中仍然存在问题廓线[8]。文献[12]提出根据顶部电离层电子密度廓线的斜率和电子密度波动的平均偏差（mean deviation, MD）来剔除这些问题廓线。本文选择COSMIC电子密度廓线的标准为在420~490 km的斜率小于-0.02；平均偏差不超过1.5。此外，峰值高度hmF2范围应在200~500 km之间。
1.2 IRI-2016模型
IRI-2016模型输出值可从http://irimodel.org/网站获取。该模型中计算NmF2的选项包括IRI_CCIR模型和IRI_URSI模型。其中，IRI_URSI模型为推荐选项。计算峰值高度hmF2的选项包括原有的IRI_CCIR模型以及新增的IRI_AMTB模型和IRI_Shubin模型。其中，IRI_AMTB模型为推荐选项。模型的输入参数包括经纬度和时间。由于掩星观测具有水平漂移，本文选用COSMIC掩星观测的F2层峰值高度处的经纬度和发生时刻作为输入量，通过控制模型选项得到不同选项下输出的NmF2和hmF2。图 1(a)、1(b)分别为当地时间（local time, LT）2018-02-03T13:35和2008-02-08T23:46的COSMIC掩星电子密度廓线和IRI-2016模型不同选项下输出结果的比较图。可以看出，COSMIC电子密度廓线和IRI-2016模型不同选项以及选项之间的输出结果存在着明显差异。
![]()
图 1 COSMIC反演的电子密度廓线和IRI-2016模型不同选项输出值的比较Figure 1. Comparison Among Electron Density Profiles Derived from COSMIC RO and IRI-2016 Model with Different Model Options1.3 F10.7指数
F10.7指数是指10.7 cm（2 800 MHz）波长的太阳辐射通量，单位为sfu，与太阳黑子数之间有强相关性，能很好地表征太阳活动强度。图 2是F10.7指数的月均值在第24太阳活动周2008—2014年期间的变化。可以看出，2008—2010年为太阳活动低年，尤其是2008、2009年，太阳活动达到了极低的水平；之后太阳活动水平逐渐提高，在2011—2014年太阳活动维持在较高水平，为太阳活动高年。
2 结果分析
2.1 IRI-2016模型在中国区域的性能分析
为了了解IRI-2016模型在中国区域的表现，首先采用模型推荐的选项计算NmF2和hmF2，将结果与COSMIC掩星观测进行比较。图 3为2008年和2014年COSMIC掩星观测与IRI-2016模型输出的峰值密度NmF2与峰值高度hmF2的散点图。其中，红色实线代表COSMIC观测值等于IRI预测值，黑色虚线代表使用最小二乘方法得到的拟合线，R和S分别表示两种方法得到的相关系数和拟合线斜率。由图 3(a)和3(b)可以看出，尽管在太阳活动低年2008年和太阳活动高年2014年，COSMIC掩星观测和IRI模型输出的NmF2相关系数非常接近（分别为0.86和0.87），但在2014年数据对较为均匀地围绕红色实线分布，且拟合线和等值线非常接近，斜率高达0.90。而在2008年，数据对较为分散，且较多的数据对分布在等值线下侧，拟合直线与等值线之间存在较大偏差。拟合直线的斜率（0.70）表明在太阳活动低年2008年，相较于COSMIC掩星观测值，IRI-2016模型明显高估了NmF2。与峰值密度NmF2情形类似，图 3(c)和图 3(d)表明, 在2014年IRI-2016模型峰值高度hmF2输出值与COSMIC掩星观测值更为接近，拟合直线几乎与等值线重合。而在2008年，拟合直线的斜率为0.76，远小于2014年的斜率0.97，说明2008年IRI-2016模型的hmF2输出值相对于COSMIC掩星观测存在明显高估。同时，由相关系数可以看出，无论在2008年还是2014年，NmF2的相关系数均高于hmF2，表明IRI-2016模型和COSMIC掩星观测输出的NmF2一致性优于hmF2。
![]()
图 3 COSMIC掩星观测与IRI-2016模型采用推荐选项输出的NmF2与hmF2散点图Figure 3. Scatter Plots of NmF2 and hmF2 from COSMIC RO and IRI-2016 Model with Recommended Model Options为了进一步定量分析太阳活动低年和高年COSMIC掩星观测值和IRI-2016模型输出值的差异，定义了相对中误差绝对值的均值（mean absolute relative error，MARE）M，计算公式为：
$$ M=\frac{\mathop{\sum }^{}\left( \left| I-C \right|/C \right)}{N} $$ 式中，C表示由掩星数据反演的NmF2或者hmF2；I表示对应的由模型计算的NmF2或者hmF2；N为观测数据个数。
表 1给出了2008—2014年的NmF2和hmF2的MARE数值。可以看出，无论是NmF2还是hmF2，太阳活动低年（2008—2010年）的MARE均高于太阳活动高年（2011—2014年），在太阳活动极小年（2008—2009年），IRI-2016模型输出值与COSMIC掩星观测值的差异最大。
表 1 2008—2014年NmF2和hmF2的MARE/%Table 1. MARE of NmF2 and hmF2 During 2008―2014/%统计项 2008年 2009年 2010年 2011年 2012年 2013年 2014年 NmF2 31.112 2 31.481 8 27.836 7 26.7882 21.888 1 24.443 1 25.204 4 hmF2 11.689 0 11.944 2 9.747 9 8.460 8 7.071 9 8.021 6 7.344 6 因此，总体而言，对中国区域，无论是NmF2还是hmF2，以COSMIC掩星数据为参照，IRI-2016模型在太阳活动高年的表现都优于太阳活动低年。
2.2 太阳活动低年COSMIC掩星与IRI-2016模型峰值参数在中国区域的空间分布
为了进一步分析中国区域COSMIC掩星反演的峰值参数与IRI-2016模型输出值间的差异，对比了差异更为明显的太阳活动低年（2008—2010年）峰值参数的空间分布。考虑到电离层活跃水平，分析的时间段选为最具代表性的春秋季的中午时分（LT 12:00—14:00）。在本文中，季节的划分标准为：3月、4月、9月、10月为春秋季，5月—8月为夏季，11月、12月、1月和2月为冬季。图 4为2008—2010年春秋季中午（LT 12:00—14:00）由COSMIC掩星与IRI-2016模型得到的峰值参数NmF2与hmF2的均值在中国区域的空间分布。图 4(a)~4(d)显示，COSMIC掩星观测与IRI模型输出的峰值参数空间分布形态均反映出低纬地区电离层更为活跃，同时观测值与模型输出值之间存在明显差异。图 4(e)、4(g) 表明，对于NmF2，与COSMIC掩星反演结果对比，IRI模型输出值在15°N~30°N的低纬区域明显高估，相对偏差最高可达100%以上；而在30°N~55°N的中纬区域，IRI模型输出的NmF2普遍表现为低估。图 4(f)、4(h)表明，对于hmF2，情况恰恰相反：在30°N~55°N中纬区域，IRI模型的hmF2输出值普遍高于COSMIC观测值，而在15°N~30°N低纬区域，模型输出值普遍低于COSMIC观测值。由图 4（f）可见，在太阳活动低年，IRI-2016模型与COSMIC在中国区域得到的hmF2的绝对偏差一般不超过50 km。
![]()
图 4 2008—2010年春秋季LT 12：00—14：00 COSMIC掩星与IRI-2016模型的NmF2与hmF2在中国区域的空间分布Figure 4. Spatial Distributions of NmF2 and hmF2 over China Region at LT 12:00—14:00 During Equinox Months of 2008—2010 Derived from COSMIC RO and IRI-2016 Model2.3 IRI-2016模型峰值参数不同选项输出结果的比较分析
以上分析中所使用的IRI-2016模型输出值均基于模型的推荐选项得到。除了推荐选项，IRI-2016对NmF2和hmF2的计算均提供其他模型选项，而使用不同选项的模型输出结果存在差异。为了了解不同选项的模型在中国区域的表现，进一步以COSMIC掩星数据为参照，对IRI-2016模型中两种计算NmF2的选项（IRI_CCIR和IRI_URSI模型选项）和3种计算hmF2的选项（IRI_CCIR、IRI_AMTB和IRI_Shubin模型选项）的输出结果进行了比较。图 5给出了不同季节COSMIC掩星与IRI-2016模型两种不同选项得到的NmF2均值的日变化情况，其中图 5(a)~5(c)和图 5(d)~5(f)分别为太阳活动低年和高年的结果。统计中如果掩星事件发生时间与某一整点时刻相差不超过30 min，则该时刻的掩星数据和IRI模型输出结果被用于计算该整点时刻的平均值。可以看出，无论采用哪个选项，NmF2在约LT 04：00达到最小值，之后快速升高，在约LT 14：00达到全天最大值，之后逐渐下降。春秋季NmF2高于其他季节，电离层最为活跃。基于两种不同选项的IRI-2016模型输出值与COSMIC观测值均有相似的日变化趋势，但分别基于两种选项的模型输出结果之间存在一定差异。图 5(a)~5(c)表明采用两种选项的IRI-2016模型输出结果非常接近，但在中午时分，IRI_CCIR模型和IRI_URSI模型均高估了峰值密度，其中IRI_URSI的高估程度更严重。在春秋季和冬季，IRI_CCIR模型和IRI_URSI模型在LT 17：00—05：00则低估了峰值密度，其中IRI_URSI低估程度更严重。图 5(d)~5(f)与图 5(a)~5(c)的对比表明，随着太阳活动水平的增加，NmF2随之增大，但中午时分采用两种选项的IRI-2016模型输出值与COSMIC观测值之间的差异均有减小。特别是在太阳活动高年的夏季，IRI-2016模型和COSMIC观测值全天都符合得很好。而在春秋季LT 16：00—04：00，两种选项的IRI-2016模型输出值均小于COSMIC观测值，其中IRI_URSI的偏离度相对更大。
![]()
图 5 不同季节COSMIC掩星与IRI-2016模型两种不同选项下NmF2的日变化图Figure 5. Diurnal Variations of NmF2 from COSMIC RO and IRI-2016 Model Using Two Different Options for Different Seasons图 6给出了不同季节COSMIC掩星观测与IRI-2016模型3种不同选项输出的hmF2日变化情况。可以看出，3种选项对应的IRI-2016模型输出与COSMIC观测的hmF2呈现相似的日变化趋势，即在LT 06：00―07：00达到最小值，在LT 12：00达到小高峰，并在约LT 24：00LT达到全天最大值。在太阳活动低年，3种计算hmF2的模型选项中，IRI_Shubin模型输出值和COSMIC观测值最为接近，其他两种模型输出值相对于COSMIC观测值均存在高估。IRI_AMTB模型输出值相对于COSMIC观测值的高估现象更为突出，冬季最显著。与太阳活动低年相比较，在太阳活动高年，3种选项对应的hmF2模型输出值与COSMIC观测值的一致性均相对更好。但IRI_AMTB模型输出值在夜间的偏差仍较大。
![]()
图 6 不同季节COSMIC掩星与IRI-2016模型不同选项下hmF2的日变化图Figure 6. Diurnal Variations of hmF2 from COSMIC RO and IRI-2016 Model Using Different Options for Different Seasons3 讨论
本文以COSMIC掩星数据为参照，对IRI-2016模型在中国区域的电离层峰值参数NmF2和hmF2输出值随太阳活动强度、季节和昼夜的变化进行了对比分析。整体而言，模型输出值与COSMIC掩星观测值的相关性较高，但太阳活动低年的相关性低于太阳活动高年。两组数据拟合直线的斜率表明，IRI-2016模型输出值与COSMIC掩星观测值的一致性在太阳活动低年低于太阳活动高年。
文献[13]使用位于磁赤道附近的Jicamarca站垂测仪数据与IRI模型输出值进行验证分析，结果发现模型输出值与垂测仪观测值在太阳活动高年的相关性高于太阳活动低年，这与本文的研究结果一致。这主要是由于太阳活动低年，特别是2008年、2009年为太阳活动非寻常极小年，太阳活动水平极低且持续时间较长，IRI模型输出值会产生较大误差，甚至出现不合理的输出结果[3]。
太阳活动低年COSMIC掩星与IRI-2016模型峰值参数在中国区域的空间分布比较情况显示，虽然两者的空间分布形态基本一致，但仍存在明显差异。在15°N~30°N范围内，IRI-2016模型往往高估NmF2而低估hmF2。在30°N~55°N范围内，IRI-2016模型则存在低估NmF2而高估hmF2的现象。COSMIC掩星与IRI-2016模型之间的差异一方面是由于IRI-2016模型本身的误差，另一方面则是由于COSMIC掩星反演结果存在的误差。在球对称假设的前提下，COSMIC掩星反演的电子密度实际上是掩星切点附近整个掩星区域的平均值。因此，在电子密度梯度较大的区域，COSMIC掩星反演的电子密度与真实大气状态相比可能存在较大差异[14]。中国区域低纬度地区位于赤道电离层异常的驼峰区域附近，电子密度梯度很大，导致COSMIC反演的电子密度廓线相较于实际值可能存在较大偏差。
文献[15]使用中国区域位于同一经度带上的4个垂测仪台站（漠河站、北京站、武汉站和三亚站）数据对COSMIC掩星峰值参数进行了比较验证，发现在漠河站两者符合程度较高，但在武汉站和北京站中午时分COSMIC NmF2观测值高于垂测仪观测值，在三亚站中午时分COSMIC NmF2观测值低于垂测仪观测值。反演存在的误差可能使得在中纬度的武汉站（31.0°N，114.6°E）COSMIC反演的NmF2高于IRI-2016模型值，在低纬度的三亚站（18.3°N，109.6°E）COSMIC反演的NmF2低于IRI-2016模型值。文献[11]使用这4个垂测仪台站数据对IRI-2016模型计算hmF2的不同选项输出值结果进行比较分析，发现IRI-2016模型推荐的IRI_AMTB模型选项在4个台站均高估hmF2。本文结果也表明IRI-2016模型的hmF2输出值在30°N~55°N高于COSMIC掩星观测值。对于15°N~30°N范围，IRI-2016模型的hmF2输出值相对于掩星资料存在低估，原因可能是由于该范围位于赤道电离层异常驼峰附近，掩星反演的hmF2高于实际值[11]。
4 结语
本文对2008—2014年IRI-2016模型在中国区域的电离层峰值密度NmF2和峰值高度hmF2输出值与COSMIC掩星数据进行了对比，并对IRI-2016模型中计算这两个峰值参数的不同选项进行了比较分析。结论如下：
1）IRI-2016模型输出值与COSMIC掩星观测数据的一致性在太阳活动高年比太阳活动低年更好。模型NmF2输出值与COSMIC掩星观测值的相关性高于hmF2输出值与COSMIC掩星观测值的相关性。
2）在太阳活动低年春秋季的中午时段，与COSMIC掩星数据相比，IRI-2016模型在30°N~55°N对NmF2和hmF2分别存在低估和高估现象；在15°N~30°N区域内则恰恰相反。
3）以COSMIC掩星观测值为参照，IRI_CCIR模型和IRI_URSI模型均对NmF2在中午时分存在高估现象，该现象在太阳活动低年更为显著。而在春秋季的LT 17：00—04：00时间段，采用这两个模型选项得到的IRI-2016模型输出值均低于COSMIC观测值。
4）在太阳活动低年，无论任何季节，IRI_CCIR模型和IRI_AMTB模型的hmF2输出值均高于COSMIC掩星观测值，其中IRI_AMTB模型输出值相对于掩星观测的偏差更大，尤其是在冬季。IRI_Shubin模型的hmF2输出值与COSMIC观测值最为接近，即使在太阳活动高年，两者的偏差也小于其他两种模型。这表明COSMIC掩星资料在IRI_Shubin模型改进hmF2输出值中发挥了重要作用。
5）根据本文的研究结果，在中国区域采用IRI-2016计算NmF2和hmF2时，分别推荐使用IRI_CCIR和IRI_Shubin模型，而非官方推荐的IRI_URSI和IRI_AMTB模型。
-
![]()
图 1 全色-多光谱影像同名点映射误差源示意图
Figure 1. Error Sources of Panchromatic-Multispectral Corresponding Point Mapping
![]()
图 5 广东组影像2补偿前后的假彩色融合影像细节对比
Figure 5. Detail Comparison in False Color of Fused Image 2 in Guangdong Group Before and After Compensation
![]()
图 6 广东组影像1三角级数项补偿残差图
Figure 6. Residual Plots of Image 1 in Guangdong Group Before and After Trigonometric Series Term Compensation
![]()
图 7 青海组影像1三角级数项补偿残差图
Figure 7. Residual Plots of Image 1 in Qinghai Group Before and After Trigonometric Series Term Compensation
表 1 实验影像信息说明
Table 1 Detail Information of Test Images
组别 影像 命名 覆盖地区 应用模式 广东组 影像一 GF701_006650_E113.3_N22.5_20210113111729 广东深圳 求解 影像二 GF701_006807_E114.2_N22.7_20210123111416 广东中山 求解 影像三 GF701_006807_E114.3_N22.9_20210123111413 广东深圳 求解 影像四 GF701_007241_E113.9_N22.7_20210221112456 广东东莞 复用 影像五 GF701_007241_E113.9_N22.9_20210221112453 广东惠州 复用 影像六 GF701_008136_E113.9_N22.7_20210421112548 广东深圳 复用 青海组 影像一 GF7_DLC_E101.2_N37.6_20220108_L1A0000677449 青海门源 求解 影像二 GF7_DLC_E101.2_N37.8_20211130_L1A0000639036 青海门源 求解 影像三 GF7_DLC_E101.3_N37.8_20220108_L1A0000677448 青海海北 复用 
下载: 导出CSV
表 2 实验影像不同模型残差补偿精度
Table 2 Compensation Accuracy of Test Images with Different Strategies
组别 影像 统计项 中误差/像素 $ {P}_{\nabla xy\le 0.25} $
/%QNR x(垂轨) y(沿轨) 广东组 影像1 初始 0.622 0.329 0.4 0.846 求解线性项后 0.119 0.202 65.8 0.863 求解分段项后 0.077 0.081 97.7 0.892 影像2 初始 0.627 0.327 0.2 0.811 求解线性项后 0.113 0.201 67.1 0.822 求解分段项后 0.067 0.071 98.6 0.846 影像3 初始 0.628 0.323 0.2 0.804 求解线性项后 0.116 0.200 66.3 0.804 求解分段项后 0.075 0.068 98.6 0.830 青海组 影像1 初始 0.345 0.170 44.5 0.983 求解线性项后 0.097 0.131 91.6 0.986 求解分段项后 0.038 0.099 97.9 0.986 影像2 初始 0.356 0.347 38.4 0.975 求解线性项后 0.121 0.307 75.7 0.977 求解分段项后 0.095 0.278 83.8 0.979
下载: 导出CSV
表 3 复用模式下分段项前后的配准精度
Table 3 Registration Accuracy of Test Images Before and After Piecewise Polynomial Term Compensation
组别 影像 纠正前后 中误差/像素 $ {P}_{\nabla xy\le 0.25} $
/%x(垂轨) y(沿轨) 广东组 影像4 复用纠正前 0.619 0.335 0.2 复用纠正后 0.073 0.073 98.3 影像5 复用纠正前 0.612 0.328 0.4 复用纠正后 0.078 0.072 97.8 影像6 复用纠正前 0.604 0.343 0.5 复用纠正后 0.082 0.075 98.3 青海组 影像7 复用纠正前 0.347 0.288 4.2 复用纠正后 0.044 0.098 98.5
下载: 导出CSV
-
[1] 张永军, 张祖勋, 龚健雅. 天空地多源遥感数据的广义摄影测量学[J]. 测绘学报, 2021, 50(1): 1-11. https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB202101001.htm Zhang Yongjun, Zhang Zuxun, Gong Jianya. Generalized Photogrammetry of Spaceborne, Airborne and Terrestrial Multi-source Remote Sensing Datasets[J]. Acta Geodaetica et Cartographica Sinica, 2021, 50(1): 1-11. https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB202101001.htm
[2] 张永军, 万一, 史文中, 等. 多源卫星影像的摄影测量遥感智能处理技术框架与初步实践[J]. 测绘学报, 2021, 50(8): 1068-1083. https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB202108009.htm Zhang Yongjun, Wan Yi, Shi Wenzhong, et al. Technical Framework and Preliminary Practices of Photogrammetric Remote Sensing Intelligent Processing of Multi-source Satellite Images[J]. Acta Geodaetica et Cartographica Sinica, 2021, 50(8): 1068-1083 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB202108009.htm
[3] 翁永玲, 田庆久. 遥感数据融合方法分析与评价综述[J]. 遥感信息, 2003, 18(3): 49-54. https://www.cnki.com.cn/Article/CJFDTOTAL-YGXX200303017.htm Weng Yongling, Tian Qingjiu. Analysis and Evaluation of Method on Remote Sensing Data Fusion[J]. Remote Sensing Information, 2003, 18(3): 49-54 https://www.cnki.com.cn/Article/CJFDTOTAL-YGXX200303017.htm
[4] 汤志强. 航天线阵CCD传感器严格几何模型的分析与验证[D]. 郑州: 信息工程大学, 2011. Tang Zhiqiang. Analysis and Validation for Rigorous Geometric Model of Space Linear Array CCD Sensors[D]. Zhengzhou: Information Engineering University, 2011.
[5] 雷蓉. 星载线阵传感器在轨几何定标的理论与算法研究[D]. 郑州: 信息工程大学, 2011. Lei Rong. Study on Theory and Algorithm of the In-flight Geometric Calibration of Spaceborne Linear Array Sensor[D]. Zhengzhou: Information Engineering University, 2011.
[6] Zheng M T, Zhang Y J, Zhu J F, et al. Self-calibration Adjustment of CBERS-02B Long-Strip Imagery[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(7): 3847-3854. doi: 10.1109/TGRS.2014.2385706
[7] 孙韬, 龙辉, 赵冬, 等. 基于五谱段合一多光谱相机影像的卫星颤振检测和补偿[J]. 光学学报, 2014, 34(7): 284-290. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB201407042.htm Sun Tao, Long Hui, Zhao Dong, et al. Detection and Compensation of Satelliteflutter Based on Image from Multispectral Camera with Five Spectral Combinations[J]. Acta Optica Sinica, 2014, 34(7): 284-290 https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB201407042.htm
[8] 王密, 朱映, 范城城. 高分辨率光学卫星影像平台震颤几何精度影响分析与处理研究综述[J]. 武汉大学学报(信息科学版), 2018, 43(12): 1899-1908. doi: 10.13203/j.whugis20180240 Wang Mi, Zhu Ying, Fan Chengcheng. Development of Platform Jitter Geometric Analysis and Processing for High-Resolution Optical Satellite Imagery[J]. Geomatics and Information Science of Wuhan University, 2018, 43(12): 1899-1908 doi: 10.13203/j.whugis20180240
[9] Grodecki J, Dial G. Block Adjustment of High-Resolution Satellite Images Described by Rational Polynomials[J]. Photogrammetric Engineering and Remote Sensing, 2003, 69(1): 59-68.
[10] Tadono T. Initial Results of Calibration and Validation for PRISM and AVNIR-2[J]. Asian Journal of Geoinformatics, 2006, 6(4): 11-20.
[11] 许丽娜, 肖奇, 何鲁晓. 考虑人类视觉特征的融合图像评价方法[J]. 武汉大学学报(信息科学版), 2019, 44(4): 546-554. doi: 10.13203/j.whugis20170168 Xu Lina, Xiao Qi, He Luxiao. Fused Image Quality Assessment Based on Human Visual Characteristics[J]. Geomatics and Information Science of Wuhan University, 2019, 44(4): 546-554 doi: 10.13203/j.whugis20170168
[12] 梁丽娟, 黄万里, 张容焱, 等. Sentinel-2卫星影像融合方法与质量评价分析[J]. 遥感技术与应用, 2019, 34(3): 612-621. https://www.cnki.com.cn/Article/CJFDTOTAL-YGJS201903018.htm Liang Lijuan, Huang Wanli, Zhang Rongyan, et al. Sentinel-2 Satellite Image Fusion Method and Quality Evaluation Analysis[J]. Remote Sensing Technology and Application, 2019, 34(3): 612-621 https://www.cnki.com.cn/Article/CJFDTOTAL-YGJS201903018.htm
[13] 曹海翊, 张新伟, 赵晨光, 等. 高分七号卫星总体设计与技术创新[J]. 中国空间科学技术, 2020, 40(5): 1-9. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKJ202005002.htm Cao Haiyi, Zhang Xinwei, Zhao Chenguang, et al. System Design and Key Technolongies of the GF-7 Satellite[J]. Chinese Space Science and Technology, 2020, 40(5): 1-9 https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKJ202005002.htm
[14] 王智均, 李德仁, 李清泉. 利用小波变换对影像进行融合的研究[J]. 武汉测绘科技大学学报, 2000, 25(2): 137-142. https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200002009.htm Wang Zhijun, Li Deren, Li Qingquan. Image Fusion with Wavelets Analysis[J]. Journal of Wuhan Technical University of Surveying and Mapping, 2000, 25(2): 137-142. https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200002009.htm
[15] 李平湘, 王智均. 高空间分辨率影像和多光谱影像融合的研究[J]. 武汉大学学报(信息科学版), 2003, 28(1): 14-17. http://ch.whu.edu.cn/article/id/4732 Li Pingxiang, Wang Zhijun. Investigation of Image Fusion Between High Resolution Image and Multi-spectral Image[J]. Geomatics and Information Science of Wuhan University, 2003, 28(1): 14-17. http://ch.whu.edu.cn/article/id/4732
[16] 李均力, 贾永红. 基于分辨率退化模型的全色和多光谱遥感影像融合方法[J]. 武汉大学学报(信息科学版), 2004, 29(4): 302-306. http://ch.whu.edu.cn/article/id/4618 Li Junli, Jia Yonghong. A New Fusion Method Based on Resolution Degradation Model for Panchromatic and Multi-spectral Images[J]. Geomatics and Information Science of Wuhan University, 2004, 29(4): 302-306 http://ch.whu.edu.cn/article/id/4618
[17] 陈雪洋, 袁超. ZY-102C卫星影像融合方法评价[J]. 测绘与空间地理信息, 2013, 36(2): 50-53. https://www.cnki.com.cn/Article/CJFDTOTAL-DBCH201302015.htm Chen Xueyang, Yuan Chao. Data Fusion Evaluation of ZY-1 02C Satellite Images[J]. Geomatics and Spatial Information Technology, 2013, 36(2): 50-53 https://www.cnki.com.cn/Article/CJFDTOTAL-DBCH201302015.htm
[18] Xie G Q, Wang M, Zhang Z Q, et al. Near Real-Time Automatic Sub-pixel Registration of Panchromatic and Multispectral Images for Pan-Sharpening[J]. Remote Sensing, 2021, 13(18): 3674.
[19] 李德仁. 我国第一颗民用三线阵立体测图卫星: 资源三号测绘卫星[J]. 测绘学报, 2012, 41(3): 317-322. https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201503009.htm Li Deren. China?s First Civilian Three-Line-Array Stereo Mapping Satellite: ZY-3[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(3): 317-322 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201503009.htm
[20] 童小华, 叶真, 刘世杰. 高分辨率卫星颤振探测补偿的关键技术方法与应用[J]. 测绘学报, 2017, 46(10): 1500-1508. https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201710033.htm Tong Xiaohua, Ye Zhen, Liu Shijie. Essential Technology and Application of Jitter Detection and Compensation for High Resolution Satellites[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(10): 1500-1508 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201710033.htm
[21] Sun J M, Shen Z H, Wang Y A, et al. LoFTR: Detector-Free Local Feature Matching with Transformers[C]//IEEE/CVF Conference on Computer Vision and Pattern Recognition, Nashville, USA, 2021.
-
期刊类型引用(2)
1. 陆浩楠，叶世榕，张庆兰，夏朋飞，鄂盛龙. 利用COSMIC-1/2 IonPrf产品分析中国区域近15年电离层时空变化特征. 武汉大学学报(信息科学版). 2024(05): 765-774 . 
百度学术
2. 唐丝语，周春红，黄智，袁洪. 国际参考电离层模型输出参数峰值高度性能. 空间科学学报. 2021(06): 898-904 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 562
- HTML全文浏览量: 90
- PDF下载量: 135
- 被引次数: 4
