wolfhece 2.2.28__py3-none-any.whl → 2.2.30__py3-none-any.whl

wolfhece/report/compare_arrays.py

@@ -4,13 +4,17 @@ import numpy.ma as ma
 from pathlib import Path
 import matplotlib.pyplot as plt
 from enum import Enum
-from scipy.ndimage import label, sum_labels
+from scipy.ndimage import label, sum_labels, find_objects
 import pymupdf as pdf
+import wx
+from tqdm import tqdm
+from matplotlib import use, get_backend
 
 from .common import A4_rect, rect_cm, list_to_html, list_to_html_aligned, get_rect_from_text
 from .common import inches2cm, pts2cm, cm2pts, cm2inches, DefaultLayoutA4, NamedTemporaryFile, pt2inches, TemporaryDirectory
 from ..wolf_array import WolfArray, header_wolf, vector, zone, Zones, wolfvertex as wv, wolfpalette
 from ..PyTranslate import _
+from .pdf import PDFViewer
 
 class ArrayDifferenceLayout(DefaultLayoutA4):
     """
@@ -98,8 +102,16 @@ class ArrayDifference():
 
         self._background = 'IGN'
 
+        self._contour = None
+        self._external_border = None
+
     @property
     def contour(self) -> vector:
+        """ Get the contour of the difference part. """
+
+        if self._contour is not None and isinstance(self._contour, vector):
+            return self._contour
+
         ret = self.reference.suxsuy_contour(abs=True)
         ret = ret[2]
 
@@ -111,8 +123,11 @@ class ArrayDifference():
     @property
     def external_border(self) -> vector:
         """
-        Get the bounds of the reference array.
+        Get the bounds of the difference part.
         """
+        if self._external_border is not None and isinstance(self._external_border, vector):
+            return self._external_border
+
         ret = vector(name=_("External border"))
         (xmin, xmax), (ymin, ymax) = self.reference.get_bounds()
         ret.add_vertex(wv(xmin, ymin))
@@ -194,6 +209,13 @@ class ArrayDifference():
                                         Walonmap= True,
                                         cat = 'IMAGERIE/ORTHO_2022_ETE',
                                         )
+        else:
+            self.reference.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
+                                        first_mask_data=False, with_legend=False,
+                                        update_palette= False,
+                                        Walonmap= False,
+                                        )
+
 
         self.reference.array.mask[:,:] = old_mask
 
@@ -225,11 +247,23 @@ class ArrayDifference():
         new = WolfArray(srcheader=h)
         new.array.mask[:,:] = True
 
-        new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
-                                       first_mask_data=False, with_legend=False,
-                                       update_palette= False,
-                                       Walonmap= True,
-                                       cat = 'IMAGERIE/ORTHO_2022_ETE')
+        if self._background.upper() == 'IGN' or self._background.upper() == 'NGI':
+            new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
+                                        first_mask_data=False, with_legend=False,
+                                        update_palette= False,
+                                        IGN= True,
+                                        cat = 'orthoimage_coverage')
+        elif self._background.upper() == 'WALONMAP':
+            new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
+                                        first_mask_data=False, with_legend=False,
+                                        update_palette= False,
+                                        Walonmap= True,
+                                        cat = 'IMAGERIE/ORTHO_2022_ETE')
+        else:
+            new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
+                                        first_mask_data=False, with_legend=False,
+                                        update_palette= False,
+                                        Walonmap= False)
 
         self.external_border.plot_matplotlib(ax=ax)
 
@@ -295,7 +329,8 @@ class ArrayDifference():
 
         fig, ax = figax
 
-        ax.hist(self.reference.array.compressed(), density=density, alpha = alpha, **kwargs)
+        data = self.reference.array.compressed()
+        ax.hist(data, bins = min(100, int(len(data)/4)), density=density, alpha = alpha, **kwargs)
         # ax.set_xlabel("Value")
         # ax.set_ylabel("Frequency")
         return fig, ax
@@ -309,7 +344,8 @@ class ArrayDifference():
 
         fig, ax = figax
 
-        ax.hist(self.to_compare.array.compressed(), density=density, alpha = alpha, **kwargs)
+        data = self.to_compare.array.compressed()
+        ax.hist(data, bins= min(100, int(len(data)/4)), density=density, alpha = alpha, **kwargs)
         # ax.set_xlabel("Value")
         # ax.set_ylabel("Frequency")
         return fig, ax
@@ -349,7 +385,7 @@ class ArrayDifference():
         fig, ax = figax
 
         difference_data = self.difference.array.compressed()
-        ax.hist(difference_data, density=density, alpha=alpha, **kwargs)
+        ax.hist(difference_data, bins= min(100, int(len(difference_data)/4)), density=density, alpha=alpha, **kwargs)
 
         # ax.set_xlabel("Value")
         # ax.set_ylabel("Frequency")
@@ -396,7 +432,6 @@ class ArrayDifference():
                 if 'Histogram' in key:
                     fig.tight_layout()
 
-
                 # convert canvas to PNG and insert it into the PDF
                 temp_file = NamedTemporaryFile(delete=False, suffix='.png')
                 fig.savefig(temp_file, format='png', bbox_inches='tight', dpi=self._dpi)
@@ -419,6 +454,7 @@ class ArrayDifference():
     def create_report(self, output_file: str | Path = None) -> Path:
         """ Create a page report for the array difference. """
 
+        from time import sleep
         if output_file is None:
             output_file = Path(f"array_difference_{self.index}.pdf")
 
@@ -429,6 +465,7 @@ class ArrayDifference():
         layout.create_report()
         self._complete_report(layout)
         layout.save_report(output_file)
+        sleep(0.2)  # Ensure the file is saved before returning
 
         return output_file
 
@@ -459,6 +496,10 @@ class CompareArrays:
 
         self.difference_parts:dict[int, ArrayDifference] = {}
 
+        self._pdf_path = None
+
+        self._background = 'IGN'
+
     @property
     def difference(self) -> WolfArray:
 
@@ -467,6 +508,24 @@ class CompareArrays:
 
         return self.array_to_compare - self.array_reference
 
+    def get_zones(self):
+        """
+        Get a Zones object containing the differences.
+        """
+
+        ret_zones = Zones()
+        exterior = zone(name=_("External border"))
+        contours = zone(name=_("Contours"))
+
+        ret_zones.add_zone(exterior, forceparent=True)
+        ret_zones.add_zone(contours, forceparent=True)
+
+        for diff in self.difference_parts.values():
+            exterior.add_vector(diff.external_border, forceparent=True)
+            contours.add_vector(diff.contour, forceparent=True)
+
+        return ret_zones
+
     def plot_position(self, figax:tuple[plt.Figure, plt.Axes]=None) -> tuple[plt.Figure, plt.Axes]:
         """
         Plot the reference array with a background.
