[发明专利]一种SAR图像变化检测差异图生成方法

权利要求书

1.一种SAR图像变化检测差异图生成方法，包括以下步骤：

步骤1.通过合成孔径雷达获取两幅不同时相同样大小的SAR图像，对这两幅SAR图像进行预处理，预处理不改变图像的大小，获得预处理后的适合变化检测使用的SAR图像X¹和X²；

步骤2.分别对图像X¹和图像X²进行全图逐像素遍历，选取第i个位置的像素点X_i作为待恢复像素点，i＝1,2,…,m×n，这里X_i泛指两幅图像X¹和X²第i个位置的像素点，其中m，n分别是步骤1中获取的两幅不同时相同样大小的SAR图像的长度和宽度，选取以X_i为中心的正方形W_i^r作为搜索窗口，选取以X_i为中心的正方形图块作为X_i的相似性测量图块，的大小应小于搜索窗口W_i^r的大小；

步骤3.在所述两幅图像X¹和X²中，对于以X_i为中心的搜索窗口W_i^r中其他像素点X_j,X_j∈W_i^r，选取以X_j为中心的正方形图块作为像素点X_j的相似性测量图块，正方形图块与正方形图块的大小相等，计算以像素点X_i为中心的图块和以X_j为中心的图块这两个相似性测量图块对应的矩阵对应位置处的像素点灰度比值r_i,j的概率密度分布函数p(r_i,j)；

步骤4.对概率密度分布函数p(r_i,j)进行高斯核函数映射得到以像素点X_i为中心的图块和以X_j为中心的图块这两个相似性测量图块对应的矩阵对应位置处的像素点比值相似度R(i,j)；

步骤5.在两幅图像X¹和X²中，分别计算待恢复像素点X_i与其搜索窗口W_i^r中所有其他像素点X_j的相似性权值ω(i,j)；

步骤6.对于像素点X_i，用相似性权值ω(i,j)与搜索窗口W_i^r中第j位置的像素点X_j的灰度值加权相乘，对加权乘积求和，得到两幅图像X¹和X²第i个位置的待恢复像素点X_i的基于比值相似度的非局部均值

步骤7.使用步骤6中得到的基于比值相似度的非局部均值NL[X_i]，得到非局部均值比差异图中对于像素点X_i的像素点，XiNLMRD=1-min{NL[Xi1]NL[Xi2],NL[Xi2]NL[Xi1]}]]>

其中，为SAR图像X¹和X²的差异图中第i位置像素点，与分别表示SAR图像X¹与X²的第i位置像素点，与表示SAR图像X¹和X²的第i位置像素点灰度值的非局部均值。