remove build/retrieve natura raster, directly use shapefile

Snakefile

@@ -165,28 +165,11 @@ if config['enable'].get('retrieve_cutout', True):
         run: move(input[0], output[0])
 
 
-if config['enable'].get('build_natura_raster', False):
-    rule build_natura_raster:
-        input:
-            natura="data/bundle/natura/Natura2000_end2015.shp",
-            cutouts=expand("cutouts/{cutouts}.nc", **config['atlite'])
-        output: "resources/natura.tiff"
-        log: "logs/build_natura_raster.log"
-        script: "scripts/build_natura_raster.py"
-
-
-if config['enable'].get('retrieve_natura_raster', True):
-    rule retrieve_natura_raster:
-        input: HTTP.remote("zenodo.org/record/4706686/files/natura.tiff", keep_local=True, static=True)
-        output: "resources/natura.tiff"
-        run: move(input[0], output[0])
-
-
 rule build_renewable_profiles:
     input:
         base_network="networks/base.nc",
         corine="data/bundle/corine/g250_clc06_V18_5.tif",
-        natura="resources/natura.tiff",
+        natura="data/bundle/natura/Natura2000_end2015.shp",
         gebco=lambda w: ("data/bundle/GEBCO_2014_2D.nc"
                          if "max_depth" in config["renewable"][w.technology].keys()
                          else []),

config.default.yaml

@@ -33,8 +33,6 @@ enable:
   retrieve_databundle: true
   build_cutout: false
   retrieve_cutout: true
-  build_natura_raster: false
-  retrieve_natura_raster: true
   custom_busmap: false
 
 electricity:

config.tutorial.yaml

@@ -33,8 +33,6 @@ enable:
   retrieve_databundle: true
   build_cutout: false
   retrieve_cutout: true
-  build_natura_raster: false
-  retrieve_natura_raster: true
   custom_busmap: false
 
 electricity:

doc/configtables/toplevel.csv

@@ -12,6 +12,4 @@ enable,,,
 -- retrieve_databundle,bool,"{true, false}","Switch to retrieve databundle from zenodo via the rule :mod:`retrieve_databundle` or whether to keep a custom databundle located in the corresponding folder."
 -- build_cutout,bool,"{true, false}","Switch to enable the building of cutouts via the rule :mod:`build_cutout`."
 -- retrieve_cutout,bool,"{true, false}","Switch to enable the retrieval of cutouts from zenodo with :mod:`retrieve_cutout`."
--- build_natura_raster,bool,"{true, false}","Switch to enable the creation of the raster ``natura.tiff`` via the rule :mod:`build_natura_raster`."
--- retrieve_natura_raster,bool,"{true, false}","Switch to enable the retrieval of ``natura.tiff`` from zenodo with :mod:`retrieve_natura_raster`."
 -- custom_busmap,bool,"{true, false}","Switch to enable the use of custom busmaps in rule :mod:`cluster_network`. If activated the rule looks for provided busmaps at ``data/custom_busmap_elec_s{simpl}_{clusters}.csv`` which should have the same format as ``resources/busmap_elec_s{simpl}_{clusters}.csv``, i.e. the index should contain the buses of ``networks/elec_s{simpl}.nc``."

doc/preparation.rst

@@ -27,7 +27,6 @@ With these and the externally extracted ENTSO-E online map topology
 Then the process continues by calculating conventional power plant capacities, potentials, and per-unit availability time series for variable renewable energy carriers and hydro power plants with the following rules:
 
 - :mod:`build_powerplants` for today's thermal power plant capacities using `powerplantmatching <https://github.com/FRESNA/powerplantmatching>`_ allocating these to the closest substation for each powerplant,
-- :mod:`build_natura_raster` for rasterising NATURA2000 natural protection areas,
 - :mod:`build_renewable_profiles` for the hourly capacity factors and installation potentials constrained by land-use in each substation's Voronoi cell for PV, onshore and offshore wind, and
 - :mod:`build_hydro_profile` for the hourly per-unit hydro power availability time series.
 
@@ -41,7 +40,6 @@ together into a detailed PyPSA network stored in ``networks/elec.nc``.
    preparation/build_shapes
    preparation/build_load_data
    preparation/build_cutout
