Commit cb23b18a authored by Chris Jewell's avatar Chris Jewell
Browse files

Fix geopackage geometry joining issue

parent 023cbc5f
"""Summarises posterior distribution into a geopackage"""
import pickle as pkl
import numpy as np
import xarray
import pandas as pd
import geopandas as gp
......@@ -31,7 +32,7 @@ def summary_geopackage(input_files, output_file, config):
# Load and filter geopackage
geo = gp.read_file(config["base_geopackage"], layer=config["base_layer"])
geo = geo[geo["lad19cd"].isin(data.coords["location"])]
geo = geo[geo["lad19cd"].isin(np.array(data.coords["location"]))]
geo = geo.sort_values(by="lad19cd")
# Dump data into the geopackage
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment