Hexagonal tiling of arbitrary geoJSON

[EllAchE]

I/my company have a use case to fully cover a polygon of arbitrary shape with an optimal (i.e. minimal) number of fixed radius circles. I don't believe this is supported in current version of turf or other libraries (would love to know if it is), however we have put together a quick solution that takes advantage of bounding boxes and existing hexagonal tiling logic to achieve this.

Here's a visual of what we have right now https://replit.com/@kashifGhafoor/geoHexCoords-Visualization.

The issue is that this solution isn't fully optimized, doesn't capture some edge cases and runs into issues when the bounding box becomes too large. We also might want to support variable radii etc.

I wanted to know if anyone has taken a look into something like this or might be interested! I'm happy to discuss my thoughts/approach to solve the problem, and possibly in the future would be open to contributing such a solution.