基于变化检测的SAR图像稀疏去噪方法.pdf

(19)中华人民共和国国家知识产权局*CN103473755A*(12)发明专利申请(10)申请公布号(10)申请公布号CNCN103473755103473755A(43)申请公布日2013.12.25(21)申请号201310404907.X(22)申请日2013.09.07(71)申请人西安电子科技大学地址710071陕西省西安市太白南路2号(72)发明人王桂婷焦李成接道伟杨淑媛马文萍马晶晶田小林王爽公茂果(74)专利代理机构陕西电子工业专利中心61205代理人田文英王品华(51)Int.Cl.G06T5/50(2006.01)权权利要求书4页利要求书4页说明书9页说明书9页附图4页附图4页(54)发明名称基于变化检测的SAR图像稀疏去噪方法(57)摘要本发明公开了一种基于变化检测的SAR图像稀疏去噪方法，其实现步骤为：（1）输入图像；（2）得到变化区域图像；（3）设定变化类；（4）得到两分类二值图；（5）得到两区域图像；（6）不变化区域去噪；（7）估计噪声标准差；（8）获得合成图；（9）稀疏字典去噪；（10）输出结果。本发明是针对SAR图像变化检测的去噪预处理，具有对变化区域一致性保持较好，能很好的抑制不变化区域的噪声，同时保留SAR图像特别是对后期变化检测影响较大的变化区域的纹理、轮廓、边缘细节信息，消除图像去噪后产生的划痕，减少因噪声抑制产生的伪信息，提高变化检测的精度的优点。CN103473755ACN103475ACN103473755A权利要求书1/4页1.一种基于变化检测的SAR图像稀疏去噪方法，包括如下步骤：（1）输入图像：1a）选取两幅卫星对同一地点、不同时间产生的已配准的大小相同的单极化SAR图像；1b）读入第一幅单极化SAR图像矩阵中的时相图像1；1c）读入第二幅单极化SAR图像矩阵中的时相图像2；（2）得到变化区域图像：2a）对两时相SAR图像作差，并将作差运算结果取绝对值，得到一幅差异图像；2b）对差异图像进行窗口大小为3×3模板的均值滤波，得到滤波后的差异图像；2c）按照下式计算，得到滤波后的差异图像中任意一个像素点的平均能量，其中，E表示任意一个像素点的平均能量，q表示以像素点为中心，大小为3×3像素的邻域内的像素序号，q=1，2，…，9，B表示以像素点为中心，大小为3×3像素的邻域，gq表示邻域B内的第q个像素的灰度值；2d）用每个差异图像中的像素点的平均能量，替代每个差异图像中的像素点的灰度值，得到一幅平均能量图；2e）计算平均能量图的灰度均值，将该均值设定为阈值，将大于等于阈值的灰度值赋为1，小于阈值的灰度值赋为0，得到一幅二值化的平均能量图；2f）对二值化的平均能量图，用3×3的结构元素进行形态学腐蚀操作，得到变化检测图像，该变化检测图像中像素值为1的区域表示检测出的两时相图像的变化区域；2g）将时相图像1和时相图像2中与两时相图像变化区域对应的空间位置之外的像素灰度值分别设置为0，得到时相图像1的初始变化区域图像和时相图像2的初始变化区域图像；（3）设定变化类：3a）用时相图像1的初始变化区域图像减去时相图像2的初始变化区域图像，得到一个初始变化区域差异图像；3b）将初始变化区域差异图像中像素值小于等于设定的阈值T1的像素点，设定为非变化类；将初始变化区域差异图像中像素值大于阈值T1的像素点设定为变化类；3c）分别将与初始变化区域差异图像中设定为变化类的像素点的空间位置对应的时相图像1初始变化区域和时相图像2初始变化区域中的像素值设置为1；（4）得到两分类二值图：4a）对初始变化区域差异图像中每一个属于非变化类的像素点，在以该像素点为中心的3×3大小的邻域内，统计大于0的像素个数；4b）比较邻域内大于0的像素个数与阈值T2、T3的大小，进行以下操作，得到时相图像1的两分类二值图和时相图像2的两分类二值图：若邻域内大于0的像素个数大于等于阈值T2，则将与满足此条件的非变化类的像素空间位置对应的时相图像1的初始变化区域图像的像素点的值设置为1；若邻域内大于0的像素个数小于阈值T3，则将与满足此条件的非变化类的像素空间位2CN103473755A权利要求书2/4页置对应的时相图像2的初始变化区域图像的像素点的值设置为1；若邻域内大于0的像素个数小于阈值T2并大于等于阈值T3，则将与满足此条件的非变化类的像素空间位置对应的时相图像1的初始变化区域图像和时相图像2的初始变化区域图像的像素点的值均设置为1；将时相图像1的初始变化区域图像和时相图像2的初始变化区域图像中没有标记的像素点的值设置为0；(5)获得两区域图像：5a）将时相图像1的两分类二值图与时相图像1的图像的空间位置对应的像素值逐像素点相乘，得到时相图像1的变化区域图像；用时相图像1减去时相1图像的变化区域图像，得到时相图像1的不变