本公開涉及信號處理技術(shù)領域，具體地，涉及一種星載合成孔徑雷達圖像的模糊抑制方法和裝置。

背景技術(shù)：

星載合成孔徑雷達是一種可以對地球進行觀測的高分辨率遙感探測雷達。由于星載合成孔徑雷達衛(wèi)星可以不受雨霧、暗夜等能見度低的氣象條件的限制，因此，可以通過它實現(xiàn)全天時、全天候、高分辨率、寬幅對地觀測。正因為這樣，星仔合成孔徑雷達在農(nóng)業(yè)、海洋、災害監(jiān)測、3d繪圖等領域具有廣闊的應用前景。

星載合成孔徑雷達具有聚束模式、條帶模式和滑動聚束模式三種工作模式，其中，條帶模式是它最主要的工作模式。然而，星載合成孔徑雷達本身的工作機理決定了雷達圖像中可能存在方位模糊的現(xiàn)象，這將會明顯降低星載合成孔徑雷達圖像的質(zhì)量，后續(xù)的與星載合成孔徑雷達圖像相關的應用也會受到一定程度的影響。

相關技術(shù)中，通常采用理想濾波器法和維納濾波器法等方法對星載合成孔徑雷達圖像進行模糊抑制，但這些方法均難以將模糊信號抑制干凈，甚至會損失圖像分辨率和信噪比。

技術(shù)實現(xiàn)要素：

本公開的目的是針對現(xiàn)有技術(shù)中的方位模糊信號抑制不徹底的問題，提供一種星載合成孔徑雷達圖像的模糊抑制方法和裝置。

為了實現(xiàn)上述目的，本公開提供一種星載合成孔徑雷達圖像的模糊抑制方法，所述方法包括：

對星載合成孔徑雷達圖像進行方位模糊重建，確定出所述星載合成孔徑雷達圖像的模糊圖像；

根據(jù)與所述星載合成孔徑雷達圖像中的數(shù)據(jù)點對應的所述模糊圖像的數(shù)據(jù)點，判斷所述星載合成孔徑雷達圖像中、滿足所述預定抑制條件的數(shù)據(jù)點是否滿足預定抑制條件；

當所述星載合成孔徑雷達圖像中的數(shù)據(jù)點滿足所述預定抑制條件時，利用所述模糊圖像，對所述星載合成孔徑雷達圖像中的數(shù)據(jù)點進行模糊抑制。

可選地，所述模糊圖像包括負頻率方位模糊圖像；

所述對星載合成孔徑雷達圖像進行方位模糊重建，確定出所述星載合成孔徑雷達圖像的模糊圖像，包括：

對所述星載合成孔徑雷達圖像s0(τ,η)按列做快速傅里葉變換，得到圖像s1(τ,fη)，其中，j為列索引，nr為圖像矩陣總列數(shù)，fs為信號采樣率，i為行索引，na為圖像矩陣總行數(shù)，fp為脈沖重復頻率，τ為距離向時間，η為方位向時間，fη為方位向頻率；

對所述圖像s1(τ,fη)按列進行方位相位補償和天線方向圖補償，每一列數(shù)據(jù)乘以傳遞函數(shù)h1，得到圖像s2(τ,fη)：

s2(τ,fη)＝h1·s1(τ,fη)

其中，h1為傳遞函數(shù)，f0為信號載頻，r0為斜距，c為光速，rref為場景中心斜距，λ為波長，vr為等效速度，fηc為多普勒中心頻率，d(fη-fp)為徙動因子，wa(fη-fηc-fp)為衛(wèi)星方位向天線方向圖的頻域表達式，wa(τ)為方向圖的時域表達式；

對所述圖像s2(τ,fη)按行做快速傅里葉變換，得到圖像s3(fτ,fη)；

對所述圖像s3(fτ,fη)按行進行一致模糊距離徙動校正和方位向調(diào)制，得到圖像s4(fτ,fη)；

對所述圖像s4(fτ,fη)按行做快速逆傅里葉變換，得到圖像s5(τ,fη)；

對所述圖像s5(τ,fη)按行進行補余模糊距離徙動校正，得到圖像s6(τ,fη)；

對所述圖像s6(τ,fη)按行做快速傅里葉變換，得到圖像s7(fτ,fη)；

對所述圖像s7(fτ,fη)按行進行脈沖壓縮和二次距離壓縮誤差校正，得到圖像s8(fτ,fη)；

