wolfhece 2.2.38__py3-none-any.whl → 2.2.40__py3-none-any.whl

wolfhece/PyDraw.py

@@ -94,6 +94,7 @@ try:
     from .images_tiles import ImagesTiles
     from .PyWMS import Alaro_Navigator, get_Alaro_legend
     from .PyPictures import PictureCollection
+    from .irm_alaro import IRM_Alaro, GribFiles, _convert_col2date_str
 
 except ImportError as e:
     print(e)
@@ -1817,7 +1818,7 @@ class Precomputed_DEM_DTM_Dialog(wx.Dialog):
         self._header.origx = float(int(onzoom[0]))
         self._header.origy = float(int(onzoom[2]))
 
-        res = self._resolution.GetSelection()
+        res = self._resolution.GetStringSelection()
         if res == '50cm':
             res = 0.5
         elif res in ['1m', '2m', '5m', '10m']:
@@ -1964,6 +1965,7 @@ class WolfMapViewer(wx.Frame):
     active_drowning: Drowning_victim_Viewer
     active_dike : DikeWolf
     active_picturecollection: PictureCollection
+    active_alaro: IRM_Alaro
 
     active_fig: MplFigViewer
 
@@ -2103,6 +2105,8 @@ class WolfMapViewer(wx.Frame):
         self.menudrowning = None
         self.menudike = None
         self.menupicturecollections = None
+        self.menuqdfidf = None
+        self.menualaro = None
 
         self.alaro_navigator = None
 
@@ -2152,7 +2156,7 @@ class WolfMapViewer(wx.Frame):
         # SIMULATION Hydrologique
 
         self.menu_hydrology = wx.Menu()
-        hydrol = self.menu_hydrology.Append(wx.ID_ANY, _('Open Hydrological model'), _('Hydrological simulation'))
+        hydrol = self.menu_hydrology.Append(wx.ID_ANY, _('Create/Open Hydrological model'), _('Hydrological simulation'))
         self.filemenu.Append(wx.ID_ANY,_('Hydrology'),self.menu_hydrology)
 
         self.filemenu.AppendSeparator()
@@ -2182,6 +2186,7 @@ class WolfMapViewer(wx.Frame):
         createbcmanager2D = self.menucreateobj.Append(wx.ID_ANY, _('Create BC manager Wolf2D...'), _('New BC manager 2D'))
         createpartsystem = self.menucreateobj.Append(wx.ID_ANY, _('Create particle system...'), _('Create a particle system - Lagrangian view'))
         create_acceptability = self.menucreateobj.Append(wx.ID_ANY, _('Create acceptability manager...'), _('Create acceptability manager'))
+        create_inbe = self.menucreateobj.Append(wx.ID_ANY, _('Create INBE manager...'), _('Create INBE manager'))
         createdrowning = self.menucreateobj.Append(wx.ID_ANY, _('Create a drowning...'),_('Create a drowning'))
         createdike = self.menucreateobj.Append(wx.ID_ANY, _('Create dike...'), _('New dike'))
 
@@ -2403,6 +2408,8 @@ class WolfMapViewer(wx.Frame):
         self.active_laz = None
         self.active_dike = None
         self.active_picturecollection = None
+        self.active_qdfidf = None
+        self.active_alaro = None
 
         self.active_fig = None
         self.active_drowning = None
@@ -2526,6 +2533,8 @@ class WolfMapViewer(wx.Frame):
         self._distances = Zones(mapviewer=self, idx=_('Distances/Areas'))
         self._distances.add_zone(zone(name='memory distances', parent=self._distances))  # distances memory
 
+        self.menu_alaro_forecasts()
+
         # self._wintab = Wintab(self.GetHandle())
 
         # if self._wintab:
@@ -2813,6 +2822,144 @@ class WolfMapViewer(wx.Frame):
             self.Bind(wx.EVT_MENU, self.action_pictcollections, allvisible)
             self.Bind(wx.EVT_MENU, self.action_pictcollections, extract)
 
+    def menu_qdfidf(self):
+        """ Menu for QDF/IDF files """
+
+        if self.menuqdfidf is None:
+            self.menuqdfidf = wx.Menu()
+            self.menubar.Append(self.menuqdfidf, _('QDF/IDF'))
+
+            download = self.menuqdfidf.Append(wx.ID_ANY, _('Download data'), _('Download all QDF/IDF data for Belgium and preprocess it !'))
+            load = self.menuqdfidf.Append(wx.ID_ANY, _('Load data'), _('Load a QDF/IDF for Belgium'))
+            pick = self.menuqdfidf.Append(wx.ID_ANY, _('Pick municipality'), _('Pick a municipality and show its data'))
+            show_tables = self.menuqdfidf.Append(wx.ID_ANY, _('Show tables'), _('Show tables on viewer'))
+            show_plots = self.menuqdfidf.Append(wx.ID_ANY, _('Show plots'), _('Show plots on viewer'))
+            scale = self.menuqdfidf.Append(wx.ID_ANY, _('Scale images'), _('Scale the QDF/IDF images'))
+
+            self.Bind(wx.EVT_MENU, self.action_qdfidf, download)
+            self.Bind(wx.EVT_MENU, self.action_qdfidf, load)
+            self.Bind(wx.EVT_MENU, self.action_qdfidf, pick)
+            self.Bind(wx.EVT_MENU, self.action_qdfidf, show_tables)
+            self.Bind(wx.EVT_MENU, self.action_qdfidf, show_plots)
+            self.Bind(wx.EVT_MENU, self.action_qdfidf, scale)
+
+    def action_qdfidf(self, event: wx.Event):
+        """ Action for QDF/IDF files """
+
+        item = event.GetEventObject().FindItemById(event.GetId())
+        itemlabel = item.ItemLabel
+
+        if self.active_qdfidf is None and itemlabel != _('Load data') and itemlabel != _('Download data'):
+            logging.warning(_('No active QDF/IDF -- Please load data first'))
+            return
+
+        if itemlabel == _('Load data'):
+
+            from .irm_qdf import QDF_Hydrology_Draw
+
+            dlg_dir = wx.DirDialog(self, _('Choose directory with QDF/IDF files'), style=wx.DD_DEFAULT_STYLE)
+            if dlg_dir.ShowModal() == wx.ID_CANCEL:
+                dlg_dir.Destroy()
+                return
+            dirpath = dlg_dir.GetPath()
+            dlg_dir.Destroy()
+
+            pgbar = wx.ProgressDialog(_('Loading QDF/IDF data'), _('Loading data...'), maximum=100, parent=self, style=wx.PD_APP_MODAL | wx.PD_AUTO_HIDE)
+            pgbar.Pulse(_('Loading QDF/IDF data... (estimated time 10 min.)'))
+
+            self.active_qdfidf = QDF_Hydrology_Draw(dirpath, idx = 'QDF/IDF', mapviewer = self)
+
+            pgbar.Update(100, _('QDF/IDF data loaded'))
+            pgbar.Destroy()
+
+            self.treelist.SetFocus()
+
+        elif itemlabel == _('Download data'):
+
+            from .irm_qdf import QDF_Hydrology_Draw
+
+            dlg_dir = wx.DirDialog(self, _('Choose an empty directory to store QDF/IDF files'), style=wx.DD_DEFAULT_STYLE | wx.DD_NEW_DIR_BUTTON)
+            if dlg_dir.ShowModal() == wx.ID_CANCEL:
+                dlg_dir.Destroy()
+                return
+
+            dirpath = dlg_dir.GetPath()
+            dlg_dir.Destroy()
+
+            # Test if the directory is empty
+            if os.listdir(dirpath):
+                logging.error(_('The directory {} is not empty. Please choose an empty directory.').format(dirpath))
+                wx.MessageBox(_('The directory {} is not empty. Please choose an empty directory.').format(dirpath), _('Error'), wx.OK | wx.ICON_ERROR)
+                return
+
+            pgbar = wx.ProgressDialog(_('Downloading QDF/IDF data'), _('Downloading data...'), maximum=100, parent=self, style=wx.PD_APP_MODAL | wx.PD_AUTO_HIDE)
+            pgbar.Pulse(_('Downloading QDF/IDF data...'))
+
+            try:
+                self.active_qdfidf = QDF_Hydrology_Draw(dirpath, idx = 'QDF/IDF', mapviewer = self)
+            except Exception as e:
+                logging.error(_('Error downloading QDF/IDF data: {}').format(e))
+                wx.MessageBox(_('Error downloading QDF/IDF data: {}').format(e), _('Error'), wx.OK | wx.ICON_ERROR)
+
+            pgbar.Update(100, _('QDF/IDF data loaded and processed'))
+            pgbar.Destroy()
+
+            self.treelist.SetFocus()
+
+        elif itemlabel == _('Pick municipality'):
+            self.action = 'pick municipality'
+            logging.info(_('Pick a municipality from the QDF/IDF data'))
+
+        elif itemlabel == _('Show tables'):
+
+            if self.active_qdfidf is not None:
+
+                pgbar = wx.ProgressDialog(_('Loading QDF/IDF data'), _('Loading data...'), maximum=100, parent=self, style=wx.PD_APP_MODAL | wx.PD_AUTO_HIDE)
+                pgbar.Pulse(_('Loading QDF/IDF data...'))
+
+                self.active_qdfidf.show_plot = False
+                self.active_qdfidf.show_table = True
+                logging.info(_('Show tables for QDF/IDF data'))
+                self.Refresh()
+
+                pgbar.Update(100, _('QDF/IDF data loaded'))
+                pgbar.Destroy()
+                self.treelist.SetFocus()
+
+
+        elif itemlabel == _('Show plots'):
+
+            if self.active_qdfidf is not None:
+
+                pgbar = wx.ProgressDialog(_('Loading QDF/IDF data'), _('Loading data...'), maximum=100, parent=self, style=wx.PD_APP_MODAL | wx.PD_AUTO_HIDE)
+                pgbar.Pulse(_('Loading QDF/IDF data...'))
+
+                self.active_qdfidf.show_table = False
+                self.active_qdfidf.show_plot = True
+                logging.info(_('Show tables for QDF/IDF data'))
+                self.Refresh()
+
+                pgbar.Update(100, _('QDF/IDF data loaded'))
+                pgbar.Destroy()
+                self.treelist.SetFocus()
+
+
+        elif itemlabel == _('Scale images'):
+            if self.active_qdfidf is not None:
+                scalefactor = wx.GetTextFromUser(_('Enter the scale factor (default is 1.0)'), _('Scale factor'), '1.0', self)
+                if scalefactor == '':
+                    scalefactor = '1.0'
+                try:
+                    scalefactor = float(scalefactor)
+                except ValueError:
+                    logging.error(_('Invalid scale factor: {}').format(scalefactor))
+                    wx.MessageBox(_('Invalid scale factor: {}').format(scalefactor), _('Error'), wx.OK | wx.ICON_ERROR)
+                    return
+
+            self.active_qdfidf.scale_images(scalefactor)
+
+            self.Refresh()
+
     def action_pictcollections(self, event: wx.Event):
         """ Action for picture collections """
 
