SAR图像变化检测

精品文档 1欢迎下载 题题 目目 SAR 图像变化检测 学学 院院 电子工程学院 专专 业业 智能科学与技术 学学 号号 精品文档 2欢迎下载 学生姓名学生姓名 教师姓名教师姓名 SARSAR 图像变化检测图像变化检测 1 1 引言引言 图像检测 是利用计算机对图像进行处理 分析和理解 以识 别各种不同模式的目标和对像的技术 地理学中指将遥感图像 进行分类的技术 图形刺激作用于感觉器官 人们辨认出它 是经验过的某一图形的过程 也叫图像再认 在图像识别中 既 要有当时进入感官的信息 也要有记忆中存储的信息 只有通过 存储的信息与当前的信息进行比较的加工过程 才能实现对图 像的再认 图像识别是人工智能的一个重要领域 2 2 算法原理简介 算法原理简介 首先对两幅实验用SAR 图像进行滤波处理 消除噪声 尽可能 精品文档 3欢迎下载 的减少噪点对变化检测结果的影响 然后利用图像比值法生成 初步的新像素矩阵 假设P1 i j 是第一幅实验用SAR 图像对应 像素点 i j 的灰度 P2 i j 是第二幅实验用图对应像素 点 i j 的灰度 对应两幅图像 变化的部分像素差异较大 所 以比值P1 i j P2 i j 距离1 就比较远 而对应不变的点 像 素虽然不至于完全一致 但是比值会在1 的附近 这个特性就 为此方法的实现奠定了良好的基础 然后对上一步得到的像素矩阵进行处理 例如进行取反操作 这个要看实际需要 乘以某一恰当的系数将像素值的差距 拉大 便可得到一幅SAR 图像变化图像 最后将实验所得的结 果跟参考进行比对 得到差异 并统计漏检和误检的像素点数 量 3 3 具体算法步骤 具体算法步骤 步骤 1 读入变化前后的实验用用 SAR 图像 步骤步骤 2 2 将两幅图像的像素矩阵相比 得到新的像素矩阵 步骤步骤 3 3 对上步所得像素矩阵进行取反 将系数等简单的处理 得到变化检测的图像 步骤步骤 4 4 对上步得到的新图像以及参考进行模糊 C 均值分类 对 参考和新图像进行比较 步骤步骤 5 5 对参考和实验结果进行相减操作 得到误差图像 精品文档 4欢迎下载 步骤步骤 6 6 统计漏测点和误测点的数量并输出 步骤步骤 7 7 结束 4 4 算法程序算法程序 function Untitled2 clc clear all close all I1 imread 1999 04 bmp I2 imread 1999 05 bmp I3 imread 1999 05 bmp I1 rgb2gray I1 I2 rgb2gray I2 I3 rgb2gray I3 I1 medfilt2 I1 3 3 I2 medfilt2 I2 3 3 I3 medfilt2 I3 3 3 精品文档 5欢迎下载 N L size I1 I I1 I2 256 I imcomplement I IM I I11 fcm2 I IM I3 I22 fcm2 I3 k 0 T J cuowu I11 I22 N L k T J disp sprintf d d T J disp sprintf d d N L k k k N L 100 disp sprintf 2 4f k I4 I11 I22 精品文档 6欢迎下载 figure 1 subplot 221 imshow I1 title before change subplot 223 imshow I2 title after change subplot 222 imshow I3 title reference subplot 224 imshow I title result figure 2 subplot 121 imshow I11 title result fcm2 精品文档 7欢迎下载 subplot 122 imshow I22 title reference fcm2 figure 3 imshow I4 title function IMMM fcm2 IM C IM IX2 maxX maxY size IM IM double IM IMM cat 2 IM IM 2 cc1 8 cc2 230 IMM cat 3 IM IM ttFcm 0 精品文档 8欢迎下载 while ttFcm 30 ttFcm ttFcm 1 c1 repmat cc1 maxX maxY c2 repmat cc2 maxX maxY c cat 3 c1 c2 ree repmat 0 000001 maxX maxY ree1 cat 3 ree ree distance IMM c distance distance distance ree1 daoshu 1 distance daoshu2 daoshu 1 daoshu 2 u distance1 distance 1 daoshu2 u1 1 distance1 distance2 distance 2 daoshu2 u2 1 distance2 z 精品文档 9欢迎下载 ccc1 sum sum u1 u1 IM sum sum u1 u1 ccc2 sum sum u2 u2 IM sum sum u2 u2 tmpMatrix abs cc1 ccc1 cc1 abs cc2 ccc2 cc2 pp cat 3 u1 u2 for i 1 maxX for j 1 maxY if max pp i j u1 i j IX2 i j 1 elseif max pp i j u2 i j IX2 i j 2 end end end if max tmpMatrix 0 0001 break else 精品文档 10欢迎下载 cc1 ccc1 cc2 ccc2 end end for i 1 maxX for j 1 maxY if IX2 i j 2 IMMM i j 255 else IMMM i j 0 end end end end function T J cuowu I11 I22 N L T 0 J 0 精品文档 11欢迎下载 for i 1 N for j 1 L if I11 i j 0 elseif I11 i j 255 end end end end end 5 5 实验结果实验结果 精品文档 12欢迎下载 精品文档 13欢迎下载 6 6 实验结果分析实验结果分析 结果漏检 44545 个像素点 误检 49 个像素点 结果图像共有 90601 个像素 其中与参考误差像素个数为 44594 此方法与参考图像的误差率 49 2202 个百分点 7 7 参考资料参考资料 1 陈传波 金先级 数字图像处理 北京 机械工业出版社 2 李金基 焦李成 张向荣 杨咚咚 基于两时相图像联合分类的SAR 图