如何在阈值处理和图像标记后为每个路灯提取houghpeaks？

Question

我正在做matlab图像处理的最后一年项目。我的项目是在静态图像中识别街道中哪些灯“关闭”。

这些是我正在遵循的步骤序列：

我已经进行了阈值处理和图像标记，即在阈值图像中标记每个光源，依此类推。相应的代码是

I = imread('st.jpg');
imshow(I);
G = rgb2gray(I);
imshow(G);
BW = im2bw(G,0.7);
imshow(BW);
f = fspecial('average',3);
I1 = filter2(f,BW,'same');
imshow(I1);
[Lw, nw] = bwlabel(I1);
imtool(Lw,[]);

s = strel('disk',5);
O = imopen(Lw,s);
figure, imshow(O,[]);

现在我需要创建每个光源的数据库，并将测试图像（关闭一个/多个光源的图像）与参考图像进行比较，以确定哪个关闭并显示它。
因为每个光源的houghpeaks是不同的。我想用霍夫变换为每个光源提取houghpeaks 我已经搜索了这么多，但我得到了任何东西。如果有人帮助我，这将对我的项目有很大的帮助 谢谢。