@@ -3008,22 +3155,355 @@ class WolfMapViewer(wx.Frame):
 
             self.menuwolf2d.Bind(wx.EVT_MENU, self.Onmenuwolf2d)
 
+    def menu_alaro_forecasts(self):
+        """ Menu for Alaro forecasts """
+        if self.menualaro is None:
+
+            self.menualaro = wx.Menu()
+            self.menubar.Append(self.menualaro, _('Alaro forecasts'))
+            self.menualaro_ftp = wx.Menu()
+            self.menualaro.AppendSubMenu(self.menualaro_ftp, _('FTP server'))
+
+            self.menualaro.Bind(wx.EVT_MENU, self.Onmenualaro)
+
+            self.alaro_show_grid = self.menualaro.Append(wx.ID_ANY, _("Show grid"), _("Show the Alaro grid on the viewer"))
+            self.alaro_list_runs = self.menualaro_ftp.Append(wx.ID_ANY, _("List runs on FTP server"), _("List available Alaro runs on the FTP server"))
+            self.alaro_download_runs = self.menualaro_ftp.Append(wx.ID_ANY, _("Download runs"), _("Download Alaro runs"))
+            self.alaro_download_runs_only_rain = self.menualaro_ftp.Append(wx.ID_ANY, _("Download runs - Only Rain and Temperature"), _("Download Alaro runs"))
+            self.alaro_load_run = self.menualaro.Append(wx.ID_ANY, _("Load run"), _("Load an available Alaro run"))
+            self.alaro_show_run = self.menualaro.Append(wx.ID_ANY, _("Add forecast(s) as arry(s)"), _("Show Alaro forecasts on the viewer as WolfArray - Total precipitation [mm]"))
+            self.alaro_show_all = self.menualaro.Append(wx.ID_ANY, _("Add all forecasts as arrays"), _("Show all available Alaro runs on the viewer as WolfArrays - Total precipitation [mm]"))
+            self.alaro_plot_xy = self.menualaro.Append(wx.ID_ANY, _("Plot for XY"), _("Plot the next selected point"))
+
+            self.menualaro_videos = wx.Menu()
+            self.menualaro.AppendSubMenu(self.menualaro_videos, _('Videos'))
+
+            self.alaro_videos_totprec = self.menualaro_videos.Append(wx.ID_ANY, _("Total precip. [mm] - One run"), _("Create a video of the Total Precipitation forecasts for a specific run"))
+            self.alaro_videos_totprec_all = self.menualaro_videos.Append(wx.ID_ANY, _("Total precip. [mm] - Multiple runs"), _("Create videos of the Total Precipitation forecasts for all runs"))
+            self.alaro_videos_comparison = self.menualaro_videos.Append(wx.ID_ANY, _("Rain intensity [mm/h] - Comparison") , _("Create a video comparing the Rain Intensity forecasts of different runs"))
+
+    def Onmenualaro(self, event: wx.Event):
+        """ Action for Alaro forecasts """
+
+        from datetime import timezone as tz
+
+        item = event.GetEventObject().FindItemById(event.GetId())
+        itemlabel = item.ItemLabel
+
+        if self.active_alaro is None:
+            self.active_alaro = IRM_Alaro()
+
+        if itemlabel == _("List runs on FTP server"):
+
+            runs = self.active_alaro.run_dates_str
+
+            for run in runs:
+                logging.info(f'Available Alaro run: {run}')
+
+            dlg = wx.MessageDialog(self, _('Available Alaro runs:\n\n') + '\n'.join(runs), _('Alaro runs'), wx.OK | wx.ICON_INFORMATION)
+            dlg.ShowModal()
+
+        elif itemlabel == _("Rain intensity [mm/h] - Comparison"):
+
+            data_dir = self.active_alaro.data_directory
+            dates = self.active_alaro.list_run_dates_cached()
+
+            dates_str = [datetime.strptime(date, '%Y%m%d%H').replace(tzinfo=tz.utc) for date in dates]
+            dates_str = [date.strftime('%Y-%m-%d %H+00') for date in dates_str]
+
+            dlg = wx.MultiChoiceDialog(self, _('Choose one or multiple Alaro runs to treat'), _('Choose Alaro runs'), choices=dates_str)
+            ret = dlg.ShowModal()
+            if ret == wx.ID_CANCEL:
+                dlg.Destroy()
+                return
+            choices = dlg.GetSelections()
+            dlg.Destroy()
+
+            if len(choices) == 0:
+                logging.warning(_('No Alaro run selected'))
+                return
+
+            dates = [dates[choice] for choice in choices]
+
+            dlg = wx.FileDialog(self, _('Choose output video filename'), wildcard='mp4 (*.mp4)|*.mp4', defaultDir=str(data_dir), defaultFile='Alaro_RainIntensity_Comparison.mp4', style=wx.FD_SAVE)
+            ret = dlg.ShowModal()
+            if ret == wx.ID_OK:
+                output_dir = dlg.GetPath()
+            dlg.Destroy()
+
+            self.active_alaro.reset_gdf()
+
+            pgbar = wx.ProgressDialog(_('Loading Alaro runs'), _('Loading data...'), maximum=100, parent=self, style=wx.PD_APP_MODAL | wx.PD_AUTO_HIDE)
+
+            pgbar.Pulse(_('Creating video...'))
+            outfile = self.active_alaro.video_gradient_cumulated_rain_compare(Path(output_dir), run_dates = dates)
+            pgbar.Update(100, _('Video created'))
+            pgbar.Destroy()
+
+            logging.info(f'Videos: {dates}')
+
+            # Propose to lauch the video
+            dlg = wx.MessageDialog(self, _('Video created: {}\n\nDo you want to launch it ?').format(outfile), _('Video created'), wx.YES_NO | wx.ICON_QUESTION)
+            ret = dlg.ShowModal()
+            if ret == wx.ID_YES:
+                if sys.platform == "win32":
+                    os.startfile(outfile)
+                # elif sys.platform == "darwin":
+                #     subprocess.call(["open", outfile])
+                # elif sys.platform == "linux":
+                #     subprocess.call(["xdg-open", outfile])
+            dlg.Destroy()
+
+        elif itemlabel == _("Total precip. [mm] - One run"):
+
+            data_dir = self.active_alaro.data_directory
+            dates = self.active_alaro.list_run_dates_cached()
+
+            dates_str = [datetime.strptime(date, '%Y%m%d%H').replace(tzinfo=tz.utc) for date in dates]
+            dates_str = [date.strftime('%Y-%m-%d %H+00') for date in dates_str]
+
+            dlg = wx.SingleChoiceDialog(self, _('Choose one Alaro run to treat'), _('Choose Alaro run'), choices=dates_str)
+            ret = dlg.ShowModal()
+            if ret == wx.ID_CANCEL:
+                dlg.Destroy()
+                return
+            choice = dlg.GetSelection()
+            dlg.Destroy()
+
+            if choice == wx.NOT_FOUND:
+                logging.warning(_('No Alaro run selected'))
+                return
+
+            dates = [dates[choice]]
+
+            dlg = wx.FileDialog(self, _('Choose output video filename'), wildcard='mp4 (*.mp4)|*.mp4', defaultDir=str(data_dir), defaultFile=f"Alaro_cumulated_rain_{dates[0]}.mp4", style=wx.FD_SAVE)
+            ret = dlg.ShowModal()
+            if ret == wx.ID_OK:
+                output_dir = dlg.GetPath()
+            dlg.Destroy()
+
+            self.active_alaro.reset_gdf()
+
+            pgbar = wx.ProgressDialog(_('Loading Alaro runs'), _('Loading data...'), maximum=100, parent=self, style=wx.PD_APP_MODAL | wx.PD_AUTO_HIDE)
+
+            pgbar.Pulse(_('Creating video...'))
+            outfile = self.active_alaro.video_cumulated_rain(dates[0], Path(output_dir))
+            pgbar.Update(100, _('Video created'))
+            pgbar.Destroy()
+
+            # Propose to lauch the video
+            dlg = wx.MessageDialog(self, _('Video created: {}\n\nDo you want to launch it ?').format(outfile), _('Video created'), wx.YES_NO | wx.ICON_QUESTION)
+            ret = dlg.ShowModal()
+            if ret == wx.ID_YES:
+                if sys.platform == "win32":
+                    os.startfile(outfile)
+                # elif sys.platform == "darwin":
+                #     subprocess.call(["open", outfile])
+                # elif sys.platform == "linux":
+                #     subprocess.call(["xdg-open", outfile])
+            dlg.Destroy()
+
+        elif itemlabel == _("Total precip. [mm] - Multiple runs"):
+
+            data_dir = self.active_alaro.data_directory
+            dates = self.active_alaro.list_run_dates_cached()
+
+            dates_str = [datetime.strptime(date, '%Y%m%d%H').replace(tzinfo=tz.utc) for date in dates]
+            dates_str = [date.strftime('%Y-%m-%d %H+00') for date in dates_str]
+
+            dlg = wx.MultiChoiceDialog(self, _('Choose multiple Alaro runs to treat'), _('Choose Alaro runs'), choices=dates_str)
+            ret = dlg.ShowModal()
+            if ret == wx.ID_CANCEL:
+                dlg.Destroy()
+                return
+            choices = dlg.GetSelections()
+            dlg.Destroy()
+
+            if len(choices) == 0:
+                logging.warning(_('No Alaro run selected'))
+                return
+
+            dates = [dates[choice] for choice in choices]
+
+            dlg = wx.DirDialog(self, _('Choose output video directory'), defaultPath=str(data_dir), style=wx.DD_DEFAULT_STYLE | wx.DD_DIR_MUST_EXIST)
+            ret = dlg.ShowModal()
+            if ret == wx.ID_OK:
+                output_dir = dlg.GetPath()
+            dlg.Destroy()
+
+            self.active_alaro.reset_gdf()
+
+            pgbar = wx.ProgressDialog(_('Loading Alaro runs'), _('Loading data...'), maximum=100, parent=self, style=wx.PD_APP_MODAL | wx.PD_AUTO_HIDE)
+
+            pgbar.Pulse(_('Creating video...'))
+            outdir = self.active_alaro.videos_cumulated_rain_allforecasts(output_dir, run_dates=dates)
+            pgbar.Update(100, _('Video created'))
+            pgbar.Destroy()
+
+        elif itemlabel == _("Download runs"):
+
+            pgbar = wx.ProgressDialog(_('Downloading Alaro runs'), _('Downloading data...'), maximum=100, parent=self, style=wx.PD_APP_MODAL | wx.PD_AUTO_HIDE)
+            pgbar.Pulse(_('Downloading Alaro runs...'))
+            ret = self.active_alaro.download_all_available_files()
+            pgbar.Update(100, _('Alaro runs downloaded'))
+            pgbar.Destroy()
+
+            for r in ret:
+                logging.info(f'Downloaded Alaro run: {r}')
+
+        elif itemlabel == _("Download runs - Only Rain and Temperature"):
+
+            pgbar = wx.ProgressDialog(_('Downloading Alaro runs'), _('Downloading data...'), maximum=100, parent=self, style=wx.PD_APP_MODAL | wx.PD_AUTO_HIDE)
+            pgbar.Pulse(_('Downloading Alaro runs...'))
+            ret = self.active_alaro.download_TotalPrecipitations_available_files()
+            pgbar.Update(100, _('Alaro runs downloaded'))
+            pgbar.Destroy()
+
+            for r in ret:
+                logging.info(f'Downloaded Alaro run: {r}')
+
+        elif itemlabel == _("Load run"):
+
+            data_dir = self.active_alaro.data_directory
+            dates = self.active_alaro.list_run_dates_cached()
+
+            dates_str = [datetime.strptime(date, '%Y%m%d%H').replace(tzinfo=tz.utc) for date in dates]
+            dates_str = [date.strftime('%Y-%m-%d %H+00') for date in dates_str]
+
+            dlg = wx.MultiChoiceDialog(self, _('Choose one or multiple Alaro runs to load'), _('Choose Alaro runs'), choices=dates_str)
+            ret = dlg.ShowModal()
+            if ret == wx.ID_CANCEL:
+                dlg.Destroy()
+                return
+            choices = dlg.GetSelections()
+            dlg.Destroy()
+
+            if len(choices) == 0:
+                logging.warning(_('No Alaro run selected'))
+                return
+
+            dates = [dates[choice] for choice in choices]
+
+            self.active_alaro.reset_gdf()
+
+            pgbar = wx.ProgressDialog(_('Loading Alaro runs'), _('Loading data...'), maximum=100, parent=self, style=wx.PD_APP_MODAL | wx.PD_AUTO_HIDE)
+            pgbar.Pulse(_('Loading Alaro runs...'))
+            self.active_alaro.load_grib_data_to_gdf(GribFiles.FILE_TotPrecip, dates)
+            pgbar.Update(100, _('Alaro runs loaded'))
+            pgbar.Destroy()
+            logging.info(f'Loaded Alaro run: {dates}')
+
+        elif itemlabel == _("Add forecast(s) as array(s)"):
+
+            if self.active_alaro._gdf is None:
+                logging.warning(_('No Alaro run loaded -- Please load a run first'))
+                return
+
+            columns = self.active_alaro.get_forecast_columns()
+
+            columns_str = [_convert_col2date_str(col) for col in columns]
+
+            dlg = wx.MultiChoiceDialog(self, _('Choose one or multiple Alaro forecast to show'), _('Choose Alaro forecasts'), choices=columns_str)
+            ret = dlg.ShowModal()
+            if ret == wx.ID_CANCEL:
+                dlg.Destroy()
+                return
+
+            choices = dlg.GetSelections()
+            dlg.Destroy()
+
+            if len(choices) == 0:
+                logging.warning(_('No Alaro forecast selected'))
+                return
+
+            columns = [columns[choice] for choice in choices]
+
+            pgbar = wx.ProgressDialog(_('Loading Alaro forecasts'), _('Loading data...'), maximum=100, parent=self, style=wx.PD_APP_MODAL | wx.PD_AUTO_HIDE)
+            pgbar.Pulse(_('Loading Alaro forecasts...'))
+            arrays = self.active_alaro.forecasts_to_arrays(columns)
+            pgbar.Update(100, _('Alaro forecasts loaded'))
+            pgbar.Destroy()
+
+            if len(arrays) == 0:
+                logging.warning(_('No Alaro forecast to show'))
+                return
+
+            for array in arrays:
+                array.array *= 1000.0  # Convert from m to mm
+                self.add_object('array', newobj=array, id=array.idx)
+
+            self.Refresh()
+
+            logging.info(_('Alaro forecasts added to the viewer'))
+
+        elif itemlabel == _("Show grid"):
+
+            if self.active_alaro._zones is None:
+                self.active_alaro._prepare_Zones_from_grib(GribFiles.FILE_TotPrecip, self.active_alaro.run_dates[0])
+
+            self.active_alaro._zones.prep_listogl()
+            self.add_object('vector', newobj=self.active_alaro._zones, id='Alaro_grid')
+            self.Refresh()
+
+        elif itemlabel == _("Add all forecasts as arrays"):
+            if self.active_alaro._gdf is None:
+                logging.warning(_('No Alaro run loaded -- Please load a run first'))
+                return
+
+            pgbar = wx.ProgressDialog(_('Loading Alaro forecasts'), _('Loading data...'), maximum=100, parent=self, style=wx.PD_APP_MODAL | wx.PD_AUTO_HIDE)
+            pgbar.Pulse(_('Loading Alaro forecasts...'))
+            arrays = self.active_alaro.forecasts_to_arrays()
+            pgbar.Update(100, _('Alaro forecasts loaded'))
+            pgbar.Destroy()
+
+            if len(arrays) == 0:
+                logging.warning(_('No Alaro forecast to show'))
+                return
+
+            for array in arrays:
+                array.array *= 1000.0  # Convert from m to mm
+                self.add_object('array', newobj=array, id=array.idx, ToCheck= False)
+
+            self.Refresh()
+            logging.info(_('Alaro forecasts added to the viewer'))
+
+        elif itemlabel == _("Plot for XY"):
+            if self.active_alaro._gdf is None:
+                logging.warning(_('No Alaro run loaded -- Please load a run first'))
+                return
+
+            self.start_action('plot alaro xy', _('Plot Alaro forecasts for a specific point -- Right click to select points'))
+
     def menu_walous(self):
 
         if self.menuwalous is None:
             self.menuwalous = wx.Menu()
 
