Add weights parameter to density_from_points()

[Copilot]

Adds optional weights parameter to density_from_points() to support weighted density maps (e.g., weighting points by their volume in voxels or mm³).

Changes

• Add weights: Union[np.ndarray, None] = None parameter, passed through to da.histogramdd()
• Validate weights shape matches number of points
• Convert weights to dask array with compatible chunking for histogramdd

Usage

points = np.array([[5.0, 5.0, 5.0], [6.0, 6.0, 6.0]])
volumes = np.array([100.0, 250.0])  # e.g., nvoxels or mm³
weighted_density = density_from_points(points, ref_img, weights=volumes)

Tests

Added 7 tests covering weighted density computation, accumulation, voxel space mode, and error handling.

Original prompt

in the density_from_points() func, add an option to provide an optional weights array, passed along to histogramdd.. e.g. could use this to weight points by their size in nvoxels, (or mm^3)

---

💬 We'd love your input! Share your thoughts on Copilot coding agent in our 2 minute survey.