Problème Template matching matlab

Problème Template matching matlab - Divers - Programmation

Marsh Posté le 13-07-2009 à 18:41:14    

Bonjour,
 
Je suis actuellement en stage au Canada et mon maitre de stage souhaiterait que je réalise un programme matlab permettant de tracker un objet sur une video et ensuite déterminer sa vitesse et tracer sa trajectoire (dendrite en photo ci-jointe). Dans mon école d'ingénieur, la formation a matlab est loin d'être une priorité et donc j'ai un peu de mal à me lancer dessus.
 
En regardant sur un forum, j'ai trouvé ce programme qui avait été posté et qui correspondrait bien voire très bien à ce que je veux.
 

Code :
  1. % this code is able to track an object within a video, to plot the
  2. % trajectory and to calculate the velocity of the object
  4. % initialisations
  6. template_rows = 1;
  7. template_columns = 1;
  8. template_xmin = 1;
  9. template_ymin = 1;
  13. % importation of the video
  15. video = aviread('number 2.avi');
  18. % video_info is used to store the information about the video
  19. % (number of frames, format of frames, number of frames per second...)
  21. video_info = aviinfo('number 2.avi');
  22. NumFrames = video_info.NumFrames;
  23. FPS = video_info.FramesPerSecond;
  26. % the crop function is used here to extract the ultrasound image from all
  27. % frames
  28. % for example sometimes the frames have a Title and some comments that the
  29. % program doesn't need for the processing
  30. % this enables to reduce the running time of the program
  31. % the user selects a rectangle containing the ultrasound image
  33. [crop,rect_image2] = imcrop(video(1).cdata);
  36. % the loop extract the ultrasound image part of each frames
  38. for j=1:NumFrames
  39.     video(j).cdata = imcrop(video(j).cdata,rect_image2);
  40.     video(j).cdata = rgb2gray(video(j).cdata);
  41. end
  44. image1 = video(1).cdata;
  45. [frame_rows,frame_columns] = size(image1);
  48. % the crop function is used here to select the object that he is interested
  49. % in finding the trajectory and velocity of
  51. [template,rect_image1] = imcrop(image1);
  53. % the code below is saving informations about the rectangle selected by the
  54. % user (size, position ...)
  56. size_template = size(template);
  57. template_rows = size_template(1);
  58. template_columns = size_template(2);
  59. template_width = template_columns - 1;
  60. template_height = template_rows - 1;
  61. if rect_image1(1)<1
  62.     rect_image1(1)=1;
  63. end
  64. if rect_image1(2)<1
  65.     rect_image1(2)=1;
  66. end
  69. % the coordinates of the object are saved
  71. x_coordinates(1,1) = round(rect_image1(1))+((template_columns-1)/2);
  72. y_coordinates(1,1) = round(rect_image1(2))+((template_rows-1)/2);
  75. % the loop below computes the cross-correlation and finds the coordinates of
  76. % the object in all frames of the video
  78. for k=2:NumFrames
  80.     image2 = video(k).cdata;
  82.     corr = normxcorr2(template,image2);
  84.     % the program using the function max_corr
  85.     [corr_max_row,corr_max_column]=max_corr(corr);
  87.     % this part is computing the updating of the template
  89.     template_xmin = corr_max_column - template_width;
  90.     template_ymin = corr_max_row - template_height;
  91.     rect_template = [template_xmin, template_ymin, template_width, template_height];
  92.     template = imcrop(image2,rect_template);
  95.     % the coordinates of the object are saved
  97.     x_coordinates(1,k) = corr_max_column - ((template_columns-1)/2);
  98.     y_coordinates(1,k) = corr_max_row - ((template_rows-1)/2);
  100. end
  103. for i=1:NumFrames
  104.     y_coordinates(1,i)=frame_rows - y_coordinates(1,i);
  105. end
  108. % the program asks the user to enter the frame length scale
  110. Scale = input('Dimension of the frames (in mm):\n');
  113. % the pixel length and the time between two consecutive frames are
  114. % calculated
  116. pixel_length = Scale/frame_columns;
  117. deltaT = 1/FPS;
  119. % the loop below calculates the velocity of the object for each frames
  121. for k=2:NumFrames-1
  122.     offset_x = sqrt((x_coordinates(1,k+1)-x_coordinates(1,k-1))^2);
  123.     offset_y = sqrt((y_coordinates(1,k+1)-y_coordinates(1,k-1))^2);
  124.     velocity(1,k-1) = (sqrt(offset_x^2 + offset_y^2))/(2*deltaT);
  125. end
  126.        
  127. % this corresponds to the plotting of the trajectory 
  128.    
  129. subplot(2,1,1), plot(x_coordinates,y_coordinates,'o');
  130. axis([0 frame_columns 0 frame_rows]);
  131. xlabel('Pixels in the x axis');
  132. ylabel('Pixels in the y axis');
  133. title('trajectory of the object')
  134. subplot (2,1,2), plot(velocity);
  135. title ('velocity of the object');
  138. % the code below displays the velocity table
  140. velocity(:,:)


 
 
 Mais lorsque je tente de le lancer sur ma video, il y a un problème avec la fonction max_corr:
 
??? Undefined function or method 'max_corr' for input arguments of type 'double'.
 
Error in ==> vitesse at 85
[corr_max_row,corr_max_column]=max_corr(corr);
 
Cette fonction devait etre realiser par un programme indépendant de celui-ci. Quelqu'un de familier avec le template matching pourrait il m'expliquer ce que réalsie cette fonction?
 
merci beaucoup pour votre aide

Reply

Marsh Posté le 13-07-2009 à 18:41:14   

Reply

Marsh Posté le 13-07-2009 à 20:01:42    

max_corr effectue tres certainement un correletaion entre ta frame et le pattern a cherché et cherche le maximum de cette correlation dans l'image et renvoit les coordonnées.

Reply

Marsh Posté le 13-07-2009 à 23:18:56    

Mais cela correspond au numéro de la ligne et de la colonne du coefficient maximum de la matrice corr?

Reply

Marsh Posté le 14-07-2009 à 10:30:05    

oui

Reply

Sujets relatifs:

Leave a Replay

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