无数据
Advanced Search
Your Location：
Home >
Browse articles >
Compensation Model of GF-7 Panchromatic and Multispectral Image Registration Error
| Updated：2023-06-29
|
Compensation Model of GF-7 Panchromatic and Multispectral Image Registration Error
- Geomatics and Information Science of Wuhan University
- Affiliations：
1.School of Remote Sensing and Information Engineering, Wuhan University, Wuhan 430079, China
2.Land and Resources Technical Center of Guangdong Province, Guangzhou 510062, China
- Author bio：
- Funds：the Foundation for Innovative Research Groups of the Natural Science Foundation of Hubei Province(2020CFA003)
DOI：10.13203/j.whugis20220714
CLC： P236Published：05 July 2023，
Received：04 November 2022，
Scan for full text
Cite this article
张永军,王梦欣,万一，等.高分七号卫星影像融合中的全色-多光谱配准误差补偿模型[J].武汉大学学报（信息科学版）,2023,48(7):1029-1038.
ZHANG Yongjun,WANG Mengxin,WAN Yi,et al.Compensation Model of GF-7 Panchromatic and Multispectral Image Registration Error[J].Geomatics and Information Science of Wuhan University,2023,48(7):1029-1038.
张永军,王梦欣,万一，等.高分七号卫星影像融合中的全色-多光谱配准误差补偿模型[J].武汉大学学报（信息科学版）,2023,48(7):1029-1038. DOI：10.13203/j.whugis20220714
ZHANG Yongjun,WANG Mengxin,WAN Yi,et al.Compensation Model of GF-7 Panchromatic and Multispectral Image Registration Error[J].Geomatics and Information Science of Wuhan University,2023,48(7):1029-1038. DOI：10.13203/j.whugis20220714
Compensation Model of GF-7 Panchromatic and Multispectral Image Registration Error
高分七号卫星影像融合中的全色-多光谱配准误差补偿模型 Compensation Model of GF-7 Panchromatic and Multispectral Image Registration Error 1 first-author 张永军，博士，教授，主要从事航空航天摄影测量、实景三维重建、遥感解译与空间知识图谱相关研究。zhangyj@whu.edu.cn 张永军，博士，教授，主要从事航空航天摄影测量、实景三维重建、遥感解译与空间知识图谱相关研究。zhangyj@whu.edu.cn zhangyj@whu.edu.cn 1 1 2 2 2 1 武汉大学遥感信息工程学院，湖北 武汉，430079 武汉大学遥感信息工程学院，湖北 武汉，430079 School of Remote Sensing and Information Engineering, Wuhan University, Wuhan 430079, China School of Remote Sensing and Information Engineering, Wuhan University, Wuhan 430079, China 广东省国土资源技术中心，广东 广州，510062 广东省国土资源技术中心，广东 广州，510062 Land and Resources Technical Center of Guangdong Province, Guangzhou 510062, China Land and Resources Technical Center of Guangdong Province, Guangzhou 510062, China 高分七号是中国第一颗面向1∶10 000立体测图的遥感卫星，其后视全色影像与多光谱影像通过融合可以得到最高0.65 m分辨率的多光谱融合影像。但在实际生产过程中，发现融合影像存在局部晕边现象。针对这一问题，分析了高分七号卫星影像的误差来源，并提出了全色-多光谱配准误差补偿模型，包含用于补偿线性误差的线性项、用于补偿子线阵电荷耦合器件安置误差和镜头畸变的分段四次多项式以及用于补偿周期性误差的三角级数项。根据该模型设计了求解模式和复用模式两种融合方案。利用中国广东地区和青海地区的高分七号影像数据进行实验验证，结果表明所提模型可以使融合配准精度由最大0.6多光谱像素提升至约0.07多光谱像素，平面误差小于0.25倍多光谱像素的区域占比由不足50%提升到约98%，完全满足融合影像生产要求。 Objectives 2 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 2 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 2 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 2 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. 高分七号 影像配准 影像融合 畸变检校 误差补偿 多光谱影像 GF-7 image registration image fusion distortion calibration error compensation multispectral image 影像融合是光学遥感影像处理中的重要过程 ［ 1 1 - 2 2 1-2 ］ 。光学遥感卫星一般以高分辨率全色影像与低分辨率多光谱影像相组合的形式获取数据，在地面处理中进行影像融合 ［ 3 3 ］ ，既获得了高分辨率的多波段图像，又节约了宝贵的卫星通信带宽、星上存储空间和信号处理功耗。在影像融合之前，一般需要先进行全色影像与多光谱影像的配准，称为融合配准。这是因为受全色光学传感器与多光谱光学传感器的电荷耦合器件（charged coupled device， CCD）安置精度和镜头畸变 ［ 4 4 - 5 5 4-5 ］ 、外方位元素误差 ［ 6 6 ］ 、卫星平台颤振 ［ 7 7 - 8 8 7-8 ］ 等因素的影响 ［ 9 9 ］ ，全色影像与多光谱影像之间往往存在非线性、不连续的几何错位。如不能进行有效的几何补偿 ［ 10 10 ］ ，融合影像往往会出现晕边现象，导致图像锐度下降，严重影响融合影像的可检测性、可分辨性和可量测性 ［ 11 11 - 12 12 11-12 ］ 。 高分七号是中国首颗民用亚米级高分辨率光学传输型立体测绘卫星，其后视传感器可以获得约0.65 m分辨率的全色影像和约2.6 m分辨率的多光谱影像 ［ 13 13 ］ 。在实验和生产中，发现部分高分七号卫星影像产品存在全色/多光谱影像配准误差，导致了融合产品的晕边现象。影像融合是遥感领域的重要研究方向，目前常用的融合方法包括主成分变换、高通滤波、小波变换等 ［ 14 14 - 17 17 14-17 ］ 。这些方法都是在全色影像和多光谱影像准确配准的条件下使用的，因此实际生产中需要首先解决配准问题。目前，常用的商业软件在对一级影像产品（即通过有理多项式系数（rational polynomial coefficient， RPC）表达地理信息的卫星影像）进行融合误差补偿时，一般采用小面元方法 ［ 18 18 ］ 。但基于小面元校正的融合配准误差补偿模型参数过多，需要依赖均匀分布的匹配点才能稳定求解，不仅限制了运算效率，还难以应对高分七号后视全色和多光谱相机子线阵CCD错位等情况。 本文以高分七号卫星影像为研究对象，基于多光谱影像和全色影像的成像几何规律，分析融合配准误差的产生机理，提出了适用于整幅影像的全色-多光谱同名像点映射误差补偿，以此进行融合配准误差的定量计算与精确补偿，最大限度地提高融合影像产品质量。利用中国广东地区和青海地区的高分七号一级数据验证了融合配准误差检测与补偿的精度，同时提出了分段项复用方法，为高效率、高精度地生产融合影像提供参考。 1 高分七号融合映射误差分析和补偿模型 1 为了设计合理的融合映射误差补偿模型，从成像机理的角度分析高分七号卫星成像误差源及其对融合配准精度的影响。 1.1 全色-多光谱影像的双像几何关系 2 线阵推扫式高分辨率卫星影像的每个扫描行在同一时刻拍摄，物方点与像方点的严密共线关系为 ［ 19 19 ］ ： <math display="block" id="M1"><mi mathvariant="bold-italic">P</mi><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi>S</mi></mrow></msub><mo>+</mo><mi>λ</mi><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mfenced separators="|"><mrow><mi mathvariant="bold-italic">d</mi><mo>+</mo><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mi mathvariant="bold-italic">p</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076936&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076930&type= 42.84133148 4.82600021 式中， <math id="M2"><mi mathvariant="bold-italic">P</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076982&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076993&type= 2.70933342 2.79399991 是物方点在WGS-84坐标系下的坐标向量； <math id="M3"><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi>S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076991&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076989&type= 3.89466691 3.72533321 是扫描行成像时刻的卫星平台中心在WGS-84坐标系下的坐标向量； <math id="M4"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077008&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076992&type= 1.77800000 2.79399991 是尺度参数； <math id="M5"><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi mathvariant="normal">J</mi><mn mathvariant="normal">2000</mn></mrow><mrow><mi mathvariant="normal">W</mi><mi mathvariant="normal">G</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">‑</mi><mn mathvariant="normal">84</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076954&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076953&type= 18.88066673 3.97933316 是成像时刻J2000天球坐标系相对WGS-84地球坐标系的转换矩阵； <math id="M6"><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">d</mi><mi mathvariant="normal">y</mi></mrow><mrow><mi mathvariant="normal">J</mi><mn mathvariant="normal">2000</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076956&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076967&type= 16.42533302 3.97933316 是卫星本体坐标系与J2000坐标系之间的转换矩阵，一般通过定姿设备（星敏感器和陀螺仪）获取； <math id="M7"><mi mathvariant="bold-italic">d</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077032&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077030&type= 2.45533323 2.87866688 是卫星平台中心与相机成像中心之间的偏移线元素； <math id="M8"><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mo>=</mo><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">m</mi></mrow><mrow><mi mathvariant="normal">b</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">d</mi><mi mathvariant="normal">y</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076961&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076975&type= 15.57866764 3.97933316 是卫星相机姿态与卫星本体坐标系之间的转换矩阵，一般由安装角确定，通过实验室标定获取，并通过星上在轨标定优化； <math id="M9"><mi mathvariant="bold-italic">p</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077159&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077157&type= 2.37066650 3.64066648 为物方点 <math id="M10"><mi mathvariant="bold-italic">P</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076982&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076993&type= 2.70933342 2.79399991 对应的像方坐标。 线阵卫星影像融合时需要确定全色影像与多光谱影像的位置映射关系，在产品生产中，一般采用RPC模型进行映射，但为了分析映射误差来源，本文给出了基于严格成像模型的全色-多光谱影像同名像点映射关系： <math display="block" id="M11"><mi mathvariant="bold-italic">P</mi><mo>=</mo><munder><mrow><munder><mrow><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi>S</mi></mrow></msub><mo>+</mo><mi>λ</mi><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mfenced separators="|"><mrow><mi mathvariant="bold-italic">d</mi><mo>+</mo><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mi mathvariant="bold-italic">p</mi></mrow></mfenced></mrow><mo>︸</mo></munder></mrow><mrow><mi mathvariant="normal">全色</mi><mi mathvariant="normal">影像</mi><mi mathvariant="normal">投影</mi><mi mathvariant="normal">关系</mi></mrow></munder><mo>=</mo><munder><mrow><munder><mrow><msubsup><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup><mo>+</mo><mi>λ</mi><mi>'</mi><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup><mfenced separators="|"><mrow><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">'</mi><mo>+</mo><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">'</mi></mrow></mfenced></mrow><mo>︸</mo></munder></mrow><mrow><mi mathvariant="normal">多光</mi><mi mathvariant="normal">谱影</mi><mi mathvariant="normal">像投</mi><mi mathvariant="normal">影关</mi><mi mathvariant="normal">系</mi></mrow></munder></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076986&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076985&type= 86.78334045 9.90600014 式中， <math id="M12"><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi>S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076991&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076989&type= 3.89466691 3.72533321 、 <math id="M13"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077008&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45076992&type= 1.77800000 2.79399991 、 <math id="M14"><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077006&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077001&type= 3.97933316 3.72533321 、 <math id="M15"><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077146&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077130&type= 3.97933316 3.72533321 、 <math id="M16"><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077153&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077151&type= 3.97933316 3.72533321 和 <math id="M17"><mi mathvariant="bold-italic">d</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077032&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077030&type= 2.45533323 2.87866688 为全色影像成像参数； <math id="M18"><mi mathvariant="bold-italic">p</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077159&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077157&type= 2.37066650 3.64066648 为全色影像像点的无畸变坐标； <math id="M19"><msubsup><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077037&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077020&type= 3.97933316 3.97933316 、 <math id="M20"><mi>λ</mi><mi>'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077089&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077079&type= 2.53999996 3.30200005 、 <math id="M21"><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077056&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077042&type= 3.97933316 3.97933316 、 <math id="M22"><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077150&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077148&type= 3.97933316 3.97933316 、 <math id="M23"><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077132&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077155&type= 3.97933316 3.97933316 和 <math id="M24"><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077052&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077050&type= 3.55599999 3.30200005 为多光谱影像成像参数； <math id="M25"><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077077&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077075&type= 3.64066648 4.14866638 为多光谱影像上的同名像点的无畸变坐标。调整 式（2） 式（2） ，获得 <math id="M26"><mi mathvariant="bold-italic">p</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077159&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077157&type= 2.37066650 3.64066648 与 <math id="M27"><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077077&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077075&type= 3.64066648 4.14866638 的函数关系为： <math display="block" id="M28"><mtable columnalign="left"><mtr><mtd><mtext> </mtext><msup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msup><mo>=</mo><munder><mrow><munder><mrow><mfrac><mrow><mi>λ</mi></mrow><mrow><msup><mrow><mi>λ</mi></mrow><mrow><mi>'</mi></mrow></msup></mrow></mfrac><mfenced separators="|"><mrow><msup><mrow><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><msup><mrow><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><msup><mrow><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow></mfenced><mo>·</mo><mfenced separators="|"><mrow><mi mathvariant="bold-italic">p</mi><mo>+</mo><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><mi mathvariant="bold-italic">d</mi></mrow></mfenced></mrow><mo>︸</mo></munder></mrow><mrow><mi mathvariant="normal">第一</mi><mi mathvariant="normal">项</mi></mrow></munder><mo>+</mo></mtd></mtr><mtr><mtd><munder><mrow><munder><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow><mi>λ</mi></mrow><mrow><mi>'</mi></mrow></msup></mrow></mfrac><msup><mrow><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><msup><mrow><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><msup><mrow><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mo>·</mo><mfenced separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi>S</mi></mrow></msub><mo>-</mo><msubsup><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></mrow></mfenced></mrow><mo>︸</mo></munder><mo>-</mo><munder><mrow><munder><mrow><msup><mrow><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mi mathvariant="bold-italic">d</mi><mi mathvariant="bold-italic">'</mi></mrow><mo>︸</mo></munder></mrow><mrow><mi mathvariant="normal">第三</mi><mi mathvariant="normal">项</mi></mrow></munder></mrow><mrow><mi mathvariant="normal">第二</mi><mi mathvariant="normal">项</mi></mrow></munder></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077083&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077067&type= 77.30065918 28.87133408 式中，第二项中的 <math id="M29"><mfenced separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi>S</mi></mrow></msub><mo>-</mo><msubsup><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077087&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077085&type= 15.66333294 3.97933316 是全色影像成像时刻和多光谱影像成像时刻的卫星位移，误差约1×10 －2 m水平； <math id="M30"><mi>λ</mi><mi>'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077089&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077079&type= 2.53999996 3.30200005 是多光谱影像的尺度参数，根据卫星航高，其数量级在1×10 6 水平； <math id="M31"><mrow><mrow><mfenced separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi>S</mi></mrow></msub><mo>-</mo><msubsup><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></mrow></mfenced></mrow><mo>/</mo><mrow><msup><mrow><mi>λ</mi></mrow><mrow><mi>'</mi></mrow></msup></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077099&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077098&type= 19.21933365 4.48733330 的误差水平大约在1×10 －8 m水平，相对于一般为1×10 －6 ~1×10 －5 m水平的卫星相机CCD像元大小， <math id="M32"><mrow><mrow><mfenced separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi>S</mi></mrow></msub><mo>-</mo><msubsup><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></mrow></mfenced></mrow><mo>/</mo><mrow><msup><mrow><mi>λ</mi></mrow><mrow><mi>'</mi></mrow></msup></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077099&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077098&type= 19.21933365 4.48733330 可以忽略不计，因此可以将第二项近似看作常量；第三项均由卫星传感器安置参数组成，也可以视为常量。因此，误差分析主要针对第一项进行。定义 式（3） 式（3） 中第一项的比例因子和旋转因子为： <math display="block" id="M33"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mi>K</mi><mo>=</mo><mrow><mrow><mi>λ</mi></mrow><mo>/</mo><mrow><msup><mrow><mi>λ</mi></mrow><mrow><mi>'</mi></mrow></msup></mrow></mrow></mtd></mtr><mtr><mtd><mi mathvariant="bold-italic">R</mi><mo>=</mo><msup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msup><msubsup><mrow/><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><msup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msup><msubsup><mrow/><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><msup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msup><msubsup><mrow/><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077103&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077101&type= 43.51866913 9.31333351 式中， <math id="M34"><mi>K</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077096&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077094&type= 3.04800010 2.79399991 是全色-多光谱影像成像尺度参数之比，为常数。经测试，高分七号全色-多光谱相机的主光轴夹角约为0.5°，扫描线夹角约为0.02°，由于对同一位置的成像时间差异不到0.5 s，因此旋转因子 <math id="M35"><mi mathvariant="bold-italic">R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077124&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077112&type= 2.87866688 2.79399991 中，地球自转差异极小， <math id="M36"><msup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msup><msubsup><mrow/><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077118&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077117&type= 11.26066685 3.97933316 近似为单位阵，卫星姿态变化极小， <math id="M37"><msup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msup><msubsup><mrow/><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077120&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077105&type= 11.26066685 3.97933316 近似为单位阵。而由于主光轴夹角和扫描线夹角极小，安置矩阵 <math id="M38"><msup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msup><msubsup><mrow/><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077108&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077121&type= 11.26066685 3.97933316 近似为单位阵。 综上所述，旋转因子 <math id="M39"><mi mathvariant="bold-italic">R</mi><mo>≈</mo><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077122&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077110&type= 9.31333351 2.79399991 ，近似为一个单位矩阵，在融合配准误差分析过程中，可以忽略 <math id="M40"><mi mathvariant="bold-italic">R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077124&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077112&type= 2.87866688 2.79399991 本身，仅考虑其误差的影响 。 综上所述， 式（3） 式（3） 可简化为： <math display="block" id="M41"><msup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msup><mo>=</mo><mi>K</mi><mo>·</mo><mi mathvariant="bold-italic">R</mi><mo>·</mo><mfenced separators="|"><mrow><mi mathvariant="bold-italic">p</mi><mo>+</mo><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><mi mathvariant="bold-italic">d</mi></mrow></mfenced><mo>+</mo><mi mathvariant="bold-italic">D</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077128&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077126&type= 42.07933044 4.82600021 式中， <math id="M42"><mi mathvariant="bold-italic">D</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077144&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077143&type= 3.04800010 2.79399991 是 式（3） 式（3） 中第二项和第三项之和，为常量。由此，可以得到理想条件下全色像点与多光谱同名像点的空间关系，但当卫星传感器的内外方位元素误差不可忽视时， 式（3） 式（3） 中的等式则不再成立，需要通过几何补偿，才能实现全色-多光谱影像的精确配准。 1.2 全色-多光谱同名像点映射误差补偿模型 2 根据 式（5） 式（5） ，全色-多光谱同名像点映射误差的来源主要有（见 图1 图1 ）：（1）角元素误差，包括卫星姿态矩阵 <math id="M43"><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077146&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077130&type= 3.97933316 3.72533321 和 <math id="M44"><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077150&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077148&type= 3.97933316 3.97933316 的误差，传感器安置矩阵 <math id="M45"><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077153&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077151&type= 3.97933316 3.72533321 和 <math id="M46"><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077132&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077155&type= 3.97933316 3.97933316 的误差。（2）内定向误差，即无畸变像坐标 <math id="M47"><mi mathvariant="bold-italic">p</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077159&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077157&type= 2.37066650 3.64066648 和 <math id="M48"><msup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077135&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077161&type= 3.13266683 3.89466691 与其对应的有畸变坐标（图像实际量测坐标）的差异。 在卫星影像融合配准中，构建全色-多光谱同名像点映射误差补偿模型的目标是求解其中一个传感器（一般为多光谱传感器）影像像点坐标改正值，从而在姿态误差和内定向误差存在的情况下配平等 式（3） 式（3） ，保证融合解算过程中获得准确的全色-多光谱同名像点。定义多光谱影像的像点坐标改正值为 <math id="M49"><mi mathvariant="bold-italic">δ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077183&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077202&type= 2.20133328 2.79399991 ，定义全色影像和多光谱影像的有畸变坐标 <math id="M50"><msub><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077175&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077173&type= 3.55599999 4.57200003 和 <math id="M51"><msubsup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077179&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077177&type= 3.55599999 4.82600021 ，内定向误差 <math id="M52"><mo>∇</mo><mi mathvariant="bold-italic">p</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077168&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077166&type= 4.74133301 3.64066648 和 <math id="M53"><mo>∇</mo><msup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077172&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077170&type= 5.50333309 3.89466691 ，定义旋转因子的误差中的非颤振项为 <math id="M54"><mo>∇</mo><msub><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077191&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077189&type= 11.00666618 4.57200003 ，颤振项为 <math id="M55"><mo>∇</mo><msub><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077196&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077193&type= 11.00666618 4.57200003 ，根据 式（5） 式（5） ，误差改正值为： <math display="block" id="M56"><mfenced separators="|"><mrow><msubsup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msubsup><mo>+</mo><mo>∇</mo><msup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msup></mrow></mfenced><mo>+</mo><mi mathvariant="bold-italic">δ</mi><mo>=</mo><mi>K</mi><mo>·</mo><mfenced separators="|"><mrow><mi mathvariant="bold-italic">R</mi><mo>+</mo><mo>∇</mo><msub><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><mo>∇</mo><msub><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfenced><mo>·</mo><mfenced separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><mo>∇</mo><mi mathvariant="bold-italic">p</mi><mo>+</mo><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><mi mathvariant="bold-italic">d</mi></mrow></mfenced><mo>+</mo><mi mathvariant="bold-italic">D</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077201&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077199&type= 109.89733887 5.84200001 要实现全色-多光谱同名像点映射误差的补偿，则需要构建 <math id="M57"><mi mathvariant="bold-italic">δ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077183&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077202&type= 2.20133328 2.79399991 的拟合函数模型，在不依赖外部控制数据的情况下求解。 旋转因子的误差项 <math id="M58"><mo>∇</mo><msub><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077187&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077184&type= 11.00666618 4.57200003 包含常量误差、随时间缓慢变化的卫星定姿误差，误差项 <math id="M59"><mo>∇</mo><msub><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077221&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077219&type= 11.00666618 4.57200003 包含随时间周期性变化的颤振误差，因此在拟合过程中，需要将时间作为自变量。但在卫星影像融合产品生产中，一般采用RPC模型进行定向，已无法获取每个扫描行对应的时间戳。针对这一问题，考虑到卫星影像扫描周期的稳定性，采用影像的行号（ l ）替代时间项，作为旋转因子误差的拟合自变量。 卫星影像的内定向误差主要包括内方位元素误差、镜头畸变、子线阵CCD安置错位3种误差源。其中内方位元素误差引起的像点偏移与探元指向角呈线性关系，镜头畸变引起的像点偏移与探元指向角呈多项式关系，而子线阵CCD安置错位引起的像点偏移与探元指向角构成分段线性函数的关系。由于卫星传感器线阵CCD的几何特性，探元指向角与影像列号（ s ）近似呈线性关系，因此卫星影像的内定向误差使用影像的列号 s 作为自变量进行拟合，拟合函数为一个线性函数（内方位元素误差）、一个多项式（镜头畸变）和一个分段线性函数（子线阵CCD安置错位）之和，在模型设计中，将其合并为分段多项式进行拟合，分段个数为子线阵CCD个数。 由于内定向误差一般仅为数个像元大小，为1×10 －6 ~1×10 －5 m水平，因此忽略内定向误差与旋转因子误差的乘积项，考虑到旋转因子 <math id="M60"><mi mathvariant="bold-italic">R</mi><mo>≈</mo><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077208&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077223&type= 9.31333351 2.79399991 ， 式（6） 式（6） 可以转化为： <math display="block" id="M61"><mi mathvariant="bold-italic">δ</mi><mo>=</mo><munder><mrow><munder><mrow><mfenced separators="|"><mrow><mi>K</mi><mo>·</mo><mo>∇</mo><mi mathvariant="bold-italic">p</mi><mo>-</mo><mo>∇</mo><msup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mi mathvariant="bold-italic">'</mi></mrow></msup></mrow></mfenced></mrow><mo>︸</mo></munder></mrow><mrow><mi mathvariant="normal">内定</mi><mi mathvariant="normal">向误</mi><mi mathvariant="normal">差</mi></mrow></munder><mo>+</mo><munder><mrow><munder><mrow><mi>K</mi><mo>·</mo><mfenced separators="|"><mrow><mo>∇</mo><msub><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><mo>∇</mo><msub><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfenced><mo>·</mo><mfenced separators="|"><mrow><mi mathvariant="bold-italic">p</mi><mo>+</mo><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><mi mathvariant="bold-italic">d</mi></mrow></mfenced></mrow><mo>︸</mo></munder></mrow><mrow><mi mathvariant="normal">旋转</mi><mi mathvariant="normal">因子</mi><mi mathvariant="normal">误差</mi></mrow></munder><mo>+</mo><mi mathvariant="normal">常数</mi><mi mathvariant="normal">项</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077227&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077226&type= 105.32532501 10.83733368 内定向误差项 <math id="M62"><mfenced separators="|"><mrow><mi>K</mi><mo>·</mo><mo>∇</mo><mi mathvariant="bold-italic">p</mi><mo>-</mo><mo>∇</mo><msup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077214&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077236&type= 21.59000015 4.82600021 由全色影像内定向误差乘以一个尺度比值，与多光谱影像内定向误差相减得到。但由于高分七号卫星的后视全色传感器和后视多光谱传感器均采用3个等长子线阵CCD构造，且扫描线方向近似平行，全色影像和多光谱影像同名像点的列号近似呈线性关系，因此内定向误差项在拟合中可以合并为一个分段多项式，并采用多光谱影像列号 s 作为唯一的自变量。定义内定向误差 <math id="M63"><mfenced separators="|"><mrow><mi>K</mi><mo>·</mo><mo>∇</mo><mi mathvariant="bold-italic">p</mi><mo>-</mo><mo>∇</mo><msup><mrow><mi mathvariant="bold-italic">p</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077214&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077236&type= 21.59000015 4.82600021 引起的高分七号多光谱影像像点坐标改正值为 <math id="M64"><mfenced separators="|"><mrow><mo>∇</mo><msubsup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>,</mo><mo>∇</mo><msubsup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077241&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077238&type= 15.40933418 4.31799984 ，其拟合方式为： <math display="block" id="M65"><mo>∇</mo><msubsup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>←</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi></mrow></mfenced><mo>=</mo><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>s</mi><mo>+</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mo>+</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup><mo>+</mo><mo>⋯</mo><mo>+</mo><msub><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>s</mi><mo>+</mo><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mo>+</mo><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup><mo>+</mo><mo>⋯</mo><mo>+</mo><msub><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>s</mi><mo>+</mo><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mo>+</mo><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup><mo>+</mo><mo>⋯</mo><mo>+</mo><msub><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077250&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077248&type= 85.34400177 13.88533401 <math display="block" id="M66"><mo>∇</mo><msubsup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>←</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>l</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi></mrow></mfenced><mo>=</mo><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><msub><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>s</mi><mo>+</mo><msub><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msub><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mo>+</mo><msub><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup><mo>+</mo><mo>⋯</mo><mo>+</mo><msub><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>e</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><msub><mrow><mi>e</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>s</mi><mo>+</mo><msub><mrow><mi>e</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msub><mrow><mi>e</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mo>+</mo><msub><mrow><mi>e</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup><mo>+</mo><mo>⋯</mo><mo>+</mo><msub><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>f</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><msub><mrow><mi>f</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>s</mi><mo>+</mo><msub><mrow><mi>f</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msub><mrow><mi>f</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mo>+</mo><msub><mrow><mi>f</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup><mo>+</mo><mo>⋯</mo><mo>+</mo><msub><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077245&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077253&type= 86.61399841 14.64733315 式中， <math id="M67"><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077257&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077263&type= 3.04800010 3.72533321 ， <math id="M68"><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077260&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077258&type= 3.04800010 3.72533321 ，…， <math id="M69"><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077278&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077269&type= 3.04800010 3.72533321 ， <math id="M70"><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077274&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077273&type= 3.04800010 3.72533321 ， <math id="M71"><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077283&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077276&type= 3.04800010 3.72533321 ，…， <math id="M72"><msub><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077296&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077293&type= 3.04800010 3.72533321 ， <math id="M73"><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077289&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077288&type= 2.70933342 3.72533321 ， <math id="M74"><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077308&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077291&type= 2.70933342 3.72533321 ，…， <math id="M75"><msub><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077312&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077310&type= 2.70933342 3.72533321 ， <math id="M76"><msub><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077303&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077301&type= 3.30200005 3.72533321 ， <math id="M77"><msub><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077317&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077316&type= 3.30200005 3.72533321 ，…， <math id="M78"><msub><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077320&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077307&type= 3.30200005 3.72533321 ， <math id="M79"><msub><mrow><mi>e</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077326&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077322&type= 2.70933342 3.72533321 ， <math id="M80"><msub><mrow><mi>e</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077341&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077339&type= 2.70933342 3.72533321 ，…， <math id="M81"><msub><mrow><mi>e</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077345&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077343&type= 2.70933342 3.72533321 ， <math id="M82"><msub><mrow><mi>f</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077334&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077346&type= 2.96333337 4.57200003 ， <math id="M83"><msub><mrow><mi>f</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077350&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077337&type= 2.96333337 4.57200003 ，…， <math id="M84"><msub><mrow><mi>f</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077352&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077355&type= 2.96333337 4.57200003 为分段模型参数； <math id="M85"><mi>s</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077515&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077513&type= 1.26999998 2.79399991 、 <math id="M86"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077714&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077712&type= 1.26999998 2.79399991 分别为多光谱影像列号和行号； <math id="M87"><msub><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077383&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077381&type= 2.45533323 3.72533321 、 <math id="M88"><msub><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077388&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077385&type= 2.45533323 3.72533321 、 <math id="M89"><msub><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077392&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077390&type= 2.45533323 3.72533321 为按列号的分段区间范围； <math id="M90"><msub><mrow><mi>C</mi></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077361&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077360&type= 8.21266651 3.97933316 、 <math id="M91"><msub><mrow><mi>C</mi></mrow><mrow><mi>l</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077398&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077397&type= 8.12800026 3.97933316 分别为影像列方向和行方向的分段模型； <math id="M92"><mo>∇</mo><msubsup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077404&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077401&type= 5.24933338 3.38666677 、 <math id="M93"><mo>∇</mo><msubsup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077409&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077407&type= 5.16466665 3.38666677 分别为列方向和行方向的分段模型补偿值。 由于推扫式成像周期的稳定性，高分七号全色影像与多光谱影像同名像点行号也近似呈线性关系。旋转因子误差的非颤振项在一景影像范围内引起的像点偏移可以通过像方线性模型近似进行拟合。定义误差项 <math id="M94"><mi>K</mi><mo>·</mo><mo>∇</mo><msub><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>·</mo><mfenced separators="|"><mrow><mi mathvariant="bold-italic">p</mi><mo>+</mo><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><mi mathvariant="bold-italic">d</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077416&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077412&type= 35.05199814 4.91066647 与 式（7） 式（7） 中它的常数项引起的像点偏移为 <math id="M95"><mfenced separators="|"><mrow><mo>∇</mo><msubsup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>,</mo><mo>∇</mo><msubsup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077365&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077419&type= 15.40933418 4.31799984 ，其拟合模型为线性模型，以多光谱影像行列号 <math id="M96"><mi>s</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077515&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077513&type= 1.26999998 2.79399991 和 <math id="M97"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077714&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077712&type= 1.26999998 2.79399991 为自变量： <math display="block" id="M98"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mo>∇</mo><msubsup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>←</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>s</mi></mrow></msub><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mi>s</mi><mo>+</mo><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mi>l</mi></mtd></mtr><mtr><mtd><mo>∇</mo><msubsup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>←</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>l</mi></mrow></msub><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub><mo>+</mo><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">5</mn></mrow></msub><mi>s</mi><mo>+</mo><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">6</mn></mrow></msub><mi>l</mi></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077433&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077431&type= 53.42466354 9.65200043 式中， <math id="M99"><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077437&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077435&type= 4.23333359 3.72533321 ， <math id="M100"><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077446&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077439&type= 4.23333359 3.72533321 ，…， <math id="M101"><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">6</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077458&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077449&type= 4.23333359 3.72533321 为线性模型参数； <math id="M102"><msub><mrow><mi>A</mi></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi><mo>,</mo><mi>l</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077454&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077460&type= 11.76866722 3.97933316 、 <math id="M103"><msub><mrow><mi>A</mi></mrow><mrow><mi>l</mi></mrow></msub><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077462&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077457&type= 12.19200039 4.48733330 分别为列方向和行方向线性模型； <math id="M104"><mo>∇</mo><msubsup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077477&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077463&type= 5.24933338 3.38666677 、 <math id="M105"><mo>∇</mo><msubsup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077481&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077466&type= 5.16466665 3.38666677 分别为列方向和行方向线性补偿值。 定义旋转因子的颤振项 <math id="M106"><mi>K</mi><mo>·</mo><mfenced separators="|"><mrow><mi mathvariant="bold-italic">p</mi><mo>+</mo><msubsup><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><mi mathvariant="bold-italic">d</mi></mrow></mfenced><mo>·</mo><mo>∇</mo><msub><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077484&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077482&type= 35.47533035 4.91066647 引起的像点偏移为 <math id="M107"><mfenced separators="|"><mrow><mo>∇</mo><msubsup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>,</mo><mo>∇</mo><msubsup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077485&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077489&type= 15.40933418 4.31799984 ，其拟合模型为以多光谱影像行号 <math id="M108"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077714&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077712&type= 1.26999998 2.79399991 为自变量的三角级数与多光谱影像坐标列号 s 线性项的乘积： <math display="block" id="M109"><mtable columnalign="left"><mtr><mtd><mo>∇</mo><msubsup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>←</mo><msub><mrow><mi>J</mi></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi><mo>,</mo><mi>l</mi></mrow></mfenced><mo>=</mo></mtd></mtr><mtr><mtd><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><msub><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>s</mi></mrow></mfenced><mo>·</mo><msub><mo>∑</mo><mrow><mi>i</mi></mrow></msub><mrow><msub><mrow><mi>T</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mfenced separators="|"><mrow><msub><mrow><mi>p</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>·</mo><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><mo>·</mo><mi>l</mi><mo>+</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfenced></mrow></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077496&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077494&type= 60.53666687 10.83733368 <math display="block" id="M110"><mtable columnalign="left"><mtr><mtd><mo>∇</mo><msubsup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>←</mo><msub><mrow><mi>J</mi></mrow><mrow><mi>l</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi><mo>,</mo><mi>l</mi></mrow></mfenced><mo>=</mo></mtd></mtr><mtr><mtd><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><msub><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mi>s</mi></mrow></mfenced><mo>·</mo><msub><mo>∑</mo><mrow><mi>i</mi></mrow></msub><mrow><msub><mrow><mi>S</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mfenced separators="|"><mrow><msub><mrow><mi>q</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>·</mo><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><mo>·</mo><mi>l</mi><mo>+</mo><msub><mrow><mi>ψ</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfenced></mrow></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077508&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077498&type= 59.60533142 10.75266743 式中， <math id="M111"><msub><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077522&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077518&type= 3.72533321 3.72533321 、 <math id="M112"><msub><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077512&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077510&type= 3.72533321 3.72533321 为影像列号 <math id="M113"><mi>s</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077515&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077513&type= 1.26999998 2.79399991 的系数； <math id="M114"><msub><mrow><mi>T</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077535&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077517&type= 3.30200005 3.72533321 、 <math id="M115"><msub><mrow><mi>S</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077537&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077529&type= 2.87866688 3.72533321 为第 i 个三角函数模型的幅值； <math id="M116"><msub><mrow><mi>p</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077541&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077548&type= 3.04800010 4.57200003 、 <math id="M117"><msub><mrow><mi>q</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077554&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077551&type= 2.62466669 4.57200003 为第 i 个三角函数模型的频率值； <math id="M118"><msub><mrow><mi>φ</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077542&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077556&type= 3.13266683 4.65666676 、 <math id="M119"><msub><mrow><mi>ψ</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077561&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077545&type= 2.87866688 4.48733330 为第 i 个三角函数模型的相位值； <math id="M120"><msub><mrow><mi>J</mi></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi><mo>,</mo><mi>l</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077564&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077546&type= 10.49866676 3.97933316 、 <math id="M121"><msub><mrow><mi>J</mi></mrow><mrow><mi>l</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi><mo>,</mo><mi>l</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077567&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077578&type= 10.41399956 3.97933316 为三角级数补偿模型； <math id="M122"><mo>∇</mo><msubsup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077570&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077580&type= 5.24933338 3.38666677 、 <math id="M123"><mo>∇</mo><msubsup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077585&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077583&type= 5.16466665 3.38666677 为颤振项补偿值。 结合上述误差项，高分七号卫星影像的全色-多光谱同名像点映射误差拟合模型为： <math display="block" id="M124"><mi mathvariant="bold-italic">δ</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mo>∇</mo><msubsup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>+</mo><mo>∇</mo><msubsup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>+</mo><mo>∇</mo><msubsup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup></mtd></mtr><mtr><mtd><mo>∇</mo><msubsup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>+</mo><mo>∇</mo><msubsup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup><mo>+</mo><mo>∇</mo><msubsup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><msup><mrow/><mrow><mi mathvariant="normal">'</mi></mrow></msup></mrow></msubsup></mtd></mtr></mtable></mrow></mfenced><mo>←</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>C</mi></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi></mrow></mfenced><mo>+</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>s</mi></mrow></msub><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mrow><mi>J</mi></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi><mo>,</mo><mi>l</mi></mrow></mfenced></mtd></mtr><mtr><mtd><msub><mrow><mi>C</mi></mrow><mrow><mi>l</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi></mrow></mfenced><mo>+</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>l</mi></mrow></msub><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mrow><mi>J</mi></mrow><mrow><mi>l</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi><mo>,</mo><mi>l</mi></mrow></mfenced></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077574&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077573&type= 82.97332764 10.92199993 在实际应用中，采用四次分段多项式拟合内定向误差的影响，而颤振项的拟合则需要根据不同时期高分七号卫星影像产品的具体误差特性。 1.3 误差补偿模型的求解 2 高质量融合影像的生产需要保证全色-多光谱同名点映射精度达到1个全色像素以内，即0.25多光谱像素，因此需要获得高精度的稠密匹配点进行配准误差补偿模型的解算。定义全色、多光谱影像之间的匹配点为 <math id="M125"><mfenced open="{" close="}" separators="|"><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">p</mi><mo>,</mo><mi mathvariant="bold-italic">p</mi><mi mathvariant="normal">'</mi></mrow></mfenced></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077588&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077587&type= 13.12333298 5.67266655 ，将全色影像上的像点 <math id="M126"><mi mathvariant="bold-italic">p</mi><mo>=</mo><mfenced separators="|"><mrow><mi>s</mi><mo>,</mo><mi>l</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077591&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077576&type= 14.81666565 4.06400013 通过“反投影-重投影”过程（即在航天飞机雷达地形测绘任务（shuttle radar topography mission， SRTM）等高程数据的辅助下利用全色影像RPC参数模型反投影到物方空间，再通过多光谱影像RPC参数模型投影到其多光谱像方空间，见 式（3） 式（3） ）投影到多光谱影像像平面得到像点 <math id="M127"><mi mathvariant="bold-italic">p</mi><mi mathvariant="normal">'</mi><mi mathvariant="normal">'</mi><mo>=</mo><mfenced separators="|"><mrow><mi>s</mi><mi mathvariant="normal">'</mi><mi mathvariant="normal">'</mi><mo>,</mo><mi>l</mi><mi mathvariant="normal">'</mi><mi mathvariant="normal">'</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077596&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077593&type= 18.96533203 4.91066647 ，从而得到映射偏差的观测值： <math display="block" id="M128"><mi mathvariant="bold-italic">δ</mi><mo>=</mo><mi mathvariant="bold-italic">p</mi><mi mathvariant="normal">'</mi><mi mathvariant="normal">'</mi><mo>-</mo><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077599&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077598&type= 18.28800011 4.14866638 在这一过程中，SRTM数据的高程误差影响可以忽略，因为高分七号全色相机和多光谱相机的光线夹角仅为约0.094°，而SRTM的标称精度为16 m，对反投影-重投影过程的影响约为0.01多光谱像素。 根据 式（12） 式（12） 即可得到观测方程，对融合配准模型进行求解。但在模型中，线性项部分和分段多项式部分存在一定的相关性，而三角级数的求解又需要在排除其他误差的干扰后，结合频率域分析等手段进行，因此本文采用分别求解的方式获得配准模型参数的估值，其步骤如下所示。 1） 求解线性项，将 <math id="M129"><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">'</mi><mo>=</mo><mfenced separators="|"><mrow><mi>s</mi><mi>'</mi><mo>,</mo><mi>l</mi><mi>'</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077612&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077610&type= 17.27199936 4.91066647 和 <math id="M130"><mi mathvariant="bold-italic">p</mi><mi>'</mi><mi>'</mi><mo>=</mo><mfenced separators="|"><mrow><mi>s</mi><mi>'</mi><mi>'</mi><mo>,</mo><mi>l</mi><mi>'</mi><mi>'</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077603&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077614&type= 18.96533203 4.91066647 的坐标差异作为误差进行求解： <math display="block" id="M131"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mi>s</mi><mi>'</mi><mi>'</mi><mo>-</mo><msup><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>s</mi></mrow></msub><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mi>l</mi><mi>'</mi><mi>'</mi><mo>-</mo><msup><mrow><mi>l</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>l</mi></mrow></msub><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077606&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077616&type= 28.02466965 9.99066734 2）定义步骤1）中求解得到的线性项模型为 <math id="M132"><msubsup><mrow><mi>A</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>*</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077619&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077618&type= 12.53066730 4.65666676 、 <math id="M133"><msubsup><mrow><mi>A</mi></mrow><mrow><mi>l</mi></mrow><mrow><mi>*</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077626&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077620&type= 12.53066730 4.65666676 ，从误差中减去线性项模型部分，剩余部分用于求解分段项参数： <math display="block" id="M134"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mi>s</mi><mi>'</mi><mi>'</mi><mo>-</mo><msup><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>-</mo><msubsup><mrow><mi>A</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>*</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi></mrow></mfenced></mtd></mtr><mtr><mtd><mi>l</mi><mi>'</mi><mi>'</mi><mo>-</mo><msup><mrow><mi>l</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>-</mo><msubsup><mrow><mi>A</mi></mrow><mrow><mi>l</mi></mrow><mrow><mi>*</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>l</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi></mrow></mfenced></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077638&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077628&type= 41.14799881 10.75266743 3）定义步骤2）中求解的分段项模型为 <math id="M135"><msubsup><mrow><mi>C</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>*</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077642&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077630&type= 9.14400005 4.65666676 、 <math id="M136"><msubsup><mrow><mi>C</mi></mrow><mrow><mi>l</mi></mrow><mrow><mi>*</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077645&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077632&type= 9.14400005 4.65666676 ，再从上述误差中减去分段项模型部分，剩余部分用于进行周期性误差的分析求解。 <math display="block" id="M137"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mi>s</mi><mi>'</mi><mi>'</mi><mo>-</mo><msup><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>-</mo><msubsup><mrow><mi>A</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>*</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo><mo>-</mo><msubsup><mrow><mi>C</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>*</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>J</mi></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi><mo>,</mo><mi>l</mi></mrow></mfenced></mtd></mtr><mtr><mtd><mi>l</mi><mi>'</mi><mi>'</mi><mo>-</mo><msup><mrow><mi>l</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>-</mo><msubsup><mrow><mi>A</mi></mrow><mrow><mi>l</mi></mrow><mrow><mi>*</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo><mo>-</mo><msubsup><mrow><mi>C</mi></mrow><mrow><mi>l</mi></mrow><mrow><mi>*</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>J</mi></mrow><mrow><mi>l</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi><mo>,</mo><mi>l</mi></mrow></mfenced></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077635&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077647&type= 56.81133270 10.75266743 其中，三角级数项的求解方法参考文献［ 20 20 ］。 在映射误差补偿模型（12）中，分段项和颤振项具有较多的参数，需要在同名点均匀覆盖整张影像的情况下进行求解。当影像存在大量水面或云层覆盖时，无法获得均匀的同名点，则无法稳定求解分段项和颤振项。但是，分段项对应的误差源是卫星相机内参数误差、镜头畸变、子线阵CCD安置误差，这一类误差可以长时间保持不变，因此获取时间相近的影像可以复用同一个分段项。在复用分段项的情况下，只需要通过匹配稀疏同名点，求解线性项即可，求解时的观测值为坐标差异与分段项预测值之差： <math display="block" id="M138"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mi>s</mi><mi>'</mi><mi>'</mi><mo>-</mo><msup><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>-</mo><msub><mrow><mover accent="true"><mi>C</mi><mo>̃</mo></mover></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>s</mi></mrow></msub><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mi>l</mi><mi>'</mi><mi>'</mi><mo>-</mo><msup><mrow><mi>l</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>-</mo><msub><mrow><mover accent="true"><mi>C</mi><mo>̃</mo></mover></mrow><mrow><mi>l</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>l</mi></mrow></msub><mo stretchy="false">(</mo><mi>s</mi><mo>,</mo><mi>l</mi><mo stretchy="false">)</mo></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077650&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077648&type= 41.40200043 12.02266598 式中， <math id="M139"><msub><mrow><mover accent="true"><mi>C</mi><mo>̃</mo></mover></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077656&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077654&type= 8.55133343 4.74133301 、 <math id="M140"><msub><mrow><mover accent="true"><mi>C</mi><mo>̃</mo></mover></mrow><mrow><mi>l</mi></mrow></msub><mfenced separators="|"><mrow><mi>s</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077659&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077668&type= 8.46666718 4.74133301 为通过相近时相的其他影像求解得到的分段项。 2 高分七号卫星影像融合配准补偿实验 1 为了验证本文融合配准模型的有效性，采用了中国广东地区6景和青海地区3景高分七号一级影像进行实验分析。广东地区影像拍摄时间为2021年1月、2月和4月，位置分布如 图2 图2 （a）所示，影像文件的前缀为“GF701”。青海地区影像拍摄时间均为2022年1月，位置分布如 图2 图2 （b）所示，影像文件的前缀均为“GF7_DLC”。影像详细信息如 表1 表1 所示。 组别 影像 命名 覆盖地区 应用模式 6 广东组 影像一 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 广东深圳 复用 3 青海组 影像一 GF7_DLC_E101.2_N37.6_20220108_L1A0000677449 青海门源 求解 影像二 GF7_DLC_E101.2_N37.8_20211130_L1A0000639036 青海门源 求解 影像三 GF7_DLC_E101.3_N37.8_20220108_L1A0000677448 青海海北 复用 2.1 融合配准精度及评价指标 2 首先对高分七号全色影像与多光谱影像进行匹配，采用§1.3所述的过程对补偿模型参数进行求解。本文在实验中采用了局部图像特征匹配 （local feature transformer，LoFTR）算法进行匹配，LoFTR算法首先在粗粒度上建立图像特征的检测、描述和匹配，然后在精粒度上细化亚像 素级别的密集匹配 ［ 21 21 ］ ，采用预训练的室外场景匹配模型参数，具有极强的弱纹理影像匹配能力。 在高分七号全色/多光谱影像进行匹配任务中，匹配点密度可以达到每100像素至少一个匹配点，且分布均匀。 融合配准的精度目标一般是达到全色影像的1 像素水平（针对高分七号，相当于0.25多光谱像素），从而消除可见的晕边，而匹配算法往往只能达到0.3~0.5多光谱像素的精度水平。因此实验通过将影像划分为规则格网，对其中匹配点对应的配准误差用鲁棒性递增求解方法进行加权平均，得到更高精度的配准误差值。为了验证格网配准误差的估算精度，本文将多光谱影像进行亚像素平移（如沿行方向和列方向平移0.2~3.8像素）和重采样后，考察平移对格网配准误差估算的影响，发现格网误差的变化（观测值）与实际平移的大小（真值）差异保持在0.1像素以内，这一结果间接证明格网配准误差的量测精度达到了0.1像素水平，能满足融合配准要求。 定义格网中心和配准偏差组成的一组观测值为 <math id="M141"><mfenced separators="|"><mrow><msup><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>,</mo><msup><mrow><mi>l</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>,</mo><mo>∇</mo><mi>s</mi><mo>,</mo><mo>∇</mo><mi>l</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077667&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077678&type= 20.57399940 3.97933316 ，其中 <math id="M142"><msup><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077698&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077680&type= 1.77800000 3.04800010 、 <math id="M143"><msup><mrow><mi>l</mi></mrow><mrow><mi>'</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077701&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077685&type= 1.86266661 3.04800010 为格网中心坐标， <math id="M144"><mo>∇</mo><mi>s</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077689&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077702&type= 4.06400013 3.13266683 、 <math id="M145"><mo>∇</mo><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077705&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077691&type= 4.06400013 3.13266683 为格网中匹配点配准误差平均值。为了评价配准模型中3个子项对融合配准的影响，§2.2统计了对 式（14） 式（14） 和 式（15） 式（15） 对应的线性项和分段项先后进行拟合以后，垂轨方向和沿轨方向的配准残差中误差和 <math id="M146"><msub><mrow><mi>P</mi></mrow><mrow><mo>∇</mo><mi>x</mi><mi>y</mi><mo>≤</mo><mn mathvariant="normal">0.25</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077814&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077812&type= 12.27666664 3.72533321 （格网配准误差小于0.25倍多光谱像素的比例），同时为了定量评价融合配准操作对融合影像质量的影响，还统计了求解两组误差后得到的融合影像的无参考图像质量参数（quality with no reference，QNR）评价指标。QNR指标是影像融合中的一种无参考数据的质量评价指标，它包含两项，一项是光谱偏差参考值，一项是几何偏差参考值，通过计算和整合两种参考值，得到融合质量的客观评价指标。§2.3分析了两组影像映射偏差在 <math id="M147"><mi>l</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077714&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077712&type= 1.26999998 2.79399991 方向上的分布情况，发现广州组影像可能存在颤振误差，对这些影像经过线性项和分段项拟合后的残差进行了三角级数拟合（参考 式（16） 式（16） ），并统计分析。§2.4测试了利用求解出来的分段项校正同一轨道其他影像映射误差的效果（参考 式（17） 式（17） ）。 2.2 基于线性项和分段项的融合配准补偿 2 选取广东地区1~3号影像和青海地区1~2号影像进行融合配准模型求解实验。上述影像的匹配点数均达到1 000万个以上。 表2 表2 给出了两组影像的格网偏移值的初值、线性项补偿残差、分段项补偿残差的统计结果。 图3 图3 和 图4 图4 分别给出了广东组影像1和青海组影像1的补偿残差图。 2 组别 2 影像 2 统计项 2 中误差/像素 2 <math id="M148"><msub><mrow><mi>P</mi></mrow><mrow><mo>∇</mo><mi>x</mi><mi>y</mi><mo>≤</mo><mn>0.25</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077804&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077790&type= 10.58333302 2.87866688 /% 2 QNR x (垂轨) y (沿轨) 9 广东组 3 影像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 3 影像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 影像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 6 青海组 3 影像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 3 影像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 广东组的3景影像在线性项补偿之前，普遍具有列方向约0.60 像素、行方向约0.30 像素的融合配准中误差， <math id="M149"><msub><mrow><mi>P</mi></mrow><mrow><mo>∇</mo><mi>x</mi><mi>y</mi><mo>≤</mo><mn mathvariant="normal">0.25</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077814&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077812&type= 12.27666664 3.72533321 低于1%。经过线性项补偿后，精度提升明显，列方向中误差均低于0.12 像素，行方向中误差均达到约0.20 像素，一定程度上消除了系统性偏移，使得 <math id="M150"><msub><mrow><mi>P</mi></mrow><mrow><mo>∇</mo><mi>x</mi><mi>y</mi><mo>≤</mo><mn mathvariant="normal">0.25</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077814&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077812&type= 12.27666664 3.72533321 大幅提升到约70%水平，但是，结合图 3 3 （b）、 3 3 （e）给出的线性项补偿后残差分布，可以发现仍然存在明显的非线性断裂误差。经过分段项补偿后，有效消除了非线性误差，中误差均达到0.06~0.08像素， <math id="M151"><msub><mrow><mi>P</mi></mrow><mrow><mo>∇</mo><mi>x</mi><mi>y</mi><mo>≤</mo><mn mathvariant="normal">0.25</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077814&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077812&type= 12.27666664 3.72533321 提升到97%以上水平，满足高质量融合的要求。青海组的影像1在通过线性项和分段项进行补偿后，有效消除了非线性误差，列方向和行方向的中误差均达到0.1像素以内， <math id="M152"><msub><mrow><mi>P</mi></mrow><mrow><mo>∇</mo><mi>x</mi><mi>y</mi><mo>≤</mo><mn mathvariant="normal">0.25</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077814&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077812&type= 12.27666664 3.72533321 提升到97%以上，满足高质量融合的要求。而影像2含有一定的云覆盖，校正后 <math id="M153"><msub><mrow><mi>P</mi></mrow><mrow><mo>∇</mo><mi>x</mi><mi>y</mi><mo>≤</mo><mn mathvariant="normal">0.25</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077814&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077812&type= 12.27666664 3.72533321 仅提升到83.8%，这一结果并不能反映实际的精度，要得到准确的评价结果需要先排除含云区域。经检查，去除云雪覆盖率为50%以上的区块后， <math id="M154"><msub><mrow><mi>P</mi></mrow><mrow><mo>∇</mo><mi>x</mi><mi>y</mi><mo>≤</mo><mn mathvariant="normal">0.25</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077814&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077812&type= 12.27666664 3.72533321 达到90%以上。 线性项和分段项补偿前后的融合影像的部分截图如 图5 图5 所示，补偿前的融合影像图 5 5 （a） ~ 5 5 （d）中，在地物颜色突变的位置，往往存在可见的晕边现象，其中图 5 5 （a）、 5 5 （b）为补偿前融合影像建筑物区域显示效果，补偿前影像的建筑物边缘呈现多彩颜色，以致无法准确判断地物边缘； 图 5 5 （c）、 5 5 （d）为在补偿前融合影像水域的显示效果，补偿前，影像的水体边缘不明确，产生了明显的位移。而在补偿后的融合影像5（e）~5（h）中，晕边现象得到显著的缓解，建筑物、水域不再存在模糊的多重边缘，提高了影像可解释性和可量测性。 2.3 基于三角级数项的融合配准补偿 2 实验结果表明，补偿线性项和分段项后，高分七号影像的配准精度已基本满足不出现融合晕边这一要求，但是由 图6 图6 和 图7 图7 中分段项补偿后残差与行方向的误差分布曲线可知，广东组影像的融合配准误差呈现一定的周期性。通过8个三角函数组成的三角级数进行拟合，发现颤振引起的周期性误差振幅均不到0.1 像素，拟合后的配准精度有一定提升，列方向的残差中误差由0.07 像素下降为0.06像素，行方向的残差中误差由0.08 像素下降为0.07 像素。而青海组影像的残差没有沿行方向呈现出周期性。实验结果表明高分七号的颤振误差基本没有对融合质量造成明显影响，可以在生产中忽略不计。 2.4 分段项的复用性验证 2 融合补偿模型中分段项的主要作用是补偿由子线阵CCD安置误差和镜头畸变引起的配准误差，由于CCD安置误差和镜头畸变误差可以在较长时期内保持不变，因此分段项参数理论上在一段时间内具有复用性。为验证补偿模型分段项参数的复用性，将广东组影像一的分段项参数直接对影像4、影像5、影像6进行补偿，将青海组影像1的分段参数对影像3进行补偿。补偿的方法是先用分段项参数进行补偿，然后利用分段项的残差求解线性项，最终得到补偿模型。分段项参数复用模式下的补偿精度如 表3 表3 所示，补偿后的广东组影像中误差基本达到0.07像素水平，青海组影像3中误差达到0.1 像素以内，4组影像的 2 组别 2 影像 2 纠正前后 2 中误差/像素 2 <math id="M155"><msub><mrow><mi>P</mi></mrow><mrow><mo>∇</mo><mi>x</mi><mi>y</mi><mo>≤</mo><mn>0.25</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077804&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077790&type= 10.58333302 2.87866688 / % x (垂轨) y (沿轨) 6 广东组 2 影像4 复用纠正前 0.619 0.335 0.2 复用纠正后 0.073 0.073 98.3 2 影像5 复用纠正前 0.612 0.328 0.4 复用纠正后 0.078 0.072 97.8 2 影像6 复用纠正前 0.604 0.343 0.5 复用纠正后 0.082 0.075 98.3 2 青海组 2 影像7 复用纠正前 0.347 0.288 4.2 复用纠正后 0.044 0.098 98.5 <math id="M156"><msub><mrow><mi>P</mi></mrow><mrow><mo>∇</mo><mi>x</mi><mi>y</mi><mo>≤</mo><mn mathvariant="normal">0.25</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077814&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077812&type= 12.27666664 3.72533321 均提升到98%水平，能满足融合要求。这一实验结果表明分段项参数在一定时间内具有较好的复用性，可以在影像含云、水等不利于匹配的情况下进行复用。 3 结 语 1 本文从高分七号卫星影像的成像机理出发，对影像融合过程中内外方位元素误差和影像内部畸变对配准精度的影响进行分析和建模，提出了由线性模型、四次分段多项式模型、三角级数模型组成的融合映射误差补偿模型。本文对两个地区的高分七号影像进行融合配准模型拟合，实验结果表明：线性模型可以有效消除内外方位元素误差引起的线性像方误差，将垂轨方向的中误差降至0.1像 素左右，沿轨方向中误差降至0.2 像素以下；四次分段多项式模型可明显消除子线阵CCD安置误差引起的像方配准残差不连续性，将垂轨方向中误差降至0.1 像素以下，沿轨方向的中误差从约0.2 像素降至0.1 像素以内。通过线性项和分段项的结合，使得 <math id="M157"><msub><mrow><mi>P</mi></mrow><mrow><mo>∇</mo><mi>x</mi><mi>y</mi><mo>≤</mo><mn mathvariant="normal">0.25</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077814&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077812&type= 12.27666664 3.72533321 提升至约98%水平，基本达到“不出现融合晕边”的要求。而颤振模型的求解和补偿有效抑制了周期性的配准误差，达到采用高分七号一级数据生产融合影像的理论极限精度。 本文测试了分段多项式参数的复用，即采用从其他（同轨或异轨）影像中求解的分段多项式模型参数和本影像中求解的线性项参数进行补偿。实验结果表明，复用模式下融合配准的精度和 <math id="M158"><msub><mrow><mi>P</mi></mrow><mrow><mo>∇</mo><mi>x</mi><mi>y</mi><mo>≤</mo><mn mathvariant="normal">0.25</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077814&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45077812&type= 12.27666664 3.72533321 基本达到求解模式的同一水平，证明了分段多项式主要补偿内部镜头畸变和子线阵CCD安置误差，在不改变零级数据生产相关参数的情况下，基本保持不变。因此，可考虑通过如下流程进行高分七号影像融合生产：（1）从整轨或同一批次获取的所有高分七号影像中，选出含云、水较少且纹理清晰（即沙漠、冰雪等不利于匹配的区域较少）的影像，通过LoFTR等先进的匹配算法进行全影像稠密匹配，然后求解本文提出的融合配准模型；（2）对其他影像，复用已求解的融合配准模型的分段项参数，然后仅匹配少量同名点求解线性项，通过复用模式进行补偿。这一方案可以灵活调整，兼顾精度和效率。 参考文献 Zhang Yongjun ， Zhang Zuxun ， Gong Jianya . Genera‑lized Photogrammetry of Spaceborne， Airborne and Terrestrial Multi-source Remote Sensing Datasets ［J］. Acta Geodaetica et Cartographica Sinica ， 2021 ， 50 （ 1 ）： 1 - 11 . Genera‑lized Photogrammetry of Spaceborne， Airborne and Terrestrial Multi-source Remote Sensing Datasets 张永军 ， 张祖勋 ， 龚健雅 . 天空地多源遥感数据的广义摄影测量学 ［J］. 测绘学报 ， 2021 ， 50 （ 1 ）： 1 - 11 . 天空地多源遥感数据的广义摄影测量学 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 . Technical Framework and Preliminary Practices of Photogrammetric Remote Sensing Intelligent Processing of Multi-source Satellite Images 张永军 ， 万一 ， 史文中 ， 等 . 多源卫星影像的摄影测量遥感智能处理技术框架与初步实践 ［J］. 测绘学报 ， 2021 ， 50 （ 8 ）： 1068 - 1083 . 多源卫星影像的摄影测量遥感智能处理技术框架与初步实践 Weng Yongling ， Tian Qingjiu . Analysis and Evaluation of Method on Remote Sensing Data Fusion ［J］. Remote Sensing Information ， 2003 ， 18 （ 3 ）： 49 - 54 . Analysis and Evaluation of Method on Remote Sensing Data Fusion 翁永玲 ， 田庆久 . 遥感数据融合方法分析与评价综述 ［J］. 遥感信息 ， 2003 ， 18 （ 3 ）： 49 - 54 . 遥感数据融合方法分析与评价综述 Tang Zhiqiang . Analysis and Validation for Rigorous Geometric Model of Space Linear Array CCD Sensors ［D］. Zhengzhou ： Information Engineering University ， 2011 . Analysis and Validation for Rigorous Geometric Model of Space Linear Array CCD Sensors 汤志强 . 航天线阵CCD传感器严格几何模型的分析与验证 ［D］. 郑州 ： 信息工程大学 ， 2011 . 航天线阵CCD传感器严格几何模型的分析与验证 Lei Rong . Study on Theory and Algorithm of the In-flight Geometric Calibration of Spaceborne Linear Array Sensor ［D］. Zhengzhou ： Information Engineering University ， 2011 . Study on Theory and Algorithm of the In-flight Geometric Calibration of Spaceborne Linear Array Sensor 雷蓉 . 星载线阵传感器在轨几何定标的理论与算法研究 ［D］. 郑州 ： 信息工程大学 ， 2011 . 星载线阵传感器在轨几何定标的理论与算法研究 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 . Self-calibration Adjustment of CBERS-02B Long-Strip Imagery 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 . Detection and Compensation of Satelliteflutter Based on Image from Multispectral Camera with Five Spectral Combinations 孙韬 ， 龙辉 ， 赵冬 ， 等 . 基于五谱段合一多光谱相机影像的卫星颤振检测和补偿 ［J］. 光学学报 ， 2014 ， 34 （ 7 ）： 284 - 290 . 基于五谱段合一多光谱相机影像的卫星颤振检测和补偿 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 . Development of Platform Jitter Geometric Analysis and Processing for High-Resolution Optical Satellite Imagery 王密 ， 朱映 ， 范城城 . 高分辨率光学卫星影像平台震颤几何精度影响分析与处理研究综述 ［J］. 武汉大学学报（信息科学版） ， 2018 ， 43 （ 12 ）： 1899 - 1908 . 高分辨率光学卫星影像平台震颤几何精度影响分析与处理研究综述 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 . Block Adjustment of High-Resolution Satellite Images Described by Rational Polynomials Tadono T . Initial Results of Calibration and Validation for PRISM and AVNIR-2 ［J］. Asian Journal of Geoinformatics ， 2006 ， 6 （ 4 ）： 11 - 20 . Initial Results of Calibration and Validation for PRISM and AVNIR-2 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 . Fused Image Quality Assessment Based on Human Visual Characteristics 许丽娜 ， 肖奇 ， 何鲁晓 . 考虑人类视觉特征的融合图像评价方法 ［J］. 武汉大学学报（信息科学版） ， 2019 ， 44 （ 4 ）： 546 - 554 . 考虑人类视觉特征的融合图像评价方法 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 . Sentinel-2 Satellite Image Fusion Method and Quality Evaluation Analysis 梁丽娟 ， 黄万里 ， 张容焱 ， 等 . Sentinel-2卫星影像融合方法与质量评价分析 ［J］. 遥感技术与应用 ， 2019 ， 34 （ 3 ）： 612 - 621 . Sentinel-2卫星影像融合方法与质量评价分析 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 . System Design and Key Technolongies of the GF-7 Satellite 曹海翊 ， 张新伟 ， 赵晨光 ， 等 . 高分七号卫星总体设计与技术创新 ［J］. 中国空间科学技术 ， 2020 ， 40 （ 5 ）： 1 - 9 . 高分七号卫星总体设计与技术创新 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 . Image Fusion with Wavelets Analysis 王智均 ， 李德仁 ， 李清泉 . 利用小波变换对影像进行融合的研究 ［J］. 武汉测绘科技大学学报 ， 2000 ， 25 （ 2 ）： 137 - 142 . 利用小波变换对影像进行融合的研究 Li Pingxiang ， Wang Zhijun . Investigation of Image Fusion Between High Resolution Image and Multi-spectral Image ［J］. Geomatics and Information Scien‑ce of Wuhan University ， 2003 ， 28 （ 1 ）： 14 - 17 . Investigation of Image Fusion Between High Resolution Image and Multi-spectral Image 李平湘 ， 王智均 . 高空间分辨率影像和多光谱影像融合的研究 ［J］. 武汉大学学报（信息科学版） ， 2003 ， 28 （ 1 ）： 14 - 17 . 高空间分辨率影像和多光谱影像融合的研究 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 . A New Fusion Method Based on Resolution Degradation Model for Panchromatic and Multi-spectral Images 李均力 ， 贾永红 . 基于分辨率退化模型的全色和多光谱遥感影像融合方法 ［J］. 武汉大学学报（信息科学版） ， 2004 ， 29 （ 4 ）： 302 - 306 . 基于分辨率退化模型的全色和多光谱遥感影像融合方法 Chen Xueyang ， Yuan Chao . Data Fusion Evaluation of ZY-1 02C Satellite Images ［J］. Geomatics and Spatial Information Technology ， 2013 ， 36 （ 2 ）： 50 - 53 . Data Fusion Evaluation of ZY-1 02C Satellite Images 陈雪洋 ， 袁超 . ZY-102C卫星影像融合方法评价 ［J］. 测绘与空间地理信息 ， 2013 ， 36 （ 2 ）： 50 - 53 . ZY-102C卫星影像融合方法评价 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 . Near Real-Time Automatic Sub-pixel Registration of Panchromatic and Multispectral Images for Pan-Sharpening 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 . China's First Civilian Three-Line-Array Stereo Mapping Satellite： ZY-3 李德仁 . 我国第一颗民用三线阵立体测图卫星：资源三号测绘卫星 ［J］. 测绘学报 ， 2012 ， 41 （ 3 ）： 317 - 322 . 我国第一颗民用三线阵立体测图卫星：资源三号测绘卫星 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 . Essential Technology and Application of Jitter Detection and Compensation for High Resolution Satellites 童小华 ， 叶真 ， 刘世杰 . 高分辨率卫星颤振探测补偿的关键技术方法与应用 ［J］. 测绘学报 ， 2017 ， 46 （ 10 ）： 1500 - 1508 . 高分辨率卫星颤振探测补偿的关键技术方法与应用 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 . LoFTR： Detector-Free Local Feature Matching with Transformers 10.13203/j.whugis20220714 P236 A 1671-8860（2023）07-1029-10 张永军,王梦欣,万一，等.高分七号卫星影像融合中的全色-多光谱配准误差补偿模型[J].武汉大学学报（信息科学版）,2023,48(7):1029-1038.DOI: 10.13203/j.whugis20220714 http://dx.doi.org/10.13203/j.whugis20220714 张永军,王梦欣,万一，等.高分七号卫星影像融合中的全色-多光谱配准误差补偿模型[J].武汉大学学报（信息科学版）,2023,48(7):1029-1038.DOI:10.13203/j.whugis20220714 ZHANG Yongjun,WANG Mengxin,WAN Yi,et al.Compensation Model of GF-7 Panchromatic and Multispectral Image Registration Error[J].Geomatics and Information Science of Wuhan University,2023,48(7):1029-1038.DOI: 10.13203/j.whugis20220714 http://dx.doi.org/10.13203/j.whugis20220714 ZHANG Yongjun,WANG Mengxin,WAN Yi,et al.Compensation Model of GF-7 Panchromatic and Multispectral Image Registration Error[J].Geomatics and Information Science of Wuhan University,2023,48(7):1029-1038.DOI:10.13203/j.whugis20220714 国家自然科学基金（42030102）；湖北省自然科学基金（2020CFA003）；“龙计划”五期合作项目（58817） ；广东省自然资源厅科技项目（GDZRZYKJ2023006）。 2022-11-04 2023-07-05 2023-06-29 39072151 武汉大学学报（信息科学版） 7 48
FULL TEXT LINK
0
Views
0
Downloads
0
CSCD
Alert me when the article has been cited
文章被引用时，请邮件提醒。
Submit
Tools
Download
- H-PDFDownload
- L-PDFDownload
- WORDDownload
- XMLDownload
- Meta-XMLDownload
Export Citation
- BibTeXDownload
- EndnoteDownload
- RefMan Download
- NoteExpressDownload
- NoteFirstDownload
Add to favorites
Add to my album
Publicity Resources
Related Articles
Remote Sensing Image Fusion Based on Low-Level Visual Features and PAPCNN in NSST Domain
An Improved Deep Novel Target Detection Method for Mars Rover Multispectral Imagery
Related Author
No data
Related Institution
Faculty of Geomatics, East China University of Technology
Key Laboratory of Marine Environmental Survey Technology and Application, Ministry of Natural Resources
Key Laboratory of Mine Environmental Monitoring and Improving Around Poyang Lake, Ministry of Natural Resources
State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University
School of Land Science and Technology, China University of Geosciences
0
Bulk Citation
Export to
TOP