-
Notifications
You must be signed in to change notification settings - Fork 3
/
log_srch_allFrames.m
105 lines (92 loc) · 2.24 KB
/
log_srch_allFrames.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
clc
clear all
close all
tstart=cputime;
f_ref(1:300,1:300)=0;
Im=imread('E:\foreman_10frames\f001.pgm');
f_ref(9:288,9:288)= Im(9:288,9:288);
srcFiles = dir('E:\foreman_10frames\*.pgm');
f_p=zeros(300,3000);
X=zeros(35,350);
Y=zeros(35,350);
SAD=zeros(1,9);
for frameNo=1:9
filename = strcat('E:\foreman_10frames\',srcFiles(frameNo+1).name);
f_2(1:300,1:300)=0;
Im2= imread(filename);
f_2(9:288,9:288)= Im2(9:288,9:288);
f_pre(1:300,1:300)=0;
s=1;
X_motion= zeros(22,22);
Y_motion=zeros(22,22);
for i=9:8:288
t=1;
for j=9:8:288
img_abs=[0 0 0 0 0];
img_16=f_ref(i-4:i+7+4,j-4:j+7+4);
img_8=f_2(i:i+7,j:j+7);
flag = 1;
I=1;
Rc=5;
Cc=5;
step_size=4;
while flag
r=[Rc,Rc-step_size,Rc,Rc,Rc+step_size];
c=[Cc,Cc,Cc-step_size,Cc+step_size,Cc];
for g=1:5
if r(g)<=0 || r(g)>=10
img_abs(g)=255*64*255;
elseif c(g)<=0 || c(g)>=10
img_abs(g)=255*64*255;
elseif g==I && g-1>0
img_abs(g)=255*64*255;
else
img_abs(g)=sum(sum((img_16(r(g):r(g)+7,c(g):c(g)+7)- img_8).^2));
end
end
[M,I] = min(img_abs);
switch (I)
case 1
step_size=step_size/2;
case 2
Rc=Rc-step_size;
case 3
Cc=Cc-step_size;
case 4
Cc=Cc+step_size;
case 5
Rc=Rc+step_size;
end
if step_size<1
flag=0;
break;
else
continue;
end
end
f_pre(i:i+7,j:j+7)=img_16(Rc:Rc+7,Cc:Cc+7);
X_motion(s,t)= Rc-5;
Y_motion(s,t)= Cc-5;
t=t+1;
end
s=s+1;
end
f_p(1:300, 1+(300*frameNo):300*(frameNo+1))=f_pre;
X(1:35, 1+(35*(frameNo-1)):35*frameNo)=X_motion;
Y(1:35, 1+(35*(frameNo-1)):35*frameNo)=Y_motion;
residu=abs(f_2-f_pre);
SAD(frameNo)=(sum(sum((residu).^2)))/90000;
figure,imshow(uint8(residu));
title('reduced residue after Logarithmic Search Operation');
figure,imshowpair(f_2,f_ref,'diff');
title('actual residue or difference between frames');
f_ref=f_2;
end
telapsed=cputime-tstart;
Frame=[1 2 3 4 5 6 7 8 9];
figure,LogSearch=plot(Frame,SAD);
title('Sum of Absolute Difference [SAD] Vs Frames Plot');
ylabel('SAD found in Logarithmic Search');
xlabel('Frame number');
display('time elapsed in search');
display(telapsed);