對所述圖像s8(fτ,fη)按行做快速逆傅里葉變換，得到圖像s9(τ,fη)；

對所述圖像s9(τ,fη)補償殘余誤差和進行方位散焦處理，得到圖像s10(τ,fη)；

對所述圖像s10(τ,fη)按列做快速逆傅里葉變換，得到所述負頻率方位模糊圖像。

可選地，所述模糊圖像包括正頻率方位模糊圖像；

所述對星載合成孔徑雷達圖像進行方位模糊重建，確定出所述星載合成孔徑雷達圖像的模糊圖像，包括：

對所述星載合成孔徑雷達圖像s0(τ,η)按列做快速傅里葉變換，得到圖像s11(τ,fη)，其中，j為列索引，nr為圖像矩陣總列數(shù)，fs為信號采樣率，i為行索引，na為圖像矩陣總行數(shù)，fp為脈沖重復頻率，τ為距離向時間η為方位向時間，fη為方位向頻率；

對所述圖像s11(τ,fη)按列進行方位相位補償和天線方向圖補償，得到圖像s12(τ,fη)：

s12(τ,fη)＝h1·s11(τ,fη)

其中，h1為傳遞函數(shù)，f0為信號載頻，r0為斜距，c為光速，rref為場景中心斜距，λ為波長，vr為等效速度，fηc為多普勒中心頻率，d(fη+fp)為徙動因子，wa(fη-fηc+fp)為衛(wèi)星方位向天線方向圖的頻域表達式，wa(τ)為方向圖的時域表達式；

對所述圖像s12(τ,fη)按行做快速傅里葉變換，得到圖像s13(fτ,fη)；

對所述圖像s13(fτ,fη)按行進行一致模糊距離徙動校正和方位向調(diào)制，得到圖像s14(fτ,fη)；

對所述圖像s14(fτ,fη)按行做快速逆傅里葉變換，得到圖像s15(τ,fη)；

對所述圖像s15(τ,fη)按行進行補余模糊距離徙動校正，得到圖像s16(τ,fη)；

對所述圖像s16(τ,fη)按行做快速傅里葉變換，得到圖像s17(fτ,fη)；

對所述圖像s17(fτ,fη)按行進行脈沖壓縮和二次距離壓縮誤差校正，得到圖像s18(fτ,fη)；

對所述圖像s18(fτ,fη)按行做快速逆傅里葉變換，得到圖像s19(τ,fη)；

對所述圖像s19(τ,fη)補償殘余誤差和進行方位散焦處理，得到圖像s20(τ,fη)；

對所述圖像s20(τ,fη)按列做快速逆傅里葉變換，得到所述正頻率方位模糊圖像。

可選地，所述根據(jù)與所述星載合成孔徑雷達圖像中的數(shù)據(jù)點對應的所述模糊圖像的數(shù)據(jù)點，判斷所述星載合成孔徑雷達圖像中的數(shù)據(jù)點是否滿足預定抑制條件，包括：

判斷所述模糊圖像的數(shù)據(jù)點的絕對值是否大于所述星載合成孔徑雷達圖像中對應的數(shù)據(jù)點的絕對值的預定倍數(shù)，其中，所述預定倍數(shù)大于0且小于1；

判斷所述方位向時間是否滿足預設的閾值條件；

當所述模糊圖像的數(shù)據(jù)點的絕對值大于所述星載合成孔徑雷達圖像中對應的數(shù)據(jù)點的絕對值的預定倍數(shù)、且所述方位向時間滿足所述閾值條件時，確定所述星載合成孔徑雷達圖像中的數(shù)據(jù)點滿足所述預定抑制條件。

可選地，所述模糊圖像包括負頻率方位模糊圖像；

所述閾值條件包括：

其中，η為所述方位向時間；fp為脈沖重復頻率；fr為多普勒調(diào)頻率；和/或

所述模糊圖像包括正頻率方位模糊圖像；

所述閾值條件包括：

其中，η為所述方位向時間；na為圖像矩陣總行數(shù)；fp為脈沖重復頻率；fr為多普勒調(diào)頻率。

可選地，所述當所述星載合成孔徑雷達圖像中的數(shù)據(jù)點滿足所述預定抑制條件時，利用所述模糊圖像，對所述星載合成孔徑雷達圖像中的數(shù)據(jù)點進行模糊抑制，包括：

