File Exchange

image thumbnail

splitFV - split a mesh

version (2.84 KB) by Sven
Splits a 2D or 3D mesh defined by faces and vertices into separately connected mesh pieces.


Updated 24 Jun 2014

View License

FVOUT = SPLITFV(F,V) separates disconnected pieces inside a patch defined by faces (F) and
vertices (V). FVOUT is a structure array with fields "faces" and "vertices". Each element of
this array indicates a separately connected patch.
FVOUT = SPLITFV(FV) takes in FV as a structure with fields "faces" and "vertices"

For example:
fullpatch.vertices = [2 4; 2 8; 8 4; 8 0; 0 4; 2 6; 2 2; 4 2; 4 0; 5 2; 5 0];
fullpatch.faces = [1 2 3; 1 3 4; 5 6 1; 7 8 9; 11 10 4];
figure, subplot(2,1,1), patch(fullpatch,'facecolor','r'), title('Unsplit mesh');
splitpatch = splitFV(fullpatch);
colours = lines(length(splitpatch));
subplot(2,1,2), hold on, title('Split mesh');
for i=1:length(splitpatch)

Note: faces and vertices should be defined such that faces sharing a coincident vertex reference
the same vertex number, rather than having a separate vertice defined for each face (yet at the
same vertex location). In other words, running the following command: size(unique(v,'rows') ==
size(v) should return TRUE. An explicit test for this has not been included in this function so
as to allow for the deliberate splitting of a mesh at a given location by simply duplicating
those vertices.

See also PATCH

Cite As

Sven (2020). splitFV - split a mesh (, MATLAB Central File Exchange. Retrieved .

Comments and Ratings (13)

sana3 sal

Hello Sir, Can you please help me understanding the code?




Christopher, good spot and I agree completely - the xyzSub variable in the previous was not doing as intended. It would work in most cases simply because meshes are usually over-connected so the bug wouldn't come up. I've updated as per your change which should now be correct and a little faster.


This implementation is slightly flawed. It does not work on the face matrix
5 9 10
5 4 9
4 8 9
4 3 8
3 7 8
3 2 7
2 6 7
2 1 6

It erroneously splits this mesh into two pieces.

It seems that the lines
[availFaceSub, xyzSub] = find(ismember(f(availFaceInds,:), openVertices));
fSets(availFaceInds(availFaceSub)) = currentSet;
openVertices = f(availFaceInds(availFaceSub),xyzSub);

should be replaced by
[availFaceSub, ~] = find(ismember(f(availFaceInds,:), openVertices));
fSets(availFaceInds(availFaceSub)) = currentSet;
openVertices = f(availFaceInds(availFaceSub),:);

xyzSub is unneeded and f(availFaceInds(availFaceSub),xyzSub) is probably not working as the author intended (I'm not sure what the author intended here)---xyzSub is basically replicating particular columns of f(availFaceInds(availFaceSub),:). Thus, if xyzSub doesn't contain 1, 2, or 3, there may be some vertices missing from the openVertices list.



Thanks for the images Tom :)

This did exactly what I wanted, straight out of the box, without editing a single line of code.

That's my kind of FEX submission!

NB For users wondering about this, it works in 3D too.

Sven, would you like some JPEGs of it decomposing different vortices in a fluid flow?


Implemented Christopher's amendment to fix bug in following connectivity matrix.

Updated icon (thanks Tom Clark!)

MATLAB Release Compatibility
Created with R2010a
Compatible with any release
Platform Compatibility
Windows macOS Linux