@@ -487,14 +546,27 @@ class CompareArrays:
         new = WolfArray(srcheader=h)
         new.array.mask[:,:] = True
 
+        if self._background.upper() == 'IGN' or self._background.upper() == 'NGI':
 
-        new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
-                                    first_mask_data=False, with_legend=False,
-                                    update_palette= False,
-                                    IGN= True,
-                                    cat = 'orthoimage_coverage',
-                                    )
-
+            new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
+                                        first_mask_data=False, with_legend=False,
+                                        update_palette= False,
+                                        IGN= True,
+                                        cat = 'orthoimage_coverage',
+                                        )
+        elif self._background.upper() == 'WALONMAP':
+            new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
+                                        first_mask_data=False, with_legend=False,
+                                        update_palette= False,
+                                        Walonmap= True,
+                                        cat = 'IMAGERIE/ORTHO_2022_ETE',
+                                        )
+        else:
+            new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
+                                        first_mask_data=False, with_legend=False,
+                                        update_palette= False,
+                                        Walonmap= False,
+                                        )
         return fig, ax
 
     def plot_cartoweb(self, figax:tuple[plt.Figure, plt.Axes]=None) -> tuple[plt.Figure, plt.Axes]:
@@ -517,14 +589,20 @@ class CompareArrays:
         new = WolfArray(srcheader=h)
         new.array.mask[:,:] = True
 
-
-        new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
-                                    first_mask_data=False, with_legend=False,
-                                    update_palette= False,
-                                    Cartoweb= True,
-                                    cat = 'overlay',
-                                    )
-
+        if self._background.upper() == 'IGN' or self._background.upper() == 'NGI':
+            new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
+                                        first_mask_data=False, with_legend=False,
+                                        update_palette= False,
+                                        Cartoweb= True,
+                                        cat = 'overlay',
+                                        )
+        else:
+            new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
+                                        first_mask_data=False, with_legend=False,
+                                        update_palette= False,
+                                        Cartoweb= False,
+                                        cat = 'overlay',
+                                        )
         return fig, ax
 
     def plot_topo_grey(self, figax:tuple[plt.Figure, plt.Axes]=None) -> tuple[plt.Figure, plt.Axes]:
@@ -547,14 +625,20 @@ class CompareArrays:
         new = WolfArray(srcheader=h)
         new.array.mask[:,:] = True
 
-
-        new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
-                                    first_mask_data=False, with_legend=False,
-                                    update_palette= False,
-                                    Cartoweb= True,
-                                    cat = 'topo_grey',
-                                    )
-
+        if self._background.upper() == 'IGN' or self._background.upper() == 'NGI':
+            new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
+                                        first_mask_data=False, with_legend=False,
+                                        update_palette= False,
+                                        Cartoweb= True,
+                                        cat = 'topo_grey',
+                                        )
+        else:
+            new.plot_matplotlib(figax=figax, figsize = self.default_size_arrays,
+                                        first_mask_data=False, with_legend=False,
+                                        update_palette= False,
+                                        Cartoweb= False,
+                                        cat = 'topo_grey',
+                                        )
         return fig, ax
 
     def plot_reference(self, figax:tuple[plt.Figure, plt.Axes]=None) -> tuple[plt.Figure, plt.Axes]:
@@ -568,7 +652,7 @@ class CompareArrays:
 
         self.array_reference.plot_matplotlib(figax=figax, figsize = self.default_size_arrays, first_mask_data=False, with_legend=True, update_palette= False)
 
-        for diff in self.difference_parts.values():
+        for diff in tqdm(self.difference_parts.values(), desc="Plotting external borders"):
             diff.external_border.plot_matplotlib(ax=ax)
 
         return fig, ax
@@ -584,7 +668,7 @@ class CompareArrays:
 
         self.array_to_compare.plot_matplotlib(figax=figax, figsize = self.default_size_arrays, first_mask_data=False, with_legend=True, update_palette= False)
 
-        for diff in self.difference_parts.values():
+        for diff in tqdm(self.difference_parts.values(), desc="Plotting contours"):
             diff.contour.plot_matplotlib(ax=ax)
 
         return fig, ax
@@ -606,11 +690,22 @@ class CompareArrays:
         diff.plot_matplotlib(figax=figax, figsize = self.default_size_arrays, first_mask_data=False, with_legend=True, update_palette= False)
         return fig, ax
 
-    def localize_differences(self, threshold: float = 0.0) -> np.ndarray:
+    def localize_differences(self, threshold: float = 0.0,
+                             ignored_patche_area:float = 1.) -> np.ndarray:
+        """ Localize the differences between the two arrays and label them.
+
+        :param threshold: The threshold value to consider a difference significant.
+        :param ignored_patche_area: The area of patches to ignore (in m²).
+        """
+
+        assert threshold >= 0, "Threshold must be a non-negative value."
 
         labeled_array = self.difference.array.data.copy()
         labeled_array[self.array_reference.array.mask] = 0
 
+        # apply threshold
+        labeled_array[np.abs(labeled_array) < threshold] = 0
+
         self.labeled_array, self.num_features = label(labeled_array)
 
         self.nb_cells = []
@@ -620,13 +715,29 @@ class CompareArrays:
 
         self.nb_cells.sort(key=lambda x: x[0], reverse=True)
 
+        # find features where nb_cells is lower than ignored_patche_area / (dx * dy)
+        ignored_patche_cells = int(ignored_patche_area / (self.array_reference.dx * self.array_reference.dy))
+        self.last_features = self.num_features
+        for idx, (nb_cell, idx_feature) in enumerate(self.nb_cells):
+            if nb_cell <= ignored_patche_cells:
+                self.last_features = idx
+                break
+
+        all_slices = find_objects(self.labeled_array)
+
+        logging.info(f"Total number of features found: {self.last_features}")
+
         # find xmin, ymin, xmax, ymax for each feature
-        for idx_feature in range(1, self.num_features + 1):
-            mask = self.labeled_array == idx_feature
-            i, j = np.where(mask)
+        for idx_feature, slices in tqdm(zip(range(1, self.num_features+1), all_slices), desc="Processing features", unit="feature"):
+            mask = self.labeled_array[slices] == idx_feature
+            nb_in_patch = np.count_nonzero(mask)
+
+            if nb_in_patch <= ignored_patche_cells:
+                logging.debug(f"Feature {idx_feature} has too few cells ({np.count_nonzero(mask)}) and will be ignored.")
+                continue
 