將所述星載合成孔徑雷達圖像中、滿足所述預定抑制條件的數(shù)據(jù)點的絕對值，與所述模糊圖像中對應的數(shù)據(jù)點的絕對值做差值，得到經(jīng)模糊抑制后的圖像。

本公開還提供一種星載合成孔徑雷達圖像的模糊抑制裝置，所述裝置包括：

方位模糊重建模塊，用于對星載合成孔徑雷達圖像進行方位模糊重建，確定出所述星載合成孔徑雷達圖像的模糊圖像；

判斷模塊，用于根據(jù)與所述星載合成孔徑雷達圖像中的數(shù)據(jù)點對應的所述模糊圖像的數(shù)據(jù)點，判斷所述星載合成孔徑雷達圖像中的數(shù)據(jù)點是否滿足預定抑制條件；

模糊抑制模塊，用于當所述星載合成孔徑雷達圖像中的數(shù)據(jù)點滿足所述預定抑制條件時，利用所述模糊圖像，對所述星載合成孔徑雷達圖像中、滿足所述預定抑制條件的數(shù)據(jù)點進行模糊抑制。

可選地，所述模糊圖像包括負頻率方位模糊圖像；

所述方位模糊重建模塊用于：

對所述星載合成孔徑雷達圖像s0(τ,η)按列做快速傅里葉變換，得到圖像s1(τ,fη)，其中，j為列索引，nr為圖像矩陣總列數(shù)，fs為信號采樣率，i為行索引，na為圖像矩陣總行數(shù)，fp為脈沖重復頻率，τ為距離向時間，η為方位向時間，fη為方位向頻率；

對所述圖像s1(τ,fη)按列進行方位相位補償和天線方向圖補償，得到圖像s2(τ,fη)：

s2(τ,fη)＝h1·s1(τ,fη)

其中，h1為傳遞函數(shù)，f0為信號載頻，r0為斜距，c為光速，rref為場景中心斜距，λ為波長，vr為等效速度，fηc為多普勒中心頻率，d(fη-fp)為徙動因子，wa(fη-fηc-fp)為衛(wèi)星方位向天線方向圖的頻域表達式，wa(τ)為方向圖的時域表達式；

對所述圖像s2(τ,fη)按行做快速傅里葉變換，得到圖像s3(fτ,fη)；

對所述圖像s3(fτ,fη)按行進行一致模糊距離徙動校正和方位向調(diào)制，得到圖像s4(fτ,fη)；

對所述圖像s4(fτ,fη)按行做快速逆傅里葉變換，得到圖像s5(τ,fη)；

對所述圖像s5(τ,fη)按行進行補余模糊距離徙動校正，得到圖像s6(τ,fη)；

對所述圖像s6(τ,fη)按行做快速傅里葉變換，得到圖像s7(fτ,fη)；

對所述圖像s7(fτ,fη)按行進行脈沖壓縮和二次距離壓縮誤差校正，得到圖像s8(fτ,fη)；

對所述圖像s8(fτ,fη)按行做快速逆傅里葉變換，得到圖像s9(τ,fη)；

對所述圖像s9(τ,fη)補償殘余誤差和進行方位散焦處理，得到圖像s10(τ,fη)；

對所述圖像s10(τ,fη)按列做快速逆傅里葉變換，得到所述負頻率方位模糊圖像。

可選地，所述模糊圖像包括正頻率方位模糊圖像；

所述方位模糊重建模塊用于：

對所述星載合成孔徑雷達圖像s0(τ,η)按列做快速傅里葉變換，得到圖像s11(τ,fη)，其中，j為列索引，nr為圖像矩陣總列數(shù)，fs為信號采樣率，i為行索引，na為圖像矩陣總行數(shù)，fp為脈沖重復頻率，τ為距離向時間，η為方位向時間，fη為方位向頻率；

對所述圖像s11(τ,fη)按列進行方位相位補償和天線方向圖補償，得到圖像s12(τ,fη)：

s12(τ,fη)＝h1·s11(τ,fη)

其中，h1為傳遞函數(shù)，f0為信號載頻，r0為斜距，c為光速，rref為場景中心斜距，λ為波長，vr為等效速度，fηc為多普勒中心頻率，d(fη+fp)為徙動因子，wa(fη-fηc+fp)為衛(wèi)星方位向天線方向圖的頻域表達式，wa(τ)為方向圖的時域表達式；

