Function euclidean_cluster_roots 

pub fn euclidean_cluster_roots(
    x: &[f32],
    y: &[f32],
    z: &[f32],
    cluster_tolerance: f32,
) -> SpatialResult<Vec<u32>>

Connected-component roots via uniform-grid union-find (minimum index per component).