[发明专利]基于独立模型法平差的卫星影像颜色相对校正方法

权利要求书

1.基于独立模型法平差的卫星影像颜色相对校正方法，其特征在于，包括如下步骤：

步骤1，对影像对的重叠区所有像素进行不变概率计算，将所有不变概率高于Hρ的像素对作为候选辐射特征点；

步骤2，对该影像对的每一张影像的每一个波段进行辐射特征点抽稀与均匀化，在该波段每一个灰度值上选取不变概率最大的Hs个候选辐射特征点，将两张影像的不同波段上选取的辐射特征点的集合作为最终的辐射特征点；

步骤3，利用正射影像仿射变换参数，对每一个辐射特征点计算平面坐标(X，Y)，然后基于平面坐标进行特征点串联，即设置一个几何误差容忍阈值Ht，若两个特征点的地面坐标欧式距离小于该阈值，则认为是辐射连接点；

步骤4，计算“虚拟辐射控制点”，对所有辐射连接点利用式(1)进行“虚拟辐射控制点”计算，

其中，R_Truth表示“虚拟辐射控制点”的灰度值，R_i表示辐射连接点在影像i上的观测值，(a_i，b_i)表示前一次迭代计算后影像i的线性调整模型，P_i是影像i的权值，表示影像的质量，所有影像的初始权值都设置为1，m表示该连接点的观测值个数；

步骤5，影像内线性颜色调整模型(a_i，b_i)计算：选择线性模型估计器，对每一个影像，利用该影像的辐射连接点观测值以及由步骤4计算得到的“虚拟辐射控制点”进行线性模型拟合，得到(a_i，b_i)的值；

步骤6，逐影像计算观测值中误差σ_i，如式(2)，并利用式(3)计算整体中误差σ⁰，

式(2)中，K表示影像i的有效观测值数量，Δ_j是连接点j的观测值的残差，R_i，j表示该连接点j在影像i中的观测值，中误差σ_i表示影像i的观测值总体误差的大小，式(3)中，M表示影像的数量；

比较本次迭代与前一次迭代的整体中误差，若本次迭代中误差大于前一次迭代，或者两次迭代的差值的绝对值小于阈值Hσ，则迭代终止；

步骤7，以影像中误差和整体中误差的较大值作为粗差剔除的基础，以式(4)作为阈值，逐影像进行粗差探测，当|Δ|＞H_outliers，则认为该观测值是粗差像素，

H_outliers＝t*max(σ⁰，σ_i) (4)

式中，t表示粗差的容忍度，|Δ|是指观测值残差的绝对值；

步骤8，利用影像中误差以及整体中误差，更新影像权值，如式(5)，

步骤9，重复步骤4-8，直到满足步骤6中迭代条件，或者迭代次数超过一定次数，则迭代终止。

2.如权利要求1所述的基于独立模型法平差的卫星影像颜色相对校正方法，其特征在于：步骤1中采用迭代加权多元变化检测算法对影像对的重叠区所有像素进行不变概率计算，Hρ设置为0.9。

3.如权利要求1所述的基于独立模型法平差的卫星影像颜色相对校正方法，其特征在于：步骤2中Hs设置为100。

4.如权利要求1所述的基于独立模型法平差的卫星影像颜色相对校正方法，其特征在于：步骤3中gsd为地面分辨率。

5.如权利要求1所述的基于独立模型法平差的卫星影像颜色相对校正方法，其特征在于：步骤5中的线性模型估计器包括最小二乘、Theil-Sen线性回归以及正交线性回归。

6.如权利要求1所述的基于独立模型法平差的卫星影像颜色相对校正方法，其特征在于：步骤6中Hσ＝10^-6。

7.如权利要求1所述的基于独立模型法平差的卫星影像颜色相对校正方法，其特征在于：步骤7中，前三次迭代中t设置为5，后续迭代中则设置为3。