02/19/2021 18:17:02 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:02 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:02 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:02 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:02 geoprocessing.align_and_resize_raster_stack(793) INFO 1 of 4 aligned: aligned_dem.tif 02/19/2021 18:17:02 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:02 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:02 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:02 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:02 geoprocessing.align_and_resize_raster_stack(793) INFO 2 of 4 aligned: aligned_lulc.tif 02/19/2021 18:17:03 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:03 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:03 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:03 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:03 geoprocessing.align_and_resize_raster_stack(793) INFO 3 of 4 aligned: aligned_erosivity.tif 02/19/2021 18:17:03 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:03 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:03 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:03 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:03 geoprocessing.align_and_resize_raster_stack(793) INFO 4 of 4 aligned: aligned_erodibility.tif 02/19/2021 18:17:03 geoprocessing.align_and_resize_raster_stack(797) INFO aligned all 4 rasters. 02/19/2021 18:17:03 Task._call(1234) INFO 100.0% complete 02/19/2021 18:17:04 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:04 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:04 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:04 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:05 Task._call(1234) INFO 100.0% complete 02/19/2021 18:17:05 Task._call(1234) INFO Calculating average aspect 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 0.00% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 5.61% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 11.21% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 16.82% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 22.42% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 28.03% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 33.63% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 39.24% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 42.45% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 48.06% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 53.66% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 59.27% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 64.87% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 70.48% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 76.09% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 81.69% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 84.91% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 86.90% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 88.89% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 90.89% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 92.88% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 94.87% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 96.86% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 98.86% complete 02/19/2021 18:17:05 Task._call(1234) INFO Average aspect 100.00% complete 02/19/2021 18:17:06 Task._call(1234) INFO 100.0% complete 02/19/2021 18:17:06 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:06 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:06 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:06 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:06 Task._call(1234) INFO filter out incomplete divergent streams 02/19/2021 18:17:07 Task._call(1234) INFO 100.0% complete 02/19/2021 18:17:07 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:07 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:07 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:07 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:07 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:07 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:07 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:07 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:07 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:07 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:07 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:07 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:07 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:07 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:07 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:07 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:08 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:08 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:08 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:08 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:09 sdr._calculate_bar_factor(1091) INFO 100.0% complete 02/19/2021 18:17:09 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:09 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:09 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:09 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:10 sdr._calculate_bar_factor(1091) INFO 100.0% complete 02/19/2021 18:17:10 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:10 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:10 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:10 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:10 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:10 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:10 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:10 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:10 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:10 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:10 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:10 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:11 Task._call(1234) INFO 100.0% complete 02/19/2021 18:17:11 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:11 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:11 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:11 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:12 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:12 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:12 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:12 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:12 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:12 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:12 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:12 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:12 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:12 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:12 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:12 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:12 Task._call(1234) INFO Calculate sediment deposition 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 0.00% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 0.00% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 0.00% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 0.00% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 0.00% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 0.00% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 0.00% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 0.00% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 0.00% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 5.61% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 11.21% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 16.82% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 22.42% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 28.03% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 33.63% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 39.24% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 0.00% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 11.21% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 22.42% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 33.63% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 44.84% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 56.05% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 67.26% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 78.47% complete 02/19/2021 18:17:13 Task._call(1234) INFO Sediment deposition 100% complete 02/19/2021 18:17:14 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:14 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:14 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:14 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:14 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:14 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:14 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:14 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:15 Task._call(1234) INFO 100.0% complete 02/19/2021 18:17:15 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:15 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:15 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:15 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:15 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:15 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:15 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:15 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:15 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:15 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:15 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:15 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:16 geoprocessing.raster_calculator(391) INFO starting stats_worker 02/19/2021 18:17:16 geoprocessing.raster_calculator(397) INFO started stats_worker 02/19/2021 18:17:16 geoprocessing.raster_calculator(477) INFO 100.0% complete 02/19/2021 18:17:16 geoprocessing.raster_calculator(480) INFO Waiting for raster stats worker result. 02/19/2021 18:17:16 geoprocessing.align_and_resize_raster_stack(793) INFO 1 of 1 aligned: clipped_rasterdz5bmpae.tif 02/19/2021 18:17:16 geoprocessing.align_and_resize_raster_stack(797) INFO aligned all 1 rasters. 02/19/2021 18:17:16 geoprocessing.zonal_statistics(1244) INFO build a lookup of aggregate field value to FID 02/19/2021 18:17:16 geoprocessing.zonal_statistics(1251) INFO creating disjoint polygon set 02/19/2021 18:17:16 geoprocessing.calculate_disjoint_polygon_set(2183) INFO build shapely polygon list 02/19/2021 18:17:16 Task.add_task(779) ERROR Something went wrong when adding task generate report (30), terminating taskgraph. Traceback (most recent call last): File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\taskgraph\Task.py", line 747, in add_task new_task._call() File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\taskgraph\Task.py", line 1234, in _call payload = self._func(*self._args, **self._kwargs) File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\natcap\invest\sdr\sdr.py", line 1369, in _generate_report 'usle_tot': pygeoprocessing.zonal_statistics( File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\pygeoprocessing\geoprocessing.py", line 1252, in zonal_statistics disjoint_fid_sets = calculate_disjoint_polygon_set( File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\pygeoprocessing\geoprocessing.py", line 2187, in calculate_disjoint_polygon_set bounding_box = shapely.prepared.prep(shapely.geometry.box(*bounding_box)) File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\shapely\geometry\geo.py", line 59, in box return Polygon(coords) File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\shapely\geometry\polygon.py", line 243, in __init__ ret = geos_polygon_from_py(shell, holes) File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\shapely\geometry\polygon.py", line 509, in geos_polygon_from_py ret = geos_linearring_from_py(shell) File "shapely\speedups\_speedups.pyx", line 408, in shapely.speedups._speedups.geos_linearring_from_py ValueError: GEOSGeom_createLinearRing_r returned a NULL pointer Traceback (most recent call last): File "C:\Users\ligia\Desktop\testesdr3.py", line 39, in natcap.invest.sdr.sdr.execute(args) File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\natcap\invest\sdr\sdr.py", line 708, in execute _ = task_graph.add_task( File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\taskgraph\Task.py", line 747, in add_task new_task._call() File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\taskgraph\Task.py", line 1234, in _call payload = self._func(*self._args, **self._kwargs) File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\natcap\invest\sdr\sdr.py", line 1369, in _generate_report 'usle_tot': pygeoprocessing.zonal_statistics( File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\pygeoprocessing\geoprocessing.py", line 1252, in zonal_statistics disjoint_fid_sets = calculate_disjoint_polygon_set( File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\pygeoprocessing\geoprocessing.py", line 2187, in calculate_disjoint_polygon_set bounding_box = shapely.prepared.prep(shapely.geometry.box(*bounding_box)) File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\shapely\geometry\geo.py", line 59, in box return Polygon(coords) File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\shapely\geometry\polygon.py", line 243, in __init__ ret = geos_polygon_from_py(shell, holes) File "C:\Users\ligia\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.8_qbz5n2kfra8p0\LocalCache\local-packages\Python38\site-packages\shapely\geometry\polygon.py", line 509, in geos_polygon_from_py ret = geos_linearring_from_py(shell) File "shapely\speedups\_speedups.pyx", line 408, in shapely.speedups._speedups.geos_linearring_from_py ValueError: GEOSGeom_createLinearRing_r returned a NULL pointer