2022-05-19 14:13:12,069 (natcap.invest.utils) utils.prepare_workspace(162) Level 100 Writing log messages to /Users/joaosantos/Documents/sdr_workspace/InVEST-Sediment-Delivery-Ratio-log-2022-05-19--14_13_12.txt 2022-05-19 14:13:12,070 (natcap.invest.ui.model) model._logged_target(1643) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.sdr.sdr 3.10.2: biophysical_table_path /Users/joaosantos/Desktop/SDR/biophysical_table_Gura.csv dem_path /Users/joaosantos/Desktop/SDR/DEM_gura.tif drainage_path erodibility_path /Users/joaosantos/Desktop/SDR/erodibility_gura.tif erosivity_path /Users/joaosantos/Desktop/SDR/erosivity_gura.tif ic_0_param 0.5 k_param 2 l_max 200 lulc_path /Users/joaosantos/Desktop/SDR/land_use_gura.tif n_workers -1 results_suffix sdr_max 0.8 threshold_flow_accumulation 1000 watersheds_path /Users/joaosantos/Desktop/SDR/watershed_gura.shp workspace_dir /Users/joaosantos/Documents/sdr_workspace 2022-05-19 14:13:12,298 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:12,300 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:12,300 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:12,333 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:12,333 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:12,363 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 4 aligned: aligned_dem.tif 2022-05-19 14:13:12,510 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:12,510 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:12,510 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:12,540 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:12,541 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:12,571 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 2 of 4 aligned: aligned_lulc.tif 2022-05-19 14:13:12,719 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:12,720 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:12,720 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:12,756 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:12,756 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:12,809 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 3 of 4 aligned: aligned_erosivity.tif 2022-05-19 14:13:12,958 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:12,958 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:12,958 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:12,990 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:12,990 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:13,022 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 4 of 4 aligned: aligned_erodibility.tif 2022-05-19 14:13:13,023 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 4 rasters. 2022-05-19 14:13:13,213 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete 2022-05-19 14:13:13,311 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:13,312 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:13,312 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:13,358 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:13,358 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:13,934 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-05-19 14:13:13,938 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Calculating average aspect 2022-05-19 14:13:13,983 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 0.00% complete 2022-05-19 14:13:13,989 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 5.61% complete 2022-05-19 14:13:13,996 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 11.21% complete 2022-05-19 14:13:14,003 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 16.82% complete 2022-05-19 14:13:14,009 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 22.42% complete 2022-05-19 14:13:14,012 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 28.03% complete 2022-05-19 14:13:14,018 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 33.63% complete 2022-05-19 14:13:14,024 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 39.24% complete 2022-05-19 14:13:14,030 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 42.45% complete 2022-05-19 14:13:14,036 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 48.06% complete 2022-05-19 14:13:14,043 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 53.66% complete 2022-05-19 14:13:14,053 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 59.27% complete 2022-05-19 14:13:14,063 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 64.87% complete 2022-05-19 14:13:14,072 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 70.48% complete 2022-05-19 14:13:14,080 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 76.09% complete 2022-05-19 14:13:14,089 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 81.69% complete 2022-05-19 14:13:14,094 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 84.91% complete 2022-05-19 14:13:14,096 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 86.90% complete 2022-05-19 14:13:14,098 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 88.89% complete 2022-05-19 14:13:14,101 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 90.89% complete 2022-05-19 14:13:14,104 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 92.88% complete 2022-05-19 14:13:14,108 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 94.87% complete 2022-05-19 14:13:14,113 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 96.86% complete 2022-05-19 14:13:14,115 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 98.86% complete 2022-05-19 14:13:14,116 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 100.00% complete 2022-05-19 14:13:14,164 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating target flow accum raster layer 2022-05-19 14:13:14,173 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG creating visited raster layer 2022-05-19 14:13:14,198 (pygeoprocessing.routing.routing) Task._call(1093) DEBUG starting search 2022-05-19 14:13:14,740 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-05-19 14:13:14,758 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:14,759 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:14,759 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:14,936 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:14,936 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:15,154 (pygeoprocessing.routing.routing) Task._call(1093) INFO filter out incomplete divergent streams 2022-05-19 14:13:15,164 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-05-19 14:13:15,177 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:15,178 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:15,178 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:15,219 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:15,220 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:15,250 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:15,251 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:15,251 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:15,273 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:15,274 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:15,310 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:15,311 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:15,311 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:15,347 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:15,348 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:15,395 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:15,395 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:15,396 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:15,470 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:15,470 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:15,535 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:15,536 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:15,536 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:15,589 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:15,589 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:15,654 (natcap.invest.sdr.sdr) sdr._calculate_bar_factor(1101) DEBUG doing flow accumulation mfd on /Users/joaosantos/Documents/sdr_workspace/intermediate_outputs/w_threshold.tif 2022-05-19 14:13:15,654 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1104) DEBUG creating target flow accum raster layer 2022-05-19 14:13:15,663 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1104) DEBUG creating visited raster layer 2022-05-19 14:13:15,672 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1104) DEBUG starting search 2022-05-19 14:13:16,251 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1104) INFO 100.0% complete 2022-05-19 14:13:16,258 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:16,258 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:16,258 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:16,367 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:16,367 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:16,416 (natcap.invest.sdr.sdr) sdr._calculate_bar_factor(1101) DEBUG doing flow accumulation mfd on /Users/joaosantos/Documents/sdr_workspace/intermediate_outputs/slope_threshold.tif 2022-05-19 14:13:16,416 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1104) DEBUG creating target flow accum raster layer 2022-05-19 14:13:16,424 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1104) DEBUG creating visited raster layer 2022-05-19 14:13:16,434 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1104) DEBUG starting search 2022-05-19 14:13:17,011 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1104) INFO 100.0% complete 2022-05-19 14:13:17,018 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:17,019 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:17,019 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:17,128 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:17,129 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:17,197 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:17,197 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:17,197 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:17,299 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:17,299 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:17,366 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:17,366 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:17,366 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:17,411 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:17,411 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:17,890 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-05-19 14:13:17,902 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:17,903 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:17,903 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:17,990 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:17,991 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:18,049 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:18,050 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:18,050 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:18,094 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:18,094 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:18,154 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:18,155 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:18,155 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:18,213 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:18,213 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:18,277 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:18,278 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:18,278 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:18,333 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:18,333 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:18,389 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Calculate sediment deposition 2022-05-19 14:13:18,409 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 0.00% complete 2022-05-19 14:13:18,450 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 0.00% complete 2022-05-19 14:13:18,481 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 0.00% complete 2022-05-19 14:13:18,516 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 0.00% complete 2022-05-19 14:13:18,537 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 0.00% complete 2022-05-19 14:13:18,539 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 0.00% complete 2022-05-19 14:13:18,565 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 0.00% complete 2022-05-19 14:13:18,600 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 0.00% complete 2022-05-19 14:13:18,621 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 0.00% complete 2022-05-19 14:13:18,640 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 5.61% complete 2022-05-19 14:13:18,658 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 11.21% complete 2022-05-19 14:13:18,706 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 16.82% complete 2022-05-19 14:13:18,780 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 22.42% complete 2022-05-19 14:13:18,848 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 28.03% complete 2022-05-19 14:13:18,905 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 33.63% complete 2022-05-19 14:13:18,941 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 39.24% complete 2022-05-19 14:13:18,948 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 0.00% complete 2022-05-19 14:13:18,950 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 11.21% complete 2022-05-19 14:13:18,952 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 22.42% complete 2022-05-19 14:13:18,959 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 33.63% complete 2022-05-19 14:13:18,985 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 44.84% complete 2022-05-19 14:13:19,002 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 56.05% complete 2022-05-19 14:13:19,008 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 67.26% complete 2022-05-19 14:13:19,010 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 78.47% complete 2022-05-19 14:13:19,011 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 100% complete 2022-05-19 14:13:19,139 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:19,140 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:19,140 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:19,224 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:19,225 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:19,286 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:19,287 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:19,287 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:19,320 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:19,321 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:19,790 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2022-05-19 14:13:19,804 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:19,805 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:19,805 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:19,887 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:19,888 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:19,953 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:19,953 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:19,953 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:20,037 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:20,037 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:20,097 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:20,097 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:20,098 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:20,139 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:20,139 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:20,204 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:20,204 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:20,204 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:20,315 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:20,315 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:20,384 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-05-19 14:13:20,384 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 25806 2022-05-19 14:13:20,385 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-05-19 14:13:20,452 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-05-19 14:13:20,453 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-05-19 14:13:20,569 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-05-19 14:13:20,569 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-05-19 14:13:20,570 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-05-19 14:13:20,570 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-05-19 14:13:20,571 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-05-19 14:13:20,573 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-05-19 14:13:20,573 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2022-05-19 14:13:20,573 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-05-19 14:13:20,574 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2022-05-19 14:13:20,574 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2022-05-19 14:13:20,579 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-05-19 14:13:20,579 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on watershed_results_sdr.shp 2022-05-19 14:13:20,582 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 watershed_results_sdr.shp 2022-05-19 14:13:20,609 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 watershed_results_sdr.shp 2022-05-19 14:13:20,650 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 1 2022-05-19 14:13:20,650 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 2.4895066e+05 1.5000000e+01 0.0000000e+00 9.9418970e+06 0.0000000e+00 -1.5000000e+01] for ('/Users/joaosantos/Documents/sdr_workspace/usle.tif', 1) 2022-05-19 14:13:20,651 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 1 2022-05-19 14:13:20,651 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_sdr.shp 2022-05-19 14:13:20,759 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-05-19 14:13:20,759 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-05-19 14:13:20,760 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-05-19 14:13:20,760 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-05-19 14:13:20,761 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-05-19 14:13:20,762 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-05-19 14:13:20,762 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2022-05-19 14:13:20,762 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-05-19 14:13:20,762 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2022-05-19 14:13:20,763 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2022-05-19 14:13:20,767 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-05-19 14:13:20,768 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on watershed_results_sdr.shp 2022-05-19 14:13:20,770 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 watershed_results_sdr.shp 2022-05-19 14:13:20,797 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 watershed_results_sdr.shp 2022-05-19 14:13:20,837 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 1 2022-05-19 14:13:20,837 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 2.4895066e+05 1.5000000e+01 0.0000000e+00 9.9418970e+06 0.0000000e+00 -1.5000000e+01] for ('/Users/joaosantos/Documents/sdr_workspace/sed_export.tif', 1) 2022-05-19 14:13:20,837 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 1 2022-05-19 14:13:20,837 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_sdr.shp 2022-05-19 14:13:20,975 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-05-19 14:13:20,975 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-05-19 14:13:20,976 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-05-19 14:13:20,976 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-05-19 14:13:20,977 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-05-19 14:13:20,977 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-05-19 14:13:20,978 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2022-05-19 14:13:20,978 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-05-19 14:13:20,978 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2022-05-19 14:13:20,978 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2022-05-19 14:13:20,983 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-05-19 14:13:20,983 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on watershed_results_sdr.shp 2022-05-19 14:13:20,985 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 watershed_results_sdr.shp 2022-05-19 14:13:21,013 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 watershed_results_sdr.shp 2022-05-19 14:13:21,054 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 1 2022-05-19 14:13:21,054 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 2.4895066e+05 1.5000000e+01 0.0000000e+00 9.9418970e+06 0.0000000e+00 -1.5000000e+01] for ('/Users/joaosantos/Documents/sdr_workspace/sed_retention.tif', 1) 2022-05-19 14:13:21,054 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 1 2022-05-19 14:13:21,054 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_sdr.shp 2022-05-19 14:13:21,159 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-05-19 14:13:21,159 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-05-19 14:13:21,160 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-05-19 14:13:21,160 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-05-19 14:13:21,161 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-05-19 14:13:21,161 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-05-19 14:13:21,162 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on watershed_results_sdr.shp 2022-05-19 14:13:21,162 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-05-19 14:13:21,162 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr.shp 2022-05-19 14:13:21,162 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on watershed_results_sdr.shp 2022-05-19 14:13:21,166 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-05-19 14:13:21,167 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on watershed_results_sdr.shp 2022-05-19 14:13:21,169 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 watershed_results_sdr.shp 2022-05-19 14:13:21,197 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 watershed_results_sdr.shp 2022-05-19 14:13:21,236 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 1 2022-05-19 14:13:21,236 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 2.4895066e+05 1.5000000e+01 0.0000000e+00 9.9418970e+06 0.0000000e+00 -1.5000000e+01] for ('/Users/joaosantos/Documents/sdr_workspace/sed_deposition.tif', 1) 2022-05-19 14:13:21,237 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 1 2022-05-19 14:13:21,237 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for watershed_results_sdr.shp 2022-05-19 14:13:21,242 (natcap.invest.utils) utils.prepare_workspace(170) INFO Elapsed time: 9.17s 2022-05-19 14:13:21,242 (natcap.invest.utils) utils.prepare_workspace(173) INFO Execution finished