How to segment/partition point cloud into upper and lower layers?

Question

Di Zhu 2017-4-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/333846-how-to-segment-partition-point-cloud-into-upper-and-lower-layers

评论： Image Analyst 2022-3-7

This problem has been troubling me for long. I have attached a dataset of unorganized point cloud. It is the partial surface of a cartilage, and it's clear that it has upper and lower surfaces. Now I'd like to use Matlab to divide the point cloud into upper and lower parts, any idea how I could do this? I thought about local search method grid by grid then separating them from the mid-plane. But this method seems very limited and doesn't work well for dataset with a severe slope.

Alternatively, I could look to segment a contour slice of it into upper and lower part. I just need a rough segmentation, the left and right end points don't matter much. After segmentation, I'd do a B-Spline surface approximation.

Any ideas? Thanks in advance!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Image Analyst 2017-4-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/333846-how-to-segment-partition-point-cloud-into-upper-and-lower-layers#answer_261893

在 MATLAB Online 中打开

Assuming you have (x,y,z) coordinates and want to find points above and below some specified z value:

zDivider = -41; % Dividing line between upper and lower points.
lowerPointsIndexes = z < zDivider; % Logical index.
% Now extract out x, y, and z for the lower part
xLower = x(lowerPointsIndexes);
yLower = y(lowerPointsIndexes);
zLower = z(lowerPointsIndexes);
% Now do for the upper points:
upperPointsIndexes = z >= zDivider; % Logical index.
% Now extract out x, y, and z for the lower part
xUpper = x(upperPointsIndexes);
yUpper = y(upperPointsIndexes);
zUpper = z(upperPointsIndexes);

4 个评论
显示 2更早的评论隐藏 2更早的评论

stefano chiappini 2022-3-7

example crown slicing.jpg

Hi, i have a similar question. I have got a point cloud crown tree and i need slicing along z direction with with 10 cm thickness. Is there a way to create a list ? my aim is to compute the volume with alpha shape concern each singular slice.

Thank you so much

Image Analyst 2022-3-7

Not sure what you have. Do you have a 3-D volumetric image where the voxel value is 1 if a tree occupies that voxel? Do you just want to get the top most z value for every (x,y) coordinate? And then maybe fill in some "holes" where the top most z value is way down near the middle or bottom of the tree? If so, attach your variable in a .mat file.

请先登录，再进行评论。