1）读取图像‘image.png’并将其转换为灰度图像I； 2）使用edge函数对I进行边缘点检测，检测使用‘Canny’算子，输出边缘图像BW（BW是二值图像，边缘点为白色，非边缘点为黑色）； 3）显示边缘图像4）使用hough函数对BW进行Hough变换得到变换域H,theta以及rho； 5）使用imshow显示变换域H（注：显示之前先使用rescale函数将H的数值范围缩放到[0,1]之间，再使用imadjust函数进行灰度调整），显示时imshow的’XData’设置为theta，以设置图像x轴的范围，’YData’设置为rho，以设置图像y轴的范围，且参数‘InitialMagnification’设置为‘fit’以缩放整个图像以适合窗口。使用xlabel和ylabel函数设置x轴和y轴的标签，通过axis on设置坐标轴的显示，axis normal自动调节坐标轴的纵横比。 6）使用houghpeaks找出前5个极大值并将极大值的坐标（即theta和rho的索引）输出为数组P； 7）通过数组P的坐标信息找到对应的（theta，rho）对，并在5）中得到的H图像中将极大值通过plot函数标示出来，极大值使用红色空心圆标示。

这是一段MATLAB代码，可以实现你提出的需求：

% 读取图像并转换为灰度图像
I = rgb2gray(imread('image.png'));

% 对灰度图像进行边缘检测
BW = edge(I, 'Canny');

% 显示边缘图像
imshow(BW);

% 对边缘图像进行Hough变换
[H,theta,rho] = hough(BW);

% 显示变换域H
H = rescale(H, 0, 1); % 缩放数值范围到[0,1]
H = imadjust(H); % 灰度调整
imshow(theta,rho,H,[],'InitialMagnification','fit');
xlabel('\theta'), ylabel('\rho');
axis on, axis normal;

% 找到前5个极大值
P = houghpeaks(H, 5, 'Threshold', ceil(0.3*max(H(:))));

% 在变换域H中标示极大值
hold on;
plot(theta(P(:,2)),rho(P(:,1)),'rs','LineWidth',2);
hold off;

这段代码首先读取图像并将其转换为灰度图像，然后使用edge函数对灰度图像进行边缘检测，输出二值图像BW。接着，使用hough函数对BW进行Hough变换，得到变换域H，theta和rho。将H缩放到[0,1]范围内，再进行灰度调整，然后使用imshow函数将H显示出来，设置XData和YData参数，以及xlabel和ylabel函数，控制坐标轴的显示。接着，使用houghpeaks函数找到前5个极大值，并将它们的坐标输出为数组P。最后，在变换域H中标示这些极大值，使用plot函数，标记为红色空心圆。