-            imin, imax = i.min(), i.max()
-            jmin, jmax = j.min(), j.max()
+            imin, imax = slices[0].start, slices[0].stop - 1
+            jmin, jmax = slices[1].start, slices[1].stop - 1
 
             imin = int(max(imin - 1, 0))
             imax = int(min(imax + 1, self.labeled_array.shape[0] - 1))
@@ -635,17 +746,25 @@ class CompareArrays:
 
             ref_crop = self.array_reference.crop(imin, jmin, imax-imin+1, jmax-jmin+1)
             to_compare_crop = self.array_to_compare.crop(imin, jmin, imax-imin+1, jmax-jmin+1)
+            label_crop = self.labeled_array[imin:imax+1, jmin:jmax+1].copy()
 
-            ref_crop.array.mask[:,:] = ~mask[imin:imax+1, jmin:jmax+1]
-            ref_crop.set_nullvalue_in_mask()
 
-            to_compare_crop.array.mask[:,:] = ~mask[imin:imax+1, jmin:jmax+1]
+            to_compare_crop.array.mask[:,:] = ref_crop.array.mask[:,:] = self.labeled_array[imin:imax+1, jmin:jmax+1] != idx_feature
+
+            ref_crop.set_nullvalue_in_mask()
             to_compare_crop.set_nullvalue_in_mask()
+            label_crop[label_crop != idx_feature] = 0
 
-            ref_crop.count()
-            to_compare_crop.nbnotnull = ref_crop.nbnotnull
+            ref_crop.nbnotnull = nb_in_patch
+            to_compare_crop.nbnotnull = nb_in_patch
 
-            self.difference_parts[idx_feature] = ArrayDifference(ref_crop, to_compare_crop, idx_feature, self.labeled_array[imin:imax+1, jmin:jmax+1].copy())
+            self.difference_parts[idx_feature] = ArrayDifference(ref_crop, to_compare_crop, idx_feature, label_crop)
+
+        assert self.last_features == len(self.difference_parts), \
+            f"Last feature index {self.last_features} does not match the number of differences found"
+
+        self.num_features = self.last_features
+        logging.info(f"Number of features after filtering: {self.num_features}")
 
         return self.labeled_array
 
@@ -721,9 +840,10 @@ class CompareArrays:
         surf = self.array_reference.dx * self.array_reference.dy
 
         # Extract the number of cells for each feature
-        nb_cells = [item[0] * surf for item in self.nb_cells]
+        nb_cells = [item[0] * surf for item in self.nb_cells[:self.last_features]]
 
-        ax.hist(nb_cells, density=density, alpha=alpha, **kwargs)
+        if len(nb_cells) > 0:
+            ax.hist(nb_cells, bins= min(100, int(len(nb_cells)/4)), density=density, alpha=alpha, **kwargs)
 
         ax.set_title(_("Histogram of surface in each feature [m²]"))
 
@@ -748,13 +868,22 @@ class CompareArrays:
 
         fig, ax = figax
 
-        # Calculate the difference between the two arrays
-        diff = self.difference
+        # # Calculate the difference between the two arrays
+        # diff = self.difference
 
-        volumes = np.bincount(self.labeled_array.ravel(), weights=diff.array.ravel()) * self.array_reference.dx * self.array_reference.dy
+        volumes = []
+        for idx in tqdm(self.nb_cells[:self.last_features], desc="Calculating volumes"):
+            # Get the feature index
+            feature_index = idx[1]
+            part = self.difference_parts[feature_index]
+            # Create a mask for the feature
+            mask = part.label == feature_index
+            # Calculate the volume for this feature
+            volumes.append(np.ma.sum(part.difference.array[mask]) * self.array_reference.dx * self.array_reference.dy)
 
         # Create a histogram of the differences
-        ax.hist(volumes, density=density, alpha=alpha, **kwargs)
+        if len(volumes) > 0:
+            ax.hist(volumes, bins= min(100, int(len(volumes)/4)), density=density, alpha=alpha, **kwargs)
 
         ax.set_title(_("Histogram of net volumes [m³]"))
 
@@ -809,7 +938,9 @@ class CompareArrays:
 
 
 
-    def create_report(self, output_file: str | Path = None, append_all_differences: bool = True) -> None:
+    def create_report(self, output_file: str | Path = None,
+                      append_all_differences: bool = True,
+                      nb_max_differences:int = -1) -> None:
         """ Create a page report for the array comparison. """
 
         if output_file is None:
@@ -822,6 +953,14 @@ class CompareArrays:
         layout.create_report()
         self._complete_report(layout)
 
+        if nb_max_differences < 0:
+            nb_max_differences = len(self.difference_parts)
+        elif nb_max_differences > len(self.difference_parts):
+            logging.warning(f"Requested {nb_max_differences} differences, but only {len(self.difference_parts)} are available. Using all available differences.")
+        elif nb_max_differences < len(self.difference_parts):
+            logging.info(f"Limiting to {nb_max_differences} differences.")
+
+        features_to_treat = [feature[1] for feature in self.nb_cells[:nb_max_differences]]
 
         with TemporaryDirectory() as temp_dir:
 
@@ -833,12 +972,83 @@ class CompareArrays:
 
             if append_all_differences:
                 # Add each difference report to the main layout
-                all_pdfs.extend([diff.create_report(Path(temp_dir) / f"array_difference_{idx}.pdf") for idx, diff in self.difference_parts.items()])
+                all_pdfs.extend([self.difference_parts[idx].create_report(Path(temp_dir) / f"array_difference_{idx}.pdf") for idx in tqdm(features_to_treat, desc="Creating individual difference reports")])
 
-            for pdf_file in all_pdfs:
+            for pdf_file in tqdm(all_pdfs, desc="Compiling PDFs"):
                 layout._doc.insert_file(pdf_file)
 
         # create a TOC
         layout._doc.set_toc(layout._doc.get_toc())
 
-        layout.save_report(output_file)
+        layout.save_report(output_file)
+        self._pdf_path = output_file
+
+    @property
+    def pdf_path(self) -> Path:
+        """ Return the path to the generated PDF report. """
+        if hasattr(self, '_pdf_path'):
+            return self._pdf_path
+        else:
+            raise AttributeError("PDF path not set. Please create the report first.")
+
+class CompareArrays_wx(PDFViewer):
+
+    def __init__(self, reference: WolfArray | str | Path,
+                 to_compare: WolfArray | str | Path,
+                 ignored_patche_area:float = 2.0,
+                 nb_max_patches:int = 10,
+                 threshold: float = 0.01,
+                 dpi=200, **kwargs):
+        """ Initialize the Simple Simulation GPU Report Viewer for comparison """
+
+        super(CompareArrays_wx, self).__init__(None, **kwargs)
+
+        use('agg')
+
+        if isinstance(reference, WolfArray) and isinstance(to_compare, WolfArray):
+            if np.any(reference.array.mask != to_compare.array.mask):
+                logging.warning("The masks of the two arrays are not identical. This may lead to unexpected results.")
+                dlg = wx.MessageDialog(self,
+                                       _("The masks of the two arrays are not identical.\nThis may lead to unexpected results.\n\nWe will use the reference mask for the comparison."),
+                                       _("Warning"),
+                                       wx.OK | wx.ICON_WARNING)
+                dlg.ShowModal()
+                dlg.Destroy()
+                to_compare = WolfArray(mold=to_compare)
+                to_compare.array.mask[:,:] = reference.array.mask[:,:]
+
+        self._report = CompareArrays(reference, to_compare)
+
+        self._report._dpi = dpi
+        self._report.localize_differences(threshold=threshold,
+                                          ignored_patche_area=ignored_patche_area)
+        self._report.create_report(nb_max_differences=nb_max_patches)
+
+        # Load the PDF into the viewer
+        if self._report.pdf_path is None:
+            logging.error("No report created. Cannot load PDF.")
+            return
+
+        self.load_pdf(self._report.pdf_path)
+        self.viewer.SetZoom(-1)  # Fit to width
+
+        # Set the title of the frame
+        self.SetTitle("Simple Simulation GPU Comparison Report")
+
+        self.Bind(wx.EVT_CLOSE, self.on_close)
+
+        use('wxagg')
+
+    def on_close(self, event):
+        """ Handle the close event to clean up resources """
+        self.viewer.pdfdoc.pdfdoc.close()
+        self.Destroy()
+
+    def get_zones(self) -> Zones:
+        """
+        Get the zones from the report.
+        """
+        ret = self._report.get_zones()
+        ret.prep_listogl()
+
+        return ret

wolfhece/report/simplesimgpu.py

@@ -11,6 +11,7 @@ from tempfile import NamedTemporaryFile
 from datetime import datetime as dt
 
 import matplotlib.pyplot as plt