-            self.menuwalous_crop = self.menuwalous.Append(wx.ID_ANY, _("Crop on active array"), _("Crop active array"))
-            self.menuwalous_cropscreen = self.menuwalous.Append(wx.ID_ANY, _("Crop on screen"), _("Crop screen"))
-            self.menuwalous_map = self.menuwalous.Append(wx.ID_ANY, _("Map active array"), _("Map active array"))
-            self.menuwalous_legend = self.menuwalous.Append(wx.ID_ANY, _("Legend"), _("Legend"))
-            self._walous_filepath = None
+            self._uts_menu = wx.Menu()
+            self._ocs_menu = wx.Menu()
+
+            self._uts_crop = self._uts_menu.Append(wx.ID_ANY, _("Crop on active array"), _("Crop active array"))
+            self._uts_cropscreen = self._uts_menu.Append(wx.ID_ANY, _("Crop on screen"), _("Crop screen"))
+            self._uts_map = self._uts_menu.Append(wx.ID_ANY, _("Map active array (WAL_UTS -> Manning)"), _("Map Walous UTS active array to Manning's n"))
+            self._uts_legend = self._uts_menu.Append(wx.ID_ANY, _("Legend"), _("Legend"))
+
+            self._ocs_crop = self._ocs_menu.Append(wx.ID_ANY, _("Crop on active array"), _("Crop active array"))
+            self._ocs_cropscreen = self._ocs_menu.Append(wx.ID_ANY, _("Crop on screen"), _("Crop screen"))
+            self._ocs_map = self._ocs_menu.Append(wx.ID_ANY, _("Map active array (WAL_OCS -> Hydrology)"), _("Map Walous OCS active array to Hydrology's landuse classification"))
+            self._ocs_map = self._ocs_menu.Append(wx.ID_ANY, _("Map active array (WAL_OCS -> Manning)"), _("Map Walous OCS active array to Manning's n"))
+            self._ocs_legend = self._ocs_menu.Append(wx.ID_ANY, _("Legend"), _("Legend"))
+
+            self._walous_UTS_filepath = None
+            self._walous_OCS_filepath = None
             self._walous_layer = None
             self._walous_map = None
 
