2023-09-25 15:13:25,688 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [C:\Users\Eri\Documents\DOUTORADO\Projeto Ipanema Emissão\SHAPES\CONQUISTA\2016_InVEST\25-09 (2)\InVEST-natcap.invest.sdr.sdr-log-2023-09-25--15_13_25.txt] 2023-09-25 15:13:25,688 (__main__) cli.main(458) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.sdr.sdr 3.13.0: biophysical_table_path C:\Users\Eri\Documents\DOUTORADO\Projeto Ipanema Emissão\SHAPES\CONQUISTA\Dados_Invest_2016\biophysical_table_Gura2016.csv dem_path C:\Users\Eri\Documents\DOUTORADO\Projeto Ipanema Emissão\SHAPES\CONQUISTA\DEM_SRTM_30m_23S.tif drainage_path erodibility_path C:\Users\Eri\Documents\DOUTORADO\Projeto Ipanema Emissão\SHAPES\CONQUISTA\K_factor_23S.tif erosivity_path C:\Users\Eri\Documents\DOUTORADO\Projeto Ipanema Emissão\SHAPES\CONQUISTA\R_factor.tif ic_0_param 0.5 k_param 2 l_max 122 lulc_path C:\Users\Eri\Documents\DOUTORADO\Projeto Ipanema Emissão\SHAPES\CONQUISTA\Dados_Invest_2016\Watershead_2016_OK.tif n_workers -1 results_suffix sdr_max 0.8 threshold_flow_accumulation 1000 watersheds_path C:\Users\Eri\Documents\DOUTORADO\Projeto Ipanema Emissão\SHAPES\CONQUISTA\Dados_Invest_2016\Waterhead_OK_23s.shp workspace_dir C:\Users\Eri\Documents\DOUTORADO\Projeto Ipanema Emissão\SHAPES\CONQUISTA\2016_InVEST\25-09 (2) 2023-09-25 15:13:26,079 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:26,085 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:26,090 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:26,091 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:26,096 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 4 aligned: aligned_dem.tif 2023-09-25 15:13:26,152 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:26,154 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:26,155 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:26,155 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:26,157 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 2 of 4 aligned: aligned_lulc.tif 2023-09-25 15:13:26,270 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:26,272 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:26,274 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:26,274 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:26,279 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 3 of 4 aligned: aligned_erosivity.tif 2023-09-25 15:13:26,369 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:26,371 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:26,374 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:26,374 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:26,382 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 4 of 4 aligned: aligned_erodibility.tif 2023-09-25 15:13:26,382 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 4 rasters. 2023-09-25 15:13:26,454 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete 2023-09-25 15:13:26,594 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:26,598 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:26,601 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:26,601 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:26,709 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-09-25 15:13:26,727 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Calculating average aspect 2023-09-25 15:13:26,743 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 0.00% complete 2023-09-25 15:13:26,749 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 66.32% complete 2023-09-25 15:13:26,753 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 100.00% complete 2023-09-25 15:13:26,835 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-09-25 15:13:26,878 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:26,903 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:26,924 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:26,924 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:26,994 (pygeoprocessing.routing.routing) Task._call(1093) INFO filter out incomplete divergent streams 2023-09-25 15:13:26,997 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-09-25 15:13:27,060 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:27,065 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:27,073 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:27,073 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:27,086 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:27,089 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:27,090 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:27,090 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:27,128 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:27,134 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:27,137 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:27,137 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:27,204 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:27,206 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:27,211 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:27,211 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:27,253 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:27,254 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:27,258 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:27,259 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:27,386 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1407) INFO 100.0% complete 2023-09-25 15:13:27,402 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:27,403 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:27,412 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:27,412 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:27,528 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1407) INFO 100.0% complete 2023-09-25 15:13:27,542 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:27,545 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:27,553 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:27,554 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:27,656 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:27,658 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:27,665 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:27,665 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:27,720 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:27,927 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:27,934 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:27,934 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:28,019 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-09-25 15:13:28,050 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:28,052 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:28,056 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:28,056 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:28,088 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:28,090 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:28,092 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:28,092 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:28,131 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:28,131 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:28,137 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:28,138 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:28,174 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:28,177 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:28,179 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:28,179 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:28,206 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Calculate sediment deposition 2023-09-25 15:13:28,244 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 0.00% complete 2023-09-25 15:13:28,260 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 0.00% complete 2023-09-25 15:13:28,261 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 100% complete 2023-09-25 15:13:28,308 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:28,311 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:28,319 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:28,320 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:28,370 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:28,372 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:28,377 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:28,377 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:28,423 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-25 15:13:28,424 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-25 15:13:28,428 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-25 15:13:28,428 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-25 15:13:28,489 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 15828860.8684 of field Shape_Area of feature 0 not successfully written. Possibly due to too larger number with respect to field width 2023-09-25 15:13:28,534 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 1 aligned: clipped_raster.tif 2023-09-25 15:13:28,534 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 1 rasters. 2023-09-25 15:13:28,546 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1524) INFO build a lookup of aggregate field value to FID 2023-09-25 15:13:28,553 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1531) INFO creating disjoint polygon set 2023-09-25 15:13:28,554 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2554) INFO build shapely polygon list 2023-09-25 15:13:28,585 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2580) INFO build shapely rtree index 2023-09-25 15:13:28,604 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2593) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:28,604 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2597) INFO build poly intersection lookup 2023-09-25 15:13:28,620 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2623) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:28,621 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2655) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:28,627 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1551) INFO processing 5 disjoint polygon sets 2023-09-25 15:13:28,632 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,635 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,641 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,649 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,650 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,654 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,656 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,656 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,659 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,660 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,660 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,664 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,669 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,669 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,671 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,672 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1756) INFO all done processing polygon sets for watershed_results_sdr.shp 2023-09-25 15:13:28,703 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 1 aligned: clipped_raster.tif 2023-09-25 15:13:28,703 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 1 rasters. 2023-09-25 15:13:28,712 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1524) INFO build a lookup of aggregate field value to FID 2023-09-25 15:13:28,715 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1531) INFO creating disjoint polygon set 2023-09-25 15:13:28,716 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2554) INFO build shapely polygon list 2023-09-25 15:13:28,728 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2580) INFO build shapely rtree index 2023-09-25 15:13:28,732 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2593) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:28,733 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2597) INFO build poly intersection lookup 2023-09-25 15:13:28,741 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2623) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:28,741 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2655) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:28,753 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1551) INFO processing 5 disjoint polygon sets 2023-09-25 15:13:28,754 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,754 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,758 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,760 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,761 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,769 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,772 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,772 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,775 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,777 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,777 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,779 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,783 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,785 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,790 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,791 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1756) INFO all done processing polygon sets for watershed_results_sdr.shp 2023-09-25 15:13:28,828 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 1 aligned: clipped_raster.tif 2023-09-25 15:13:28,830 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 1 rasters. 2023-09-25 15:13:28,841 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1524) INFO build a lookup of aggregate field value to FID 2023-09-25 15:13:28,843 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1531) INFO creating disjoint polygon set 2023-09-25 15:13:28,845 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2554) INFO build shapely polygon list 2023-09-25 15:13:28,861 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2580) INFO build shapely rtree index 2023-09-25 15:13:28,862 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2593) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:28,862 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2597) INFO build poly intersection lookup 2023-09-25 15:13:28,873 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2623) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:28,874 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2655) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:28,880 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1551) INFO processing 5 disjoint polygon sets 2023-09-25 15:13:28,884 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,885 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,890 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,893 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,893 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,899 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,903 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,904 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,906 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,909 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,909 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,911 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,914 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,915 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,918 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_sdr.shp 2023-09-25 15:13:28,921 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1756) INFO all done processing polygon sets for watershed_results_sdr.shp 2023-09-25 15:13:28,952 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 1 aligned: clipped_raster.tif 2023-09-25 15:13:28,952 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 1 rasters. 2023-09-25 15:13:28,962 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1524) INFO build a lookup of aggregate field value to FID 2023-09-25 15:13:28,965 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1531) INFO creating disjoint polygon set 2023-09-25 15:13:28,967 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2554) INFO build shapely polygon list 2023-09-25 15:13:28,979 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2580) INFO build shapely rtree index 2023-09-25 15:13:28,984 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2593) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:28,984 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2597) INFO build poly intersection lookup 2023-09-25 15:13:28,990 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2623) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:28,990 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2655) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:28,995 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1551) INFO processing 5 disjoint polygon sets 2023-09-25 15:13:28,997 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:28,998 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,004 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,007 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:29,007 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,011 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,014 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:29,016 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,021 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,023 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:29,023 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,025 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,026 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:29,026 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,029 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,033 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1756) INFO all done processing polygon sets for watershed_results_sdr.shp 2023-09-25 15:13:29,067 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 1 aligned: clipped_raster.tif 2023-09-25 15:13:29,067 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 1 rasters. 2023-09-25 15:13:29,076 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1524) INFO build a lookup of aggregate field value to FID 2023-09-25 15:13:29,077 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1531) INFO creating disjoint polygon set 2023-09-25 15:13:29,078 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2554) INFO build shapely polygon list 2023-09-25 15:13:29,090 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2580) INFO build shapely rtree index 2023-09-25 15:13:29,091 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2593) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:29,091 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2597) INFO build poly intersection lookup 2023-09-25 15:13:29,095 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2623) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:29,096 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2655) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2023-09-25 15:13:29,106 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1551) INFO processing 5 disjoint polygon sets 2023-09-25 15:13:29,106 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 1 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:29,107 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 1 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,110 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 1 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,113 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 2 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:29,113 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 2 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,119 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 2 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,124 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 3 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:29,124 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 3 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,126 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 3 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,128 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 4 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:29,128 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 4 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,133 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 4 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,136 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 5 of 5 100.0% processed on watershed_results_sdr.shp 2023-09-25 15:13:29,136 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 5 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,139 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 5 of 5 watershed_results_sdr.shp 2023-09-25 15:13:29,140 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1756) INFO all done processing polygon sets for watershed_results_sdr.shp 2023-09-25 15:13:29,146 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 15828860.8684 of field Shape_Area of feature 0 not successfully written. Possibly due to too larger number with respect to field width 2023-09-25 15:13:29,170 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 3.48s 2023-09-25 15:13:29,170 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished