2021-12-03 11:21:59,201 (natcap.invest.utils) utils.prepare_workspace(158) INFO Writing log messages to C:\Users\felix\Documents\annual_water_yield_workspace\InVEST-Annual-Water-Yield-log-2021-12-03--11_21_59.txt 2021-12-03 11:21:59,223 (natcap.invest.ui.model) model._logged_target(1644) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.annual_water_yield 3.9.2: biophysical_table_path C:/Users/felix/OneDrive/Documents/Uni/Masterarbeit/inVEST/Annual Water Yield/input Data/LU classes in floodplains csv.csv depth_to_root_rest_layer_path C:/Users/felix/OneDrive/Documents/Uni/Masterarbeit/inVEST/Annual Water Yield/input Data/BDRICM M 250m clipped reprojected.tif do_scarcity_and_valuation False eto_path C:/Users/felix/OneDrive/Documents/Uni/Masterarbeit/inVEST/Annual Water Yield/input Data/Zuschnitt Evapotranspiration reprojected.tif lulc_path C:/Users/felix/OneDrive/Documents/Uni/Masterarbeit/inVEST/Annual Water Yield/input Data/LULC try nodata.tif n_workers -1 pawc_path C:/Users/felix/OneDrive/Documents/Uni/Masterarbeit/inVEST/Annual Water Yield/input Data/weighted average clipped reprojected.tif precipitation_path C:/Users/felix/OneDrive/Documents/Uni/Masterarbeit/inVEST/Annual Water Yield/input Data/Precipitation Zuschnitt reprojected.tif results_suffix seasonality_constant 5 sub_watersheds_path watersheds_path C:/Users/felix/OneDrive/Documents/Uni/Masterarbeit/inVEST/Annual Water Yield/input Data/watersheds_renamed.gpkg workspace_dir C:\Users\felix\Documents\annual_water_yield_workspace 2021-12-03 11:21:59,271 (natcap.invest.annual_water_yield) annual_water_yield.execute(269) INFO Validating arguments 2021-12-03 11:21:59,282 (natcap.invest.validation) validation._wrapped_validate_func(920) DEBUG Using ARG_SPEC for validation 2021-12-03 11:21:59,291 (natcap.invest.validation) validation._wrapped_validate_func(924) INFO Starting whole-model validation with ARGS_SPEC 2021-12-03 11:21:59,297 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-12-03 11:21:59,319 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-12-03 11:21:59,319 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-12-03 11:21:59,321 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-12-03 11:21:59,321 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-12-03 11:21:59,342 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-12-03 11:21:59,342 (natcap.invest.validation) validation.validate(802) DEBUG Provided key do_scarcity_and_valuation does not exist in ARGS_SPEC 2021-12-03 11:21:59,342 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-12-03 11:21:59,363 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-12-03 11:21:59,363 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-12-03 11:21:59,382 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-12-03 11:21:59,382 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-12-03 11:21:59,397 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-12-03 11:21:59,398 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-12-03 11:21:59,403 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-12-03 11:21:59,403 (natcap.invest.validation) validation.timeout(633) DEBUG Starting file checking thread with timeout=5 2021-12-03 11:21:59,425 (natcap.invest.validation) validation.timeout(645) DEBUG File checking thread completed. 2021-12-03 11:21:59,472 (natcap.invest.validation) validation._wrapped_validate_func(969) DEBUG Validation warnings: [] 2021-12-03 11:22:00,304 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 5 aligned: eto.tif 2021-12-03 11:22:01,032 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 2 of 5 aligned: precip.tif 2021-12-03 11:22:01,455 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 3 of 5 aligned: depth_to_root_rest_layer.tif 2021-12-03 11:22:03,070 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 4 of 5 aligned: pawc.tif 2021-12-03 11:22:03,768 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 5 of 5 aligned: clipped_lulc.tif 2021-12-03 11:22:03,769 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 5 rasters. 2021-12-03 11:22:03,804 (natcap.invest.annual_water_yield) annual_water_yield.execute(405) DEBUG bio_lucodes {512, 131, 132, 133, 9.0, 141, 142, 411, 412, 311, 312, 313, 321, 322, 324, 331, 332, 333, 335, 211, 221, 222, 231, 111, 112, 242, 243, 244, 121, 122, 123, 124, 511} 2021-12-03 11:22:03,804 (natcap.invest.annual_water_yield) annual_water_yield.execute(453) INFO Reclassifying temp_Kc raster 2021-12-03 11:22:03,831 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-12-03 11:22:03,833 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 2848 2021-12-03 11:22:03,833 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-12-03 11:22:04,306 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-12-03 11:22:04,310 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-12-03 11:22:04,554 (natcap.invest.annual_water_yield) annual_water_yield.execute(465) INFO Reclassifying tmp_root raster 2021-12-03 11:22:04,573 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-12-03 11:22:04,574 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 2848 2021-12-03 11:22:04,574 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-12-03 11:22:05,008 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-12-03 11:22:05,008 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-12-03 11:22:05,319 (natcap.invest.annual_water_yield) annual_water_yield.execute(479) INFO Reclassifying tmp_veg raster 2021-12-03 11:22:05,335 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-12-03 11:22:05,336 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 2848 2021-12-03 11:22:05,336 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-12-03 11:22:05,748 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-12-03 11:22:05,749 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-12-03 11:22:05,996 (natcap.invest.annual_water_yield) annual_water_yield.execute(492) INFO Calculate PET from Ref Evap times Kc 2021-12-03 11:22:06,029 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-12-03 11:22:06,030 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 2848 2021-12-03 11:22:06,031 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-12-03 11:22:06,551 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-12-03 11:22:06,552 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-12-03 11:22:06,936 (natcap.invest.annual_water_yield) annual_water_yield.execute(510) DEBUG Performing fractp operation 2021-12-03 11:22:06,974 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-12-03 11:22:06,975 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 2848 2021-12-03 11:22:06,975 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-12-03 11:22:08,809 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-12-03 11:22:08,809 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-12-03 11:22:09,217 (natcap.invest.annual_water_yield) annual_water_yield.execute(523) INFO Performing wyield operation 2021-12-03 11:22:09,235 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-12-03 11:22:09,236 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 2848 2021-12-03 11:22:09,236 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-12-03 11:22:09,756 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-12-03 11:22:09,756 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-12-03 11:22:10,128 (natcap.invest.annual_water_yield) annual_water_yield.execute(536) DEBUG Performing aet operation 2021-12-03 11:22:10,145 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2021-12-03 11:22:10,146 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 2848 2021-12-03 11:22:10,146 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2021-12-03 11:22:10,639 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2021-12-03 11:22:10,639 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2021-12-03 11:22:11,898 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2021-12-03 11:22:11,898 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2021-12-03 11:22:11,904 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2021-12-03 11:22:11,938 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2021-12-03 11:22:11,942 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2021-12-03 11:22:11,974 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2021-12-03 11:22:11,993 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watersheds_renamed.gpkg 2021-12-03 11:22:11,994 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2021-12-03 11:22:12,289 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watersheds_renamed.gpkg 2021-12-03 11:22:12,290 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watersheds_renamed.gpkg 2021-12-03 11:22:12,304 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 6 disjoint polygon sets 2021-12-03 11:22:12,309 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:12,333 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:12,714 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:13,009 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:13,014 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:13,575 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:13,756 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:13,761 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:14,032 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:14,163 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:14,168 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:14,451 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:14,588 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:14,593 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:15,456 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:15,682 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:15,687 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:16,011 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:16,166 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 14 2021-12-03 11:22:16,167 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 4.3686516e+05 1.0000833e+02 0.0000000e+00 5.5464440e+06 0.0000000e+00 -1.0000833e+02] for ('C:\\Users\\felix\\Documents\\annual_water_yield_workspace\\intermediate\\precip.tif', 1) 2021-12-03 11:22:16,168 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 14 2021-12-03 11:22:16,168 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watersheds_renamed.gpkg 2021-12-03 11:22:17,069 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2021-12-03 11:22:17,070 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2021-12-03 11:22:17,071 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2021-12-03 11:22:17,077 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2021-12-03 11:22:17,082 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2021-12-03 11:22:17,113 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2021-12-03 11:22:17,119 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watersheds_renamed.gpkg 2021-12-03 11:22:17,120 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2021-12-03 11:22:17,417 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watersheds_renamed.gpkg 2021-12-03 11:22:17,418 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watersheds_renamed.gpkg 2021-12-03 11:22:17,431 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 6 disjoint polygon sets 2021-12-03 11:22:17,435 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:17,442 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:17,843 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:18,175 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:18,181 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:18,831 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:19,107 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:19,114 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:19,493 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:19,689 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:19,696 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:20,074 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:20,304 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:20,311 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:21,240 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:21,492 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:21,499 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:21,933 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:22,095 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 14 2021-12-03 11:22:22,095 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 4.3686516e+05 1.0000833e+02 0.0000000e+00 5.5464440e+06 0.0000000e+00 -1.0000833e+02] for ('C:\\Users\\felix\\Documents\\annual_water_yield_workspace\\intermediate\\pet.tif', 1) 2021-12-03 11:22:22,096 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 14 2021-12-03 11:22:22,096 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watersheds_renamed.gpkg 2021-12-03 11:22:23,277 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2021-12-03 11:22:23,278 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2021-12-03 11:22:23,280 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2021-12-03 11:22:23,288 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2021-12-03 11:22:23,293 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2021-12-03 11:22:23,321 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2021-12-03 11:22:23,330 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watersheds_renamed.gpkg 2021-12-03 11:22:23,331 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2021-12-03 11:22:23,538 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watersheds_renamed.gpkg 2021-12-03 11:22:23,539 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watersheds_renamed.gpkg 2021-12-03 11:22:23,551 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 6 disjoint polygon sets 2021-12-03 11:22:23,553 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:23,560 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:23,998 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:24,424 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:24,430 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:25,116 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:25,363 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:25,370 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:25,742 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:25,948 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:25,955 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:26,328 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:26,474 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:26,481 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:27,526 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:27,729 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:27,736 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:28,173 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:28,378 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 14 2021-12-03 11:22:28,378 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 4.3686516e+05 1.0000833e+02 0.0000000e+00 5.5464440e+06 0.0000000e+00 -1.0000833e+02] for ('C:\\Users\\felix\\Documents\\annual_water_yield_workspace\\output\\per_pixel\\aet.tif', 1) 2021-12-03 11:22:28,379 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 14 2021-12-03 11:22:28,379 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watersheds_renamed.gpkg 2021-12-03 11:22:29,332 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2021-12-03 11:22:29,332 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2021-12-03 11:22:29,334 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2021-12-03 11:22:29,343 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2021-12-03 11:22:29,350 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2021-12-03 11:22:29,421 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2021-12-03 11:22:29,429 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watersheds_renamed.gpkg 2021-12-03 11:22:29,430 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2021-12-03 11:22:29,761 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watersheds_renamed.gpkg 2021-12-03 11:22:29,762 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watersheds_renamed.gpkg 2021-12-03 11:22:29,775 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 6 disjoint polygon sets 2021-12-03 11:22:29,779 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:29,786 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:30,242 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:30,554 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:30,562 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:31,155 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:31,325 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:31,331 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:31,645 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:31,847 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:31,854 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:32,232 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:32,413 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:32,419 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:33,291 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:33,492 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 6 100.0% processed on watersheds_renamed.gpkg 2021-12-03 11:22:33,498 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:33,806 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 6 watersheds_renamed.gpkg 2021-12-03 11:22:33,935 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 14 2021-12-03 11:22:33,935 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 4.3686516e+05 1.0000833e+02 0.0000000e+00 5.5464440e+06 0.0000000e+00 -1.0000833e+02] for ('C:\\Users\\felix\\Documents\\annual_water_yield_workspace\\output\\per_pixel\\wyield.tif', 1) 2021-12-03 11:22:33,936 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 14 2021-12-03 11:22:33,936 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watersheds_renamed.gpkg 2021-12-03 11:22:34,055 (taskgraph.Task) Task.add_task(706) ERROR Something went wrong when adding task create_ws_id_vector_output (12), terminating taskgraph. Traceback (most recent call last): File "c:\program files\python38\lib\site-packages\taskgraph\Task.py", line 674, in add_task new_task._call() File "c:\program files\python38\lib\site-packages\taskgraph\Task.py", line 1093, in _call payload = self._func(*self._args, **self._kwargs) File "c:\program files\python38\lib\site-packages\natcap\invest\annual_water_yield.py", line 675, in create_vector_output _add_zonal_stats_dict_to_shape( File "c:\program files\python38\lib\site-packages\natcap\invest\annual_water_yield.py", line 1096, in _add_zonal_stats_dict_to_shape if stats_map[feature_fid]['count'] > 0: KeyError: 0 2021-12-03 11:22:34,063 (natcap.invest.utils) utils.prepare_workspace(163) ERROR Exception while executing Annual-Water-Yield Traceback (most recent call last): File "c:\program files\python38\lib\site-packages\natcap\invest\utils.py", line 161, in prepare_workspace yield File "c:\program files\python38\lib\site-packages\natcap\invest\ui\model.py", line 1650, in _logged_target return self.target(args=args) File "c:\program files\python38\lib\site-packages\natcap\invest\annual_water_yield.py", line 597, in execute create_output_vector_task = graph.add_task( File "c:\program files\python38\lib\site-packages\taskgraph\Task.py", line 674, in add_task new_task._call() File "c:\program files\python38\lib\site-packages\taskgraph\Task.py", line 1093, in _call payload = self._func(*self._args, **self._kwargs) File "c:\program files\python38\lib\site-packages\natcap\invest\annual_water_yield.py", line 675, in create_vector_output _add_zonal_stats_dict_to_shape( File "c:\program files\python38\lib\site-packages\natcap\invest\annual_water_yield.py", line 1096, in _add_zonal_stats_dict_to_shape if stats_map[feature_fid]['count'] > 0: KeyError: 0 2021-12-03 11:22:34,067 (natcap.invest.utils) utils.prepare_workspace(166) INFO Elapsed time: 34.86s 2021-12-03 11:22:34,067 (natcap.invest.utils) utils.prepare_workspace(169) INFO Execution finished