I can't run Urban Flood Risk Mitigation

Software Support
1

Hello, it’s my first time using Invest,
I saved Parameter and data (.tgz file that I upload)

I am trying to run the urban flood risk mitigation model, but whenever I try to run this errors appears:

03/02/2025 17:48:21 pygeoprocessing.geoprocessing ERROR exception encountered in raster_calculator
Traceback (most recent call last):
File “natcap\invest\urban_flood_risk_mitigation.py”, line 892, in _lu_to_cn_op
File “numpy_core\fromnumeric.py”, line 455, in choose
File “numpy_core\fromnumeric.py”, line 57, in _wrapfunc
ValueError: invalid entry in choice array

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File “pygeoprocessing\geoprocessing.py”, line 475, in raster_calculator
File “natcap\invest\urban_flood_risk_mitigation.py”, line 899, in _lu_to_cn_op
ValueError: invalid entry in choice array
Check that the Soil Group raster does not contain values other than (1, 2, 3, 4)
03/02/2025 17:48:21 pygeoprocessing.geoprocessing INFO Waiting for raster stats worker result.
03/02/2025 17:48:21 taskgraph.Task ERROR Something went wrong when adding task create Curve Number raster (1), terminating taskgraph.
Traceback (most recent call last):
File “natcap\invest\urban_flood_risk_mitigation.py”, line 892, in _lu_to_cn_op
File “numpy_core\fromnumeric.py”, line 455, in choose
File “numpy_core\fromnumeric.py”, line 57, in _wrapfunc
ValueError: invalid entry in choice array

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File “taskgraph\Task.py”, line 674, in add_task
File “taskgraph\Task.py”, line 1093, in _call
File “src\pygeoprocessing\geoprocessing_core.pyx”, line 77, in pygeoprocessing.geoprocessing_core.gdal_use_exceptions.wrapper
File “src\pygeoprocessing\geoprocessing_core.pyx”, line 78, in pygeoprocessing.geoprocessing_core.gdal_use_exceptions.wrapper
File “pygeoprocessing\geoprocessing.py”, line 475, in raster_calculator
File “natcap\invest\urban_flood_risk_mitigation.py”, line 899, in _lu_to_cn_op
ValueError: invalid entry in choice array
Check that the Soil Group raster does not contain values other than (1, 2, 3, 4)
03/02/2025 17:48:21 natcap.invest.utils ERROR Exception while executing natcap.invest.urban_flood_risk_mitigation
Traceback (most recent call last):
File “natcap\invest\urban_flood_risk_mitigation.py”, line 892, in _lu_to_cn_op
File “numpy_core\fromnumeric.py”, line 455, in choose
File “numpy_core\fromnumeric.py”, line 57, in _wrapfunc
ValueError: invalid entry in choice array

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File “natcap\invest\utils.py”, line 186, in prepare_workspace
File “invest\cli.py”, line 472, in main
File “natcap\invest\urban_flood_risk_mitigation.py”, line 325, in execute
File “taskgraph\Task.py”, line 674, in add_task
File “taskgraph\Task.py”, line 1093, in _call
File “src\pygeoprocessing\geoprocessing_core.pyx”, line 77, in pygeoprocessing.geoprocessing_core.gdal_use_exceptions.wrapper
File “src\pygeoprocessing\geoprocessing_core.pyx”, line 78, in pygeoprocessing.geoprocessing_core.gdal_use_exceptions.wrapper
File “pygeoprocessing\geoprocessing.py”, line 475, in raster_calculator
File “natcap\invest\urban_flood_risk_mitigation.py”, line 899, in _lu_to_cn_op
ValueError: invalid entry in choice array
Check that the Soil Group raster does not contain values other than (1, 2, 3, 4)
03/02/2025 17:48:21 natcap.invest.utils INFO Elapsed time: 0.15s
03/02/2025 17:48:21 natcap.invest.utils INFO Execution finished; version: 3.14.3
Could you help me, please?

2

Hi @YohanaMora -

The log file includes this line:

Check that the Soil Group raster does not contain values other than (1, 2, 3, 4)

As the User Guide also notes, the Soil Group raster may only contain values 1-4. Have you checked that this is the case for your soil group input raster?

Sometimes soil group data contains values like 14, indicating dual soil groups A/D etc. These need to be converted into one of the valid values 1-4. The User Guide provides some guidance about what to consider when converting these values.

