[SciPy-Dev] Contribution: Estimating the overlap of a ray with volume-cells (voxels) for e.g. tomography

[Sebastian A. Mueller]

Dear scipy,

I wrote a package to estimate the distances a ray overlaps with 
volume-cells in three dimensions. This is needed to estimate 
system-matrices in tomography. The interface is only one function. It 
takes the ray's support and direction vectors, and the bin-edges for the 
space partitioning. It returns the overlaps of the ray and the voxels. 
The space-partitioning must be rectangular, but it does not need to be 
equally spaced. In combination with scipy.sparse, it can efficiently 
create system-matrices for iterative, tomographic reconstructions. The 
algorithm uses oct-tree-traversal for efficiency.

On pypi:
https://pypi.org/project/ray-voxel-overlap/0.0.2/

On github:
https://github.com/cherenkov-plenoscope/ray_voxel_overlap

I would like to contribute this to scipy. It addresses multiple fields 
in science such as medicine, biology, microscopy, physics, astronomy, 
geo-earth, image-processing.
At its core, it is pure linear algebra. My package contains a pure 
python-implementation for development and testing, and a (much faster) 
C-implementation wrapped in cython.

Do you think such a tool belongs into scipy? If not, where else should 
it go?
This is my first attempt to contribute to scipy.

Cheers
Sebastian