Skip to content
This repository has been archived by the owner on Jan 9, 2024. It is now read-only.

Curve Number raster creation: handle null values in either landcover or soils datasets #13

Open
gassc opened this issue Oct 28, 2020 · 0 comments

Comments

@gassc
Copy link
Member

gassc commented Oct 28, 2020

For the tool that supports the creation of a curve number, have a way to handle null values in soils hydro classes and landcover classes. It's most likely we'll encounter soil layers without hydrologic classifications. Possible solutions:

  • determine and build in the logic for an appropriate fallback
  • expose a user-provided parameter in the toolbox for providing a fallback
  • ensure that the lookup table itself can be used to handle nulls, e.g., with a blank soil hydro value in the lookup row, the tool will fall back to the provided curve number for the landcover class when the soil is null:
land cover class soil hydro curve number
2 45

This could also work to handle conditions where landcover is null but soil is known, or both are null.

That last approach is probably the easiest to implement. The change would happen around this line:

cn_raster_component = Con((landcover_raster_obj == rec.utc) & (soils_raster == rec.soil), rec.cn, 0)

Using the ArcPy Spatial Analyst / map algebra toolbox, an IsNull function following or proceeding that line could do the work.

Also: create docs around the Curve Number raster creation process, describing how that lookup table is made and adapted for other landcover datasets.

Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant