Suggestion: Adding support for unstructured grids about scikit-fmm HOT 1 CLOSED

Hi, Thanks for the message.

Scikit-fmm could be adapted to work on non-uniform grids. The basics of the algorithm would be unchanged but the quadratic equation solved to update each point would be different and a data structure to hold the grid would be needed. Code-wise we would create a new sub-class of the BaseMarcher class that has new updatePointOrderOne and updatePointOrderTwo functions. A few other things, like initializing the initially frozen points and the narrow band, may also need to change.

Having the spatial gradients pre-defined at each point would probably not make it much easier. The gradient has to be calculated in a specific upwind way for the algorithm to work.

The Sethian book has the formulation for the fast marching method on unstructured grids which we could implement.

Several other people have asked about this in the past. Some people have been able to map their data to a regular grid for fast-marching as a work around. I do not have any plans to add unstructured grids any time soon but I will keep it in mind. If you want to start developing the feature I can help.

from scikit-fmm.