2023-05-05 15:51:41,682 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [/Users/kevinli/Documents/PdS/SedimentDeliveryRatio/InVEST-natcap.invest.sdr.sdr-log-2023-05-05--15_51_41.txt] 2023-05-05 15:51:41,682 (__main__) cli.main(458) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.sdr.sdr 3.13.0: biophysical_table_path /Users/kevinli/Documents/Data/PdS_InVEST/BiophysicalTables/biophys_SDR.csv dem_path /Users/kevinli/Documents/Data/PdS_InVEST/DEM.tif drainage_path erodibility_path /Users/kevinli/Documents/Data/PdS_InVEST/erod_ca.tif erosivity_path /Users/kevinli/Documents/Data/PdS_InVEST/erosivity.tif ic_0_param 0.5 k_param 2 l_max 122 lulc_path /Users/kevinli/Documents/Data/PdS_InVEST/CDL_reproj/CDL_2017_reproj.tif n_workers -1 results_suffix sdr_PdS_2017 sdr_max 0.8 threshold_flow_accumulation 1000 watersheds_path /Users/kevinli/Documents/Data/Watersheds/site_huc10.shp workspace_dir /Users/kevinli/Documents/PdS/SedimentDeliveryRatio 2023-05-05 15:51:41,990 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-05-05 15:51:41,990 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-05-05 15:51:42,018 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-05-05 15:51:42,019 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-05-05 15:51:42,073 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 4 aligned: aligned_dem_sdr_PdS_2017.tif 2023-05-05 15:51:42,220 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-05-05 15:51:42,220 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-05-05 15:51:42,239 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-05-05 15:51:42,239 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-05-05 15:51:42,278 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 2 of 4 aligned: aligned_lulc_sdr_PdS_2017.tif 2023-05-05 15:51:42,423 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-05-05 15:51:42,423 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-05-05 15:51:42,439 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-05-05 15:51:42,439 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-05-05 15:51:42,480 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 3 of 4 aligned: aligned_erosivity_sdr_PdS_2017.tif 2023-05-05 15:51:42,610 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-05-05 15:51:42,610 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-05-05 15:51:42,624 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-05-05 15:51:42,624 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-05-05 15:51:42,664 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 4 of 4 aligned: aligned_erodibility_sdr_PdS_2017.tif 2023-05-05 15:51:42,665 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 4 rasters. 2023-05-05 15:51:43,126 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete 2023-05-05 15:51:43,336 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-05-05 15:51:43,337 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-05-05 15:51:43,373 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-05-05 15:51:43,373 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-05-05 15:51:44,029 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-05-05 15:51:44,032 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Calculating average aspect 2023-05-05 15:51:44,069 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 0.00% complete 2023-05-05 15:51:44,072 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 3.89% complete 2023-05-05 15:51:44,075 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 7.79% complete 2023-05-05 15:51:44,078 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 11.68% complete 2023-05-05 15:51:44,083 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 15.58% complete 2023-05-05 15:51:44,086 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 19.47% complete 2023-05-05 15:51:44,090 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 23.37% complete 2023-05-05 15:51:44,091 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 23.70% complete 2023-05-05 15:51:44,096 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 27.60% complete 2023-05-05 15:51:44,101 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 31.49% complete 2023-05-05 15:51:44,107 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 35.39% complete 2023-05-05 15:51:44,113 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 39.28% complete 2023-05-05 15:51:44,119 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 43.18% complete 2023-05-05 15:51:44,125 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 47.07% complete 2023-05-05 15:51:44,125 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 47.41% complete 2023-05-05 15:51:44,130 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 51.30% complete 2023-05-05 15:51:44,137 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 55.20% complete 2023-05-05 15:51:44,143 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 59.09% complete 2023-05-05 15:51:44,150 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 62.99% complete 2023-05-05 15:51:44,156 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 66.88% complete 2023-05-05 15:51:44,160 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 70.78% complete 2023-05-05 15:51:44,160 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 71.11% complete 2023-05-05 15:51:44,166 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 75.01% complete 2023-05-05 15:51:44,172 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 78.90% complete 2023-05-05 15:51:44,179 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 82.80% complete 2023-05-05 15:51:44,184 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 86.69% complete 2023-05-05 15:51:44,189 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 90.59% complete 2023-05-05 15:51:44,192 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 94.48% complete 2023-05-05 15:51:44,192 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 94.81% complete 2023-05-05 15:51:44,193 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 95.67% complete 2023-05-05 15:51:44,195 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 96.52% complete 2023-05-05 15:51:44,196 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 97.37% complete 2023-05-05 15:51:44,196 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 98.22% complete 2023-05-05 15:51:44,197 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 99.07% complete 2023-05-05 15:51:44,198 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 99.93% complete 2023-05-05 15:51:44,199 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Average aspect 100.00% complete 2023-05-05 15:51:44,914 (pygeoprocessing.routing.routing) Task._call(1093) INFO 100.0% complete 2023-05-05 15:51:44,931 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-05-05 15:51:44,931 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-05-05 15:51:44,934 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(571) ERROR exception encountered in raster_calculator Traceback (most recent call last): File "pygeoprocessing/geoprocessing.py", line 517, in raster_calculator File "natcap/invest/sdr/sdr.py", line 1116, in ls_factor_function IndexError: index 4 is out of bounds for axis 0 with size 4 2023-05-05 15:51:44,935 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(586) INFO Waiting for raster stats worker result. 2023-05-05 15:51:44,935 (pygeoprocessing.geoprocessing_core) threading.run(953) WARNING No valid pixels were received, sending None. 2023-05-05 15:51:44,936 (taskgraph.Task) Task.add_task(706) ERROR Something went wrong when adding task ls factor calculation (7), terminating taskgraph. Traceback (most recent call last): File "taskgraph/Task.py", line 674, in add_task File "taskgraph/Task.py", line 1093, in _call File "natcap/invest/sdr/sdr.py", line 1131, in _calculate_ls_factor File "pygeoprocessing/geoprocessing.py", line 517, in raster_calculator File "natcap/invest/sdr/sdr.py", line 1116, in ls_factor_function IndexError: index 4 is out of bounds for axis 0 with size 4 2023-05-05 15:51:44,936 (natcap.invest.utils) utils.prepare_workspace(166) ERROR Exception while executing natcap.invest.sdr.sdr Traceback (most recent call last): File "natcap/invest/utils.py", line 164, in prepare_workspace File "invest/cli.py", line 469, in main File "natcap/invest/sdr/sdr.py", line 634, in execute File "taskgraph/Task.py", line 674, in add_task File "taskgraph/Task.py", line 1093, in _call File "natcap/invest/sdr/sdr.py", line 1131, in _calculate_ls_factor File "pygeoprocessing/geoprocessing.py", line 517, in raster_calculator File "natcap/invest/sdr/sdr.py", line 1116, in ls_factor_function IndexError: index 4 is out of bounds for axis 0 with size 4 2023-05-05 15:51:44,936 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 3.25s 2023-05-05 15:51:44,936 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished