leveecmd 2.8.0__py3-none-any.whl

leveecmd/__init__.py

@@ -0,0 +1 @@
+""" Beschrijving van leveecmd """

leveecmd/__main__.py

@@ -0,0 +1,3 @@
+if __name__ == "__main__":  # pragma: no cover
+    from leveecmd.app.main import app
+    app()

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/buffer_ref_line.py

@@ -0,0 +1,39 @@
+from __future__ import annotations
+import geopandas as gpd
+from shapely.geometry import Polygon, LineString
+import os
+from typing import List, TYPE_CHECKING
+import warnings
+
+from leveecmd.cmd_components.live_controller import LiveController, DummyLiveController
+from pathlib import Path
+if TYPE_CHECKING:
+    from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.tool import InputAHNLangsKering
+
+
+def buffer_ref_line(
+        path_to_ref_line: str, buffer_distance: float, dir_to_ref_line_buffer: str, path_to_ref_line_buffer: str):
+    # Read ref file (and anticipate measure warning)
+    with warnings.catch_warnings():
+        warnings.filterwarnings("ignore", message="Measured \\(M\\) geometry types are not supported.*")
+        lines = gpd.read_file(path_to_ref_line)
+
+    # Create buffer
+    polygons: List[Polygon] = []
+    for line in lines.geometry:
+        line: LineString
+        polygons.append(line.buffer(distance=buffer_distance))
+    polygon_gdf = gpd.GeoDataFrame(geometry=polygons, crs=lines.crs)
+    os.makedirs(dir_to_ref_line_buffer, exist_ok=True)
+    polygon_gdf.to_file(Path(path_to_ref_line_buffer))
+
+
+def step_buffer_ref_line(
+        tool_input: InputAHNLangsKering, live_controller: LiveController = DummyLiveController(), **_
+):  # pragma: no cover
+    buffer_ref_line(
+        path_to_ref_line=tool_input.path_to_ref_line,
+        buffer_distance=tool_input.buffer_distance,
+        dir_to_ref_line_buffer=tool_input.dir_to_ref_line_buffer,
+        path_to_ref_line_buffer=tool_input.path_to_ref_line_buffer)
+    live_controller.print(f"Finished buffering ref line")

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/clip_rasters.py

@@ -0,0 +1,65 @@
+from __future__ import annotations
+from typing import TYPE_CHECKING
+import geopandas as gpd
+from rasterio.mask import mask
+from shapely.geometry import mapping
+import os
+from shapely.geometry import box
+import os.path
+import rasterio
+from leveecmd.cmd_components.live_controller import DummyLiveController, LiveController
+if TYPE_CHECKING:
+    from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.tool import InputAHNLangsKering
+
+
+def clip_rasters(path_to_ref_line_buffer: str, export_dir: str, required_raster_names: list[str], ahn_tiles_dir: str):
+
+    # Load the polygon shapefile
+    polygon = gpd.read_file(path_to_ref_line_buffer)
+    geoms = [mapping(geom) for geom in polygon.geometry]
+
+    # Create clip folder
+    clips_dir = os.path.join(export_dir, "clips")
+    os.makedirs(clips_dir, exist_ok=True)
+
+    # Create clips
+    for raster_file_name in required_raster_names:
+        print(f"Now clipping raster: {raster_file_name}.")
+
+        file_path = os.path.join(ahn_tiles_dir, raster_file_name)
+        if not os.path.isfile(file_path):  # pragma: no cover
+            continue
+        if not raster_file_name.endswith(".tif"):  # pragma: no cover
+            continue
+
+        # Open the raster file
+        with rasterio.open(file_path) as src:
+
+            raster_bounds = src.bounds
+            raster_box = box(*raster_bounds)
+            if not polygon.intersects(raster_box).any():  # pragma: no cover
+                print(f"{raster_file_name} does not overlap polygon and is skipped")
+                continue
+
+            out_image, out_transform = mask(src, geoms, crop=True)
+            out_meta = src.meta.copy()
+
+        # Update metadata
+        out_meta.update({
+            "driver": "GTiff",
+            "height": out_image.shape[1],
+            "width": out_image.shape[2],
+            "transform": out_transform
+        })
+
+        # Save the clipped raster
+        with rasterio.open(os.path.join(clips_dir, f"clip_{raster_file_name}"), "w", **out_meta) as dest:
+            dest.write(out_image)
+        print(f" -> Finished clipping raster {raster_file_name}.")
+
+
+def step_clip_rasters(tool_input: InputAHNLangsKering, live_controller: LiveController = DummyLiveController(), **_):  # pragma: no cover
+    clip_rasters(
+        path_to_ref_line_buffer=tool_input.path_to_ref_line_buffer, export_dir=tool_input.export_dir,
+        required_raster_names=tool_input.required_raster_names, ahn_tiles_dir=tool_input.ahn_tiles_dir)
+    live_controller.print(f" -> Finished clipping all rasters.")

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/collect_raster_names_to_download.py