對所述圖像s12(τ,fη)按行做快速傅里葉變換，得到圖像s13(fτ,fη)；

對所述圖像s13(fτ,fη)按行進行一致模糊距離徙動校正和方位向調(diào)制，得到圖像s14(fτ,fη)；

對所述圖像s14(fτ,fη)按行做快速逆傅里葉變換，得到圖像s15(τ,fη)；

對所述圖像s15(τ,fη)按行進行補余模糊距離徙動校正，得到圖像s16(τ,fη)；

對所述圖像s16(τ,fη)按行做快速傅里葉變換，得到圖像s17(fτ,fη)；

對所述圖像s17(fτ,fη)按行進行脈沖壓縮和二次距離壓縮誤差校正，得到圖像s18(fτ,fη)；

對所述圖像s18(fτ,fη)按行做快速逆傅里葉變換，得到圖像s19(τ,fη)；

對所述圖像s19(τ,fη)補償殘余誤差和進行方位散焦處理，得到圖像s20(τ,fη)；

對所述圖像s20(τ,fη)按列做快速逆傅里葉變換，得到所述正頻率方位模糊圖像。

可選地，所述判斷模塊包括：

第一判斷子模塊，用于判斷所述模糊圖像的數(shù)據(jù)點的絕對值是否大于所述星載合成孔徑雷達圖像中對應的數(shù)據(jù)點的絕對值的預定倍數(shù)，其中，所述預定倍數(shù)大于0且小于1；

第二判斷子模塊，用于判斷所述方位向時間是否滿足預設的閾值條件；

確定子模塊，用于當所述模糊圖像的數(shù)據(jù)點的絕對值大于所述星載合成孔徑雷達圖像中對應的數(shù)據(jù)點的絕對值的預定倍數(shù)、且所述方位向時間滿足所述閾值條件時，確定所述星載合成孔徑雷達圖像中的數(shù)據(jù)點滿足所述預定抑制條件。

可選地，所述模糊圖像包括負頻率方位模糊圖像；

所述閾值條件包括：

其中，η為所述方位向時間；fp為脈沖重復頻率；fr為多普勒調(diào)頻率；和/或

所述模糊圖像包括正頻率方位模糊圖像；

所述閾值條件包括：

其中，η為所述方位向時間；na為圖像矩陣總行數(shù)；fp為脈沖重復頻率；fr為方位向時間。

可選地，所述模糊抑制模塊用于將所述星載合成孔徑雷達圖像中、滿足所述預定抑制條件的數(shù)據(jù)點的絕對值，與所述模糊圖像中對應的數(shù)據(jù)點的絕對值做差值，得到經(jīng)模糊抑制后的圖像。

通過上述技術(shù)方案，對星載合成孔徑雷達圖像進行方位模糊重建，得到模糊圖像，并利用模糊圖像，對星載合成孔徑雷達圖像中、滿足預定抑制條件的數(shù)據(jù)點進行模糊抑制。這樣，能夠有效去除星載合成孔徑雷達圖像中的方位模糊，提高星載合成孔徑雷達圖像的清晰度。

本公開的其他特征和優(yōu)點將在隨后的具體實施方式部分予以詳細說明。

附圖說明

附圖是用來提供對本公開的進一步理解，并且構(gòu)成說明書的一部分，與下面的具體實施方式一起用于解釋本公開，但并不構(gòu)成對本公開的限制。在附圖中：

圖1是根據(jù)一示例性實施例示出的一種星載合成孔徑雷達圖像的模糊抑制方法的流程圖。

圖2是根據(jù)另一示例性實施例示出的一種星載合成孔徑雷達圖像的模糊抑制方法的流程圖。

圖3a是根據(jù)一示例性實施例示出的一種星載合成孔徑雷達圖像的模糊抑制裝置的框圖。

圖3b是根據(jù)另一示例性實施例示出的一種星載合成孔徑雷達圖像的模糊抑制裝置中的判斷模塊的框圖。

具體實施方式

以下結(jié)合附圖對本公開的具體實施方式進行詳細說明。應當理解的是，此處所描述的具體實施方式僅用于說明和解釋本公開，并不用于限制本公開。

請參考圖1，圖1是一示例性實施例提供的星載合成孔徑雷達圖像的模糊抑制方法的流程圖，如圖1所示，該方法包括s11、s12以及s13。下面依次對各個步驟進行說明。

