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
|
function H=doog1(sig,r,th,N);
% H=doog1(sig,r,th,N);
% by Serge Belongie
no_pts=N; % no. of points in x,y grid
[x,y]=meshgrid(-(N/2)+1/2:(N/2)-1/2,-(N/2)+1/2:(N/2)-1/2);
phi=pi*th/180;
sigy=sig;
sigx=r*sig;
R=[cos(phi) -sin(phi); sin(phi) cos(phi)];
C=R*diag([sigx,sigy])*R';
X=[x(:) y(:)];
Gb=gaussian(X,[0 0]',C);
Gb=reshape(Gb,N,N);
m=R*[0 sig]';
a=1;
b=-1;
% make odd-symmetric filter
Ga=gaussian(X,m/2,C);
Ga=reshape(Ga,N,N);
Gb=rot90(Ga,2);
H=a*Ga+b*Gb;
|
