2024-01-12 08:37:10,994 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [C:\Users\Jaye Bee Lopena\Documents\sediment_thesis\2nd-run\2011\results_try_new_version\InVEST-natcap.invest.sdr.sdr-log-2024-01-12--08_37_10.txt] 2024-01-12 08:37:10,994 (__main__) cli.main(459) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.sdr.sdr 3.14.1: biophysical_table_path C:\Users\Jaye Bee Lopena\Documents\sediment_thesis\2nd-run\2011\2011_biophysical_table.csv dem_path C:\Users\Jaye Bee Lopena\Documents\sediment_thesis\2nd-run\2011\DEM\NEWDEM\dem_10m.tif drainage_path erodibility_path C:\Users\Jaye Bee Lopena\Documents\sediment_thesis\2nd-run\2011\erodability\resample_K10.tif erosivity_path C:\Users\Jaye Bee Lopena\Documents\sediment_thesis\2nd-run\2011\erosivity(cell10_edited_final_2011)\resample_Rfinal.tif ic_0_param 0.5 k_param 2 l_max 122 lulc_path C:\Users\Jaye Bee Lopena\Documents\sediment_thesis\2nd-run\2011\2011lulc\LULC_2011_PolygonToRaster1.tif n_workers -1 results_suffix sdr_max 0.8 threshold_flow_accumulation 1000 watersheds_path C:\Users\Jaye Bee Lopena\Documents\sediment_thesis\2nd-run\2011\watershed_2011\2011watershed.shp workspace_dir C:\Users\Jaye Bee Lopena\Documents\sediment_thesis\2nd-run\2011\results_try_new_version 2024-01-12 08:37:13,180 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:13,210 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:13,276 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:13,276 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:13,376 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 1 of 4 aligned: aligned_dem.tif 2024-01-12 08:37:14,255 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:14,267 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:14,280 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:14,280 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:14,347 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 2 of 4 aligned: aligned_lulc.tif 2024-01-12 08:37:14,989 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:15,000 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:15,042 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:15,042 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:15,158 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 3 of 4 aligned: aligned_erosivity.tif 2024-01-12 08:37:15,689 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:15,694 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:15,726 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:15,726 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:15,777 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 4 of 4 aligned: aligned_erodibility.tif 2024-01-12 08:37:15,777 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1154) INFO aligned all 4 rasters. 2024-01-12 08:37:16,046 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:16,048 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:16,170 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:16,170 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:16,474 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:16,479 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:16,556 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:16,556 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:16,872 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:16,879 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:16,923 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:16,923 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:17,222 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:17,238 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:17,321 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:17,321 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:17,622 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:17,624 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:17,688 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:17,688 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:18,569 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete 2024-01-12 08:37:19,689 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:19,705 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:19,770 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:19,770 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:20,935 (pygeoprocessing.routing.routing) Task._call(1093) INFO Flow dir MFD 100.0% complete 2024-01-12 08:37:22,000 (pygeoprocessing.routing.routing) Task._call(1093) INFO Flow accum MFD 100.0% complete 2024-01-12 08:37:22,284 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:22,295 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:22,671 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:22,671 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:23,436 (pygeoprocessing.routing.routing) Task._call(1093) INFO Extract streams MFD: filter out incomplete divergent streams 2024-01-12 08:37:23,468 (pygeoprocessing.routing.routing) Task._call(1093) INFO Extract streams MFD: 100.0% complete 2024-01-12 08:37:23,629 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:23,637 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:23,737 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:23,737 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:23,834 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:23,845 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:23,884 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:23,884 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:24,155 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:24,157 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:24,247 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:24,247 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:24,498 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:24,512 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:24,716 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:24,716 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:25,047 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:25,057 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:25,147 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:25,147 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:26,558 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1361) INFO Flow accum MFD 100.0% complete 2024-01-12 08:37:26,598 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:26,605 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:26,859 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:26,859 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:28,405 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1361) INFO Flow accum MFD 100.0% complete 2024-01-12 08:37:28,445 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:28,465 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:28,717 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:28,717 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:29,276 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:29,281 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:29,580 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:29,580 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:30,102 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:30,114 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:30,225 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:30,225 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:31,387 (pygeoprocessing.routing.routing) Task._call(1093) INFO Dist to channel MFD 100.0% complete 2024-01-12 08:37:31,901 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:31,917 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:32,153 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:32,153 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:32,665 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:32,692 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:32,793 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:32,793 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:33,257 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:33,269 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:33,403 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:33,403 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:33,724 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:33,740 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:33,838 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:33,838 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:34,457 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Calculate sediment deposition 2024-01-12 08:37:35,899 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 100% complete 2024-01-12 08:37:36,280 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:36,291 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:36,394 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:36,394 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:36,941 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:36,957 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:37,100 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:37,100 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:37,487 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-01-12 08:37:37,498 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-01-12 08:37:37,616 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-01-12 08:37:37,616 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-01-12 08:37:37,970 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 20.4316999999999993 of field Area_Ha of feature 0 not successfully written. Possibly due to too larger number with respect to field width 2024-01-12 08:37:38,020 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 366.315999999999974 of field Area_Ha of feature 2 not successfully written. Possibly due to too larger number with respect to field width 2024-01-12 08:37:38,030 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 11.9393999999999991 of field Area_Ha of feature 3 not successfully written. Possibly due to too larger number with respect to field width 2024-01-12 08:37:38,053 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 5318.89999999999964 of field Area_Ha of feature 4 not successfully written. Possibly due to too larger number with respect to field width 2024-01-12 08:37:38,058 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 13.7946000000000009 of field Area_Ha of feature 5 not successfully written. Possibly due to too larger number with respect to field width 2024-01-12 08:37:38,735 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector 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. 2024-01-12 08:37:39,050 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection 2024-01-12 08:37:39,124 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector 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. 2024-01-12 08:37:39,142 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg 2024-01-12 08:37:39,239 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector 2024-01-12 08:37:39,365 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets 2024-01-12 08:37:39,382 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list 2024-01-12 08:37:39,478 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index 2024-01-12 08:37:39,650 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:39,650 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup 2024-01-12 08:37:39,703 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:39,707 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:40,151 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1 2024-01-12 08:37:40,162 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 6 2024-01-12 08:37:40,379 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 1 of 6 2024-01-12 08:37:40,403 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 2 of 6 2024-01-12 08:37:40,454 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 3 of 6 2024-01-12 08:37:40,481 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 4 of 6 2024-01-12 08:37:40,522 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 5 of 6 2024-01-12 08:37:40,557 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_sdr.shp 2024-01-12 08:37:40,638 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector 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. 2024-01-12 08:37:40,873 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection 2024-01-12 08:37:40,900 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector 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. 2024-01-12 08:37:40,912 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg 2024-01-12 08:37:40,983 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector 2024-01-12 08:37:41,392 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets 2024-01-12 08:37:41,401 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list 2024-01-12 08:37:41,412 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index 2024-01-12 08:37:41,414 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:41,416 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup 2024-01-12 08:37:41,425 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:41,425 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:41,688 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1 2024-01-12 08:37:41,707 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 6 2024-01-12 08:37:41,750 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 1 of 6 2024-01-12 08:37:41,771 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 2 of 6 2024-01-12 08:37:41,803 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 3 of 6 2024-01-12 08:37:41,857 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 4 of 6 2024-01-12 08:37:41,906 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 5 of 6 2024-01-12 08:37:41,951 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_sdr.shp 2024-01-12 08:37:42,036 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector 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. 2024-01-12 08:37:42,386 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection 2024-01-12 08:37:42,449 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector 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. 2024-01-12 08:37:42,484 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg 2024-01-12 08:37:42,616 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector 2024-01-12 08:37:42,795 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets 2024-01-12 08:37:42,807 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list 2024-01-12 08:37:42,824 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index 2024-01-12 08:37:42,830 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:42,832 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup 2024-01-12 08:37:42,849 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:42,850 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:43,249 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1 2024-01-12 08:37:43,261 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 6 2024-01-12 08:37:43,299 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 1 of 6 2024-01-12 08:37:43,335 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 2 of 6 2024-01-12 08:37:43,370 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 3 of 6 2024-01-12 08:37:43,405 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 4 of 6 2024-01-12 08:37:43,452 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 5 of 6 2024-01-12 08:37:43,488 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_sdr.shp 2024-01-12 08:37:43,552 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector 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. 2024-01-12 08:37:43,839 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection 2024-01-12 08:37:43,867 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector 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. 2024-01-12 08:37:43,886 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg 2024-01-12 08:37:44,080 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector 2024-01-12 08:37:44,356 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets 2024-01-12 08:37:44,368 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list 2024-01-12 08:37:44,388 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index 2024-01-12 08:37:44,394 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:44,396 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup 2024-01-12 08:37:44,409 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:44,410 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:44,671 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1 2024-01-12 08:37:44,687 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 6 2024-01-12 08:37:44,732 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 1 of 6 2024-01-12 08:37:44,783 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 2 of 6 2024-01-12 08:37:44,821 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 3 of 6 2024-01-12 08:37:44,864 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 4 of 6 2024-01-12 08:37:44,919 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 5 of 6 2024-01-12 08:37:44,965 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_sdr.shp 2024-01-12 08:37:45,066 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector 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. 2024-01-12 08:37:45,579 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection 2024-01-12 08:37:45,588 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector 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. 2024-01-12 08:37:45,596 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg 2024-01-12 08:37:45,648 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector 2024-01-12 08:37:45,812 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets 2024-01-12 08:37:45,821 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list 2024-01-12 08:37:45,831 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index 2024-01-12 08:37:45,833 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:45,833 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup 2024-01-12 08:37:45,845 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:45,847 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2024-01-12 08:37:46,148 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1 2024-01-12 08:37:46,166 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 6 2024-01-12 08:37:46,198 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 1 of 6 2024-01-12 08:37:46,235 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 2 of 6 2024-01-12 08:37:46,282 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 3 of 6 2024-01-12 08:37:46,320 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 4 of 6 2024-01-12 08:37:46,368 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 5 of 6 2024-01-12 08:37:46,408 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_sdr.shp 2024-01-12 08:37:46,464 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 20.4316999999999993 of field Area_Ha of feature 0 not successfully written. Possibly due to too larger number with respect to field width 2024-01-12 08:37:46,488 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 366.315999999999974 of field Area_Ha of feature 2 not successfully written. Possibly due to too larger number with respect to field width 2024-01-12 08:37:46,502 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 11.9393999999999991 of field Area_Ha of feature 3 not successfully written. Possibly due to too larger number with respect to field width 2024-01-12 08:37:46,518 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 5318.89999999999964 of field Area_Ha of feature 4 not successfully written. Possibly due to too larger number with respect to field width 2024-01-12 08:37:46,524 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 13.7946000000000009 of field Area_Ha of feature 5 not successfully written. Possibly due to too larger number with respect to field width 2024-01-12 08:37:46,801 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 35.78s 2024-01-12 08:37:46,801 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished