Рассчитать углы от извлеченных вен листьев - PullRequest
Новая
       5

Рассчитать углы от извлеченных вен листьев

0 голосов
/ 01 ноября 2019

Я пытаюсь вычислить углы, которые существуют между жилками листа в matlab. Я подумал о том, чтобы обнаружить линии (вены) с помощью Hough, а затем вычислить углы. Но результат, который я получаю, когда пытаюсь обнаружить строки, немного странен: enter image description here

Вот код: 

[H,theta,rho] = hough(vein);
P = houghpeaks(H,4);
x = theta(P(:,2));
y = rho(P(:,1));
%lines = houghlines(BW,theta,rho,P,'FillGap',2,'MinLength',4);
lines = houghlines(BW,theta,rho,P);


figure, imshow(vein), hold on
max_len = 0;
for k = 1:length(lines)
   xy = [lines(k).point1; lines(k).point2];
   plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','green');

   % Plot beginnings and ends of lines
   plot(xy(1,1),xy(1,2),'x','LineWidth',2,'Color','yellow');
   plot(xy(2,1),xy(2,2),'x','LineWidth',2,'Color','red');

   % Determine the endpoints of the longest line segment
   len = norm(lines(k).point1 - lines(k).point2);
   if ( len > max_len)
      max_len = len;
      xy_long = xy;
   end
end
% highlight the longest line segment
plot(xy_long(:,1),xy_long(:,2),'LineWidth',2,'Color','red');

Это мой оригинализображение вены: enter image description here

Есть идеи?

1 Ответ

1 голос
/ 02 ноября 2019

Полагаю, вы просто перепутали названия матриц.

В вашем коде есть vein матрица и BW матрица.
Вы используете:

  • вена здесь: [H,theta,rho] = hough(vein);
  • BW здесь: lines = houghlines(BW,theta,rho,P);

Я добавил строку BW = vein;, и результат намного меньше странно : 

vein = imbinarize(rgb2gray(imread('https://i.stack.imgur.com/cuwb0.jpg'))); %Read image, and convert to binary
BW = vein; %Copy vein to BW

[H,theta,rho] = hough(vein);
P = houghpeaks(H,4);
x = theta(P(:,2));
y = rho(P(:,1));
%lines = houghlines(BW,theta,rho,P,'FillGap',2,'MinLength',4);
lines = houghlines(BW,theta,rho,P);


figure, imshow(vein), hold on
max_len = 0;
for k = 1:length(lines)
   xy = [lines(k).point1; lines(k).point2];
   plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','green');

   % Plot beginnings and ends of lines
   plot(xy(1,1),xy(1,2),'x','LineWidth',2,'Color','yellow');
   plot(xy(2,1),xy(2,2),'x','LineWidth',2,'Color','red');

   % Determine the endpoints of the longest line segment
   len = norm(lines(k).point1 - lines(k).point2);
   if ( len > max_len)
      max_len = len;
      xy_long = xy;
   end
end
% highlight the longest line segment
plot(xy_long(:,1),xy_long(:,2),'LineWidth',2,'Color','red');

Результат:

enter image description here

...