2022-12-14 16:48:15,909 (natcap.invest.utils) utils.prepare_workspace(162) Level 100 Writing log messages to G:/study/ecosystem service data/NDR/output\InVEST-Nutrient-Delivery-Ratio-log-2022-12-14--16_48_09.txt 2022-12-14 16:48:15,914 (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 G:/study/ecosystem service data/NDR/Biophysical Table.csv calc_n True calc_p True dem_path G:/study/ecosystem service data/SDR/input data/DEM_Pro.tif k_param 2 lulc_path G:/study/ecosystem service data/SDR/LUCC_pro/LUCC2005_Pro.tif n_workers -1 results_suffix 2005 runoff_proxy_path G:/study/ecosystem service data/NDR/input/pre2005.tif subsurface_critical_length_n 200 subsurface_eff_n 0.75 threshold_flow_accumulation 1000 watersheds_path G:/study/ecosystem service data/SDR/input data/YZ3_Project.shp workspace_dir G:/study/ecosystem service data/NDR/output 2022-12-14 16:48:17,144 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1741) INFO starting reprojection 2022-12-14 16:48:17,151 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer YZ3_Project of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2022-12-14 16:48:17,222 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1783) INFO reprojection 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:19,098 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-12-14 16:48:19,098 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-12-14 16:48:19,104 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-12-14 16:48:19,127 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-12-14 16:48:19,141 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-12-14 16:48:19,243 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-12-14 16:48:19,347 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:19,348 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-12-14 16:48:19,478 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:19,479 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:19,498 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 6 disjoint polygon sets 2022-12-14 16:48:19,515 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:19,521 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:19,741 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:19,989 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:19,993 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:20,185 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:20,325 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:20,329 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:20,533 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:20,664 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:20,668 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:20,887 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:21,011 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:21,016 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:21,202 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:21,302 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:21,306 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:21,486 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:21,583 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 76 2022-12-14 16:48:21,584 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-2.5689375e+06 1.0000000e+03 0.0000000e+00 4.4704890e+06 0.0000000e+00 -1.0000000e+03] for ('G:/study/ecosystem service data/NDR/output\\n_subsurface_export_2005.tif', 1) 2022-12-14 16:48:21,584 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 76 2022-12-14 16:48:21,584 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_2005.gpkg 2022-12-14 16:48:24,880 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-12-14 16:48:24,881 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-12-14 16:48:24,885 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-12-14 16:48:24,894 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-12-14 16:48:24,908 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-12-14 16:48:24,931 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-12-14 16:48:24,949 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:24,951 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-12-14 16:48:25,088 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:25,091 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:25,132 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 6 disjoint polygon sets 2022-12-14 16:48:25,148 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:25,154 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:25,354 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:25,659 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:25,664 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:25,854 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:25,969 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:25,973 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:26,180 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:26,310 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:26,314 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:26,522 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:26,680 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:26,686 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:26,904 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:27,030 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:27,037 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:27,230 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:27,331 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 76 2022-12-14 16:48:27,332 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-2.5689375e+06 1.0000000e+03 0.0000000e+00 4.4704890e+06 0.0000000e+00 -1.0000000e+03] for ('G:/study/ecosystem service data/NDR/output\\n_total_export_2005.tif', 1) 2022-12-14 16:48:27,333 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 76 2022-12-14 16:48:27,333 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_2005.gpkg 2022-12-14 16:48:28,634 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-12-14 16:48:28,634 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-12-14 16:48:28,642 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-12-14 16:48:28,654 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-12-14 16:48:28,671 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-12-14 16:48:28,706 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-12-14 16:48:28,753 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:28,756 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-12-14 16:48:28,947 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:28,949 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:28,980 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 6 disjoint polygon sets 2022-12-14 16:48:28,996 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:29,002 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:29,212 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:29,477 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:29,482 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:29,682 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:29,809 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:29,813 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:30,027 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:30,167 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:30,172 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:30,379 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:30,507 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:30,511 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:30,700 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:30,804 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:30,809 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:30,993 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:31,098 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 76 2022-12-14 16:48:31,099 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-2.5689375e+06 1.0000000e+03 0.0000000e+00 4.4704890e+06 0.0000000e+00 -1.0000000e+03] for ('G:/study/ecosystem service data/NDR/output\\intermediate_outputs\\sub_load_n_2005.tif', 1) 2022-12-14 16:48:31,099 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 76 2022-12-14 16:48:31,100 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_2005.gpkg 2022-12-14 16:48:32,733 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-12-14 16:48:32,734 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-12-14 16:48:32,739 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-12-14 16:48:32,748 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-12-14 16:48:32,766 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-12-14 16:48:32,792 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-12-14 16:48:32,810 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:32,811 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-12-14 16:48:32,940 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:32,942 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:33,064 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 6 disjoint polygon sets 2022-12-14 16:48:33,080 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:33,085 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:33,284 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:33,685 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:33,689 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:33,879 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:33,999 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:34,003 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:34,210 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:34,346 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:34,349 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:34,558 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:34,688 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:34,691 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:34,881 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:34,988 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:34,992 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:35,185 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:35,281 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 76 2022-12-14 16:48:35,282 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-2.5689375e+06 1.0000000e+03 0.0000000e+00 4.4704890e+06 0.0000000e+00 -1.0000000e+03] for ('G:/study/ecosystem service data/NDR/output\\n_surface_export_2005.tif', 1) 2022-12-14 16:48:35,282 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 76 2022-12-14 16:48:35,283 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_2005.gpkg 2022-12-14 16:48:38,372 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-12-14 16:48:38,373 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-12-14 16:48:38,377 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-12-14 16:48:38,387 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-12-14 16:48:38,401 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-12-14 16:48:38,427 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-12-14 16:48:38,446 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:38,447 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-12-14 16:48:38,580 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:38,582 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:38,626 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 6 disjoint polygon sets 2022-12-14 16:48:38,642 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:38,647 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:38,850 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:39,236 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:39,240 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:39,428 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:39,545 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:39,549 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:39,761 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:39,897 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:39,902 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:40,109 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:40,231 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:40,236 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:40,419 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:40,522 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:40,526 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:40,717 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:40,838 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 76 2022-12-14 16:48:40,839 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-2.5689375e+06 1.0000000e+03 0.0000000e+00 4.4704890e+06 0.0000000e+00 -1.0000000e+03] for ('G:/study/ecosystem service data/NDR/output\\intermediate_outputs\\surface_load_n_2005.tif', 1) 2022-12-14 16:48:40,839 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 76 2022-12-14 16:48:40,840 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_2005.gpkg 2022-12-14 16:48:44,519 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-12-14 16:48:44,520 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-12-14 16:48:44,525 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-12-14 16:48:44,536 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-12-14 16:48:44,551 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-12-14 16:48:44,578 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-12-14 16:48:44,597 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:44,601 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-12-14 16:48:44,743 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:44,745 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:44,795 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 6 disjoint polygon sets 2022-12-14 16:48:44,810 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:44,815 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:45,012 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:45,301 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:45,306 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:45,492 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:45,611 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:45,617 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:45,835 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:45,971 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:45,975 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:46,183 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:46,311 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:46,316 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:46,501 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:46,607 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:46,611 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:46,795 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:46,904 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 76 2022-12-14 16:48:46,904 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-2.5689375e+06 1.0000000e+03 0.0000000e+00 4.4704890e+06 0.0000000e+00 -1.0000000e+03] for ('G:/study/ecosystem service data/NDR/output\\p_surface_export_2005.tif', 1) 2022-12-14 16:48:46,905 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 76 2022-12-14 16:48:46,905 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_2005.gpkg 2022-12-14 16:48:50,067 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-12-14 16:48:50,068 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-12-14 16:48:50,072 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-12-14 16:48:50,083 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-12-14 16:48:50,096 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-12-14 16:48:50,126 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-12-14 16:48:50,147 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:50,148 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-12-14 16:48:50,279 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:50,281 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:50,456 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 6 disjoint polygon sets 2022-12-14 16:48:50,473 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:50,478 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:50,684 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:50,977 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:50,982 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:51,175 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:51,293 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:51,297 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:51,512 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:51,653 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:51,657 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:51,882 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:52,008 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:52,012 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:52,196 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:52,301 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 6 100.0% processed on watershed_results_ndr_2005.gpkg 2022-12-14 16:48:52,306 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:52,494 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 6 watershed_results_ndr_2005.gpkg 2022-12-14 16:48:52,591 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 76 2022-12-14 16:48:52,591 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [-2.5689375e+06 1.0000000e+03 0.0000000e+00 4.4704890e+06 0.0000000e+00 -1.0000000e+03] for ('G:/study/ecosystem service data/NDR/output\\intermediate_outputs\\surface_load_p_2005.tif', 1) 2022-12-14 16:48:52,593 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 76 2022-12-14 16:48:52,593 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_2005.gpkg 2022-12-14 16:48:52,661 (natcap.invest.ndr.ndr) ndr.execute(730) INFO Writing summaries to output shapefile 2022-12-14 16:48:52,765 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer YZ3_Project of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2022-12-14 16:48:57,762 (natcap.invest.ndr.ndr) ndr.execute(734) INFO NDR complete! 2022-12-14 16:48:57,763 (natcap.invest.ndr.ndr) ndr.execute(735) INFO _ _ ____ ____ 2022-12-14 16:48:57,763 (natcap.invest.ndr.ndr) ndr.execute(736) INFO | \ |"| | _"\U | _"\ u 2022-12-14 16:48:57,764 (natcap.invest.ndr.ndr) ndr.execute(737) INFO <| \| |>/| | | |\| |_) |/ 2022-12-14 16:48:57,764 (natcap.invest.ndr.ndr) ndr.execute(738) INFO U| |\ |uU| |_| |\| _ < 2022-12-14 16:48:57,765 (natcap.invest.ndr.ndr) ndr.execute(739) INFO |_| \_| |____/ u|_| \_\ 2022-12-14 16:48:57,766 (natcap.invest.ndr.ndr) ndr.execute(740) INFO || \\,-.|||_ // \\_ 2022-12-14 16:48:57,767 (natcap.invest.ndr.ndr) ndr.execute(741) INFO (_") (_/(__)_) (__) (__) 2022-12-14 16:48:57,770 (natcap.invest.utils) utils.prepare_workspace(170) INFO Elapsed time: 41.86s 2022-12-14 16:48:57,770 (natcap.invest.utils) utils.prepare_workspace(173) INFO Execution finished