@@ -0,0 +1,34 @@
+from __future__ import annotations
+from typing import List, TYPE_CHECKING
+import os
+from leveecmd.cmd_components.live_controller import DummyLiveController, LiveController
+if TYPE_CHECKING:
+    from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.tool import InputAHNLangsKering
+
+
+def collect_raster_names_to_download(required_raster_names: list[str], ahn_tiles_dir: str):
+    # TODO Later Should Middel: Wanneer AHN6 komt, zal deze functie foutief denken dat die al in de storage-folder
+    #  staat. In werkelijkheid staat daar AHN5. Maar er is geen manier om het versienummer to controleren op dit moment.
+
+    rasters_to_download: List[str] = []
+    for raster_file_name in required_raster_names:
+        path_to_ahn_tile = os.path.join(ahn_tiles_dir, raster_file_name)
+        if not os.path.exists(path_to_ahn_tile):
+            rasters_to_download.append(raster_file_name)
+            continue
+    return rasters_to_download
+
+
+def step_collect_raster_names_to_download(
+        tool_input: InputAHNLangsKering, live_controller: LiveController = DummyLiveController(), **_):  # pragma: no cover
+    tool_input.rasters_to_download = collect_raster_names_to_download(
+        required_raster_names=tool_input.required_raster_names, ahn_tiles_dir=tool_input.ahn_tiles_dir)
+
+    # Feedback in terminal
+    if tool_input.rasters_to_download.__len__() == 0:
+        live_controller.print(
+            msg=f"Rasters al aanwezig in storage directory ({tool_input.ahn_tiles_dir}).", pinned=True)
+    else:
+        live_controller.print(
+            msg=f"Volgende rasters worden gedownload: {', '.join(tool_input.rasters_to_download)}.\n"
+                f"Dit gebeurt naar de directory {tool_input.ahn_tiles_dir}.", pinned=True)

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/collect_required_raster_names.py

@@ -0,0 +1,28 @@
+from __future__ import annotations
+import geopandas as gpd
+from typing import TYPE_CHECKING
+from leveecmd.cmd_components.live_controller import DummyLiveController, LiveController
+from importlib.resources import files
+if TYPE_CHECKING:
+    from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.tool import InputAHNLangsKering
+
+
+PATH_TO_KAARTBLAD_INDEX_SHP = (
+    files('leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.files.dependencies.kaart_blad_index')
+    .joinpath('kaart_blad_index.shp'))
+PARTIAL_URL_TO_PDOK = "https://service.pdok.nl/rws/ahn/atom/downloads/dtm_05m/"
+
+
+def collect_required_raster_names(path_to_ref_line_buffer: str):
+    gdf_kaartblad_index = gpd.read_file(PATH_TO_KAARTBLAD_INDEX_SHP)
+    gdf_ref_line_buffer = gpd.read_file(path_to_ref_line_buffer)
+    pg_ref_line_buffer = gdf_ref_line_buffer.geometry.iloc[0]
+    gdf_overlapping: gpd.GeoDataFrame = gdf_kaartblad_index[gdf_kaartblad_index.geometry.intersects(pg_ref_line_buffer)]
+    return gdf_overlapping['file_name'].values.tolist()
+
+
+def step_collect_required_raster_names(
+        tool_input: InputAHNLangsKering, live_controller: LiveController = DummyLiveController(), **_):  # pragma: no cover
+    tool_input.required_raster_names = collect_required_raster_names(
+        path_to_ref_line_buffer=tool_input.path_to_ref_line_buffer)
+    live_controller.print(msg=f"Benodigde rasters: {', '.join(tool_input.required_raster_names)}", pinned=True)

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/download_rasters.py

@@ -0,0 +1,60 @@
+from __future__ import annotations
+from typing import TYPE_CHECKING
+import requests
+import os
+import sys
+from leveecmd.cmd_components.live_controller import DummyLiveController, LiveController
+from leveecmd.utils.other import BColors
+if TYPE_CHECKING:
+    from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.tool import InputAHNLangsKering
+
+
+def download_raster(raster_file_name: str, ahn_tiles_dir: str):
+    url = f"https://service.pdok.nl/rws/ahn/atom/downloads/dtm_05m/{raster_file_name}"
+
+    # Download raster
+    export_path = os.path.join(ahn_tiles_dir, raster_file_name)
+    try:
+        response = requests.get(url)
+    except requests.exceptions.RequestException:  # pragma: no cover
+        print(f"{BColors.WARNING}Unable to connect to PDOK. Are you sure you are connected to the internet?\n"
+              f"Tried connection url was {url}.{BColors.ENDC}")
+        sys.exit()
+    with open(export_path, "wb") as f:
+        f.write(response.content)
+
+
+def step_download_raster(
+        tool_input: InputAHNLangsKering,
+        raster_file_name: str,
+        live_controller: LiveController = DummyLiveController(),
+        **_):  # pragma: no cover
+    live_controller.print(msg=f"Now downloading {raster_file_name}.")
+    download_raster(raster_file_name=raster_file_name, ahn_tiles_dir=tool_input.ahn_tiles_dir)
+    live_controller.print(msg=f" -> Downloaded raster {raster_file_name}.")
+
+
+def download_rasters(rasters_to_download: list[str], ahn_tiles_dir: str):
+
+    if rasters_to_download.__len__() == 0:  # pragma: no cover
+        return
+
+    print(f"Now downloading the following rasters:")
+    print(", ".join(rasters_to_download))
+    print(f"Note: It may take several minutes per raster.")
+
+    for raster_file_name in rasters_to_download:
+        download_raster(raster_file_name=raster_file_name, ahn_tiles_dir=ahn_tiles_dir)
+
+
+def step_download_rasters(tool_input: InputAHNLangsKering, live_controller: LiveController = DummyLiveController(), **_):  # pragma: no cover
+
+    if tool_input.rasters_to_download.__len__() == 0:
+        return
+
+    live_controller.print(f"Now downloading the following rasters:")
+    live_controller.print(", ".join(tool_input.rasters_to_download))
+    live_controller.print(f"Note: It may take several minutes per raster.")
+
+    for raster_file_name in tool_input.rasters_to_download:
+        step_download_raster(tool_input=tool_input, raster_file_name=raster_file_name, live_controller=live_controller)

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/files/dependencies/dijktrajecten/dijktrajecten.cst

@@ -0,0 +1 @@
+ISO-8859-1

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/files/dependencies/dijktrajecten/dijktrajecten.prj

@@ -0,0 +1 @@
+PROJCS["Amersfoort / RD New", GEOGCS["Amersfoort", DATUM["Amersfoort", SPHEROID["Bessel 1841", 6377397.155, 299.1528128, AUTHORITY["EPSG","7004"]], TOWGS84[565.2369, 50.0087, 465.658, -0.40685733032239757, -0.3507326765425626, 1.8703473836067956, 4.0812], AUTHORITY["EPSG","6289"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], AXIS["Geodetic latitude", NORTH], AXIS["Geodetic longitude", EAST], AUTHORITY["EPSG","4289"]], PROJECTION["Oblique_Stereographic", AUTHORITY["EPSG","9809"]], PARAMETER["central_meridian", 5.387638888888891], PARAMETER["latitude_of_origin", 52.15616055555556], PARAMETER["scale_factor", 0.9999079], PARAMETER["false_easting", 155000.0], PARAMETER["false_northing", 463000.0], UNIT["m", 1.0], AXIS["Easting", EAST], AXIS["Northing", NORTH], AUTHORITY["EPSG","28992"]]

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/files/dependencies/kaart_blad_index/kaart_blad_index.cpg

@@ -0,0 +1 @@
+UTF-8

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/files/dependencies/kaart_blad_index/kaart_blad_index.prj

@@ -0,0 +1 @@
+PROJCS["RD_New",GEOGCS["GCS_Amersfoort",DATUM["D_Amersfoort",SPHEROID["Bessel_1841",6377397.155,299.1528128]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Double_Stereographic"],PARAMETER["False_Easting",155000.0],PARAMETER["False_Northing",463000.0],PARAMETER["Central_Meridian",5.38763888888889],PARAMETER["Scale_Factor",0.9999079],PARAMETER["Latitude_Of_Origin",52.1561605555556],UNIT["Meter",1.0]]

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/files/examples/ref_line/mline224.cpg

@@ -0,0 +1 @@
+UTF-8

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/files/examples/ref_line/mline224.prj

@@ -0,0 +1 @@
+PROJCS["RD_New",GEOGCS["GCS_Amersfoort",DATUM["D_Amersfoort",SPHEROID["Bessel_1841",6377397.155,299.1528128]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Double_Stereographic"],PARAMETER["False_Easting",155000.0],PARAMETER["False_Northing",463000.0],PARAMETER["Central_Meridian",5.38763888888889],PARAMETER["Scale_Factor",0.9999079],PARAMETER["Latitude_Of_Origin",52.1561605555556],UNIT["Meter",1.0]]

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/input_form.py

@@ -0,0 +1,197 @@
+from __future__ import annotations
+from InquirerPy import inquirer
+from leveecmd.cmd_components.live_controller import LiveController
+from leveecmd.utils.other import BColors
+import os
+import shutil
+from typing import Optional
+from typing import TYPE_CHECKING
+from geopandas import GeoDataFrame, read_file
+from shapely.geometry import LineString, MultiLineString
+from importlib.resources import files
+from importlib.resources.readers import MultiplexedPath
+from pathlib import Path
+import importlib.resources
+from leveecmd.utils.other import clear_terminal
+if TYPE_CHECKING:
+    from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.tool import InputAHNLangsKering
+
+
+def input_ref_line(
+        live_controller: LiveController,
+) -> Optional[str]:
+
+    shp_file_path: Optional[str] = None
+    shp_file_is_valid = False
+    live_controller.pause()
+    while shp_file_is_valid is False:
+        shp_file_path: str = inquirer.text(
+            message="Specificeer het volledige bestandspad naar de referentielijn (shapefile). ",
+        ).execute()
+
+        shp_file_path = shp_file_path.replace('"', '')
+
+        if not shp_file_path.endswith(".shp"):
+            print(BColors.WARNING, f"Het bestand met een shapefile zijn, eindigend op .shp. Jouw invoer eindigt op "
+                                   f"extensie .{shp_file_path.split(sep='.')[-1]}.", BColors.ENDC)
+            continue
+
+        shp_file_is_valid = True
+        live_controller.print(msg=f"Bestandspad naar opgegeven referentielijn is '{shp_file_path}'.", pinned=True)
+    live_controller.resume()
+
+    return shp_file_path
+
+
+def input_normtraject_number(
+        tool_input: InputAHNLangsKering,
+        live_controller: LiveController,
+) -> Optional[str]:
+
+    shp_file_path: Optional[str] = None
+    traject_number_is_valid = False
+    live_controller.pause()
+    while traject_number_is_valid is False:
+        traject_number: str = inquirer.text(
+            message="Specificeer het nummer van het normtraject. ",
+        ).execute()
+
+        # Gather traject numbers
+        with importlib.resources.path(
+                package='leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.files.dependencies.dijktrajecten',
+                resource='dijktrajecten.shp') as shp_path:
+            gdf: GeoDataFrame = read_file(shp_path)
+        possible_traject_numbers = gdf['TRAJECT_ID'].values.tolist()
+        possible_traject_numbers.sort()
+
+        # Check if number exists
+        if traject_number not in possible_traject_numbers:
+            print(BColors.WARNING,
+                  f"Het opgegeven trajectnummer '{traject_number}' is niet beschikbaar. De mogelijke opties zijn "
+                  f"{', '.join(possible_traject_numbers)}. ", BColors.ENDC)
+            continue
+
+        # Export as shape (for user to verify, plus in line with other input options)
+        geom: LineString | MultiLineString = gdf[gdf['TRAJECT_ID'] == traject_number]['geometry'].iloc[0]
+        if isinstance(geom, LineString):
+            ls: LineString = geom
+        else:
+            raise NotImplementedError(f"Type {type(geom)} is not yet supported.")
+        gdf = GeoDataFrame(geometry=[ls], crs="EPSG:28992")
+        file_name = f"mline{traject_number}.shp"
+        export_dir_shp = os.path.join(tool_input.export_dir, "ref_line")
+        os.makedirs(export_dir_shp, exist_ok=True)
+        shp_file_path = os.path.join(export_dir_shp, file_name)
+        gdf.to_file(Path(os.path.join(export_dir_shp, file_name)))
+
+        traject_number_is_valid = True
+        live_controller.print(msg=f"Opgegeven trajectnummer is '{traject_number}'.", pinned=True)
+        live_controller.print(msg=f"Kopie van referentielijn opgeslagen onder '{shp_file_path}'.", pinned=True)
+
+    live_controller.resume()
+
+    return shp_file_path
+
+
+def input_export_dir(
+        live_controller: LiveController,
+):
+
+    export_dir_path: Optional[str] = None
+    export_dir_is_valid = False
+    live_controller.pause()
+    while export_dir_is_valid is False:
+        export_dir_path: str = inquirer.text(
+            message="Specificeer het volledige bestandspad naar de exporteer locatie.",
+        ).execute()
+
+        export_dir_path = export_dir_path.replace('"', '')
+
+        if not os.path.exists(export_dir_path):
+            print(BColors.WARNING, f"De opgegeven map bestaat niet.", BColors.ENDC)
+            continue
+        if not os.path.isdir(export_dir_path):
+            print(BColors.WARNING, f"De opgegeven locatie is geen map.", BColors.ENDC)
+            continue
+        if os.listdir(export_dir_path):
+            print(BColors.WARNING, f"De opgegeven map is niet leeg. ", BColors.ENDC)
+            continue
+
+        export_dir_is_valid = True
+        live_controller.print(msg=f"Opgegeven export directory is '{export_dir_path}'.", pinned=True)
+
+    live_controller.resume()
+
+    return export_dir_path
+
+
+def get_path_from_multiplexed_path(multiplexed: MultiplexedPath) -> Path:
+    # noinspection PyProtectedMember
+    # noinspection PyUnresolvedReferences
+    paths: list[Path] = multiplexed._paths
+    assert paths.__len__() == 1
+    return paths[0]
+
+
+def input_ref_line_example_or_custom(
+        tool_input: InputAHNLangsKering,
+        live_controller: LiveController,
+):
+
+    live_controller.pause()
+    choice = inquirer.select(
+        message="Wil je een voorbeeld referentielijn gebruiken? De mline van een normtraject? "
+                "Of een eigen shapefile met referentielijn?",
+        choices=["Voorbeeld referentielijn", "Normtraject", "Eigen referentielijn"],
+        default="Voorbeeld referentielijn",
+    ).execute()
+    live_controller.resume()
+
+    live_controller.print(msg=f"Type gekozen referentielijn is '{choice}'.", pinned=True)
+    if choice == "Voorbeeld referentielijn":
+        shp_dir = files("leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.files.examples.ref_line")
+        assert isinstance(shp_dir, MultiplexedPath)
+        shp_dir = get_path_from_multiplexed_path(shp_dir)
+        assert isinstance(shp_dir, Path), f"Expected type Path from shp_dir, got type {type(shp_dir)}."
+        new_dir = os.path.join(tool_input.export_dir, "ref_line")
+        shutil.copytree(src=shp_dir, dst=new_dir)
+        tool_input.path_to_ref_line = os.path.join(new_dir, "mline224.shp")
+    elif choice == "Normtraject":
+        tool_input.path_to_ref_line = input_normtraject_number(tool_input=tool_input, live_controller=live_controller)
+    elif choice == "Eigen referentielijn":
+        tool_input.path_to_ref_line = input_ref_line(live_controller=live_controller)
+    else:
+        raise ValueError()
+
+
+def input_buffer_distance(
+        tool_input: InputAHNLangsKering,
+        live_controller: LiveController,
+):
+
+    buffer_distance_is_valid = False
+    live_controller.pause()
+    while buffer_distance_is_valid is False:
+        buffer_afstand: str = inquirer.text(
+            message="Specificeer de buffer afstand vanaf referentie lijn in hele meters. ",
+        ).execute()
+
+        try:
+            tool_input.buffer_distance = int(buffer_afstand)
+            buffer_distance_is_valid = True
+            live_controller.print(msg=f"Gekozen buffer distance is '{tool_input.buffer_distance} meter'.", pinned=True)
+        except ValueError:
+            print(BColors.WARNING, f"De invoer '{buffer_afstand}' is geen geldig geheel getal.", BColors.ENDC)
+
+    live_controller.resume()
+
+
+def input_form(
+        tool_input: InputAHNLangsKering,
+        live_controller: LiveController,
+        **_
+):
+    clear_terminal()
+    tool_input.export_dir = input_export_dir(live_controller=live_controller)
+    input_ref_line_example_or_custom(tool_input=tool_input, live_controller=live_controller)
+    input_buffer_distance(tool_input=tool_input, live_controller=live_controller)

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/merge_rasters.py

@@ -0,0 +1,56 @@
+from __future__ import annotations
+from typing import TYPE_CHECKING
+import os
+import os.path
+import rasterio
+from rasterio.merge import merge
+from pathlib import Path
+import numpy as np
+from leveecmd.cmd_components.live_controller import DummyLiveController, LiveController
+import shutil
+if TYPE_CHECKING:
+    from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.tool import InputAHNLangsKering
+
+
+def merge_clips(export_dir: str):
+
+    # Folder containing the clipped rasters
+    clips_dir = Path(os.path.join(export_dir, "clips"))
+
+    # List all TIF files
+    clip_files = list(clips_dir.glob("*.tif"))
+
+    # Open all raster files
+    src_files_to_mosaic = [rasterio.open(fp) for fp in clip_files]
+
+    # Merge the rasters
+    mosaic, out_transform = merge(src_files_to_mosaic, nodata=np.nan)
+
+    # Copy metadata from one of the input files
+    out_meta = src_files_to_mosaic[0].meta.copy()
+    out_meta.update({
+        "height": mosaic.shape[1],
+        "width": mosaic.shape[2],
+        "transform": out_transform,
+        "nodata": np.nan
+    })
+
+    # Save the merged raster
+    dir_merged = Path(os.path.join(export_dir, "merged"))
+    os.makedirs(dir_merged, exist_ok=True)
+    merged_file_path = os.path.join(dir_merged, "merged_output.tif")
+    with rasterio.open(merged_file_path, "w", **out_meta) as dest:
+        dest.write(mosaic)
+
+    # Remove redundant files (now the app is done)
+    for raster in src_files_to_mosaic:
+        raster.close()
+    shutil.rmtree(clips_dir)
+
+    return merged_file_path
+
+
+def step_merge_clips(tool_input: InputAHNLangsKering, live_controller: LiveController = DummyLiveController(), **_):  # pragma: no cover
+    merged_file_path = merge_clips(export_dir=tool_input.export_dir)
+    live_controller.print(msg=f"Exported to {merged_file_path}", pinned=True)
+    live_controller.print(f"Merge completed")

leveecmd/app/cat_ahn/tool_download_ahn_langs_kering/test_logic.py

@@ -0,0 +1,42 @@
+from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.buffer_ref_line import buffer_ref_line
+from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.collect_required_raster_names import (
+    collect_required_raster_names)
+from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.collect_raster_names_to_download import (
+    collect_raster_names_to_download)
+from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.download_rasters import download_rasters
+from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.clip_rasters import clip_rasters
+from leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.merge_rasters import merge_clips
+from importlib.resources import files
+import os
+import pytest
+
+
+@pytest.mark.slow
+def test_ahn_langs_kering_logic(tmpdir):
+
+    # TODO: Require internet for this unit test, otherwise end here. 
+
+    # Buffer
+    dir_to_ref_line_buffer = os.path.join(tmpdir, "ref_line_buffer")
+    path_to_ref_line_buffer = os.path.join(dir_to_ref_line_buffer, "ref_line_buffer.shp")
+    path_to_ref_line = str(files(
+        "leveecmd.app.cat_ahn.tool_download_ahn_langs_kering.files.examples.ref_line").joinpath("mline224.shp"))
+    buffer_ref_line(
+        path_to_ref_line=path_to_ref_line, buffer_distance=50, dir_to_ref_line_buffer=dir_to_ref_line_buffer,
+        path_to_ref_line_buffer=path_to_ref_line_buffer)
+
+    # Required rasters
+    required_raster_names = collect_required_raster_names(path_to_ref_line_buffer=path_to_ref_line_buffer)
+    ahn_tiles_dir = os.path.join(tmpdir, "ahn_tiles")
+    os.makedirs(ahn_tiles_dir, exist_ok=True)
+    rasters_to_download = collect_raster_names_to_download(
+        required_raster_names=required_raster_names, ahn_tiles_dir=ahn_tiles_dir)
+    download_rasters(rasters_to_download=rasters_to_download, ahn_tiles_dir=ahn_tiles_dir)
+
+    # Clip rasters
+    clip_rasters(
+        path_to_ref_line_buffer=path_to_ref_line_buffer, export_dir=tmpdir,
+        required_raster_names=required_raster_names, ahn_tiles_dir=ahn_tiles_dir)
+
+    # Merge clips
+    merge_clips(export_dir=tmpdir)