DelineateIT error: "The datasets' intersection is empty (i.e., not all the datasets touch each other)"

NatCap Software Support
1

Hello,

I received my first error trying to use the DelineateIT tool with a DEM (1m res. - 29439x31843 cells).

I already checked the coordinate system of both files (DEM and outlet point) and they are both correct.

Below you can find the error log I received.

Thank you for your help.
Vinc.

2020-02-20 12:21:05,183 utils.prepare_workspace(109) INFO Writing log messages to H:\01_GISDATA\ES\DelineateIT\StudyArea_1\InVEST-DelineateIT-Watershed-Delineation-log-2020-02-20–12_21_05.txt
2020-02-20 12:21:05,186 model._logged_target(1633) Level 100 Starting model with parameters:
Arguments for InVEST natcap.invest.routing.delineateit 3.7.0:
dem_uri H:/01_GISDATA/ES/DelineateIT/StudyArea_1/DEM_1M_2018/DEM1M_2018.ovr
flow_threshold 4000
n_workers -1
outlet_shapefile_uri H:/01_GISDATA/ES/DelineateIT/StudyArea_1/Outlet_points/p2.shp
results_suffix ST1
snap_distance 35
workspace_dir H:\01_GISDATA\ES\DelineateIT\StudyArea_1

2020-02-20 12:21:05,188 delineateit.execute(60) INFO creating directory H:\01_GISDATA\ES\DelineateIT\StudyArea_1
2020-02-20 12:21:05,193 model._logged_target(1639) ERROR Exception while executing <function execute at 0x0F1A7230>
Traceback (most recent call last):
File “Z:\opt\atlassian\pipelines\agent\build\env\lib\site-packages\natcap\invest\ui\model.py”, line 1636, in _logged_target
File “Z:\opt\atlassian\pipelines\agent\build\env\lib\site-packages\natcap\invest\routing\delineateit.py”, line 79, in execute
File “Z:\opt\atlassian\pipelines\agent\build\env\lib\site-packages\natcap\invest\pygeoprocessing_0_3_3\routing\routing.py”, line 401, in delineate_watershed
File “Z:\opt\atlassian\pipelines\agent\build\env\lib\site-packages\natcap\invest\pygeoprocessing_0_3_3\geoprocessing.py”, line 3238, in tile_dataset_uri
File “Z:\opt\atlassian\pipelines\agent\build\env\lib\site-packages\natcap\invest\pygeoprocessing_0_3_3\geoprocessing.py”, line 2283, in vectorize_datasets
File “Z:\opt\atlassian\pipelines\agent\build\env\lib\site-packages\natcap\invest\pygeoprocessing_0_3_3\geoprocessing.py”, line 2035, in align_dataset_list
Exception: The datasets’ intersection is empty (i.e., not all the datasets touch each other).
2020-02-20 12:21:05,194 model._logged_target(1642) INFO Execution finished
2020-02-20 12:21:05,198 utils.prepare_workspace(115) INFO Elapsed time: 0.01s

2

Hi @Vincenzo -

Are the DEM raster and outlet shapefile in the same projected coordinate system? If not, that could be a cause of the error. Also, it looks like you’re entering an .ovr file for the DEM, try using the .tif since I’m not sure how InVEST handles .ovr.

~ Stacie

3

The .tif format works better.

Thank you

4

You might also like the new InVEST version 3.8.0, which has some nice performance updates for DelineateIT. And which won’t even let us make the mistake of using the .ovr !

5

Thank you Dave,

The 3.7 it allowed me to upload the pyramid file (.ovr)… The 3.8 does not.