clear?%工作區清零
clc?%命令行清零
Path1=‘F:\代做\變化檢測\變化后.jpg‘;
Path2=‘F:\代做\變化檢測\變化前.jpg‘;
Data1=imread（Path1）;%讀取圖像的信息
figure（1）;imshow（Data1）;%畫出變化前圖片
Data2=imread（Path2）;%讀取圖像的信息
figure（2）;imshow（Data2）;%畫出變化后圖片

[rowcol]=size（Data1）;
data1=zeros（rowcol）;
data2=zeros（rowcol）;
data=zeros（rowcol）;
DATA=zeros（rowcol）;
for?i=1:1:row*col
data1（i）=Data1（i）;
data2（i）=Data2（i）;
data（i）=data1（i）-data2（i）;
end

for?i=1:row*col
????if?abs（data（i））<60?%注意這是進行判斷有效值范圍，如果有效范圍是-1到1，則改成max（bz）>-1即可
????????DATA（i）=0;
%?????elseif?（abs（data（i））>=20）&&（abs（data（i））<40）
%???????