How to get a defocus image from single image to which already applied wavelet transform and canny edge detector ?

1 次查看（过去 30 天）
Haseeb Hassan
Haseeb Hassan 2018-4-18
编辑： Haseeb Hassan 2018-4-18
Hello everyone. I have a color image, trying to get first the wavelets transform and then apply canny edge detctor
separately on the color image.After getting wavelet transform of an image and edges through canny detector, then
want to get the focus image from the two previous applied operations.The algorthm should be look like this.
1-colar Image--->Wavelet Transform ,
2-color Image-->Canny Edge Detector ,
3- Wavelet Transform Image & canny edges Image----> Defocus Image.
Please see in the image what actually i want to do.
clear all ; close all
x = imread('football.jpg');
figure; imshow(x);
%%Applying db2 for getting Avg Aprroximation,Horizontal,Vertical and Diagonal Details
[xar,xhr,xvr,xdr] = dwt2(x(:,:,1),'db2');
[xag,xhg,xvg,xdg] = dwt2(x(:,:,2),'db2');
[xab,xhb,xvb,xdb] = dwt2(x(:,:,3),'db2');
xa(:,:,1) = xar ; xa(:,:,2) = xag ; xa(:,:,3) = xab ;
xh(:,:,1) = xhr ; xh(:,:,2) = xhg ; xh(:,:,3) = xhb ;
xv(:,:,1) = xvr ; xv(:,:,2) = xvg ; xv(:,:,3) = xvb ;
xd(:,:,1) = xdr ; xd(:,:,2) = xdg ; xd(:,:,3) = xdb ;
%%figure, imshow(xa/255) ;
%%figure, imshow(xh) ;
%%figure, imshow(xv) ;
%%figure, imshow(xd) ;
%%To get the window in standarized Form
X1 = [xa*0.003 log10(xv)*0.3 ; log10(xh)*0.3 log10(xd)*0.3] ;
figure ; imshow(X1)
%%Agian applying db2 to the subsampled Avg Approximationation Image, Vertical,horizontal details and diagonal
details.
[xaar,xhhr,xvvr,xddr] = dwt2(xa(:,:,1),'db2');
[xaag,xhhg,xvvg,xddg] = dwt2(xa(:,:,2),'db2');
[xaab,xhhb,xvvb,xddb] = dwt2(xa(:,:,3),'db2');
xaa(:,:,1) = xaar ; xaa(:,:,2) = xaag ; xaa(:,:,3) = xaab ;
xhh(:,:,1) = xhhr ; xhh(:,:,2) = xhhg ; xhh(:,:,3) = xhhb ;
xvv(:,:,1) = xvvr ; xaa(:,:,2) = xvvg ; xvv(:,:,3) = xvvb ;
xdd(:,:,1) = xddr ; xdd(:,:,2) = xddg ; xdd(:,:,3) = xddb ;
%%figure,imshow(xaa/255);
%%figure,imshow(xhh);
%%figure,imshow(xvv);
%%figure,imshow(xdd);
X11 = [ xaa*0.001 log10(xvv)*0.3 ; log10(xhh)*0.3 log10(xdd)*0.3 ] ;
%%figure; imshow(X11)
[r,c,s] = size(xv);
figure ; imshow([X11(1:r,1:c,:) xv*0.05 ; xh*0.05 xd*0.05 ])
%%Jacobian Edges Applied(https://www.mathworks.com/matlabcentral/fileexchange/28114-fast-edges-of-a-color-image--
actual-color--not-converting-to-grayscale-)
%get color edges and normalize magnitude
C = coloredges(x);
C = C / max(C(:));
%get grayscale edges and normalize magnitude
G_image = single(rgb2gray(x)) / 255;
G = sqrt(imfilter(G_image, fspecial('sobel')').^2 + imfilter(G_image, fspecial('sobel')).^2);
G = G / max(G(:));
%show comparison
figure, imshow(uint8(255 * cat(3, C, G, G)))

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Object Analysis 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by