[发明专利]基于迭代最大类间方差的SAR图像变化检测方法

摘要

本发明公开了一种基于迭代最大类间方差的SAR图像变化检测方法，主要解决现有技术对于变化区域占整体区域的小部分情况下检测会产生过多虚警点的问题，其实现方案是：1.获取不同时相同一地区的SAR图像并进行预处理；2.生成对数比差异图；3.迭代使用最大类间方差法遍寻差异图查找变化元素与非变化元素相对均匀的区域；4.将这一区域的分割阈值应用到全图，得到分割结果图；5.对分割结果图进行中值滤波，得到最终变化检测结果图。本发明解决了直方图为单峰时最大类间方差法虚警过多的问题，提高了运行效率较高，可用于快速SAR图像变化检测，尤其是对于小面积变化区域的检测。

主权项

1.一种基于迭代最大类间方差的SAR图像变化检测方法，包括如下：(1)通过合成孔径雷达获取在不同时相拍摄的同一地区的两幅SAR图像，对这两幅SAR图像依次进行辐射校正、几何校正、图像配准的预处理，获得预处理后的适合变化检测使用的两幅SAR图像I1和I2；(2)计算第一幅SAR图像I1在坐标(i,j)处的灰度值I1(i,j)和第二幅SAR图像I2在坐标(i,j)处的灰度值I2(i,j)的比值，对比值运算的结果再进行一步对数运算，得到坐标(i,j)处的对数比值DI(i,j)：DI(i,j)＝log(I1(i,j)/I2(i,j))，遍历两幅SAR图像所有的像素点，得到对数比差异图：DI＝{DI(i,j)}；(3)使用最大类间方差法对对数比差异图DI进行初分割，得到初步分割阈值T1；(4)用阈值T1对对数比差异图DI进行二值分割得到初步分割图SI1；(5)设置大小为w*w的窗口对初步分割图SI1进行滑窗，统计窗口内零元素与非零元素的个数，直至小窗口内零元素与非零元素个数差小于等于3时，停止滑窗，得到满足条件的窗口中心元素(p,q)，其中w*w窗口大小取值为待处理图像大小的0.5％‑1.5％；(6)取对数比差异图DI中中心像素为(p,q)，窗口大小为w*w的区域P，使用最大类间方差法获取区域P的分割阈值T，用阈值T对对数比差异图DI进行二值分割，得到新的分割图SI；(7)计算区域P的分割阈值T与对数比差异图DI的预分割阈值T1的差值，若差值小于10‑4则认为阈值收敛，执行步骤(8)，否则，将区域P的分割阈值T的值赋给对数比差异图预分割阈值T1，将新的分割图SI的值赋给初步分割图SI1，返回步骤(5)；(8)对分割图SI进行3*3中值滤波，消除点状噪声对变化检测结果的影响，生成变化检测结果图。