Commit 4c0e7240 authored by Chris Jewell's avatar Chris Jewell
Browse files

Merge branch 'add_inferencedata_format' into 'master'

Fix geopackage geometry joining issue

See merge request !31
parents da9aa2d7 cb23b18a
"""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
......
Markdown is supported
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