If it looks like your soil group raster only does have values 1-4, then it would be useful if you’d upload that file (or point us to it somewhere on cloud storage) and I’ll check it out.

~ Stacie

3

Thanks.

Hello, I changed the groupsoil raster but now the problem is this:
03/03/2025 11:47:46 pygeoprocessing.geoprocessing INFO poly intersection lookup approximately 2.2% complete on reprojected_aoi
03/03/2025 11:47:51 taskgraph.Task ERROR A taskgraph _task_executor failed on Task zonal_statistics over runoff_retention raster (9). Terminating taskgraph.
multiprocessing.pool.RemoteTraceback:
“”"
Traceback (most recent call last):
File “multiprocessing\pool.py”, line 125, in worker
File “src\pygeoprocessing\geoprocessing_core.pyx”, line 77, in pygeoprocessing.geoprocessing_core.gdal_use_exceptions.wrapper
File “src\pygeoprocessing\geoprocessing_core.pyx”, line 78, in pygeoprocessing.geoprocessing_core.gdal_use_exceptions.wrapper
File “pygeoprocessing\geoprocessing.py”, line 1803, in zonal_statistics
File “src\pygeoprocessing\geoprocessing_core.pyx”, line 77, in pygeoprocessing.geoprocessing_core.gdal_use_exceptions.wrapper
File “src\pygeoprocessing\geoprocessing_core.pyx”, line 78, in pygeoprocessing.geoprocessing_core.gdal_use_exceptions.wrapper
File “pygeoprocessing\geoprocessing.py”, line 2881, in calculate_disjoint_polygon_set
File “shapely\prepared.py”, line 57, in intersects
File “shapely\geometry\base.py”, line 719, in intersects
File “shapely\decorators.py”, line 77, in wrapped
File “shapely\predicates.py”, line 778, in intersects
shapely.errors.GEOSException: bad allocation
“”"

The above exception was the direct cause of the following exception:

Traceback (most recent call last):
File “taskgraph\Task.py”, line 461, in _task_executor
File “taskgraph\Task.py”, line 1090, in _call
File “multiprocessing\pool.py”, line 774, in get
shapely.errors.GEOSException: bad allocation
03/03/2025 11:47:51 natcap.invest.utils ERROR Exception while executing natcap.invest.urban_flood_risk_mitigation
multiprocessing.pool.RemoteTraceback:
“”"
Traceback (most recent call last):
File “multiprocessing\pool.py”, line 125, in worker
File “src\pygeoprocessing\geoprocessing_core.pyx”, line 77, in pygeoprocessing.geoprocessing_core.gdal_use_exceptions.wrapper
File “src\pygeoprocessing\geoprocessing_core.pyx”, line 78, in pygeoprocessing.geoprocessing_core.gdal_use_exceptions.wrapper
File “pygeoprocessing\geoprocessing.py”, line 1803, in zonal_statistics
File “src\pygeoprocessing\geoprocessing_core.pyx”, line 77, in pygeoprocessing.geoprocessing_core.gdal_use_exceptions.wrapper
File “src\pygeoprocessing\geoprocessing_core.pyx”, line 78, in pygeoprocessing.geoprocessing_core.gdal_use_exceptions.wrapper
File “pygeoprocessing\geoprocessing.py”, line 2881, in calculate_disjoint_polygon_set
File “shapely\prepared.py”, line 57, in intersects
File “shapely\geometry\base.py”, line 719, in intersects
File “shapely\decorators.py”, line 77, in wrapped
File “shapely\predicates.py”, line 778, in intersects
shapely.errors.GEOSException: bad allocation
“”"

The above exception was the direct cause of the following exception:

Traceback (most recent call last):
File “natcap\invest\utils.py”, line 186, in prepare_workspace
File “invest\cli.py”, line 472, in main
File “natcap\invest\urban_flood_risk_mitigation.py”, line 497, in execute
File “taskgraph\Task.py”, line 1299, in get
File “taskgraph\Task.py”, line 1272, in join
File “taskgraph\Task.py”, line 461, in _task_executor
File “taskgraph\Task.py”, line 1090, in _call
File “multiprocessing\pool.py”, line 774, in get
shapely.errors.GEOSException: bad allocation

I had changed to raster pixels for 300m X 300 m but the problem persists. ¿Can be my PC the problem?

This link contains the folder with parameter and data that I am using now:

Thank you so much for your help