clc;close?all;clear?all;
p1=imread（‘3.jpg‘）%用于讀取圖片文件中的數據
p2=imread（‘4.jpg‘）%讀入圖片并分別賦值到p1?p2
dif1=imabsdiff（p1p2）%對p1p2做差值，找出不同點，并賦值dif1
dif=im2bw（dif10.1）%將dif1轉化為2值圖像
I=dif
imshow（I）
title（‘binary?image‘）;
[mn]=size（I）;
txt=get（gca‘title‘）;%獲得矩陣大小，gca返回屬性值
set（txt‘fontsize‘16）;
L=bwlabel（I）;
stats=regionprops（L‘all‘）;
set（gcf‘color‘‘w‘）;
set（gca‘units‘‘pixels‘‘Visible‘‘off‘）;
q=get（gca‘position‘）;
q（1）=0;%設置左邊距離值為零
q（2）=0;%設置右邊距離值為零
set（gca‘position‘q）;
for?i=1:length（stats）
????hold?on;
????rectangle（‘position‘stats（i）.BoundingBox‘edgecolor‘‘r‘‘linewidth‘2）;
????temp?=?stats（i）.Centroid;
????drawnow;
end??????????????????????????????????????%標記區域，并畫紅色矩形框
frame=getframe（gcf[00nm]）;
im=frame2im（frame）;??????????????????????%將圖片轉變成jpg格式
dif4=im2uint8（dif）;?%把dif4轉換為無符號8位整形
rgb=cat（3difdifdif）;?%將dif變為3維矩陣
rgb=uint8（rgb）;
rgb=rgb*255;
h=imabsdiff（imrgb）;%將im變成jpg格式數組，使得可以與p1p2相加
figure
subplot（121）imshow（h+p1）
subplot（122）imshow（h+p2）%顯示最終圖片

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????32114??2009-11-21?00:03??找茬?程序\1.jpg

?????文件??????26840??2011-09-04?00:04??找茬?程序\2.jpg

?????文件??????23510??2011-08-31?12:33??找茬?程序\3.jpg

?????文件??????20304??2011-08-31?12:34??找茬?程序\4.JPG

?????文件???????1170??2011-09-06?20:31??找茬?程序\main.asv

?????文件???????1170??2011-09-06?20:32??找茬?程序\main.m

????..A.SH.?????58880??2011-09-06?20:31??找茬?程序\Thumbs.db

?????目錄??????????0??2011-09-06?20:31??找茬?程序

-----------?---------??----------?-----??----

???????????????163988????????????????????8