Package: geohexgrid
A tiny Python 3.9+ library for making geographic flat-top hexagonal grids like QGIS’screate grid function.
That’s it.
Not designed for making other kinds of grids ordiscrete global grid systems.Here ahexagonal grid, orhex gridfor short, is a finite subset of a hexagonal tiling.
A hexagonal tiling is a covering of the plane with regular hexagons in which exactly three hexagons meet at each vertex.
(For more details, seethe Wikipedia article on hexagonal tilings.)
Thecircumradiusof a hex grid is the circumradius of any one of its hexagons, that is, the radius of a circle circumscribing any one of the hexagons.
This library favours the word ‘grid’ over ‘tiling’, because ‘grid’ is used more often in geographic context, the main context of this library.The two main features of this library areMaking a flat-top hexagonal grid of given circumradius that minimally covers a GeoDataFrame of features, where distance units come from the GeoDataFrame’s coordinate reference system (CRS), e.g. no units for no CRS, metres for the New Zealand Transverse Mercator (NZTM) CRS, and decimal degrees for the WGS84 CRS.By default, hex grids made with a common CRS and circumradis share an origin and thus have equal hexagons (and hexagon IDs) where they overlap.
In other words, the grids share a single (infinite) hexagonal tiling of the plane, which is useful when reconciling multiple grids across different geographic areas.The main non-feature of this library isMaking any other kind of grid, e.g. ones with pointy-top hexagons, squares, triangles, kisrhombilles, Penrose tiles…Here’s an typical example.importgeopandasasgpdimportgeohexgridasghg# Load New Zealand territorial authorities projected in EPSG 2193 (NZTM)nz=gpd.read_file(DATA_DIR/"nz_tas.gpkg")# Cover it minimally with hexagons of circumradius 10 kilometresgrid=ghg.make_grid_from_gdf(nz,R=10_000)# Plotbase=nz.plot(color="black",figsize=(20,20),aspect="equal")grid.plot(ax=base,color="white",edgecolor="red",alpha=0.5)But why hexagons?!
Becausehexagons are the bestagons.
More seriously, no one grid type works best for all geographic applications.
MRCagney, this library’s funder, often works with isochrones, which favour simple convex equal area grids with equidistant neighbour cells, that is, hex gids.AuthorsAlex Raichev (2014-09), maintainerInstallationInstall from PyPI, e.g. viapoetry add geohexgrid.ExamplesSee the Jupyter notebook atnotebooks/examples.ipynb.NotesThis project’s development status is Alpha.
Alex uses this project for work and changes it breakingly when it suits his needs.This project uses semantic versioning.Thanks toMRCagneyfor periodically funding this project.Red Blog Games has agreat write up of hexagonal gridsfor computer games.Alex wanted to chose a shorter name for this package, such as ‘hexgrid’, ‘geohex’, or ‘hexcover’, but those were already taken or too close to taken on PyPI.Changes2.0.0, 2023-11-14Refactored for simpler architecture, gapless grids, and a ~15x speed up in the main functiongrid_from_gdf.1.1.0, 2023-10-27Added theclipoption to the functiongrid_from_gdf.Updated dependencies.Re-ordered functions.Changed the cell ID separotor to a comma.1.0.0, 2022-08-15First release.