Function grid_bounds 

pub fn grid_bounds(
    x: &[f32],
    y: &[f32],
    z: &[f32],
    cell_size: f32,
) -> SpatialResult<([f32; 3], [u32; 3])>

Returns grid origin (min corner) and cell counts for cell size radius.