我正在使用3D图像,我想实现Gabor功能。为此,我需要生成具有不同比例和方向(也许以3D,两个角度)的Gabor滤波器组,然后将其与我的图像融合。我的问题是,是否有可能实施3D gabor滤镜,然后将其与3D图像融合。还是我要在2D切片中逐片进行。如果是这样,将不胜感激在matlab中实现它的任何建议。
先感谢您
使其成为2D。Gabor滤镜会自动将图像更改为灰度,因此无论如何3D变得毫无用处。
这是一个用于matlab的gabor过滤器,我使用了它,它运行良好。
function gaborArray = gaborFilterBank(u,v,m,n)
% GABORFILTERBANK generates a custum Gabor filter bank.
% It creates a u by v array, whose elements are m by n matries;
% each matrix being a 2-D Gabor filter.
%
%
% Inputs:
% u : No. of scales (usually set to 5)
% v : No. of orientations (usually set to 8)
% m : No. of rows in a 2-D Gabor filter (an odd integer number usually set to 39)
% n : No. of columns in a 2-D Gabor filter (an odd integer number usually set to 39)
%
% Output:
% gaborArray: A u by v array, element of which are m by n
% matries; each matrix being a 2-D Gabor filter
%
%
% Sample use:
%
% gaborArray = gaborFilterBank(5,8,39,39);
%
%
% Details can be found in:
%
% M. Haghighat, S. Zonouz, M. Abdel-Mottaleb, "Identification Using
% Encrypted Biometrics," Computer Analysis of Images and Patterns,
% Springer Berlin Heidelberg, pp. 440-448, 2013.
%
%
% (C) Mohammad Haghighat, University of Miami
% [email protected]
% I WILL APPRECIATE IF YOU CITE OUR PAPER IN YOUR WORK.
if (nargin ~= 4) % Check correct number of arguments
error('There should be four inputs.')
end
%% Create Gabor filters
% Create u*v gabor filters each being an m*n matrix
gaborArray = cell(u,v);
fmax = 0.25;
gama = sqrt(2);
eta = sqrt(2);
for i = 1:u
fu = fmax/((sqrt(2))^(i-1));
alpha = fu/gama;
beta = fu/eta;
for j = 1:v
tetav = ((j-1)/v)*pi;
gFilter = zeros(m,n);
for x = 1:m
for y = 1:n
xprime = (x-((m+1)/2))*cos(tetav)+(y-((n+1)/2))*sin(tetav);
yprime = -(x-((m+1)/2))*sin(tetav)+(y-((n+1)/2))*cos(tetav);
gFilter(x,y) = (fu^2/(pi*gama*eta))*exp(-((alpha^2)*(xprime^2)+(beta^2)*(yprime^2)))*exp(1i*2*pi*fu*xprime);
end
end
gaborArray{i,j} = gFilter;
end
end
%% Show Gabor filters
% Show magnitudes of Gabor filters:
figure('NumberTitle','Off','Name','Magnitudes of Gabor filters');
for i = 1:u
for j = 1:v
subplot(u,v,(i-1)*v+j);
imshow(abs(gaborArray{i,j}),[]);
end
end
% Show real parts of Gabor filters:
figure('NumberTitle','Off','Name','Real parts of Gabor filters');
for i = 1:u
for j = 1:v
subplot(u,v,(i-1)*v+j);
imshow(real(gaborArray{i,j}),[]);
end
end
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句