在s11中，對星載合成孔徑雷達圖像進行方位模糊重建，確定出星載合成孔徑雷達圖像的模糊圖像。

本公開中，對星載合成孔徑雷達圖像進行方位模糊重建，重建模糊的幅度響應和相位歷史，得到與星載合成孔徑雷達圖像對應的模糊圖像。

此外，模糊圖像包括負頻率方位模糊圖像和/或正頻率方位模糊圖像。

首先，對確定負頻率方位模糊圖像的方法進行說明：

對星載合成孔徑雷達圖像s0(τ,η)按列做快速傅里葉變換，得到圖像s1(τ,fη)：

s1(τ,fη)＝ffta[s0(τ,η)]

其中，ffta為按列做快速傅里葉變換，j為列索引，nr為圖像矩陣總列數(shù)，fs為信號采樣率，i為行索引，na為圖像矩陣總行數(shù)，fp為脈沖重復頻率τ為距離向時間，η為方位向時間，fη為方位向頻率；

對圖像s1(τ,fη)按列進行方位相位補償和天線方向圖補償，每一列數(shù)據(jù)乘以傳遞函數(shù)h1，得到圖像s2(τ,fη)：

s2(τ,fη)＝h1·s1(τ,fη)

其中，f0為信號載頻，r0為斜距，c為光速，rref為場景中心斜距，λ為波長，vr為等效速度，fηc為多普勒中心頻率，d(fη-fp)為徙動因子，wa(fη-fηc-fp)為衛(wèi)星方位向天線方向圖的頻域表達式，wa(τ)為方向圖的時域表達式；

對圖像s2(τ,fη)按行做快速傅里葉變換，得到圖像s3(fτ,fη)：

s3(fτ,fη)＝fftr[s2(τ,fη)]

其中，fftr為按行做快速傅里葉變換；

對圖像s3(fτ,fη)按行進行一致模糊距離徙動校正和方位向調(diào)制，每一行數(shù)據(jù)乘以傳遞函數(shù)h2，得到圖像s4(fτ,fη)：

其中，q1(fη)為第一散焦因子，kr為信號調(diào)頻率，d(fηc)為徙動因子，fτ為距離向頻率；

對圖像s4(fτ,fη)按行做快速逆傅里葉變換，得到圖像s5(τ,fη)：

s5(τ,fη)＝ifftr[s4(fτ,fη)]

其中，ifftr為按行做快速逆傅里葉變換；

對圖像s5(τ,fη)按行進行補余模糊距離徙動校正，每一行數(shù)據(jù)乘以校正函數(shù)h3，得到圖像s6(τ,fη)：

其中，knew為信號調(diào)頻率，α為線性變標因子；

對圖像s6(τ,fη)按行做快速傅里葉變換，得到圖像s7(fτ,fη)：

s7(fτ,fη)＝fftr[s6(τ,fη)]

對圖像s7(fτ,fη)按行進行脈沖壓縮和二次距離壓縮誤差校正，每一行數(shù)據(jù)乘以校正函數(shù)h4，得到圖像s8(fτ,fη)：

其中，q2(fη)為第二散焦因子，d(fη)為徙動因子，km(fη)為距離向調(diào)頻率，km(fη-fp)為距離向調(diào)頻率；

對圖像s8(fτ,fη)按行做快速逆傅里葉變換，得到圖像s9(τ,fη)：

s9(τ,fη)＝ifftr[s8(fτ,fη)]

對圖像s9(τ,fη)補償殘余誤差和進行方位散焦處理，矩陣中每一個點均需要乘以對應的補償函數(shù)h5，得到圖像s10(τ,fη)：

對圖像s10(τ,fη)按列做快速逆傅里葉變換，得到負頻率方位模糊圖像：

s-1(τ,η)＝iffta[s10(τ,fη)]

其中，s-1(τ,η)為所述負頻率方位模糊圖像，iffta為按列做快速逆傅里葉變換。

其次，對確定正頻率方位模糊圖像的方法進行說明：

正頻率方位模糊圖像也可以按照上述步驟進行，只需要將上述公式中所有的-fp替換成+fp，即可得到方位模糊重建之后的正頻率方位模糊圖像s+1(τ,η)，本實施例中不做贅述。

在s12中，根據(jù)與星載合成孔徑雷達圖像中的數(shù)據(jù)點對應的模糊圖像的數(shù)據(jù)點，判斷星載合成孔徑雷達圖像中的數(shù)據(jù)點是否滿足預定抑制條件。

