mrpro.algorithms.dcf.dcf_voronoi.Voronoi
- class mrpro.algorithms.dcf.dcf_voronoi.Voronoi(points, furthest_site=False, incremental=False, qhull_options=None)
Bases:
_QhullUser
Voronoi diagrams in N dimensions.
Added in version 0.12.0.
- Parameters:
points (ndarray of floats, shape (npoints, ndim)) – Coordinates of points to construct a Voronoi diagram from
furthest_site (bool, optional) – Whether to compute a furthest-site Voronoi diagram. Default: False
incremental (bool, optional) – Allow adding new points incrementally. This takes up some additional resources.
qhull_options (str, optional) – Additional options to pass to Qhull. See Qhull manual for details. (Default: “Qbb Qc Qz Qx” for ndim > 4 and “Qbb Qc Qz” otherwise. Incremental mode omits “Qz”.)
- points
Coordinates of input points.
- Type:
ndarray of double, shape (npoints, ndim)
- vertices
Coordinates of the Voronoi vertices.
- Type:
ndarray of double, shape (nvertices, ndim)
- ridge_points
Indices of the points between which each Voronoi ridge lies.
- Type:
ndarray of ints, shape
(nridges, 2)
- ridge_vertices
Indices of the Voronoi vertices forming each Voronoi ridge.
- Type:
list of list of ints, shape
(nridges, *)
- regions
Indices of the Voronoi vertices forming each Voronoi region. -1 indicates vertex outside the Voronoi diagram. When qhull option “Qz” was specified, an empty sublist represents the Voronoi region for a point at infinity that was added internally.
- Type:
list of list of ints, shape
(nregions, *)
- point_region
Index of the Voronoi region for each input point. If qhull option “Qc” was not specified, the list will contain -1 for points that are not associated with a Voronoi region. If qhull option “Qz” was specified, there will be one less element than the number of regions because an extra point at infinity is added internally to facilitate computation.
- Type:
array of ints, shape (npoints)
- furthest_site
True if this was a furthest site triangulation and False if not.
Added in version 1.4.0.
- Raises:
QhullError – Raised when Qhull encounters an error condition, such as geometrical degeneracy when options to resolve are not enabled.
ValueError – Raised if an incompatible array is given as input.
Notes
The Voronoi diagram is computed using the Qhull library.
Examples
Voronoi diagram for a set of point:
>>> import numpy as np >>> points = np.array([[0, 0], [0, 1], [0, 2], [1, 0], [1, 1], [1, 2], ... [2, 0], [2, 1], [2, 2]]) >>> from scipy.spatial import Voronoi, voronoi_plot_2d >>> vor = Voronoi(points)
Plot it:
>>> import matplotlib.pyplot as plt >>> fig = voronoi_plot_2d(vor) >>> plt.show()
The Voronoi vertices:
>>> vor.vertices array([[0.5, 0.5], [0.5, 1.5], [1.5, 0.5], [1.5, 1.5]])
There is a single finite Voronoi region, and four finite Voronoi ridges:
>>> vor.regions [[], [-1, 0], [-1, 1], [1, -1, 0], [3, -1, 2], [-1, 3], [-1, 2], [0, 1, 3, 2], [2, -1, 0], [3, -1, 1]] >>> vor.ridge_vertices [[-1, 0], [-1, 0], [-1, 1], [-1, 1], [0, 1], [-1, 3], [-1, 2], [2, 3], [-1, 3], [-1, 2], [1, 3], [0, 2]]
The ridges are perpendicular between lines drawn between the following input points:
>>> vor.ridge_points array([[0, 3], [0, 1], [2, 5], [2, 1], [1, 4], [7, 8], [7, 6], [7, 4], [8, 5], [6, 3], [4, 5], [4, 3]], dtype=int32)
- __init__(points, furthest_site=False, incremental=False, qhull_options=None)
Methods
__init__
(points[, furthest_site, ...])add_points
(points[, restart])Process a set of additional new points.
close
()Finish incremental processing.
Attributes
ridge_dict
- add_points(points, restart=False)
Process a set of additional new points.
- Parameters:
points (ndarray) – New points to add. The dimensionality should match that of the initial points.
restart (bool, optional) – Whether to restart processing from scratch, rather than adding points incrementally.
- Raises:
QhullError – Raised when Qhull encounters an error condition, such as geometrical degeneracy when options to resolve are not enabled.
See also
Notes
You need to specify
incremental=True
when constructing the object to be able to add points incrementally. Incremental addition of points is also not possible after close has been called.
- close()
Finish incremental processing.
Call this to free resources taken up by Qhull, when using the incremental mode. After calling this, adding more points is no longer possible.