Generate hexagons over areas not entirely within a polygon?

[marklit]

I want to generate H3s for the coastal areas as well as anything inland. With h3.LatLngPoly it only wants to generate hexagons that are entirely in land and miss areas like Miami for example. Any idea how I could do this?

Here's the code that I have atm:

$ wget https://raw.githubusercontent.com/scdoshi/us-geojson/master/geojson/nation/US.geojson

import geopandas as gpd
import h3

print(h3.versions()) # {'c': '4.0.0', 'python': '4.0.0b5'}

usa = gpd.read_file('US.geojson')
polygons = list(list(x.geoms) for x in usa.geometry)[0]

polygons.sort(key=lambda x: x.area)
coords = polygons[-1].exterior.coords[:]
coords = list(c[:2][::-1] for c in coords)

poly = h3.LatLngPoly(coords)

h3_list = set(h3.h3shape_to_cells(poly, res=3))

[grim7reaper]

I think it's now possible to do that with the H3 C library, but it's not yet exposed in h3-py IIUC.

Meanwhile you could try h3ronpy, with ContainmentMode=IntersectsBoundary (or Cover) it should do the trick. Cf. example here.

If it's a one-shot thing, you could also use h3o-cli:

h3o-cli geomToCells -r 3 -m intersects-boundary  < US.geojson | h3o-cli cellToPolygon > H3-coverage.geojson

[marklit]

Thanks for that.