在本公開中，根據(jù)與星載合成孔徑雷達圖像的數(shù)據(jù)點對應的模糊圖像的數(shù)據(jù)點，判斷星載合成孔徑雷達圖像的數(shù)據(jù)點是否滿足預定抑制條件，也就是判斷該數(shù)據(jù)點是否為模糊點，如果該數(shù)據(jù)點是模糊點，則對該數(shù)據(jù)點進行模糊抑制。

圖2是根據(jù)另一示例性實施例示出的一種星載合成孔徑雷達圖像的模糊抑制方法的流程圖。如圖2所示，上述s12可以包括以下步驟。

在s121中，判斷模糊圖像的數(shù)據(jù)點的絕對值是否大于星載合成孔徑雷達圖像中對應的數(shù)據(jù)點的絕對值的預定倍數(shù)。

在本公開中，該預定倍數(shù)大于0且小于1，它可以例如是0.2。

示例地，當模糊圖像為負頻率方位模糊圖像時，當模糊圖像滿足以下等式(1)時，可以判定模糊圖像的數(shù)據(jù)點的絕對值大于星載合成孔徑雷達圖像中對應的數(shù)據(jù)點的絕對值的預定倍數(shù)：

|s-1(τ,η)|＞m×|s0(τ,η)|(1)

其中，s-1(τ,η)為所述負頻率方位模糊圖像，s0(τ,η)為星載合成孔徑雷達圖像，m為所述預定倍數(shù)。

又示例地，當模糊圖像為正頻率方位模糊圖像時，當模糊圖像滿足以下等式(2)時可以判定模糊圖像的數(shù)據(jù)點的絕對值大于星載合成孔徑雷達圖像中對應的數(shù)據(jù)點的絕對值的預定倍數(shù)：

|s+1(τ,η)|＞m×|s0(τ,η)|(2)

其中，s+1(τ,η)為所述正頻率方位模糊圖像，s0(τ,η)為星載合成孔徑雷達圖像，m為所述預定倍數(shù)。

在s122中，判斷方位向時間是否滿足預設的閾值條件。

在一種實施方式中，當模糊圖像為負頻率方位模糊圖像時，該閾值條件可以是其中，η為所述方位向時間，fp為脈沖重復頻率，fr為多普勒調(diào)頻率。

在另一種實施方式中，當模糊圖像為正頻率方位模糊圖像時，該閾值條件可以是其中，η為所述方位向時間，na為圖像矩陣總行數(shù)，fp為脈沖重復頻率，fr為多普勒調(diào)頻率。

另外，需要說明的是，上述s122可以在上述s121之前執(zhí)行，也可以與s121同時執(zhí)行，還可以在s121之后執(zhí)行，s121和s122之間的執(zhí)行順序在此不作具體限定。

在s123中，當模糊圖像的數(shù)據(jù)點的絕對值大于星載合成孔徑雷達圖像中對應的數(shù)據(jù)點的絕對值的預定倍數(shù)、且方位向時間滿足閾值條件時，確定星載合成孔徑雷達圖像中的數(shù)據(jù)點滿足預定抑制條件。

示例地，當模糊圖像為負頻率方位模糊圖像時，當模糊圖像滿足以下等式(3)時，可以確定星載合成孔徑雷達圖像中的數(shù)據(jù)點滿足預定抑制條件：

|s-1(τ,η)|＞m×|s0(τ,η)|，且

當星載合成孔徑雷達圖像中的數(shù)據(jù)點滿足預定抑制條件時，可以用邏輯值1來表示；當星載合成孔徑雷達圖像中的數(shù)據(jù)點不滿足預定抑制條件時，可以用邏輯值0來表示。此外，還可以將該判定結(jié)果保存在負模糊判定矩陣中。

示例地，可以通過以下等式(4)來建立負模糊判定矩陣：

其中，b-1(τ,η)為所述負模糊判定矩陣。

又示例地，當模糊圖像為正頻率方位模糊圖像時，當模糊圖像滿足以下等式(5)時，可以確定星載合成孔徑雷達圖像中的數(shù)據(jù)點滿足預定抑制條件：

|s+1(τ,η)|＞m×|s0(τ,η)|，且

同樣地，可以通過以下等式(6)來建立正模糊判定矩陣：

其中，b+1(τ,η)為所述正模糊判定矩陣。

