function?[?img_resnChange]?=?changedetect（?ABmethodthreshold）
%?變化檢測函數
%?輸入：變化前后圖像A、B，以及檢測方法名method對應閾值threshold
%?輸出：結果差異二值圖像img_res，差異像素數nChange
img1=rgb2gray（A）;
img2=rgb2gray（B）;
[rowcol]?=?size（img1）;
img3=ones（rowcol）;
%?比值法
if?strcmp（method‘ratio‘）
????nChange?=?0;
????for?i=1:row
????????for?j=1:col
???????????if（img1（ij）==0&&img2（ij）~=0）
???????????????if（abs（img1（ij）/img2（ij）-1）>threshold）%閾值
???????????????????img3（ij）=0;
???????????????????nChange?=?nChange+1;
???????????????end
???????????elseif（img2（ij）==0&&img1（ij）~=0）
???????????????if（abs（img2（ij）/img1（ij）-1）>threshold）%閾值
???????????????????img3（ij）=0;
???????????????????nChange?=?nChange+1;
???????????????end
???????????else
???????????????if（abs（img1（ij）/img2（ij）-1）>threshold）%閾值
???????????????????img3（ij）=0;
???????????????????nChange?=?nChange+1;
???????????????end
???????????end??
????????end
????end
end
%?歸一化法
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（ij）=（img1（ij）-min1）/（max1-min1）;
????????????img2（ij）=（img2（ij）-min2）/（max2-min2）;
????????end
????end?
????for?i=1:row
????????for?j=1:col
????????????if（abs（img1（ij）-img2（ij））>threshold）
????????????????img3（ij）=0;
????????????????nChange?=?nChange+1;
????????????end
????????end
????end
end
%?相關系數法
if?strcmp（method‘corrcoef‘）
????nChange?=?0;
????temp1=zeros（33）;
????temp2=zeros（33）;
????img1=double（img1）;
????img2=double（img2）;
????img3=double（img3）;
????%?計算5*5大小的窗口的相關系數
????for?i=2:row-1
????????for?j=2:col-1
????????????temp1=[img1（i-1j-1）img1（i-1j）img1（i-1j+1）;img1（ij-1）img1（ij）img1（ij+1）;img1（i+1j-1）img1（i+1j）img1（i+1j+1）];
????????????temp2=[img2（i-1j-1）img2（i-1j）img2（i-1j+1）;img2（ij-1）img2（ij）img2（ij+1）;img2（i+1j-1）img2（i+1j）img2（i+1j+1）];
????????????corr=corrcoef（temp1temp2）;
???????????if（abs（corr（12））???????????????img3（ij）=0;
???????????????nChange?=?nChange+1;
???????????end
????????end
????end
end
img_res?=?img3;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2014-12-06?00:28??遙感圖像變化檢測\
?????文件?????4826954??2014-11-03?12:22??遙感圖像變化檢測\0804_1w.bmp
?????文件?????4826954??2014-11-03?12:22??遙感圖像變化檢測\0902_1w.bmp
?????文件????????2339??2014-12-03?19:57??遙感圖像變化檢測\changedetect.m
?????文件???????11382??2014-12-03?19:36??遙感圖像變化檢測\ImgChangeDetect.fig
?????文件????????7974??2014-12-03?19:56??遙感圖像變化檢測\ImgChangeDetect.m