一种卫星遥感影像物理匀色方法与流程

本发明涉及航天遥感卫星领域，尤指一种卫星遥感影像物理匀色方法。

背景技术：

为了保证卫星传感器既能获取高空间分辨率卫星影像，同时能够具备宽覆盖对地观测能力，人们经常采用多相机拼接的方法实现该目的。例如，gf-1卫星搭载的宽幅盖wfv传感器采用4台相机拼接成像技术，如图1所示，实现了高分辨率与大视场相结合的研制总目标，其能够获取空间分辨率为16m，幅宽优于800km的卫星影像。然而，由于各相机相关参数(相机观测角度、相机光谱响应函数等)有所不同，导致不同相机获取的卫星影像重叠区域出现了较为明显的辐射不一致性问题，即存在色差现象。因此若要使同一传感器不同相机获取的影像进行拼接处理后不存在色差问题，需要对其结果进行影像匀色处理。

目前，已经有许多研究学者针对该问题提出了相应的解决方法。例如，王建忠提出将小波变换用于影像匀色处理。朱述龙等人在对小波变换匀色方法和基于重叠影像的匀色方法的优劣性进行分析后，提出了一种拼接缝消除的强制改正方法。易尧华等人在对常用的影像色调平衡方法进行总结后，提出了一种适用于大型影像数据库中影像匀色方法。周延刚提出了一种基于色度空间变换和动态宽度的遥感影像拼接缝消除方法，并利用航空影像验证了算法的可行性。李德仁等人利用wallis滤波方法实现了对多幅影像的匀色处理。zomet等人通过对影像接边处的对比度进行平滑处理，消除了影像拼接缝处色差。葛仕明等人提出的基于梯度场的匀色方法，在消除接边处色差的同时，能够较好的处理影像几何错位问题。孙黎明等人利用直方图配准算法对cbers-02bccd影像拼接产生的色差进行匀色处理。周丽雅等人在对传统的wallis匀光方法进行改进后，提出了一种反差一致性保持的影像匀光算法，该方法通过对影像的合理分块及自动统计标准均值 和方差来实现影像的匀光处理。田金炎等人提出了一种基于wallis与距离权重增强的航空影像匀色方法。

虽然采用上述匀色方法均可以实现拼接后影像接边处无明显色差的目视效果，但是这些方法都会破坏影像原有的辐射信息，且拼接结果不具备任何物理意义，这将无法开展进一步的遥感定量化产品生产及相关研究工作。因此，从gf-1卫星宽覆盖wfv传感器影像匀色过程中遇到的实际问题出发，本发明考虑利用wfv影像定标系数将原始dn值数据统一转换成具有辐射物理意义的表观辐亮度信息，从而完成影像匀色处理。但是由于官方采用的反射率基法获取各相机影像定标系数的过程是相互独立的，这将无法保证wfv影像重叠区域具有较好的辐射一致性。

技术实现要素：

针对现有算法无法满足对卫星影像进行匀色处理时不破坏影像原始辐射信息的要求，本发明提出了一种基于交叉辐射定标的影像物理匀色方法，该方法在不改变原始影像辐射信息的同时，可以有效的降低相邻相机影像重叠区域辐射不一致性问题，实现影像物理匀色处理。

本发明提出了一种卫星遥感影像物理匀色方法，其特征在于，包括以下步骤：

1).利用互交叉辐射定标方法获取某一台wfv相机绝对辐射定标系数；

2).采用统计分析方法建立相邻相机影像的dn值相关模型；

3).采用仿真模拟的方法建立相邻相机影像的表观辐亮度相关模型

4).利用自交叉辐射定标方法获取其它wfv相机绝对辐射定标系数。

进一步，该方法还包括步骤5)利用wfv相机绝对辐射定标系数实现wfv影像物理匀色。

进一步，其特征在于，步骤1)具体为：采用目前常用的交叉定标方法，以辐射性能较高的modis影像为参考基准，将敦煌辐射校正场作为传递平台，获取某一台wfv相机绝对辐射定标系数。由于该过 程是跨平台的辐射性能传递，因此也可称为互交叉定标。

进一步，其特征在于，步骤2)具体为：选取不同时相、不同区域、同一轨道的相邻相机影像，对影像重叠区域中光谱特性均一的地物影像dn值进行统计，建立相邻相机影像dn值之间的相关关系模型。

进一步，其特征在于，步骤3)具体为：利用多种典型地物光谱实测曲线，结合影像重叠区域观测几何特性，基于6s辐射传输模型进行仿真模拟，建立相邻相机影像表观辐亮度之间的相关关系模型。

进一步，其特征在于，步骤4)具体为：利用自交叉辐射定标方法将已完成互交叉辐射定标的相机辐射特性传递给其它相机，获取其它相机各自的绝对辐射定标系数。

进一步，其特征在于，步骤5)具体为：利用上述各相机绝对辐射定标系数，将原始wfv影像dn值信息转换为表观辐亮度信息，进而实现影像物理匀色处理。

本发明的基于交叉辐射定标的影像物理匀色方法可以有效的实现同一传感器不同相机间影像物理匀色处理，除了不破坏原始影像辐射信息外，还有许多优势，如匀色方法简单、处理速度快、不受季节影响、等等。

附图说明

图1wfv传感器成像示意图；

图2卫星遥感物理匀色流程图；

图3gf-1卫星wfv4相机互交叉辐射定标结果(a为第一波段结果；b为第二波段结果；c为第三波段结果；d为第四波段结果)；

图4wfv3与wfv4相机影像dn值相关模型图(a为第一波段结果；b为第二波段结果；c为第三波段结果；d为第四波段结果)；

图5wfv3与wfv4相机影像表观辐亮度相关模型图(a为第一波段结果；b为第二波段结果；c为第三波段结果；d为第四波段结果)；

图6wfv3与wfv4相机影像重叠区域表观辐亮度差异的绝对值图(a为第一波段结果；b为第二波段结果；c为第三波段结果；d为第四波段结果)；

图7wfv1与wfv2相机匀色前后对比图(a为原始影像；b为利 用官方首次定标系数匀色结果；c为利用官方最新定标系数匀色结果；d为本发明物理匀色结果)；

图8wfv2与wfv3相机匀色前后对比图(a为原始影像；b为利用官方首次定标系数匀色结果；c为利用官方最新定标系数匀色结果；d为本发明物理匀色结果)；

图9wfv3与wfv4相机匀色前后对比图(a为原始影像；b为利用官方首次定标系数匀色结果；c为利用官方最新定标系数匀色结果；d为本发明物理匀色结果)。

具体实施方式

如图2所示本发明的一种卫星遥感影像物理匀色方法包括传统的交叉辐射定标、相邻相机影像dn相关性统计、相邻相机影像表观辐亮度相关性仿真、自交叉辐射定标以及物理匀色。本发明以wfv3和wfv4为例，阐述影像物理匀色方法过程。

1)基于modis影像的某一台wfv相机交叉辐射定标，以wfv4相机为例。首先收集敦煌地区同步或准同步的影像对数据，并利用数据筛选限制条件获取有效的影像对。其中数据筛选条件包括：a)modis与wfv影像在敦煌辐射校正场的过境时间差小于1h；b)卫星过境时，试验区上方不能有云的遮挡；c)wfv影像应该覆盖敦煌辐射校正场的中心附近的区域；d)敦煌辐射校正场不应该位于modis影像左右边缘处。本发明共得到了自卫星发射后一年左右时间内敦煌辐射校正场地区的7景有效影像对数据。然后利用传统的交叉辐射定标方法获取wfv4相机绝对辐射定标系数，如图3所示。其中wfv4相机定标公式如公式1所示。

l^lwfv4＝aⁱdn4_to_l4·dn^lwfv4+b^ldn4_to_l4(1)

式中，dnⁱwfv4和lⁱwfv4分别代表wfv4影像第i波段的dn值和表观辐亮度，aⁱdn4_to_l4和bⁱdn4_to_l4为第i波段的定标系数(i＝1，2，3，4，下同)。

