基于多尺度融合的布匹瑕疵检测技术（代码完整，数据齐全）

%%%%%%%VERSION 2 %%ANOTHER DESCRIBTION OF GABOR FILTER %The Gabor filter is basically a Gaussian (with variances sx and sy along x and y-axes respectively) %modulated by a complex sinusoid (with centre frequencies U and V along x and y-axes respectively) %described by the following equation %% % -1 x' ^ y' ^ %%% G(x,y,theta,f) = exp ([----{(----) 2+(----) 2}])*cos(2*pi*f*x'); % 2 sx' sy' %%% x' = x*cos(theta)+y*sin(theta); %%% y' = y*cos(theta)-x*sin(theta); %% Describtion : %% I : Input image %% Sx & Sy : Variances along x and y-axes respectively %% f : The frequency of the sinusoidal function正弦函数的频率 %% theta : The orientation of Gabor filter滤波器的方向 %% G : The output filter as described above输出滤波器 %% gabout : The output filtered image过滤的图像 %% Author : Ahmad poursaberi e-mail : a.poursaberi@ece.ut.ac.ir %% Faulty of Engineering, Electrical&Computer Department,Tehran %% University,Iran,June 2004 function gabout=gaborfilter2(I,Sx,Sy,f,theta); if isa(I,'double')~=1 I = double(I); end for x = -fix(Sx):fix(Sx) %fix:朝零方向取整 for y = -fix(Sy):fix(Sy) xPrime = x * cos(theta) + y * sin(theta);%将 Gabor(x,y)中的x,y坐标以theta角度进行方向坐标变换: yPrime = y * cos(theta) - x * sin(theta); G(fix(Sx)+x+1,fix(Sy)+y+1) = exp(-.5*((xPrime/Sx)^2+(yPrime/Sy)^2)) * cos(2*pi*f*xPrime)/(2*pi); end end Imgabout = conv2(I,double(imag(G)),'same'); Regabout = conv2(I,double(real(G)),'same'); gabout = sqrt(Imgabout.*Imgabout + Regabout.*Regabout); %G %size(G) %figure;imshow(G,[]); %figure;g=log(abs(fftshift(fft2(double(G)))));imshow(g,[]); %figure;imshow(gabout,[]); %figure;g=log(abs(fftshift(fft2(double(gabout)))));imshow(g,[]);