[发明专利]一种航空遥感影像的最小二乘区域网匀色方法

摘要

本发明公开了一种航空遥感影像的最小二乘区域网匀色方法，该方法借鉴光束法区域网空中三角测量的思想，首先统计单像正射纠正影像的相邻影像叠加区域的像素色调信息，再构建每幅影像像素色调信息补偿参数的方程组，利用最小二乘平差方法，解算每幅影像的像素色调信息补偿参数，最后根据解算得到的色调补偿参数调整每幅影像的色调信息。本发明方法能够对大测区以及测区内影像内容差异较大的影像进行色调差异消除处理，在消除色差的过程中还能够方便地指定测区内部任意位置、任意数量的影像作为色调控制影像来控制整个测区的基本色调。

主权项

1.一种航空遥感影像的最小二乘区域网匀色方法，其特征在于，依次包括以下步骤：步骤1，确认测区内所有待镶嵌的单像正射纠正影像在镶嵌过程中的接边关系，并获得q对接边影像对；步骤2，对步骤1获得的q对接边影像对，按照地理坐标对每对接边影像对中的两张单像正射纠正影像进行叠加，叠加后排除叠加区域内单像正射纠正影像的背景像素和纹理差异较大的像素，然后，分别统计q对接边影像对的叠加区域内未被排除的像素数和包括像素灰度均值和像素灰度方差的像素色调信息；步骤3，计算测区内所有待镶嵌的单像正射纠正影像的像素色调信息补偿参数，该步骤进一步包括以下子步骤：3-1、利用步骤2统计得到的像素色调信息构建矩阵方程：A_BX＝L_B，其中，X＝[x₀ x₁ x₂ ... x_p-1]^T，x₀、x₁、x₂、x₃、...x_p-1分别为测区内所有待镶嵌的单像正射纠正影像的像素色调信息补偿参数，p为所有待镶嵌影像的张数；A_B为q×p矩阵，且AB=1-10...01-10...0...1-10...;]]>LB=a01-a00a11-a10a21-a20...aq-11-aq-10T,]]>为第i对接边影像对中的第一张影像x在叠加区域内的像素灰度均值或像素灰度方差，为第i对接边影像对中的第二张影像y在叠加区域内的像素灰度均值或像素灰度方差，且和同为像素灰度均值或像素灰度方差；3-2、将p张单像正射纠正影像中色调较好的t张影像设置为色调控制影像，并构建矩阵方程：A_CX＝L_C，其中，X＝[x₀ x₁ x₂ ... x_p-1]^T，x₀、x₁、x₂、x₃、...x_p-1分别为p张单像正射纠正影像的像素色调信息补偿参数；A_c[i，j]＝c_i，j，i∈[0，t-1]，j∈[0，p-1]，i表示在t张色调控制影像中的第i张影像，当j＝i时，c_i，j＝1，其它情况下，c_i，j＝0；L_C为t行零常数的列向量，L_C＝[0 0…0]^T；3-3、令A=ABAC,]]>L=LBLC,]]>构建矩阵方程AX＝L，则p张单像正射纠正影像的像素色调信息补偿参数X＝(A^TA)^-1A^TL；3-4、构建权矩阵P=p00···00p1000······00···0pq-100E]]>修正步骤3-3得到的像素色调信息补偿参数，则p张单像正射纠正影像的色调补偿参数X＝(A^TPA)^-1A^TPL，其中，r∈[0，q-1]，t_r为叠加区域内的有效像素数，E为(p-q)阶单位矩阵；步骤4，根据步骤3所得的像素色调补偿参数，调整每张单向正射纠正影像的像素色调信息。