点云-半径搜索法-Radius Search

• 核心作用
  • 在于通过设定一个==空间范围（半径）==寻找点的邻域点集合，从而支持对局部区域的分析和操作。
• 因为空间半径不会随着密度变化而改变点云输出的结果，处理密度变化大的点云时很重要。

应用场景

1. 稀疏点检测：当点云密度不均匀时，半径搜索可以检测稀疏区域。
2. 局部平滑：利用半径内点的均值更新点坐标，去除噪声。
3. 体素滤波 (Voxel Filter)：基于半径选取中心点，实现点云降采样。

场景中的作用

1. 邻域点搜索：
   • 用于寻找查询点附近的所有点，可以用于密度估计、曲率计算、法向量估计等。
2. 点云降噪与分割：
   • 半径搜索可以根据邻域密度识别孤立点（噪声）或区域边界。
3. 稀疏与下采样：
   • 在点云简化过程中，通过设置固定半径找到邻域点，以均匀选取代表性点。
4. 点云配准与滤波：
   • 配准算法（如ICP）中，可以用半径搜索代替KNN搜索，找到距离目标点云一定范围内的点。

计算

半径搜索基于一个查询点 ( p )，在给定的半径 ( r ) 内找到所有满足条件的点 ( p_i )。具体计算如下：

1. 距离定义：

   • 同样使用欧几里得距离：$d(p,p_i)=\sqrt{(x-x_i{)}^2+(y-y_i{)}^2+(z-z_i{)}^2}$
   • 满足$d(p,p_i)\le r$ 的点 ( p_i ) 被视为邻居点。

2. 搜索方式：

   • Brute Force：
     • 对所有点计算距离，保留距离小于 ( r ) 的点，复杂度为 ( O(n) )。
   • 加速方法：
     • 使用 kd-tree、Octree 等空间划分结构，搜索效率提升至 ( O(\log n) )。

3. 半径范围结果：

   • 搜索结果为一个点集合，表示所有邻域点的索引或坐标。

实现框架

实现框架实现：PCL、Open3D、scikit-learn

import open3d as o3d
import numpy as np# 生成随机点云
points = np.random.rand(1000, 3)  # 1000个3D点
query_point = points[0]  # 查询点
radius = 0.1  # 搜索半径# 创建点云对象
pcd = o3d.geometry.PointCloud()
pcd.points = o3d.utility.Vector3dVector(points)# 构建kd-tree
pcd_tree = o3d.geometry.KDTreeFlann(pcd)# 半径搜索
indices = pcd_tree.search_radius_vector_3d(query_point, radius)[1]print("查询点：", query_point)
print("半径邻域点索引：", indices)
print("邻域点数量：", len(indices))