2023-03-17 21:00:50,955 (natcap.invest.utils) utils.prepare_workspace(162) Level 100 Writing log messages to F:/乱七八糟\InVEST-Nutrient-Delivery-Ratio-log-2023-03-17--21_00_50.txt 2023-03-17 21:00:50,957 (natcap.invest.ui.model) model._logged_target(1643) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.ndr.ndr 3.10.2: biophysical_table_path D:/GEP项目/生物物理参数表/水质净化生物物理参数表.csv calc_n True calc_p True dem_path D:/GEP项目/陕西省30m精度DEM/填洼_包括陕西省缓冲区dem1.tif k_param 2 lulc_path D:/GEP项目/DLBM 土地利用数据/陕西省土地利用.tif n_workers -1 results_suffix 6 runoff_proxy_path D:/GEP项目/陕西省气象数据2020/陕西省2020年降雨量重采样30m.tif subsurface_critical_length_n 50 subsurface_eff_n 0 threshold_flow_accumulation 1000000 watersheds_path D:/GEP项目/陕西省流域/陕西省流域合并后.shp workspace_dir F:/乱七八糟 2023-03-17 21:00:50,997 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1741) INFO starting reprojection 2023-03-17 21:00:51,001 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer 陕西省流域合并后 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. 2023-03-17 21:00:51,007 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1783) INFO reprojection 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:00:54,877 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:00:54,877 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:00:54,877 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:00:56,020 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:00:56,021 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:00:57,361 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 3 aligned: aligned_dem_6.tif 2023-03-17 21:00:59,193 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:00:59,194 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:00:59,194 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:00:59,651 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:00:59,651 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:00:59,982 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 2 of 3 aligned: aligned_lulc_6.tif 2023-03-17 21:01:03,523 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:01:03,524 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:01:03,524 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:01:04,447 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:01:04,448 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:01:05,582 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 3 of 3 aligned: aligned_runoff_proxy_6.tif 2023-03-17 21:01:05,583 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 3 rasters. 2023-03-17 21:01:14,633 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete 2023-03-17 21:01:25,035 (pygeoprocessing.routing.routing) Task._call(1093) INFO 66.3% complete 2023-03-17 21:01:30,585 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-03-17 21:01:30,593 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating target flow accum raster layer 2023-03-17 21:01:30,603 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating visited raster layer 2023-03-17 21:01:30,614 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG starting search 2023-03-17 21:01:41,002 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 27.3% complete 2023-03-17 21:01:48,042 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-03-17 21:01:53,253 (pygeoprocessing.routing.routing) Task._call(1093) INFO filter out incomplete divergent streams 2023-03-17 21:01:53,571 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-03-17 21:01:56,354 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:01:56,355 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:01:56,355 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:01:57,568 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:01:57,569 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:01:59,718 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:01:59,719 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:01:59,719 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:02:00,767 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:02:00,767 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:02:01,888 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating target flow accum raster layer 2023-03-17 21:02:01,898 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating visited raster layer 2023-03-17 21:02:01,912 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG starting search 2023-03-17 21:02:12,005 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 23.6% complete 2023-03-17 21:02:21,004 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-03-17 21:02:21,032 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:02:21,032 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:02:21,032 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:02:24,854 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:02:24,854 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:02:25,803 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:02:25,804 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:02:25,804 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:02:28,846 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:02:28,847 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:02:30,140 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:02:30,141 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:02:30,142 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:02:31,362 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:02:31,362 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:02:43,006 (pygeoprocessing.routing.routing) Task._call(1093) INFO 66.3% complete 2023-03-17 21:02:50,152 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-03-17 21:03:01,017 (pygeoprocessing.routing.routing) Task._call(1093) INFO 71.4% complete 2023-03-17 21:03:06,421 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-03-17 21:03:06,482 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:03:06,483 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:03:06,483 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:03:09,020 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:03:09,021 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:03:09,264 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:03:09,265 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:03:09,265 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:03:12,128 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:03:12,128 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:03:13,109 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:03:13,110 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:03:13,110 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:03:15,058 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:03:15,058 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:03:15,854 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:03:15,856 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:03:15,856 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:03:17,698 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:03:17,699 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:03:18,962 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:03:18,963 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:03:18,963 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:03:20,637 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:03:20,638 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:03:21,918 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:03:21,919 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:03:21,919 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:03:23,254 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:03:23,254 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:03:24,066 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:03:24,067 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:03:24,067 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:03:25,298 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:03:25,298 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:03:26,116 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:03:26,117 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:03:26,117 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:03:29,091 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:03:29,092 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:03:51,410 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:03:51,410 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:03:51,411 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:03:53,153 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:03:53,153 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:03:54,257 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:03:54,257 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:03:54,257 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:03:55,993 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:03:55,994 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:03:57,164 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:03:57,165 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:03:57,165 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:03:58,935 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:03:58,935 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:03:59,578 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:03:59,579 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:03:59,579 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:04:01,286 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:04:01,287 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:04:01,934 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:04:01,935 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:04:01,935 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:04:03,196 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:04:03,196 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:04:03,880 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:04:03,881 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:04:03,881 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:04:05,491 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:04:05,491 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:04:08,597 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-17 21:04:08,597 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2023-03-17 21:04:08,599 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2023-03-17 21:04:08,600 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2023-03-17 21:04:08,604 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2023-03-17 21:04:08,609 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2023-03-17 21:04:08,612 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:08,612 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2023-03-17 21:04:08,630 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:08,631 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:08,638 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2023-03-17 21:04:08,641 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:08,665 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:09,409 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:10,253 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:10,269 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:11,041 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:11,545 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:11,561 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:12,376 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:12,991 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:13,006 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:13,771 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:14,290 (pygeoprocessing.geoprocessing) geoprocessing.(1306) INFO zonal stats approximately 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:14,291 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:14,306 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:15,072 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:15,597 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 2 of 30 2023-03-17 21:04:15,597 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-1.5874833e+04 8.6883499e+01 0.0000000e+00 4.3814780e+06 0.0000000e+00 -8.6883499e+01] for ('F:/乱七八糟\\n_subsurface_export_6.tif', 1) 2023-03-17 21:04:15,598 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 30 2023-03-17 21:04:15,598 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_6.gpkg 2023-03-17 21:04:18,276 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-17 21:04:18,276 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2023-03-17 21:04:18,278 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2023-03-17 21:04:18,280 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2023-03-17 21:04:18,284 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2023-03-17 21:04:18,287 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2023-03-17 21:04:18,290 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:18,290 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2023-03-17 21:04:18,310 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:18,310 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:18,320 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2023-03-17 21:04:18,322 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:18,347 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:19,127 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:20,213 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:20,228 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:21,002 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:21,522 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:21,537 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:22,367 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:22,973 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:22,990 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:23,752 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:24,276 (pygeoprocessing.geoprocessing) geoprocessing.(1306) INFO zonal stats approximately 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:24,277 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:24,292 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:25,050 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:25,542 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 2 of 30 2023-03-17 21:04:25,543 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-1.5874833e+04 8.6883499e+01 0.0000000e+00 4.3814780e+06 0.0000000e+00 -8.6883499e+01] for ('F:/乱七八糟\\n_total_export_6.tif', 1) 2023-03-17 21:04:25,543 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 30 2023-03-17 21:04:25,543 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_6.gpkg 2023-03-17 21:04:27,578 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-17 21:04:27,579 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2023-03-17 21:04:27,580 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2023-03-17 21:04:27,581 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2023-03-17 21:04:27,584 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2023-03-17 21:04:27,587 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2023-03-17 21:04:27,590 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:27,590 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2023-03-17 21:04:27,608 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:27,608 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:27,616 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2023-03-17 21:04:27,618 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:27,641 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:28,381 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:29,212 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:29,228 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:29,982 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:30,499 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:30,515 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:31,328 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:31,942 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:31,958 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:32,733 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:33,260 (pygeoprocessing.geoprocessing) geoprocessing.(1306) INFO zonal stats approximately 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:33,260 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:33,277 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:34,049 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:34,545 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 2 of 30 2023-03-17 21:04:34,545 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-1.5874833e+04 8.6883499e+01 0.0000000e+00 4.3814780e+06 0.0000000e+00 -8.6883499e+01] for ('F:/乱七八糟\\intermediate_outputs\\sub_load_n_6.tif', 1) 2023-03-17 21:04:34,546 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 30 2023-03-17 21:04:34,546 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_6.gpkg 2023-03-17 21:04:37,178 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-17 21:04:37,178 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2023-03-17 21:04:37,180 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2023-03-17 21:04:37,181 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2023-03-17 21:04:37,184 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2023-03-17 21:04:37,188 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2023-03-17 21:04:37,190 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:37,190 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2023-03-17 21:04:37,209 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:37,209 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:37,218 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2023-03-17 21:04:37,221 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:37,243 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:37,988 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:38,989 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:39,005 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:39,759 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:40,265 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:40,282 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:41,096 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:41,701 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:41,722 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:42,489 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:43,005 (pygeoprocessing.geoprocessing) geoprocessing.(1306) INFO zonal stats approximately 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:43,006 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:43,022 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:43,797 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:44,284 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 2 of 30 2023-03-17 21:04:44,285 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-1.5874833e+04 8.6883499e+01 0.0000000e+00 4.3814780e+06 0.0000000e+00 -8.6883499e+01] for ('F:/乱七八糟\\n_surface_export_6.tif', 1) 2023-03-17 21:04:44,285 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 30 2023-03-17 21:04:44,285 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_6.gpkg 2023-03-17 21:04:47,127 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-17 21:04:47,127 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2023-03-17 21:04:47,129 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2023-03-17 21:04:47,130 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2023-03-17 21:04:47,133 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2023-03-17 21:04:47,136 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2023-03-17 21:04:47,139 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:47,139 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2023-03-17 21:04:47,157 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:47,157 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:47,167 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2023-03-17 21:04:47,169 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:47,194 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:47,935 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:48,978 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:48,994 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:49,750 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:50,288 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:50,304 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:51,118 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:51,757 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:51,775 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:52,590 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:53,155 (pygeoprocessing.geoprocessing) geoprocessing.(1306) INFO zonal stats approximately 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:04:53,156 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:04:53,172 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:53,933 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:04:54,434 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 2 of 30 2023-03-17 21:04:54,435 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-1.5874833e+04 8.6883499e+01 0.0000000e+00 4.3814780e+06 0.0000000e+00 -8.6883499e+01] for ('F:/乱七八糟\\intermediate_outputs\\surface_load_n_6.tif', 1) 2023-03-17 21:04:54,435 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 30 2023-03-17 21:04:54,435 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_6.gpkg 2023-03-17 21:04:54,545 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:04:54,546 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:04:54,546 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:04:56,533 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:04:56,534 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:04:57,353 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:04:57,354 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:04:57,354 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:04:59,351 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:04:59,351 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:05:00,716 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:05:00,717 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:05:00,717 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:05:01,898 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:05:01,898 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:05:03,141 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:05:03,142 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:05:03,142 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:05:04,410 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:05:04,410 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:05:05,265 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:05:05,266 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:05:05,266 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:05:06,604 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:05:06,605 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:05:07,434 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:05:07,435 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:05:07,435 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:05:10,432 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:05:10,432 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:05:32,887 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:05:32,887 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:05:32,887 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:05:34,619 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:05:34,619 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:05:35,728 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-03-17 21:05:35,728 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 3836 2023-03-17 21:05:35,728 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-03-17 21:05:37,552 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-03-17 21:05:37,552 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-03-17 21:05:41,120 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-17 21:05:41,120 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2023-03-17 21:05:41,121 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2023-03-17 21:05:41,123 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2023-03-17 21:05:41,127 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2023-03-17 21:05:41,130 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2023-03-17 21:05:41,133 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:05:41,133 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2023-03-17 21:05:41,152 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:05:41,153 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:05:41,162 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2023-03-17 21:05:41,165 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:05:41,187 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:41,962 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:42,954 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:05:42,970 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:43,727 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:44,239 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:05:44,255 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:45,075 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:45,699 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:05:45,715 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:46,483 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:47,004 (pygeoprocessing.geoprocessing) geoprocessing.(1306) INFO zonal stats approximately 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:05:47,004 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:05:47,020 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:47,793 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:48,289 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 2 of 30 2023-03-17 21:05:48,289 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-1.5874833e+04 8.6883499e+01 0.0000000e+00 4.3814780e+06 0.0000000e+00 -8.6883499e+01] for ('F:/乱七八糟\\p_surface_export_6.tif', 1) 2023-03-17 21:05:48,290 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 30 2023-03-17 21:05:48,290 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_6.gpkg 2023-03-17 21:05:51,118 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2023-03-17 21:05:51,118 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2023-03-17 21:05:51,120 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2023-03-17 21:05:51,121 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2023-03-17 21:05:51,124 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2023-03-17 21:05:51,127 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2023-03-17 21:05:51,130 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:05:51,130 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2023-03-17 21:05:51,147 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:05:51,148 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:05:51,158 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 5 disjoint polygon sets 2023-03-17 21:05:51,160 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:05:51,184 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:51,928 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:52,972 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:05:52,991 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:53,751 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:54,294 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:05:54,310 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:55,121 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:55,741 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:05:55,758 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:56,521 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:57,049 (pygeoprocessing.geoprocessing) geoprocessing.(1306) INFO zonal stats approximately 100.0% complete on watershed_results_ndr_6.gpkg 2023-03-17 21:05:57,050 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_ndr_6.gpkg 2023-03-17 21:05:57,066 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:57,823 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_ndr_6.gpkg 2023-03-17 21:05:58,326 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 2 of 30 2023-03-17 21:05:58,326 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-1.5874833e+04 8.6883499e+01 0.0000000e+00 4.3814780e+06 0.0000000e+00 -8.6883499e+01] for ('F:/乱七八糟\\intermediate_outputs\\surface_load_p_6.tif', 1) 2023-03-17 21:05:58,327 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 30 2023-03-17 21:05:58,327 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_6.gpkg 2023-03-17 21:05:58,414 (natcap.invest.ndr.ndr) ndr.execute(730) INFO Writing summaries to output shapefile 2023-03-17 21:05:58,437 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer 陕西省流域合并后 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. 2023-03-17 21:05:58,488 (natcap.invest.ndr.ndr) ndr.execute(734) INFO NDR complete! 2023-03-17 21:05:58,488 (natcap.invest.ndr.ndr) ndr.execute(735) INFO _ _ ____ ____ 2023-03-17 21:05:58,488 (natcap.invest.ndr.ndr) ndr.execute(736) INFO | \ |"| | _"\U | _"\ u 2023-03-17 21:05:58,488 (natcap.invest.ndr.ndr) ndr.execute(737) INFO <| \| |>/| | | |\| |_) |/ 2023-03-17 21:05:58,488 (natcap.invest.ndr.ndr) ndr.execute(738) INFO U| |\ |uU| |_| |\| _ < 2023-03-17 21:05:58,488 (natcap.invest.ndr.ndr) ndr.execute(739) INFO |_| \_| |____/ u|_| \_\ 2023-03-17 21:05:58,489 (natcap.invest.ndr.ndr) ndr.execute(740) INFO || \\,-.|||_ // \\_ 2023-03-17 21:05:58,489 (natcap.invest.ndr.ndr) ndr.execute(741) INFO (_") (_/(__)_) (__) (__) 2023-03-17 21:05:58,490 (natcap.invest.utils) utils.prepare_workspace(170) INFO Elapsed time: 5m 7.529999999999973s 2023-03-17 21:05:58,490 (natcap.invest.utils) utils.prepare_workspace(173) INFO Execution finished