2023-03-30 20:09:39,705 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [C:\Users\我的小新\Desktop\论文实验记录\生态系统服务\水源涵养\result\1990\InVEST-natcap.invest.annual_water_yield-log-2023-03-30--20_09_39.txt] 2023-03-30 20:09:39,706 (__main__) cli.main(475) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.annual_water_yield 3.12.1: biophysical_table_path C:\Users\我的小新\Desktop\论文实验记录\生态系统服务\水源涵养\biophysical_table_gura.csv demand_table_path depth_to_root_rest_layer_path C:\Users\我的小新\Desktop\论文实验记录\生态系统服务\水源涵养\根系深度数据\cut_after\genxi_mmA.tif eto_path C:\Users\我的小新\Desktop\论文实验记录\生态系统服务\水源涵养\蒸发量数据\数据预处理-年均数据\年均蒸散量\mean_19901_A.tif lulc_path C:\Users\我的小新\Desktop\论文实验记录\内蒙古\重采样_1KM\1990_1km1.tif n_workers -1 pawc_path C:\Users\我的小新\Desktop\论文实验记录\生态系统服务\水源涵养\Plant Available Water Content（PAWC）\土壤数据——北京大学地理数据平台——HWSD的中国土壤数据集\pawcfile.tif\PAWC_A.tif precipitation_path C:\Users\我的小新\Desktop\论文实验记录\生态系统服务\水源涵养\降水数据\2、数据预处理—年度降水数据\算平均数据\premean_19901_A.tif results_suffix 1990 seasonality_constant 4.6 sub_watersheds_path valuation_table_path watersheds_path C:\Users\我的小新\Desktop\论文实验记录\生态系统服务\水源涵养\水域\nm_watersheds.shp workspace_dir C:\Users\我的小新\Desktop\论文实验记录\生态系统服务\水源涵养\result\1990 2023-03-30 20:09:39,706 (natcap.invest.annual_water_yield) annual_water_yield.execute(325) INFO Validating arguments 2023-03-30 20:09:39,706 (natcap.invest.validation) validation._wrapped_validate_func(1030) INFO Starting whole-model validation with ARGS_SPEC 2023-03-30 20:09:41,632 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 1 of 5 aligned: eto_1990.tif 2023-03-30 20:09:41,982 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 2 of 5 aligned: precip_1990.tif 2023-03-30 20:09:42,329 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 3 of 5 aligned: depth_to_root_rest_layer_1990.tif 2023-03-30 20:09:42,607 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 4 of 5 aligned: pawc_1990.tif 2023-03-30 20:09:42,803 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 5 of 5 aligned: clipped_lulc_1990.tif 2023-03-30 20:09:42,803 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(829) INFO aligned all 5 rasters. 2023-03-30 20:09:42,858 (natcap.invest.annual_water_yield) annual_water_yield.execute(509) INFO Reclassifying temp_Kc raster 2023-03-30 20:09:42,885 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:09:42,887 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:09:43,034 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:09:43,035 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:09:43,128 (natcap.invest.annual_water_yield) annual_water_yield.execute(521) INFO Reclassifying tmp_root raster 2023-03-30 20:09:43,151 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:09:43,152 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:09:43,323 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:09:43,323 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:09:43,421 (natcap.invest.annual_water_yield) annual_water_yield.execute(535) INFO Reclassifying tmp_veg raster 2023-03-30 20:09:43,451 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:09:43,452 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:09:43,633 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:09:43,633 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:09:43,727 (natcap.invest.annual_water_yield) annual_water_yield.execute(548) INFO Calculate PET from Ref Evap times Kc 2023-03-30 20:09:43,764 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:09:43,765 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:09:44,023 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:09:44,024 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:09:44,347 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:09:44,348 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:09:45,163 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:09:45,164 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:09:45,307 (natcap.invest.annual_water_yield) annual_water_yield.execute(579) INFO Performing wyield operation 2023-03-30 20:09:45,333 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:09:45,334 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:09:45,535 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:09:45,535 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:09:45,688 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:09:45,689 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:09:45,878 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:09:45,878 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:09:46,858 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-30 20:09:46,858 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(829) INFO aligned all 1 rasters. 2023-03-30 20:09:47,256 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1310) INFO build a lookup of aggregate field value to FID 2023-03-30 20:09:47,261 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1317) INFO creating disjoint polygon set 2023-03-30 20:09:47,263 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2314) INFO build shapely polygon list 2023-03-30 20:09:47,563 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2339) INFO build shapely rtree index 2023-03-30 20:09:47,579 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2352) INFO poly feature lookup 100.0% complete on watershed_results_wyield_1990.shp 2023-03-30 20:09:47,579 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2356) INFO build poly intersection lookup 2023-03-30 20:09:47,614 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2384) INFO poly intersection feature lookup 100.0% complete on watershed_results_wyield_1990.shp 2023-03-30 20:09:47,615 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2418) INFO maximal subset build 100.0% complete on watershed_results_wyield_1990.shp 2023-03-30 20:09:47,630 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO processing 3 disjoint polygon sets 2023-03-30 20:09:47,634 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 1 of 3 100.0% processed on watershed_results_wyield_1990.shp 2023-03-30 20:09:47,649 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 1 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:47,879 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 1 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:47,986 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 2 of 3 100.0% processed on watershed_results_wyield_1990.shp 2023-03-30 20:09:47,990 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 2 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:48,116 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 2 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:48,151 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 3 of 3 100.0% processed on watershed_results_wyield_1990.shp 2023-03-30 20:09:48,157 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 3 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:48,299 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 3 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:48,350 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1546) INFO all done processing polygon sets for watershed_results_wyield_1990.shp 2023-03-30 20:09:48,951 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-30 20:09:48,953 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(829) INFO aligned all 1 rasters. 2023-03-30 20:09:49,231 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1310) INFO build a lookup of aggregate field value to FID 2023-03-30 20:09:49,236 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1317) INFO creating disjoint polygon set 2023-03-30 20:09:49,240 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2314) INFO build shapely polygon list 2023-03-30 20:09:49,513 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2339) INFO build shapely rtree index 2023-03-30 20:09:49,515 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2352) INFO poly feature lookup 100.0% complete on watershed_results_wyield_1990.shp 2023-03-30 20:09:49,515 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2356) INFO build poly intersection lookup 2023-03-30 20:09:49,535 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2384) INFO poly intersection feature lookup 100.0% complete on watershed_results_wyield_1990.shp 2023-03-30 20:09:49,536 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2418) INFO maximal subset build 100.0% complete on watershed_results_wyield_1990.shp 2023-03-30 20:09:49,550 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO processing 3 disjoint polygon sets 2023-03-30 20:09:49,555 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 1 of 3 100.0% processed on watershed_results_wyield_1990.shp 2023-03-30 20:09:49,568 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 1 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:49,781 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 1 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:49,880 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 2 of 3 100.0% processed on watershed_results_wyield_1990.shp 2023-03-30 20:09:49,882 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 2 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:49,999 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 2 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:50,033 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 3 of 3 100.0% processed on watershed_results_wyield_1990.shp 2023-03-30 20:09:50,040 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 3 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:50,173 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 3 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:50,215 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1546) INFO all done processing polygon sets for watershed_results_wyield_1990.shp 2023-03-30 20:09:50,774 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-30 20:09:50,774 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(829) INFO aligned all 1 rasters. 2023-03-30 20:09:51,017 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1310) INFO build a lookup of aggregate field value to FID 2023-03-30 20:09:51,023 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1317) INFO creating disjoint polygon set 2023-03-30 20:09:51,024 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2314) INFO build shapely polygon list 2023-03-30 20:09:51,273 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2339) INFO build shapely rtree index 2023-03-30 20:09:51,274 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2352) INFO poly feature lookup 100.0% complete on watershed_results_wyield_1990.shp 2023-03-30 20:09:51,274 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2356) INFO build poly intersection lookup 2023-03-30 20:09:51,293 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2384) INFO poly intersection feature lookup 100.0% complete on watershed_results_wyield_1990.shp 2023-03-30 20:09:51,293 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2418) INFO maximal subset build 100.0% complete on watershed_results_wyield_1990.shp 2023-03-30 20:09:51,309 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO processing 3 disjoint polygon sets 2023-03-30 20:09:51,311 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 1 of 3 100.0% processed on watershed_results_wyield_1990.shp 2023-03-30 20:09:51,326 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 1 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:51,541 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 1 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:51,650 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 2 of 3 100.0% processed on watershed_results_wyield_1990.shp 2023-03-30 20:09:51,654 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 2 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:51,782 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 2 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:51,816 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 3 of 3 100.0% processed on watershed_results_wyield_1990.shp 2023-03-30 20:09:51,823 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 3 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:51,962 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 3 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:52,010 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1546) INFO all done processing polygon sets for watershed_results_wyield_1990.shp 2023-03-30 20:09:52,594 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-30 20:09:52,595 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(829) INFO aligned all 1 rasters. 2023-03-30 20:09:52,849 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1310) INFO build a lookup of aggregate field value to FID 2023-03-30 20:09:52,853 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1317) INFO creating disjoint polygon set 2023-03-30 20:09:52,854 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2314) INFO build shapely polygon list 2023-03-30 20:09:53,099 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2339) INFO build shapely rtree index 2023-03-30 20:09:53,101 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2352) INFO poly feature lookup 100.0% complete on watershed_results_wyield_1990.shp 2023-03-30 20:09:53,102 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2356) INFO build poly intersection lookup 2023-03-30 20:09:53,119 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2384) INFO poly intersection feature lookup 100.0% complete on watershed_results_wyield_1990.shp 2023-03-30 20:09:53,120 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2418) INFO maximal subset build 100.0% complete on watershed_results_wyield_1990.shp 2023-03-30 20:09:53,132 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO processing 3 disjoint polygon sets 2023-03-30 20:09:53,138 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 1 of 3 100.0% processed on watershed_results_wyield_1990.shp 2023-03-30 20:09:53,151 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 1 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:53,361 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 1 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:53,451 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 2 of 3 100.0% processed on watershed_results_wyield_1990.shp 2023-03-30 20:09:53,455 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 2 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:53,572 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 2 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:53,601 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 3 of 3 100.0% processed on watershed_results_wyield_1990.shp 2023-03-30 20:09:53,604 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 3 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:53,745 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 3 of 3 watershed_results_wyield_1990.shp 2023-03-30 20:09:53,798 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1546) INFO all done processing polygon sets for watershed_results_wyield_1990.shp 2023-03-30 20:09:55,346 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 15.64s 2023-03-30 20:09:55,347 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished