[發(fā)明專利]一種基于空域的自適應(yīng)MTFC遙感圖像復(fù)原方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201410086318.6
• 申請(qǐng)日: 2014-03-10
• 公開(kāi)（公告）號(hào): CN103871033A
• 公開(kāi)（公告）日: 2014-06-18
• 發(fā)明（設(shè)計(jì)）人: 何紅艷;齊文雯;周楠;周虎;黃巧林;岳春宇;邢坤;鮑云飛;梅強(qiáng)
• 申請(qǐng)（專利權(quán)）人: 北京空間機(jī)電研究所
• 主分類號(hào): G06T5/00
• 分類號(hào): G06T5/00
• 代理公司: 中國(guó)航天科技專利中心 11009
• 代理人: 陳鵬

1.一種基于空域的自適應(yīng)MTFC遙感圖像復(fù)原方法，其特征在于步驟如下：

（1）在地面通過(guò)實(shí)測(cè)得到遙感器系統(tǒng)從零頻到奈奎斯特頻率范圍內(nèi)的MTF值；

（2）根據(jù)成像需求確定理想的MTF提升曲線，并與步驟（1）的結(jié)果進(jìn)行比較，得到遙感器系統(tǒng)從零頻到奈奎斯特頻率范圍內(nèi)不同頻率處MTF的補(bǔ)償值|H(e^jw)|；

（3）利用步驟（2）的結(jié)果，計(jì)算得到一維卷積系數(shù)h(k)，|H(e^jw)|=A(w,k)*h(k)，其中h(k)是中心對(duì)稱的，h(k)的長(zhǎng)度為M；w為頻率，k取

（4）利用步驟（3）得到的一維卷積系數(shù)h(k)構(gòu)造二維卷積核h(i,j)，h(i,j)為一個(gè)M*M的正方形矩陣且滿足中心對(duì)稱和軸對(duì)稱，其中i,j的取值范圍分別為h(i,j)中不同的卷積系數(shù)的個(gè)數(shù)為象限的二維頻域響應(yīng)函數(shù)表示為：

|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)為一個(gè)長(zhǎng)度的函數(shù)，其中t_y∈[0,t_x]，|H(f_x,f_y)|為二維頻率處的MTF補(bǔ)償值，k(t_x,t_y)為卷積系數(shù)的重復(fù)次數(shù)；

（5）實(shí)測(cè)遙感器系統(tǒng)在不同成像參數(shù)下不同灰度下的SNR，利用擬合的方式得到遙感器系統(tǒng)在不同成像參數(shù)下不同灰度下的SNR曲線；

（6）根據(jù)步驟（5）得到的SNR曲線構(gòu)造抑噪查找表系數(shù)M(i)，M(i)=10^((y(i)-y_max)/20)，i在[0,2^N-1]內(nèi)取整，N為量化位數(shù)，i為某一確定成像參數(shù)下的某一DN值，y(i)為i與對(duì)應(yīng)的SNR，y_max為同一個(gè)確定的成像參數(shù)下的最大SNR；

（7）利用步驟（4）得到的二維卷積核h(i,j)對(duì)待復(fù)原的圖像進(jìn)行卷積，然后進(jìn)行噪聲抑制后輸出復(fù)原后的圖像；所述復(fù)原后的圖像中的每一個(gè)像元的灰度值i₂滿足條件i₂=i₀+k(i₁-i₀)，k=K||i₁-i₀||M([i₀])，其中i₀為像元的原始灰度值，i₁為為像元經(jīng)過(guò)二維卷積后的灰度值，[]表示取整函數(shù)，M([i₀])通過(guò)查找步驟（6）得到的抑噪查找表系數(shù)獲得，K為大于0的整數(shù)，是可變的抑噪閾值參數(shù)，根據(jù)當(dāng)前像素補(bǔ)償前的SNR而定，K的取值應(yīng)保證補(bǔ)償后該像素DN值的均勻景物SNR下降不超過(guò)設(shè)定值。