【图像检测】基于计算机视觉实现遥感图像变换检测含Matlab源码（基于matlab的遥感图像处理）

【图像检测】基于计算机视觉实现遥感图像变换检测含Matlab源码（基于matlab的遥感图像处理）

1 简介

变化监测是从不同时期的遥感图像中,定量地分析和确定地物变化的特征的过程,其中涉及到变化的类型、分布状况及变化信息的描述,即需要确定变化前后的地物类型,界限和分析变化的属性.本文采用在遥感软件中进行图像配准后再在MATLAB下使用差值法进行变化监测,能够监测出较好的变化效果.

2 部分代码

function [ img_res,nChange] = changedetect( A,B,method,threshold)% 变化检测函数% 输入：变化前后图像A、B，以及检测方法名method,对应阈值threshold% 输出：结果差异二值图像img_res，差异像素数nChangeimg1=rgb2gray(A);img2=rgb2gray(B);[row,col] = size(img1);img3=ones(row,col);% 比值法if strcmp(method,'ratio') nChange = 0; for i=1:row for j=1:col if(img1(i,j)==0&&img2(i,j)~=0) if(abs(img1(i,j)/img2(i,j)-1)>threshold)%阈值 img3(i,j)=0; nChange = nChange+1; end elseif(img2(i,j)==0&&img1(i,j)~=0) if(abs(img2(i,j)/img1(i,j)-1)>threshold)%阈值 img3(i,j)=0; nChange = nChange+1; end else if(abs(img1(i,j)/img2(i,j)-1)>threshold)%阈值 img3(i,j)=0; nChange = nChange+1; end end end endend% 归一化法if strcmp(method,'norm') nChange = 0; img1=double(img1); img2=double(img2); img3=double(img3); max1=max(max(img1)); max2=max(max(img2)); min1=min(min(img1)); min2=min(min(img2)); for i=1:row for j=1:col img1(i,j)=(img1(i,j)-min1)/(max1-min1); img2(i,j)=(img2(i,j)-min2)/(max2-min2); end end for i=1:row for j=1:col if(abs(img1(i,j)-img2(i,j))>threshold) img3(i,j)=0; nChange = nChange+1; end end endend% 相关系数法if strcmp(method,'corrcoef') nChange = 0; temp1=zeros(3,3); temp2=zeros(3,3); img1=double(img1); img2=double(img2); img3=double(img3); % 计算5*5大小的窗口的相关系数 for i=2:row-1 for j=2:col-1 temp1=[img1(i-1,j-1),img1(i-1,j),img1(i-1,j+1);img1(i,j-1),img1(i,j),img1(i,j+1);img1(i+1,j-1),img1(i+1,j),img1(i+1,j+1)]; temp2=[img2(i-1,j-1),img2(i-1,j),img2(i-1,j+1);img2(i,j-1),img2(i,j),img2(i,j+1);img2(i+1,j-1),img2(i+1,j),img2(i+1,j+1)]; corr=corrcoef(temp1,temp2); if(abs(corr(1,2))

3 仿真结果

4 参考文献

[1]张孟君, 李晨钊, 舒红,等. 基于MATLAB的遥感图像变化检测[J]. 国土资源遥感, 2005(4):5.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。