python.exe .\seasonalwateryield.py 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(392) INFO prepare and test inputs for common errors 03/12/2021 07:42:32 (taskgraph.Task) Task._execution_monitor(824) INFO taskgraph execution status: tasks added: 0 tasks complete: 0 (0.0%) tasks waiting for a free worker: 0 (qsize: 0) tasks executing (0): graph is open 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(458) INFO Building file registry 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(464) INFO Checking that the AOI is not the output aggregate vector 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(473) INFO Aligning and clipping dataset list 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(575) INFO quick flow 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(589) INFO loading number of monthly events 03/12/2021 07:42:32 (taskgraph.Task) Task.add_task(726) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (n_events for month 0 (5)) as a dependent task to this one (n_events for month 1 (6)). 03/12/2021 07:42:32 (taskgraph.Task) Task.add_task(726) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (n_events for month 1 (6)) as a dependent task to this one (n_events for month 8 (12)). 03/12/2021 07:42:32 (taskgraph.Task) Task.add_task(726) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (n_events for month 3 (7)) as a dependent task to this one (n_events for month 10 (13)). 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(655) INFO calculate quick flow for month 1 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(655) INFO calculate quick flow for month 2 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(655) INFO calculate quick flow for month 3 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(655) INFO calculate quick flow for month 4 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(655) INFO calculate quick flow for month 5 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(655) INFO calculate quick flow for month 6 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(655) INFO calculate quick flow for month 7 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(655) INFO calculate quick flow for month 8 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(655) INFO calculate quick flow for month 9 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(655) INFO calculate quick flow for month 10 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(655) INFO calculate quick flow for month 11 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(655) INFO calculate quick flow for month 12 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(682) INFO calculate local recharge 03/12/2021 07:42:32 (taskgraph.Task) Task.add_task(726) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (classify kc month 0 (29)) as a dependent task to this one (classify kc month 1 (30)). 03/12/2021 07:42:32 (taskgraph.Task) Task.add_task(726) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (classify kc month 4 (32)) as a dependent task to this one (classify kc month 5 (33)). 03/12/2021 07:42:32 (taskgraph.Task) Task.add_task(726) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (classify kc month 6 (33)) as a dependent task to this one (classify kc month 7 (34)). 03/12/2021 07:42:32 (taskgraph.Task) Task.add_task(726) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (classify kc month 2 (30)) as a dependent task to this one (classify kc month 9 (35)). 03/12/2021 07:42:32 (taskgraph.Task) Task.add_task(726) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (classify kc month 3 (31)) as a dependent task to this one (classify kc month 10 (35)). 03/12/2021 07:42:32 (taskgraph.Task) Task.add_task(726) INFO A task was created that had an identical args signature sans target paths, but a different target_path_list of the same length. To avoid recomputation, dynamically adding previous Task (classify kc month 10 (35)) as a dependent task to this one (classify kc month 11 (35)). 03/12/2021 07:42:32 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._execute(752) INFO calculate L_sum 03/12/2021 07:42:33 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 1 of 28 aligned: prcp_a0_t01.tif 03/12/2021 07:42:33 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 2 of 28 aligned: prcp_a1_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 3 of 28 aligned: prcp_a2_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 4 of 28 aligned: prcp_a3_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 5 of 28 aligned: prcp_a4_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 6 of 28 aligned: prcp_a5_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 7 of 28 aligned: prcp_a6_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 8 of 28 aligned: prcp_a7_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 9 of 28 aligned: prcp_a8_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 10 of 28 aligned: prcp_a9_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 11 of 28 aligned: prcp_a10_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 12 of 28 aligned: prcp_a11_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 13 of 28 aligned: soil_group_aligned_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 14 of 28 aligned: et0_a0_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 15 of 28 aligned: et0_a1_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 16 of 28 aligned: et0_a2_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 17 of 28 aligned: et0_a3_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 18 of 28 aligned: et0_a4_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 19 of 28 aligned: et0_a5_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 20 of 28 aligned: et0_a6_t01.tif 03/12/2021 07:42:34 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 21 of 28 aligned: et0_a7_t01.tif 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 22 of 28 aligned: et0_a8_t01.tif 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 23 of 28 aligned: et0_a9_t01.tif 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 24 of 28 aligned: et0_a10_t01.tif 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 25 of 28 aligned: et0_a11_t01.tif 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 26 of 28 aligned: lulc_aligned_t01.tif 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 27 of 28 aligned: dem_aligned_t01.tif 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 28 of 28 aligned: cz_aligned_t01.tif 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(797) INFO aligned all 28 rasters. 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:35 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:36 (taskgraph.Task) Task.is_precalculated(1379) INFO not precalculated (classify kc month 1 (30)), Task hash exists, but there are these mismatches: Recorded path not in target path list d:\development\invest-sample-data\seasonal_water_yield\output\cache_dir\kc_0_t01.tif 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:36 (taskgraph.Task) Task.is_precalculated(1379) INFO not precalculated (classify kc month 5 (33)), Task hash exists, but there are these mismatches: Recorded path not in target path list d:\development\invest-sample-data\seasonal_water_yield\output\cache_dir\kc_4_t01.tif 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:36 (taskgraph.Task) Task.is_precalculated(1379) INFO not precalculated (n_events for month 1 (6)), Task hash exists, but there are these mismatches: Recorded path not in target path list d:\development\invest-sample-data\seasonal_water_yield\output\cache_dir\n_events0_t01.tif 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:36 (taskgraph.Task) Task.is_precalculated(1379) INFO not precalculated (classify kc month 7 (34)), Task hash exists, but there are these mismatches: Recorded path not in target path list d:\development\invest-sample-data\seasonal_water_yield\output\cache_dir\kc_6_t01.tif 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:36 (taskgraph.Task) Task.is_precalculated(1379) INFO not precalculated (n_events for month 8 (12)), Task hash exists, but there are these mismatches: Recorded path not in target path list d:\development\invest-sample-data\seasonal_water_yield\output\cache_dir\n_events1_t01.tif 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:36 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:37 (pygeoprocessing.routing.routing) pool.worker(125) INFO 100.0% complete 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:37 (taskgraph.Task) Task._execution_monitor(824) INFO taskgraph execution status: tasks added: 49 tasks complete: 19 (38.8%) tasks waiting for a free worker: 6 (qsize: 6) tasks executing (1): graph is closed flow dir mfd (2): executing for 0.20s 03/12/2021 07:42:37 (pygeoprocessing.routing.routing) pool.worker(125) INFO 100.0% complete 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:37 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:38 (taskgraph.Task) Task.is_precalculated(1379) INFO not precalculated (classify kc month 10 (35)), Task hash exists, but there are these mismatches: Recorded path not in target path list d:\development\invest-sample-data\seasonal_water_yield\output\cache_dir\kc_3_t01.tif 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:38 (taskgraph.Task) Task.is_precalculated(1379) INFO not precalculated (classify kc month 9 (35)), Task hash exists, but there are these mismatches: Recorded path not in target path list d:\development\invest-sample-data\seasonal_water_yield\output\cache_dir\kc_2_t01.tif 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:38 (pygeoprocessing.routing.routing) pool.worker(125) DEBUG creating target flow accum raster layer 03/12/2021 07:42:38 (pygeoprocessing.routing.routing) pool.worker(125) DEBUG creating visited raster layer 03/12/2021 07:42:38 (pygeoprocessing.routing.routing) pool.worker(125) DEBUG starting search 03/12/2021 07:42:38 (pygeoprocessing.routing.routing) pool.worker(125) INFO 100.0% complete 03/12/2021 07:42:38 (taskgraph.Task) Task.is_precalculated(1379) INFO not precalculated (classify kc month 11 (35)), Task hash exists, but there are these mismatches: Recorded path not in target path list d:\development\invest-sample-data\seasonal_water_yield\output\cache_dir\kc_10_t01.tif 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:38 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:39 (pygeoprocessing.routing.routing) pool.worker(125) INFO filter out incomplete divergent streams 03/12/2021 07:42:39 (pygeoprocessing.routing.routing) pool.worker(125) INFO 100.0% complete 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:39 (taskgraph.Task) Task.is_precalculated(1379) INFO not precalculated (n_events for month 10 (13)), Task hash exists, but there are these mismatches: Recorded path not in target path list d:\development\invest-sample-data\seasonal_water_yield\output\cache_dir\n_events3_t01.tif 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:39 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:40 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:41 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:42 (taskgraph.Task) Task._execution_monitor(824) INFO taskgraph execution status: tasks added: 49 tasks complete: 43 (87.8%) tasks waiting for a free worker: 0 (qsize: 0) tasks executing (1): graph is closed calculate QFi (28): executing for 0.01s 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:42 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:46 (pygeoprocessing.routing.routing) pool.worker(125) DEBUG creating target flow accum raster layer 03/12/2021 07:42:46 (pygeoprocessing.routing.routing) pool.worker(125) DEBUG creating visited raster layer 03/12/2021 07:42:46 (pygeoprocessing.routing.routing) pool.worker(125) DEBUG starting search 03/12/2021 07:42:46 (pygeoprocessing.routing.routing) pool.worker(125) INFO 100.0% complete 03/12/2021 07:42:46 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(391) INFO starting stats_worker 03/12/2021 07:42:46 (pygeoprocessing.geoprocessing_core) threading.run(892) DEBUG stats worker PID: 1956 03/12/2021 07:42:46 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(397) INFO started stats_worker 03/12/2021 07:42:46 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(477) INFO 100.0% complete 03/12/2021 07:42:46 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 03/12/2021 07:42:47 (taskgraph.Task) Task._execution_monitor(824) INFO taskgraph execution status: tasks added: 49 tasks complete: 47 (95.9%) tasks waiting for a free worker: 1 (qsize: 1) tasks executing (1): graph is closed calculate B_sum (39): executing for 0.53s 03/12/2021 07:42:48 (natcap.invest.seasonal_water_yield.seasonal_water_yield) seasonal_water_yield._aggregate_recharge(1129) WARNING D:\development\invest-sample-data\Seasonal_Water_Yield\output\aggregated_results_swy_t01.shp exists, deleting and writing new output 03/12/2021 07:42:48 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1184) DEBUG > 03/12/2021 07:42:48 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(793) INFO 1 of 1 aligned: clipped_raster9yxhxn0v.tif 03/12/2021 07:42:48 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(797) INFO aligned all 1 rasters. 03/12/2021 07:42:48 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1244) INFO build a lookup of aggregate field value to FID 03/12/2021 07:42:48 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1251) INFO creating disjoint polygon set 03/12/2021 07:42:48 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2183) INFO build shapely polygon list 03/12/2021 07:42:48 (taskgraph.Task) Task._task_executor(533) ERROR A taskgraph _task_executor failed on Task aggregate recharge (37). Terminating taskgraph. multiprocessing.pool.RemoteTraceback: """ Traceback (most recent call last): File "C:\Program Files\WindowsApps\PythonSoftwareFoundation.Python.3.9_3.9.752.0_x64__qbz5n2kfra8p0\lib\multiprocessing\pool.py", line 125, in worker result = (True, func(*args, **kwds)) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\natcap\invest\seasonal_water_yield\seasonal_water_yield.py", line 1147, in _aggregate_recharge aggregate_stats = pygeoprocessing.zonal_statistics( File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\pygeoprocessing\geoprocessing.py", line 1252, in zonal_statistics disjoint_fid_sets = calculate_disjoint_polygon_set( File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\pygeoprocessing\geoprocessing.py", line 2187, in calculate_disjoint_polygon_set bounding_box = shapely.prepared.prep(shapely.geometry.box(*bounding_box)) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\shapely\geometry\geo.py", line 59, in box return Polygon(coords) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\shapely\geometry\polygon.py", line 243, in __init__ ret = geos_polygon_from_py(shell, holes) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\shapely\geometry\polygon.py", line 509, in geos_polygon_from_py ret = geos_linearring_from_py(shell) File "shapely\speedups\_speedups.pyx", line 408, in shapely.speedups._speedups.geos_linearring_from_py ValueError: GEOSGeom_createLinearRing_r returned a NULL pointer """ The above exception was the direct cause of the following exception: Traceback (most recent call last): File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\taskgraph\Task.py", line 528, in _task_executor task._call() File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\taskgraph\Task.py", line 1231, in _call payload = result.get() File "C:\Program Files\WindowsApps\PythonSoftwareFoundation.Python.3.9_3.9.752.0_x64__qbz5n2kfra8p0\lib\multiprocessing\pool.py", line 771, in get raise self._value ValueError: GEOSGeom_createLinearRing_r returned a NULL pointer 03/12/2021 07:42:48 (taskgraph.Task) Task.join(875) ERROR Exception raised when joining task Task object 1743790913376: {'exception_object': ValueError('GEOSGeom_createLinearRing_r returned a NULL pointer'), 'ignore_directories': True, 'ignore_path_list': [], 'priority': 0, 'self._reexecution_info': {'args_clean': ['d:\\development\\invest-sample-data\\seasonal_water_yield\\watershed_gura.shp', 'd:\\development\\invest-sample-data\\seasonal_water_yield\\output\\l_t01.tif', 'd:\\development\\invest-sample-data\\seasonal_water_yield\\output\\vri_t01.tif', 'in_target_path_list'], 'file_stat_list': [('d:\\development\\invest-sample-data\\seasonal_water_yield\\watershed_gura.shp', 'sizetimestamp', '51084::1614259719.653430::d:\\development\\invest-sample-data\\seasonal_water_yield\\watershed_gura.shp'), ('d:\\development\\invest-sample-data\\seasonal_water_yield\\output\\l_t01.tif', 'sizetimestamp', '754186::1615531366.298816::d:\\development\\invest-sample-data\\seasonal_water_yield\\output\\l_t01.tif'), ('d:\\development\\invest-sample-data\\seasonal_water_yield\\output\\vri_t01.tif', 'sizetimestamp', '617017::1615531367.028808::d:\\development\\invest-sample-data\\seasonal_water_yield\\output\\vri_t01.tif')], 'func_name': '_aggregate_recharge', 'kwargs_clean': {}, 'other_arguments': [[None, None, None, 'd:\\development\\in_target_path_list'], {}], 'source_code_hash': 'd9762ddb96d23e2812811764a62eddf4ba58c502'}, 'self._result': None, 'target_path_list': ['d:\\development\\invest-sample-data\\seasonal_water_yield\\output\\aggregated_results_swy_t01.shp'], 'task_id_hash': '0ef9cb8a75583d3476768599ac4b6243c8496985', 'task_name': 'aggregate recharge (37)', 'task_reexecution_hash': 'a3977cf89fb83172527fbecf6fbb6c6b9876dd47'}. It's possible that this task did not cause the exception, rather another exception terminated the task_graph. Check the log to see if there are other exceptions. multiprocessing.pool.RemoteTraceback: """ Traceback (most recent call last): File "C:\Program Files\WindowsApps\PythonSoftwareFoundation.Python.3.9_3.9.752.0_x64__qbz5n2kfra8p0\lib\multiprocessing\pool.py", line 125, in worker result = (True, func(*args, **kwds)) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\natcap\invest\seasonal_water_yield\seasonal_water_yield.py", line 1147, in _aggregate_recharge aggregate_stats = pygeoprocessing.zonal_statistics( File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\pygeoprocessing\geoprocessing.py", line 1252, in zonal_statistics disjoint_fid_sets = calculate_disjoint_polygon_set( File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\pygeoprocessing\geoprocessing.py", line 2187, in calculate_disjoint_polygon_set bounding_box = shapely.prepared.prep(shapely.geometry.box(*bounding_box)) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\shapely\geometry\geo.py", line 59, in box return Polygon(coords) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\shapely\geometry\polygon.py", line 243, in __init__ ret = geos_polygon_from_py(shell, holes) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\shapely\geometry\polygon.py", line 509, in geos_polygon_from_py ret = geos_linearring_from_py(shell) File "shapely\speedups\_speedups.pyx", line 408, in shapely.speedups._speedups.geos_linearring_from_py ValueError: GEOSGeom_createLinearRing_r returned a NULL pointer """ The above exception was the direct cause of the following exception: Traceback (most recent call last): File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\taskgraph\Task.py", line 859, in join timedout = not task.join(timeout) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\taskgraph\Task.py", line 1399, in join raise self.exception_object File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\taskgraph\Task.py", line 528, in _task_executor task._call() File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\taskgraph\Task.py", line 1231, in _call payload = result.get() File "C:\Program Files\WindowsApps\PythonSoftwareFoundation.Python.3.9_3.9.752.0_x64__qbz5n2kfra8p0\lib\multiprocessing\pool.py", line 771, in get raise self._value ValueError: GEOSGeom_createLinearRing_r returned a NULL pointer multiprocessing.pool.RemoteTraceback: """ Traceback (most recent call last): File "C:\Program Files\WindowsApps\PythonSoftwareFoundation.Python.3.9_3.9.752.0_x64__qbz5n2kfra8p0\lib\multiprocessing\pool.py", line 125, in worker result = (True, func(*args, **kwds)) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\natcap\invest\seasonal_water_yield\seasonal_water_yield.py", line 1147, in _aggregate_recharge aggregate_stats = pygeoprocessing.zonal_statistics( File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\pygeoprocessing\geoprocessing.py", line 1252, in zonal_statistics disjoint_fid_sets = calculate_disjoint_polygon_set( File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\pygeoprocessing\geoprocessing.py", line 2187, in calculate_disjoint_polygon_set bounding_box = shapely.prepared.prep(shapely.geometry.box(*bounding_box)) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\shapely\geometry\geo.py", line 59, in box return Polygon(coords) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\shapely\geometry\polygon.py", line 243, in __init__ ret = geos_polygon_from_py(shell, holes) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\shapely\geometry\polygon.py", line 509, in geos_polygon_from_py ret = geos_linearring_from_py(shell) File "shapely\speedups\_speedups.pyx", line 408, in shapely.speedups._speedups.geos_linearring_from_py ValueError: GEOSGeom_createLinearRing_r returned a NULL pointer """ The above exception was the direct cause of the following exception: Traceback (most recent call last): File "D:\development\seasonalwateryield.py", line 48, in natcap.invest.seasonal_water_yield.seasonal_water_yield.execute(args) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\natcap\invest\seasonal_water_yield\seasonal_water_yield.py", line 377, in execute _execute(args) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\natcap\invest\seasonal_water_yield\seasonal_water_yield.py", line 786, in _execute task_graph.join() File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\taskgraph\Task.py", line 859, in join timedout = not task.join(timeout) File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\taskgraph\Task.py", line 1399, in join raise self.exception_object File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\taskgraph\Task.py", line 528, in _task_executor task._call() File "C:\Users\hydrologis\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.9_qbz5n2kfra8p0\LocalCache\local-packages\Python39\site-packages\taskgraph\Task.py", line 1231, in _call payload = result.get() File "C:\Program Files\WindowsApps\PythonSoftwareFoundation.Python.3.9_3.9.752.0_x64__qbz5n2kfra8p0\lib\multiprocessing\pool.py", line 771, in get raise self._value ValueError: GEOSGeom_createLinearRing_r returned a NULL pointer PS D:\development>