GB∕T 40523-2021 1：25000 1：50000光学测绘卫星几何检校技术规范.pdf

书 书 书犐 犆犛 犃 ? ? ? ? ? ? ? ? ? ? ?犌犅犜 ?犜 犲 犮 犺 狀 狅 犾 狅 犵 狔狊 狆 犲 犮 犻 犳 犻 犮 犪 狋 犻 狅 狀犳 狅 狉犵 犲 狅犿犲 狋 狉 犻 犮犮 犪 犾 犻 犫 狉 犪 狋 犻 狅 狀狅 犳 狅 狆 狋 犻 犮 犪 犾犿犪 狆 狆 犻 狀 犵狊 犪 狋 犲 犾 犾 犻 狋 犲 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?目次前言范围规范性引用文件术语和定义总则 检校内容 精度要求 检校流程检校条件 卫星数据获取条件 检校控制数据条件检校项目和方法 检校项目 检校方法检校精度验证 检校后影像产品生产 检校精度验证检校结果及检校周期 检校结果 检校周期质量控制和记录要求附录（资料性附录）人工标志布设样例附录（资料性附录）外定向检校方法附录（资料性附录）基于指向角的内定向几何检校方法附录（资料性附录）基于畸变参数的内定向几何检校方法附录（资料性附录）检校结果示例 犌犅犜 前言本标准按照 给出的规则起草。本标准由中华人民共和国自然资源部提出。本标准由全国地理信息标准化技术委员会（ ）归口。本标准起草单位：自然资源部国土卫星遥感应用中心。本标准主要起草人：唐新明、祝小勇、张过、范大昭、刘小强、付兴科、王霞、王甄铭、窦显辉、李傲。犌犅犜 光学测绘卫星几何检校技术规范范围本标准规定了 、 光学测绘卫星几何检校的总体要求、检校条件、检校项目和方法、检校精度验证、检校结果及检校周期。本标准适用于 、 光学测绘卫星在轨几何检校。规范性引用文件下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件，仅注日期的版本适用于本文件。凡是不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。 光学遥感测绘卫星影像产品 全球定位系统实时动态测量（）技术规范 光学遥感测绘卫星影像产品生产技术规范术语和定义下列术语和定义适用于本文件。 卫星传感器几何检校狊 犪 狋 犲 犾 犾 犻 狋 犲狊 犲 狀 狊 狅 狉 狊犵 犲 狅犿犲 狋 狉 犻 犮犮 犪 犾 犻 犫 狉 犪 狋 犻 狅 狀采用卫星地面检校设施或满足几何精度要求的已有地理信息数据（包括数字正射影像和数字高程模型等） ，定期对测绘卫星传感器的内外定向参数进行解算，建立传感器严密几何模型的过程。 卫星传感器内定向参数犻 狀 狋 犲 狉 犻 狅 狉狅 狉 犻 犲 狀 狋 犪 狋 犻 狅 狀狆 犪 狉 犪犿犲 狋 犲 狉 狊狅 犳狊 犪 狋 犲 犾 犾 犻 狋 犲狊 犲 狀 狊 狅 狉 狊描述卫星影像像点在像方空间坐标系中的位置所需要的参数。注：卫星传感器内定向参数包括像主点坐标、主距、畸变参数等。 卫星传感器外定向参数犲 狓 狋 犲 狉 犻 狅 狉狅 狉 犻 犲 狀 狋 犪 狋 犻 狅 狀狆 犪 狉 犪犿犲 狋 犲 狉 狊狅 犳狊 犪 狋 犲 犾 犾 犻 狋 犲狊 犲 狀 狊 狅 狉 狊描述卫星传感器摄影光束在物方空间坐标系中的位置和方位的基本参数及其随时间的变化量。注：卫星传感器外定向参数包括卫星轨道测量数据、摄影中心与轨道测量的偏移量、卫星姿态数据、偏置矩阵等。 偏置矩阵犪 犾 犻 犵 狀犿犲 狀 狋犿犪 狋 狉 犻 狓卫星传感器在像方空间坐标系中的实际观测方位与卫星姿态测量方位之间夹角所构成的正交旋转矩阵。 指向角狆 狅 犻 狀 狋 犻 狀 犵犳 犻 犲 犾 犱犪 狀 犵 犾 犲卫星传感器成像器件的探元在像方空间坐标系中的位置矢量与犡、犢轴的夹角。 人工标志犪 狉 狋 犻 犳 犻 犮 犻 犪 犾狋 犪 狉 犵 犲 狋人工预先布设的、在卫星遥感影像上易于识别的特征目标。犌犅犜 注：改写 ，定义 。总则 检校内容检校内容包括监测卫星传感器成像几何精度，定期解算探元指向角、偏置矩阵、主点、主距等参数及其随时间的变化量，减少成像过程中存在的系统误差及其随时间的累积误差，建立传感器严密几何模型。 精度要求光学测绘卫星的内定向检校精度应采用指向角角度转换为像元表示，中误差应优于 像元；外定向检校精度应采用卫星姿态角转换为像元表示，中误差应优于像元。 检校流程几何检校流程见图，具体如下：）检校任务设计。根据卫星传感器几何精度监测数据，确定检校任务，进行检校任务设计，主要包括明确检校范围，制定成像计划，设计检校和验证方法等。）卫星影像获取。根据检校任务设计，针对检校区域制定成像计划，获取卫星影像。）检校控制数据获取。利用高分辨率地面参考数据（、）和待检校卫星影像，采用高精度影像匹配算法获取密集均匀分布的匹配控制数据；或同时开展人工标志布设工作，对标志影像提取像点坐标，获取人工标志控制数据。）检校参数解算。根据获得的检校控制数据，结合精密定轨数据、精密定姿数据、成像时间数据等，解算内、外定向参数，建立严密成像模型。）检校精度验证。利用地面控制数据，对成像模型进行几何精度验证。）检校参数提供。分析参数适用范围和变化趋势，编制精度分析报告，将传感器像元指向角、传感器偏置矩阵等检校参数提供生产使用。图几何检校在轨工作流程犌犅犜 不合格项检校条件 卫星数据获取条件 检校区域位置宜满足以下条件：）检校区域宜大于被检星载光学传感器幅宽；）所在地区每月晴天宜不少于 ；）地形宜选用平地；）地物宜具有多样性、变化缓慢，高层建筑少；）交通便利。 人工标志布设场地宜满足以下条件：）场地宜大于实际布设面积；场地应平整，相对高差宜小于 ，坡度宜小于 ；标志布设场地应视线开阔，无建（构）筑物遮挡；）气象条件宜选择白天、无雨、风力小于级、阵风不大于级。 卫星拍摄条件要求如下：）卫星检校区域无云；）检校成像时，卫星侧摆角度应小于 。 检校控制数据条件 检校控制数据获取方式检校控制数据可采用人工标志获取或自动匹配获取。 自动匹配检校控制数据获取自动匹配获取检校控制数据要求如下：）精度要求：用于检校的数字正射影像分辨率应优于 ，平面精度应优于 ；用于检校的数字高程模型格网尺寸应优于，高程精度应优于。）控制数据像方坐标确定：通过待检校影像与数字正射影像匹配获得，采用高精度匹配算法和误匹配点剔除算法，匹配精度应优于 像元，误匹配率应小于 。）控制数据物方坐标确定：通过对数字正射影像和数字高程模型进行内插方式获得。 人工标志检校控制数据获取人工标志获取检校控制数据获取要求如下：）形状与大小：人工标志的形状与大小应根据实际需要设计，并应能在卫星影像上精确识别和定位。人工标志形状设计参见附录中图 、图 ；人工标志的大小设计沿轨向、垂轨向应大于像元。）材料：人工标志的材料应具备色彩均匀、无光泽、耐高低温、强度高、朗伯特性好、重量轻等特点，应采用高反射率和低反射率兼顾的特制材料，低反射率应低于，高反射率应高于 。）标志布设：人工标志应在卫星成像拍摄前完成布设，标志数量应根据检校方法复杂度确定，每片对应人工标志数量应大于个，标志分布应覆盖每片全视场。）物方三维坐标测定：人工标志坐标测量按 执行，并应分别在成像前后进行，测量要求采用国家大地 坐标系，高程采用大地高，平面精度应优于 ，高程精度应优于 ，测量后应进行记录归档。犌犅犜 检校项目和方法 检校项目检校项目包括：）根据检校控制数据计算光学测绘卫星影像的外定向参数；）根据检校控制数据计算光学测绘卫星影像的内定向参数；）根据检校控制数据对检校参数精度进行分析和验证。 检校方法检校方法包括：）运用空间后方交会、空间前方交会、区域网平差等原理，以共线方程为基础，解求被检传感器的内、外定向参数；）光学测绘卫星检校外定向检校方法参见附录，内定向宜采用指向角检校法或畸变参数检校法，具体参见附录、附录。检校精度验证 检校后影像产品生产检校后按 进行影像产品生产，产品应满足 的要求。 检校精度验证产品精度验证应满足以下要求：）产品精度利用解算得到的像方坐标和控制数据真实像方坐标按误差计算公式进行计算获得。）产品精度主要评价指标包括沿线阵方向精度犿狓（单位为像元） 、沿轨道方向精度犿狔（单位为像元）和像平面精度犿狓狔（单位为像元） 。其中犿狓、犿狔、犿狓狔分别按公式（）公式（）计算：犿狓狀犻（狓犻狓犻）狀槡（）犿狔狀犻（狔犻狔犻）狀槡（）犿狓狔犿狓犿槡狔（）式中：狓犻，狔犻 由几何模型解算的像方坐标，单位为像元；狓犻，狔犻 由控制点量测的真实像方坐标，单位为像元；狀 检测点数，不少于 个，单位为个。检校结果及检校周期 检校结果检校结果应满足以下要求：犌犅犜 ）检校满足检校精度要求的可提供给数据生产使用；）检校结果不满足检校精度要求的应进行分析，并在报告中注明不满足精度的项目，确有必要时重新检校。 检校周期光学测绘卫星传感器几何检校周期宜不超过个月。检校条件容许时，宜适当缩短检校周期。质量控制和记录要求 质量控制要求如下：）几何检校质量控制参照测绘产品质量控制的有关规定执行；）质量检查应覆盖几何检校关键环节；）检校质量主要通过控制数据质量、沿线阵和轨道方向检校后残差分布进行控制；）几何检校作业过程应严格按流程做好记录。 记录主要包括：）内、外定向输出记录表，记录表样例参见附录的 、 ；）几何检校残差表，记录表样例参见 。）人工标志量测记录，如控制点略图、观测手簿、计算结果等；）检查报告；）技术总结。犌犅犜 附录犃（资料性附录）人工标志布设样例图 、图 给出了人工标志布设样例图。图犃 人工标志布设样例一图犃 人工标志布设样例二犌犅犜 附录犅（资料性附录）外定向检校方法犅 利用检校控制点解求偏置矩阵，补偿载荷安装系统误差。犅 外定向的解求方法宜按摄影测量中后方交会方法进行。偏置矩阵可表示为犚犝，见公式（ ） 。犚犝犪犪犪犫犫犫犮犮犮熿燀燄燅犚犚犚 熿燀燄燅 熿燀燄燅 熿燀燄燅（ ）式中：犚犚犚 传感器在俯仰、横滚、偏航三个方向上的补偿矩阵； 传感器在俯仰方向上的补偿角度； 传感器在横滚方向上的补偿角度； 传感器在偏航方向上的补偿角度。犅 卫星传感器的严密几何模型可按公式（ ）变换。 狔（ ） 狓（ ）熿燀燄燅犿犚犝犚 犡犡犢犢犣犣熿燀燄燅（ ）式中：狓，狔 传感器探元在卫星本体坐标系下的指向角；犿 比例尺分母；犚 卫星本体坐标系到地球固定地面参考系的旋转矩阵；犡，犢，犣 地面点坐标；犡，犢，犣 卫星在地球固定地面参考系下的位置。犅 公式（ ）可变换为公式（ ） 。犳狓犪犡犫犢犮犣犪犡犫犢犮犣 狓（ ）犳狔犪犡犫犢犮犣犪犡犫犢犮犣 狔（ ）（ ）犅 误差方程见公式（ ） 。犞犅狓犾，犠（ ）式中：犅犳狓犳狓犳狓犳狔犳狔犳狔烄烆烌烎，犾犳狓犳狔熿燀燄燅，狓，（），犠为权矩阵。解算狓（犅犅）犅犾，可计算出偏置矩阵三个角元素，（），求解出偏置矩阵。犌犅犜 附录犆（资料性附录）基于指向角的内定向几何检校方法犆 将计算得到的偏置矩阵引入严密成像几何模型，利用控制点进行解算，获取每个探元在卫星本体坐标系下的指向角。犆 主要步骤如下：）求解各控制点对应的探元指向角针对传感器第犻个探元，用图像第犻列配准获得的多个控制点求解在卫星本体坐标系下的指向角，将严密成像几何模型转换为公式（ ） 、公式（ ） ：犡犡犢犢犣犣熿燀燄燅犿犚 犚犝 狓（ ） 狔（ ）熿燀燄燅 （ ）犳狓犪 狓（ ）犪 狔（ ）犪犮 狓（ ）犮 狔（ ）犮犡犡犣犣犳狔犫 狓（ ）犫 狔（ ）犫犮 狓（ ）犮 狔（ ）犮犢犢犣犣（ ）误差方程公式见（ ） ：犞犅狓犾，犠（ ）式中：犅犳狓 狓（ ）（）犳狓 狔（ ）（）犳狔 狓（ ）（）犳狔 狔（ ）（）熿燀燄燅，犾犳狓犳狔烄烆烌烎，狓 狓（ ） 狔（ ）（），犠为权矩阵。利用解算方法狓（犅犅）犅犾，可计算出犻探元 狓（ ） 狔（ ）（），得到该探元指向角。）求解各探元指向角对各配准列犻探元进行探元指向角解求，得到若干探元的探元指向角：（ （狓） （狔） ）， （ （狓） （狔） ）， ， （ （狓） （狔） ）犻。）探元指向角处理对解求获取的所有探元指向角进行多项式拟合，再计算所有探元的探元指向角。多项式拟合宜取不超过次多项式，见公式（ ） 。 狓（ ）犪犪狊犪狊犪狊犪狊犪狊 狔（ ）犫犫狊犫狊犫狊犫狊犫狊（ ）式中：狊 探元编号，对应图像列号。犌犅犜 附录犇（资料性附录）基于畸变参数的内定向几何检校方法犇 将计算得到的偏置矩阵引入严密成像几何模型，利用检校控制点解求相机畸变参数。犇 主要步骤如下：）建立畸变模型严密成像几何模型见公式（ ） ：犡犢犣熿燀燄燅 犡犢犣熿燀燄燅犿犚 犚 犚犝犚 狓狓狔狔犳熿燀燄燅（ ）注：因为焦距犳与缩放比例因子存在很强的相关性，故只解求缩放比例因子，不解求犳改正量。狓狓犽狉犽狉（）狓狆狉狓（）狆狓狔狔 狔狔犽狉犽狉（）狔狆狉狔（）狆狓狔狊狔狔建立畸变模型见公式（ ）公式（ ） （狓为沿轨向，狔为垂轨向） ：犉狓（）狓犽狉犽狉（）狓狆狉狓（）狆狓狔狔 狓犉狔（）狔犽狉犽狉（）狔狆狉狔（）狆狓狔狊狔狔狔（ ）狉狓狓（）狔狔（）槡（ ）狓狓狓狔狔狔（ ）式中：（狓，狔） 像平面坐标；狓，狔 主点偏移量；狓，狔 主点位置；犽，犽 径向畸变系数；狆，狆 偏心畸变系数； 在焦平面内的旋转量；狊狔 缩放比例因子。）求解畸变参数对同一相机不同求解不同的主点偏移量（狓，狔） 、旋转量（）和缩放比例因子（狊狔） 。对不同相机求取不同的径向畸变系数（犽，犽） 、偏心畸变系数（狆，狆） 。将严密成像几何模型公式（ ）转换为公式（ ） ：犡犡犢犢犣犣熿燀燄燅犿犚 犚 犚犝犚 狓狓狔狔犳熿燀燄燅（ ）公式（ ）可变换为犚 犚 犡犡犢犢犣犣熿燀燄燅犿犚犝犚 狓狓狔狔犳熿燀燄燅（ ）犚 犪犪犪犫犫犫犮犮犮熿燀燄燅犌犅犜 其中，令犡犡犢犢犣犣熿燀燄燅 犿犚犝犪狓狓（）犪狔狔（）犪犳犫狓狓（）犫狔狔（）犫犳犮狓狓（）犮狔狔（）犮犳熿燀燄燅则公式（ ）可变换为犚犝犡犡犢犢犣犣熿燀燄燅 犿犪狓狓（）犪狔狔（）犪犳犫狓狓（）犫狔狔（）犫犳犮狓狓（）犮狔狔（）犮犳熿燀燄燅（ ）设犚犝犪狌犪狌犪狌犫狌犫狌犫狌犮狌犮狌犮狌熿燀燄燅则犚犝犪狌犫狌犮狌犪狌犫狌犮狌犪狌犫狌犮狌熿燀燄燅设狓狓狔狔狕 狕熿燀燄燅 犪狓狓（）犪狔狔（）犪犳犫狓狓（）犫狔狔（）犫犳犮狓狓（）犮狔狔（）犮犳熿燀燄燅则公式（ ）可变换为狓狓狕 狕犪狌犡犡犫狌犢犢犮狌犣犣犪狌犡犡犫狌犢犢犮狌犣犣狔狔狕 狕犪狌犡犡犫狌犢犢犮狌犣犣犪狌犡犡犫狌犢犢犮狌犣犣（ ）令犳狓犪狌犡犡犫狌犢犢犮狌犣犣犪狌犡犡犫狌犢犢犮狌犣犣狓狓狕 狕犳狔犪狌犡犡犫狌犢犢犮狌犣犣犪狌犡犡犫狌犢犢犮狌犣犣狔狔狕 狕犡犢犣熿燀燄燅犚犝犡犡犢犢犣犣熿燀燄燅 列出误差方程犞犅犡犔，犠犅犳狓（ ）犳狓（ ）犳狓（ ）犳狓（ ） 狓犳狓（ ） 狔犳狓（ ）狊狔犳狓（ ）犳狓（ ）犽犳狓（ ）犽犳狓（ ）狆犳狓（ ）狆犳狔（ ）犳狔（ ）犳狔（ ）犳狔（ ） 狓犳狔（ ） 狔犳狔（ ）狊狔犳狔（ ）犳狔（ ）犽犳狔（ ）犽犳狔（ ）狆犳狔（ ）狆熿燀燄燅式中：犠 权矩阵。利用解算原理犡（犅犅）犅犔，解出主点偏移量（狓，狔） 、旋转量（） 、缩放比例因子（狊狔） 、径向畸变系数（犽，犽）和偏心畸变系数（狆，狆） 。 犌犅犜 附录犈（资料性附录）检校结果示例犈 外定向检校结果记录示例见表 。表犈 外定向检校参数示例姿态补偿角度检校角度 俯仰 横滚 偏航 犈 内定向检校结果记录示例见表 。表犈 内定向检校参数示例探元序号角度狓 角度狔 犈 几何检校残差表示例见表 。表犈 几何检校残差表示例序号像点位置像素检校残差像素狓狔狓狔 犌犅犜