s1-denoise 0.0.0.dev0__tar.gz

s1_denoise-0.0.0.dev0/.gitignore

@@ -0,0 +1,40 @@
+# Python
+__pycache__/
+.venv/
+dist/
+
+# Editor / tool config
+.cache/
+.config/
+.vscode/
+
+# Archives & logs
+*.tar
+*.xz
+*.gz
+*.log
+*.out
+*.csv
+
+# Test data
+aux-files/
+test/batch*
+test/input*
+test/aux-files
+test/fonts
+test/*.png
+
+# Documentation
+*.pdf
+
+# Sentinel-1 data & products
+*.SAFE
+job-orders/*
+
+# Pipeline output
+output/
+compare-noise-plots/
+
+# Documentation site
+site/
+documentation*.zip

s1_denoise-0.0.0.dev0/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (C) Aresys S.r.l. <info@aresys.it>
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

s1_denoise-0.0.0.dev0/PKG-INFO

@@ -0,0 +1,147 @@
+Metadata-Version: 2.5
+Name: s1-denoise
+Version: 0.0.0.dev0
+Summary: Noise vector reconstruction for Level 1 Sentinel-1 products
+Project-URL: Homepage, https://github.com/aresys-srl/s1-denoise
+Project-URL: Documentation, https://opensource.aresys.it/s1-denoise/
+Project-URL: Repository, https://github.com/aresys-srl/s1-denoise
+Project-URL: Issues, https://github.com/aresys-srl/s1-denoise/issues
+Project-URL: Changelog, https://opensource.aresys.it/s1-denoise/API/changelog/
+Author-email: "Aresys S.R.L." <info@aresys.it>
+License: MIT License
+        
+        Copyright (C) Aresys S.r.l. <info@aresys.it>
+        
+        Permission is hereby granted, free of charge, to any person obtaining a copy
+        of this software and associated documentation files (the "Software"), to deal
+        in the Software without restriction, including without limitation the rights
+        to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+        copies of the Software, and to permit persons to whom the Software is
+        furnished to do so, subject to the following conditions:
+        
+        The above copyright notice and this permission notice shall be included in all
+        copies or substantial portions of the Software.
+        
+        THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+        IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+        FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+        AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+        LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+        OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+        SOFTWARE.
+License-File: LICENSE
+Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: Python :: 3
+Requires-Python: >=3.11
+Requires-Dist: click
+Requires-Dist: colorlog
+Requires-Dist: lxml
+Requires-Dist: matplotlib
+Requires-Dist: numpy
+Requires-Dist: perseo-core
+Requires-Dist: pillow
+Requires-Dist: rasterio
+Requires-Dist: scipy
+Requires-Dist: sentinel1decoder
+Requires-Dist: xmltodict
+Description-Content-Type: text/markdown
+
+# Sentinel-1 Level 1 Noise Processor
+
+#### Thermal noise vector reconstruction for Sentinel-1 Level 1 products
+
+The S1 L1N Processor reconstructs the thermal noise range vectors for Sentinel-1 IW and EW products (both SLC and GRD) and writes them to output XML files that follow the standard Sentinel-1 annotation format. The reconstructed vectors can then be compared against the nominal noise annotations embedded in the SAFE product to assess noise estimation quality.
+
+For each processed channel the processor:
+
+1. reads acquisition geometry and calibration data from auxiliary files (PP1, CAL, INS);
+2. models the thermal noise using range spread loss, elevation antenna pattern, and calibration LUT;
+3. optionally refines the noise power correction using an associated Level 0 Noise (L0N) product;
+4. discards noise pulses whose power deviates from the expected noise level in AUX_INS by more than −3/+1 dB;
+5. for GRD products, handles sub-swath separation, per-subswath noise reconstruction and cross-subswath merging;
+6. writes the output noise XML files and, if nominal noise annotations are present in the SAFE product, reports the Mean Absolute Error (MAE) and maximum absolute error in dB.
+
+📖 **Documentation:** <https://opensource.aresys.it/s1-denoise/>
+
+---
+
+## Installation
+
+**Requirements:** Python ≥ 3.11, pip.
+
+Install the package from PyPI into a virtual environment of your choice:
+
+```bash
+pip install s1-denoise
+```
+
+Alternatively, install the wheel attached to a [GitHub release](https://github.com/aresys-srl/s1-denoise/releases):
+
+```bash
+pip install s1_denoise-<version>-py3-none-any.whl
+```
+
+Both options install the `denoise` command into the active Python environment.
+
+---
+
+## Usage
+
+The processor is driven by a **Job Order XML file** that specifies input products, auxiliary files, and the output path.
+
+```bash
+denoise "/path/to/job/order.xml"
+```
+
+The processor reads each L1 SAFE product listed in the Job Order, generates the noise profiles, and writes the output XML files to the configured output directory. If nominal noise annotations are found inside the SAFE product, the MAE against them is logged automatically.
+
+For all available options:
+
+```bash
+denoise --help
+```
+
+To print the installed version:
+
+```bash
+denoise --version
+```
+
+### Job Order structure
+
+The Job Order is an XML file that configures:
+
+| Field | Description |
+|---|---|
+| Input L1 SAFE products | One or more Sentinel-1 SAFE folders to process |
+| Input L0N products (optional) | Associated Level 0 Noise products for refined power correction |
+| Auxiliary files | AUX_PP1, AUX_CAL, AUX_INS product paths |
+| Output path | Directory where output XML files are written |
+| Log level | Verbosity of the processor log |
+
+A Job Order template is provided in [`templates/JobOrder.xml`](templates/JobOrder.xml). See the [Job Order documentation](https://opensource.aresys.it/s1-denoise/user-guide/job-order/) for the supported file types.
+
+---
+
+## Visualisation
+
+To visually compare the original noise annotation with the reconstructed one:
+
+```bash
+./scripts/compare_noise.py <noise_xml_filepath> <denoise_xml_filepath> --colorbar 0.01
+```
+
+The `--colorbar` argument sets the symmetric boundary (in dB) for the 2D difference plot colorbar.
+
+---
+
+## Documentation
+
+The documentation is built with [Zensical](https://zensical.org/). To build it locally:
+
+```bash
+pip install zensical mkdocstrings-python
+python tools/build_doc.py
+```
+
+The site is generated in the `site/` folder. Use `zensical serve` for a live preview while editing the pages in `docs/`.

s1_denoise-0.0.0.dev0/README.md

@@ -0,0 +1,99 @@
+# Sentinel-1 Level 1 Noise Processor
+
+#### Thermal noise vector reconstruction for Sentinel-1 Level 1 products
+
+The S1 L1N Processor reconstructs the thermal noise range vectors for Sentinel-1 IW and EW products (both SLC and GRD) and writes them to output XML files that follow the standard Sentinel-1 annotation format. The reconstructed vectors can then be compared against the nominal noise annotations embedded in the SAFE product to assess noise estimation quality.
+
+For each processed channel the processor:
+
+1. reads acquisition geometry and calibration data from auxiliary files (PP1, CAL, INS);
+2. models the thermal noise using range spread loss, elevation antenna pattern, and calibration LUT;
+3. optionally refines the noise power correction using an associated Level 0 Noise (L0N) product;
+4. discards noise pulses whose power deviates from the expected noise level in AUX_INS by more than −3/+1 dB;
+5. for GRD products, handles sub-swath separation, per-subswath noise reconstruction and cross-subswath merging;
+6. writes the output noise XML files and, if nominal noise annotations are present in the SAFE product, reports the Mean Absolute Error (MAE) and maximum absolute error in dB.
+
+📖 **Documentation:** <https://opensource.aresys.it/s1-denoise/>
+
+---
+
+## Installation
+
+**Requirements:** Python ≥ 3.11, pip.
+
+Install the package from PyPI into a virtual environment of your choice:
+
+```bash
+pip install s1-denoise
+```
+
+Alternatively, install the wheel attached to a [GitHub release](https://github.com/aresys-srl/s1-denoise/releases):
+
+```bash
+pip install s1_denoise-<version>-py3-none-any.whl
+```
+
+Both options install the `denoise` command into the active Python environment.
+
+---
+
+## Usage
+
+The processor is driven by a **Job Order XML file** that specifies input products, auxiliary files, and the output path.
+
+```bash
+denoise "/path/to/job/order.xml"
+```
+
+The processor reads each L1 SAFE product listed in the Job Order, generates the noise profiles, and writes the output XML files to the configured output directory. If nominal noise annotations are found inside the SAFE product, the MAE against them is logged automatically.
+
+For all available options:
+
+```bash
+denoise --help
+```
+
+To print the installed version:
+
+```bash
+denoise --version
+```
+
+### Job Order structure
+
+The Job Order is an XML file that configures:
+
+| Field | Description |
+|---|---|
+| Input L1 SAFE products | One or more Sentinel-1 SAFE folders to process |
+| Input L0N products (optional) | Associated Level 0 Noise products for refined power correction |
+| Auxiliary files | AUX_PP1, AUX_CAL, AUX_INS product paths |
+| Output path | Directory where output XML files are written |
+| Log level | Verbosity of the processor log |
+
+A Job Order template is provided in [`templates/JobOrder.xml`](templates/JobOrder.xml). See the [Job Order documentation](https://opensource.aresys.it/s1-denoise/user-guide/job-order/) for the supported file types.
+
+---
+
+## Visualisation
+
+To visually compare the original noise annotation with the reconstructed one:
+
+```bash
+./scripts/compare_noise.py <noise_xml_filepath> <denoise_xml_filepath> --colorbar 0.01
+```
+
+The `--colorbar` argument sets the symmetric boundary (in dB) for the 2D difference plot colorbar.
+
+---
+
+## Documentation
+
+The documentation is built with [Zensical](https://zensical.org/). To build it locally:
+
+```bash
+pip install zensical mkdocstrings-python
+python tools/build_doc.py
+```
+
+The site is generated in the `site/` folder. Use `zensical serve` for a live preview while editing the pages in `docs/`.

s1_denoise-0.0.0.dev0/pyproject.toml

@@ -0,0 +1,56 @@
+[project]
+name = "s1-denoise"
+description = "Noise vector reconstruction for Level 1 Sentinel-1 products"
+authors = [
+  { name="Aresys S.R.L.", email="info@aresys.it" }
+]
+readme = "README.md"
+requires-python = ">=3.11"
+classifiers = [
+    "Programming Language :: Python :: 3",
+    "Operating System :: OS Independent",
+]
+dynamic = ["version"]
+license = { file = "LICENSE" }
+dependencies = [
+    "click",
+    "colorlog",
+    "lxml",
+    "matplotlib",
+    "numpy",
+    "Pillow",
+    "rasterio",
+    "scipy",
+    "sentinel1decoder",
+    "xmltodict",
+    "perseo-core"
+]
+
+[project.urls]
+Homepage = "https://github.com/aresys-srl/s1-denoise"
+Documentation = "https://opensource.aresys.it/s1-denoise/"
+Repository = "https://github.com/aresys-srl/s1-denoise"
+Issues = "https://github.com/aresys-srl/s1-denoise/issues"
+Changelog = "https://opensource.aresys.it/s1-denoise/API/changelog/"
+
+[project.scripts]
+denoise = "s1_denoise.__main__:main"
+
+[build-system]
+requires = ["hatchling"]
+build-backend = "hatchling.build"
+
+[tool.hatch.version]
+path = "src/s1_denoise/__init__.py"
+
+[tool.hatch.build.targets.wheel]
+packages = ["src/s1_denoise"]
+
+[tool.hatch.build.targets.sdist]
+include = [
+    "src/s1_denoise/**/*.py",
+    "README.md",
+    "LICENSE",
+    "pyproject.toml",
+    "scripts/compare_noise.py"
+]

s1_denoise-0.0.0.dev0/scripts/compare_noise.py

@@ -0,0 +1,230 @@
+#!.venv/bin/python3
+
+"""
+Compare a nominal and a reconstructed Sentinel-1 noise XML file and plot their differences.
+
+Prints the mean and maximum absolute errors in dB and saves to ``compare-noise-plots/`` a 2D map of the difference
+and range profiles at a few azimuth times.
+"""
+
+import numpy as np
+import argparse
+import sys
+import matplotlib.pyplot as plt
+from pathlib import Path
+from math import nan
+
+from s1_denoise.error_calculation import compare_noise, DN_VALUES
+from s1_denoise.utils import linear_to_db
+
+OUTPUT_DIR = Path("compare-noise-plots")
+
+
+def plot_2d(common_pixels, common_times_sec, diff_db, out_name, colorbar, pol, subswath, prod_type):
+    """
+    Plot the 2D map of the noise difference and save it as ``whole_raster.png``.
+
+    Parameters
+    ----------
+    common_pixels : numpy.ndarray
+        Range pixel axis.
+    common_times_sec : numpy.ndarray
+        Azimuth time axis [s] from acquisition start.
+    diff_db : numpy.ndarray
+        Reconstructed minus nominal noise [dB], shape ``(len(common_times_sec), len(common_pixels))``.
+    out_name : str
+        Product name, shown as subtitle.
+    colorbar : float or None
+        Symmetric color scale limit [dB]; if None, the maximum absolute difference is used.
+    pol : str
+        Polarization.
+    subswath : str
+        Subswath, shown in the title for non-GRD products.
+    prod_type : str
+        Product type (``"GRD"`` or ``"SLC"``).
+    """
+
+    plt.figure(figsize=(9, 7))
+    if colorbar is not None:
+        vmax = colorbar
+    else:
+        vmax = np.nanmax(np.abs(diff_db))  
+    vmin = -vmax # symmetric range
+
+    plt.pcolormesh(
+        common_pixels,
+        common_times_sec,
+        diff_db,
+        shading="auto",
+        cmap="RdBu",
+        vmax=vmax,
+        vmin=vmin
+    )
+    plt.colorbar(label="Noise difference [dB]")
+    plt.xlabel("Range pixel")
+    plt.ylabel("Azimuth time from acquisition start [s]")
+    if prod_type == "GRD":
+        plt.title(f"Denoising results ({pol})", weight="bold", pad=25, fontsize=14)
+    else:
+        plt.title(f"Denoising results ({pol}, {subswath})", weight="bold", pad=25, fontsize=14)
+    plt.gca().text(0.5, 1.02, f"{out_name}", transform=plt.gca().transAxes, ha="center", fontsize=10, fontfamily="monospace")
+
+    plt.tight_layout()
+    out_path = OUTPUT_DIR / "whole_raster.png"
+    plt.savefig(out_path, dpi=300)
+    plt.close()
+    print(f"Saved plot to {out_path}")
+
+
+def plot_cut_profiles(common_pixels, f_noise, f_denoise, cut_times_sec, product_name, dn, pol, subswath, prod_type):
+    """
+    Plot the 1D noise range profiles and their difference at the given azimuth times.
+
+    One figure per azimuth time is saved as ``cut_profile_<t>s.png``: the top panel shows nominal and reconstructed noise
+    in dB, the bottom panel their difference, where points with a jump larger than 0.01 dB between adjacent pixels
+    (subswath transitions) are masked.
+
+    Parameters
+    ----------
+    common_pixels : numpy.ndarray
+        Range pixel axis.
+    f_noise : scipy.interpolate.RegularGridInterpolator
+        Interpolator of the nominal noise over (azimuth time [s], range pixel).
+    f_denoise : scipy.interpolate.RegularGridInterpolator
+        Interpolator of the reconstructed noise over (azimuth time [s], range pixel).
+    cut_times_sec : numpy.ndarray
+        Azimuth times [s] of the cuts.
+    product_name : str
+        Product name, shown as title.
+    dn : float
+        Digital number scaling factor; noise values are divided by ``dn**2``.
+    pol : str
+        Polarization.
+    subswath : str
+        Subswath, shown in the title for SLC products.
+    prod_type : str
+        Product type (``"GRD"`` or ``"SLC"``).
+    """
+
+    for t_sec in cut_times_sec:
+        # Query points: shape (num_pixels, 2) = (time, pixel)
+        points = np.column_stack([np.full_like(common_pixels, t_sec), common_pixels])
+
+        # Evaluate interpolators
+        noise_profile = f_noise(points) / (dn**2)
+        denoise_profile = f_denoise(points) / (dn**2)
+
+        # Avoid zero values for log10
+        noise_profile[noise_profile == 0] = nan
+        denoise_profile[denoise_profile == 0] = nan
+
+        noise_db = linear_to_db(noise_profile)
+        denoise_db = linear_to_db(denoise_profile)
+        diff_db = denoise_db - noise_db
+
+        # Plot
+        fig, axs = plt.subplots(2, 1, figsize=(8, 9), sharex=True)
+        fig.suptitle(f"{product_name}", fontfamily = "monospace", fontweight="bold")
+
+        # Top panel: original vs denoised
+        axs[0].plot(common_pixels, noise_db, label="noise", linestyle="--")
+        axs[0].plot(common_pixels, denoise_db, label="denoise", alpha=0.7)
+        axs[0].set_ylabel("Noise power [dB]")
+        if prod_type == "SLC":
+            axs[0].set_title(f"Profiles ({pol}, {subswath}) at azimuth ≈ {t_sec:.2f} s")
+        else:
+            axs[0].set_title(f"Profiles ({pol}) at azimuth ≈ {t_sec:.2f} s")
+        axs[0].legend()
+        axs[0].grid(True)
+
+        # Bottom panel: difference (remove swath jumps)
+        deriv = (diff_db[1:] - diff_db[:-1])
+        deriv = np.insert(deriv, 0, 0.0)
+        mask = abs(deriv) > 0.01
+
+        # print(max(abs(deriv)))
+
+        diff_db[mask] = np.nan
+
+        axs[1].plot(common_pixels, diff_db, color="red")
+        axs[1].set_xlabel("Range pixel")
+        axs[1].set_ylabel("Noise difference [dB]")
+        if prod_type == "SLC":
+            axs[1].set_title("Difference")
+        else:
+            axs[1].set_title("Difference (filtering subswath jumps)")
+        axs[1].grid(True)
+
+        plt.tight_layout()
+        out_path = OUTPUT_DIR / f"cut_profile_{int(t_sec)}s.png"
+        plt.savefig(out_path, dpi=300)
+        plt.close()
+        print(f"Saved plot to {out_path}")
+
+
+def main():
+    """
+    Run the command line interface.
+
+    Parses the two noise XML files, checks that subswath, product type and polarization (from the file names) match,
+    computes the errors and generates the plots.
+
+    Raises
+    ------
+    SystemExit
+        If the two files refer to different subswaths, product types or polarizations (exit code 128).
+    """
+
+    parser = argparse.ArgumentParser(
+        description="Compare two Sentinel-1 noise XML files (GRD or SLC)"
+    )
+
+    parser.add_argument("noise_file", type=Path, help="Nominal noise XML")
+    parser.add_argument("denoise_file", type=Path, help="Candidate denoise XML")
+    parser.add_argument("--colorbar", type=float, default=None, help="Colorbar symmetric limit value for plotting")
+
+    args = parser.parse_args()
+
+    OUTPUT_DIR.mkdir(parents=True, exist_ok=True)
+
+    product_name = args.noise_file.parent.parent.parent.stem
+
+    parts_noise = args.noise_file.stem.split("-")
+    parts_denoise = args.denoise_file.stem.split("-")
+    
+    subswath = parts_noise[2].upper()
+    prod_type = parts_noise[3].upper()
+    pol = parts_noise[4].upper()
+
+    subswath_d = parts_denoise[2].upper()
+    prod_type_d = parts_denoise[3].upper()
+    pol_d = parts_denoise[4].upper()
+
+    if subswath != subswath_d:
+        print(f"Incompatible subswaths: {subswath} for noise file, {subswath_d} for denoise file")
+        sys.exit(128)
+    elif prod_type != prod_type_d:
+        print(f"Incompatible product type: {prod_type} for noise file, {prod_type_d} for denoise file")
+        sys.exit(128)
+    elif pol != pol_d:
+        print(f"Incompatible polarization: {pol} for noise file, {pol_d} for denoise file")
+        sys.exit(128)
+
+    mode = product_name[4:6]
+    dn = DN_VALUES.get((mode, prod_type))
+
+    print(f"Mode: {mode}\nProduct type: {prod_type}")
+
+    common_pixels, common_times_sec, diff_db, f1, f2, mean_error, max_error = compare_noise(args.noise_file, args.denoise_file, dn)
+    print(f"Absolute mean error on denoising: {mean_error:.2e} dB")
+    print(f"Maximum absolute error on denoising: {max_error:.2e} dB")
+
+    plot_2d(common_pixels, common_times_sec, diff_db, product_name, args.colorbar, pol, subswath, prod_type)
+
+    # azimuth cut times in seconds
+    cut_times_sec = np.linspace(1, max(common_times_sec), 6, endpoint=False)[1:]  
+    plot_cut_profiles(common_pixels, f1, f2, cut_times_sec, product_name, dn, pol, subswath, prod_type)
+
+
+if __name__ == "__main__":
+    main()

s1_denoise-0.0.0.dev0/src/s1_denoise/__init__.py

@@ -0,0 +1,18 @@
+# SPDX-FileCopyrightText: Aresys S.r.l. <info@aresys.it>
+# SPDX-License-Identifier: MIT
+
+"""
+Configuration and constants module for the S1 L1N processor.
+
+Defines the package version, the processor name used in log records, the reference slant range ``R_REF``
+(in meters) and the noise pulse outlier detection bounds ``NOISE_LOWER_BOUND`` / ``NOISE_UPPER_BOUND`` (in dB).
+"""
+
+__version__ = '0.0.0.dev0'
+
+PROCESSOR_NAME = 's1_l1n_processor'
+
+# CONSTANTS
+R_REF = 8e5  # [m]
+
+NOISE_LOWER_BOUND, NOISE_UPPER_BOUND = 3.0, 1.0  # [dB] bounds for noise pulse outlier detection

s1_denoise-0.0.0.dev0/src/s1_denoise/__main__.py

@@ -0,0 +1,36 @@
+# SPDX-FileCopyrightText: Aresys S.r.l. <info@aresys.it>
+# SPDX-License-Identifier: MIT
+
+"""Entry point module for the denoising command-line interface."""
+
+import sys
+import click
+
+from s1_denoise.cli import denoise_cli
+
+def main():
+    """
+    Run the denoising CLI, mapping Click errors to exit code 128.
+
+    The command is invoked with ``standalone_mode=False`` so that Click usage errors and user aborts are
+    reported here and terminate the processor with exit code 128 instead of Click's default exit codes.
+
+    Raises
+    ------
+    SystemExit
+        With code 128 on a Click usage error or abort; exit codes raised by the CLI itself are propagated.
+    """
+    # standalone_mode=False lets click errors (e.g. invalid arguments) propagate,
+    # so that the processor exits with code 128 also in these cases
+    try:
+        denoise_cli(standalone_mode=False)
+    except click.ClickException as e:
+        e.show()
+        sys.exit(128)
+    except click.Abort:
+        click.echo("Aborted!", err=True)
+        sys.exit(128)
+
+
+if __name__ == "__main__":
+    main()