In a particular class of detectors, our data comes out as pairs of points in two dimensions, and we want to string these points into lines.
The data is noisy, and is binned in one direction but not in the other. We can't guarantee a hit in every bin even when each detector element is working, so there may be skips.
Our current analysis chain looks like
- Adjust hits for the calibration of individual detector elements
- Find clusters
- Rough fit lines to the clusters
- Connect up clusters into longer line-like structures
- ...
This question concerns step (3).
We've been using a Hough transform for that step and it works well, but as we try to scale up from the test-bed to simulation of a full-scale project it becomes unacceptably slow.
I'm looking for a faster way.
For those who might care the actual use case here is a Liquid Argon Time-Projection Chamber