具体实施方式

为了使本申请实施例的目的、技术方案和优点更加清楚，下面将结合附图，对本申请实施例中的技术方案进行描述。

在本申请实施例的描述中，“示例性的”、“例如”或者“举例来说”等词用于表示作例子、例证或说明。本申请实施例中被描述为“示例性的”、“例如”或者“举例来说”的任何实施例或设计方案不应被解释为比其它实施例或设计方案更优选或更具优势。确切而言，使用“示例性的”、“例如”或者“举例来说”等词旨在以具体方式呈现相关概念。

在本申请实施例的描述中，术语“和/或”，仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，单独存在B，同时存在A和B这三种情况。另外，除非另有说明，术语“多个”的含义是指两个或两个以上。例如，多个系统是指两个或两个以上的系统，多个屏幕终端是指两个或两个以上的屏幕终端。

此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。术语“包括”、“包含”、“具有”及它们的变形都意味着“包括但不限于”，除非是以其他方式另外特别强调。

图1是本申请实施例提供的一种遥感卫星影像辐射校正方法的流程图。如图所示，该方法包括如下的步骤S1-步骤S5。

在步骤S1中，获取遥感卫星影像。

本实施例中，遥感卫星影像可以是未被校正的0级影像，遥感卫星影像中包括非渐晕影像块和渐晕影像块。

在一个示例中，可以按照渐晕区域的宽度，将遥感卫星影像中的渐晕区域和非渐晕区域划分为多个影像块，从而获得非渐晕影像块和渐晕影像块。示例性的，如图2a所示的遥感卫星影像图，该卫星影像图中的黑色条是渐晕区域，除黑色条之外的其他部分是非渐晕区域。可选地，非渐晕影像块和渐晕影像块可以是100×100大小的影像块。可以理解的，遥感卫星影像可以包括至少一个非渐晕影像块和至少一个渐晕影像块。

在步骤S2中，根据非渐晕影像块的整体均值和渐晕影像块的列均值，确定岭回归模型。

本实施例中，可以先计算每个非渐晕影像块中所有像素的整体均值，以及计算每个渐晕影像块中每列像素的列均值，然后对整体均值和列均值进行岭回归拟合，获得岭回归模型的参数，从而获得非渐晕影像块和渐晕影像块之间的岭回归模型。其中，列均值为所述岭回归函数的输入量，整体均值为所述岭回归函数的输出量。可以理解的，渐晕影像块中每列的列均值对应一个岭回归模型的参数。

在一个示例中，在获得岭回归模型时，可以根据非渐晕影像块的灰度共生矩阵，在至少一个非渐晕影像块中确定出纹理均匀的第一非渐晕影像块，将第一非渐晕影像块的整体均值用于岭回归拟合。具体地，在获得每个非渐晕影像块的灰度共生矩阵之后，根据灰度共生矩阵计算每个非渐晕影像块的逆差矩和熵值，然后根据预先设定的逆差矩阈值和熵阈值，在至少一个非渐晕影像块中筛选出第一非渐晕影像块。可选地，可以将逆差矩大于逆差矩阈值和熵小于熵阈值的非渐晕影像块视为纹理均匀的第一非渐晕影像块。

在步骤S3中，利用岭回归模型校正渐晕影像块的像素值，以获得第一遥感卫星影像。

本实施例中，在获得岭回归模型之后，将渐晕影像块中列均值对应的像素值输入岭回归模型中进行渐晕校正，获得岭回归模型输出的校正后的像素值。然后将渐晕影像块列均值对应的该像素值调整为岭回归模型输出的像素值，从而获得渐晕影像块对应的经过校正后的第一影像块。重复这个过程，直至渐晕影像块中列均值对应的所有像素值均被渐晕校正。可以理解的，第一遥感卫星影像是经过渐晕校正获得，因此，第一遥感卫星影像包括没有经过渐晕校正的非渐晕影像块和渐晕校正获得的第一影像块。在一个示例中，还可以将渐晕影像块中列均值对应的一列像素值全部输入岭回归模型中进行渐晕校正，对应的，岭回归模型也输出一列像素值，然后将列均值对应的一列像素值调整为岭回归模型也输出一列像素值，获得经过校正后的第一影像块。

