0

二进制图像

我已经骨架化了二进制图像和连接信息。我想以连接点为中心画圆，并想找到圆和二值图像的交点。我写了以下代码：

 BW = imread('circles.png');
 imshow(BW);
 BW2 = bwmorph(BW,'remove');
 figure, imshow(BW2)
 BW3 = bwmorph(BW,'skel',Inf);
 figure, imshow(BW3)
 BW3t = bwmorph(BW3,'thin');
 figure, imshow(BW3t)


 [rj, cj, re, ce] = findendsjunctions(BW3t, 1);
 hold on 
 plot(cj(1),rj(1),'ob')
 hold on
 circle([cj(1),rj(1)],4,50,':r');

findendsjunctions.m 和相关文件 show.m 可以从这里下载：http ://www.csse.uwa.edu.au/~pk/research/matlabfns/LineSegments/findendsjunctions.m和这里http://www.csse。 uwa.edu.au/~pk/research/matlabfns/Misc/show.m。circle.m 可以从这里下载：http: //www.mathworks.co.uk/matlabcentral/fileexchange/2876-draw-a-circle/content/circle.m

我想知道圆圈是否与它周围的 2、3 或 4 条血管相交（在图像中标记为星形）。即使圆形横穿单个血管多次，但输出应该是每个血管的一个交点。

请建议我如何找到圆形和二元血管的交集。

谢谢

4

1 回答 1

0

我找到了圆和二进制图像的交点以及 3 个点的坐标（在我的问题提供的图像中标记为星形）。我已经更改了函数 circle.m （在我上面的问题中提到）以输出圆的圆周的所有 X 和 Y 坐标，然后我编写了以下 matlab 代码：

   [H, X, Y]=circle([cj(1),rj(1)],4,50,':r');
   c = improfile(BW3t,X,Y)
   x=1:length(c)
   figure
   plot(x, c,'r')

峰值的数量表示圆形切割二值图像的位置

   [maxtab, mintab]=peakdet(c, 1)
   [pks,locs] = findpeaks(c)
   pt1=[X(locs(1)) Y(locs(1))]
   pt2=[X(locs(2)) Y(locs(2))]
   pt3=[X(locs(3)) Y(locs(3))]

   hold on
   plot(pt1(1),pt1(2),'om','LineWidth',2)
   hold on
   plot(pt2(1),pt2(2),'og','LineWidth',2)
   hold on
   plot(pt3(1),pt3(2),'ob','LineWidth',2)

pt1, pt2 pt3 是圆切割二值图像的三个点

于 2013-06-02T22:49:26.343 回答