2022-07-29 21:46:25,761 (natcap.invest.utils) utils.prepare_workspace(162) Level 100 Writing log messages to [F:\sanxia\water_purification\2000\InVEST-natcap.invest.ndr.ndr-log-2022-07-29--21_46_25.txt] 2022-07-29 21:46:25,761 (__main__) cli.main(475) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.ndr.ndr 3.11.0: biophysical_table_path F:\sanxia\water_purification\BiophysicalTable.csv calc_n True calc_p False dem_path F:\sanxia\water_purification\Dem_fill.tif k_param 2 lulc_path F:\sanxia\lucc\c2000.tif n_workers -1 results_suffix 2000 runoff_proxy_path F:\sanxia\water_purification\copy_rpi_2000.tif subsurface_critical_length_n 200 subsurface_eff_n 0.8 threshold_flow_accumulation 1000 watersheds_path F:\sanxia\water_purification\sanxia.shp workspace_dir F:\sanxia\water_purification\2000 2022-07-29 21:46:25,956 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1741) INFO starting reprojection 2022-07-29 21:46:25,960 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1783) INFO reprojection 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:46:31,102 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 60.0% complete F:\sanxia\water_purification\2000\intermediate_outputs\cache_dir\aligned_dem_2000.tif 2022-07-29 21:46:34,032 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 100.0% complete F:\sanxia\water_purification\2000\intermediate_outputs\cache_dir\aligned_dem_2000.tif 2022-07-29 21:46:40,938 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:46:40,939 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:46:44,160 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:46:44,161 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:46:46,389 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 3 aligned: aligned_dem_2000.tif 2022-07-29 21:46:51,471 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 88.0% complete F:\sanxia\water_purification\2000\intermediate_outputs\cache_dir\aligned_lulc_2000.tif 2022-07-29 21:46:51,662 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 100.0% complete F:\sanxia\water_purification\2000\intermediate_outputs\cache_dir\aligned_lulc_2000.tif 2022-07-29 21:46:54,768 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:46:54,769 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:46:56,309 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:46:56,309 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:46:57,179 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 2 of 3 aligned: aligned_lulc_2000.tif 2022-07-29 21:47:02,217 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 83.0% complete F:\sanxia\water_purification\2000\intermediate_outputs\cache_dir\aligned_runoff_proxy_2000.tif 2022-07-29 21:47:03,159 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 100.0% complete F:\sanxia\water_purification\2000\intermediate_outputs\cache_dir\aligned_runoff_proxy_2000.tif 2022-07-29 21:47:08,442 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:47:08,443 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:47:11,322 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:47:11,323 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:47:14,819 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 3 of 3 aligned: aligned_runoff_proxy_2000.tif 2022-07-29 21:47:14,820 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 3 rasters. 2022-07-29 21:47:25,000 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): 26554368 of 119713280 pixels complete 2022-07-29 21:47:36,003 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): 106206464 of 119713280 pixels complete 2022-07-29 21:47:37,101 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete 2022-07-29 21:47:48,021 (pygeoprocessing.routing.routing) Task._call(1093) INFO 6.3% complete 2022-07-29 21:47:59,035 (pygeoprocessing.routing.routing) Task._call(1093) INFO 22.2% complete 2022-07-29 21:48:10,000 (pygeoprocessing.routing.routing) Task._call(1093) INFO 57.0% complete 2022-07-29 21:48:15,636 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-07-29 21:48:26,010 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 8.5% complete 2022-07-29 21:48:37,000 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 16.8% complete 2022-07-29 21:48:48,014 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 22.1% complete 2022-07-29 21:48:55,205 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-07-29 21:49:06,000 (pygeoprocessing.routing.routing) Task._call(1093) INFO 19.0% complete 2022-07-29 21:49:17,000 (pygeoprocessing.routing.routing) Task._call(1093) INFO 69.7% complete 2022-07-29 21:49:20,773 (pygeoprocessing.routing.routing) Task._call(1093) INFO filter out incomplete divergent streams 2022-07-29 21:49:21,928 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-07-29 21:49:29,711 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:49:29,712 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:49:34,711 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 94.9% complete 2022-07-29 21:49:34,991 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:49:34,991 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:49:42,030 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:49:42,031 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:49:45,866 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:49:45,867 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:50:00,000 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 7.6% complete 2022-07-29 21:50:11,000 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 14.8% complete 2022-07-29 21:50:22,000 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 20.7% complete 2022-07-29 21:50:33,001 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 25.8% complete 2022-07-29 21:50:34,567 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-07-29 21:50:34,698 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:50:34,700 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:50:39,699 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 28.5% complete 2022-07-29 21:50:44,705 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 52.2% complete 2022-07-29 21:50:49,705 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 80.6% complete 2022-07-29 21:50:52,897 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:50:52,897 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:50:55,601 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:50:55,602 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:51:00,603 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 36.8% complete 2022-07-29 21:51:05,608 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 73.4% complete 2022-07-29 21:51:09,374 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:51:09,374 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:51:13,333 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:51:13,334 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:51:18,035 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:51:18,036 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:51:33,138 (pygeoprocessing.routing.routing) Task._call(1093) INFO 12.7% complete 2022-07-29 21:51:44,025 (pygeoprocessing.routing.routing) Task._call(1093) INFO 25.3% complete 2022-07-29 21:51:55,096 (pygeoprocessing.routing.routing) Task._call(1093) INFO 44.4% complete 2022-07-29 21:52:06,040 (pygeoprocessing.routing.routing) Task._call(1093) INFO 88.7% complete 2022-07-29 21:52:07,939 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-07-29 21:52:19,085 (pygeoprocessing.routing.routing) Task._call(1093) INFO 15.8% complete 2022-07-29 21:52:30,010 (pygeoprocessing.routing.routing) Task._call(1093) INFO 34.9% complete 2022-07-29 21:52:41,025 (pygeoprocessing.routing.routing) Task._call(1093) INFO 72.9% complete 2022-07-29 21:52:45,500 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-07-29 21:52:45,705 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:52:45,706 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:52:50,706 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 50.0% complete 2022-07-29 21:52:55,486 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:52:55,487 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:52:56,434 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:52:56,435 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:53:01,442 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 35.6% complete 2022-07-29 21:53:06,444 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 72.8% complete 2022-07-29 21:53:10,720 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:53:10,721 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:53:14,555 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:53:14,556 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:53:18,281 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:53:18,281 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:53:21,125 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:53:21,126 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:53:26,129 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 61.6% complete 2022-07-29 21:53:29,519 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:53:29,519 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:53:33,681 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:53:33,682 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:53:38,681 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 92.2% complete 2022-07-29 21:53:39,142 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:53:39,143 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:53:43,213 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:53:43,214 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:53:45,910 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:53:45,910 (pygeoprocessing.geoprocessing_core) threading.run(946) WARNING No valid pixels were received, sending None. 2022-07-29 21:53:45,910 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:53:48,186 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:53:48,187 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:53:50,618 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:53:50,618 (pygeoprocessing.geoprocessing_core) threading.run(946) WARNING No valid pixels were received, sending None. 2022-07-29 21:53:50,619 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:53:53,001 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:53:53,002 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:53:58,004 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 45.2% complete 2022-07-29 21:54:03,009 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 84.5% complete 2022-07-29 21:54:04,985 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:54:04,985 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:54:49,235 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:54:49,236 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:54:54,237 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 60.5% complete 2022-07-29 21:54:57,608 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:54:57,608 (pygeoprocessing.geoprocessing_core) threading.run(946) WARNING No valid pixels were received, sending None. 2022-07-29 21:54:57,608 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:55:00,212 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:55:00,214 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:55:05,210 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 83.9% complete 2022-07-29 21:55:06,194 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:55:06,194 (pygeoprocessing.geoprocessing_core) threading.run(946) WARNING No valid pixels were received, sending None. 2022-07-29 21:55:06,195 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:55:08,805 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:55:08,806 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:55:13,804 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 85.2% complete 2022-07-29 21:55:14,629 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:55:14,629 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:55:17,774 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:55:17,775 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:55:22,775 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 50.6% complete 2022-07-29 21:55:27,473 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:55:27,473 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:55:31,451 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:55:31,452 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:55:36,452 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 83.4% complete 2022-07-29 21:55:37,448 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:55:37,448 (pygeoprocessing.geoprocessing_core) threading.run(946) WARNING No valid pixels were received, sending None. 2022-07-29 21:55:37,449 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:55:39,849 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-07-29 21:55:39,850 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-07-29 21:55:44,848 (pygeoprocessing.geoprocessing) geoprocessing.(477) INFO 73.9% complete 2022-07-29 21:55:46,534 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-07-29 21:55:46,535 (pygeoprocessing.geoprocessing_core) threading.run(946) WARNING No valid pixels were received, sending None. 2022-07-29 21:55:46,535 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-07-29 21:55:55,567 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-07-29 21:55:55,568 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-07-29 21:55:55,570 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-07-29 21:55:55,571 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-07-29 21:55:55,578 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-07-29 21:55:55,589 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-07-29 21:55:55,602 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:55:55,602 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-07-29 21:55:55,604 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:55:55,604 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:55:55,616 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-07-29 21:55:55,617 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on watershed_results_ndr_2000.gpkg 2022-07-29 21:55:56,019 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 watershed_results_ndr_2000.gpkg 2022-07-29 21:55:58,828 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 watershed_results_ndr_2000.gpkg 2022-07-29 21:56:03,547 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_2000.gpkg 2022-07-29 21:56:10,115 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-07-29 21:56:10,116 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-07-29 21:56:10,120 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-07-29 21:56:10,121 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-07-29 21:56:10,134 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-07-29 21:56:10,136 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-07-29 21:56:10,139 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:10,140 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-07-29 21:56:10,140 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:10,141 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:10,164 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-07-29 21:56:10,165 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:10,321 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 watershed_results_ndr_2000.gpkg 2022-07-29 21:56:13,161 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 watershed_results_ndr_2000.gpkg 2022-07-29 21:56:17,688 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_2000.gpkg 2022-07-29 21:56:25,657 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-07-29 21:56:25,658 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-07-29 21:56:25,660 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-07-29 21:56:25,660 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-07-29 21:56:25,667 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-07-29 21:56:25,668 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-07-29 21:56:25,670 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:25,670 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-07-29 21:56:25,670 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:25,671 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:25,682 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-07-29 21:56:25,683 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:25,766 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 watershed_results_ndr_2000.gpkg 2022-07-29 21:56:28,619 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 watershed_results_ndr_2000.gpkg 2022-07-29 21:56:33,183 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_2000.gpkg 2022-07-29 21:56:40,108 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-07-29 21:56:40,109 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-07-29 21:56:40,112 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-07-29 21:56:40,113 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-07-29 21:56:40,124 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-07-29 21:56:40,126 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-07-29 21:56:40,129 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:40,129 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-07-29 21:56:40,130 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:40,130 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:40,150 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-07-29 21:56:40,151 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:40,253 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 watershed_results_ndr_2000.gpkg 2022-07-29 21:56:43,327 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 watershed_results_ndr_2000.gpkg 2022-07-29 21:56:47,411 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_2000.gpkg 2022-07-29 21:56:52,757 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 88.0% complete F:\sanxia\water_purification\2000\intermediate_outputs\cache_dir\tmpsq4bi4tm\clipped_raster.tif 2022-07-29 21:56:53,150 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 100.0% complete F:\sanxia\water_purification\2000\intermediate_outputs\cache_dir\tmpsq4bi4tm\clipped_raster.tif 2022-07-29 21:56:56,536 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-07-29 21:56:56,536 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-07-29 21:56:56,538 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-07-29 21:56:56,539 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-07-29 21:56:56,546 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-07-29 21:56:56,548 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-07-29 21:56:56,549 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:56,549 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-07-29 21:56:56,550 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:56,550 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:56,575 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-07-29 21:56:56,576 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on watershed_results_ndr_2000.gpkg 2022-07-29 21:56:56,661 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 watershed_results_ndr_2000.gpkg 2022-07-29 21:56:59,453 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 watershed_results_ndr_2000.gpkg 2022-07-29 21:57:04,165 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_2000.gpkg 2022-07-29 21:57:04,356 (natcap.invest.ndr.ndr) ndr.execute(732) INFO Writing summaries to output shapefile 2022-07-29 21:57:04,581 (natcap.invest.ndr.ndr) ndr.execute(736) INFO NDR complete! 2022-07-29 21:57:04,581 (natcap.invest.ndr.ndr) ndr.execute(737) INFO _ _ ____ ____ 2022-07-29 21:57:04,582 (natcap.invest.ndr.ndr) ndr.execute(738) INFO | \ |"| | _"\U | _"\ u 2022-07-29 21:57:04,582 (natcap.invest.ndr.ndr) ndr.execute(739) INFO <| \| |>/| | | |\| |_) |/ 2022-07-29 21:57:04,582 (natcap.invest.ndr.ndr) ndr.execute(740) INFO U| |\ |uU| |_| |\| _ < 2022-07-29 21:57:04,582 (natcap.invest.ndr.ndr) ndr.execute(741) INFO |_| \_| |____/ u|_| \_\ 2022-07-29 21:57:04,582 (natcap.invest.ndr.ndr) ndr.execute(742) INFO || \\,-.|||_ // \\_ 2022-07-29 21:57:04,582 (natcap.invest.ndr.ndr) ndr.execute(743) INFO (_") (_/(__)_) (__) (__) 2022-07-29 21:57:04,582 (natcap.invest.utils) utils.prepare_workspace(170) INFO Elapsed time: 10m 38.82000000000005s 2022-07-29 21:57:04,583 (natcap.invest.utils) utils.prepare_workspace(173) INFO Execution finished