Bio-medical Ontologies Maintenance and Change Management

2
D (Q, T)
= ∑∑ fϕ
fϕ
f ϕ
Q
d ( Q,
T ) , where d = ( f − f
Multimedia Medical Databases 113
T
)
2
(4.4)
The algorithm in pseudo-cod for detecting color texture with Gabor method is:
**Procedure Gabor_filter
int scale = 3; int orientation = 4; double Ul = 0.1;
double Uh = 0.4; int flag = 0; int side = 60;
//the main function
*function GaborFilteredImg(feature[][], Matrix
img,Matrix img_i, int side, double Ul, double Uh, int
scale, int orientation, int flag)
begin
int border = side;int height = img.getNRows();
int width = img.getNCols();
int xs = height+2.0*border;
int ys = width+2.0*border;
Matrix IMG (xs,ys); Matrix IMG_imag (xs,ys);
Matrix F_real (xs,ys);
Matrix F_imag(xs,ys);
Matrix G_real (xs,ys); Matrix G_imag (xs,ys);
Matrix Gr (2*side+1, 2*side+1);
Matrix Gi (2*side+1, 2*side+1);
Matrix Tmp_1 (xs,ys); Matrix Tmp_2 (xs,ys);
Matrix F_1 (xs,ys); Matrix F_2 (xs,ys);
//The Fourier transform of the image matrix
*Function Matrix FFT2D(F_real, F_imag, IMG,
IMG_imag);
for s = 0; scale; do
for n = 0; orientation do
*Function Gabor(Gr, Gi, s+1, n+1, Ul, Uh, scale,
orientation, flag);
//The Fourier transform of the Gabor filter
*Function Matrix FFT2D(G_real, G_imag, F_1, F_2);
for i=0;xs do
for j=0;ys do
//The product of the two Fourier transforms
Tmp_1.set(i,j, G_real[i,j]* F_real[i,j]);
Tmp_2.set(i,j, G_imag[i,j]* F_imag[i,j]);
IMG.set(i,j , Tmp_1[i,j] - Tmp_2.[i,j]);
Tmp_1.set(i,j, G_real[i,j] * F_imag[i,j]);
Tmp_2.set(i,j,G_imag[i,j] * F_real[i,j]);
IMG_imag.set(i,j,Tmp_1.[i,j] + Tmp_2[i,j]);
End;
End;