The proposal looks really good. One question I have is whether the hexagonal grid is essential -- would this work on a rectangular grid? Assuming that you replicate their model on a hexagonal grid, you might look around for suggestions on implementing a hexagonal grid in numpy and running a 6-neighbor kernel with correlate2d.
The proposal looks really good. One question I have is whether the hexagonal grid is essential -- would this work on a rectangular grid? Assuming that you replicate their model on a hexagonal grid, you might look around for suggestions on implementing a hexagonal grid in numpy and running a 6-neighbor kernel with correlate2d.