GlacioHack · belletva · Jun 17, 2026 · Apr 9, 2026 · Apr 9, 2026 · Apr 10, 2026
diff --git a/doc/source/_workflows/outputs_accuracy/report.html b/doc/source/_workflows/outputs_accuracy/report.html
@@ -2,7 +2,7 @@
 <head><meta charset='UTF-8'><title>Qualify elevation results</title></head>
 <body>
 <h1>Accuracy assessment report — xDEM</h1>
-<p>xDEM version: 0.2.3.dev26</p><p>Date: 16/06/2026 14:38:58</p><p>Computing time: 12.87 seconds</p><h2>Elevation inputs</h2>
+<p>xDEM version: 0.2.3.dev26</p><p>Date: 16/06/2026 18:46:15</p><p>Computing time: 10.76 seconds</p><h2>Elevation inputs</h2>
 <img src='plots/inputs.png' alt='Image PNG' style='width: 100%; height: auto;'>
 <h2>Masked elevation data</h2>
 <img src='plots/masked_elev_map.png' alt='Image PNG' style='width: 100%; height: auto;'>
@@ -16,7 +16,7 @@ <h2>Information about inputs</h2>
 </table>
 </div>
 <h2>Preprocessed elevation data</h2>
-<img src='plots/preprocessed_to_be_aligned_elev_map.png' alt='Image PNG' style='width: 100%; height: auto;'>
+<img src='plots/preprocessed_to_be_aligned_elev_map.png' alt='Image PNG' style='width: 100%; height: auto; justify-content: center'>
 <div style='clear: both; margin-bottom: 30px;'>
 <h2>Coregistration user configuration</h2>
 <table border='1' cellspacing='0' cellpadding='5'>
@@ -52,9 +52,9 @@ <h2>Statistics</h2>
 <tr><td>Valid count</td><td>1133607</td><td>1145221</td><td>1133607</td></tr>
 </table>
 <h2>Elevation differences</h2>
-<img src='plots/diff_elev_diff_coreg_map.png' alt='Image PNG' style='width: 100%; height: auto'>
+<img src='plots/diff_elev_diff_coreg_map.png' alt='Image PNG' style='width: 100%; height: auto; justify-content: center'>
 <h2>Differences histogram</h2>
-<img src='plots/elev_diff_histo.png' alt='Image PNG' style='width: 100%; height: auto'>
+<img src='plots/elev_diff_histo.png' alt='Image PNG' style='width: 100%; height: auto; justify-content: center'>
 
          </body>
          </html>
diff --git a/doc/source/_workflows/outputs_accuracy/report.pdf b/doc/source/_workflows/outputs_accuracy/report.pdf
diff --git a/doc/source/_workflows/outputs_topo/report.html b/doc/source/_workflows/outputs_topo/report.html
@@ -2,7 +2,7 @@
 <head><meta charset='UTF-8'><title>Topographic summary results</title></head>
 <body>
 <h1>Topography summary report — xDEM</h1>
-<p>xDEM version: 0.2.3.dev26</p><p>Date: 16/06/2026 14:39:12</p><p>Computing time: 7.20 seconds</p><h2>Elevation input</h2>
+<p>xDEM version: 0.2.3.dev26</p><p>Date: 16/06/2026 18:46:25</p><p>Computing time: 6.03 seconds</p><h2>Elevation input</h2>
 <img src='plots/elev_map.png' alt='Image PNG' style='width: 100%; height: auto;'>
 <h2>Masked elevation data</h2>
 <img src='plots/masked_elev_map.png' alt='Image PNG' style='width: 100%; height: auto;'>

diff --git a/doc/source/_workflows/outputs_topo/report.pdf b/doc/source/_workflows/outputs_topo/report.pdf
diff --git a/doc/source/cli_accuracy.md b/doc/source/cli_accuracy.md
@@ -301,19 +301,16 @@ Tree of outputs for level 1:
 - root
   ├─ tables
   │   ├─ [aligned_elev_stats.csv] (if coregistration)
-  │   ├─ [diff_elev_before_coreg_stats.csv] (if coregistration)
   │   ├─ [diff_elev_after_coreg_stats.csv] (if coregistration)
-  │   ├─ reference_elev_stats.csv
-  │   └─ to_be_aligned_elev_stats.csv
+  │   ├─ [diff_elev_before_coreg_stats.csv] (if coregistration)
+  │   ├─ [diff_elev_without_coreg_stats.csv] (if no coregistration)
 - ├─ plots
-  │   ├─ inputs.png
-  │   ├─ [masked_elev_map.png] (if `path_to_mask` is given in input)
-  │   ├─ [preprocessed_to_be_aligned_elev_map.png or preprocessed_reference_elev_map.png] (if sampling_grid)
-  │   ├─ [diff_elev_before_coreg_map.png] (if coregistration)
-  │   ├─ [diff_elev_after_coreg_map.png] (if coregistration)
-  │   ├─ [diff_elev_before_after_hist.png] (if coregistration)
+  │   ├─ [diff_elev_diff_coreg_map.png] (if coregistration)
   │   ├─ [diff_elev_without_coreg_map.png] (if no coregistration)
-  │   └─ [elev_diff_histo.png] (if coregistration))
+  │   ├─ [elev_diff_histo.png] (if coregistration))
+- │   ├─ inputs.png
+  │   ├─ [masked_elev_map.png] (if `path_to_mask` is given in input)
+  │   └─ [preprocessed_to_be_aligned_elev_map.png or preprocessed_reference_elev_map.png] (if sampling_grid)
   ├─ [rasters
   │   └─ [aligned_elev.tif] (if coregistration)
   ├─ report.html
@@ -327,28 +324,28 @@ Tree of outputs for level 2:
 - root
   ├─ tables
   │   ├─ [aligned_elev_stats.csv] (if coregistration)
-  │   ├─ [diff_elev_before_coreg_stats.csv] (if coregistration)
   │   ├─ [diff_elev_after_coreg_stats.csv] (if coregistration)
+  │   ├─ [diff_elev_before_coreg_stats.csv] (if coregistration)
+  │   ├─ [diff_elev_without_coreg_stats.csv] (if no coregistration)
   │   ├─ reference_elev_stats.csv
   │   └─ to_be_aligned_elev_stats.csv
   ├─ plots
-  │   ├─ inputs.png
-  │   ├─ [masked_elev_map.png] (if `path_to_mask` is given in input)
-  │   ├─ [preprocessed_to_be_aligned_elev_map.png or preprocessed_reference_elev_map.png] (if sampling_grid)
-  │   ├─ [diff_elev_before_coreg_map.png] (if coregistration)
   │   ├─ [diff_elev_after_coreg_map.png] (if coregistration)
+  │   ├─ [diff_elev_before_coreg_map.png] (if coregistration)
   │   ├─ [diff_elev_coreg_tba_map.png]  (if coregistration)
-  │   ├─ [diff_elev_before_after_hist.png] (if coregistration)
   │   ├─ [diff_elev_without_coreg_map.png] (if no coregistration)
-  │   └─ [elev_diff_histo.png] (if coregistration))
+  │   ├─ [elev_diff_histo.png] (if coregistration))
+- │   ├─ inputs.png
+  │   ├─ [masked_elev_map.png] (if `path_to_mask` is given in input)
+  │   └─ [preprocessed_to_be_aligned_elev_map.png or preprocessed_reference_elev_map.png] (if sampling_grid)
   ├─ rasters
-  │   ├─ [reference_elev_reprojected.tif] (if grid resampling)
-  │   ├─ [to_be_aligned_elev_reprojected.tif] (if grid resampling)
   │   ├─ [aligned_elev.tif] (if coregistration)
-  │   ├─ [diff_elev_before_coreg_map.tif] (if coregistration)
   │   ├─ [diff_elev_after_coreg_map.tif] (if coregistration)
+  │   ├─ [diff_elev_before_coreg_map.tif] (if coregistration)
   │   ├─ [diff_elev_coreg_tba_map.tif] (if coregistration)
-  │   └─ [diff_elev_without_coreg_map.tif] (if no coregistration)
+  │   ├─ [diff_elev_without_coreg_map.tif] (if no coregistration)
+  │   ├─ [reference_elev_reprojected.tif] (if grid resampling)
+  │   └─ [to_be_aligned_elev_reprojected.tif] (if grid resampling)
   ├─ report.html
   ├─ [report.pdf] (if `generate_pdf` if `True`)
   └─ used_config.yaml

diff --git a/tests/test_workflows/test_accuracy.py b/tests/test_workflows/test_accuracy.py
@@ -125,7 +125,12 @@ def test_run(get_accuracy_inputs_test, tmp_path, level, generated_pdf):
 
     assert Path(tmp_path / "tables").joinpath("aligned_elev_stats.csv").exists()
 
