Discover MakerZone

MATLAB and Simulink resources for Arduino, LEGO, and Raspberry Pi

Learn more

Discover what MATLAB® can do for your career.

Opportunities for recent engineering grads.

Apply Today

New to MATLAB?

Remove the spurious edge of skeleton

Asked by Hong


on 26 Sep 2013
Latest activity Answered by Hong


on 30 Sep 2013
Accepted Answer by Teja Muppirala

Teja Muppirala

Recently, I've had an object to get it's skeleton,and then simplified the skeleton

but my problem is how to remove the edge of the skeleton for simplified the skeleton

here is my few steps:

First,I've got the skeleton from my object with bwmorph(ima,'skel',Inf)

here is the object which has been filled and the skeleton below:

filling the object:

the skeleton:

Second,for the skeleton, I've mark the branch points and end points

to mark the branch points and end points,I use bwmorph(skel, 'branchpoints') and

bwmorph(skel, 'endpoints');

here is the figure that I marked the points below:

the reds are branch points

the blues are end points

and the last step,

from the above figure, I want to remove the edges between branch point and end point (here I call the edge is spurious edge)

here is the resulting figure that I want

my problem is how can I remove the spurious edges

I've referenced this article

and the suggestions are set the threshold of the distance

but in my situation,the lengths of the spurious edges may have too long or too short...





No products are associated with this question.

2 Answers

Answer by Teja Muppirala

Teja Muppirala

on 26 Sep 2013
Accepted answer

I think there must be a more efficient way to do it, but this at least works. BWDISTGEODESIC starts at a given point, and then calculates the distance as you travel along the path. Start at a endpoint, start walking and find all pixels that are closer than the nearest branchpoint. Then remove those pixels.

skel= bwmorph(ima,'skel',Inf);
B = bwmorph(skel, 'branchpoints');
E = bwmorph(skel, 'endpoints');
[y,x] = find(E);
B_loc = find(B);
Dmask = false(size(skel));
for k = 1:numel(x)
    D = bwdistgeodesic(skel,x(k),y(k));
    distanceToBranchPt = min(D(B_loc));
    Dmask(D < distanceToBranchPt) =true;
skelD = skel - Dmask;
hold all;
[y,x] = find(B); plot(x,y,'ro')


Teja Muppirala

Teja Muppirala

Answer by Hong


on 30 Sep 2013

thanks to Teja Muppirala

this solution is really help me a lot

as your method, it could worked ^^

thanks to you again :")




Contact us