Update test config; increase figure resolution

Author: sjpfenninger

tests/integration/test_config.yaml

@@ -1,10 +1,31 @@
 module_area_potentials:
-  # Options for buffering: either a string of the form "epsg:xxxx" or "UTM"/"utm"
-  # - "UTM": project each shape to the UTM zone of its centroid for buffering
-  # - "epsg:xxxx": use the specified CRS for all buffering
-  #   A good option is "epsg:8857" (WGS 84 / Equal Earth Greenwich) for global coverage
   buffer_crs: "epsg:8857"
 
+  land_cover_types:
+    POST_FLOODING: FARM
+    RAINFED_CROPLANDS: FARM
+    MOSAIC_CROPLAND: FARM
+    MOSAIC_VEGETATION: FARM
+    CLOSED_TO_OPEN_BROADLEAVED_FOREST: FOREST
+    CLOSED_BROADLEAVED_FOREST: FOREST
+    OPEN_BROADLEAVED_FOREST: FOREST
+    CLOSED_NEEDLELEAVED_FOREST: FOREST
+    OPEN_NEEDLELEAVED_FOREST: FOREST
+    CLOSED_TO_OPEN_MIXED_FOREST: FOREST
+    MOSAIC_FOREST: FOREST
+    CLOSED_TO_OPEN_REGULARLY_FLOODED_FOREST: FOREST
+    CLOSED_REGULARLY_FLOODED_FOREST: FOREST
+    MOSAIC_GRASSLAND: OTHER
+    CLOSED_TO_OPEN_SHRUBLAND: OTHER
+    CLOSED_TO_OPEN_HERBS: OTHER
+    SPARSE_VEGETATION: OTHER
+    CLOSED_TO_OPEN_REGULARLY_FLOODED_GRASSLAND: OTHER
+    BARE_AREAS: OTHER
+    ARTIFICIAL_SURFACES_AND_URBAN_AREAS: URBAN
+    WATER_BODIES: WATER
+    PERMANENT_SNOW: NOT_SUITABLE
+    NO_DATA: NOT_SUITABLE
+
   techs:
     pv_rooftop:
       initial_area: settlement_area

workflow/scripts/resample.py

@@ -107,9 +107,9 @@ def determine_pixel_areas(raster_input):
     """
     # the following is based on https://gis.stackexchange.com/a/288034/77760
     # and assumes the data to be in EPSG:4326
-    assert raster_input.rio.crs.to_epsg() == 4326, (
-        "raster_input does not have the projection EPSG:4326"
-    )
+    assert (
+        raster_input.rio.crs.to_epsg() == 4326
+    ), "raster_input does not have the projection EPSG:4326"
     resolution = raster_input.rio.resolution()[0]  # resolution in degrees
     varea_of_pixel = np.vectorize(lambda lat: _area_of_pixel(resolution, lat))
     pixel_area = varea_of_pixel(raster_input.y) * 1000**2  # convert to m^2
@@ -289,7 +289,7 @@ def resample_inputs(
 
     print("Saving image to plot path:", plot_path)
     # If needed, resample `resampled` to fit within a maximum of `max_pixels` pixels
-    max_pixels = 1000000
+    max_pixels = 5000000
     total_pixels = resampled.sizes["y"] * resampled.sizes["x"]
     if total_pixels > max_pixels:
         # Calculate the new resolution to fit within the max_pixels limit

workflow/scripts/script_utils.py

@@ -29,6 +29,6 @@ def plot_all_dataset_variables(ds, ncols=2, savefig=None):
     plt.tight_layout()
 
     if savefig:
-        plt.savefig(savefig, bbox_inches="tight")
+        plt.savefig(savefig, dpi=300, bbox_inches="tight")
 
     return fig