在s13中，當星載合成孔徑雷達圖像中的數(shù)據(jù)點滿足預定抑制條件時，利用模糊圖像，對星載合成孔徑雷達圖像中、滿足預定抑制條件的數(shù)據(jù)點進行模糊抑制。

在一種實施方式中，可以通過將星載合成孔徑雷達圖像中、滿足預定抑制條件的數(shù)據(jù)點的絕對值，與模糊圖像中對應的數(shù)據(jù)點的絕對值做差值，來得到經(jīng)模糊抑制后的圖像。

示例地，可以通過以下等式(7)來得到經(jīng)模糊抑制后的圖像：

其中，sfinal(τ,η)為所述經(jīng)模糊抑制后的圖像。

通過上述技術(shù)方案，對星載合成孔徑雷達圖像進行方位模糊重建，得到模糊圖像，并利用模糊圖像，對星載合成孔徑雷達圖像中、滿足預定抑制條件的數(shù)據(jù)點進行模糊抑制。這樣，能夠有效去除星載合成孔徑雷達圖像中的方位模糊，提高星載合成孔徑雷達圖像的清晰度。

圖3a是根據(jù)一示例性實施例示出的一種星載合成孔徑雷達圖像的模糊抑制裝置。參照圖3a，該裝置30可以包括：方位模糊重建模塊301，用于對星載合成孔徑雷達圖像進行方位模糊重建，確定出所述星載合成孔徑雷達圖像的模糊圖像；判斷模塊302，用于根據(jù)與所述星載合成孔徑雷達圖像中的數(shù)據(jù)點對應的所述模糊圖像的數(shù)據(jù)點，判斷所述星載合成孔徑雷達圖像中的數(shù)據(jù)點是否滿足預定抑制條件；模糊抑制模塊303，用于當所述星載合成孔徑雷達圖像中的數(shù)據(jù)點滿足所述預定抑制條件時，利用所述模糊圖像，對所述星載合成孔徑雷達圖像中、滿足所述預定抑制條件的數(shù)據(jù)點進行模糊抑制。

可選地，所述模糊圖像包括負頻率方位模糊圖像；所述方位模糊重建模塊301用于：對所述星載合成孔徑雷達圖像s0(τ,η)按列做快速傅里葉變換，得到圖像s1(τ,fη)，其中，j為列索引，nr為圖像矩陣總列數(shù)，fs為信號采樣率，i為行索引，na為圖像矩陣總行數(shù)，fp為脈沖重復頻率，τ為距離向時間，η為方位向時間，fη為方位向頻率；

對所述圖像s1(τ,fη)按列進行方位相位補償和天線方向圖補償，得到圖像s2(τ,fη)：

s2(τ,fη)＝h1·s1(τ,fη)

其中，h1為傳遞函數(shù)，f0為信號載頻，r0為斜距，c為光速，rref為場景中心斜距，λ為波長，vr為等效速度，fηc為多普勒中心頻率，d(fη-fp)為徙動因子，wa(fη-fηc-fp)為衛(wèi)星方位向天線方向圖的頻域表達式，wa(τ)為方向圖的時域表達式；

對所述圖像s2(τ,fη)按行做快速傅里葉變換，得到圖像s3(fτ,fη)；

對所述圖像s3(fτ,fη)按行進行一致模糊距離徙動校正和方位向調(diào)制，得到圖像s4(fτ,fη)；

對所述圖像s4(fτ,fη)按行做快速逆傅里葉變換，得到圖像s5(τ,fη)；

對所述圖像s5(τ,fη)按行進行補余模糊距離徙動校正，得到圖像s6(τ,fη)；

對所述圖像s6(τ,fη)按行做快速傅里葉變換，得到圖像s7(fτ,fη)；

對所述圖像s7(fτ,fη)按行進行脈沖壓縮和二次距離壓縮誤差校正，得到圖像s8(fτ,fη)；

對所述圖像s8(fτ,fη)按行做快速逆傅里葉變換，得到圖像s9(τ,fη)；

對所述圖像s9(τ,fη)補償殘余誤差和進行方位散焦處理，得到圖像s10(τ,fη)；

對所述圖像s10(τ,fη)按列做快速逆傅里葉變換，得到所述負頻率方位模糊圖像。

可選地，所述模糊圖像包括正頻率方位模糊圖像；

所述方位模糊重建模塊301用于：