2)建立相邻相机影像dn值相关性模型，以wfv3和wfv4相机为例。通过对大量重叠区域wfv3和wfv4影像dn值进行统计，从图4中可以看出，其dn值之间存在很高的线性相关性。因此，wfv3和wfv4 影像dn值之间的相关模型可用公式2表示。

dnⁱwfv3＝cⁱdn4_to_dn3·dnⁱwfv4+dⁱdn4_to_dn3(2)

3)建立相邻相机影像表观辐亮度相关性模型，以wfv3和wfv4相机为例。对于影像重叠区域地物而言，6s辐射传输模型的输入参数大多数是相同的，仅地表反射率和观测角度有所差异。本发明统计多组影像xml文件中卫星角度信息并结合wfv传感器成像示意图，如图1所示，可以计算出重叠区域地物的各卫星观测角度。选取了8种典型地物实测反射率数据(人工黑/白靶标、塑胶场地、水泥地、草地、裸土、花岗岩、人工草坪)，将其与wfv3和wfv4相机光谱响应函数进行卷积处理，获取各波段等效地表反射率。另外，设定太阳天顶角为30°、太阳方位角为160°、wfv3相机卫星天顶角为16.22°、wfv4相机卫星天顶角为17.78°、卫星方位角为280°、大气模式为中纬度夏季、气溶胶模型为大陆型、550nmaod为0.15、地表高程为1500m。通过对wfv3和wfv4影像表观辐亮度统计发现，两者有很高的线性相关特征，如图4所示。因此两者表观辐亮度可用公式3表示。

lⁱwfv3＝eⁱl4_to_l3·lⁱwfv4+fⁱl4_to_l3(3)

式中，lⁱwfv3和lⁱwfv4代表wfv3和wfv4影像第i波段表观辐亮度，eⁱl4_to_l3和fⁱl4_to_l3为第i波段表观辐亮度线性拟合系数。

4)自交叉辐射定标。假设wfv3相机定标公式如公式4所示。在获取了wfv4相机在轨互交叉辐射定标系数以及wfv3和wfv4影像dn值、表观辐亮度的相关模型后，将公式1、2、3带入公式4中后得到公式5，利用公式5采用最小二乘方法获取wfv3相机绝对辐射定标系数。

lⁱwfv3＝gⁱdn3_to_l3·dnⁱwfv3+hⁱdn3_to_l3(4)

式中，dnⁱwfv3和lⁱwfv3代表wfv3影像第i波段的dn值和表观辐亮度，gⁱdn3_to_l3和hⁱdn3_to_l4为待求的wfv3影像第i波段的辐射定标系数。

在同一轨道、不同时相wfv3和wfv4影像重叠区中提取了多个光 谱均一的地物作为检测点，分别利用官方首次/最新公布的定标系数和本发明获得定标结果计算其各自表观辐亮度，并统计两相机表观辐亮度差异的绝对值。从图6可以发现，与官方定标系数相比，利用本发明获得的定标结果可以大大降低相机间辐射不一致性。对于四个波段而言，与官方首次公布的定标系数相比，其辐射不一致性分别降低了81.5％、34.07％、72.34％和83.83％。与官方最新系数相比，其辐射不一致性分别降低了22.82％、81.37％、92.17％和15.95％。

利用相同的辐射传递方法获取wfv2和wfv1相机绝对辐射定标系数。

表1gf-1卫星互/自交叉辐射定标系数

5)同轨4相机影像物理匀色。随机选取一组wfv传感器影像进行物理匀色，其成像日期为2014年11月18日。在利用平面平差方法消除原始影像几何错位信息后，采用官方首次/最新公布的定标系数和本发明获得的定标结果将其分别转换为表观辐亮度，选取局部区域进行展示如图6、7、8所示。

从结果中可以看出，官方首次/最新公布的定标系数相比，利用本发明得到的定标系数可以很好的消除原始影像辐射不一致性问题。与传统匀色方法相比，本发明的匀色过程并不破坏原始影像辐射信息，从而实现卫星遥感影像物理匀色处理。