-    assert Path(tmp_path / "plots").joinpath("diff_elev_diff_coreg_map.png").exists()
+    if level == 1:
+        assert Path(tmp_path / "plots").joinpath("diff_elev_diff_coreg_map.png").exists()
+    else:
+        assert Path(tmp_path / "plots").joinpath("diff_elev_before_coreg_map.png").exists()
+        assert Path(tmp_path / "plots").joinpath("diff_elev_after_coreg_map.png").exists()
+
     assert Path(tmp_path / "plots").joinpath("elev_diff_histo.png").exists()
     assert Path(tmp_path / "plots").joinpath("masked_elev_map.png").exists()
     assert Path(tmp_path / "plots").joinpath("inputs.png").exists()

diff --git a/xdem/workflows/accuracy.py b/xdem/workflows/accuracy.py
@@ -97,7 +97,6 @@ def _load_data(self) -> tuple[float, float]:
             title_dem_right="To-be-aligned elevation",
             vmin=vmin,
             vmax=vmax,
-            cbar_title=f"Elevation ({self.reference_elev.crs.linear_units})",
         )
         if ref_mask is not None or tba_mask is not None:
             if ref_mask is not None:
@@ -115,7 +114,6 @@ def _load_data(self) -> tuple[float, float]:
                 title_dem_right="Masked terrain for to-be-aligned elevation",
                 vmin=vmin,
                 vmax=vmax,
-                cbar_title=f"Elevation ({self.reference_elev.crs.linear_units})",
             )
 
         self.dico_to_show = [
@@ -214,7 +212,6 @@ def _prepare_datas(self, vmin: float, vmax: float) -> None:
                 filename="preprocessed_to_be_aligned_elev_map",
                 vmin=vmin,
                 vmax=vmax,
-                cbar_title=f"Elevation ({self.to_be_aligned_elev.crs.linear_units})",
             )
         else:
             self.reference_elev = self.reference_elev.crop(coord_intersection)
@@ -224,7 +221,6 @@ def _prepare_datas(self, vmin: float, vmax: float) -> None:
                 filename="preprocessed_reference_elev_map",
                 vmin=vmin,
                 vmax=vmax,
-                cbar_title=f"Elevation ({self.reference_elev.crs.linear_units})",
             )
 
         if self.level > 1:
@@ -342,32 +338,66 @@ def run(self) -> None:
                 self.stats_after["median"] + 3 * self.stats_after["nmad"],
             )
 
-            self.generate_plot(
-                dem=self.diff_before,
-                title="Difference between To-be-align and Reference elevation\n(before coregistration)",
-                filename="diff_elev_diff_coreg_map",
-                dem_right=self.diff_after,
-                title_dem_right="Difference between Aligned and Reference elevation\n(after coregistration)",
-                vmin=vmin_diff,
-                vmax=vmax_diff,
-                cmap="RdBu",
-                cbar_title=f"Elevation differences ({self.diff_before.crs.linear_units})",
-            )
+            if self.level == 1:
+                self.generate_plot(
+                    dem=self.diff_before,
+                    title="Difference between To-be-aligned and Reference elevation\n(before coregistration)",
+                    filename="diff_elev_diff_coreg_map",
+                    dem_right=self.diff_after,
+                    title_dem_right="Difference between Aligned and Reference elevation\n(after coregistration)",
+                    vmin=vmin_diff,
+                    vmax=vmax_diff,
+                    cmap="RdBu",
+                )
+            else:
+                self.generate_plot_with_profiles(
+                    dem=self.diff_before,
+                    title="Difference between To-be-aligned and Reference elevation\n(before coregistration)",
+                    filename="diff_elev_before_coreg_map",
+                    vmin=vmin_diff,
+                    vmax=vmax_diff,
+                    cmap="RdBu",
+                )
+
+                self.generate_plot_with_profiles(
+                    dem=self.diff_after,
+                    title="Difference between Aligned and Reference elevation\n(after coregistration)",
+                    filename="diff_elev_after_coreg_map",
+                    vmin=vmin_diff,
+                    vmax=vmax_diff,
+                    cmap="RdBu",
+                )
 
-            if self.level > 1:
                 self.diff_coreg_tba = aligned_elev.reproject(self.to_be_aligned_elev) - self.to_be_aligned_elev
 
-                self.generate_plot(
+                self.generate_plot_with_profiles(
                     dem=self.diff_coreg_tba,
-                    title="Difference between Aligned and To-be-align elevation\n(no coregistration)",
+                    title="Difference between Aligned and To-be-aligned elevation\n(after coregistration)",
                     filename="diff_elev_coreg_tba_map",
                     cmap="RdBu",
-                    cbar_title=f"Elevation differences ({self.diff_after.crs.linear_units})",
                 )
         else:
             self.diff = self.to_be_aligned_elev - self.reference_elev
             self.stats = self.diff.get_stats(stats_keys)
             vmin, vmax = -(self.stats["median"] + 3 * self.stats["nmad"]), self.stats["median"] + 3 * self.stats["nmad"]
+            if self.level == 1:
+                self.generate_plot(
+                    self.diff,
+                    title="Difference between To-be-aligned and Reference elevation",
+                    filename="diff_elev_without_coreg_map",
+                    vmin=vmin,
+                    vmax=vmax,
+                    cmap="RdBu",
+                )
+            else:
+                self.generate_plot_with_profiles(
+                    dem=self.diff,
+                    title="Difference between To-be-aligned and Reference elevation",
+                    filename="diff_elev_without_coreg_map",
+                    vmin=vmin,
+                    vmax=vmax,
+                    cmap="RdBu",
+                )
             self.generate_plot(
                 self.diff,
                 title="Difference between To-be-align and Reference elevation",
@@ -494,26 +524,32 @@ def print_dict(title: str, dictionary: dict[str, Any]) -> str:
             div_html += "</div>\n"
             return div_html
 
+        def print_png(title: str, width: int = 100) -> str:
+            return (
+                f"<img src='plots/{title}.png' alt='Image PNG' style='width: {width}%; "
+                f"height: auto; justify-content: center'>\n"
+            )
+
         # Metadata: Inputs
         inputs_information = list_dict[0]
         html += print_dict(inputs_information[0], inputs_information[1])
 
         # Plot preprocessed data if did
         if "sampling_grid" in self.config["inputs"] and self.config["inputs"]["sampling_grid"] is not None:
             if self.config["inputs"]["sampling_grid"] == "reference_elev":
-                preprocessed_data = "plots/preprocessed_to_be_aligned_elev_map.png"
+                preprocessed_data = "preprocessed_to_be_aligned_elev_map"
             else:
-                preprocessed_data = "plots/preprocessed_reference_elev_map.png"
+                preprocessed_data = "preprocessed_reference_elev_map"
 
             html += "<h2>Preprocessed elevation data</h2>\n"
-            html += "<img src='" + preprocessed_data + "' alt='Image PNG' style='width: 100%; height: auto;'>\n"
+            html += print_png(preprocessed_data)
 
         # Metadata: Inputs
         for title, dictionary in list_dict[1:]:  # type: ignore
             html += print_dict(title, dictionary)
 
         if self.compute_coreg and self.level > 1:
-            html += "<img src='plots/diff_elev_coreg_tba_map.png' alt='Image PNG' style='width: 100%; height: auto'>\n"
+            html += print_png("diff_elev_coreg_tba_map")
 
         # Statistics table:
         if self.df_stats is not None:
@@ -530,16 +566,19 @@ def print_dict(title: str, dictionary: dict[str, Any]) -> str:
         # Coregistration: Add elevation difference plot and histograms before/after
         if self.compute_coreg:
             html += "<h2>Elevation differences</h2>\n"
-            html += "<img src='plots/diff_elev_diff_coreg_map.png' alt='Image PNG' style='width: 100%; height: auto'>\n"
+            if self.level == 1:
+                html += print_png("diff_elev_diff_coreg_map")
+            else:
+                html += print_png("diff_elev_before_coreg_map")
+                html += print_png("diff_elev_after_coreg_map")
 
             html += "<h2>Differences histogram</h2>\n"
-            html += "<img src='plots/elev_diff_histo.png' alt='Image PNG' style='width: 100%; height: auto'>\n"
+            html += print_png("elev_diff_histo")
 
         else:
             html += "<h2>Elevation differences</h2>\n"
-            html += (
-                "<img src='plots/diff_elev_without_coreg_map.png' alt='Image PNG' style='width: 100%; height: auto'>\n"
-            )
+            html += print_png("diff_elev_without_coreg_map")
+
         html += """
          </body>
          </html>