對所述星載合成孔徑雷達圖像s0(τ,η)按列做快速傅里葉變換，得到圖像s11(τ,fη)，其中，j為列索引，nr為圖像矩陣總列數(shù)，fs為信號采樣率，i為行索引，na為圖像矩陣總行數(shù)，fp為脈沖重復頻率，τ為距離向時間，η為方位向時間，fη為方位向頻率；

對所述圖像s11(τ,fη)按列進行方位相位補償和天線方向圖補償，得到圖像s12(τ,fη)：

s12(τ,fη)＝h1·s11(τ,fη)

其中，h1為傳遞函數(shù)，f0為信號載頻，r0為斜距，c為光速，rref為場景中心斜距，λ為波長，vr為等效速度，fηc為多普勒中心頻率，d(fη+fp)為徙動因子，wa(fη-fηc+fp)為衛(wèi)星方位向天線方向圖的頻域表達式，wa(τ)為方向圖的時域表達式；

對所述圖像s12(τ,fη)按行做快速傅里葉變換，得到圖像s13(fτ,fη)；

對所述圖像s13(fτ,fη)按行進行一致模糊距離徙動校正和方位向調(diào)制，得到圖像s14(fτ,fη)；

對所述圖像s14(fτ,fη)按行做快速逆傅里葉變換，得到圖像s15(τ,fη)；

對所述圖像s15(τ,fη)按行進行補余模糊距離徙動校正，得到圖像s16(τ,fη)；

對所述圖像s16(τ,fη)按行做快速傅里葉變換，得到圖像s17(fτ,fη)；

對所述圖像s17(fτ,fη)按行進行脈沖壓縮和二次距離壓縮誤差校正，得到圖像s18(fτ,fη)；

對所述圖像s18(fτ,fη)按行做快速逆傅里葉變換，得到圖像s19(τ,fη)；

對所述圖像s19(τ,fη)補償殘余誤差和進行方位散焦處理，得到圖像s20(τ,fη)；

對所述圖像s20(τ,fη)按列做快速逆傅里葉變換，得到所述正頻率方位模糊圖像。

如圖3b所示，所述判斷模塊302可以包括：第一判斷子模塊3021，用于判斷所述模糊圖像的數(shù)據(jù)點的絕對值是否大于所述星載合成孔徑雷達圖像中對應的數(shù)據(jù)點的絕對值的預定倍數(shù)，其中，所述預定倍數(shù)大于0且小于1；第二判斷子模塊3022，用于判斷所述方位向時間是否滿足預設的閾值條件；確定子模塊3023，用于當所述模糊圖像的數(shù)據(jù)點的絕對值大于所述星載合成孔徑雷達圖像中對應的數(shù)據(jù)點的絕對值的預定倍數(shù)、且所述方位向時間滿足所述閾值條件時，確定所述星載合成孔徑雷達圖像中的數(shù)據(jù)點滿足所述預定抑制條件。

可選地，所述模糊圖像包括負頻率方位模糊圖像；所述閾值條件包括：

其中，η為所述方位向時間；fp為脈沖重復頻率；fr為多普勒調(diào)頻率；和/或

所述模糊圖像包括正頻率方位模糊圖像；所述閾值條件包括：

其中，η為所述方位向時間；na為圖像矩陣總行數(shù)；fp為脈沖重復頻率；fr為多普勒調(diào)頻率。

可選地，所述模糊抑制模塊303用于將所述星載合成孔徑雷達圖像中、滿足所述預定抑制條件的數(shù)據(jù)點的絕對值，與所述模糊圖像中對應的數(shù)據(jù)點的絕對值做差值，得到經(jīng)模糊抑制后的圖像。

以上結(jié)合附圖詳細描述了本公開的優(yōu)選實施方式，但是，本公開并不限于上述實施方式中的具體細節(jié)，在本公開的技術(shù)構(gòu)思范圍內(nèi)，可以對本公開的技術(shù)方案進行多種簡單變型，這些簡單變型均屬于本公開的保護范圍。

另外需要說明的是，在上述具體實施方式中所描述的各個具體技術(shù)特征，在不矛盾的情況下，可以通過任何合適的方式進行組合。為了避免不必要的重復，本公開對各種可能的組合方式不再另行說明。

此外，本公開的各種不同的實施方式之間也可以進行任意組合，只要其不違背本公開的思想，其同樣應當視為本公開所公開的內(nèi)容。