I am newbie to image processing,I work on retinal image analysis. I use gabor filter to perform some operation on my image. I used gabor filter for 12 different orientations of an retinal image and i got 12 different output. now i want to superimpose these images on my original image(first visualized image with original image,second visualized image with original image...till the 12th imgage) and to display all the 12 images.How do i do it?
here is my code
Code: Select all
I = getimage();
I=I(:,:,2);
lambda = 8;
theta = 0;
psi = [0 pi/2];
gamma = 0.5;
bw = 1;
N = 12;
img_in = im2double(I);
%img_in(:,:,2:3) = []; % discard redundant channels, it's gray anyway
img_out = zeros(size(img_in,1), size(img_in,2), N);
imshow(img_in);
title('input image');
for n=1:N
gb = gabor_fn(bw,gamma,psi(1),lambda,theta)...
+ 1i * gabor_fn(bw,gamma,psi(2),lambda,theta);
img_out(:,:,n) = imfilter(img_in, gb, 'symmetric');
imshow(img_out(:,:,n)); //gives me 12 different visualization output
theta = 15 * n; % next orientation
end
https://www.dropbox.com/s/on1zfnugs3gtm8q/green.jpg?m
and this is the output which i got using gabor filter
https://www.dropbox.com/s/nnevl4swc7jzzm3/gab.jpg
and this is the type of image i want to get..(this is not the exact image)
https://www.dropbox.com/s/riiqd83tpvvdu8q/03_test.jpg?m