wolfhece 2.2.41__py3-none-any.whl → 2.2.43__py3-none-any.whl

wolfhece/pywalous.py

@@ -22,6 +22,125 @@ import wx.grid
 
 from .PyTranslate import _
 from .PyPalette import wolfpalette
+from .wolf_array import WolfArray, header_wolf
+
+
+"""
+Hydrology classification for Land Use in Wallonia:
+1 = forêt
+2 = prairie
+3 = culture
+4 = pavés/urbain
+5 = rivière
+6 = plan d'eau
+"""
+
+HYDROLOGY_LANDUSE_FR = {"forêt" : 1.,
+                        "prairie" : 2.,
+                        "culture" : 3.,
+                        "pavés/urbain" : 4.,
+                        "rivière" : 5.,
+                        "plan d'eau" : 6.,
+                     }
+
+HYDROLOGY_LANDUSE_EN = {"forest" : 1.,
+                        "meadow" : 2.,
+                        "crop" : 3.,
+                        "paved/urban" : 4.,
+                        "river" : 5.,
+                        "water body" : 6.,
+                        }
+
+# WALOUS OCS
+# ----------
+
+"""
+Extrait de https://metawal.wallonie.be/geonetwork/srv/api/records/86462606-4a21-49c8-ab9e-564ccba681b7/attachments/DescriptionLegende_WALOUS_OCS.pdf
+
+1. Revêtement artificiel du sol: cette classe reprend tous les revêtements artificiels du sol
+de nature peu ou pas perméable (ex béton, le bitume ou les pavés). Ceci comprend le réseau routier, les trottoirs, les terrasses, les parkings (les emplacements de parking semiperméables non végétalisés sont également dans cette catégorie) et les terrains ou pistes de sport en matériaux synthétiques.
+2. Constructions artificielles hors sol : cette classe reprend tous les bâtiments et autres constructions s'élevant au-dessus du sol.
+3. Réseau ferroviaire : cette classe reprend les rails et ballasts des chemins de fer encore en activité (les RaVEL entre dans la classe 1). Cette classe se différencie de la classe 1 par sa grande perméabilité.
+4. Sols nus : cette classe inclut tout type de roche mère n'étant couverte par des végétaux supérieurs à aucun moment de l'année. Ces sols sont soit naturels (roches affleurantes, falaises, berges caillouteuses…), soit générés par l'activité humaine (extraction, sols compactés, coupe à blanc de l'année…)
+5. Eaux de surface : cette classe comprend toutes les surfaces d'eau libre, naturelles ou artificielles. Ceci inclut donc à la fois les cours d'eau (rivières, fleuves et canaux) et les plans d'eau (mares, étangs, lacs, bassins de décantation, piscines extérieures).
+6. Couvert herbacé en rotation : cette classe reprend les parcelles combinant un couvert herbacé une partie de l'année et un sol nu temporairement mis à nu. On y retrouve toutes les cultures annuelles, ainsi que les prairies temporaires succédant à une culture annuelle.
+7. Couvert herbacé continu : cette classe reprend tous les sols recouverts par de la végétation herbacée tout au long de l'année. Cette végétation peut être d'origine naturelle (landes, mégaphorbiaies, tourbières, pelouses naturelles, végétation rudérale recolonisant une friche ou une ancienne coupe à blanc…), agricole (prés et prairies non labourés dans la saison) ou artificielle (jardin, terrains de sport, parcs…)
+8. Arbres résineux : cette classe comprend tous les arbres de plus de 3 m (isolés, en haie ou en peuplement) du groupe des résineux (gymnosperme).
+9. Arbres feuillus : cette classe comprend tous les arbres de plus de 3 m (isolés, en haie ou en peuplement) du groupe des feuillus (angiosperme).
+80. Arbustes résineux : cette classe comprend tous les arbres et arbustes de moins de 3 m (isolés, en haie ou en peuplement) du groupe des résineux (gymnosperme).
+90. Arbustes feuillus : cette classe comprend tous les arbres et arbuste de moins de 3 m (isolés, en haie ou en peuplement) du groupe des feuillus (angiosperme).
+
+
+<paletteEntry color="#8a8a8a" value="1" alpha="255" label="Revêtement artificiel au sol"/>
+<paletteEntry color="#dc0f0f" value="2" alpha="255" label="Constructions artificielles hors sol"/>
+<paletteEntry color="#4e4e4e" value="3" alpha="255" label="Réseau ferroviaire"/>
+<paletteEntry color="#d0d0d0" value="4" alpha="255" label="Sols nus"/>
+<paletteEntry color="#2461f7" value="5" alpha="255" label="Eaux de surface"/>
+<paletteEntry color="#ffff73" value="6" alpha="255" label="Couvert herbacé en rotation dans l'année (ex: culture annuelle)"/>
+<paletteEntry color="#e9ffbe" value="7" alpha="255" label="Couvert herbacé toute l'année"/>
+<paletteEntry color="#003200" value="8" alpha="255" label="Résineux (> 3m)"/>
+<paletteEntry color="#007800" value="80" alpha="255" label="Résineux (≤ 3m)"/>
+<paletteEntry color="#28c828" value="9" alpha="255" label="Feuillus (> 3m)"/>
+<paletteEntry color="#b7e8b0" value="90" alpha="255" label="Feuillus (≤ 3m)"/>
+<paletteEntry color="#e5ea3f" value="0" alpha="0" label="Pas de données"/>
+"""
+
+WALOUS_OCS = { "Revêtement artificiel au sol": 1.,
+              "Constructions artificielles hors sol": 2.,
+              "Réseau ferroviaire": 3.,
+              "Sols nus": 4.,
+              "Eaux de surface": 5.,
+              "Couvert herbacé en rotation dans l'année": 6.,
+              "Couvert herbacé toute l'année": 7.,
+              "Arbres résineux (>= 3m)": 8.,
+              "Arbres feuillus (>= 3m)": 9.,
+              "Arbustes résineux (< 3m)": 80.,
+              "Arbustes feuillus (< 3m)": 90.}
+
+WALOUS_OCS_COLORMAP = {
+    0.: (229, 234, 63, 0),       # #e5ea3f Pas de données
+    1.: (138, 138, 138, 255),    # #8a8a8a Revêtement artificiel au sol
+    2.: (220, 15, 15, 255),      # #dc0f0f Constructions artificielles hors sol
+    3.: (78, 78, 78, 255),       # #4e4e4e Réseau ferroviaire
+    4.: (208, 208, 208, 255),    # #d0d0d0 Sols nus
+    5.: (36, 97, 247, 255),      # #2461f7 Eaux de surface
+    6.: (255, 255, 115, 255),    # #ffff73 Couvert herbacé en rotation
+    7.: (233, 255, 190, 255),    # #e9ffbe Couvert herbacé continu
+    8.: (0, 50, 0, 255),         # #003200 Arbres résineux (> 3m)
+    9.: (40, 200, 40, 255),      # #28c828 Arbres feuillus (> 3m)
+    80.: (0, 120, 0, 255),       # #007800 Arbustes résineux (≤ 3m)
+    90.: (183, 232, 176, 255),   # #b7e8b0 Arbustes feuillus (≤ 3m)
+}
+
+WALOUS_OCS2MANNING = {1.: 0.02, # Revêtement artificiel au sol
+                     2.: 0.025, # Constructions artificielles hors sol
+                     3.: 0.04,  # Réseau ferroviaire
+                     4.: 0.04,  # Sols nus
+                     5.: 0.033,  # Eaux de surface
+                     6.: 0.04,  # Couvert herbacé en rotation dans l'année
+                     7.: 0.03,  # Couvert herbacé toute l'année
+                     8.: 0.04, # Résineux (> 3m)
+                     9.: 0.04,  # Feuillus (> 3m)
+                     80.: 0.03, # Résineux (≤ 3m)
+                     90.: 0.03, # Feuillus (≤ 3m)
+                     }
+
+
+WALOUS_OCS2HYDROLOGY = {1.: 4., # Revêtement artificiel au sol
+                    2.: 4., # Constructions artificielles hors sol
+                    3.: 4., # Réseau ferroviaire
+                    4.: 4., # Sols nus
+                    5.: 5., # Eaux de surface
+                    6.: 3., # Couvert herbacé en rotation dans l'année
+                    7.: 2., # Couvert herbacé toute l'année
+                    8.: 1., # Résineux (> 3m)
+                    9.: 1., # Feuillus (> 3m)
+                    80.: 3., # Arbustes résineux (< 3m)
+                    90.: 3.  # Arbustes feuillus (< 3m)
+                    }
+
+# WALOUS UTS
+# ----------
 
 WALOUS_UTS_MAJ_NIV1 = {'Production primaire': 1.,
                    'Production secondaire': 2.,
@@ -85,59 +204,37 @@ WALOUS_UTS_COLORMAP_MAJ_NIV2 = {11.: (153,230,0,255),
                             66.: (206,136,102,255),
                             70.: (255,255,190,255)}
 
-WALOUS_OCS = { "Couvert herbacé en rotation dans l'année": 1.,
-              "Couvert herbacé toute l'année": 2.,
-              "Résineux (>3 m)": 3.,
-              "Feuillus (> 3m)": 4.,
-              "Résineux (<= 3m)": 5.,
-              "Feuillus (<= 3m)": 6.,
-              "Sols nus": 7.,
-              "Eaux de surface": 8.,
-              "Revêtement artificiels au sol": 9.,
-              "Réseau ferroviaire": 10.,
-              "Constructions artificielles hors sol": 11.}
-
-WALOUS_OCS_COLORMAP = {1.: (255,254,162,255), # Couvert herbacé en rotation dans l'année
-                        2.: (253,254,194,255), # Couvert herbacé toute l'année
-                        3.: (3,49,5,255), # Résineux (>3 m)
-                        4.: (55,198,61,255), # Feuillus (> 3m)
-                        5.: (18,119,23,255), # Résineux (<= 3m)
-                        6.: (185,231,179,255), # Feuillus (<= 3m)
-                        7.: (208,208,208,255), # Sols nus
-                        8.: (45,101,242,255), # Eaux de surface
-                        9.: (138,138,138,255), # Revêtement artificiels au sol
-                        10.: (78,78,78,255), # Réseau ferroviaire
-                        11.: (199,22,19,255), # Constructions artificielles hors sol
-                        }
-
-HYDROLOGY_LANDUSE_FR = {"forêt" : 1.,
-                        "prairie" : 2.,
-                        "culture" : 3.,
-                        "pavés/urbain" : 4.,
-                        "rivière" : 5.,
-                        "plan d'eau" : 6.,
-                     }
+WALOUS_UTS2MANNING_MAJ_NIV1 = {1.: 0.04, # Production primaire
+                           2.: 0.02, # Production secondaire
+                           3.: 0.02, # Production tertiaire
+                           4.: 0.03, # Réseaux de transport, Logistique et réseaux d'utilité publique
+                           5.: 0.025,# Usage résidentiel
+                           6.: 0.04, # Autres usages
+                           7.: 0.05} # Zones naturelles
 
-HYDROLOGY_LANDUSE_EN = {"forest" : 1.,
-                        "meadow" : 2.,
-                        "crop" : 3.,
-                        "paved/urban" : 4.,
-                        "river" : 5.,
-                        "water body" : 6.,
-                        }
+WALOUS_UTS2MANNING_MAJ_NIV2 = {11.: 0.04, # Agriculture
+                           12.: 0.04, # Sylviculture
+                           13.: 0.04, # Industries Extractives
+                           14.: 0.04, # Aquaculture et pêche
+                           20.: 0.03, # Production secondaire non définie
+                           21.: 0.03, # Industrie de matières premières
+                           22.: 0.03, # Industrie lourde
+                           23.: 0.03, # Industrie légère
+                           24.: 0.03, # Production d'énergie
+                           31.: 0.02, # Service commerciaux
+                           32.: 0.02, # Services financiers, spécialisés et d'information
+                           33.: 0.02, # Services publics
+                           34.: 0.02, # Services culturels, Services de loisirs et Services récréatifs
+                           41.: 0.025, # Réseaux de transport
+                           42.: 0.02, # Services Logistiques et d'entreposage
+                           43.: 0.025, # Réseau d'utilité publique
+                           51.: 0.025, # Usage résidentiel permanent
+                           52.: 0.025, # Usage résidentiel avec d'autres usages compatibles
+                           53.: 0.025, # Autres usages résidentiels
+                           62.: 0.04, # Zones abandonnées
+                           66.: 0.04, # Usage inconnu
+                           70.: 0.05} # Zones naturelles
 
-MAPPING_WALOUS_TO_HYDROLOGY = {1.: 2., # Couvert herbacé en rotation dans l'année
-                       2.: 2., # Couvert herbacé toute l'année
-                       3.: 1., # Résineux (>3 m)
-                       4.: 1., # Feuillus (> 3m)
-                       5.: 3., # Résineux (<= 3m)
-                       6.: 3., # Feuillus (<= 3m)
-                       7.: 3., # Sols nus
-                       8.: 5., # Eaux de surface
-                       9.: 4., # Revêtement artificiels au sol
-                       10.: 4., # Constructions artificielles hors sol
-                       11.: 4., # Réseau ferroviaire
-                    }
 
 def get_palette_walous_uts(which:Literal['MAJ_NIV1', 'MAJ_NIV2']) -> wolfpalette:
     """
@@ -177,7 +274,7 @@ def get_palette_walous_ocs() -> wolfpalette:
 
     return locpal
 
-def update_palette_walous_uts(which:Literal['MAJ_NIV1', 'MAJ_NIV2'], pal:wolfpalette):
+def update_palette_walous_uts(which:Literal['UTS_MAJ_NIV1', 'UTS_MAJ_NIV2'], pal:wolfpalette):
     """
     Update the palette for WALOUS MAJ_NIV1
 
@@ -185,13 +282,13 @@ def update_palette_walous_uts(which:Literal['MAJ_NIV1', 'MAJ_NIV2'], pal:wolfpal
     :return : updated palette
     """
 
-    if which == 'MAJ_NIV1':
+    if which in ['UTS_MAJ_NIV1', 'MAJ_NIV1']:
         for k, v in WALOUS_UTS_COLORMAP_MAJ_NIV1.items():
-            pal.set_values_colors(values = list(WALOUS_UTS_MAJ_NIV1.values()),
+            pal.set_values_colors(values = [val - 1e-6 for val in list(WALOUS_UTS_MAJ_NIV1.values())],
                                  colors = list(WALOUS_UTS_COLORMAP_MAJ_NIV1.values()))
-    elif which == 'MAJ_NIV2':
+    elif which in ['UTS_MAJ_NIV2', 'MAJ_NIV2']:
         for k, v in WALOUS_UTS_COLORMAP_MAJ_NIV2.items():
-            pal.set_values_colors(values = list(WALOUS_UTS_MAJ_NIV2.values()),
+            pal.set_values_colors(values = [val - 1e-6 for val in list(WALOUS_UTS_MAJ_NIV2.values())],
                                  colors = list(WALOUS_UTS_COLORMAP_MAJ_NIV2.values()))
 
     pal.interval_cst = True
@@ -208,7 +305,7 @@ def update_palette_walous_ocs(pal:wolfpalette):
     """
 
     for k, v in WALOUS_OCS_COLORMAP.items():
-        pal.set_values_colors(values = list(WALOUS_OCS_COLORMAP.keys()),
+        pal.set_values_colors(values = [val - 1e-6 for val in list(WALOUS_OCS_COLORMAP.keys())],
                              colors = list(WALOUS_OCS_COLORMAP.values()))
 
     pal.interval_cst = True
@@ -217,71 +314,6 @@ def update_palette_walous_ocs(pal:wolfpalette):
     return 0
 
 
-WALOUS2MANNING_MAJ_NIV1 = {1.: 0.04, # Production primaire
-                           2.: 0.02, # Production secondaire
-                           3.: 0.02, # Production tertiaire
-                           4.: 0.03, # Réseaux de transport, Logistique et réseaux d'utilité publique
-                           5.: 0.025,# Usage résidentiel
-                           6.: 0.04, # Autres usages
-                           7.: 0.05} # Zones naturelles
-
-WALOUS2MANNING_MAJ_NIV2 = {11.: 0.04, # Agriculture
-                           12.: 0.04, # Sylviculture
-                           13.: 0.04, # Industries Extractives
-                           14.: 0.04, # Aquaculture et pêche
-                           20.: 0.03, # Production secondaire non définie
-                           21.: 0.03, # Industrie de matières premières
-                           22.: 0.03, # Industrie lourde
-                           23.: 0.03, # Industrie légère
-                           24.: 0.03, # Production d'énergie
-                           31.: 0.02, # Service commerciaux
-                           32.: 0.02, # Services financiers, spécialisés et d'information
-                           33.: 0.02, # Services publics
-                           34.: 0.02, # Services culturels, Services de loisirs et Services récréatifs
-                           41.: 0.025, # Réseaux de transport
-                           42.: 0.02, # Services Logistiques et d'entreposage
-                           43.: 0.025, # Réseau d'utilité publique
-                           51.: 0.025, # Usage résidentiel permanent
-                           52.: 0.025, # Usage résidentiel avec d'autres usages compatibles
-                           53.: 0.025, # Autres usages résidentiels
-                           62.: 0.04, # Zones abandonnées
-                           66.: 0.04, # Usage inconnu
-                           70.: 0.05} # Zones naturelles
-
-WALOUSOCS2MANNING = {1.: 0.04, # Couvert herbacé en rotation dans l'année
-                     2.: 0.04, # Couvert herbacé toute l'année
-                     3.: 0.04, # Résineux (>3 m)
-                     4.: 0.04, # Feuillus (> 3m)
-                     5.: 0.04, # Résineux (<= 3m)
-                     6.: 0.04, # Feuillus (<= 3m)
-                     7.: 0.02, # Sols nus
-                     8.: 0.033, # Eaux de surface
-                     9.: 0.02, # Revêtement artificiels au sol
-                     10.: 0.02, # Réseau ferroviaire
-                     11.: 0.02} # Constructions artificielles hors sol
-
-"""
-Hydrology classification for Land Use in Wallonia:
-1 = forêt
-2 = prairie
-3 = culture
-4 = pavés/urbain
-5 = rivière
-6 = plan d'eau
-"""
-WALOUS2HYDROLOGY = {1.: 2, # Couvert herbacé en rotation dans l'année
-                    2.: 2, # Couvert herbacé toute l'année
-                    3.: 1, # Résineux (>3 m)
-                    4.: 1, # Feillus (> 3m)
-                    5.: 3, # Résineux (<= 3m)
-                    6.: 3, # Feuillus (<= 3m)
-                    7.: 4, # Sols nus
-                    8.: 5, # Eaux de surface
-                    9.: 4, # Revêtement artificiels au sol
-                    10.: 4, # Réseau ferroviaire
-                    11.: 4, # Constructions artificielles hors sol
-                    }
-
 class Walous_data():
     """
     La donnée Walous est liée à l'utilisation des sols en Wallonie
@@ -297,18 +329,21 @@ class Walous_data():
     def __init__(self,
                  dir_data:str = '',
                  fn:str = 'WAL_UTS__2018_L72',
+                 extension:str = '.shp',
                  bounds:Union[list[float, float, float, float],list[list[float, float], list[float, float]]] = None) -> None:
         """
         Constructor
 
         :param dir_data : directory of the data
-        :param fn : filename without extension (shp)
+        :param fn : filename without extension
+        :param extension : extension of the file (.shp, .gpkg, ...)
         :param bounds : Two ways to set spatial bounds -- [xmin, ymin, xmax, ymax] or [[xmin, xmax], [ymin, ymax]]
 
         """
 
-        self._dir = dir_data    # directory of the data
-        self._fn  = fn          # filename without extension
+        self._dir = str(dir_data)    # directory of the data
+        self._fn  = str(fn)          # filename without extension
+        self._extension = extension # extension of the file
         self._gdf = None        # geopandas dataframe
 
         if bounds is not None:
@@ -328,7 +363,7 @@ class Walous_data():
         if self._gdf is None or force:
             assert self._dir!="" and self._fn != ''
 
-            filepath = (Path(self._dir)/ self._fn).with_suffix('.shp')
+            filepath = (Path(self._dir)/ self._fn).with_suffix(self._extension)
 
             if filepath.exists():
                 if bounds is not None:
@@ -369,6 +404,8 @@ class Walous_data():
 
             chdir(detsdir)
 
+            fnout = str((Path(fnout)).with_suffix(self._extension))
+
             self._gdf.to_file(basename(fnout))
 
             chdir(curdir)
@@ -390,15 +427,15 @@ class Walous_data():
 
     def rasterize(self,
                   bounds:Union[list[float, float, float, float],list[list[float, float], list[float, float]]],
-                  layer:Literal['MAJ_NIV1','MAJ_NIV2'] ='MAJ_NIV1',
+                  layer:Literal['UTS_MAJ_NIV1','UTS_MAJ_NIV2'] = 'UTS_MAJ_NIV1',
                   fn_out:str = 'out.tif',
                   pixel_size:float = 0.5,
                   NoData_value:float = -99999.,
-                  num_type = gdal.GDT_Float32
+                  num_type = gdal.GDT_Float32,
                   ):
 
         """
-        Rasterization of polygon data to tif
+        Rasterization of polygon data to tif.
 
         :param bounds : [xmin, ymin, xmax, ymax] or [[xmin, xmax], [ymin, ymax]]
         :param layer : layer to rasterize
@@ -408,6 +445,8 @@ class Walous_data():
         :param num_type : type of the number
         """
 
+        layer = layer.upper().replace('UTS_', '')
+
         if bounds is None:
             logging.error('Bounds must be set')
             return None
@@ -438,10 +477,10 @@ class Walous_data():
         try:
 
             # Add a new column for mapping based on the desired layer
-            self._gdf['Mapping'] = np.float32(self._gdf[layer])
-
             if layer == 'MAJ_NIV2':
-                self._gdf['Mapping'] = np.float32(self._gdf['Mapping'].replace('_', ''))
+                self._gdf['Mapping'] = np.float32(self._gdf[layer].replace('_', ''))
+            else:
+                self._gdf['Mapping'] = np.float32(self._gdf[layer])
 
             source_ds:ogr.DataSource
             source_layer:ogr.Layer
@@ -485,7 +524,7 @@ class Walous_data():
 class  DlgMapWalous2Manning(wx.Dialog):
     """ Modal dialog for mapping WALOUS value to another ones """
 
-    def __init__(self, parent, title:str = _("Mapping WALOUS value to ..."), which:str = 'MAJ_NIV1'):
+    def __init__(self, parent, title:str = _("Mapping WALOUS value to ..."), which:Literal['UTS_MAJ_NIV1', 'UTS_MAJ_NIV2', 'OCS_MANNING', 'OCS_HYDROLOGY'] = 'UTS_MAJ_NIV1'):
 
         super(DlgMapWalous2Manning, self).__init__(parent, title=title, size=(450, 400))
 
@@ -495,7 +534,7 @@ class  DlgMapWalous2Manning(wx.Dialog):
 
         self._table = wx.grid.Grid(panel)
 
-        if which == 'MAJ_NIV1':
+        if which == 'UTS_MAJ_NIV1':
             self._table.CreateGrid(len(WALOUS_UTS_MAJ_NIV1), 3)
             self._table.SetColLabelValue(0, _("Name"))
             self._table.SetColLabelValue(1, _("Value - UTS"))
@@ -505,11 +544,11 @@ class  DlgMapWalous2Manning(wx.Dialog):
             for i, (k, v) in enumerate(WALOUS_UTS_MAJ_NIV1.items()):
                 self._table.SetCellValue(i, 0, str(k))
                 self._table.SetCellValue(i, 1, str(v))
-                self._table.SetCellValue(i, 2, str(WALOUS2MANNING_MAJ_NIV1[v]))
+                self._table.SetCellValue(i, 2, str(WALOUS_UTS2MANNING_MAJ_NIV1[v]))
                 self._table.SetCellAlignment(i, 1, wx.ALIGN_CENTER, wx.ALIGN_CENTER)
                 self._table.SetCellAlignment(i, 2, wx.ALIGN_CENTER, wx.ALIGN_CENTER)
 
-        elif which == 'MAJ_NIV2':
+        elif which == 'UTS_MAJ_NIV2':
             self._table.CreateGrid(len(WALOUS_UTS_MAJ_NIV2), 3)
             self._table.SetColLabelValue(0, _("Name"))
             self._table.SetColLabelValue(1, _("Value - UTS"))
@@ -519,7 +558,33 @@ class  DlgMapWalous2Manning(wx.Dialog):
             for i, (k, v) in enumerate(WALOUS_UTS_MAJ_NIV2.items()):
                 self._table.SetCellValue(i, 0, str(k))
                 self._table.SetCellValue(i, 1, str(v))
-                self._table.SetCellValue(i, 2, str(WALOUS2MANNING_MAJ_NIV2[v]))
+                self._table.SetCellValue(i, 2, str(WALOUS_UTS2MANNING_MAJ_NIV2[v]))
+                self._table.SetCellAlignment(i, 1, wx.ALIGN_CENTER, wx.ALIGN_CENTER)
+                self._table.SetCellAlignment(i, 2, wx.ALIGN_CENTER, wx.ALIGN_CENTER)
+
+        elif which == 'OCS_MANNING':
+            self._table.CreateGrid(len(WALOUS_OCS2MANNING), 3)
+            self._table.SetColLabelValue(0, _("Name"))
+            self._table.SetColLabelValue(1, _("Value - OCS"))
+            self._table.SetColLabelValue(2, _("Manning 'n'"))
+            self._table.HideRowLabels()
+            for i, (k, v) in enumerate(WALOUS_OCS.items()):
+                self._table.SetCellValue(i, 0, str(k))
+                self._table.SetCellValue(i, 1, str(v))
+                self._table.SetCellValue(i, 2, str(WALOUS_OCS2MANNING[v]))
+                self._table.SetCellAlignment(i, 1, wx.ALIGN_CENTER, wx.ALIGN_CENTER)
+                self._table.SetCellAlignment(i, 2, wx.ALIGN_CENTER, wx.ALIGN_CENTER)
+
+        elif which == 'OCS_HYDROLOGY':
+            self._table.CreateGrid(len(WALOUS_OCS2HYDROLOGY), 3)
+            self._table.SetColLabelValue(0, _("Name"))
+            self._table.SetColLabelValue(1, _("Value - OCS"))
+            self._table.SetColLabelValue(2, _("LandUse - Hydrology"))
+            self._table.HideRowLabels()
+            for i, (k, v) in enumerate(WALOUS_OCS2HYDROLOGY.items()):
+                self._table.SetCellValue(i, 0, str(k))
+                self._table.SetCellValue(i, 1, str(v))
+                self._table.SetCellValue(i, 2, str(WALOUS_OCS2HYDROLOGY[v]))
                 self._table.SetCellAlignment(i, 1, wx.ALIGN_CENTER, wx.ALIGN_CENTER)
                 self._table.SetCellAlignment(i, 2, wx.ALIGN_CENTER, wx.ALIGN_CENTER)
 
@@ -575,11 +640,11 @@ class  DlgMapWalousOCS2Manning(DlgMapWalous2Manning):
 
     def __init__(self, parent, title:str = _("Mapping WALOUS value to ...")):
 
-        super(DlgMapWalousOCS2Manning, self).__init__(parent, title= title)
+        super(DlgMapWalousOCS2Manning, self).__init__(parent, title= title, which='OCS_MANNING')
 
         self._table.SetColLabelValue(2, _("Land Use - OCS"))
 
-        for i, (k, v) in enumerate(WALOUSOCS2MANNING.items()):
+        for i, (k, v) in enumerate(WALOUS_OCS2MANNING.items()):
             self._table.SetColLabelValue(1, _("Value - OCS"))
             self._table.SetCellValue(i, 2, str(v))
 
@@ -593,11 +658,11 @@ class  DlgMapWalous2Hydrology(DlgMapWalous2Manning):
 
     def __init__(self, parent, title:str = _("Mapping WALOUS value to ...")):
 
-        super(DlgMapWalous2Hydrology, self).__init__(parent, title= title)
+        super(DlgMapWalous2Hydrology, self).__init__(parent, title= title, which='OCS_HYDROLOGY')
 
         self._table.SetColLabelValue(2, _("Land Use - OCS"))
 
-        for i, (k, v) in enumerate(WALOUS2HYDROLOGY.items()):
+        for i, (k, v) in enumerate(WALOUS_OCS2HYDROLOGY.items()):
             self._table.SetColLabelValue(1, _("Value - OCS"))
             self._table.SetCellValue(i, 2, str(v))
 
@@ -692,3 +757,143 @@ class Walous_OCS_Legend(wx.Dialog):
 
     def close(self):
         self.Destroy()
+
+
+def get_array_WALOUS_OCS(fname:Path, spatial_res:float, bounds:Union[list[float, float, float, float],list[list[float, float], list[float, float]]]) -> Union[WolfArray, None]:
+    """ Get WALOUS OCS as WolfArray
+
+    :param fname: path to the WALOUS OCS file
+    :param spatial_res: desired spatial resolution (in m)
+    :param bounds: [xmin, ymin, xmax, ymax] or [[xmin, xmax], [ymin, ymax]]
+    :return: WolfArray or None if error
+    """
+
+    if not fname.exists():
+        logging.error(_('WALOUS OCS file not found: {}').format(fname))
+        return None
+
+    # Convert bounds to [[xmin, xmax], [ymin, ymax]] if needed
+    if len(bounds)==4:
+        bounds = [[bounds[0], bounds[2]], [bounds[1], bounds[3]]]
+
+    # Convert bvounds to list to be able to modify it
+    bounds = [list(bounds[0]), list(bounds[1])]
+
+    bounds_copy = [list(bounds[0]), list(bounds[1])]
+
+    header_OCS = header_wolf.read_header(fname)
+    if header_OCS.dx != spatial_res:
+        # Adapt bounds to ensure that the rebin will be correct.
+        # If spatial_res is a multiple of header_OCS.dx, no need to change bounds.
+        # If not, change the bounds to ensure that the crop will include data in all footprints.
+
+        if (bounds[0][0] - header_OCS.origx) % header_OCS.dx != 0:
+            bounds[0][0] = header_OCS.origx + ((bounds[0][0] - header_OCS.origx) // header_OCS.dx) * header_OCS.dx
+        if (bounds[0][1] - bounds[0][0]) % header_OCS.dx != 0:
+            bounds[0][1] = bounds[0][0] + ((bounds[0][1] - bounds[0][0]) // header_OCS.dx + 1) * header_OCS.dx
+        if (bounds[1][0] - header_OCS.origy) % header_OCS.dy != 0:
+            bounds[1][0] = header_OCS.origy + ((bounds[1][0] - header_OCS.origy) // header_OCS.dy) * header_OCS.dy
+        if (bounds[1][1] - bounds[1][0]) % header_OCS.dy != 0:
+            bounds[1][1] = bounds[1][0] + ((bounds[1][1] - bounds[1][0]) // header_OCS.dy + 1) * header_OCS.dy
+
+    locwalous = WolfArray(fname=fname, crop = [bounds[0][0], bounds[0][1], bounds[1][0], bounds[1][1]])
+
+    if locwalous.dx != spatial_res:
+        locwalous.rebin(spatial_res / locwalous.dx, operation='min')
+        logging.info(_('Rebin to {} m because original data are {} m').format(spatial_res, locwalous.dx))
+        locwalous = WolfArray(mold=locwalous, crop = bounds_copy)
+
+    return locwalous
+
+def get_array_WALOUS_OCS_2023_10m(spatial_res:float, bounds:Union[list[float, float, float, float],list[list[float, float], list[float, float]]]) -> Union[WolfArray, None]:
+    """ Get WALOUS OCS 2023 as WolfArray
+
+    :param spatial_res: desired spatial resolution (in m)
+    :param bounds: [xmin, ymin, xmax, ymax] or [[xmin, xmax], [ymin, ymax]]
+    :return: WolfArray or None if error
+    """
+
+    from wolfhece.pydownloader import toys_dataset
+
+    fname = toys_dataset(dir='Walous_OCS', file='WALOUS_2023_lbt72_10m.tif')
+
+    return get_array_WALOUS_OCS(fname, spatial_res, bounds)
+
+def get_array_WALOUS_OCS_2023_4m(spatial_res:float, bounds:Union[list[float, float, float, float],list[list[float, float], list[float, float]]]) -> Union[WolfArray, None]:
+    """ Get WALOUS OCS 2023 as WolfArray
+
+    :param spatial_res: desired spatial resolution (in m)
+    :param bounds: [xmin, ymin, xmax, ymax] or [[xmin, xmax], [ymin, ymax]]
+    :return: WolfArray or None if error
+    """
+
+    from wolfhece.pydownloader import toys_dataset
+
+    fname = toys_dataset(dir='Walous_OCS', file='WALOUS_2023_lbt72_4m.tif')
+
+    return get_array_WALOUS_OCS(fname, spatial_res, bounds)
+
+def get_array_WALOUS_OCS_2020_10m(spatial_res:float, bounds:Union[list[float, float, float, float],list[list[float, float], list[float, float]]]) -> Union[WolfArray, None]:
+    """ Get WALOUS OCS 2020 as WolfArray
+
+    :param spatial_res: desired spatial resolution (in m)
+    :param bounds: [xmin, ymin, xmax, ymax] or [[xmin, xmax], [ymin, ymax]]
+    :return: WolfArray or None if error
+    """
+
+    from wolfhece.pydownloader import toys_dataset
+
+    fname = toys_dataset(dir='Walous_OCS', file='WALOUS_2020_lbt72_10m.tif')
+
+    return get_array_WALOUS_OCS(fname, spatial_res, bounds)
+
+def get_array_WALOUS_OCS_2020_4m(spatial_res:float, bounds:Union[list[float, float, float, float],list[list[float, float], list[float, float]]]) -> Union[WolfArray, None]:
+    """ Get WALOUS OCS 2020 as WolfArray
+
+    :param spatial_res: desired spatial resolution (in m)
+    :param bounds: [xmin, ymin, xmax, ymax] or [[xmin, xmax], [ymin, ymax]]
+    :return: WolfArray or None if error
+    """
+
+    from wolfhece.pydownloader import toys_dataset
+
+    fname = toys_dataset(dir='Walous_OCS', file='WALOUS_2020_lbt72_4m.tif')
+
+    return get_array_WALOUS_OCS(fname, spatial_res, bounds)
+
+def get_array_WALOUS_UTS(fname:Path, fnout:Path, spatial_res:float, bounds:Union[list[float, float, float, float],list[list[float, float], list[float, float]]], which:Literal['UTS_MAJ_NIV1', 'UTS_MAJ_NIV2'] = 'UTS_MAJ_NIV1') -> Union[WolfArray, None]:
+    """ Get WALOUS UTS as WolfArray
+
+    :param fname: path to the WALOUS UTS file
+    :param spatial_res: desired spatial resolution (in m)
+    :param bounds: [xmin, ymin, xmax, ymax] or [[xmin, xmax], [ymin, ymax]]
+    :param which: which level to use ('UTS_MAJ_NIV1' or 'UTS_MAJ_NIV2')
+    :return: WolfArray or None if error
+    """
+
+    fname = Path(fname)
+    fnout = Path(fnout)
+
+    if not fname.exists():
+        logging.error(_('WALOUS UTS file not found: {}').format(fname))
+        return None
+
+    # Convert bounds to [[xmin, xmax], [ymin, ymax]] if needed
+    if len(bounds)==4:
+        bounds = [[bounds[0], bounds[2]], [bounds[1], bounds[3]]]
+
+    # Convert bvounds to list to be able to modify it
+    bounds = [list(bounds[0]), list(bounds[1])]
+
+    bounds_copy = [list(bounds[0]), list(bounds[1])]
+
+    locwalous = Walous_data(fname.parent, fname.name, extension=fname.suffix, bounds=bounds)
+    ret = locwalous.rasterize(bounds=bounds, layer= which, fn_out=fnout.with_suffix('.tif'), pixel_size=spatial_res)
+
+    ret = Path(ret)
+
+    if ret.exists():
+        locwalous_array = WolfArray(fname=ret)
+        return locwalous_array
+    else:
+        return None