Hello, everyone. When processing the annual water supply volume using the “invest” model version 3.14.1, I encountered the following error report. I would like to ask everyone how to solve this issue. Thank you all.
2024-04-03 00:29:47,648 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [E:\water4—3\InVEST-natcap.invest.annual_water_yield-log-2024-04-03–00_29_47.txt]
2024-04-03 00:29:47,649 (main) cli.main(459) Level 100 Starting model with parameters:
Arguments for InVEST natcap.invest.annual_water_yield 3.14.1:
biophysical_table_path E:\invest water yield\biophysical_table_gura1.csv
demand_table_path
depth_to_root_rest_layer_path E:\invest water yield\Processed\2020 water\dtb2020\rrld2020.tif
eto_path E:\invest water yield\Processed\2020 water\pet2020\etp2020.tif
lulc_path E:\invest water yield\Processed\2020 water\lucc2020\lucc_2020(42).tif
n_workers -1
pawc_path E:\invest water yield\Processed\2020 water\hwsd2020\pawc2020.tif
precipitation_path E:\invest water yield\Processed\Pre2020\Pre2020.tif
results_suffix 2020
seasonality_constant 0.532
sub_watersheds_path
valuation_table_path
watersheds_path E:\invest water yield\Processed\Pre2020\MAwsid_A.shp
workspace_dir E:\water4—3
2024-04-03 00:29:47,722 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 1 of 5 aligned: eto_2020.tif
2024-04-03 00:29:47,738 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 2 of 5 aligned: precip_2020.tif
2024-04-03 00:29:47,754 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 3 of 5 aligned: depth_to_root_rest_layer_2020.tif
2024-04-03 00:29:47,769 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 4 of 5 aligned: pawc_2020.tif
2024-04-03 00:29:47,784 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 5 of 5 aligned: clipped_lulc_2020.tif
2024-04-03 00:29:47,784 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1154) INFO aligned all 5 rasters.
2024-04-03 00:29:47,802 (natcap.invest.annual_water_yield) annual_water_yield.execute(696) INFO Reclassifying temp_Kc raster
2024-04-03 00:29:47,810 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2024-04-03 00:29:47,810 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-1 (stats_worker), started daemon 22428)>
2024-04-03 00:29:47,813 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2024-04-03 00:29:47,813 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2024-04-03 00:29:47,819 (natcap.invest.annual_water_yield) annual_water_yield.execute(708) INFO Reclassifying tmp_root raster
2024-04-03 00:29:47,825 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2024-04-03 00:29:47,825 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-2 (stats_worker), started daemon 28956)>
2024-04-03 00:29:47,828 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2024-04-03 00:29:47,828 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2024-04-03 00:29:47,834 (natcap.invest.annual_water_yield) annual_water_yield.execute(722) INFO Reclassifying tmp_veg raster
2024-04-03 00:29:47,842 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2024-04-03 00:29:47,842 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-3 (stats_worker), started daemon 30184)>
2024-04-03 00:29:47,844 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2024-04-03 00:29:47,844 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2024-04-03 00:29:47,850 (natcap.invest.annual_water_yield) annual_water_yield.execute(735) INFO Calculate PET from Ref Evap times Kc
2024-04-03 00:29:47,863 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2024-04-03 00:29:47,863 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-4 (stats_worker), started daemon 43296)>
2024-04-03 00:29:47,865 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2024-04-03 00:29:47,865 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2024-04-03 00:29:47,887 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2024-04-03 00:29:47,888 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-5 (stats_worker), started daemon 38648)>
2024-04-03 00:29:47,895 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2024-04-03 00:29:47,895 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2024-04-03 00:29:47,904 (natcap.invest.annual_water_yield) annual_water_yield.execute(766) INFO Performing wyield operation
2024-04-03 00:29:47,918 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2024-04-03 00:29:47,918 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-6 (stats_worker), started daemon 40228)>
2024-04-03 00:29:47,920 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2024-04-03 00:29:47,922 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2024-04-03 00:29:47,943 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2024-04-03 00:29:47,943 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-7 (stats_worker), started daemon 9828)>
2024-04-03 00:29:47,945 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2024-04-03 00:29:47,945 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2024-04-03 00:29:47,961 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field HYBAS_ID of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:47,961 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field NEXT_SINK of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:47,961 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field MAIN_BAS of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:47,961 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 38344733252.0999985 of field Shape_Area of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,029 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection
2024-04-03 00:29:48,030 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg
2024-04-03 00:29:48,042 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector
2024-04-03 00:29:48,049 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets
2024-04-03 00:29:48,049 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list
2024-04-03 00:29:48,050 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index
2024-04-03 00:29:48,050 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_wyield_2020.shp
2024-04-03 00:29:48,050 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup
2024-04-03 00:29:48,051 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_wyield_2020.shp
2024-04-03 00:29:48,051 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_wyield_2020.shp
2024-04-03 00:29:48,055 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1
2024-04-03 00:29:48,056 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 1
2024-04-03 00:29:48,058 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_wyield_2020.shp
2024-04-03 00:29:48,099 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection
2024-04-03 00:29:48,101 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg
2024-04-03 00:29:48,114 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector
2024-04-03 00:29:48,119 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets
2024-04-03 00:29:48,120 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list
2024-04-03 00:29:48,121 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index
2024-04-03 00:29:48,121 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_wyield_2020.shp
2024-04-03 00:29:48,121 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup
2024-04-03 00:29:48,121 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_wyield_2020.shp
2024-04-03 00:29:48,121 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_wyield_2020.shp
2024-04-03 00:29:48,125 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1
2024-04-03 00:29:48,128 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 1
2024-04-03 00:29:48,130 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_wyield_2020.shp
2024-04-03 00:29:48,173 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection
2024-04-03 00:29:48,175 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg
2024-04-03 00:29:48,190 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector
2024-04-03 00:29:48,196 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets
2024-04-03 00:29:48,196 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list
2024-04-03 00:29:48,197 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index
2024-04-03 00:29:48,197 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_wyield_2020.shp
2024-04-03 00:29:48,198 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup
2024-04-03 00:29:48,198 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_wyield_2020.shp
2024-04-03 00:29:48,198 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_wyield_2020.shp
2024-04-03 00:29:48,202 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1
2024-04-03 00:29:48,204 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 1
2024-04-03 00:29:48,206 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_wyield_2020.shp
2024-04-03 00:29:48,253 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection
2024-04-03 00:29:48,255 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg
2024-04-03 00:29:48,267 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector
2024-04-03 00:29:48,272 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets
2024-04-03 00:29:48,273 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list
2024-04-03 00:29:48,273 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index
2024-04-03 00:29:48,273 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_wyield_2020.shp
2024-04-03 00:29:48,273 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup
2024-04-03 00:29:48,273 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_wyield_2020.shp
2024-04-03 00:29:48,273 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_wyield_2020.shp
2024-04-03 00:29:48,278 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1
2024-04-03 00:29:48,280 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 1
2024-04-03 00:29:48,282 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_wyield_2020.shp
2024-04-03 00:29:48,292 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field HYBAS_ID of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,293 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field NEXT_SINK of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,293 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field MAIN_BAS of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,293 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 38344733252.0999985 of field Shape_Area of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,295 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field HYBAS_ID of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,295 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field NEXT_SINK of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,295 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field MAIN_BAS of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,295 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 38344733252.0999985 of field Shape_Area of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,296 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field HYBAS_ID of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,296 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field NEXT_SINK of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,296 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field MAIN_BAS of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,296 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 38344733252.0999985 of field Shape_Area of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,297 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field HYBAS_ID of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,297 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field NEXT_SINK of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,297 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field MAIN_BAS of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,297 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 38344733252.0999985 of field Shape_Area of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,297 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field HYBAS_ID of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,298 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field NEXT_SINK of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,298 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 4060050290 of field MAIN_BAS of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,298 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 38344733252.0999985 of field Shape_Area of feature 0 not successfully written. Possibly due to too larger number with respect to field width
2024-04-03 00:29:48,309 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 0.66s
2024-04-03 00:29:48,309 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished