2023-10-23 14:43:57,478 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [N:\PACIFIC\InVEST results\InVEST-natcap.invest.ndr.ndr-log-2023-10-23--14_43_57.txt] 2023-10-23 14:43:57,478 (__main__) cli.main(459) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.ndr.ndr 3.14.0: biophysical_table_path N:\PACIFIC\InVEST workspace\biophysical_table_benta.csv calc_n False calc_p True dem_path N:\PACIFIC\InVEST workspace\dtm_8_buff.tif k_param 2 lulc_path N:\PACIFIC\GIS\lcm2015gb25m_1.tif n_workers -1 results_suffix _benta runoff_proxy_path N:\PACIFIC\InVEST workspace\HadUK-Grid_MeanAnnualRainfall_1990_2017.tif subsurface_critical_length_n subsurface_eff_n threshold_flow_accumulation 100 watersheds_path N:\PACIFIC\InVEST workspace\WFD_Management_Catchments_Cycle_1_square_repair.shp workspace_dir N:\PACIFIC\InVEST results 2023-10-23 14:43:57,628 (taskgraph.Task) Task.is_precalculated(1252) INFO not precalculated (create target vector (0)), Task hash exists, but there are these mismatches: Modified times don't match cached: (1698065105662017280.000000) actual: (1698065513948575488.000000) 2023-10-23 14:43:57,944 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2037) INFO starting reprojection 2023-10-23 14:43:57,979 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer WFD_Management_Catchments_Cycle_1_square_repair 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-10-23 14:43:58,178 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2077) INFO reprojection 100.0% complete on watershed_results_ndr_benta.gpkg 2023-10-23 14:44:04,502 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 3 aligned: aligned_dem_benta.tif 2023-10-23 14:44:08,188 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 2 of 3 aligned: aligned_lulc_benta.tif 2023-10-23 14:44:10,841 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 3 of 3 aligned: aligned_runoff_proxy_benta.tif 2023-10-23 14:44:10,858 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 3 rasters. 2023-10-23 14:44:56,003 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): 23103459 of 42178224 pixels complete 2023-10-23 14:45:05,078 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete 2023-10-23 14:45:16,122 (pygeoprocessing.routing.routing) Task._call(1093) INFO 43.3% complete 2023-10-23 14:45:27,017 (pygeoprocessing.routing.routing) Task._call(1093) INFO 59.0% complete 2023-10-23 14:45:38,362 (pygeoprocessing.routing.routing) Task._call(1093) INFO 74.7% complete 2023-10-23 14:45:49,089 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-10-23 14:46:00,051 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 12.1% complete 2023-10-23 14:46:11,010 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 19.9% complete 2023-10-23 14:46:22,022 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 29.3% complete 2023-10-23 14:46:25,776 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-10-23 14:46:37,238 (pygeoprocessing.routing.routing) Task._call(1093) INFO 70.8% complete 2023-10-23 14:46:42,676 (pygeoprocessing.routing.routing) Task._call(1093) INFO filter out incomplete divergent streams 2023-10-23 14:46:43,027 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-10-23 14:46:48,049 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:46:48,052 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:46:49,416 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:46:49,423 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:46:53,670 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:46:53,677 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:46:54,670 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:46:54,672 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:47:06,132 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 11.2% complete 2023-10-23 14:47:17,060 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 18.8% complete 2023-10-23 14:47:28,038 (pygeoprocessing.routing.routing) Task._call(1093) INFO mfd flow accum 28.5% complete 2023-10-23 14:47:32,558 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-10-23 14:47:33,560 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:47:33,572 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:47:38,548 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO s_bar_benta.tif 74.3% complete 2023-10-23 14:47:39,883 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:47:39,888 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:47:42,770 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:47:42,777 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:47:45,506 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:47:45,511 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:47:48,258 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:47:48,275 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:47:49,343 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:47:49,359 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:48:02,239 (pygeoprocessing.routing.routing) Task._call(1093) INFO 47.2% complete 2023-10-23 14:48:13,268 (pygeoprocessing.routing.routing) Task._call(1093) INFO 62.9% complete 2023-10-23 14:48:24,094 (pygeoprocessing.routing.routing) Task._call(1093) INFO 82.6% complete 2023-10-23 14:48:27,397 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-10-23 14:48:38,174 (pygeoprocessing.routing.routing) Task._call(1093) INFO 51.1% complete 2023-10-23 14:48:49,302 (pygeoprocessing.routing.routing) Task._call(1093) INFO 74.7% complete 2023-10-23 14:48:56,843 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-10-23 14:48:57,194 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:48:57,200 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:48:59,846 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:48:59,855 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:49:00,897 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:49:00,903 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:49:04,718 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:49:04,733 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:49:07,003 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:49:07,017 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:49:09,791 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:49:09,806 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:49:10,791 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:49:10,809 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:49:12,292 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:49:12,297 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:49:13,511 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:49:13,527 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:49:14,361 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:49:14,364 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:49:15,561 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:49:15,563 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:49:16,429 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:49:16,432 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:49:17,446 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:49:17,449 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:49:18,298 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:49:18,302 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:49:20,000 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:49:20,013 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:49:23,187 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:49:23,202 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:50:10,383 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:50:10,386 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:50:13,221 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:50:13,234 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:50:15,521 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-23 14:50:15,525 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-23 14:50:17,025 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-23 14:50:17,027 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-23 14:50:21,010 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 1 aligned: clipped_raster.tif 2023-10-23 14:50:21,013 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 1 rasters. 2023-10-23 14:50:21,030 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1524) INFO build a lookup of aggregate field value to FID 2023-10-23 14:50:21,228 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1531) INFO creating disjoint polygon set 2023-10-23 14:50:21,247 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2554) INFO build shapely polygon list 2023-10-23 14:50:21,295 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2580) INFO build shapely rtree index 2023-10-23 14:50:21,312 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2593) INFO poly feature lookup 100.0% complete on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:21,315 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2597) INFO build poly intersection lookup 2023-10-23 14:50:21,794 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2623) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:21,798 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2655) INFO maximal subset build 100.0% complete on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:21,878 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1551) INFO processing 7 disjoint polygon sets 2023-10-23 14:50:21,898 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 1 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:21,979 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 1 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:22,964 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 1 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:23,647 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 2 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:23,669 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 2 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:24,581 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 2 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:24,831 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 3 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:24,853 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 3 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:25,615 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 3 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:25,816 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 4 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:25,843 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 4 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:26,618 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 4 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:26,800 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 5 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:26,822 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 5 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:27,534 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 5 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:27,702 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1553) INFO zonal stats approximately 85.7% complete on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:27,711 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 6 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:27,724 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 6 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:28,401 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 6 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:28,536 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 7 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:28,557 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 7 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:29,185 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 7 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:29,320 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1756) INFO all done processing polygon sets for watershed_results_ndr_benta.gpkg 2023-10-23 14:50:32,289 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 1 aligned: clipped_raster.tif 2023-10-23 14:50:32,293 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 1 rasters. 2023-10-23 14:50:32,308 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1524) INFO build a lookup of aggregate field value to FID 2023-10-23 14:50:32,325 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1531) INFO creating disjoint polygon set 2023-10-23 14:50:32,348 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2554) INFO build shapely polygon list 2023-10-23 14:50:32,407 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2580) INFO build shapely rtree index 2023-10-23 14:50:32,424 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2593) INFO poly feature lookup 100.0% complete on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:32,427 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2597) INFO build poly intersection lookup 2023-10-23 14:50:32,940 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2623) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:32,948 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2655) INFO maximal subset build 100.0% complete on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:33,007 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1551) INFO processing 7 disjoint polygon sets 2023-10-23 14:50:33,042 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 1 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:33,078 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 1 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:34,007 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 1 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:34,574 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 2 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:34,602 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 2 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:35,542 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 2 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:35,777 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 3 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:35,801 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 3 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:36,609 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 3 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:36,810 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 4 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:36,831 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 4 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:37,594 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 4 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:37,778 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 5 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:37,798 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 5 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:38,511 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 5 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:38,662 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1553) INFO zonal stats approximately 85.7% complete on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:38,681 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 6 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:38,700 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 6 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:39,379 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 6 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:39,528 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 7 of 7 100.0% processed on watershed_results_ndr_benta.gpkg 2023-10-23 14:50:39,548 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 7 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:40,214 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 7 of 7 watershed_results_ndr_benta.gpkg 2023-10-23 14:50:40,398 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1756) INFO all done processing polygon sets for watershed_results_ndr_benta.gpkg 2023-10-23 14:50:40,847 (natcap.invest.ndr.ndr) ndr.execute(984) INFO Writing summaries to output shapefile 2023-10-23 14:50:40,982 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer WFD_Management_Catchments_Cycle_1_square_repair 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-10-23 14:50:46,152 (natcap.invest.ndr.ndr) ndr.execute(988) INFO NDR complete! 2023-10-23 14:50:46,169 (natcap.invest.ndr.ndr) ndr.execute(989) INFO _ _ ____ ____ 2023-10-23 14:50:46,172 (natcap.invest.ndr.ndr) ndr.execute(990) INFO | \ |"| | _"\U | _"\ u 2023-10-23 14:50:46,174 (natcap.invest.ndr.ndr) ndr.execute(991) INFO <| \| |>/| | | |\| |_) |/ 2023-10-23 14:50:46,177 (natcap.invest.ndr.ndr) ndr.execute(992) INFO U| |\ |uU| |_| |\| _ < 2023-10-23 14:50:46,179 (natcap.invest.ndr.ndr) ndr.execute(993) INFO |_| \_| |____/ u|_| \_\ 2023-10-23 14:50:46,181 (natcap.invest.ndr.ndr) ndr.execute(994) INFO || \\,-.|||_ // \\_ 2023-10-23 14:50:46,183 (natcap.invest.ndr.ndr) ndr.execute(995) INFO (_") (_/(__)_) (__) (__) 2023-10-23 14:50:46,184 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 6m 48.70999999999998s 2023-10-23 14:50:46,186 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished