in data/geospatial/crs-gen.py [0:0]
def check_crs_crs(name, expected_crs):
expected_crs = pyproj.CRS(expected_crs)
file = parquet.ParquetFile(HERE / name, arrow_extensions_enabled=True)
ext_type = file.schema_arrow.field(1).type
actual_crs = pyproj.CRS(ext_type.crs)
if actual_crs != expected_crs:
raise ValueError(f"Expected '{expected_crs}' crs but got '{actual_crs}'")