+            self.menuwalous.AppendSubMenu(self._uts_menu, _('UTS'))
+            self.menuwalous.AppendSubMenu(self._ocs_menu, _('OCS'))
             self.menubar.Append(self.menuwalous, _('Walous'))
 
-            self.menuwalous.Bind(wx.EVT_MENU, self.Onmenuwalous)
+            self._uts_menu.Bind(wx.EVT_MENU, self.Onmenuwalous_uts)
+            self._ocs_menu.Bind(wx.EVT_MENU, self.Onmenuwalous_ocs)
 
     def menu_drowning(self):
 
@@ -3103,21 +3583,22 @@ class WolfMapViewer(wx.Frame):
 
             return [xmin, ymin, xmax, ymax]
 
-    def get_bounds(self):
+    def get_bounds(self, gridsize:float = None) -> tuple:
         """
         Retourne les limites de la zone d'affichage, voir aussi get_canvas_bounds
 
         :return: ([xmin, xmax], [ymin, ymax])
         """
-        xmin, ymin, xmax, ymax = self.get_canvas_bounds()
+        xmin, ymin, xmax, ymax = self.get_canvas_bounds(gridsize=gridsize)
+
         return ([xmin, xmax], [ymin, ymax])
 
-    def get_bounds_as_polygon(self) -> vector:
+    def get_bounds_as_polygon(self, gridsize:float = None) -> vector:
         """
         Retourne les limites de la zone d'affichage sous forme de polygone
         :return: vector
         """
