2022-07-15 12:53:39,944 (natcap.invest.utils) utils.prepare_workspace(162) Level 100 Writing log messages to [C:/Users/yanyan/Desktop/InVEST NDR/new version\InVEST-Nutrient-Delivery-Ratio-log-2022-07-15--12_53_39.txt] 2022-07-15 12:53:39,946 (natcap.invest.ui.model) model._logged_target(1626) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.ndr.ndr 3.11.0: biophysical_table_path C:/Users/yanyan/Desktop/InVEST NDR/NRD biophysical table.csv calc_n True calc_p True dem_path C:/Users/yanyan/Desktop/DEM/DEM 50N/WuhanDEM_50N_30m.tif k_param 2 lulc_path C:/Users/yanyan/Desktop/final calculation/2015Y LULC projected raster new/2015_WGS1984 50N_30m.tif n_workers -1 results_suffix runoff_proxy_path C:/Users/yanyan/Desktop/InVEST NDR/Nutrient Runoff Proxy/avarage annual pre.tif subsurface_critical_length_n 200 subsurface_eff_n 0.8 threshold_flow_accumulation 20000 watersheds_path C:/Users/yanyan/Desktop/DEM/whole DEM to watershed/Wuhan_watershed_50N.shp workspace_dir C:/Users/yanyan/Desktop/InVEST NDR/new version 2022-07-15 12:53:40,035 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1741) INFO starting reprojection 2022-07-15 12:53:40,040 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer Wuhan_watershed_50N of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2022-07-15 12:53:40,045 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1783) INFO reprojection 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:53:41,753 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:53:41,754 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:53:41,754 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:53:42,166 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:53:42,167 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:53:42,597 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 3 aligned: aligned_dem.tif 2022-07-15 12:53:43,425 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:53:43,426 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:53:43,426 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:53:43,618 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:53:43,618 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:53:43,721 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 2 of 3 aligned: aligned_lulc.tif 2022-07-15 12:53:45,116 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:53:45,117 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:53:45,117 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:53:45,486 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:53:45,486 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:53:45,896 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 3 of 3 aligned: aligned_runoff_proxy.tif 2022-07-15 12:53:45,897 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 3 rasters. 2022-07-15 12:53:56,002 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): 12626432 of 23118831 pixels complete 2022-07-15 12:54:07,004 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): 20658688 of 23118831 pixels complete 2022-07-15 12:54:08,459 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete 2022-07-15 12:54:19,272 (pygeoprocessing.routing.routing) Task._call(1093) INFO 49.6% complete 2022-07-15 12:54:30,869 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-07-15 12:54:30,884 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating target flow accum raster layer 2022-07-15 12:54:30,897 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating visited raster layer 2022-07-15 12:54:30,924 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG starting search 2022-07-15 12:54:41,025 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 31.4% complete 2022-07-15 12:54:44,916 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-07-15 12:54:50,540 (pygeoprocessing.routing.routing) Task._call(1093) INFO filter out incomplete divergent streams 2022-07-15 12:54:50,754 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-07-15 12:54:52,722 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:54:52,723 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:54:52,723 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:54:53,695 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:54:53,696 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:54:55,191 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:54:55,193 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:54:55,193 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:54:56,019 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:54:56,020 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:54:56,794 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating target flow accum raster layer 2022-07-15 12:54:56,823 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating visited raster layer 2022-07-15 12:54:56,846 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG starting search 2022-07-15 12:55:07,170 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-07-15 12:55:07,204 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:07,205 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:07,205 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:08,937 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:08,937 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:09,525 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:09,527 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:09,527 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:10,727 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:10,727 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:11,341 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:11,342 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:11,342 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:11,754 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:11,754 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:19,197 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-07-15 12:55:25,765 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-07-15 12:55:25,802 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:25,803 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:25,803 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:26,796 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:26,796 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:26,895 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:26,897 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:26,897 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:28,043 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:28,044 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:28,556 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:28,556 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:28,556 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:28,970 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:28,970 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:29,277 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:29,278 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:29,278 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:29,940 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:29,941 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:30,379 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:30,380 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:30,380 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:30,881 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:30,881 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:31,308 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:31,309 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:31,309 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:31,667 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:31,667 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:32,029 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:32,063 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:32,063 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:32,676 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:32,676 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:33,008 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:33,009 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:33,010 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:34,030 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:34,030 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:43,170 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:43,171 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:43,172 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:43,873 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:43,873 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:44,389 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:44,390 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:44,390 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:44,994 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:44,994 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:45,552 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:45,553 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:45,553 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:46,099 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:46,100 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:46,368 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:46,369 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:46,370 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:47,074 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:47,075 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:47,535 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:47,536 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:47,536 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:48,093 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:48,094 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:48,367 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:55:48,368 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:55:48,368 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:55:48,974 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:55:48,975 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:55:50,507 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-07-15 12:55:50,507 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-07-15 12:55:50,512 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-07-15 12:55:50,513 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-07-15 12:55:50,518 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-07-15 12:55:50,521 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-07-15 12:55:50,523 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:55:50,523 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-07-15 12:55:50,532 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:55:50,533 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:55:50,546 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2022-07-15 12:55:50,548 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:50,553 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:50,837 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:51,144 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:51,148 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:51,466 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:51,690 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:51,694 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:51,998 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:52,203 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:52,208 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:52,515 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:52,712 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:52,716 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:53,007 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:53,206 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 4 of 16 2022-07-15 12:55:53,206 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.8218986e+05 3.0000000e+01 0.0000000e+00 3.4729610e+06 0.0000000e+00 -3.0000000e+01] for ('C:/Users/yanyan/Desktop/InVEST NDR/new version\\n_subsurface_export.tif', 1) 2022-07-15 12:55:53,207 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 16 2022-07-15 12:55:53,207 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr.gpkg 2022-07-15 12:55:54,371 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-07-15 12:55:54,372 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-07-15 12:55:54,381 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-07-15 12:55:54,383 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-07-15 12:55:54,407 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-07-15 12:55:54,411 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-07-15 12:55:54,414 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:55:54,415 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-07-15 12:55:54,430 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:55:54,431 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:55:54,450 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2022-07-15 12:55:54,452 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:54,460 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:54,873 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:55,306 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:55,310 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:55,620 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:55,827 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:55,832 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:56,141 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:56,345 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:56,349 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:56,650 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:56,851 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:56,855 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:57,152 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:57,352 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 4 of 16 2022-07-15 12:55:57,353 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.8218986e+05 3.0000000e+01 0.0000000e+00 3.4729610e+06 0.0000000e+00 -3.0000000e+01] for ('C:/Users/yanyan/Desktop/InVEST NDR/new version\\n_total_export.tif', 1) 2022-07-15 12:55:57,354 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 16 2022-07-15 12:55:57,354 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr.gpkg 2022-07-15 12:55:58,238 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-07-15 12:55:58,239 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-07-15 12:55:58,246 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-07-15 12:55:58,247 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-07-15 12:55:58,254 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-07-15 12:55:58,256 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-07-15 12:55:58,259 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:55:58,259 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-07-15 12:55:58,268 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:55:58,269 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:55:58,282 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2022-07-15 12:55:58,284 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:58,288 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:58,596 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:58,917 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:58,922 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:59,242 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:59,461 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:59,465 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:59,771 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:55:59,975 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:55:59,979 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:00,289 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:00,492 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:00,496 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:00,801 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:01,029 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 4 of 16 2022-07-15 12:56:01,029 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.8218986e+05 3.0000000e+01 0.0000000e+00 3.4729610e+06 0.0000000e+00 -3.0000000e+01] for ('C:/Users/yanyan/Desktop/InVEST NDR/new version\\intermediate_outputs\\sub_load_n.tif', 1) 2022-07-15 12:56:01,031 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 16 2022-07-15 12:56:01,031 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr.gpkg 2022-07-15 12:56:02,403 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-07-15 12:56:02,404 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-07-15 12:56:02,409 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-07-15 12:56:02,410 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-07-15 12:56:02,415 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-07-15 12:56:02,417 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-07-15 12:56:02,419 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:56:02,420 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-07-15 12:56:02,429 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:56:02,429 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:56:02,441 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2022-07-15 12:56:02,443 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:02,448 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:02,733 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:03,134 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:03,138 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:03,447 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:03,653 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:03,657 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:03,964 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:04,168 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:04,172 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:04,481 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:04,680 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:04,684 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:04,977 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:05,181 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 4 of 16 2022-07-15 12:56:05,181 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.8218986e+05 3.0000000e+01 0.0000000e+00 3.4729610e+06 0.0000000e+00 -3.0000000e+01] for ('C:/Users/yanyan/Desktop/InVEST NDR/new version\\n_surface_export.tif', 1) 2022-07-15 12:56:05,182 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 16 2022-07-15 12:56:05,182 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr.gpkg 2022-07-15 12:56:06,147 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-07-15 12:56:06,147 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-07-15 12:56:06,155 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-07-15 12:56:06,156 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-07-15 12:56:06,160 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-07-15 12:56:06,162 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-07-15 12:56:06,164 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:56:06,165 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-07-15 12:56:06,173 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:56:06,173 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:56:06,185 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2022-07-15 12:56:06,188 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:06,192 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:06,565 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:06,995 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:07,000 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:07,328 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:07,567 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:07,572 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:07,900 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:08,104 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:08,108 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:08,413 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:08,614 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:08,618 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:08,910 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:09,111 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 4 of 16 2022-07-15 12:56:09,111 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.8218986e+05 3.0000000e+01 0.0000000e+00 3.4729610e+06 0.0000000e+00 -3.0000000e+01] for ('C:/Users/yanyan/Desktop/InVEST NDR/new version\\intermediate_outputs\\surface_load_n.tif', 1) 2022-07-15 12:56:09,113 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 16 2022-07-15 12:56:09,113 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr.gpkg 2022-07-15 12:56:09,144 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:56:09,145 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:56:09,145 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:56:09,700 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:56:09,700 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:56:10,040 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:56:10,041 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:56:10,041 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:56:10,852 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:56:10,853 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:56:11,294 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:56:11,295 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:56:11,295 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:56:11,630 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:56:11,630 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:56:12,127 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:56:12,128 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:56:12,128 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:56:12,917 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:56:12,918 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:56:13,391 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:56:13,392 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:56:13,393 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:56:13,903 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:56:13,903 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:56:14,368 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:56:14,370 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:56:14,370 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:56:16,185 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:56:16,187 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:56:33,309 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:56:33,310 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:56:33,311 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:56:34,727 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:56:34,728 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:56:35,708 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-15 12:56:35,710 (pygeoprocessing.geoprocessing_core) threading.run(946) DEBUG stats worker PID: 29104 2022-07-15 12:56:35,710 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-15 12:56:36,916 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-15 12:56:36,917 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-15 12:56:39,965 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-07-15 12:56:39,966 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-07-15 12:56:39,977 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-07-15 12:56:39,978 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-07-15 12:56:39,985 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-07-15 12:56:39,988 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-07-15 12:56:39,991 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:56:39,991 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-07-15 12:56:40,002 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:56:40,003 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:56:40,023 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2022-07-15 12:56:40,028 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:40,042 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:40,546 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:41,439 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:41,449 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:42,042 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:42,441 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:42,446 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:43,094 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:43,498 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:43,509 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:44,163 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:44,582 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:44,594 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:45,114 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:45,541 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 4 of 16 2022-07-15 12:56:45,542 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.8218986e+05 3.0000000e+01 0.0000000e+00 3.4729610e+06 0.0000000e+00 -3.0000000e+01] for ('C:/Users/yanyan/Desktop/InVEST NDR/new version\\p_surface_export.tif', 1) 2022-07-15 12:56:45,543 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 16 2022-07-15 12:56:45,543 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr.gpkg 2022-07-15 12:56:47,123 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-07-15 12:56:47,124 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-07-15 12:56:47,131 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-07-15 12:56:47,132 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-07-15 12:56:47,138 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-07-15 12:56:47,141 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-07-15 12:56:47,143 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:56:47,144 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-07-15 12:56:47,153 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:56:47,153 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr.gpkg 2022-07-15 12:56:47,166 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2022-07-15 12:56:47,170 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:47,175 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:47,483 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:47,859 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:47,864 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:48,177 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:48,386 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:48,390 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:48,696 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:48,901 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:48,905 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:49,216 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:49,423 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr.gpkg 2022-07-15 12:56:49,426 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:49,725 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr.gpkg 2022-07-15 12:56:49,925 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 4 of 16 2022-07-15 12:56:49,925 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.8218986e+05 3.0000000e+01 0.0000000e+00 3.4729610e+06 0.0000000e+00 -3.0000000e+01] for ('C:/Users/yanyan/Desktop/InVEST NDR/new version\\intermediate_outputs\\surface_load_p.tif', 1) 2022-07-15 12:56:49,926 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 16 2022-07-15 12:56:49,926 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr.gpkg 2022-07-15 12:56:49,941 (natcap.invest.ndr.ndr) ndr.execute(732) INFO Writing summaries to output shapefile 2022-07-15 12:56:49,991 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer Wuhan_watershed_50N of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2022-07-15 12:56:50,033 (natcap.invest.ndr.ndr) ndr.execute(736) INFO NDR complete! 2022-07-15 12:56:50,034 (natcap.invest.ndr.ndr) ndr.execute(737) INFO _ _ ____ ____ 2022-07-15 12:56:50,034 (natcap.invest.ndr.ndr) ndr.execute(738) INFO | \ |"| | _"\U | _"\ u 2022-07-15 12:56:50,035 (natcap.invest.ndr.ndr) ndr.execute(739) INFO <| \| |>/| | | |\| |_) |/ 2022-07-15 12:56:50,035 (natcap.invest.ndr.ndr) ndr.execute(740) INFO U| |\ |uU| |_| |\| _ < 2022-07-15 12:56:50,035 (natcap.invest.ndr.ndr) ndr.execute(741) INFO |_| \_| |____/ u|_| \_\ 2022-07-15 12:56:50,035 (natcap.invest.ndr.ndr) ndr.execute(742) INFO || \\,-.|||_ // \\_ 2022-07-15 12:56:50,035 (natcap.invest.ndr.ndr) ndr.execute(743) INFO (_") (_/(__)_) (__) (__) 2022-07-15 12:56:50,036 (natcap.invest.utils) utils.prepare_workspace(170) INFO Elapsed time: 3m 10.090000000000003s 2022-07-15 12:56:50,037 (natcap.invest.utils) utils.prepare_workspace(173) INFO Execution finished