I’ve got a vexing GIS analysis problem which is conceptually fairly straightforward, but in reality more challenging to figure out….
Basically, I need to calculate the VOLUME of floodwater within a 100 year floodplai, using ArcGIS 10.0 (all extensions available).
My input datasets are:
LiDAR DEM (1m pixel or 3m pixel) of bare ground (raster)
100 year floodplain (output of USACE HEC-RAS analysis), Polygon Shapefile, but attribute table has no useful content.
So far I have overlaid the floodplain polygons over the LiDAR raster, [using Extract by Mask tool from Spatial Analyst], so that I have the ground elevations are along the outside edge of the floodplain (think of the top of the bathtub). Key concept is that the top of the floodplain gradually falls in elevation as you move from the headwaters down through the watershed.
Then, (this seems to be the impossible part) I need to figure out how I can take the elevations along the outside edge pixels of the raster floodplain (output of previous Extract by Mask) and extend them to the equivalent elevation on the other side of the floodplain, and “raise” all the pixels in between the 2 outside edges of the floodplain to the “top level” of this 100-year floodplain.
Once I have this “top elevation” of the floodplain water, I should be able to use the CUT AND FILL tools in ArcGIS to simply subtract the LiDAR ground elevations from the “top floodplain elevation” and compute a volume between these two rasters.
I’ve been researching this for days and have yet to find a solution.
Perhaps I’m thinking of this completely wrong?
I would have expected to find others who have dealt with this problem, but while there are similar problems, none appear to offer a solution.