2022-05-12 22:19:39,393 (natcap.invest.utils) utils.prepare_workspace(162) Level 100 Writing log messages to /Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/InVEST-Workspaces/NDR/CJ-NDR-5.12.22/InVEST-Nutrient-Delivery-Ratio-log-2022-05-12--22_19_39.txt 2022-05-12 22:19:39,393 (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 /Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/Data_Inputs/NDR_INVEST_INPUTS/biophysical_table_NDR_v3.csv calc_n True calc_p False dem_path /Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/Data_Inputs/NDR_INVEST_INPUTS/DEM-EPSG-26986.tif k_param 2 lulc_path /Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/Data_Inputs/SCENARIOS/LULC-scenario-files/cran-wetland/LULC-cran-wetland-100-0.tif n_workers -1 results_suffix _cran_wetland_100_0 runoff_proxy_path /Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/Data_Inputs/NDR_INVEST_INPUTS/precipitation-EPSG-26986.tif subsurface_critical_length_n 150 subsurface_eff_n .05 threshold_flow_accumulation 1000 watersheds_path /Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/Data_Inputs/Helper_Data_Folders/Watersheds/NDR-calibration-watershed/WARRIV_SUBW_16JUL10.shp workspace_dir /Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/InVEST-Workspaces/NDR/CJ-NDR-5.12.22 2022-05-12 22:19:39,413 (natcap.invest.utils) utils.build_lookup_from_csv(570) WARNING Empty or NaN values were found in the table: /Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/Data_Inputs/NDR_INVEST_INPUTS/biophysical_table_NDR_v3.csv. 2022-05-12 22:19:39,521 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1741) INFO starting reprojection 2022-05-12 22:19:39,537 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer WARRIV_SUBW_16JUL10 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-05-12 22:19:39,539 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1783) INFO reprojection 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:39,816 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-05-12 22:19:39,817 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-05-12 22:19:39,819 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-05-12 22:19:39,820 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-05-12 22:19:39,825 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-05-12 22:19:39,852 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-05-12 22:19:39,860 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:39,860 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-05-12 22:19:39,889 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:39,890 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:39,900 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 7 disjoint polygon sets 2022-05-12 22:19:39,906 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:39,909 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:39,940 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:39,985 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:39,986 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:39,995 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,020 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,038 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,065 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,076 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,077 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,087 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,093 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,094 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,103 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,110 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,110 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,119 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,128 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 7 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,129 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 7 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,152 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 7 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,158 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 50 2022-05-12 22:19:40,178 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 2.6227238e+05 2.5679913e+01 0.0000000e+00 8.4849338e+05 0.0000000e+00 -2.5679913e+01] for ('/Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/InVEST-Workspaces/NDR/CJ-NDR-5.12.22/n_subsurface_export_cran_wetland_100_0.tif', 1) 2022-05-12 22:19:40,180 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 50 2022-05-12 22:19:40,198 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,255 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-05-12 22:19:40,256 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-05-12 22:19:40,257 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-05-12 22:19:40,259 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-05-12 22:19:40,262 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-05-12 22:19:40,266 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-05-12 22:19:40,271 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,271 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-05-12 22:19:40,288 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,288 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,293 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 7 disjoint polygon sets 2022-05-12 22:19:40,297 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,297 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,307 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,323 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,323 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,331 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,338 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,338 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,345 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,352 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,352 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,361 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,368 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,369 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,394 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,400 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,401 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,408 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,414 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 7 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,414 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 7 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,421 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 7 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,425 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 50 2022-05-12 22:19:40,426 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 2.6227238e+05 2.5679913e+01 0.0000000e+00 8.4849338e+05 0.0000000e+00 -2.5679913e+01] for ('/Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/InVEST-Workspaces/NDR/CJ-NDR-5.12.22/n_total_export_cran_wetland_100_0.tif', 1) 2022-05-12 22:19:40,426 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 50 2022-05-12 22:19:40,426 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,492 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-05-12 22:19:40,492 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-05-12 22:19:40,494 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-05-12 22:19:40,495 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-05-12 22:19:40,499 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-05-12 22:19:40,503 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-05-12 22:19:40,509 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,509 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-05-12 22:19:40,529 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,530 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,536 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 7 disjoint polygon sets 2022-05-12 22:19:40,539 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,539 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,552 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,568 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,569 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,577 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,588 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,589 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,598 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,607 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,608 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,619 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,627 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,628 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,639 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,647 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,648 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,657 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,664 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 7 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,665 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 7 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,673 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 7 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,679 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 50 2022-05-12 22:19:40,679 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 2.6227238e+05 2.5679913e+01 0.0000000e+00 8.4849338e+05 0.0000000e+00 -2.5679913e+01] for ('/Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/InVEST-Workspaces/NDR/CJ-NDR-5.12.22/intermediate_outputs/sub_load_n_cran_wetland_100_0.tif', 1) 2022-05-12 22:19:40,680 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 50 2022-05-12 22:19:40,680 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,745 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-05-12 22:19:40,746 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-05-12 22:19:40,749 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-05-12 22:19:40,750 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-05-12 22:19:40,754 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-05-12 22:19:40,759 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-05-12 22:19:40,765 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,765 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-05-12 22:19:40,790 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,791 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,798 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 7 disjoint polygon sets 2022-05-12 22:19:40,803 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,804 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,815 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,837 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,838 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,846 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,858 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,858 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,868 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,878 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,879 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,888 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,895 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,896 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,904 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,911 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,911 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,919 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,924 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 7 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,924 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 7 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,930 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 7 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,934 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 50 2022-05-12 22:19:40,935 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 2.6227238e+05 2.5679913e+01 0.0000000e+00 8.4849338e+05 0.0000000e+00 -2.5679913e+01] for ('/Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/InVEST-Workspaces/NDR/CJ-NDR-5.12.22/n_surface_export_cran_wetland_100_0.tif', 1) 2022-05-12 22:19:40,935 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 50 2022-05-12 22:19:40,935 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:40,986 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-05-12 22:19:40,986 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-05-12 22:19:40,988 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-05-12 22:19:40,989 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-05-12 22:19:40,993 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-05-12 22:19:40,996 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-05-12 22:19:41,003 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,003 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-05-12 22:19:41,022 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,026 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,050 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 7 disjoint polygon sets 2022-05-12 22:19:41,053 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,059 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,089 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,104 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 2 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,105 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 2 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,113 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 2 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,122 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 3 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,123 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 3 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,131 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 3 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,138 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 4 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,138 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 4 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,147 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 4 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,172 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 5 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,173 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 5 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,183 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 5 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,189 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 6 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,190 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 6 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,197 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 6 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,202 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 7 of 7 100.0% processed on watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,203 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 7 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,209 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 7 of 7 watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,213 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 50 2022-05-12 22:19:41,213 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 2.6227238e+05 2.5679913e+01 0.0000000e+00 8.4849338e+05 0.0000000e+00 -2.5679913e+01] for ('/Users/carolinejaffe/Dropbox (MIT)/UROP_Share_Folder_Dropbox/InVEST-Workspaces/NDR/CJ-NDR-5.12.22/intermediate_outputs/surface_load_n_cran_wetland_100_0.tif', 1) 2022-05-12 22:19:41,214 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 50 2022-05-12 22:19:41,214 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_ndr_cran_wetland_100_0.gpkg 2022-05-12 22:19:41,220 (natcap.invest.ndr.ndr) ndr.execute(730) INFO Writing summaries to output shapefile 2022-05-12 22:19:41,289 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer WARRIV_SUBW_16JUL10 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-05-12 22:19:41,306 (natcap.invest.ndr.ndr) ndr.execute(734) INFO NDR complete! 2022-05-12 22:19:41,307 (natcap.invest.ndr.ndr) ndr.execute(735) INFO _ _ ____ ____ 2022-05-12 22:19:41,307 (natcap.invest.ndr.ndr) ndr.execute(736) INFO | \ |"| | _"\U | _"\ u 2022-05-12 22:19:41,308 (natcap.invest.ndr.ndr) ndr.execute(737) INFO <| \| |>/| | | |\| |_) |/ 2022-05-12 22:19:41,308 (natcap.invest.ndr.ndr) ndr.execute(738) INFO U| |\ |uU| |_| |\| _ < 2022-05-12 22:19:41,308 (natcap.invest.ndr.ndr) ndr.execute(739) INFO |_| \_| |____/ u|_| \_\ 2022-05-12 22:19:41,309 (natcap.invest.ndr.ndr) ndr.execute(740) INFO || \\,-.|||_ // \\_ 2022-05-12 22:19:41,309 (natcap.invest.ndr.ndr) ndr.execute(741) INFO (_") (_/(__)_) (__) (__) 2022-05-12 22:19:41,310 (natcap.invest.utils) utils.prepare_workspace(170) INFO Elapsed time: 1.92s 2022-05-12 22:19:41,310 (natcap.invest.utils) utils.prepare_workspace(173) INFO Execution finished