-        xmin, ymin, xmax, ymax = self.get_canvas_bounds()
+        xmin, ymin, xmax, ymax = self.get_canvas_bounds(gridsize=gridsize)
         poly = vector()
         poly.add_vertex(wolfvertex(xmin, ymin))
         poly.add_vertex(wolfvertex(xmax, ymin))
@@ -3126,8 +3607,169 @@ class WolfMapViewer(wx.Frame):
         poly.force_to_close()
         return poly
 
-    def Onmenuwalous(self, event: wx.MenuEvent):
+    def Onmenuwalous_ocs(self, event: wx.MenuEvent):
+        """ Handle the Walous OCS menu events """
+        id = event.GetId()
+        item = self.menubar.FindItemById(event.GetId())
+
+        if item is None:
+            return
+
+        itemlabel = item.ItemLabel
+
+        if itemlabel in [_("Crop on active array"), _("Crop on screen")]:
+
+            if itemlabel == _("Crop on screen"):
+
+                bounds = self.get_canvas_bounds(gridsize=1.)
+
+                def_outdrir = ''
+                spatial_res = 1.
+
+                if self.active_array is not None:
+                    spatial_res = self.active_array.dx
+
+                dlg = wx.TextEntryDialog(None,_("Spatial resolution [m] ?"), value = str(spatial_res))
+
+                dlg.ShowModal()
+                try:
+                    spatial_res = float(dlg.GetValue())
+                    dlg.Destroy()
+                except:
+                    dlg.Destroy()
+                    logging.warning(_("Bad value -- Rety"))
+                    return
+
+            else:
+
+                if self.active_array is None:
+                    logging.warning(_('No active array -- Please activate data first'))
+                    return
+
+                bounds = self.active_array.get_bounds()
+                def_outdrir = Path(self.active_array.filename).parent
+                spatial_res = self.active_array.dx
+
+            from .pywalous import update_palette_walous_ocs
+
+            if self._walous_OCS_filepath is None:
+                dlg = wx.FileDialog(self, _("Choose the Walous OCS Tif file"), wildcard="Tif file (*.tif)|*.tif|all (*.*)|*.*", style=wx.FD_OPEN)
+                if dlg.ShowModal() == wx.ID_CANCEL:
+                    dlg.Destroy()
+                    return
+
+                self._walous_OCS_filepath = Path(dlg.GetPath())
+                dlg.Destroy()
+
+
+            dlg = wx.FileDialog(self, _("Choose the output file"), wildcard="Geotif (*.tif)|*.tif|all (*.*)|*.*", style=wx.FD_SAVE, defaultDir=str(def_outdrir))
+            if dlg.ShowModal() == wx.ID_CANCEL:
+                dlg.Destroy()
+                return
+
+            output = Path(dlg.GetPath())
+            dlg.Destroy()
+
+            locwalous = WolfArray(fname=self._walous_OCS_filepath,
+                                  crop = bounds)
+
+            if Path(output).exists():
+                logging.info(_('File {} created').format(output))
+            else:
+                logging.error(_('File {} not created').format(output))
+                return
+
+            dlg = wx.MessageDialog(self, _('Do you want to load the created file ?'), _('Load file'), wx.YES_NO | wx.ICON_QUESTION)
+            ret = dlg.ShowModal()
+
+            if ret == wx.ID_CANCEL:
+                dlg.Destroy()
+                return
+
+            elif ret == wx.ID_YES:
+                walousarray = WolfArray(fname=output)
+                update_palette_walous_ocs(walousarray.mypal)
+                self.add_object('array', newobj=walousarray, id = 'walous_ocs_crop')
+                dlg.Destroy()
+
+        elif itemlabel == _("Legend"):
+            from .pywalous import Walous_OCS_Legend
+            newlegend = Walous_OCS_Legend(self)
+            newlegend.Show()
+
+        elif itemlabel == _("Map active array (WAL_OCS -> Hydrology)"):
+
+            from .pywalous import DlgMapWalous2Hydrology
+
+            if self.active_array is None:
+                logging.warning(_('No active array -- Please activate data first'))
+                return
+
+            vals = self.active_array.get_unique_values()
+
+            if self._walous_layer is None:
+
+                if vals[0] > 11:
+                    logging.error(_('You have values greater than 11 -- Please check your data'))
+                    return
+
+            dlg = DlgMapWalous2Hydrology(self,)
+
+            ret = dlg.ShowModal()
+
+            if ret == wx.ID_CANCEL:
+                dlg.Destroy()
+                return
+
+            mapvals = dlg.get_mapping()
+            dlg.Destroy()
+
+            if mapvals == -1:
+                logging.error(_('Bad values -- retry'))
+                return
+
+            self.active_array.map_values(mapvals)
+
+            self.active_array.reset_plot()
+
+        elif itemlabel == _("Map active array (WAL_OCS -> Manning)"):
+
+            from .pywalous import DlgMapWalousOCS2Manning
+
+            if self.active_array is None:
+                logging.warning(_('No active array -- Please activate data first'))
+                return
+
+            vals = self.active_array.get_unique_values()
+
+            if self._walous_layer is None:
+
+                if vals[0] > 11:
+                    logging.error(_('You have values greater than 11 -- Please check your data'))
+                    return
+
+            dlg = DlgMapWalousOCS2Manning(self,)
+
+            ret = dlg.ShowModal()
+
+            if ret == wx.ID_CANCEL:
+                dlg.Destroy()
+                return
+
+            mapvals = dlg.get_mapping()
+            dlg.Destroy()
+
+            if mapvals == -1:
+                logging.error(_('Bad values -- retry'))
+                return
+
+            self.active_array.map_values(mapvals)
 
+            self.active_array.reset_plot()
+
+
+    def Onmenuwalous_uts(self, event: wx.MenuEvent):
+        """ Handle the Walous UTS menu events """
         id = event.GetId()
         item = self.menubar.FindItemById(event.GetId())
 
@@ -3170,15 +3812,15 @@ class WolfMapViewer(wx.Frame):
                 def_outdrir = Path(self.active_array.filename).parent
                 spatial_res = self.active_array.dx
 
-            from .pywalous import Walous_data, WALOUS2MANNING_MAJ_NIV1, WALOUS2MANNING_MAJ_NIV2, update_palette_walous
+            from .pywalous import Walous_data, WALOUS2MANNING_MAJ_NIV1, WALOUS2MANNING_MAJ_NIV2, update_palette_walous_uts
 
-            if self._walous_filepath is None:
+            if self._walous_UTS_filepath is None:
                 dlg = wx.FileDialog(self, _("Choose the Walous shape file"), wildcard="Shapefile (*.shp)|*.shp|all (*.*)|*.*", style=wx.FD_OPEN)
                 if dlg.ShowModal() == wx.ID_CANCEL:
                     dlg.Destroy()
                     return
 
-                self._walous_filepath = Path(dlg.GetPath())
+                self._walous_UTS_filepath = Path(dlg.GetPath())
                 dlg.Destroy()
 
 
