       RGB = imread('radarData22222-S3.png');
       I  = rgb2gray(RGB); % convert to intensity
       BW = edge(I,'canny'); % extract edges
       [H,T,R] = hough(BW,'RhoResolution',0.5,'Theta',-90:0.5:89.5);
 
       figure;
       % display the original image
       subplot(2,1,1);
       imshow(RGB);
       title('RadarData22222Scan3.png');
 
       % display the hough matrix
       subplot(2,1,2);
       imshow(imadjust(rescale(H)),'XData',T,'YData',R,...
              'InitialMagnification','fit');
       title('Hough transform of RadarData22222Scan3.png');
       xlabel('\theta'), ylabel('\rho');
       axis on, axis normal, hold on;
       colormap(gca,hot);
