2023-03-30 20:16:04,072 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [C:\Users\我的小新\Desktop\论文实验记录\生态系统服务\水源涵养\result\1995\InVEST-natcap.invest.annual_water_yield-log-2023-03-30--20_16_04.txt] 2023-03-30 20:16:04,073 (__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_19951_A.tif lulc_path C:\Users\我的小新\Desktop\论文实验记录\内蒙古\重采样_1KM\1995_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_19951_A.tif results_suffix 1995 seasonality_constant 4.6 sub_watersheds_path valuation_table_path watersheds_path C:\Users\我的小新\Desktop\论文实验记录\生态系统服务\水源涵养\水域\nm_watersheds.shp workspace_dir C:\Users\我的小新\Desktop\论文实验记录\生态系统服务\水源涵养\result\1995 2023-03-30 20:16:04,073 (natcap.invest.annual_water_yield) annual_water_yield.execute(325) INFO Validating arguments 2023-03-30 20:16:04,074 (natcap.invest.validation) validation._wrapped_validate_func(1030) INFO Starting whole-model validation with ARGS_SPEC 2023-03-30 20:16:06,181 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 1 of 5 aligned: eto_1995.tif 2023-03-30 20:16:06,663 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 2 of 5 aligned: precip_1995.tif 2023-03-30 20:16:07,082 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 3 of 5 aligned: depth_to_root_rest_layer_1995.tif 2023-03-30 20:16:07,414 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 4 of 5 aligned: pawc_1995.tif 2023-03-30 20:16:07,564 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 5 of 5 aligned: clipped_lulc_1995.tif 2023-03-30 20:16:07,564 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(829) INFO aligned all 5 rasters. 2023-03-30 20:16:07,612 (natcap.invest.annual_water_yield) annual_water_yield.execute(509) INFO Reclassifying temp_Kc raster 2023-03-30 20:16:07,645 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:16:07,646 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:16:07,794 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:16:07,794 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:16:07,887 (natcap.invest.annual_water_yield) annual_water_yield.execute(521) INFO Reclassifying tmp_root raster 2023-03-30 20:16:07,918 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:16:07,919 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:16:08,067 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:16:08,067 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:16:08,168 (natcap.invest.annual_water_yield) annual_water_yield.execute(535) INFO Reclassifying tmp_veg raster 2023-03-30 20:16:08,201 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:16:08,202 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:16:08,350 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:16:08,351 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:16:08,449 (natcap.invest.annual_water_yield) annual_water_yield.execute(548) INFO Calculate PET from Ref Evap times Kc 2023-03-30 20:16:08,480 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:16:08,481 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:16:08,693 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:16:08,693 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:16:09,010 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:16:09,011 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:16:09,847 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:16:09,847 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:16:10,000 (natcap.invest.annual_water_yield) annual_water_yield.execute(579) INFO Performing wyield operation 2023-03-30 20:16:10,029 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:16:10,030 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:16:10,251 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:16:10,251 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:16:10,416 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-30 20:16:10,417 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-30 20:16:10,605 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(483) INFO 100.0% complete 2023-03-30 20:16:10,605 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(486) INFO Waiting for raster stats worker result. 2023-03-30 20:16:11,629 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-30 20:16:11,629 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(829) INFO aligned all 1 rasters. 2023-03-30 20:16:11,880 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1310) INFO build a lookup of aggregate field value to FID 2023-03-30 20:16:11,884 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1317) INFO creating disjoint polygon set 2023-03-30 20:16:11,886 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2314) INFO build shapely polygon list 2023-03-30 20:16:12,160 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2339) INFO build shapely rtree index 2023-03-30 20:16:12,162 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2352) INFO poly feature lookup 100.0% complete on watershed_results_wyield_1995.shp 2023-03-30 20:16:12,162 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2356) INFO build poly intersection lookup 2023-03-30 20:16:12,182 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2384) INFO poly intersection feature lookup 100.0% complete on watershed_results_wyield_1995.shp 2023-03-30 20:16:12,183 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2418) INFO maximal subset build 100.0% complete on watershed_results_wyield_1995.shp 2023-03-30 20:16:12,195 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO processing 3 disjoint polygon sets 2023-03-30 20:16:12,198 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 1 of 3 100.0% processed on watershed_results_wyield_1995.shp 2023-03-30 20:16:12,212 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 1 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:12,518 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 1 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:12,657 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 2 of 3 100.0% processed on watershed_results_wyield_1995.shp 2023-03-30 20:16:12,659 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 2 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:12,789 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 2 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:12,821 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 3 of 3 100.0% processed on watershed_results_wyield_1995.shp 2023-03-30 20:16:12,827 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 3 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:12,967 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 3 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:13,015 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1546) INFO all done processing polygon sets for watershed_results_wyield_1995.shp 2023-03-30 20:16:13,637 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-30 20:16:13,638 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(829) INFO aligned all 1 rasters. 2023-03-30 20:16:13,889 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1310) INFO build a lookup of aggregate field value to FID 2023-03-30 20:16:13,894 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1317) INFO creating disjoint polygon set 2023-03-30 20:16:13,897 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2314) INFO build shapely polygon list 2023-03-30 20:16:14,183 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2339) INFO build shapely rtree index 2023-03-30 20:16:14,184 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2352) INFO poly feature lookup 100.0% complete on watershed_results_wyield_1995.shp 2023-03-30 20:16:14,185 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2356) INFO build poly intersection lookup 2023-03-30 20:16:14,212 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2384) INFO poly intersection feature lookup 100.0% complete on watershed_results_wyield_1995.shp 2023-03-30 20:16:14,212 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2418) INFO maximal subset build 100.0% complete on watershed_results_wyield_1995.shp 2023-03-30 20:16:14,231 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO processing 3 disjoint polygon sets 2023-03-30 20:16:14,235 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 1 of 3 100.0% processed on watershed_results_wyield_1995.shp 2023-03-30 20:16:14,252 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 1 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:14,511 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 1 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:14,627 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 2 of 3 100.0% processed on watershed_results_wyield_1995.shp 2023-03-30 20:16:14,629 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 2 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:14,755 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 2 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:14,784 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 3 of 3 100.0% processed on watershed_results_wyield_1995.shp 2023-03-30 20:16:14,788 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 3 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:14,924 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 3 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:14,965 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1546) INFO all done processing polygon sets for watershed_results_wyield_1995.shp 2023-03-30 20:16:15,531 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-30 20:16:15,531 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(829) INFO aligned all 1 rasters. 2023-03-30 20:16:15,774 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1310) INFO build a lookup of aggregate field value to FID 2023-03-30 20:16:15,778 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1317) INFO creating disjoint polygon set 2023-03-30 20:16:15,779 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2314) INFO build shapely polygon list 2023-03-30 20:16:16,030 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2339) INFO build shapely rtree index 2023-03-30 20:16:16,032 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2352) INFO poly feature lookup 100.0% complete on watershed_results_wyield_1995.shp 2023-03-30 20:16:16,032 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2356) INFO build poly intersection lookup 2023-03-30 20:16:16,052 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2384) INFO poly intersection feature lookup 100.0% complete on watershed_results_wyield_1995.shp 2023-03-30 20:16:16,053 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2418) INFO maximal subset build 100.0% complete on watershed_results_wyield_1995.shp 2023-03-30 20:16:16,067 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO processing 3 disjoint polygon sets 2023-03-30 20:16:16,076 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 1 of 3 100.0% processed on watershed_results_wyield_1995.shp 2023-03-30 20:16:16,094 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 1 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:16,347 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 1 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:16,447 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 2 of 3 100.0% processed on watershed_results_wyield_1995.shp 2023-03-30 20:16:16,450 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 2 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:16,575 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 2 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:16,603 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 3 of 3 100.0% processed on watershed_results_wyield_1995.shp 2023-03-30 20:16:16,606 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 3 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:16,753 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 3 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:16,812 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1546) INFO all done processing polygon sets for watershed_results_wyield_1995.shp 2023-03-30 20:16:17,444 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(825) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-30 20:16:17,444 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(829) INFO aligned all 1 rasters. 2023-03-30 20:16:17,715 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1310) INFO build a lookup of aggregate field value to FID 2023-03-30 20:16:17,720 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1317) INFO creating disjoint polygon set 2023-03-30 20:16:17,721 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2314) INFO build shapely polygon list 2023-03-30 20:16:17,972 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2339) INFO build shapely rtree index 2023-03-30 20:16:17,974 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2352) INFO poly feature lookup 100.0% complete on watershed_results_wyield_1995.shp 2023-03-30 20:16:17,974 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2356) INFO build poly intersection lookup 2023-03-30 20:16:18,001 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2384) INFO poly intersection feature lookup 100.0% complete on watershed_results_wyield_1995.shp 2023-03-30 20:16:18,001 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2418) INFO maximal subset build 100.0% complete on watershed_results_wyield_1995.shp 2023-03-30 20:16:18,013 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO processing 3 disjoint polygon sets 2023-03-30 20:16:18,016 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 1 of 3 100.0% processed on watershed_results_wyield_1995.shp 2023-03-30 20:16:18,031 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 1 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:18,317 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 1 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:18,423 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 2 of 3 100.0% processed on watershed_results_wyield_1995.shp 2023-03-30 20:16:18,426 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 2 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:18,554 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 2 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:18,598 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1371) INFO disjoint polygon set 3 of 3 100.0% processed on watershed_results_wyield_1995.shp 2023-03-30 20:16:18,603 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1379) INFO rasterizing disjoint polygon set 3 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:18,757 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1397) INFO summarizing rasterized disjoint polygon set 3 of 3 watershed_results_wyield_1995.shp 2023-03-30 20:16:18,816 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1546) INFO all done processing polygon sets for watershed_results_wyield_1995.shp 2023-03-30 20:16:20,782 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 16.71s 2023-03-30 20:16:20,782 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished