Automatically determine downstream grid cell ID

TaylorJordan — Sat, 22 Jul 2023 19:29:01 GMT

Hello,

I'll start by explaining a little about the project I'm working on. I'm trying to analyze a watershed (~60 sq. km) using a lumped rainfall-runoff model. To improve the resolution of the watershed, I've created a grid of polygons with a cell size of 100mX100m (each with an elevation assigned to it). The grid cells act as subcatchments in the model and each subcatchment needs to have an outlet assigned to it. With the size of my watershed and number of grid cells, I need the outlet IDs to be automatically determined as opposed to doing it all myself.

I'm sure there's a way for ArcGIS Pro to automatically tell me the name/ID of the neighboring cell that receives runoff. The Flow Direction tool doesn't quite tell me what I need, since it just determines which direction each cell drains. Does anyone know of a way to do this? Ideally, the final result would be a column in the Attribute Table for my polygon grid called "Outlet", which would have the name of the neighboring cell that that cell drains to.

Hopefully that all makes sense. I can't seem to find a solution to this in the rainfall-runoff model I'm using, so I'm hoping to process it in ArcGIS first and then import into my model. That way, the outlet will already be defined for each cell to use in the model. For reference, I'm using ArcGIS Pro 3.0.2.

I appreciate any help.

Thanks!

Taylor

Re: Automatically determine downstream grid cell ID

JohannesLindner — Mon, 24 Jul 2023 17:40:00 GMT

Hmmm, this works, but I really hope I missed an easier solution.

# extract your watershed from the larger digital elevation model ws_dem = arcpy.sa.ExtractByMask("DEM", "Cutter") # resample to 100x100 ws_dem_resampled = arcpy.management.Resample(ws_dem, "ws_dem_resampled", "100 100") # fill and flow dir ws_dem_filled = arcpy.sa.Fill(ws_dem_resampled) flow_dir = arcpy.sa.FlowDirection(ws_dem_filled) # create an id raster (numbered sequentially) cell_points = arcpy.conversion.RasterToPoint(ws_dem_filled, "cell_points", "Value") id_raster = arcpy.Raster(arcpy.conversion.PointToRaster(cell_points, "pointid", "id_raster", cellsize=flow_dir)) # calculate downstream id directions = { 1: (1, 0), 2: (1, 1), 4: (0, 1), 8: (-1, 1), 16: (-1, 0), 32: (-1, -1), 64: (0, -1), 128: (1, -1) } ds_id_raster = arcpy.Raster(id_raster.getRasterInfo()) for x, y in ds_id_raster: d = flow_dir[x, y] i = id_raster[x,y ] try: dy, dx = directions[d] i_ = id_raster[x+dx, y+dy] except: continue ds_id_raster[x, y] = i_ ds_id_raster.save("ds_id_raster") # convert raster to polygons cells = arcpy.conversion.RasterToPolygon(id_raster, "Cells", "NO_SIMPLIFY") # extract values into the points and join arcpy.sa.ExtractMultiValuesToPoints(cell_points, [(id_raster, "Cell"), (ds_id_raster, "CellDS"), (flow_dir, "FlowDir"), (ws_dem_filled, "Elevation")]) arcpy.management.JoinField(cells, "gridcode", cell_points, "pointid", ["Cell", "CellDS", "FlowDir", "Elevation"])

This script will:

resample your dem

calculate flow direction

Give a unique id to each raster cell

Use the flow direction to get the downstream cell id

Extract all this information into a polygon grid

topic Automatically determine downstream grid cell ID in Geoprocessing Questions

Automatically determine downstream grid cell ID

Re: Automatically determine downstream grid cell ID