A Python library for making geographic flat-top hexagon grids like QGIS's `create grid` function
Project description
A tiny Python 3.9+ library for making geographic flat-top hexagonal grids like QGIS’s create grid function. That’s it. Not designed for making other kinds of grids or discrete global grid systems.
Here a hexagonal grid, or hex grid for 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, see the Wikipedia article on hexagonal tilings.) The circumradius of 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 are
Making 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 is
Making any other kind of grid, e.g. ones with pointy-top hexagons, squares, triangles, kisrhombilles, Penrose tiles…
Here’s an typical example.
import geopandas as gpd
import geohexgrid as ghg
# 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 kilometres
grid = ghg.make_grid_from_gdf(nz, R=10_000)
# Plot
base = nz.plot(color="black", figsize=(20, 20), aspect="equal")
grid.plot(ax=base, color="white", edgecolor="red", alpha=0.5)
But why hexagons?! Because hexagons 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.
Installation
Install from PyPI, e.g. via poetry add geohexgrid.
Examples
See the Jupyter notebook at notebooks/examples.ipynb.
Notes
This 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 to MRCagney for periodically funding this project.
Red Blog Games has a great write up of hexagonal grids for 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.
Changes
2.0.0, 2023-11-14
Refactored for simpler architecture, gapless grids, and a ~15x speed up in the main function grid_from_gdf.
1.1.0, 2023-10-27
Added the clip option to the function grid_from_gdf.
Updated dependencies.
Re-ordered functions.
Changed the cell ID separotor to a comma.
1.0.0, 2022-08-15
First release.
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for geohexgrid-2.0.0-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 1c22eef70e1ac759c3753391eaccb1d44ec2e810b2e85d0edcb92f9e314a75e6 |
|
MD5 | 5a3331e21d07299c8195a419b9b3ff75 |
|
BLAKE2b-256 | 3df2ed93b6b3bb0aa9d4454a96d54ef7dc115d23eaff311f9cf3aaf86c2ca7d2 |