A=imread（‘tire.tif‘）;?%打開原圖
[hw]=size（A）;???%取圖像的長寬，用在imtransform中
figure;
imshow（A）;title（‘前一幀‘）???????%輸出前一幀
tform=maketform（‘affine‘[1?0?5;?0?1?5;?0?0?1]‘）;
B=imtransform（Atform‘XData‘[1?w]‘YData‘[1?h]）;??????%得到處理素材
figure;
imshow（B）;title（‘當前幀‘）??????%輸出當前幀


m=0;?????????????????%建立變量m
n=0;?????????????????%建立變量n
[mn]=?size（A）;??????%把A中圖像行數賦給m，列數賦給n
Hproj=zeros（m1）;????%建立m*1的一維數組Hproj
Hproj1=zeros（m1）;???%建立m*1的一維數組Hproj1
Vproj=zeros（1n）;????%建立1*n的一維數組Vproj???
Vproj1=zeros（1n）;???%建立1*n的一維數組Vproj1
for?k=1:m?????
????Hproj（k）?=?sum（A（k:））;??
????Hproj1（k）?=?sum（B（k:））;
?end;????????????????%每行求和，賦給對應的Hproj數組中的元素
figure;plot（Hproj）;????xlabel（‘前一幀行投影‘）;%得到前一幀行投影???????
figure;plot（Hproj1）;?xlabel（‘當前幀行投影‘）;?%得到當前幀行投影
for?v=1:n
????Vproj（v）?=?sum（A（:v））;??
????Vproj1（v）?=?sum（B（:v））;?
?end;???????????????%每行求和，賦給對應的Hproj數組中的元素
figure;plot（Vproj）;??xlabel（‘前一幀列投影‘）;???????%輸出一維曲線
figure;plot（Vproj1）;?xlabel（‘當前幀列投影‘）;
cc=zeros（160）;?
for?i=1:60
????for?j=1:120
??????cc（i）=cc（i）+（Hproj（i+j）-Hproj1（30+j））^2;?????%行相關計算
????end?;
end?;
figure;plot（cc‘g‘）;???xlabel（‘垂直方向‘）?;
[simin]=min（cc）;?
dy=31-imin;?????????????????????%得到垂直偏移量
cr=zeros（160）;?
for?d=1:60
????for?b=1:120
??????cr（d）=cr（d）+（Vproj（d+b）-Vproj1（30+b））^2;?????%列相關計算
????end?;
end?;
figure;plot（cr‘r‘）;???xlabel（‘水平方向‘）?;??
[gdmin]=min（cr）;
dx=31-dmin;???????????????????????%得到水平偏移量
tform=maketform（‘affine‘[1?0?-dx;?0?1?-dy;?0?0?1]‘）;???%運動矢量補償
C=imtransform（Btform‘XData‘[1?w]‘YData‘[1?h]）;
figure;???????????????????????????????
imshow（C）;title（‘穩像后的當前幀‘）;??????%處理結果輸出
K=imsubtract（AB）;
figure;
imshow（K）;?title（‘穩像前兩幀差‘）???????%輸出穩像前兩幀差
H=imsubtract（AC）;
figure;
imshow（H）;title（‘穩像后兩幀差‘）;????????%輸出穩像后兩幀差

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2026??2012-06-05?07:49??wancheng.m