radius_search: Radius Search of Points

In rTLS: Processing of Point Clouds Derived from Terrestrial Laser Scanning

Radius Search of Points

Description

Adapted radius searching of points based on RcppHNSW

Usage

radius_search(
  query,
  ref,
  radius,
  max_neighbour,
  distance = "euclidean",
  same = FALSE,
  threads = 1L,
  verbose = FALSE,
  progress = FALSE,
  ...
)

Arguments

query	A data.table containing the set of query points where each row represent a point and each column a given coordinate.
ref	A numeric containing the set of reference points where each row represent a point and each column a given coordinate.
radius	A numeric describing maximum euclidean distance form the each query points in which a point can be consider a neighbor.
max_neighbour	An integer specifying the maximum number of ref points to look around to consider for a given radius.
distance	Type of distance to calculate. "euclidean" as default. Look hnsw_knn for more options.
same	Logic. If TRUE, it delete neighbors with distance of 0, useful when the k search is based on the same query.
threads	An integer specifying the number of threads to use for parallel processing. Experiment to see what works best for your data on your hardware.
verbose	If TRUE, log messages to the console.
progress	If TRUE, log a progress bar when verbose = TRUE. Tracking progress could cause a small overhead.
...	Arguments passed to hnsw_build and hnsw_search.

Details

This function is based on hnswlib C++ library (Malkov & Yashunin 2016) and its bindings for R (RcppHNSW; Melville 2020) for a fast estimation of neighbors points. It is adapted to simplify the workflow within rTLS. If you use this function, please consider cite the C++ library and RcppHNSW package.

Value

A data.table with three columns describing the indices of the query and ref points and the distances.

Author(s)

J. Antonio Guzmán Q.

References

Malkov, Y. A., & Yashunin, D. A. (2016). Efficient and robust approximate nearest neighbor search using Hierarchical Navigable Small World graphs. arXiv preprint arXiv:1603.09320.

See Also

radius_search

Examples

#Point cloud
data("pc_tree")


#Radius search of 1
radius_search(pc_tree, pc_tree, radius = 1, max_neighbour = 100)

rTLS documentation built on Dec. 28, 2022, 1:28 a.m.