+import matplotlib as mpl
 import seaborn as sns
 
 import pymupdf as pdf
@@ -30,6 +31,7 @@ class SimpleSimGPU_Report():
     def __init__(self, sim:SimpleSimulation | Path | str, **kwargs):
         """ Initialize the Simple Simulation GPU Report Viewer """
 
+        self._summary = {}
         self._doc = None
 
         if isinstance(sim, Path):
@@ -38,6 +40,7 @@ class SimpleSimGPU_Report():
             except Exception as e:
                 logging.error(f"Failed to load simulation from path {sim}: {e}")
                 self._sim = None
+                self._summary['errors'] = e
                 return
         elif isinstance(sim, str):
             try:
@@ -45,6 +48,7 @@ class SimpleSimGPU_Report():
             except Exception as e:
                 logging.error(f"Failed to load simulation from string path {sim}: {e}")
                 self._sim = None
+                self._summary['errors'] = e
                 return
         elif not isinstance(sim, SimpleSimulation):
             try:
@@ -52,12 +56,12 @@ class SimpleSimGPU_Report():
             except Exception as e:
                 logging.error(f"Failed to set simulation: {e}")
                 self._sim = None
+                self._summary['errors'] = e
                 return
         else:
             logging.error("Invalid type for simulation. Must be SimpleSimulation, Path, or str.")
             return
 
-        self._summary = {}
         self._summary['warnings'] = self._summary_warnings()
         self._summary['errors'] = self._summary_errors()
 
@@ -272,9 +276,13 @@ class SimpleSimGPU_Report():
         # Plot the histogram of waterdepth adn add it to the PDF
         fig, ax = plt.subplots(figsize=(8, 6))
         # Plot the histogram of water depth
-        ax.hist(sim.h[sim.h > 0.], bins=100, density=True)
-        # ax.set_title('Histogram of Water Depth')
-        ax.set_xlim(0, np.max(sim.h[sim.h > 0.]) * 1.25)  # Set xlim to 110% of max value
+
+        h_min = np.min(sim.h[sim.nap == 1])
+        h_max = np.max(sim.h[sim.nap == 1])
+
+        if h_max > h_min:
+            ax.hist(sim.h[sim.h > 0.], bins=100, density=True)
+            ax.set_xlim(0, h_max)  # Set xlim to 110% of max value
         ax.set_xlabel('Water Depth [m]')
         ax.set_ylabel('Frequency')
 
@@ -301,7 +309,7 @@ class SimpleSimGPU_Report():
         ax.hist(sim.manning[sim.nap == 1], bins=100, density = True)
         # ax.set_title('Histogram of Manning Coefficient')
         ax.set_xlabel('Manning [$\\frac {s} {m^{1/3}} $]')
-        ax.set_xlim(0, np.max(sim.manning[sim.nap == 1]) * 1.25)  # Set xlim to 110% of max value
+        ax.set_xlim(0, np.max(sim.manning[sim.nap == 1]) * 1.1)  # Set xlim to 110% of max value
         ax.set_ylabel('Frequency')
 
         # set font size of the labels
@@ -771,15 +779,21 @@ class SimpleSimGPU_Report():
 
 class SimpleSimGPU_Report_wx(PDFViewer):
 
-    def __init__(self, sim:SimpleSimulation | Path | str, **kwargs):
+    def __init__(self, sim:SimpleSimulation | Path | str, show:bool=False, **kwargs):
         """ Initialize the Simple Simulation GPU Report Viewer """
 
+        mpl.use('Agg')  # Use a non-interactive backend for matplotlib
+
         super(SimpleSimGPU_Report_wx, self).__init__(None, **kwargs)
 
         self._report = SimpleSimGPU_Report(sim, **kwargs)
 
         if self._report._sim is None:
             logging.error("No simulation data available to create report.")
+            dlg = wx.MessageDialog(self, "No simulation data available to create report.\n\nPlease check the errors in the logs.",
+                                  "Error", wx.OK | wx.ICON_ERROR)
+            dlg.ShowModal()
+            dlg.Destroy()
             return
 
         self._report.create_report()
@@ -798,6 +812,11 @@ class SimpleSimGPU_Report_wx(PDFViewer):
 
         self.Bind(wx.EVT_CLOSE, self.on_close)
 
+        if show:
+            self.Show()
+
+        mpl.use('WxAgg')  # Reset matplotlib to use the WxAgg backend for other plots
+
     def on_close(self, event):
         """ Handle the close event of the frame """