[GIS] How to extend a (river water level) raster, transfering its values to the added area

bufferfloodqgisrasterrivers

Exploring the raster tools in QGIS and related software, I found raster buffer functions to surround existing rasters. The values of the raster buffer are restricted to be distances from the interior raster or constant however. I would rather like to derive the values of the added cells from the boundary of the interior raster…

Data: 1. river flood level raster; 2. polygon at the side of 1.; 3. bended river axis as polyline; 4. DEM.

Aim: I'd like to laterally extend the river flood level raster, which is restricted between dikes, constraining the flood plain. Some dikes are to be removed and additional flood retention area (delineated by polygon 2.) will be flooded.

It's not about modelling the flow or change of water level. – I just want to transfer the given flood levels at the boundary of the raster to the neighboring area, though keeping the overal water level inclination along the main flow direction of the river.

The laterally extended flood water level will be used to approximate inundation depths by subtracting DEM values in a second step. I don't need very precise water level results, because the study is about sedimentation and soil contamination mainly.

Any idea how to laterally extend the (flood level) raster values by use of QGIS & friends?

Best Answer

Original dataset showing flood depth and polygon I've created along it:

enter image description here

Let's produce points on the boundary of flood plain using flood level raster called RL, raster calculator plus some tools:

<code>enter image description here</code>

Next - assign proximity areas to above points and convert those areas into corresponding flood level:

enter image description here

Final step using raster calculator to define flood level where it is above DEM. If you'll get to here, you'll have no trouble to do it.

Result below compares flood LEVEL in the river and polygon (I dropped original elevation inside polygon by 10 m). As one can see it picks the trend, however every single imperfection in original flood level get exaggerated in output:

enter image description here

Related Question