Délimitation « Simple » de la plaine inondable en utilisant uniquement un MNT LiDAR et une couche de cours d'eau

splatt_wfp

Bonjour,

Je n'ai jamais utilisé les outils ArcHydro auparavant et j'espère que quelqu'un pourra m'aider avec une délimitation basique de la plaine inondable. J'ai un MNT LiDAR et une couche de cours d'eau dérivée avec des valeurs Z. J'espère cartographier l'étendue de la plaine inondable à partir de ces lignes centrales de cours d'eau en fonction d'une élévation de 3 m du niveau de l'eau. Rien de plus compliqué que cela. D'après ce que j'ai lu jusqu'à présent, je dois "Flood Out WSE le long des cours d'eau par étapes incrémentielles"... Mais je ne suis pas sûr des outils ArcHydro à utiliser pour cela.

Quelqu'un peut-il aider ?

Merci,

Steve

 

   

DeanDjokic

Steve,

Depending on the quality of your DEM and streams, this could be an "easy" exercise.  If your DEM and the flow direction match the placement and elevation of your streams, then you can directly use the HAND (Height Above Nearest Drainage) approach.  You will need the following rasters (I suggest that you use .tif format for your rasters):

1. DEM (i'll call it DEM.tif).  Should be in correct horizontal/vertical units - e.g. meters.
2. Flow direction raster derived from the DEM (FDR.tif) that "flows" into your streams.
3. Stream raster (STR.tif).  This can have value of 1 where stream cells are and nodata elsewhere.  The elevation of the streams (thalweg) should be reflected in the DEM.

The steps are then (these are all core Spatial Analyst functions - you do not need Arc Hydro):

1. Use "Flow Distance" (Hydrology toolset) to calculate the HAND (HAND.tif) surface (that is the "output raster".  Specify inputs stream raster (STR), surface raster (DEM) and flow direction raster (FDR). Keep all other parameters at their default values.
2. Once you have the HAND surface, to get the depth and extent of flooding for a specific depth (D) in the same vertical units as the DEM, use the simple conditional statement (use CON SA function from “Conditional” toolset) and specify:
   1. HAND.tif as your conditional raster.
   2. Expression as: Value is less than or equal to 3 (D is your flooding depth – in your case it is 3).
   3. 3 as your true raster (D is your flooding depth – in your case it is 3).
   4. Leave false raster entry empty (meaning all values outside of your flood extent will be nodata).
   5. Specify output raster (I like to put name of the depth in the output name – e.g. fe3m.tif – for flood extent of 3m, but it can be anything).

That’s it.  In the output raster, where you have values, area is flooded and where it is nodata, it is “dry”.  If you want to calculate depth of flooding for each cell, just use the MINUS function and use fe3m.tif and HAND.tif as inputs (in that order).

Hope this works for you.