find peaks of hough transform and extend the hough lines in matlab

main.m
% Load the image and perform edge detection
img = imread('image.png');
bw = edge(img);

% Compute the Hough transform
[H,theta,rho] = hough(bw);

% Find the peaks in the Hough space
peaks = houghpeaks(H, 10); % finds the top 10 peaks

% Plot the Hough space with the peaks overlaid
figure, imshow(imadjust(mat2gray(H)),[],'XData',theta,'YData',rho,...
              'InitialMagnification','fit');
xlabel('\theta'), ylabel('\rho');
axis on, axis normal, hold on;
plot(theta(peaks(:,2)),rho(peaks(:,1)),'s','color','white');
523 chars
17 lines

main.m
% Find the lines corresponding to the Hough peaks
lines = houghlines(bw, theta, rho, peaks);

% Plot the original image with the extended Hough lines overlaid
figure, imshow(img), hold on
max_len = 0;
for k = 1:length(lines)
   xy = [lines(k).point1; lines(k).point2];
   
   % Extend the lines to the boundary of the image
   plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','green');
   
   % Compute the length of the line segment
   len = norm(lines(k).point1 - lines(k).point2);
   
   % Keep track of the longest line found
   if ( len > max_len)
      max_len = len;
      xy_long = xy;
   end
end

600 chars
23 lines