- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
SAR图像变化检测
题 目 SAR图像变化检测
学 院 电子工程学院
专 业 智能科学与技术
学 号
学生姓名
教师姓名
SAR图像变化检测
引言
图像检测,是利用计算机对图像进行处理、分析和理解,以识别各种不同模式的目标和对像的技术。地理学中指将遥感图像进行分类的技术。 图形刺激作用于感觉器官,人们辨认出它是经验过的某一图形的过程,也叫图像再认。在图像识别中,既要有当时进入感官的信息,也要有记忆中存储的信息。只有通过存储的信息与当前的信息进行比较的加工过程,才能实现对图像的再认。图像识别是人工智能的一个重要领域。
2.算法原理简介
首先对两幅实验用SAR 图像进行滤波处理,消除噪声,尽可能的减少噪点对变化检测结果的影响。然后利用图像比值法生成初步的新像素矩阵。假设P1(i,j)是第一幅实验用SAR 图像对应像素点(i,j)的灰度,P2(i,j)是第二幅实验用图对应像素点(i,j)的灰度。对应两幅图像,变化的部分像素差异较大,所以比值P1(i,j)/P2(i,j)距离1 就比较远;而对应不变的点,像素虽然不至于完全一致,但是比值会在1 的附近。这个特性就为此方法的实现奠定了良好的基础。
然后对上一步得到的像素矩阵进行处理,例如进行取反操作(这个要看实际需要),乘以某一恰当的系数将像素值的差距拉大。便可得到一幅SAR 图像变化图像。最后将实验所得的结果跟参考进行比对,得到差异,并统计漏检和误检的像素点数量。
3.具体算法步骤:
步骤1:读入变化前后的实验用用SAR图像
步骤2:将两幅图像的像素矩阵相比,得到新的像素矩阵
步骤3:对上步所得像素矩阵进行取反,将系数等简单的处理,得到变化检测的图像
步骤4: 对上步得到的新图像以及参考进行模糊C均值分类,对参考和新图像进行比较
步骤5: 对参考和实验结果进行相减操作,得到误差图像
步骤6: 统计漏测点和误测点的数量并输出
步骤7: 结束
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]);
[N,L]=size(I1);
I=(I1./I2).*256;%±è?μ′|àí?°à-′ó2??à
I=imcomplement(I);%è?·′2ù×÷
IM=I;
I11=fcm2(I);
IM=I3;
I22=fcm2(I3);
k=0;
[T,J]=cuowu(I11,I22,N,L);
k=T+J;
disp(sprintf(?á1????ì%d??????μ?£??ó?ì%d??????μ?,T,J));
disp(sprintf(?á1?í???12óD%d??????£????Dó?2????ó2???????êy?a%d,N*L,k));
k=k/(N*L)*100;
disp(sprintf(′?·?·¨ó?2???í???μ??ó2??ê%2.4f??°ù·?μ?,k));
I4=I11-I22; %?ó2?í???
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);
subplot(122)
imshow(I22);
title(reference_fcm2);
figure(3)
imshow(I4);
title(·?ààoó?á1?ó?2???μ??ó2?);
function IMMM=fcm2(IM) %·?ààêy
文档评论(0)