|
The MFD algorithm in ArcMap 10.6.1 is illustrated here wtih two simple rasters.
Flow direction was calculated with:
arcpy.gp.FlowDirection_sa("demf", "littledir.crf", "FORCE", "", "MFD")
Because flow was forced off the edge, the only values shown from:
arcpy.gp.FlowAccumulation_sa("littledir.crf", "littleacc.crf", "", "FLOAT", "MFD")
are portions of the unit area of the center cell.
1) While the digonal flows might be expected to be reduced by .7071, the algorithm makes them inexplicibly tiny.
2) There is an inexplicably large weight given to steeper drops.
3) The floatingpoint elevation values are rounded to three decimal places.
This is why the ESRI MFD flows differ so from my smootheed flows in these
Antarctic Dry Valleys figures.
back to the classic stability page
back to the main fuzziness page
back to the new stability page
on which Harvey implements the Montgomery/Dietrich fuzzy flowaccumulation algorithm in python
within the SHALSTAB slope stability calculator.