一种遥感影像沿轨方向局部畸变的快速检测及纠正方法与流程

1.本技术涉及遥感影像检测技术领域，尤其涉及一种遥感影像沿轨方向局部畸变的快速检测及纠正方法。


背景技术：

2.高景一号卫星电荷耦合(charge coupled device，ccd)器件焦面拼接相机采用4片全色多光谱五色tdiccd进行反射镜拼接成一条直线。根据相机内各ccd器件各谱段各探元的精确真实位置、每行影像的成像时间数据、积分时间数据、姿态数据、轨道数据和分景编目信息数据，进行传感器校正，构建每个探元的观测矢量，统一解决几何畸变、几何拼接、成像时间归一化、波段配准等问题，实现每个探元分辨率的统一，从而保证卫星影像的内部几何精度。若上述任何一个因素出现较大误差或错误，均可导致遥感影像几何模型错误，进而影像内部几何精度变差。
3.高景一号卫星在轨运行中，偶发“成像时间数据跳变”情况(发生概率不足1
‰
)，在跳变时刻前后若干时间内成像的数据经传感器校正后会发生影像局部沿轨方向畸变。因各波段分时成像，每个波段发生沿轨方向畸变的时间相同，但地理位置不同，波段合成后，则表现为数十至数百扫描行存在亚像素至数个像素级的波段错位，而其余扫描行正常。
4.现有的对于遥感影像沿轨方向畸变自动检测过程主要为采用sift算法，互信息以及相位一致性等算法获取同名特征点来进行影像匹配，但是采用这些算法计算复杂度较高，运算量较大。由于高景一号四星每天平均订购生产影像数目较多。如果采用现有算法对当日生产的所有影像数据进行计算机自动沿轨方向局部畸变检测，会消耗大量的软硬件资源，且用时长，不能保证数据交付的时效性，进而检测的效率以及时效性较差；另外，由于在一景数据中出现沿轨方向畸变的行数有数十行至数百行，其他扫描行正常，若获取的同名特征点不在出现沿轨方向畸变的行数范围内，则极有可能无法检测出沿轨方向局部畸变问题，即发生漏检。且出现沿轨方向局部畸变的影像，其不同波段畸变的地理位置和畸变的程度均不相同，因此在影像出现沿轨方向畸变的区域进行同名特征点提取极易出错或无法获取同名点，即发生错检或漏检。另外，线阵推扫成像的影像数据存在以下情况：相机各波段存在不同数量不同位置的坏像元，在原始影像上表现为坏线，需进行插值处理。如果相邻2个及2个以上像元均为坏像元，插值处理后会出现地物纹理或边缘扩散。特别是坏线落在反射特性差异较大的地物边缘上(如草地与道路的边缘、林地与建筑物的边缘等)，扩散的边缘波段合成后呈彩带状，易误检测成波段错位。影像辐射几何处理过程中进行的若干次灰度重采样，导致某些角度的纹理或边缘(如与沿轨和垂轨方向夹角均较大)呈锯齿状；水体、沙漠影像纹理较单一，农田、林地影像纹理较细碎，这些也会影响现有自动检测的控制点自动匹配精度，出现错检，进而导致检测准确性较差。


技术实现要素：

