eaf30a9b65
* rewrite mocksnakemake for parsing real Snakefile * continue add function to scripts * going through all scripts, setting new mocksnakemake * fix plotting scripts * fix build_country_flh * fix build_country_flh II * adjust config files * fix make_summary for tutorial network * create dir also for output * incorporate suggestions * consistent import of mocksnakemake * consistent import of mocksnakemake II * Update scripts/_helpers.py Co-Authored-By: euronion <42553970+euronion@users.noreply.github.com> * Update scripts/_helpers.py Co-Authored-By: euronion <42553970+euronion@users.noreply.github.com> * Update scripts/_helpers.py Co-Authored-By: euronion <42553970+euronion@users.noreply.github.com> * Update scripts/_helpers.py Co-Authored-By: euronion <42553970+euronion@users.noreply.github.com> * Update scripts/plot_network.py Co-Authored-By: euronion <42553970+euronion@users.noreply.github.com> * Update scripts/plot_network.py Co-Authored-By: euronion <42553970+euronion@users.noreply.github.com> * Update scripts/retrieve_databundle.py Co-Authored-By: euronion <42553970+euronion@users.noreply.github.com> * use pathlib for mocksnakemake * rename mocksnakemake into mock_snakemake * revert change in data * Update scripts/_helpers.py Co-Authored-By: euronion <42553970+euronion@users.noreply.github.com> * remove setting logfile in mock_snakemake, use Path in configure_logging * fix fallback path and base_dir fix return type of make_io_accessable * reformulate mock_snakemake * incorporate suggestion, fix typos * mock_snakemake: apply absolute paths again, add assertion error *.py: make hard coded io path accessable for mock_snakemake * retrieve_natura_raster: use snakemake.output for fn_out * include suggestion * Apply suggestions from code review Co-Authored-By: Jonas Hörsch <jonas.hoersch@posteo.de> * linting, add return ad end of file * Update scripts/plot_p_nom_max.py Co-Authored-By: Jonas Hörsch <jonas.hoersch@posteo.de> * Update scripts/plot_p_nom_max.py fixes #112 Co-Authored-By: Jonas Hörsch <jonas.hoersch@posteo.de> * plot_p_nom_max: small correction * config.tutorial.yaml fix snapshots end * use techs instead of technology * revert try out from previous commit, complete replacing * change clusters -> clusts in plot_p_nom_max due to wildcard constraints of clusters * change clusters -> clusts in plot_p_nom_max due to wildcard constraints of clusters II
104 lines
2.9 KiB
Python
104 lines
2.9 KiB
Python
"""
|
|
Creates Voronoi shapes for each bus representing both onshore and offshore regions.
|
|
|
|
Relevant Settings
|
|
-----------------
|
|
|
|
.. code:: yaml
|
|
|
|
countries:
|
|
|
|
.. seealso::
|
|
Documentation of the configuration file ``config.yaml`` at
|
|
:ref:`toplevel_cf`
|
|
|
|
Inputs
|
|
------
|
|
|
|
- ``resources/country_shapes.geojson``: confer :ref:`shapes`
|
|
- ``resources/offshore_shapes.geojson``: confer :ref:`shapes`
|
|
- ``networks/base.nc``: confer :ref:`base`
|
|
|
|
Outputs
|
|
-------
|
|
|
|
- ``resources/regions_onshore.geojson``:
|
|
|
|
.. image:: ../img/regions_onshore.png
|
|
:scale: 33 %
|
|
|
|
- ``resources/regions_offshore.geojson``:
|
|
|
|
.. image:: ../img/regions_offshore.png
|
|
:scale: 33 %
|
|
|
|
Description
|
|
-----------
|
|
|
|
"""
|
|
|
|
import logging
|
|
logger = logging.getLogger(__name__)
|
|
from _helpers import configure_logging
|
|
|
|
from vresutils.graph import voronoi_partition_pts
|
|
|
|
import os
|
|
|
|
import pandas as pd
|
|
import geopandas as gpd
|
|
|
|
import pypsa
|
|
|
|
if __name__ == "__main__":
|
|
if 'snakemake' not in globals():
|
|
from _helpers import mock_snakemake
|
|
snakemake = mock_snakemake('build_bus_regions')
|
|
configure_logging(snakemake)
|
|
|
|
countries = snakemake.config['countries']
|
|
|
|
n = pypsa.Network(snakemake.input.base_network)
|
|
|
|
country_shapes = gpd.read_file(snakemake.input.country_shapes).set_index('name')['geometry']
|
|
offshore_shapes = gpd.read_file(snakemake.input.offshore_shapes).set_index('name')['geometry']
|
|
|
|
onshore_regions = []
|
|
offshore_regions = []
|
|
|
|
for country in countries:
|
|
c_b = n.buses.country == country
|
|
|
|
onshore_shape = country_shapes[country]
|
|
onshore_locs = n.buses.loc[c_b & n.buses.substation_lv, ["x", "y"]]
|
|
onshore_regions.append(gpd.GeoDataFrame({
|
|
'name': onshore_locs.index,
|
|
'x': onshore_locs['x'],
|
|
'y': onshore_locs['y'],
|
|
'geometry': voronoi_partition_pts(onshore_locs.values, onshore_shape),
|
|
'country': country
|
|
}))
|
|
|
|
if country not in offshore_shapes.index: continue
|
|
offshore_shape = offshore_shapes[country]
|
|
offshore_locs = n.buses.loc[c_b & n.buses.substation_off, ["x", "y"]]
|
|
offshore_regions_c = gpd.GeoDataFrame({
|
|
'name': offshore_locs.index,
|
|
'x': offshore_locs['x'],
|
|
'y': offshore_locs['y'],
|
|
'geometry': voronoi_partition_pts(offshore_locs.values, offshore_shape),
|
|
'country': country
|
|
}, index=offshore_locs.index)
|
|
offshore_regions_c = offshore_regions_c.loc[offshore_regions_c.area > 1e-2]
|
|
offshore_regions.append(offshore_regions_c)
|
|
|
|
def save_to_geojson(s, fn):
|
|
if os.path.exists(fn):
|
|
os.unlink(fn)
|
|
schema = {**gpd.io.file.infer_schema(s), 'geometry': 'Unknown'}
|
|
s.to_file(fn, driver='GeoJSON', schema=schema)
|
|
|
|
save_to_geojson(pd.concat(onshore_regions), snakemake.output.regions_onshore)
|
|
|
|
save_to_geojson(pd.concat(offshore_regions), snakemake.output.regions_offshore)
|