-   preparation/build_natura_raster
    preparation/prepare_links_p_nom
    preparation/base_network
    preparation/build_bus_regions

doc/preparation/build_natura_raster.rst

@@ -1,39 +0,0 @@
-..
-  SPDX-FileCopyrightText: 2019-2020 The PyPSA-Eur Authors
-
-  SPDX-License-Identifier: CC-BY-4.0
-
-.. _natura:
-
-Rule ``build_natura_raster``
-===============================
-
-.. graphviz::
-    :align: center
-
-    digraph snakemake_dag {
-        graph [bgcolor=white,
-            margin=0,
-            size="8,5"
-        ];
-        node [fontname=sans,
-            fontsize=10,
-            penwidth=2,
-            shape=box,
-            style=rounded
-        ];
-        edge [color=grey,
-            penwidth=2
-        ];
-        9	 [color="0.22 0.6 0.85",
-            label=build_renewable_profiles];
-        12	 [color="0.31 0.6 0.85",
-            fillcolor=gray,
-            label=build_natura_raster,
-            style=filled];
-        12 -> 9;
-    }
-
-|
-
-.. automodule:: build_natura_raster

doc/preparation/build_renewable_profiles.rst

@@ -41,9 +41,6 @@ Rule ``build_renewable_profiles``
         8	 [color="0.00 0.6 0.85",
             label=build_shapes];
         8 -> 9;
-        12	 [color="0.31 0.6 0.85",
-            label=build_natura_raster];
-        12 -> 9;
         13	 [color="0.56 0.6 0.85",
             label=build_cutout];
         13 -> 9;

doc/preparation/retrieve.rst