@@ -3199,7 +3841,7 @@ class WolfMapViewer(wx.Frame):
 
             self._walous_layer = dlg.GetStringSelection()
 
-            locwalous = Walous_data(self._walous_filepath.parent, self._walous_filepath.name)
+            locwalous = Walous_data(self._walous_UTS_filepath.parent, self._walous_UTS_filepath.name)
             ret = locwalous.rasterize(bounds=bounds,
                                       layer=self._walous_layer,
                                       fn_out=output,
@@ -3224,20 +3866,20 @@ class WolfMapViewer(wx.Frame):
 
             elif ret == wx.ID_YES:
                 walousarray = WolfArray(fname=output)
-                update_palette_walous(self._walous_layer, walousarray.mypal)
-                self.add_object('array', newobj=walousarray, id = 'walous_crop')
+                update_palette_walous_uts(self._walous_layer, walousarray.mypal)
+                self.add_object('array', newobj=walousarray, id = 'walous_uts_crop')
                 dlg.Destroy()
 
         elif itemlabel == _("Legend"):
 
-            from .pywalous import WalousLegend
+            from .pywalous import Walous_UTS_Legend
 
-            newlegend = WalousLegend(self)
+            newlegend = Walous_UTS_Legend(self)
             newlegend.Show()
 
-        elif itemlabel == _("Map active array"):
+        elif itemlabel == _("Map active array (WAL_UTS -> Manning)"):
 
-            from .pywalous import DlgMapWalous, WALOUS2MANNING_MAJ_NIV1, WALOUS2MANNING_MAJ_NIV2
+            from .pywalous import DlgMapWalous2Manning
 
             if self.active_array is None:
                 logging.warning(_('No active array -- Please activate data first'))
@@ -3252,7 +3894,7 @@ class WolfMapViewer(wx.Frame):
                 else:
                     self._walous_layer = 'MAJ_NIV1'
 
-            dlg = DlgMapWalous(self, which=self._walous_layer)
+            dlg = DlgMapWalous2Manning(self, which=self._walous_layer)
 
             ret = dlg.ShowModal()
 
@@ -4888,6 +5530,10 @@ class WolfMapViewer(wx.Frame):
 
                         if nb_ticks_x > 10 or nb_ticks_y > 10:
                             logging.error(_('Too many ticks for the image. Please raise the ticks size in the global options.'))
+                            dlg = wx.MessageDialog(None, _('Too many ticks for the image. Please raise the ticks size in the global options.'), _('Error'), style=wx.OK)
+                            dlg.ShowModal()
+                            dlg.Destroy()
+                            wx.TheClipboard.Close()
                             return
 
                     # Création d'un graphique Matplotlib
@@ -4964,10 +5610,25 @@ class WolfMapViewer(wx.Frame):
                         if palette.values is not None:
                             bufpal = io.BytesIO()
                             palette.export_image(bufpal,'v')
-                            bufpal.seek(0)
 
-                            #lecture du buffer et conversion en image avec PIL
-                            impal = Image.open(bufpal)
+                            try:
+                                bufpal.seek(0)
+
+                                #lecture du buffer et conversion en image avec PIL
+                                impal = Image.open(bufpal)
+
+                            except Exception as e:
+
+                                text = _('Error while creating the colormap/palette image !')
+                                text += '\n'
+                                text += _('Please check if an array or a 2D result is active !')
+
+                                logging.error(text)
+                                dlg = wx.MessageDialog(None, text, _('Error'), style=wx.OK)
+                                dlg.ShowModal()
+                                dlg.Destroy()
+                                return
+
                             impal = impal.resize((int(impal.size[0]*im.size[1]*.8/impal.size[1]),int(im.size[1]*.8)))
 
                             imnew = Image.new('RGB',(im.size[0]+impal.size[0], im.size[1]), (255,255,255))
@@ -5268,6 +5929,8 @@ class WolfMapViewer(wx.Frame):
         self.menu_1to9 =self.menu_options.Append(wx.ID_ANY, _('Colors for selections 1->9'), _('Selections'))
         self.Bind(wx.EVT_MENU, self.colors1to9.change_colors, self.menu_1to9)
 
+        self.menu_qdfidf()
+
         self.Show(True)
 
     def OnChangeTitle(self, e):
@@ -5356,7 +6019,7 @@ class WolfMapViewer(wx.Frame):
 
     def setbounds(self, updatescale=True):
         """
-        Calcule les limites visibles de la fenêtrte graphique sur base des
+        Calcule les limites visibles de la fenêtre graphique sur base des
         facteurs d'échelle courants
         """
 
@@ -5393,12 +6056,14 @@ class WolfMapViewer(wx.Frame):
             self.ymin = self.mousey - self.height / 2.
             self.ymax = self.ymin + self.height
 
+
         self.mybackisloaded = False
         self.myfrontisloaded = False
 
         self.Refresh()
         self.mimicme()
 
+
     def setsizecanvas(self,width,height):
         """ Redimensionne la fenêtre graphique """
         self.canvas.SetClientSize(width, height)
@@ -5560,41 +6225,48 @@ class WolfMapViewer(wx.Frame):
         'Wet_Bulb_Poten_Temper',
         'freezing_level_zeroDegC_isotherm'],
         """
-        alaro = {'ALARO': {'10m_u_wind_component': '10_m_u__wind_component',
-                          '10m_v_wind_component': '10_m_v__wind_component',
-                          '2m_Max_temp_since_ppp': '2_m_Max_temp_since_ppp',
-                          '2m_Min_temp_since_ppp': '2_m_Min_temp_since_ppp',
-                          '2m_dewpoint_temperature': '2_m_dewpoint_temperature',
-                          '2m_temperature': '2_m_temperature',
-                          '2m_Relative_humidity': '2m_Relative_humidity',
+        # alaro = {'ALARO': {'10m_u_wind_component': '10_m_u__wind_component',
+        #                   '10m_v_wind_component': '10_m_v__wind_component',
+        #                   '2m_Max_temp_since_ppp': '2_m_Max_temp_since_ppp',
+        #                   '2m_Min_temp_since_ppp': '2_m_Min_temp_since_ppp',
+        #                   '2m_dewpoint_temperature': '2_m_dewpoint_temperature',
+        #                   '2m_temperature': '2_m_temperature',
+        #                   '2m_Relative_humidity': '2m_Relative_humidity',
+        #                   'Convective_rain': 'Convective_rain',
+        #                   'Convective_snow': 'Convective_snow',
+        #                   'Geopotential': 'Geopotential',
+        #                   'Inst_flx_Conv_Cld_Cover': 'Inst_flx_Conv_Cld_Cover',
+        #                   'Inst_flx_High_Cld_Cover': 'Inst_flx_High_Cld_Cover',
+        #                   'Inst_flx_Low_Cld_Cover': 'Inst_flx_Low_Cld_Cover',
+        #                   'Inst_flx_Medium_Cld_Cover': 'Inst_flx_Medium_Cld_Cover',
+        #                   'Inst_flx_Tot_Cld_cover': 'Inst_flx_Tot_Cld_cover',
+        #                   'Large_scale_rain': 'Large_scale_rain',
+        #                   'Large_scale_snow': 'Large_scale_snow',
+        #                   'Mean_sea_level_pressure': 'Mean_sea_level_pressure',
+        #                   'Relative_humidity': 'Relative_humidity',
+        #                   'Relative_humidity_isobaric': 'Relative_humidity_isobaric',
+        #                   'SBL_Meridian_gust': 'SBL_Meridian_gust',
+        #                   'SBL_Zonal_gust': 'SBL_Zonal_gust',
+        #                   'Specific_humidity': 'Specific_humidity',
+        #                   'Surf_Solar_radiation': 'Surf_Solar_radiation',
+        #                   'Surf_Thermal_radiation': 'Surf_Thermal_radiation',
+        #                   'Surface_CAPE': 'Surface_CAPE',
+        #                   'Surface_Temperature': 'Surface_Temperature',
+        #                   'Surface_orography': 'Surface_orography',
+        #                   'Temperature': 'Temperature',
+        #                   'Total_precipitation': 'Total_precipitation',
+        #                   'U-velocity': 'U-velocity',
+        #                   'V-velocity': 'V-velocity',
+        #                   'Vertical_velocity': 'Vertical_velocity',
+        #                   'Wet_Bulb_Poten_Temper': 'Wet_Bulb_Poten_Temper',
+        #                   'freezing_level_zeroDegC_isotherm': 'freezing_level_zeroDegC_isotherm',}}
+        alaro = {'ALARO': {'2m_temperature': '2_m_temperature',
                           'Convective_rain': 'Convective_rain',
                           'Convective_snow': 'Convective_snow',
-                          'Geopotential': 'Geopotential',
-                          'Inst_flx_Conv_Cld_Cover': 'Inst_flx_Conv_Cld_Cover',
-                          'Inst_flx_High_Cld_Cover': 'Inst_flx_High_Cld_Cover',
-                          'Inst_flx_Low_Cld_Cover': 'Inst_flx_Low_Cld_Cover',
-                          'Inst_flx_Medium_Cld_Cover': 'Inst_flx_Medium_Cld_Cover',
-                          'Inst_flx_Tot_Cld_cover': 'Inst_flx_Tot_Cld_cover',
                           'Large_scale_rain': 'Large_scale_rain',
                           'Large_scale_snow': 'Large_scale_snow',
-                          'Mean_sea_level_pressure': 'Mean_sea_level_pressure',
-                          'Relative_humidity': 'Relative_humidity',
-                          'Relative_humidity_isobaric': 'Relative_humidity_isobaric',
-                          'SBL_Meridian_gust': 'SBL_Meridian_gust',
-                          'SBL_Zonal_gust': 'SBL_Zonal_gust',
-                          'Specific_humidity': 'Specific_humidity',
-                          'Surf_Solar_radiation': 'Surf_Solar_radiation',
-                          'Surf_Thermal_radiation': 'Surf_Thermal_radiation',
-                          'Surface_CAPE': 'Surface_CAPE',
                           'Surface_Temperature': 'Surface_Temperature',
-                          'Surface_orography': 'Surface_orography',
-                          'Temperature': 'Temperature',
-                          'Total_precipitation': 'Total_precipitation',
-                          'U-velocity': 'U-velocity',
-                          'V-velocity': 'V-velocity',
-                          'Vertical_velocity': 'Vertical_velocity',
-                          'Wet_Bulb_Poten_Temper': 'Wet_Bulb_Poten_Temper',
-                          'freezing_level_zeroDegC_isotherm': 'freezing_level_zeroDegC_isotherm',}}
+                          'Total_precipitation': 'Total_precipitation',}}
 
 
         for idx, (k, item) in enumerate(orthos.items()):
@@ -9281,7 +9953,7 @@ class WolfMapViewer(wx.Frame):
 
             self.create_2D_GPU_model()
 
-        elif itemlabel == _('Open Hydrological model'):
+        elif itemlabel == _('Create/Open Hydrological model'):
 
             self.open_hydrological_model()
 
@@ -9474,6 +10146,13 @@ class WolfMapViewer(wx.Frame):
             newmanager.mapviewer = self
             newmanager.Show()
 
+        elif itemlabel == _('Create INBE manager...'):
+
+            from .insyde_be.INBE_gui import INBEGui
+            newmanager = INBEGui()
+            newmanager.mapviewer = self
+            newmanager.Show()
+
         elif itemlabel==_('Create BC manager Wolf2D...'):
 
             if self.active_array is not None:
@@ -10486,7 +11165,7 @@ class WolfMapViewer(wx.Frame):
 
             if newobj is None:
 
-                if filename.endswith('.npz'):
+                if str(filename).endswith('.npz'):
 
                     wait = wx.BusyCursor()
                     logging.info(_('Start of importing arrays from npz file'))
@@ -11197,7 +11876,7 @@ class WolfMapViewer(wx.Frame):
 
         # curdict[id.lower()] = newobj
         if filename != '':
-            newobj._filename_vector = filename.lower() # FIXME useful ??
+            newobj._filename_vector = Path(filename).name.lower() # FIXME useful ??
         newobj.checked = ToCheck
 
         if isinstance(newobj,crosssections):
@@ -11210,24 +11889,113 @@ class WolfMapViewer(wx.Frame):
 
         return 0
 
+    def replace_object(self, id: str, newobj, drawing_type: draw_type = None):
+        """ Replace an object in the list of objects of type drawing_type """
+
+        if drawing_type is None:
+            for curdict in draw_type:
+                keys = self.get_list_keys(curdict, checked_state=None)
+                if id.lower() in keys:
+                    # The object exists in the current dictionary
+                    obj = self.get_obj_from_id(id, drawing_type=curdict)
+                    obj.reset_listogl()
+                    # Searching the object in all lists
+                    if obj is not None:
+                        curlist = self._get_list(drawing_type=curdict)
+                        if obj in curlist:
+                            pos = curlist.index(obj)
+                            if isinstance(newobj, curlist[pos].__class__):
+                                # Updating the tree item
+                                self.treelist.SetItemData(self.get_treeitem_from_obj(obj), newobj)
+                                curlist[pos] = newobj
+                                newobj.idx = id.lower()
+                            else:
+                                logging.error(f'Cannot replace {id} with {newobj.idx} - Different type of object')
+                        else:
+                            logging.error(f'Object {id} not found in list')
+                    else:
+                        logging.error(f'Object {id} not found in list')
+                else:
+                    logging.error(f'Object {id} not found in dictionary {curdict}')
+        else:
+            keys = self.get_list_keys(drawing_type, checked_state=None)
+            if id.lower() in keys:
+                # The object exists in the current dictionary
+                obj = self.get_obj_from_id(id, drawing_type=drawing_type)
+                obj.reset_listogl()
+                # Searching the object in all lists
+                if obj is not None:
+                    curlist = self._get_list(drawing_type=drawing_type)
+                    if obj in curlist:
+                        pos = curlist.index(obj)
+                        if isinstance(newobj, curlist[pos].__class__):
+                            # Updating the tree item
+                            self.treelist.SetItemData(self.get_treeitem_from_obj(obj), newobj)
+                            curlist[pos] = newobj
+                            newobj.idx = id.lower()
+                        else:
+                            logging.error(f'Cannot replace {id} with {newobj.idx} - Different type of object')
+                    else:
+                        logging.error(f'Object {id} not found in list')
+                else:
+                    logging.error(f'Object {id} not found in list')
+            else:
+                logging.error(f'Object {id} not found in dictionary {drawing_type}')
+
+        obj = self.get_obj_from_id(id, drawing_type=drawing_type)
+        obj_from_tree = self.get_obj_from_treeitem(self.get_treeitem_from_id(id, drawing_type=drawing_type))
+        if obj is not None and obj_from_tree is not None:
+            if obj is obj_from_tree:
+                logging.debug(f'Object {id} replaced successfully in the list and tree item')
+            else:
+                logging.error(f'Object {id} replaced in the list but not in the tree item - {obj} != {obj_from_tree}')
+        else:
+            logging.error(f'Object {id} not found in the list or tree item after replacement')
+
     def get_obj_from_treeitem(self, treeitem):
         """ Find the object associated with treeitem """
 
         return self.treelist.GetItemData(treeitem)
 
-    def getobj_from_id(self, id: str):
+
+    def get_treeitem_from_id(self, id: str, drawing_type: draw_type = None):
+        """ Find the tree item associated with id """
+
+        obj = self.get_obj_from_id(id, drawing_type=drawing_type)
+        if obj is not None:
+            return self.get_treeitem_from_obj(obj)
+        return None
+
+    def get_treeitem_from_obj(self, obj):
+        """ Find the tree item associated with obj.
+
+        Alias for "gettreeitem".
+        """
+
+        return self.gettreeitem(obj)
+
+    def getobj_from_id(self, id: str, drawing_type: draw_type = None):
         """ Find the object associated with id """
 
-        for curdict in draw_type:
-            keys = self.get_list_keys(curdict, checked_state=None)
+        if drawing_type is None:
+            for curdict in draw_type:
+                keys = self.get_list_keys(curdict, checked_state=None)
+                if id.lower() in keys:
+                    try:
+                        idx = keys.index(id.lower())
+                        return self.get_list_objects(curdict, checked_state=None)[idx]
+                    except:
+                        return None
+        else:
+            keys = self.get_list_keys(drawing_type, checked_state=None)
             if id.lower() in keys:
                 try:
                     idx = keys.index(id.lower())
-                    return self.get_list_objects(curdict, checked_state=None)[idx]
+                    return self.get_list_objects(drawing_type, checked_state=None)[idx]
                 except:
                     return None
 
-    def get_obj_from_id(self, id: str, drawtype: draw_type):
+    def get_obj_from_id(self, id: str, drawing_type: draw_type = None):
         """ Find the object associated with id in a specifid drawtype
 
         If you want to search in all drawtypes, use getobj_from_id instead.
