[matlab] PB laplacien Canny

PB laplacien Canny [matlab] - Divers - Programmation

Marsh Posté le 02-05-2012 à 19:37:15    

Bonjour à tous,
 
Je suis entrain de tester l'algorithme de Canny pour détecter des contours dans une image. Mon souci c'est que je ne comprends pas pourquoi lorsque j'affiche l'image du laplacien mes contours n'apparaissent pas en noir et les zones uniformes en blanc. Normalement les contours sont les zéros du laplacien, je ne vois pas mon erreur pourtant mon gradient semble juste. Je n'arrive pas à obtenir les contours à l'aide du laplacien. Ai-je fais une erreur de programmation.
Pouvez vous m'aidez svp ?
 

Code :
  1. close all
  2. clear
  5. %%%recuperation de l'image
  6. I=imread('image/lena.jpg');   
  7. I=rgb2gray(I);
  8. I=double(I)/255.0;
  10. figure,imshow(I,'notruesize'),title('image originale');
  11. figure,imshow(edge(I,'canny')),title('image contour matlab');
  14. %%%%%%%%%%%%%%%%%%% filtre de Canny%%%%%%%%%%%%%%%%%%%%%
  16. % creation du filtre
  18. n=3; %taille du filtre pour le calcul du gradient n=2*tailleFenetre +1
  19. n1=5; % taille du filtre pour le laplacien du gradient
  21. tailleFenetre=(n-1)/2;
  22. tailleFenetre1=(n1-1)/2;
  23. sigma2=(n/(2*3))^2; % Nota bene : rayon=3*sigma (gaussienne)
  25. %% calul du gradient
  26. h_canny=zeros(n,n); % filtre horizontal pour le grandient en x
  27. for i1=-tailleFenetre:tailleFenetre
  28. for j1=-tailleFenetre:tailleFenetre
  29.    h_canny(tailleFenetre+i1 +1,tailleFenetre+j1+1)=(-i1/sigma2)*exp(-(i1*i1 +j1*j1)/2*sigma2);
  30. end
  31. end
  33. v_canny=zeros(n,n);% filtre vertical pour le gradient en y
  34. for i1=-tailleFenetre:tailleFenetre
  35. for j1=-tailleFenetre:tailleFenetre
  36.    v_canny(tailleFenetre+i1 +1,tailleFenetre+j1+1)=(-j1/sigma2)*exp(-(i1*i1 +j1*j1)/2*sigma2);
  37. end
  38. end
  40. %% laplacien (derivee seconde)
  41. sigma2=(n1/(2*3))^2; % Nota bene : rayon=3*sigma
  42. sigma4=sigma2*sigma2;
  44. h2_canny=zeros(n1,n1); %filtre horizontal pour le laplacien en x
  45. for i1=-tailleFenetre1:tailleFenetre1
  46. for j1=-tailleFenetre1:tailleFenetre1
  47.    h2_canny(tailleFenetre1+i1 +1,tailleFenetre1+j1+1)=((i1*i1 - sigma2)/sigma4)*exp(-(i1*i1+j1*j1)/2*sigma2);
  48. end
  49. end
  51. v2_canny=zeros(n1,n1); %filtre vertical pour le laplacien en y
  52. for i1=-tailleFenetre1:tailleFenetre1
  53. for j1=-tailleFenetre1:tailleFenetre1
  54.    v2_canny(tailleFenetre1+i1 +1,tailleFenetre1+j1+1)=((j1*j1 - sigma2)/sigma4)*exp(-(i1*i1+j1*j1)/2*sigma2);
  55. end
  56. end
  58. v2_canny(tailleFenetre1+1,tailleFenetre1+1)=0 % modification de l'element central pour que la somme des termes de la matrice soit nulle
  59. v2_canny(tailleFenetre1+1,tailleFenetre1+1)=-sum(sum((v2_canny)));
  60. h2_canny(tailleFenetre1+1,tailleFenetre1+1)=0
  61. h2_canny(tailleFenetre1+1,tailleFenetre1+1)=-sum(sum(h2_canny));
  63. %%%%%%%%%%%%%%%% calcul de la norme du gradient
  64. gradH=filter2(h_canny,I);
  65. gradV=filter2(v_canny,I);
  66. grad=sqrt(gradH.*gradH + gradV.*gradV);
  67. %grad=abs(gradH)+abs(gradV);
  68. figure,imshow(grad,[min(min(grad)) max(max(grad))],'notruesize'),title('norme gradien');
  71. %%%%%%%%%%%%%%%%%%%%%%%%%%%% calcul du laplacien %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  73. IlaplacienH=filter2(h2_canny,I);
  74. IlaplacienV=filter2(v2_canny,I);
  75. Ilaplacien=IlaplacienH+IlaplacienV ;
  76. figure,imshow(abs(Ilaplacien),[min(min(abs(Ilaplacien))) max(max(abs(Ilaplacien)))],'notruesize');
  77. contourLaplacien=(abs(Ilaplacien)<0.2);
  78. figure,imshow(contourLaplacien);

Reply

Marsh Posté le 02-05-2012 à 19:37:15   

Reply

Sujets relatifs:

Leave a Replay

Make sure you enter the(*)required information where indicate.HTML code is not allowed