[发明专利]一种基于异质性预矫正的小波域SAR图像去斑方法

权利要求书

1.一种基于异质性预矫正的小波域SAR图像去斑方法，其特征在于，所述方法包括如下步骤：

(1)小波变换，利用平稳小波变换对图像做J层分解，在每个尺度上得到一个低频子带和三个方向的高频子带；

(2)计算异质度，逐层计算每个小波系数所对应的多尺度局部变差系数MLCV，和该层的模态mode值，具体包括以下三个子步骤：

(2-1)对于第j层子带d＝LH,HL,HH表示子带方向、k表示空间位置，对于第j层子带LL表示子带方向、k表示空间位置,分别按照公式计算在实现过程中按3×3模板窗口的方式遍历图像子带；

(2-2)按照下述公式计算MLCV，

c(k)=E(wkLH)2+E(wkHL)2+E(wkHH)2E(wkLL),]]>

式中c(k)代表每个小波系数所对应的多尺度局部变差系数MLCV；

(2-3)统计c(k)的分布得到其直方图，并求出直方图中纵坐标对应的最大值，该最大值所对应的横坐标值即为模态mode；

(3)异质度分类，对小波系数进行分类，根据计算的MLCV值和模态mode值对每一层三个高频子带中的小波系数按如下公式进行异质度分类：

同质区域：0＜c(k)≤mode，

弱纹理区域：mode＜c(k)≤2mode，

强纹理区域：2mode＜c(k)≤5mode，

孤立散射点区域：c(k)＞5mode；

(4)预矫正处理，按照下列预矫正函数对每一层三个高频子带(LH,HL,HH)中的小波系数做预矫正处理：

pre-correcting(wk)={a·[11+exp[-λc(k)mode-gate]-0.5]+1}·wk,]]>

上式中：w_k表示预矫正前的小波系数，幅度因子α、指数因子λ及gate的取值范围用(α,λ,gate)表示，(α,λ,gate)在同质区域、弱纹理区域、强纹理区域的取值范围分别为(0.65-0.85,2-2.5,0.9-1.1)、(0.01-0.5,0.1-0.5,1.2-1.8)、(0.6-0.8,0.7-0.85,1.8-2.3)，孤立散射点区域的小波系数保留原有值不变，不进行预矫正处理；

(5)小波域去斑处理。