# I want to find end to end distance of this skeleton image

3 ビュー (過去 30 日間)
Adarsh Tripathi 2022 年 6 月 25 日
コメント済み: Adarsh Tripathi 2022 年 6 月 26 日
k = rgb2gray(k);
%imshow(insect);
%applyig threshold to convert the image into binary so that backgroung can
%be separated
BI = k < 160;
%imshow(BI)
BI = bwareaopen(BI,40);
%imshow(BI)
%finding the skeleton of image
k = bwskel(BI);
%imshow(k);
[a b]=size(k);
output=zeros(a,b);
for i=2:a-1
for j=2:b-1
ws=[k(i-1,j-1),k(i-1,j),k(i-1,j+1),k(i,j-1),k(i,j)...
,k(i,j+1),k(i+1,j-1),k(i+1,j),k(i+1,j+1)];
kg=sum(ws);
if(k(i,j)==1 && kg==2)
output(i,j)=1;
end
end
end
figure;
output=imdilate(output,strel('diamond',5));
imshow(output+k); i have this code with which i highlighted the end points but dont know how to go further
##### 2 件のコメント表示非表示 1 件の古いコメント
Adarsh Tripathi 2022 年 6 月 25 日

Thanks for your reply @walter Roberson Actually I attached the image before the code named insect in which end points are clearly specified i don't know how this below image is visible. And also just after sometime of posting i solved it myself but again thanks for your early reply.

サインインしてコメントする。

### 回答 (1 件)

Image Analyst 2022 年 6 月 25 日
What exactly is the end to end distance? Do you mean the Feret diameter bwferet?
Do you mean the sum of all white pixels to the neighbor pixels? Have you tried pdist2, and extract all pairs just 1 or sqrt(2) pixels apart, sum them, and divide by 2? Something like (untested):
[y, x] = find(k); % Find coordinates.
xy = [x(:), y(:)]; % Combine into matrix.
distances = pdist2(xy, xy); % Find distance of every point to every other point.
mask = distances < mean([sqrt(2), sqrt(3)]); % Find pixels that are adjacent or on nearest diagonal ONLY.
sumOfDistances = sum(distances(mask)) / 2 % Divide by 2 since we get the sum for each endpoint so it's counted twice.
##### 9 件のコメント表示非表示 8 件の古いコメント
Adarsh Tripathi 2022 年 6 月 26 日
@walter Roberson That is significant because changing of that distance will give the activity of the worm with time.

サインインしてコメントする。

R2022a

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!