CN118365611B 遥感影像内部几何畸变定量化质检方法、装置和电子设备 （航天宏图信息技术股份有限公司）

本发明提供了一种遥感影像内部几何畸变且还采用待检查遥感影像的内部几何畸变中误差和待检查遥感影像的粗差率作为待检查遥感2获取待检查遥感影像和基准参考影像的同名点，其中，所述基准对所述同名点中的所述基准参考影像的地理坐标进行转换，得到同名点中的所述待检查遥感影像的像素坐标计算各所述同名点的像基于所述同名点中的所述待检查遥感影像的像素坐标在所述待检查遥感影像中构建三角网，并根据各所述同名点的像素误差计算所述三角网中每个三角形的每条边的边误基于各所述边误差对所述三角网中的各条边采用不同颜色区根据各所述边误差计算所述待检查遥感影像的内部几何畸变中误差将所述待检查遥感影像的内部几何畸变中误差和所述待检查遥感影像的粗差率作为若所述边误差中的第一边误差小于所述内部畸变误差阈值，若所述边误差中的第二边误差不小于所述内部畸变误差阈值，且第二边误差小于2倍若所述边误差中的第三边误差不小于2倍的所述内部畸变误差阈值，则所述第三边误采用Harris特征点算子对每一个所述格网进行特征点提取，得到在所述基准参考影像中确定与各所述特征点匹配的匹配点，若所述待检查遥感影像为所述原始遥感影像，则采用RPC模型或共线方程模型对所有3若所述待检查遥感影像为所述正射遥感影像，则采用几何多项若所述初始同名点中的第一初始同名点的绝对误差或几何多项式误差大于预设的误若所述待检查遥感影像为所述原始遥感影像，则采用RPC模型对所述同名点中的所述若所述待检查遥感影像为所述正射遥感影像，则采用所述待检查根据各所述理论像素坐标和对应的同名点中的所述待检查遥感影像的像素坐标计算根据同名点像素误差计算算式计算各所述同名点像素坐标中的理论行坐标，colsrc表示与所述理论像素坐标对应的同名点中的所述待检查遥感影像的像素坐标中的列坐标，rowsrc表示与所述理论像素坐标对应的同名点中的所述采用开源库CGAL对所述同名点中的所述待检查遥感影像的像素坐标进行三角网构建，根据边误差计算算式nems-,/T—下T丽—F根据内部几何畸变中误差计算算式计算所述待检查4根据各所述边误差和所述内部畸变误差阈值确定所述边误差中不小于2倍的所述内部根据所述待检查遥感影像的内部几何畸变的质检结果确定所述待检查遥感影像是否基于所述待检查遥感影像的内部几何畸变的质检结果生成所述待检查遥感影像的内差的个数、所述边误差中不小于所述内部畸变误差阈值且小于2倍的所述内部畸变误差阈转换和计算单元，用于对所述同名点中的所述基准参考影像的地理论像素坐标和对应的同名点中的所述待检查遥感影像的像素坐标计算各所述同名点的三角网构建和边误差计算单元，用于基于所述同名点中颜色区分单元，用于基于各所述边误差对所述三角网中的各条边采用不同颜色区分，计算单元，用于根据各所述边误差计算所述待检查遥感影像的内部几何畸变中误差，设定单元，用于将所述待检查遥感影像的内部几何畸变中误差其中，所述颜色区分单元还用于：若所述边误差中的第一边5误差对应的边设置为第二颜色；若所述边误差中的第三边误差不小于2倍的所述内部畸变8.一种电子设备，包括存储器、处理器及存储9.一种计算机可读存储介质，其特征在于，所述6[0003]传统的遥感影像自动化精度质检大多都是通过自动匹配待检查遥感影像和基准参考影像的同名点(即待检查遥感影像中的中的任意一个特征点A与基准参考影像中与上述特征点A匹配的匹配点A’之间的位置坐标差)，继而根据待检查遥感影像中的所有的同名点的绝对几何精度误差自动统计得到整个化精度质检的过程中，只对待检查遥感影像的绝对几基准参考影像之间的几何精度误差)进行了质检，但是并没有对待检查遥感影像的内部几何畸变进行质检，而内部几何畸变(即待检查遥[0005]有鉴于此，本发明的目的在于提供一种遥感影像内部几何畸变定量化质检方法、无法定量计算待检查遥感影像的内部几何畸应的同名点中的所述待检查遥感影像的像素坐标计算各所述同名点的[0009]基于所述同名点中的所述待检查遥感影像的像素坐标在所述待检查遥感影像中7[0012]将所述待检查遥感影像的内部几何畸变中误差和所述待检查遥感影像的粗差率[0017]若所述待检查遥感影像为所述原始遥感影像，则采用RPC模型或共线方程模型对[0019]若所述初始同名点中的第一初始同名点的绝对误差或几何多项式误差大于预设[0021]若所述待检查遥感影像为所述原始遥感影像，则采用RPC模型对所述同名点中的影像的地理坐标在所述待检查遥感影像上的理[0023]根据各所述理论像素坐标和对应的同名点中的所述待检查遥感影像的像素坐标[0024]根据同名点像素误差计算算式计算各所述同8理论像素坐标中的理论行坐标，colsrc表示与所述理论像素坐标对应的同名点中的所述待检查遥感影像的像素坐标中的列坐标，rowsrc表示与所述理论像素坐标对应的同名点中的[0026]采用开源库CGAL对所述同名点中的所述待检查遥感影像的像素坐标进行三角网ms-,/T—T下T丽—TF[0031]若所述边误差中的第三边误差不小于2倍的所述内部畸变误差阈值，则所述第三[0032]进一步的，根据各所述边误差计算所述待检查遥感影像[0033]根据内部几何畸变中误差计算算式计算所述待[0034]根据各所述边误差和所述内部畸变误差阈值确定所述边误差中不小于2倍的所述[0037]根据所述待检查遥感影像的内部几何畸变的质检结果确定所述待检查遥感影像9[0038]基于所述待检查遥感影像的内部几何畸变的质检结果生成所述待检查遥感影像边误差的个数、所述边误差中不小于所述内部畸变误差阈值且小于2倍的所述内部畸变误差阈值的边误差的个数、所述边误差中不小于2倍的所述内部畸变误差阈值的边误差的个[0039]第二方面，本发明实施例还提供了一种遥感影像内部几[0041]转换和计算单元，用于对所述同名点中的所述基准参考影像的地理坐标进行转各所述理论像素坐标和对应的同名点中的所述待检查遥感影像的像素坐标计算各所述同影像的粗差率作为所述待检查遥感影像的内部几何畸变的质可运行指令促使所述处理器运行上述第一方面任差对三角网中的各条边采用不同颜色区分，以表征待检查遥感影像的内部几何畸变的分像的内部几何畸变中误差和待检查遥感影像的粗差率作为待检查遥感影像的内部几何畸何畸变以及无法定量计算待检查遥感影像的内部几何畸变[0050]图1为本发明实施例提供的一种遥感影像内部几何畸变定量化质检方法的流程[0051]图2为本发明实施例提供的一个待检查遥感影像内部几何畸变三角网分色表示[0052]图3为本发明实施例提供的另一个待检查遥感影像内部几何畸变三角网分色表示[0053]图4为本发明实施例提供的一种遥感影像内部几何畸变定量化质检装置的示意术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范的三角网的各条边是用不同颜色区分的，用以表征待检查遥感影像的内部几何畸变的分[0060]根据本发明实施例，提供了一种遥感影像内部几何畸变定量化质检方法的实施[0061]图1是根据本发明实施例的一种遥感影像内部几何畸变定量化质检方法的流程边误差、预设的内部畸变误差阈值和三角网中的边的总数量计算待检查遥感影像的粗差差对三角网中的各条边采用不同颜色区分，以表征待检查遥感影像的内部几何畸变的分像的内部几何畸变中误差和待检查遥感影像的粗差率作为待检查遥感影像的内部几何畸何畸变以及无法定量计算待检查遥感影像的内部几何畸变[0078](2)采用Harris特征点算子对每一个格网进行特征点提取，得到各个格网的特征[0079](3)在基准参考影像中确定与各特征点匹配的匹配点，并将各特征点和与其匹配XYXYZ项式模型)或共线方程模型对所有初始同名点坐标进行最小二乘拟合，得到每个初始同名[0096](5)若待检查遥感影像为正射遥感影像，则采用几何多项式模型对所有初始同名[0097](6)若初始同名点中的第一初始同名点的绝对误差或几何多项式误差大于预设的[0098]在本发明的一个可选实施例中，对同名点中的基准参考影像的地理坐标进行转[0101](2)若待检查遥感影像为正射遥感影像，则采用待检查遥感影像的六参数转换模型(开源库GDAL中提供了六参数转换模型)对同名点中的基准参考影像的地理坐标进行转[0103]根据同名点像素误差计算算式计算各同名点src表示与理论像素坐标对应的同名点中的待检查遥感影像的像素坐标中的列src表示与理论像素坐标对应的同名点中的待检查遥感影像的像素坐标中的行坐[0105](1)采用开源库CGAL对同名点中的待检查遥感影像的像素坐标进行三角网构建，[0109](1)若边误差中的第一边误差小于内部畸变误差阈值，则第一边误差对应的边设[0116](1)根据内部几何畸变中误差计算算式计算待[0117](2)根据各边误差和内部畸变误差阈值确定边误差中不小于2倍的内部畸变误差[0118]上述边误差中不小于2倍的内部畸变误差阈值的边误差的个数即不小于2倍的内[0119](3)根据粗差率计算算式计算待检查遥感影像的粗差率，其中，[0121](1)根据待检查遥感影像的内部几何畸变的质检结果确定待检查遥感影像是否合[0122](2)基于待检查遥感影像的内部几何畸变的质检结果生成待检查遥感影像的内部[0123]图3中示出了内部精度质检报告中的带有颜色区分的三角网的待检查遥感影像的件为待检查遥感影像的粗差率<预设的内部畸变误差阈值比例，并且待检查遥感影像的内中的列误差—B同名点的像素误差中的行误差—B同名点的像素误差中的列误差之间的对[0124]本发明的方法从视觉上可以清晰快速定位到待检查遥感影像的哪些地方存在着内部几何畸变定量化质检装置主要用于执行本发明实施例一中所提供的遥感影像内部几[0127]图4是根据本发明实施例的一种遥感影像内部几何畸变定量化质检装置的示意应的同名点中的待检查遥感影像的像素坐标计算各同名点的差对三角网中的各条边采用不同颜色区分，以表征待检查遥感影像的内部几何畸变的分像的内部几何畸变中误差和待检查遥感影像的粗差率作为待检查遥感影像的内部几何畸何畸变以及无法定量计算待检查遥感影像的内部像中确定与各特征点匹配的匹配点，并将各特征点和与其匹配的匹配点作为初始同名点，的同名点中的待检查遥感影像的像素坐标中的列坐标，rowsrc表示与理论像素坐标对应的计算算式计算三同名点的像素误差中的列误差，RMSCOL2表示构成边的另一个同名点的像素误差中的列误若边误差中的第三边误差不小于2倍的内部畸变误差阈值，则第三边误差对应的边设置为中的边的总数量；根据各边误差和内部畸变误差阈值确定边误差中不小于2倍的内部畸变误差阈值的边误差的个数；根据粗差率计算算式计算待检查遥感影像的[0141]基于待检查遥感影像的内部几何畸变的质检结果生成待检查遥感影像的内部精所述处理器601与所述存储器602之间通过总线通信，所述处理器601执行所述机器可读指方法的各步骤可以通过处理器601中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器601可以是通用处理器，包括中央处理器(CentralProcessingUnit，简称CPU)、网络处理器(NetworkProcessor，简称NP)等；还可以是数字