2025-02-25 14:30:19,331 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [E:\1ESV\SDR\SDR\InVEST-natcap.invest.sdr.sdr-log-2025-02-25--14_30_19.txt]
2025-02-25 14:30:19,331 (main) cli.main(459) Level 100 Starting model with parameters:
Arguments for InVEST natcap.invest.sdr.sdr 3.14.1:
biophysical_table_path E:\1ESV\SDR\SDR1.csv
dem_path E:\1ESV\DEM\Fill_DEM1.tif
drainage_path
erodibility_path E:\1ESV\SDR\PAWC和K\1\K1.tif
erosivity_path E:\1ESV\SDR\R\REI_2000.TIF
ic_0_param 0.5
k_param 2
l_max 122
lulc_path E:\1ESV\LUCC\【立方数据学社】CLCD_v01_2000_albert.tif
n_workers -1
results_suffix SDR2000
sdr_max 0.8
threshold_flow_accumulation 1000
watersheds_path E:\1ESV\bianjie\liuyu.shp
workspace_dir E:\1ESV\SDR\SDR
2025-02-25 14:30:22,672 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:30:22,676 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-1 (stats_worker), started daemon 14004)>
2025-02-25 14:30:23,404 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:30:23,404 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:30:24,318 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 1 of 4 aligned: aligned_dem_SDR2000.tif
2025-02-25 14:30:28,299 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:30:28,300 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-2 (stats_worker), started daemon 13376)>
2025-02-25 14:30:28,556 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:30:28,556 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:30:28,756 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 2 of 4 aligned: aligned_lulc_SDR2000.tif
2025-02-25 14:30:32,246 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:30:32,262 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-3 (stats_worker), started daemon 9364)>
2025-02-25 14:30:33,094 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:30:33,094 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:30:34,490 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 3 of 4 aligned: aligned_erosivity_SDR2000.tif
2025-02-25 14:30:37,598 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:30:37,604 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-4 (stats_worker), started daemon 15732)>
2025-02-25 14:30:38,313 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:30:38,313 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:30:39,360 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 4 of 4 aligned: aligned_erodibility_SDR2000.tif
2025-02-25 14:30:39,360 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1154) INFO aligned all 4 rasters.
2025-02-25 14:30:39,460 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:30:39,463 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-5 (stats_worker), started daemon 1888)>
2025-02-25 14:30:42,651 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:30:42,651 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:30:43,549 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:30:43,558 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-6 (stats_worker), started daemon 11740)>
2025-02-25 14:30:45,676 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:30:45,676 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:30:47,006 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:30:47,007 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-7 (stats_worker), started daemon 20172)>
2025-02-25 14:30:48,286 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:30:48,286 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:30:49,433 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:30:49,441 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-8 (stats_worker), started daemon 4680)>
2025-02-25 14:30:51,942 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:30:51,942 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:30:53,622 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:30:53,629 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-9 (stats_worker), started daemon 18804)>
2025-02-25 14:30:55,732 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:30:55,732 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:31:08,040 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): 42078976 of 59390406 pixels complete
2025-02-25 14:31:09,760 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete
2025-02-25 14:31:12,869 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:31:12,872 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-10 (stats_worker), started daemon 14112)>
2025-02-25 14:31:13,899 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:31:13,899 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:31:26,042 (pygeoprocessing.routing.routing) Task._call(1093) INFO 37.0% complete
2025-02-25 14:31:37,136 (pygeoprocessing.routing.routing) Task._call(1093) INFO 70.9% complete
2025-02-25 14:31:41,043 (pygeoprocessing.routing.routing) Task._call(1093) INFO Flow dir MFD 100.0% complete
2025-02-25 14:31:52,011 (pygeoprocessing.routing.routing) Task._call(1093) INFO Flow accum MFD 30.4% complete
2025-02-25 14:31:58,378 (pygeoprocessing.routing.routing) Task._call(1093) INFO Flow accum MFD 100.0% complete
2025-02-25 14:31:58,445 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:31:58,445 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-11 (stats_worker), started daemon 9888)>
2025-02-25 14:32:03,448 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(575) INFO ls_SDR2000.tif 76.2% complete
2025-02-25 14:32:04,844 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:32:04,844 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:32:14,808 (pygeoprocessing.routing.routing) Task._call(1093) INFO Extract streams MFD: filter out incomplete divergent streams
2025-02-25 14:32:15,032 (pygeoprocessing.routing.routing) Task._call(1093) INFO Extract streams MFD: 100.0% complete
2025-02-25 14:32:15,068 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:32:15,078 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-12 (stats_worker), started daemon 3036)>
2025-02-25 14:32:16,778 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:32:16,778 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:32:17,492 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:32:17,501 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-13 (stats_worker), started daemon 19532)>
2025-02-25 14:32:18,324 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:32:18,324 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:32:19,005 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:32:19,013 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-14 (stats_worker), started daemon 14320)>
2025-02-25 14:32:20,751 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:32:20,751 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:32:21,499 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:32:21,512 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-15 (stats_worker), started daemon 17996)>
2025-02-25 14:32:26,045 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:32:26,045 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:32:26,785 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:32:26,786 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-16 (stats_worker), started daemon 12496)>
2025-02-25 14:32:28,330 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:32:28,330 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:32:40,013 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1361) INFO Flow accum MFD 24.9% complete
2025-02-25 14:32:50,303 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1361) INFO Flow accum MFD 100.0% complete
2025-02-25 14:32:50,385 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:32:50,385 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-17 (stats_worker), started daemon 8320)>
2025-02-25 14:32:54,973 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:32:54,973 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:33:07,008 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1361) INFO Flow accum MFD 25.3% complete
2025-02-25 14:33:17,058 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1361) INFO Flow accum MFD 100.0% complete
2025-02-25 14:33:17,130 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:33:17,134 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-18 (stats_worker), started daemon 20072)>
2025-02-25 14:33:21,900 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:33:21,900 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:33:23,529 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:33:23,532 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-19 (stats_worker), started daemon 17424)>
2025-02-25 14:33:28,532 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(575) INFO d_up_SDR2000.tif 85.0% complete
2025-02-25 14:33:29,529 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:33:29,529 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:33:30,610 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:33:30,625 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-20 (stats_worker), started daemon 19688)>
2025-02-25 14:33:32,321 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:33:32,321 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:33:44,024 (pygeoprocessing.routing.routing) Task._call(1093) INFO Dist to channel MFD 61.6% complete
2025-02-25 14:33:48,893 (pygeoprocessing.routing.routing) Task._call(1093) INFO Dist to channel MFD 100.0% complete
2025-02-25 14:33:48,976 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:33:48,986 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-21 (stats_worker), started daemon 19872)>
2025-02-25 14:33:52,583 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:33:52,583 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:33:53,547 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:33:53,551 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-22 (stats_worker), started daemon 1764)>
2025-02-25 14:33:54,744 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:33:54,744 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:33:56,052 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:33:56,053 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-23 (stats_worker), started daemon 18960)>
2025-02-25 14:33:57,885 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:33:57,885 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:33:59,381 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:33:59,387 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-24 (stats_worker), started daemon 17204)>
2025-02-25 14:34:01,176 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:34:01,176 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:34:02,589 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Calculate sediment deposition
2025-02-25 14:34:08,075 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 23.66% complete
2025-02-25 14:34:14,025 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 34.66% complete
2025-02-25 14:34:20,008 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 43.79% complete
2025-02-25 14:34:26,109 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 50.50% complete
2025-02-25 14:34:32,093 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 59.41% complete
2025-02-25 14:34:38,110 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 66.67% complete
2025-02-25 14:34:44,093 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 72.50% complete
2025-02-25 14:34:50,841 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 84.82% complete
2025-02-25 14:34:53,085 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 100% complete
2025-02-25 14:34:53,600 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:34:53,608 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-25 (stats_worker), started daemon 20124)>
2025-02-25 14:34:56,309 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:34:56,309 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:34:58,437 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:34:58,440 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-26 (stats_worker), started daemon 19788)>
2025-02-25 14:35:03,423 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:35:03,423 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:35:05,683 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker
2025-02-25 14:35:05,691 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker <Thread(Thread-27 (stats_worker), started daemon 13676)>
2025-02-25 14:35:08,742 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete
2025-02-25 14:35:08,742 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result.
2025-02-25 14:35:09,024 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] One or several characters couldn't be converted correctly from UTF-8 to ISO-8859-1. This warning will not be emitted anymore.
