2021-11-10 14:24:19,439 (natcap.invest.utils) utils.prepare_workspace(158) INFO Writing log messages to /Users/Sonali 1/Downloads/test/InVEST-Urban-Flood-Risk-Mitigation-log-2021-11-10--14_24_19.txt 2021-11-10 14:24:19,440 (natcap.invest.ui.model) model._logged_target(1644) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.urban_flood_risk_mitigation 3.9.2: aoi_watersheds_path /Users/Sonali 1/Library/Mobile Documents/com~apple~CloudDocs/NTU/PhD/Modules/ES7030/InVEST Assignment/InVEST inputs - ES3202/Watersheds (NCSg)/catchmentsMOD2.shp built_infrastructure_vector_path /Users/Sonali 1/Downloads/sg_flood_current_new.shp curve_number_table_path /Users/Sonali 1/Library/Mobile Documents/com~apple~CloudDocs/NTU/PhD/Modules/ES7030/InVEST Assignment/InVEST inputs - ES3202/Biophysical table_ES3202_upload.csv infrastructure_damage_loss_table_path /Users/Sonali 1/Library/Mobile Documents/com~apple~CloudDocs/NTU/PhD/Modules/ES7030/Final Project/damage_loss_current_unique.csv lulc_path /Users/Sonali 1/Library/Mobile Documents/com~apple~CloudDocs/NTU/PhD/Modules/ES7030/InVEST Assignment/InVEST inputs - ES3202/LULC.tif n_workers -1 rainfall_depth 105 results_suffix soils_hydrological_group_raster_path /Users/Sonali 1/Library/Mobile Documents/com~apple~CloudDocs/NTU/PhD/Modules/ES7030/InVEST Assignment/InVEST inputs - ES3202/HSoilGroup.tif workspace_dir /Users/Sonali 1/Downloads/test 2021-11-10 14:24:19,441 (natcap.invest.validation) validation._wrapped_validate_func(920) DEBUG Using ARG_SPEC for validation 2021-11-10 14:24:19,441 (natcap.invest.validation) validation._wrapped_validate_func(924) INFO Starting whole-model validation with ARGS_SPEC 2021-11-10 14:24:19,442 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-11-10 14:24:19,496 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-11-10 14:24:19,496 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-11-10 14:24:19,502 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-11-10 14:24:19,502 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-11-10 14:24:19,514 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-11-10 14:24:19,514 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-11-10 14:24:19,567 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-11-10 14:24:19,570 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-11-10 14:24:19,587 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-11-10 14:24:19,587 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-11-10 14:24:19,603 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-11-10 14:24:19,603 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-11-10 14:24:19,603 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-11-10 14:24:19,666 (natcap.invest.validation) validation._wrapped_validate_func(969) DEBUG Validation warnings: [] 2021-11-10 14:24:21,975 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 2 aligned: aligned_lulc.tif 2021-11-10 14:24:23,906 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 2 of 2 aligned: aligned_soils_hydrological_group.tif 2021-11-10 14:24:23,907 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 2 rasters. 2021-11-10 14:24:23,930 (natcap.invest.utils) utils.build_lookup_from_csv(565) WARNING Empty or NaN values were found in the table: /Users/Sonali 1/Library/Mobile Documents/com~apple~CloudDocs/NTU/PhD/Modules/ES7030/InVEST Assignment/InVEST inputs - ES3202/Biophysical table_ES3202_upload.csv. 2021-11-10 14:24:23,932 (natcap.invest.utils) utils.build_lookup_from_csv(572) WARNING Encountered an entirely blank row on line(s) [nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan]. Dropping rows from table. 2021-11-10 14:24:23,987 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-11-10 14:24:24,006 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 39071 2021-11-10 14:24:24,005 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-11-10 14:24:26,908 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-11-10 14:24:26,914 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-11-10 14:24:27,536 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-11-10 14:24:27,537 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 39071 2021-11-10 14:24:27,537 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-11-10 14:24:27,947 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-11-10 14:24:27,948 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-11-10 14:24:28,502 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-11-10 14:24:28,503 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 39071 2021-11-10 14:24:28,503 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-11-10 14:24:29,041 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-11-10 14:24:29,042 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-11-10 14:24:29,655 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-11-10 14:24:29,657 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 39071 2021-11-10 14:24:29,657 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-11-10 14:24:30,018 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-11-10 14:24:30,018 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-11-10 14:24:30,487 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-11-10 14:24:30,488 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 39071 2021-11-10 14:24:30,489 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-11-10 14:24:30,841 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-11-10 14:24:30,841 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-11-10 14:24:31,276 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-11-10 14:24:31,276 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 39071 2021-11-10 14:24:31,276 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-11-10 14:24:31,598 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-11-10 14:24:31,598 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-11-10 14:24:32,073 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1741) INFO starting reprojection 2021-11-10 14:24:32,077 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer catchmentsMOD2 of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2021-11-10 14:24:32,141 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1783) INFO reprojection 100.0% complete on reprojected_aoi.gpkg 2021-11-10 14:24:33,336 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2021-11-10 14:24:33,337 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2021-11-10 14:24:33,338 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2021-11-10 14:24:33,347 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2021-11-10 14:24:33,351 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2021-11-10 14:24:33,413 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2021-11-10 14:24:33,445 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on reprojected_aoi.gpkg 2021-11-10 14:24:33,445 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2021-11-10 14:24:33,767 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on reprojected_aoi.gpkg 2021-11-10 14:24:33,768 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on reprojected_aoi.gpkg 2021-11-10 14:24:33,782 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 7 disjoint polygon sets 2021-11-10 14:24:33,806 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:33,854 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:34,247 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:34,725 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:34,733 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:35,139 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:35,415 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:35,422 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:35,933 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:36,210 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:36,218 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:36,629 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:36,906 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:36,916 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:37,306 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:37,554 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:37,565 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:37,919 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:38,130 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 7 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:38,137 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 7 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:38,482 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 7 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:38,683 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 3 of 370 2021-11-10 14:24:38,683 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 3.4491953e+05 9.0000000e+00 0.0000000e+00 1.6259433e+05 0.0000000e+00 -9.0000000e+00] for ('/Users/Sonali 1/Downloads/test/intermediate_files/Q_m3.tif', 1) 2021-11-10 14:24:38,687 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 370 2021-11-10 14:24:38,688 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for reprojected_aoi.gpkg 2021-11-10 14:24:39,807 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2021-11-10 14:24:39,807 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2021-11-10 14:24:39,809 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2021-11-10 14:24:39,817 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2021-11-10 14:24:39,822 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2021-11-10 14:24:39,864 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2021-11-10 14:24:39,892 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on reprojected_aoi.gpkg 2021-11-10 14:24:39,892 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2021-11-10 14:24:40,211 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on reprojected_aoi.gpkg 2021-11-10 14:24:40,212 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on reprojected_aoi.gpkg 2021-11-10 14:24:40,222 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 7 disjoint polygon sets 2021-11-10 14:24:40,245 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:40,296 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:40,686 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:41,145 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:41,153 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:41,536 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:41,802 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:41,809 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:42,196 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:42,460 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:42,467 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:42,845 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:43,124 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:43,133 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:43,511 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:43,762 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:43,769 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:44,116 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:44,346 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 7 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:44,353 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 7 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:44,685 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 7 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:44,875 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 3 of 370 2021-11-10 14:24:44,876 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 3.4491953e+05 9.0000000e+00 0.0000000e+00 1.6259433e+05 0.0000000e+00 -9.0000000e+00] for ('/Users/Sonali 1/Downloads/test/Runoff_retention.tif', 1) 2021-11-10 14:24:44,878 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 370 2021-11-10 14:24:44,878 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for reprojected_aoi.gpkg 2021-11-10 14:24:45,957 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2021-11-10 14:24:45,959 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2021-11-10 14:24:45,961 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2021-11-10 14:24:45,968 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2021-11-10 14:24:45,975 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2021-11-10 14:24:46,013 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2021-11-10 14:24:46,040 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on reprojected_aoi.gpkg 2021-11-10 14:24:46,040 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2021-11-10 14:24:46,372 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on reprojected_aoi.gpkg 2021-11-10 14:24:46,374 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on reprojected_aoi.gpkg 2021-11-10 14:24:46,385 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 7 disjoint polygon sets 2021-11-10 14:24:46,403 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:46,456 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:46,834 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:47,288 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:47,296 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:47,681 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:47,942 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:47,949 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:48,341 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:48,621 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:48,628 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:49,010 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:49,283 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:49,292 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:49,665 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:49,897 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:49,904 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:50,248 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:50,467 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 7 of 7 100.0% processed on reprojected_aoi.gpkg 2021-11-10 14:24:50,477 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 7 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:50,820 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 7 of 7 reprojected_aoi.gpkg 2021-11-10 14:24:51,019 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 3 of 370 2021-11-10 14:24:51,020 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 3.4491953e+05 9.0000000e+00 0.0000000e+00 1.6259433e+05 0.0000000e+00 -9.0000000e+00] for ('/Users/Sonali 1/Downloads/test/Runoff_retention_m3.tif', 1) 2021-11-10 14:24:51,021 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 370 2021-11-10 14:24:51,021 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for reprojected_aoi.gpkg 2021-11-10 14:24:51,064 (osgeo) utils._log_gdal_errors(97) ERROR [errno 1] Wrong field type for fid 2021-11-10 14:24:51,066 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1741) INFO starting reprojection 2021-11-10 14:24:51,067 (osgeo) utils._log_gdal_errors(97) ERROR [errno 1] Invalid index : 45 2021-11-10 14:24:51,069 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer sg_flood_current_new of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2021-11-10 14:24:51,070 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1783) INFO reprojection 100.0% complete on structures_reprojected.gpkg 2021-11-10 14:24:51,143 (taskgraph.Task) Task.add_task(706) ERROR Something went wrong when adding task calculate damage to infrastructure in aoi (12), terminating taskgraph. Traceback (most recent call last): File "taskgraph/Task.py", line 674, in add_task File "taskgraph/Task.py", line 1093, in _call File "natcap/invest/urban_flood_risk_mitigation.py", line 611, in _calculate_damage_to_infrastructure_in_aoi KeyError: 3038 2021-11-10 14:24:51,146 (natcap.invest.utils) utils.prepare_workspace(163) ERROR Exception while executing Urban-Flood-Risk-Mitigation Traceback (most recent call last): File "natcap/invest/utils.py", line 161, in prepare_workspace File "natcap/invest/ui/model.py", line 1650, in _logged_target File "natcap/invest/urban_flood_risk_mitigation.py", line 379, in execute File "taskgraph/Task.py", line 674, in add_task File "taskgraph/Task.py", line 1093, in _call File "natcap/invest/urban_flood_risk_mitigation.py", line 611, in _calculate_damage_to_infrastructure_in_aoi KeyError: 3038 2021-11-10 14:24:51,146 (natcap.invest.utils) utils.prepare_workspace(166) INFO Elapsed time: 31.71s 2021-11-10 14:24:51,147 (natcap.invest.utils) utils.prepare_workspace(169) INFO Execution finished