5.本技术解决的技术问题是：针对现有技术中检测效率以及准确性较低的问题，本
申请提供了一种遥感影像沿轨方向局部畸变的快速检测及纠正方法，本技术实施例所提供的方案中，一方面通过自动检测方式避免影像数目较大，人工检测所导致的检测效率以及精度不高的问题，另一方面，避免通过sift算法，互信息或相位一致性等算法获取同名特征点来进行影像自动匹配来检测过程中，由于相机各波段可能会存在坏像元、影像辐射几何处理过程中进行的若干次灰度重采样或水体、沙漠影像纹理较单一，农田、林地影像纹理较细碎等原因，影响自动匹配精度，所导致出现错检的问题。
6.第一方面，本技术实施例提供一种遥感影像沿轨方向局部畸变的快速检测及纠正方法，该方法包括：
7.获取待检测影像的成像时间数据文件和元信息文件，从元信息文件中提取待检测影像的信息，其中，所述待检测影像的信息包括所述待检测影像的景号、影像起始绝对行计数和影像结束绝对行计数；
8.根据所述景号从成像时间数据文件中获取影像起始绝对行计数至影像结束绝对行计数部分每一行所记录的积分时间，逐行判断所述每一行所对应的积分时间是否发生跳变；
9.若发生跳变，则生成预警信息、输出所述待检测影像的信息以及出现跳变的行计数和积分时间，其中，所述预警信息用于指示待检测影像会出现沿轨方向局部畸变。
10.可选地，逐行判断所述内容中每一行所对应的积分时间是否发生跳变，包括：
11.逐行判断所述内容中每一行所对应的积分时间是否小于0；或
12.根据影像起始绝对行计数至影像结束绝对行计数部分每一行所对应的积分时间计算第一平均积分时间，逐行判断所述每一行所对应的积分时间是否大于指定倍数的第一平均积分时间。
13.可选地，还包括：根据所述预警信息指示进行人工判读，获取人工判读结果；根据所述人工判读结果以及所述预警信息生成所述待检测影像所对应的沿轨方向局部畸变检测报告。
14.可选地，若人工判断结果为待检测影像出现沿轨方向局部畸变，还包括：
15.根据所述成像时间数据文件中每一行所记录的积分时间计算第二平均积分时间，确定所述成像时间数据文件中出现积分时间跳变的至少一行的每一行所对应的第一行计数、第一积分时间以及第一行时；
16.根据相邻两次出现积分时间跳变行所对应的第一行时、第一积分时间以及第二平均积分时间对所述相邻两次出现积分时间跳变行之间的第一行时进行平滑处理得到纠正后的行时；
17.根据所述纠正后的行时对所述成像时间数据文件中出现积分时间跳变行的每一行所对应的积分时间进行纠正得到纠正后的积分时间。
18.可选地，若相邻两次出现积分时间跳变行为所述成像时间数据文件中第一次和第二次出现积分时间跳变行，根据相邻两次出现积分时间跳变行所对应的第一行时、第一积分时间以及第二平均积分时间对所述相邻两次出现积分时间跳变行之间的第一行时进行平滑处理得到纠正后的行时，包括：
19.通过下式对所述相邻两次出现积分时间跳变行之间的第一行时进行平滑处理得到纠正后的行时：
20.lt
re
(i)＝lt(i)
‑
it(l1) it
mean
21.其中，lt
re
(i)表示纠正后的第i行所对应的行时；lt(i)表示纠正前的第i行所对应的行时；it(l1)表示第l1行所对应的积分时间；it
mean
表示第二平均积分时间；l1 1≤i≤l2，l1表示所述第一次出现积分时间跳变行对应的第一行计数，l2表示第二次出现积分时间跳变行对应的第一行计数。
22.可选地，若相邻两次出现积分时间跳变行为所述成像时间数据文件中第j次和第j 1次出现积分时间跳变行，j为不小于2的正整数，根据相邻两次出现积分时间跳变行所对应的第一行时、第一积分时间以及第二平均积分时间对所述相邻两次出现积分时间跳变行之间的第一行时进行平滑处理得到纠正后的行时，包括：
23.通过下式对所述相邻两次出现积分时间跳变行之间的第一行时进行平滑处理得到纠正后的行时：
24.lt
re
(i)＝lt(i)
‑
(lt(l
j
1)
‑
lt
re
(l
j
)) it
mean
25.其中，l
j
表示第j次出现积分时间跳变的行所对应的第一行计数；lt(l
j
1)表示第j次出现积分时间跳变的下一行所对应的行时；lt
re
(l
j
)表示第j次出现积分时间跳变的行所对应的纠正后的行时；j表示第j次出现积分时间跳变，2≤n；l
j
1≤i≤l
j 1
。
26.可选地，根据所述纠正后的行时对所述成像时间数据文件中出现积分时间跳变行的每一行所对应的积分时间进行纠正得到纠正后的积分时间，包括：
27.通过下式对所述每一行所对应的积分时间进行纠正得到纠正后的积分时间：
28.it
re
(i)＝lt
re
(i 1)
‑
lt
re
(i)
29.其中，it
re
(i)表示第i行纠正后的积分时间；lt
re
(i 1)表示第i 1行纠正后的行时；lt
re
(i)表示第i行纠正后的行时；i＝l
1，
l2，，l
n
，l
1，
l2，，
n
表示所述第一行计数。
30.可选地，还包括：根据所述纠正后的行时和纠正后的积分时间得到新的成像时间数据文件，根据所述新的成像时间数据文件重新生成新的遥感影像。
31.可选地，还包括：对所述新的遥感影像进行几何定位精度检测；若几何定位精度超限，则根据预设参考图像对所述新的遥感影像进行rpc精化处理得到满足预设几何定位精度要求的影像。
32.与现有技术相比，本技术实施例所提供的方案至少具有如下有益效果：
33.1、本技术实施例所提供的方案中，一方面通过自动检测方式避免影像数目较大，人工检测所导致的检测效率以及准确性不高的问题，另一方面，避免通过sift算法，互信息或相位一致性等算法获取同名特征点来进行影像自动匹配来检测过程中，由于相机各波段可能会存在坏像元、影像辐射几何处理过程中进行的若干次灰度重采样或水体、沙漠影像纹理较单一，农田、林地影像纹理较细碎等原因，影响自动匹配精度，所导致出现错检的问题。
34.2、本技术实施例所提供的方案中，通过从元信息文件中提取待检测影像的的景号、影像起始绝对行计数和影像结束绝对行计数，根据景号从成像时间数据文件中获取影像起始绝对行计数至影像结束绝对行计数部分每一行所记录的积分时间，逐行判断每一行所对应的积分时间是否发生跳变，根据积分时间跳变来确定影像可能会存在沿轨方向局部畸变，避免由于水体、沙漠影像纹理较单一，农田、林地影像纹理较细碎，影响自动匹配的精度，易漏检，导致检测准确性不高。
附图说明
35.图1为本技术实施例所提供的一种遥感影像沿轨方向局部畸变的快速检测及纠正方法的流程示意图；
36.图2为本技术实施例所提供的一种元信息文件的结构示意图；
37.图3为本技术实施例所提供的一种成像时间数据文件(.it文件)的结构示意图；
38.图4为本技术实施例所提供的一种积分时间随行计数变化的示意图；
39.图5a为本技术实施例所提供的一种纠正前影像的示意图；
40.图5b为本技术实施例所提供的一种纠正后影像的示意图；
41.图6a为本技术实施例所提供的又一种纠正前影像的示意图；
42.图6b为本技术实施例所提供的又一种纠正后影像的示意图。
具体实施方式
43.本技术实施例提供的方案中，所描述的实施例仅是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本技术保护的范围。
44.以下结合说明书附图对本技术实施例所提供的一种遥感影像沿轨方向局部畸变的快速检测及纠正方法做进一步详细的说明，该方法具体实现方式可以包括以下步骤(方法流程如图1所示)：
45.步骤101，获取待检测影像的成像时间数据文件和元信息文件，从元信息文件中提取待检测影像的信息，其中，所述待检测影像的信息包括所述待检测影像的景号、影像起始绝对行计数和影像结束绝对行计数。
46.具体的，在本技术实施例所提供的方案中，卫星在生成遥感影像所对应的0级数据时，同时会生成并保存对应的对应辅助文件，其中，辅助文件包括时间数据文件(.it文件)以及元信息文件(.xml文件)等。在元信息文件(.xml文件)中记录了多光谱和全色影像所对应的起始绝对行计数和结束绝对行计数，其中，绝对行计数是指待检测影像所对应的数据在整个条带所对应的数据中的行计数值。由于本技术主要针对的是多光谱波段间沿轨方向的畸变检测，故从元信息文件中只需提取多光谱影像所对应数据的起始绝对行计数和结束绝对行计数。
47.为了便于理解上述从元信息文件中提取待检测影像所对应数据的起始绝对行计数和结束绝对行计数过程，下面对其进行简要介绍。
48.参见图2，为本技术实施例提供的一种元信息文件(.xml文件)的结构示意图。在图2中，分别读取待检测影像所对应的景号、卫星标识、接收轨道号、成像轨道号、条带号、成像模式以及成像数据起始绝对行计数(l1)和结束绝对行计数(l2)等，即分别，即分别读取<sceneid>字段对应的景号(2904594)、<satelliteid>字段对应的卫星标识(gj1b)、<receive station id>字段对应的接收轨道号(gua)、<orbitid>字段对应的成像轨道号(18501)、<porbitid>字段对应的条带号(18493)、<datasetid>字段对应的成像模式(168597)、<scenestartlin>字段对应的起始绝对行计数(692339,2770471)、<scenestopline>字段对应的结束绝对行计数(699138，2796870)。
49.步骤102，根据所述景号从成像时间数据文件中获取影像起始绝对行计数至影像
结束绝对行计数部分每一行所记录的积分时间，逐行判断所述每一行所对应的积分时间是否发生跳变。
50.具体的，计算机设备在获取待检测影像的景号以及成像起始绝对行计数和结束绝对行计数之后，根据景号打开待检测影像所在条带所对应的成像时间数据文件(.it文件)，然后从it文件中截取待检测影像的起始绝对行计数(l
star
)至结束绝对行计数(l
stop
)之间的内容，然后逐行判断l
star
到l
stop
每一行所对应的积分时间是否发生跳变。在本技术实施例所提供的方案中，判断每一行所对应的积分时间是否发生跳变的方式有多种，下面以其中一种为例进行说明。
51.在一种可能实现的方式中，逐行判断所述内容中每一行所对应的积分时间是否发生跳变，包括：逐行判断所述内容中每一行所对应的积分时间是否小于0；或根据影像起始绝对行计数至影像结束绝对行计数部分每一行所对应的积分时间计算第一平均积分时间，逐行判断所述每一行所对应的积分时间是否大于指定倍数的第一平均积分时间。
52.具体的，在it文件中第一列记录的是绝对行计数，第三列记录的是积分时间，逐行判断l
st
到l
stop
每一行第三列所对应的积分时间是否小于0或积分时间是否大于指定倍数的第一平均积分时间，例如，指定倍数为2倍。
53.参见图3所示，为本技术实施例所提供的一种成像时间数据文件(.it文件)的结构示意图。在图3中，行计数为693151所对应的积分时间为0.001742005，而其相邻上一行的行计数为693150，其所对应的积分时间为0.000346005，相邻下一行的行计数为693152，其所对应的积分时间为0.000344992，故行计数为693151的行积分时间与相邻上一行和相邻下一行存在明显的跳变。参见图4所示，本技术实施例提供的一种积分时间随行计数变化的示意图。在图4中可以明显看出部分行存在积分时间跳变的异常值。
54.步骤103，若发生跳变，则生成预警信息、输出所述待检测影像的信息以及出现跳变的行计数和积分时间，其中，所述预警信息用于指示待检测影像会出现沿轨方向局部畸变，其中，待检测影像的信息包括卫星标识、接收轨道号、成像轨道号、条带号、景号、景起始绝对行计数、景结束绝对行计数、异常行计数及积分时间异常值。
55.具体的，若l
star
到l
stop
行中任一行存在积分时间跳变，则生成预警信息、输出所述待检测影像的信息以及出现跳变的行计数和积分时间；否则，过程结束。
56.进一步，为了提高检测的准确性，在检测出待检测影像出现沿轨方向畸变时，需指示对待检测影像进行人工判读，通过人工来进一步确定待检测影像是否真的存在沿轨方向畸变。
57.在一种可能实现的方式中，还包括：根据所述预警信息指示进行人工判读，获取人工判读结果；根据所述人工判读结果以及所述预警信息生成所述待检测影像所对应的沿轨方向局部畸变检测报告。
58.具体的，技术人员根据预警信息提取对应产品tiff图像，根据每景预警影像的异常行计数在该景的位置，快速定位影像发生沿轨方向局部畸变的位置，得到人工判读结果，根据所述人工判读结果以及所述预警信息生成所述待检测影像所对应的沿轨方向局部畸变检测报告。
59.进一步，为了正常生产遥感影像，需对检测出存在沿轨方向畸变的影像数据进行纠正，在一种可能实现的方式中，若人工判断结果为待检测影像出现沿轨方向局部畸变，还
包括：根据所述成像时间数据文件中每一行所记录的积分时间计算第二平均积分时间，确定所述成像时间数据文件中出现积分时间跳变的至少一行的每一行所对应的第一行计数、第一积分时间以及第一行时；根据相邻两次出现积分时间跳变行所对应的第一行时、第一积分时间以及第二平均积分时间对所述相邻两次出现积分时间跳变行之间的第一行时进行平滑处理得到纠正后的行时；根据所述纠正后的行时对所述成像时间数据文件中出现积分时间跳变行的每一行所对应的积分时间进行纠正得到纠正后的积分时间。
60.具体的，在本技术实施例所提供的方案中，根据it文件中l
star
行到l
stop
行所对应的内容中每一行所记录的积分时间计算第二平均积分时间，然后确定l
star
行到l
stop
行存在积分时间跳变的至少一行，以及存在积分时间跳变的至少一行所对应第一行计数(l1,l2，，l
n
)、第一积分时间(it(l1),it(l2),
…
,it(l
n
))以及第一行时(lt(l1),lt(l2),
…
,lt(l
n
))，其中，n为不小于1的正整数。根据相邻两次出现积分时间跳变行所对应的第一行时、第一积分时间以及第二平均积分时间对所述相邻两次出现积分时间跳变行之间的第一行时进行平滑处理得到纠正后的行时；根据所述纠正后的行时对所述成像时间数据文件中出现积分时间跳变行的每一行所对应的积分时间进行纠正得到纠正后的积分时间。在it文件中l
star
行到l
stop
行所对应的内容中，存在积分时间跳变的相邻两行，可以是行计数连续的两行，也可以是行计数不连续的两行，例如，相邻两次出现积分时间跳变行为693151和693152，或相邻两次出现积分时间跳变行为693151和693526。
61.进一步，若相邻两次出现积分时间跳变行为所述成像时间数据文件中第一次和第二次出现积分时间跳变的行，其中，第一次出现积分时间跳变的行是指按照l
star
行到l
stop
的顺序第一次出现积分时间跳变的行，则通过下式对第一次和第二次出现积分时间跳变行之间的第一行时进行平滑处理得到纠正后的行时：
62.lt
re
(i)＝lt(i)
‑
it(l1) it
mean
63.其中，lt
re
(i)表示纠正后的第i行所对应的行时；lt(i)表示纠正前的第i行所对应的行时；it(l1)表示第l1行所对应的积分时间；it
mean
表示第二平均积分时间；l1 1≤i≤l2，l1表示所述第一次出现积分时间跳变行对应的第一行计数，l2表示第二次出现积分时间跳变行对应的第一行计数。
64.进一步，若相邻两次出现积分时间跳变行为所述成像时间数据文件中第j次和第j 1次出现积分时间跳变行，j为不小于2的正整数，例如，第二次和第三次出现积分时间跳变行，或者第四次和第五次出现积分时间跳变行等。
65.若相邻两次出现积分时间跳变行为所述成像时间数据文件中第j次和第j 1次出现积分时间跳变行，j为不小于2的正整数，根据相邻两次出现积分时间跳变行所对应的第一行时、第一积分时间以及第二平均积分时间对所述相邻两次出现积分时间跳变行之间的第一行时进行平滑处理得到纠正后的行时，包括：
66.通过下式对所述相邻两次出现积分时间跳变行之间的第一行时进行平滑处理得到纠正后的行时：
67.lt
re
(i)＝lt(i)
‑
(lt(l
j
1)
‑
lt
re
(l
j
)) it
mean
68.其中，l
j
表示第j次出现积分时间跳变的行所对应的第一行计数；lt(l
j
1)表示第j次出现积分时间跳变的下一行所对应的行时；lt
re
(l
j
)表示第j次出现积分时间跳变的行所对应的纠正后的行时；j表示第j次出现积分时间跳变，2≤j≤n；l
j
1≤i≤l
j 1
。
69.进一步，在对it文件中积分时间进行纠正后，还需要对行时进行纠正。在一种可能实现的方式中，根据所述纠正后的行时对所述成像时间数据文件中出现积分时间跳变行的每一行所对应的积分时间进行纠正得到纠正后的积分时间，包括：
70.通过下式对所述每一行所对应的积分时间进行纠正得到纠正后的积分时间：
71.it
re
(i)＝lt
re
(i 1)
‑
lt
re
(i)
72.其中，it
re
(i)表示第i行纠正后的积分时间；lt
re
(i 1)表示第i 1行纠正后的行时；lt
re
(i)表示第i行纠正后的行时；i＝l
1，
l2,,l
n
，l
1,
l2，，
n
表示所述第一行计数。
73.进一步，在一种可能实现的方式中，还包括：根据所述纠正后的行时和纠正后的积分时间得到新的成像时间数据文件，根据所述新的成像时间数据文件重新生成新的遥感影像。
74.具体的，将新的成像时间数据文件替换生产系统原it文件，重新订购该条带影像，生产得到无畸变的新影像。
75.进一步，在一种可能实现的方式中，还包括：对所述新的遥感影像进行几何定位精度检测；若几何定位精度超限，则根据预设参考图像对所述新的遥感影像进行理性多项式系数(rational polynomial coefficients，rpc)精化处理得到满足预设几何定位精度要求的影像。
76.进一步，为了验证本技术所提供的遥感影像沿轨方向局部畸变的检测效果，下面以高景一号卫星为例对本技术实施例所提供方案的效果进行说明。
77.例如，待检验影像数据为2020年11月份高景一号四星所有产品影像数据。
78.1)验证方法
79.(1)采用本技术实施例所提供的技术、现有技术1、现有技术2、现有技术3及人工目视判读方法分别对上述验证数据进行沿轨方向局部畸变，对比分析其准确率和时效性，其中，现有技术1是指基于sifi特征遥感影像自动配准与拼接，现有技术2是指基于互信息的高性能遥感影像配准，现有技术3是指基于相位一致性的异源图像匹配。
80.(2)采用本技术实施例所提供的技术对检测出沿轨方向存在局部畸变的影像进行纠正，并与纠正前影像进行对比验证其纠正效果。
81.2)验证结果
82.(1)采用本技术实施例所提供的技术、现有技术1、现有技术2、现有技术3及人工目视判读方法对高景一号四星2020年11月份影像产品进行沿轨方向存在局部畸变检测的结果如下：表1表示基于本技术实施例所提供的技术的检测方法检测高景一号四星11月份影像产品沿轨方向存在局部畸变的输出结果；表2表示基于本技术实施例所提供的技术以及现有技术的检测方法检测高景一号四星2020年11月份影像产品沿轨方向存在局部畸变检测结果对比。
83.表1
[0084][0085]
表2
[0086]
[0087][0088]
(2)采用本技术实施例所提供的技术对上述检测出沿轨方向存在局部畸变的影响进行纠正，并根据参考影像对其进行几何定位精度检测，若影像存在精度超限情况，则对其进行rpc精化处理。以景号为2900681和2900694的影像为例，参见表3为景号为2900681影像的信息，参见表4为景号2900694影像的信息，其纠正前后对比结果如图5a和图5b、图6a和图6b所示。
[0089]
表3
[0090][0091]
表4
[0092][0093]
3)验证结果分析与结论依据上述验证结果，结论如下：
[0094]
(1)高景一号卫星地面数据处理系统每天生产约300景影像产品，所有产品如果均要通过人工目视检测沿轨方向局部畸变，每景约需要1分钟，300景共需要300分钟；若使用
现有技术1、2、3检测沿轨方向局部畸变，由于运算量均较大，每景分别需要约2分钟、1.8分钟和1.7分钟，300景总用时分别为600分钟、540分钟以及510分钟；而使用基于本技术技术的影像沿轨方向局部畸变检测方法每景约需要0.027分钟，300景检测仅需8分钟，大大节省了人力、时间和软硬件资源，保证了数据交付时效性。
[0095]
(2)人眼判读来检测卫星影像沿轨方向局部畸变问题极易发生漏检错检的情况，11月数据漏检3景，错检13景，准确率为20％，其中错检的13景中有5景是由多次灰度重采样引起的某些角度的纹理或边缘呈锯齿状导致的错检；有8景是由坏像元插值处理后出现地物纹理或边缘扩散、彩带等导致的错检。使用现有技术1、2、3检测影像沿轨方向局部畸变，分别错检4景、6景、6景，漏检3景、2景、4景，准确率分别为65％、60％、50％，以上错检的景均为弱纹理图像匹配同名点失败导致的错判，漏检的景为同名特征点不在出现局部波段错位的行数范围内导致的漏检；而基于本技术技术的检测影像沿轨方向局部畸变问题，无错检和漏检发生，准确率可达100％，准确率大幅度提升，保证了数据交付的质量。
[0096]
(3)检测出沿轨方向局部畸变的影像通过该申请技术均得到纠正，畸变消失，可顺利交付用户使用，大大降低了用户投诉率，提高用户满意度，并减少了星地资源的浪费。
[0097]
因此，本技术实施例所提供的方案中，一方面通过自动检测方式避免影像数目较大，人工检测所导致的检测效率以及准确性不高的问题，另一方面，避免通过sift算法，互信息或相位一致性等算法获取同名特征点来进行影像自动匹配来检测过程中，由于相机各波段可能会存在坏像元、影像辐射几何处理过程中进行的若干次灰度重采样或水体、沙漠影像纹理较单一，农田、林地影像纹理较细碎等原因，影响自动匹配精度，所导致出现错检的问题。
[0098]
本技术实施例所提供的方案中，通过从元信息文件中提取待检测影像的景号、影像起始绝对行计数和影像结束绝对行计数，根据景号从成像时间数据文件中获取影像起始绝对行计数至影像结束绝对行计数部分每一行所记录的积分时间，逐行判断每一行所对应的积分时间是否发生跳变，根据积分时间跳变来确定影像可能会存在沿轨方向局部畸变，避免由于水体、沙漠影像纹理较单一，农田、林地影像纹理较细碎，影响自动匹配的精度，易漏检，导致检测准确性不高。
[0099]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
[0100]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0101]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指
令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0102]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0103]
显然，本领域的技术人员可以对本技术进行各种改动和变型而不脱离本技术的精神和范围。这样，倘若本技术的这些修改和变型属于本技术权利要求及其等同技术的范围之内，则本技术也意图包含这些改动和变型在内。