[发明专利]基于变差函数的高分SAR新增建设用地提取软件

权利要求书

1.一种基于高分SAR新增建设用地检测方法，其步骤为：

第一步、数据预处理

1-1)输入不同时相的高分辨率SAR图像

1-2)检测输入图像的尺寸是否相同

若相同进行滤波处理，若不同则终止程序

1-3)进行基于结构检测的相干斑抑制滤波处理，首先进行强点目标检测，对于强点目标滤波算法保持其原值不作处理，对于非强点目标初始窗口大小为N，计算局域方差系数C_ij，设置阈值Cu，当C_ij≤Cu时，认定窗口区域为同质区域，进一步寻找最大同质区域大小为N×N矩形，然后进行区域增长，选择最终参与滤波的像素，进行均值滤波；当C_ij＞Cu时，认定窗口区域为非同质区域，进行边缘与线检测，对于微纹理区选择N×N矩形窗口计算滤波参数b进行MMSE滤波，对于边缘选择基于边缘的滤波窗口计算滤波参数b进行MMSE滤波，同样对于线选择基于线的滤波窗口计算滤波参数b后进行MMSE滤波，最终输出滤波后的影像；

第二步、改进变差函数提取新增建设用地

2-1)对滤波后的图像计算变差函数生成纹理特征图

变差函数

其中，ω为二维模板，采用0°，45°，90°，135°四个方向，在计算像素点的变差函数值时，计算像素点(x,y)与模板四个方向上间距h的点(x+h,y)的半方差值，将此值作为点(x,y)的半方差值；然后取窗口内所有像素的半方差值的中值代替均值作为变差函数值；

2-2)不同时相的纹理特征图作比值

以新时相的图像除以老时相的图像，老时相图像可能存在为0的数据点，在相除之前将所有为0的数据点赋值为0.000001，相除后的结果进行归一化至0-255之间；

2-3)确定阈值提取新增建设用地

根据经验取得阈值a＝200，当归一化后的数据点值大于a时，认为此数据点为新增建设用地的数据点；

2-4)对提取的结果进行后处理，得到最终的结果

搜索初步提取结果中所有的联通区，通过膨胀腐蚀来连接破碎的新增建设用地边界，并通过填洞获取完整的建筑区，最后利用删除去掉非新增建设用地的小面积区域；

2-5)精度评价

可以通过输入真值进行提取结果的精度评价；

检测率：DR＝TP/(TP+FN)*100

虚警率：Far＝FP/(TP+FP)*100

漏警率：MAR＝FN/(TP+FN)*100

其中TP为正确提取的像元数量，FN为真值中是新增建设用地但程序未提取的像元数量，FP为真值中是非建设用地但程序错误提取的像元数量。