I=imread（‘circuit.tif‘）;
rotI=imrotate（I33‘crop‘）;
BW=edge（rotI‘canny‘）;
[HTR]=hough（BW）;
imshow（H[]‘XData‘R‘InitialMagnification‘‘fit‘）;
xlabel（‘\theta‘）ylabel（‘\rho‘）;
axis