Avoiding 3rd degree skeletons and calculating the volume of the 2rd degree branch skeleton

Hello all,

Here attached is a .csv file containing the voxel information. I converted it to 3D binary array and used Skeleton3d.m (https://nl.mathworks.com/matlabcentral/fileexchange/43400-skeleton3d?s_tid=FX_rc1_behav) to skeletonize. The purpose of skeletonizing is to determine where the branch starts and find the area and volume of that entire segment corresponding to that branch in 3D. In that process, I notice there are many 3rd and 4th degree branches (like a branch that originates from a branch from the main stem).

Any suggestion on how to avoid the 3rd degree branch and calculate the volume of that not just the branch but the entire segment of that particular 3D structure would be helpful.

Thanks!