[发明专利]一种基于空域的自适应MTFC遥感图像复原方法

权利要求书

1.一种基于空域的自适应MTFC遥感图像复原方法，其特征在于步骤如下：

（1）在地面通过实测得到遥感器系统从零频到奈奎斯特频率范围内的MTF值；

（2）根据成像需求确定理想的MTF提升曲线，并与步骤（1）的结果进行比较，得到遥感器系统从零频到奈奎斯特频率范围内不同频率处MTF的补偿值|H(e^jw)|；

（3）利用步骤（2）的结果，计算得到一维卷积系数h(k)，|H(e^jw)|=A(w,k)*h(k)，其中h(k)是中心对称的，h(k)的长度为M；w为频率，k取

（4）利用步骤（3）得到的一维卷积系数h(k)构造二维卷积核h(i,j)，h(i,j)为一个M*M的正方形矩阵且满足中心对称和轴对称，其中i,j的取值范围分别为h(i,j)中不同的卷积系数的个数为象限的二维频域响应函数表示为：

|H(fx,fy)|=Σtx=0M-12Σty=0txx(tx,ty)×k(tx,ty)[cos(fx×tx)×cos(fy×ty)+cos(fy×tx)×cos(fx×ty)]/2]]>

式中|H(fx,fy)|=H(fx2+fy2)=Σk=0M-12h(k)·2cos(M-12-k)fx2+fy2,]]>x(t_x，t_y)为一个长度的函数，其中t_y∈[0,t_x]，|H(f_x,f_y)|为二维频率处的MTF补偿值，k(t_x,t_y)为卷积系数的重复次数；

（5）实测遥感器系统在不同成像参数下不同灰度下的SNR，利用拟合的方式得到遥感器系统在不同成像参数下不同灰度下的SNR曲线；

（6）根据步骤（5）得到的SNR曲线构造抑噪查找表系数M(i)，M(i)=10^((y(i)-y_max)/20)，i在[0,2^N-1]内取整，N为量化位数，i为某一确定成像参数下的某一DN值，y(i)为i与对应的SNR，y_max为同一个确定的成像参数下的最大SNR；

（7）利用步骤（4）得到的二维卷积核h(i,j)对待复原的图像进行卷积，然后进行噪声抑制后输出复原后的图像；所述复原后的图像中的每一个像元的灰度值i₂满足条件i₂=i₀+k(i₁-i₀)，k=K||i₁-i₀||M([i₀])，其中i₀为像元的原始灰度值，i₁为为像元经过二维卷积后的灰度值，[]表示取整函数，M([i₀])通过查找步骤（6）得到的抑噪查找表系数获得，K为大于0的整数，是可变的抑噪阈值参数，根据当前像素补偿前的SNR而定，K的取值应保证补偿后该像素DN值的均匀景物SNR下降不超过设定值。