@@ -50,30 +50,3 @@ The :ref:`tutorial` uses a smaller cutout than required for the full model (30 M
 
 .. seealso::
     For details see :mod:`build_cutout` and read the `atlite documentation <https://atlite.readthedocs.io>`_.
-
-
-Rule ``retrieve_natura_raster``
--------------------------------
-
-.. image:: https://zenodo.org/badge/DOI/10.5281/zenodo.4706686.svg
-   :target: https://doi.org/10.5281/zenodo.4706686
-
-This rule, as a substitute for :mod:`build_natura_raster`, downloads an already rasterized version (`natura.tiff <https://zenodo.org/record/4706686/files/natura.tiff>`_) of `Natura 2000 <https://en.wikipedia.org/wiki/Natura_2000>`_ natural protection areas to reduce computation times. The file is placed into the ``resources`` sub-directory.
-
-**Relevant Settings**
-
-.. code:: yaml
-
-    enable:
-        build_natura_raster:
-
-.. seealso::
-    Documentation of the configuration file ``config.yaml`` at
-    :ref:`toplevel_cf`
-
-**Outputs**
-
-- ``resources/natura.tiff``: Rasterized version of `Natura 2000 <https://en.wikipedia.org/wiki/Natura_2000>`_ natural protection areas to reduce computation times.
-
-.. seealso::
-    For details see :mod:`build_natura_raster`.

doc/release_notes.rst

@@ -62,6 +62,8 @@ Upcoming Release
 
 * New network topology extracted from the ENTSO-E interactive map.
 
+* Remove rules to build or retrieve rasterized NATURA 2000 dataset. Renewable potential calculation now directly uses the shapefiles.
+
 
 PyPSA-Eur 0.4.0 (22th September 2021)
 =====================================

doc/tutorial.rst

@@ -35,8 +35,8 @@ To run the tutorial, use this as your configuration file ``config.yaml``.
 
     .../pypsa-eur % cp config.tutorial.yaml config.yaml
 
-This configuration is set to download a reduced data set via the rules :mod:`retrieve_databundle`,
+This configuration is set to download a reduced data set via the rules :mod:`retrieve_databundle`
-:mod:`retrieve_natura_raster`, :mod:`retrieve_cutout` totalling at less than 250 MB.
+and :mod:`retrieve_cutout` totalling at less than 250 MB.
 The full set of data dependencies would consume 5.3 GB.
 For more information on the data dependencies of PyPSA-Eur, continue reading :ref:`data`.
 

scripts/build_natura_raster.py

@@ -1,92 +0,0 @@
-# SPDX-FileCopyrightText: : 2017-2020 The PyPSA-Eur Authors
-#
-# SPDX-License-Identifier: MIT
-
-"""
-Rasters the vector data of the `Natura 2000 <https://en.wikipedia.org/wiki/Natura_2000>`_ natural protection areas onto all cutout regions.
-
-Relevant Settings
------------------
-
-.. code:: yaml
-
-    renewable:
-        {technology}:
-            cutout:
-
-.. seealso::
-    Documentation of the configuration file ``config.yaml`` at
-    :ref:`renewable_cf`
-
-Inputs
-------
-
-- ``data/bundle/natura/Natura2000_end2015.shp``: `Natura 2000 <https://en.wikipedia.org/wiki/Natura_2000>`_ natural protection areas.
-
-    .. image:: ../img/natura.png
-        :scale: 33 %
-
-Outputs
--------
-
-- ``resources/natura.tiff``: Rasterized version of `Natura 2000 <https://en.wikipedia.org/wiki/Natura_2000>`_ natural protection areas to reduce computation times.
-
-    .. image:: ../img/natura.png
-        :scale: 33 %
-
-Description
------------
-
-"""
-
-import logging
-from _helpers import configure_logging, retrieve_snakemake_keys
-
-import atlite
-import geopandas as gpd
-import rasterio as rio
-from rasterio.features import geometry_mask
-from rasterio.warp import transform_bounds
-
-logger = logging.getLogger(__name__)
-
-
-def determine_cutout_xXyY(cutout_name):
-    cutout = atlite.Cutout(cutout_name)
-    assert cutout.crs.to_epsg() == 4326
-    x, X, y, Y = cutout.extent
-    dx, dy = cutout.dx, cutout.dy
-    return [x - dx/2., X + dx/2., y - dy/2., Y + dy/2.]
-
-
-def get_transform_and_shape(bounds, res):
-    left, bottom = [(b // res)* res for b in bounds[:2]]
-    right, top = [(b // res + 1) * res for b in bounds[2:]]
-    shape = int((top - bottom) // res), int((right - left) / res)
-    transform = rio.Affine(res, 0, left, 0, -res, top)
-    return transform, shape
-
-
-if __name__ == "__main__":
-    if 'snakemake' not in globals():
-        from _helpers import mock_snakemake
-        snakemake = mock_snakemake('build_natura_raster')
-    configure_logging(snakemake)
-
-    paths, config, wildcards, logs, out = retrieve_snakemake_keys(snakemake)
-
-    cutouts = paths.cutouts
-    xs, Xs, ys, Ys = zip(*(determine_cutout_xXyY(cutout) for cutout in cutouts))
-    bounds = transform_bounds(4326, 3035, min(xs), min(ys), max(Xs), max(Ys))
-    transform, out_shape = get_transform_and_shape(bounds, res=100)
-
-    # adjusted boundaries
-    shapes = gpd.read_file(paths.natura).to_crs(3035)
-    raster = ~geometry_mask(shapes.geometry, out_shape[::-1], transform)
-    raster = raster.astype(rio.uint8)
-
-    with rio.open(out[0], 'w', driver='GTiff', dtype=rio.uint8,
-                  count=1, transform=transform, crs=3035, compress='lzw',
-                  width=raster.shape[1], height=raster.shape[0]) as dst:
-        dst.write(raster, indexes=1)
-

scripts/build_renewable_profiles.py

@@ -224,7 +224,7 @@ if __name__ == '__main__':
     excluder = atlite.ExclusionContainer(crs=3035, res=100)
 
     if config['natura']:
-        excluder.add_raster(snakemake.input.natura, nodata=0, allow_no_overlap=True)
+        excluder.add_geometry(snakemake.input.natura)
 
     corine = config.get("corine", {})
     if "grid_codes" in corine:

test/config.test1.yaml

@@ -32,8 +32,6 @@ enable:
   retrieve_databundle: true
   build_cutout: false
   retrieve_cutout: true
-  build_natura_raster: false
-  retrieve_natura_raster: true
   custom_busmap: false
 
 electricity: