“I want to compare Chinese county-level population from 2000 to 2020, but the number of county units differs between years — counting by today’s county boundaries would be misaligned.”
- Download the China package (
data/countries/CHN.zip, about 1.9 GB); unzipping givesCHN_adm3_timeseries.parquetand friends — the bundledREADME.txtlists row counts and fields. - Select the year you need via
snapshot_year, then join your statistical table byunit_id. - Export a GPKG (opens directly in QGIS/ArcGIS) or keep working in Python.
import geopandas as gpd, pandas as pd
chn = gpd.read_parquet("CHN_adm3_timeseries.parquet")
b_2020 = chn[chn.snapshot_year == 2020][["unit_id", "name", "geometry"]]
pop = pd.read_csv("my_pop_2020.csv")
gdf = b_2020.merge(pop, on="unit_id", how="left")
gdf.to_file("chn_2020_with_pop.gpkg", driver="GPKG")