@@ -11237,11 +12005,11 @@ class WolfMapViewer(wx.Frame):
 
         """
 
-        keys = self.get_list_keys(drawtype, checked_state=None)
+        keys = self.get_list_keys(drawing_type, checked_state=None)
         if id.lower() in keys:
             try:
                 idx = keys.index(id.lower())
-                return self.get_list_objects(drawtype, checked_state=None)[idx]
+                return self.get_list_objects(drawing_type, checked_state=None)[idx]
             except:
                 return None
 
@@ -11304,7 +12072,9 @@ class WolfMapViewer(wx.Frame):
             return None
 
     def get_list_keys(self, drawing_type:draw_type = None, checked_state:bool=True):
-        """ Create a list of keys of type draw_type
+        """ Create a list of keys of type draw_type.
+
+        Return a list of keys (idx) in LOWER CASE of objects of type draw_type.
 
         :param drawing_type: type of object to search - If None, return all objects
         :param checked_state: if True/False, return only keys of objects that are plotted or not. None return all objects.
@@ -11321,7 +12091,9 @@ class WolfMapViewer(wx.Frame):
         return self.get_list_keys(drawing_type, checked_state)
 
     def get_list_objects(self, drawing_type:draw_type = None, checked_state:bool=True):
-        """ Create a list of objects of type draw_type
+        """ Create a list of objects of type draw_type.
+
+        Return a list of keys (idx) in LOWER CASE of objects of type draw_type.
 
         :param drawing_type: type of object to search -- If None, return all objects.
         :param checked_state: if True/False, return only objects that are plotted or not. None return all objects.
@@ -11444,7 +12216,7 @@ class WolfMapViewer(wx.Frame):
                 self.treelist.CheckItem(self.selected_treeitem, False)
                 myobj.uncheck_plot()
 
-    def removeobj_from_id(self, id:str):
+    def removeobj_from_id(self, id:str, draw_type:draw_type = None):
         """ Remove object from id """
 
         myobj = self.getobj_from_id(id)
@@ -11646,14 +12418,26 @@ class WolfMapViewer(wx.Frame):
 
         elif text == _('Rename'):
             #Modification du nom de l'objet sélectionné
+
             if self.selected_object is not None:
                 #récupération de l'id courant
+                all_ids = self.get_list_ids(checked_state=None)
+
                 label = self.selected_object.idx
+
+                all_ids.remove(label.lower())
+
                 dlg = wx.TextEntryDialog(self, message=_('Chose a new label :'), value=label)
                 ret=dlg.ShowModal()
 
                 if ret == wx.ID_OK:
                     newlab = dlg.GetValue()
+
+                    if newlab.lower() in all_ids:
+                        wx.MessageBox(_('This label already exists. Please choose another one.'), _('Error'), wx.OK | wx.ICON_ERROR)
+                        dlg.Destroy()
+                        return
+
                     #MAJ de l'id dans l'objet
                     self.selected_object.idx = newlab
                     #MAJ de l'arbre
@@ -11806,6 +12590,15 @@ class WolfMapViewer(wx.Frame):
                 elif isinstance(self.selected_object, Drowning_victim_Viewer):
                     self.selected_object.saveas()
 
+                elif isinstance(self.selected_object, crosssections):
+                    filterArray = "vecz (*.vecz)|*.vecz|SXY (*.sxy)|*.sxy"
+                    fdlg = wx.FileDialog(self, "Choose file name for Cross Sections :" + self.selected_object.idx, wildcard=filterArray,
+                                         style=wx.FD_SAVE)
+                    ret = fdlg.ShowModal()
+                    if ret == wx.ID_OK:
+                        self.selected_object.saveas(fdlg.GetPath())
+                    fdlg.Destroy()
+
         elif text == _('Properties'):
 
             myobj = self.selected_object
@@ -12426,8 +13219,9 @@ class WolfMapViewer(wx.Frame):
         """ Show images of the checked alaro layers"""
         objs = self.get_list_objects(drawing_type=draw_type.WMSFORE, checked_state=True)
         for obj in objs:
-            img = Image.open(get_Alaro_legend(obj.idx.replace('alaro ', '')))
-            img.show()
+            if obj.category == 'ALARO':
+                img = Image.open(get_Alaro_legend(obj.idx.replace('alaro ', '')))
+                img.show()
 
     def getXY(self, pospix):
 
@@ -12536,6 +13330,15 @@ class WolfMapViewer(wx.Frame):
 
             self.rightdown = (x, y)
 
+        elif self.action == 'plot alaro xy':
+
+            if self.active_alaro._gdf is None:
+                logging.warning(_('No Alaro run loaded -- Please load a run first'))
+                return
+
+            fig = self.active_alaro.plot_Rain_and_TotPrecip4XY(x, y)
+            fig.show()
+
         elif self.action == 'distance along vector':
 
             # add a vertex to the vector
@@ -12602,6 +13405,14 @@ class WolfMapViewer(wx.Frame):
                 self.active_landmap.load_texture(x,y, which='low')
             self.Refresh()
 
+        elif self.action == 'pick municipality':
+            # Pick a municipality if loaded
+            if self.active_qdfidf is None:
+                logging.warning(_('No municipality data available -- Please activate the data and retry !'))
+                return
+
+            self.active_qdfidf.pick_municipality(x, y, self.get_canvas_bounds())
+
         elif self.action == 'pick a picture':
             # Pick a picture
 
@@ -15454,6 +16265,13 @@ class WolfMapViewer(wx.Frame):
             # Dessin des images
             self._plotting(draw_type.PICTURECOLLECTION)
 
+            # Dessin des QDF/IDF
+            if self.active_qdfidf is not None:
+                self.active_qdfidf.plot(sx = self.sx, sy=self.sy,
+                                        xmin=self.xmin, ymin=self.ymin,
+                                        xmax=self.xmax, ymax=self.ymax,
+                                        size = (self.xmax - self.xmin) / 100.)
+
             # Gestion des BC (si actif)
             if self.active_bc is not None:
                 self.active_bc.plot()