示例性的，在将图2a的渐晕影像块的像素值代入岭回归模型进行渐晕校正，可以得到如图2b所示的校正后的卫星影像，可以看到，在图2b中，经过岭回归模型校正后，消除了其中的渐晕问题，让影像更加清晰可辨。

在步骤S4中，对第一遥感卫星影像进行影像配准，获得第二遥感卫星影像。

本实施例中，可以采用相位相关法获得第一遥感卫星影像的位移差，基于位移差对第一遥感卫星影像进行亚像素级影像配准，获得配准后的第二遥感卫星影像。具体地，可以根据获得的位移差对第一遥感卫星影像进行影像配准，从而获得第二遥感卫星影像。影像配准还可以采用本领域常规的配准方法，不限于本实施例中提供的相位相关法的配准方法。

在步骤S5中，确定第二遥感卫星影像中的重叠区域，并利用平差模型校正所述第二遥感卫星影像中重叠区域的像素值，获得第三遥感卫星影像。

本实施例中，可以将平差模型的参数初始化，然后将重叠区域的像素值输入平差模型，并采用梯度下降法迭代更新平差模型的参数，当平差模型迭代更新后的参数满足迭代终止条件时，获得平差优化模型输出的重叠区域的第二像素值；最后，将重叠区域的像素值调整为第二像素值，即可实现对重叠区域中的影像的色差校正，第二遥感卫星影像经过校正之后，即为第三遥感卫星影像。

在一个示例中，可以按照预设的像素值区间和每个区间的像素值个数从重叠区域的所有像素值中筛选出像素值样本，然后将像素值样本用于平差模型的训练，确定出平差模型的参数。然后再利用平差模型对重叠区域的所有像素值进行平差优化，平差优化的过程见前述的介绍。这样可以节省确定平差模型的时间，提高效率。可选地，为了保证平差模型的精确度，像素值区间可以设为30，每个区间的像素值可以设为500。

示例性的，如图3a、图3b和图3c所示的卫星影像，如图中的标记，三幅影像中都存在明显的色差现象，因此，在通过像素值样本确定出平差模型之后，使用平差模型对图3a、图3b和图3c所示的卫星影像中重叠区域的影像进行色差优化，即将三幅影像中重叠区域的像素值代入平差模型，然后将重叠区域的像素值调整为平差模型输出的对应的像素值。图3a、图3b和图3c所示三幅影像经色差优化获得的卫星影像分别如图4a、图4b和图4c所示，其中，图4a是图3a的影像经色差优化后的影像图，图4b是图3b的影像经色差优化后的影像图，图4c是图3c的影像经色差优化后的影像图。以图4a为例，从图中可以明显看出，与图3a相同标记位置处，本申请的方法可以消除影像的色差问题

在一个可选的实施方式中，在对第一遥感卫星影像进行影像配准之前，将第一遥感卫星影像中大于像素阈值的像素值调整为像素阈值，其中，像素阈值可以根据CCD图像传感器的最大接收亮度值设定。

图5是本申请实施例提供的一种遥感卫星影像辐射校正装置的示意图。如图5所示，该装置包括获取模块、第一校正模块、第一校正模块、配准模块和第二校正模块，其中，获取模块，用于获取遥感卫星影像；所述遥感卫星影像中包括非渐晕影像块和渐晕影像块；第一校正模块，用于根据所述非渐晕影像块的整体均值和所述渐晕影像块的列均值，确定岭回归模型，以及利用所述岭回归模型校正所述渐晕影像块的像素值，获得第一遥感卫星影像；配准模块，用于对所述第一遥感卫星影像进行影像配准，获得第二遥感卫星影像；第二校正模块，用于确定所述第二遥感卫星影像中的重叠区域，并利用平差模型校正所述重叠区域的像素值，获得第三遥感卫星影像。

