2024-04-17 21:47:03,670 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [E:\datasets\invset_OUTPUT0417\SDR_T4\InVEST-natcap.invest.sdr.sdr-log-2024-04-17--21_47_03.txt] 2024-04-17 21:47:03,670 (__main__) cli.main(459) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.sdr.sdr 3.14.1: biophysical_table_path E:\datasets\invest process\SOIL\BIOPHYSICTAB_ND.csv dem_path E:\datasets\predata\dem125\dem_sfq0417_Resample1.tif drainage_path erodibility_path E:\datasets\invest process\SOIL\K0406.tif erosivity_path E:\datasets\invest process\SOIL\RR04061.tif ic_0_param 0.5 k_param 2 l_max 122 lulc_path E:\datasets\LULC0415\LULC_2020NSFQ.TIF n_workers -1 results_suffix T4_2020 sdr_max 0.8 threshold_flow_accumulation 1000 watersheds_path E:\datasets\invest process\SOIL\SFQ_LYBJ0331.shp workspace_dir E:\datasets\invset_OUTPUT0417\SDR_T4 2024-04-17 21:47:04,048 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:04,049 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:04,101 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:04,101 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:04,164 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 1 of 4 aligned: aligned_dem_T4_2020.tif 2024-04-17 21:47:04,417 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:04,418 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:04,452 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:04,452 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:04,483 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 2 of 4 aligned: aligned_lulc_T4_2020.tif 2024-04-17 21:47:04,897 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:04,898 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:04,974 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:04,975 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:05,073 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 3 of 4 aligned: aligned_erosivity_T4_2020.tif 2024-04-17 21:47:05,450 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:05,450 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:05,524 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:05,524 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:05,604 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1150) INFO 4 of 4 aligned: aligned_erodibility_T4_2020.tif 2024-04-17 21:47:05,605 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1154) INFO aligned all 4 rasters. 2024-04-17 21:47:05,636 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:05,637 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:05,834 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:05,834 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:05,940 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:05,940 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:06,042 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:06,042 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:06,152 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:06,152 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:06,238 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:06,239 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:06,333 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:06,334 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:06,472 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:06,472 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:06,584 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:06,585 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:06,724 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:06,725 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:08,002 (pygeoprocessing.routing.routing) Task._call(1093) INFO (fill pits): complete 2024-04-17 21:47:08,287 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:08,287 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:08,408 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:08,408 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:11,762 (pygeoprocessing.routing.routing) Task._call(1093) INFO Flow dir MFD 100.0% complete 2024-04-17 21:47:13,322 (pygeoprocessing.routing.routing) Task._call(1093) INFO Flow accum MFD 100.0% complete 2024-04-17 21:47:13,345 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:13,345 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:13,907 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:13,907 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:14,797 (pygeoprocessing.routing.routing) Task._call(1093) INFO Extract streams MFD: filter out incomplete divergent streams 2024-04-17 21:47:14,826 (pygeoprocessing.routing.routing) Task._call(1093) INFO Extract streams MFD: 100.0% complete 2024-04-17 21:47:14,851 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:14,852 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:15,003 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:15,003 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:15,093 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:15,094 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:15,182 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:15,182 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:15,279 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:15,280 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:15,433 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:15,434 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:15,536 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:15,537 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:15,865 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:15,866 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:16,038 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:16,038 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:16,202 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:16,202 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:18,222 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1361) INFO Flow accum MFD 100.0% complete 2024-04-17 21:47:18,237 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:18,238 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:18,633 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:18,634 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:20,664 (pygeoprocessing.routing.routing) sdr._calculate_bar_factor(1361) INFO Flow accum MFD 100.0% complete 2024-04-17 21:47:20,681 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:20,681 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:21,101 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:21,101 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:21,415 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:21,416 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:21,914 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:21,914 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:22,226 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:22,226 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:22,372 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:22,373 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:24,052 (pygeoprocessing.routing.routing) Task._call(1093) INFO Dist to channel MFD 100.0% complete 2024-04-17 21:47:24,068 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:24,069 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:24,390 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:24,390 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:24,533 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:24,534 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:24,641 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:24,641 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:24,797 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:24,798 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:24,970 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:24,970 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:25,120 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:25,121 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:25,281 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:25,282 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:25,417 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Calculate sediment deposition 2024-04-17 21:47:28,318 (natcap.invest.sdr.sdr_core) Task._call(1093) INFO Sediment deposition 100% complete 2024-04-17 21:47:28,604 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:28,605 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:28,776 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:28,777 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:28,952 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:28,953 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:29,224 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:29,224 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:29,386 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(495) INFO starting stats_worker 2024-04-17 21:47:29,387 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(501) INFO started stats_worker 2024-04-17 21:47:29,572 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(580) INFO 100.0% complete 2024-04-17 21:47:29,573 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(583) INFO Waiting for raster stats worker result. 2024-04-17 21:47:29,603 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 267.278999999999996 of field AREA of feature 1 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,603 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 76.8829999999999956 of field AREA of feature 2 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,603 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 409.432000000000016 of field AREA of feature 3 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,604 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 256.163000000000011 of field AREA of feature 4 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,604 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 90.1049000000000007 of field AREA of feature 5 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,604 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 89.5742999999999938 of field AREA of feature 6 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,604 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 80.8633999999999986 of field AREA of feature 7 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,604 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 52.224499999999999 of field AREA of feature 8 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,604 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 151.818999999999988 of field AREA of feature 10 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,605 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 671.975999999999999 of field AREA of feature 11 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,605 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 116.909000000000006 of field AREA of feature 12 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,605 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 39.7006000000000014 of field AREA of feature 13 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,606 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 100.507999999999996 of field AREA of feature 17 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:29,674 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2024-04-17 21:47:29,734 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection 2024-04-17 21:47:29,736 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2024-04-17 21:47:29,738 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg 2024-04-17 21:47:29,767 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector 2024-04-17 21:47:29,944 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets 2024-04-17 21:47:29,945 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list 2024-04-17 21:47:29,946 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index 2024-04-17 21:47:29,948 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:29,948 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup 2024-04-17 21:47:29,952 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:29,953 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:30,182 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1 2024-04-17 21:47:30,185 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 5 2024-04-17 21:47:30,283 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 1 of 5 2024-04-17 21:47:30,358 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 2 of 5 2024-04-17 21:47:30,425 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 3 of 5 2024-04-17 21:47:30,484 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 4 of 5 2024-04-17 21:47:30,537 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:30,569 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2024-04-17 21:47:30,626 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection 2024-04-17 21:47:30,629 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2024-04-17 21:47:30,631 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg 2024-04-17 21:47:30,651 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector 2024-04-17 21:47:30,837 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets 2024-04-17 21:47:30,838 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list 2024-04-17 21:47:30,839 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index 2024-04-17 21:47:30,840 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:30,840 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup 2024-04-17 21:47:30,843 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:30,843 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:31,056 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1 2024-04-17 21:47:31,059 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 5 2024-04-17 21:47:31,145 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 1 of 5 2024-04-17 21:47:31,217 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 2 of 5 2024-04-17 21:47:31,285 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 3 of 5 2024-04-17 21:47:31,350 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 4 of 5 2024-04-17 21:47:31,410 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:31,442 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2024-04-17 21:47:31,494 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection 2024-04-17 21:47:31,496 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2024-04-17 21:47:31,498 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg 2024-04-17 21:47:31,520 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector 2024-04-17 21:47:31,693 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets 2024-04-17 21:47:31,694 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list 2024-04-17 21:47:31,695 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index 2024-04-17 21:47:31,696 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:31,696 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup 2024-04-17 21:47:31,699 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:31,699 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:31,909 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1 2024-04-17 21:47:31,911 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 5 2024-04-17 21:47:31,996 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 1 of 5 2024-04-17 21:47:32,067 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 2 of 5 2024-04-17 21:47:32,133 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 3 of 5 2024-04-17 21:47:32,194 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 4 of 5 2024-04-17 21:47:32,257 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:32,299 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2024-04-17 21:47:32,351 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection 2024-04-17 21:47:32,353 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2024-04-17 21:47:32,354 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg 2024-04-17 21:47:32,378 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector 2024-04-17 21:47:32,580 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets 2024-04-17 21:47:32,580 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list 2024-04-17 21:47:32,582 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index 2024-04-17 21:47:32,582 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:32,582 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup 2024-04-17 21:47:32,585 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:32,586 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:32,808 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1 2024-04-17 21:47:32,811 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 5 2024-04-17 21:47:32,902 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 1 of 5 2024-04-17 21:47:32,977 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 2 of 5 2024-04-17 21:47:33,046 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 3 of 5 2024-04-17 21:47:33,106 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 4 of 5 2024-04-17 21:47:33,166 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:33,204 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2024-04-17 21:47:33,256 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2250) INFO starting reprojection 2024-04-17 21:47:33,261 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2024-04-17 21:47:33,262 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2290) INFO reprojection 100.0% complete on reprojected.gpkg 2024-04-17 21:47:33,284 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1785) INFO Clipping rasters to their intersection with the vector 2024-04-17 21:47:33,484 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1810) INFO calculating disjoint polygon sets 2024-04-17 21:47:33,484 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2781) INFO build shapely polygon list 2024-04-17 21:47:33,485 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2807) INFO build shapely rtree index 2024-04-17 21:47:33,486 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2820) INFO poly feature lookup 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:33,486 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2824) INFO build poly intersection lookup 2024-04-17 21:47:33,489 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2850) INFO poly intersection feature lookup 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:33,490 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2882) INFO maximal subset build 100.0% complete on watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:33,704 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1846) INFO calculating stats on raster 0 of 1 2024-04-17 21:47:33,707 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 0 of 5 2024-04-17 21:47:33,796 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 1 of 5 2024-04-17 21:47:33,872 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 2 of 5 2024-04-17 21:47:33,946 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 3 of 5 2024-04-17 21:47:34,006 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1856) INFO disjoint polygon set 4 of 5 2024-04-17 21:47:34,065 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1984) INFO all done processing polygon sets for watershed_results_sdr_T4_2020.shp 2024-04-17 21:47:34,077 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 267.278999999999996 of field AREA of feature 1 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,077 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 76.8829999999999956 of field AREA of feature 2 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,077 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 409.432000000000016 of field AREA of feature 3 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,077 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 256.163000000000011 of field AREA of feature 4 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,077 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 90.1049000000000007 of field AREA of feature 5 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,078 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 89.5742999999999938 of field AREA of feature 6 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,078 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 80.8633999999999986 of field AREA of feature 7 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,078 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 52.224499999999999 of field AREA of feature 8 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,078 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 151.818999999999988 of field AREA of feature 10 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,078 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 671.975999999999999 of field AREA of feature 11 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,078 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 116.909000000000006 of field AREA of feature 12 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,078 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 39.7006000000000014 of field AREA of feature 13 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,079 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 100.507999999999996 of field AREA of feature 17 not successfully written. Possibly due to too larger number with respect to field width 2024-04-17 21:47:34,083 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 30.41s 2024-04-17 21:47:34,083 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished