2021 1-25000 1-50000光学测绘卫星几何检校技术规范

犌犅／犜４０５２３—２０２１

目次

前言…………………………Ⅰ

１范围………………………１

２规范性引用文件…………………………１

３术语和定义………………１

４总则………………………２

４．１检校内容……………２

４．２精度要求……………２

４．３检校流程……………２

５检校条件…………………３

５．１卫星数据获取条件…………………３

５．２检校控制数据条件…………………３

６检校项目和方法…………………………４

６．１检校项目……………４

６．２检校方法……………４

７检校精度验证……………４

７．１检校后影像产品生产………………４

７．２检校精度验证………………………４

８检校结果及检校周期……………………４

８．１检校结果……………４

８．２检校周期……………５

９质量控制和记录要求……………………５

附录Ａ（资料性附录）人工标志布设样例………………６

附录Ｂ（资料性附录）外定向检校方法…………………７

附录Ｃ（资料性附录）基于指向角的内定向几何检校方法……………８

附录Ｄ（资料性附录）基于畸变参数的内定向几何检校方法…………９

附录Ｅ（资料性附录）检校结果示例……………………１１

犌犅／犜４０５２３—２０２１

前言

本标准按照ＧＢ／Ｔ１．１—２００９给出的规则起草。

本标准由中华人民共和国自然资源部提出。

本标准由全国地理信息标准化技术委员会（ＳＡＣ／ＴＣ２３０）归口。

本标准起草单位：自然资源部国土卫星遥感应用中心。

本标准主要起草人：唐新明、祝小勇、张过、范大昭、刘小强、付兴科、王霞、王甄铭、窦显辉、李傲。

Ⅰ

犌犅／犜４０５２３—２０２１

１∶２５０００１∶５００００光学测绘卫星

几何检校技术规范

１范围

本标准规定了１∶２５０００、１∶５００００光学测绘卫星几何检校的总体要求、检校条件、检校项目和方

法、检校精度验证、检校结果及检校周期。

本标准适用于１∶２５０００、１∶５００００光学测绘卫星在轨几何检校。

２规范性引用文件

下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件，仅注日期的版本适用于本文

件。凡是不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。

ＧＢ／Ｔ３５６４２１∶２５０００１∶５００００光学遥感测绘卫星影像产品

ＣＨ／Ｔ２００９全球定位系统实时动态测量（ＲＴＫ）技术规范

ＣＨ／Ｔ３０１９１∶２５０００１∶５００００光学遥感测绘卫星影像产品生产技术规范

３术语和定义

下列术语和定义适用于本文件。

３．１

卫星传感器几何检校狊犪狋犲犾犾犻狋犲狊犲狀狊狅狉狊犵犲狅犿犲狋狉犻犮犮犪犾犻犫狉犪狋犻狅狀

采用卫星地面检校设施或满足几何精度要求的已有地理信息数据（包括数字正射影像和数字高程

模型等），定期对测绘卫星传感器的内外定向参数进行解算，建立传感器严密几何模型的过程。

３．２

卫星传感器内定向参数犻狀狋犲狉犻狅狉狅狉犻犲狀狋犪狋犻狅狀狆犪狉犪犿犲狋犲狉狊狅犳狊犪狋犲犾犾犻狋犲狊犲狀狊狅狉狊

描述卫星影像像点在像方空间坐标系中的位置所需要的参数。

注：卫星传感器内定向参数包括像主点坐标、主距、畸变参数等。

３．３

卫星传感器外定向参数犲狓狋犲狉犻狅狉狅狉犻犲狀狋犪狋犻狅狀狆犪狉犪犿犲狋犲狉狊狅犳狊犪狋犲犾犾犻狋犲狊犲狀狊狅狉狊

描述卫星传感器摄影光束在物方空间坐标系中的位置和方位的基本参数及其随时间的变化量。

注：卫星传感器外定向参数包括卫星轨道测量数据、摄影中心与轨道测量的偏移量、卫星姿态数据、偏置矩阵等。

３．４

偏置矩阵犪犾犻犵狀犿犲狀狋犿犪狋狉犻狓

卫星传感器在像方空间坐标系中的实际观测方位与卫星姿态测量方位之间夹角所构成的正交旋转

矩阵。

３．５

指向角狆狅犻狀狋犻狀犵犳犻犲犾犱犪狀犵犾犲

卫星传感器成像器件的探元在像方空间坐标系中的位置矢量与犡、犢轴的夹角。

３．６

人工标志犪狉狋犻犳犻犮犻犪犾狋犪狉犵犲狋

人工预先布设的、在卫星遥感影像上易于识别的特征目标。

１

犌犅／犜４０５２３—２０２１

注：改写ＧＢ／Ｔ１４９５０—２００９，定义４．１３２。

４总则

４．１检校内容

检校内容包括监测卫星传感器成像几何精度，定期解算探元指向角、偏置矩阵、主点、主距等参数及

其随时间的变化量，减少成像过程中存在的系统误差及其随时间的累积误差，建立传感器严密几何

模型。

４．２精度要求

光学测绘卫星的内定向检校精度应采用指向角角度转换为像元表示，中误差应优于０．３像元；外定

向检校精度应采用卫星姿态角转换为像元表示，中误差应优于１像元。

４．３检校流程

几何检校流程见图１，具体如下：

ａ）检校任务设计。根据卫星传感器几何精度监测数据，确定检校任务，进行检校任务设计，主要

包括明确检校范围，制定成像计划，设计检校和验证方法等。

ｂ）卫星影像获取。根据检校任务设计，针对检校区域制定成像计划，获取卫星影像。

ｃ）检校控制数据获取。利用高分辨率地面参考数据（ＤＯＭ、ＤＥＭ）和待检校卫星影像，采用高精

度影像匹配算法获取密集均匀分布的匹配控制数据；或同时开展人工标志布设工作，对标志影

像提取像点坐标，获取人工标志控制数据。

ｄ）检校参数解算。根据获得的检校控制数据，结合精密定轨数据、精密定姿数据、成像时间数据

等，解算内、外定向参数，建立严密成像模型。

ｅ）检校精度验证。利用地面控制数据，对成像模型进行几何精度验证。

ｆ）检校参数提供。分析参数适用范围和变化趋势，编制精度分析报告，将传感器ＣＣＤ像元指向

角、传感器偏置矩阵等检校参数提供生产使用。

图１几何检校在轨工作流程

２

犌犅／犜４０５２３—２０２１

５检校条件

５．１卫星数据获取条件

５．１．１检校区域位置宜满足以下条件：

ａ）检校区域宜大于被检星载光学传感器幅宽；

ｂ）所在地区每月晴天宜不少于１５ｄ；

ｃ）地形宜选用平地；

ｄ）地物宜具有多样性、变化缓慢，高层建筑少；

ｅ）交通便利。

５．１．２人工标志布设场地宜满足以下条件：

ａ）场地宜大于实际布设面积；场地应平整，相对高差宜小于３０ｃｍ，坡度宜小于５°；标志布设场地

应视线开阔，无建（构）筑物遮挡；

ｂ）气象条件宜选择白天、无雨、风力小于４级、阵风不大于７级。

５．１．３卫星拍摄条件要求如下：

ａ）卫星检校区域无云；

ｂ）检校成像时，卫星侧摆角度应小于１５°。

５．２检校控制数据条件

５．２．１检校控制数据获取方式

检校控制数据可采用人工标志获取或自动匹配获取。

５．２．２自动匹配检校控制数据获取

自动匹配获取检校控制数据要求如下：

ａ）精度要求：用于检校的数字正射影像分辨率应优于０．２ｍ，平面精度应优于０．５ｍ；用于检校的

数字高程模型格网尺寸应优于５ｍ，高程精度应优于１ｍ。

ｂ）控制数据像方坐标确定：通过待检校影像与数字正射影像匹配获得，采用高精度匹配算法和

误匹配点剔除算法，匹配精度应优于０．１像元，误匹配率应小于１０％。

ｃ）控制数据物方坐标确定：通过对数字正射影像和数字高程模型进行内插方式获得。

５．２．３人工标志检校控制数据获取

人工标志获取检校控制数据获取要求如下：

ａ）形状与大小：人工标志的形状与大小应根据实际需要设计，并应能在卫星影像上精确识别和定

位。人工标志形状设计参见附录Ａ中图Ａ．１、图Ａ．２；人工标志的大小设计沿轨向、垂轨向应

大于６像元。

ｂ）材料：人工标志的材料应具备色彩均匀、无光泽、耐高低温、强度高、朗伯特性好、重量轻等特

点，应采用高反射率和低反射率兼顾的特制材料，低反射率应低于５％，高反射率应高于６０％。

ｃ）标志布设：人工标志应在卫星成像拍摄前完成布设，标志数量应根据检校方法复杂度确定，每

片ＣＣＤ对应人工标志数量应大于２个，标志分布应覆盖每片ＣＣＤ全视场。

ｄ）物方三维坐标测定：人工标志坐标测量按ＣＨ／Ｔ２００９执行，并应分别在成像前后进行，测量要

求采用国家大地２０００坐标系，高程采用大地高，平面精度应优于０．１ｍ，高程精度应优于

０．１ｍ，测量后应进行记录归档。

３

犌犅／犜４０５２３—２０２１

６检校项目和方法

６．１检校项目

检校项目包括：

ａ）根据检校控制数据计算光学测绘卫星影像的外定向参数；

ｂ）根据检校控制数据计算光学测绘卫星影像的内定向参数；

ｃ）根据检校控制数据对检校参数精度进行分析和验证。

６．２检校方法

检校方法包括：

ａ）运用空间后方交会、空间前方交会、区域网平差等原理，以共线方程为基础，解求被检传感器的

内、外定向参数；

ｂ）光学测绘卫星检校外定向检校方法参见附录Ｂ，内定向宜采用指向角检校法或畸变参数检校

法，具体参见附录Ｃ、附录Ｄ。

７检校精度验证

７．１检校后影像产品生产

检校后按ＣＨ／Ｔ３０１９进行影像产品生产，产品应满足ＧＢ／Ｔ３５６４２的要求。

７．２检校精度验证

产品精度验证应满足以下要求：

ａ）产品精度利用解算得到的像方坐标和控制数据真实像方坐标按误差计算公式进行计算获得。

ｂ）产品精度主要评价指标包括沿线阵方向精度犿狓（单位为像元）、沿轨道方向精度犿狔（单位为

像元）和像平面精度犿狓狔（单位为像元）。其中犿狓、犿狔、犿狓狔分别按公式（１）～公式（３）计算：

狀

（狓犻－狓犻０）２

犻∑＝１

犿狓＝槡狀…………（１）

狀

（狔犻－狔犻０）２

犻∑＝１

犿狔＝槡狀…………（２）

２２

犿狓狔＝槡犿狓＋犿狔…………（３）

式中：

狓犻，狔犻———由几何模型解算的像方坐标，单位为像元；

狓犻０，狔犻０———由控制点量测的真实像方坐标，单位为像元；

狀———检测点数，不少于２０个，单位为个。

８检校结果及检校周期

８．１检校结果

检校结果应满足以下要求：

４

犌犅／犜４０５２３—２０２１

ａ）检校满足检校精度要求的可提供给数据生产使用；

ｂ）检校结果不满足检校精度要求的应进行分析，并在报告中注明不满足精度的项目，确有必要

时重新检校。

８．２检校周期

光学测绘卫星传感器几何检校周期宜不超过３个月。检校条件容许时，宜适当缩短检校周期。

９质量控制和记录要求

９．１质量控制要求如下：

ａ）几何检校质量控制参照测绘产品质量控制的有关规定执行；

ｂ）质量检查应覆盖几何检校关键环节；

ｃ）检校质量主要通过控制数据质量、沿线阵和轨道方向检校后残差分布进行控制；

ｄ）几何检校作业过程应严格按流程做好记录。

９．２记录主要包括：

ａ）内、外定向输出记录表，记录表样例参见附录Ｅ的Ｅ．１、Ｅ．２；

ｂ）几何检校残差表，记录表样例参见Ｅ．３。

ｃ）人工标志量测记录，如控制点略图、观测手簿、计算结果等；

ｄ）检查报告；

ｅ）技术总结。

５

犌犅／犜４０５２３—２０２１

附录犃

（资料性附录）

人工标志布设样例

图Ａ．１、图Ａ．２给出了人工标志布设样例图。

图犃．１人工标志布设样例一图犃．２人工标志布设样例二

６

犌犅／犜４０５２３—２０２１

附录犅

（资料性附录）

外定向检校方法

犅．１利用检校控制点解求偏置矩阵，补偿载荷安装系统误差。

犅．２外定向的解求方法宜按摄影测量中后方交会方法进行。偏置矩阵可表示为犚犝，见公式（Ｂ．１）。

熿犪１犪２犪３燄熿ｃｏｓφ０－ｓｉｎφ燄熿１００燄熿ｃｏｓκ－ｓｉｎκ０燄

犚犝＝犫１犫２犫３＝犚犚ω犚κ＝０１００ｃｏｓω－ｓｉｎωｓｉｎκｃｏｓκ０

燀犮１犮２犮３燅燀ｓｉｎφ０ｃｏｓφ燅燀０ｓｉｎωｃｏｓω燅燀００１燅

…………（Ｂ．１）

式中：

犚φ犚ω犚κ———传感器在俯仰、横滚、偏航三个方向上的补偿矩阵；

φ———传感器在俯仰方向上的补偿角度；

ω———传感器在横滚方向上的补偿角度；

κ———传感器在偏航方向上的补偿角度。

犅．３卫星传感器的严密几何模型可按公式（Ｂ．２）变换。

熿ｔａｎ（Ψ狔）燄熿犡－犡ｓ燄

１－１－１

ｔａｎ（Ψ狓）＝犿犚犝犚ｂｏｄｙ２ＣＴＳ犢－犢ｓ…………（Ｂ．２）

燀１燅燀犣－犣ｓ燅

式中：

Ψ狓，Ψ狔———传感器探元在卫星本体坐标系下的指向角；

犿———比例尺分母；

犚ｂｏｄｙ２ＣＴＳ———卫星本体坐标系到地球固定地面参考系的旋转矩阵；

犡，犢，犣———地面点坐标；

犡ｓ，犢ｓ，犣ｓ———卫星在地球固定地面参考系下的位置。

犅．４公式（Ｂ．１）可变换为公式（Ｂ．３）。

犪１犡＋犫１犢＋犮１犣

犳狓＝－ｔａｎ（Ψ狓）＝０

犪３犡＋犫３犢＋犮３犣…………（Ｂ．３）

犪２犡＋犫２犢＋犮２犣

犳狔＝－ｔａｎ（Ψ狔）＝０

犪３犡＋犫３犢＋犮３犣

犅．５误差方程见公式（Ｂ．４）。

犞＝犅狓－犾，犠…………（Ｂ．４）

式中：

烄犳狓犳狓犳狓烌

φωκ熿－犳０狓燄

犅＝，犾＝０，狓＝（φ，ω，κ），犠为权矩阵。

犳狔犳狔犳狔燀－犳狔燅

烆φωκ烎

解算狓＝（犅Ｔ犅）－１犅Ｔ犾，可计算出偏置矩阵三个角元素（φ，ω，κ），求解出偏置矩阵。

７

犌犅／犜４０５２３—２０２１

附录犆

（资料性附录）

基于指向角的内定向几何检校方法

犆．１将计算得到的偏置矩阵引入严密成像几何模型，利用控制点进行解算，获取每个探元在卫星本体

坐标系下的指向角。

犆．２主要步骤如下：

ａ）求解各控制点对应的探元指向角

针对传感器第犻个探元，用图像第犻列配准获得的多个控制点求解在卫星本体坐标系下的指

向角，将严密成像几何模型转换为公式（Ｃ．１）、公式（Ｃ．２）：

熿犡－犡ｓ燄熿ｔａｎ（Ψ狓）燄

犢－犢ｓ＝犿犚ｂｏｄｙ２ＣＴＳ犚犝ｔａｎ（Ψ狔）……（Ｃ．１）

燀犣－犣ｓ燅燀１燅ｂｏｄｙ

犪１ｔａｎ（Ψ狓）＋犪２ｔａｎ（Ψ狔）＋犪３犡犡

犳狓＝－＝０

犮１ｔａｎ（Ψ狓）＋犮２ｔａｎ（Ψ狔）＋犮３犣犣………………（Ｃ．２）

犫１ｔａｎ（Ψ狓）＋犫２ｔａｎ（Ψ狔）＋犫３犢犢

犳狔＝－＝０

犮１ｔａｎ（Ψ狓）＋犮２ｔａｎ（Ψ狔）＋犮３犣犣

误差方程公式见（Ｃ．３）：

犞＝犅狓－犾，犠…………（Ｃ．３）

式中：

熿犳狓犳狓燄

０

（ｔａｎ（Ψ狓））（ｔａｎ（Ψ狔））烄－犳狓烌

犅＝，犾＝０，狓＝（ｔａｎ（Ψ狓）ｔａｎ（Ψ狔）），犠为权矩阵。

犳狔犳狔烆－犳狔烎

燀（ｔａｎ（Ψ狓））（ｔａｎ（Ψ狔））燅

利用解算方法狓＝（犅Ｔ犅）－１犅Ｔ犾，可计算出犻探元（ｔａｎ（Ψ狓）ｔａｎ（Ψ狔）），得到该探元指向角。

ｂ）求解各探元指向角

对各配准列犻探元进行探元指向角解求，得到若干探元的探元指向角：（ｔａｎ（Ψ狓）

ｔａｎ（Ψ狔））１，（ｔａｎ（Ψ狓）ｔａｎ（Ψ狔））２，…，（ｔａｎ（Ψ狓）ｔａｎ（Ψ狔））犻。

ｃ）探元指向角处理

对解求获取的所有探元指向角进行多项式拟合，再计算所有探元的探元指向角。

多项式拟合宜取不超过５次多项式，见公式（Ｃ．４）。

ｔａｎ（Ψ狓）＝犪０＋犪１狊＋犪２狊２＋犪３狊３＋犪４狊４＋犪５狊５

２３４５………………（Ｃ．４）

ｔａｎ（Ψ狔）＝犫０＋犫１狊＋犫２狊＋犫３狊＋犫４狊＋犫５狊

式中：

狊———探元编号，对应图像列号。

８

犌犅／犜４０５２３—２０２１

附录犇

（资料性附录）

基于畸变参数的内定向几何检校方法

犇．１将计算得到的偏置矩阵引入严密成像几何模型，利用检校控制点解求相机畸变参数。

犇．２主要步骤如下：

ａ）建立畸变模型

严密成像几何模型见公式（Ｄ．１）：

熿犡燄熿犡ｓ燄熿狓－Δ狓燄

犢＝犢ｓ＋犿犚ｏｒｂｉｔ２ＣＴＳ犚ｂｏｄｙ２ｏｒｂｉｔ犚犝犚ｃａｍｅｒａ２ｂｏｄｙ狔－Δ狔……………（Ｄ．１）

燀犣燅ＷＧＳ８４燀犣ｓ燅燀犳燅

注：因为焦距犳与缩放比例因子存在很强的相关性，故只解求缩放比例因子，不解求犳改正量。

Δ狓＝Δ狓０＋（犽１狉２＋犽２狉４）狓＋狆１（狉２＋２狓２）＋２狆２狓狔＋狔ｓｉｎθ

２４２２

Δ狔＝Δ狔０＋（犽１狉＋犽２狉）狔＋狆２（狉＋２狔）＋２狆１狓狔＋狊狔狔

建立畸变模型见公式（Ｄ．２）～公式（Ｄ．４）（狓为沿轨向，狔为垂轨向）：

犉（狓）＝Δ狓０＋（犽１狉２＋犽２狉４）狓＋狆１（狉２＋２狓２）＋２狆２狓狔＋狔ｓｉｎθ－Δ狓

２４２２………（Ｄ．２）

犉（狔）＝Δ狔０＋（犽１狉＋犽２狉）狔＋狆２（狉＋２狔）＋２狆１狓狔＋狊狔狔－Δ狔

狉＝槡（狓－狓０）２＋（狔－狔０）２…………（Ｄ．３）

狓＝狓－狓０…………（Ｄ．４）

狔＝狔－狔０

式中：

（狓，狔）———像平面坐标；

Δ狓０，Δ狔０———主点偏移量；

狓０，狔０———主点位置；

犽１，犽２———径向畸变系数；

狆１，狆２———偏心畸变系数；

θ———ＣＣＤ在焦平面内的旋转量；

狊狔———缩放比例因子。

ｂ）求解畸变参数

对同一相机不同ＣＣＤ求解不同的主点偏移量（Δ狓０，Δ狔０）、旋转量（θ）和缩放比例因子（狊狔）。

对不同相机求取不同的径向畸变系数（犽１，犽２）、偏心畸变系数（狆１，狆２）。

将严密成像几何模型公式（Ｄ．１）转换为公式（Ｄ．５）：

熿犡－犡ｓ燄熿狓－Δ狓燄

犢－犢ｓ＝犿犚ｏｒｂｉｔ２ＣＴＳ犚ｂｏｄｙ２ｏｒｂｉｔ犚犝犚ｃａｍｅｒａ２ｂｏｄｙ狔－Δ狔………………（Ｄ．５）

燀犣－犣ｓ燅ＣＴＳ燀犳燅

公式（Ｄ．５）可变换为

熿犡－犡ｓ燄熿狓－Δ狓燄

犚ＣＴＳ８４２ｏｒｂｉｔ犚ｏｒｂｉｔ２ｂｏｄｙ犢－犢ｓ＝犿犚犝犚ｃａｍｅｒａ２ｂｏｄｙ狔－Δ狔………………（Ｄ．６）

燀犣－犣ｓ燅ＣＴＳ燀犳燅

熿犪１犪２犪３燄

犚ｃａｍｅｒａ２ｂｏｄｙ＝犫１犫２犫３

燀犮１犮２犮３燅

９

犌犅／犜４０５２３—２０２１

其中，令

熿犡犡燄熿犪１（狓－Δ狓）＋犪２（狔－Δ狔）＋犪３犳燄

犢犢＝犿犚犝犫１（狓－Δ狓）＋犫２（狔－Δ狔）＋犫３犳

燀犣犣燅ｂｏｄｙ燀犮１（狓－Δ狓）＋犮２（狔－Δ狔）＋犮３犳燅

则公式（Ｄ．６）可变换为

熿犡犡燄熿犪１（狓－Δ狓）＋犪２（狔－Δ狔）＋犪３犳燄

犚犝－１犢犢＝犿犫１（狓－Δ狓）＋犫２（狔－Δ狔）＋犫３犳………………（Ｄ．７）

燀犣犣燅ｂｏｄｙ燀犮１（狓－Δ狓）＋犮２（狔－Δ狔）＋犮３犳燅

设

熿犪狌１犪狌２犪狌３燄

犚犝＝犫狌１犫狌２犫狌３

燀犮狌１犮狌２犮狌３燅

则

熿犪狌１犫狌２犮狌３燄

犚犝－１＝犪狌１犫狌２犮狌３

燀犪狌１犫狌２犮狌３燅

设

熿狓狓燄熿犪１（狓－Δ狓）＋犪２（狔－Δ狔）＋犪３犳燄

狔狔＝犫１（狓－Δ狓）＋犫２（狔－Δ狔）＋犫３犳

燀狕狕燅ｂｏｄｙ燀犮１（狓－Δ狓）＋犮２（狔－Δ狔）＋犮３犳燅

则公式（Ｄ．７）可变换为

狓狓＝犪狌１犡犡＋犫狌１犢犢＋犮狌１犣犣

狕狕犪狌３犡犡＋犫狌３犢犢＋犮狌３犣犣…………（Ｄ．８）

狔狔＝犪狌２犡犡＋犫狌２犢犢＋犮狌２犣犣

狕狕犪狌３犡犡＋犫狌３犢犢＋犮狌３犣犣

令

犪狌１犡犡＋犫狌１犢犢＋犮狌１犣犣狓狓

狓＝－

犳犪狌３犡犡＋犫狌３犢犢＋犮狌３犣犣狕狕

犪狌２犡犡＋犫狌２犢犢＋犮狌２犣犣狔狔

狔＝－

犳犪狌３犡犡＋犫狌３犢犢＋犮狌３犣犣狕狕

熿犡燄熿犡犡燄

犢＝犚犝－１犢犢

燀犣燅燀犣犣燅ｂｏｄｙ

列出误差方程

犞＝犅犡－犔，犠

熿（犳狓）（犳狓）（犳狓）（犳狓）（犳狓）（犳狓）（犳狓）（犳狓）（犳狓）（犳狓）（犳狓）燄

犅＝φωκΔ狓０Δ狔０狊狔θ犽１犽２狆１狆２

（犳狔）（犳狔）（犳狔）（犳狔）（犳狔）（犳狔）（犳狔）（犳狔）（犳狔）（犳狔）（犳狔）

燀φωκΔ狓０Δ狔０狊狔θ犽１犽２狆１狆２燅

式中：

犠———权矩阵。

Ｔ－１Ｔ

利用解算原理犡＝（犅犅）犅犔，解出主点偏移量（Δ狓０，Δ狔０）、旋转量（θ）、缩放比例因子（狊狔）、径

向畸变系数（犽１，犽２）和偏心畸变系数（狆１，狆２）。

１０