可以理解的，上述校正装置仅仅是本申请实施例的一种示例，并不限值其具体划分，因此，在其他实施例中，校正装置可以包括高于或低于相比上述五个模块的模块，来实现本申请方法实施例记载的方法步骤。

在一种可能的实施方式中，所述第一校正模块具体用于：根据所述非渐晕影像块的灰度共生矩阵，确定所述非渐晕影像块为第一非渐晕影像块；对所述第一非渐晕影像块的整体均值和所述渐晕影像块的列均值进行岭回归拟合，获得所述岭回归模型。

在一种可能的实施方式中，所述第一校正模块还用于：根据所述非渐晕影像块的灰度共生矩阵，确定所述非渐晕影像块的逆差矩和熵值；根据所述非渐晕影像块的逆差矩和熵值，以及逆差矩阈值和熵阈值，确定所述非渐晕影像块为第一非渐晕影像块。

在一种可能的实施方式中，所述第一校正模块具体用于：将所述渐晕影像块中所述列均值对应的像素值输入所述岭回归模型，获得所述岭回归模型输出的像素值；将所述列像素均值对应的像素值调整为所述岭回归模型输出的像素值，获得第一影像块；所述第一遥感卫星影像包括所述非渐晕影像块和所述第一影像块。

在一种可能的实施方式中，所述配准模块还用于：在所述对所述第一遥感卫星影像进行影像配准之前，将所述第一遥感卫星影像中大于像素阈值的像素值调整为所述像素阈值。

在一种可能的实施方式中，所述第二校正模块具体用于：将所述重叠区域的像素值输入平差模型，采用梯度下降法迭代更新所述平差模型的参数；当所述平差模型的参数满足迭代终止条件时，获得所述平差优化模型输出的第二像素值；将所述重叠区域的像素值调整为所述第二像素值，获得第三遥感卫星影像。

本申请实施例还提供一种计算设备，该计算设备包括存储器和处理器。其中，存储器中存储有计算机指令，处理器读取该计算机指令时，实现上述方法实施例中的方法步骤。

可以理解的是，本申请的实施例中的处理器可以是中央处理单元(centralprocessing unit，CPU)，还可以是其他通用处理器、数字信号处理器(digital signalprocessor，DSP)、专用集成电路(application specific integrated circuit，ASIC)、现场可编程门阵列(field programmable gate array，FPGA)或者其他可编程逻辑器件、晶体管逻辑器件，硬件部件或者其任意组合。通用处理器可以是微处理器，也可以是任何常规的处理器。

本申请的实施例中的方法步骤可以通过硬件的方式来实现，也可以由处理器执行软件指令的方式来实现。软件指令可以由相应的软件模块组成，软件模块可以被存放于随机存取存储器(random access memory，RAM)、闪存、只读存储器(read-only memory，ROM)、可编程只读存储器(programmable rom，PROM)、可擦除可编程只读存储器(erasable PROM，EPROM)、电可擦除可编程只读存储器(electrically EPROM，EEPROM)、寄存器、硬盘、移动硬盘、CD-ROM或者本领域熟知的任何其它形式的存储介质中。一种示例性的存储介质耦合至处理器，从而使处理器能够从该存储介质读取信息，且可向该存储介质写入信息。当然，存储介质也可以是处理器的组成部分。处理器和存储介质可以位于ASIC中。

在上述实施例中，可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用软件实现时，可以全部或部分地以计算机程序产品的形式实现。所述计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行所述计算机程序指令时，全部或部分地产生按照本申请实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中，或者通过所述计算机可读存储介质进行传输。所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL))或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质，(例如，软盘、硬盘、磁带)、光介质(例如，DVD)、或者半导体介质(例如固态硬盘(solid state disk，SSD))等。

可以理解的是，在本申请的实施例中涉及的各种数字编号仅为描述方便进行的区分，并不用来限制本申请的实施例的范围。