如何在Open3D中从初始点云移除已提取的子点云？

要从原pcd点云中排除ceiling和floor子点云，关键是拿到这两个子点云在原云里的索引，再求原云索引的补集，以下是两种实用方法：

方法一：分割时直接保存索引（推荐） #

如果你的ceiling和floor是通过索引筛选得到的（比如用segment_plane或select_by_index），直接保存当时的索引列表即可：

1. 分割天花板时保存索引：

ceiling_indices, _ = pcd.segment_plane(distance_threshold=0.01, ransac_n=3, num_iterations=1000)
ceiling = pcd.select_by_index(ceiling_indices)

2. 同理分割地面并保存floor_indices，然后合并两个索引列表：

combined_indices = ceiling_indices + floor_indices

3. 生成原云的全部索引，筛选出不在合并列表里的索引：

all_indices = list(range(len(pcd.points)))
remaining_indices = [idx for idx in all_indices if idx not in combined_indices]
# 获取剩余点云
remaining_pcd = pcd.select_by_index(remaining_indices)

方法二：通过坐标匹配获取索引（无原始索引时用） #

如果分割时没保存索引，可以通过点坐标匹配来查找，但要注意浮点精度问题：

import numpy as np

# 将点云数据转为numpy数组
pcd_points = np.asarray(pcd.points)
ceiling_points = np.asarray(ceiling.points)
floor_points = np.asarray(floor.points)

# 合并要排除的点
removed_points = np.vstack((ceiling_points, floor_points))

# 设置精度阈值，避免浮点误差导致匹配失败
threshold = 1e-6
# 标记原云中属于ceiling或floor的点
mask = np.all(np.isclose(pcd_points[:, None], removed_points, atol=threshold), axis=2).any(axis=1)

# 取反得到剩余点的索引
remaining_indices = np.where(~mask)[0].tolist()
remaining_pcd = pcd.select_by_index(remaining_indices)

注意事项 #

• 方法一精度更高，没有坐标匹配的误差，优先推荐在分割子点云时同步保存索引。
• 方法二中的threshold需要根据你的点云精度调整，过小可能漏匹配，过大可能误匹配。

内容的提问来源于stack exchange，提问作者user16570427