TransWikia.com

Why am I getting this error when creating a DEM?

Geographic Information Systems Asked by DanG on February 13, 2021

I’m using QGIS and trying to create a DEM using this tutorial:
https://www.youtube.com/watch?v=aJm9lewzooU
I’m getting error messages when trying to run "natural neighbor" with my file of coordinates. I pasted the error log below, as well as the screenshot for the settings I used. Do you know what is causing this error? enter image description here

QGIS version: 3.16.1-Hannover
QGIS code revision: b381a90dca[![enter image description here][2]][2]
Qt version: 5.11.2
GDAL version: 3.1.4
GEOS version: 3.8.1-CAPI-1.13.3
PROJ version: Rel. 6.3.2, May 1st, 2020
Processing algorithm…
Algorithm 'Natural neighbour' starting…
Input parameters:
{ 'FIELD' : 'altitude (m)', 'METHOD' : 1, 'SHAPES' : 'delimitedtext://file:///C:/Users/EVG/Desktop/latlong.csv?type=csv&maxFields=10000&detectTypes=yes&xField=longitude&yField=latitude&crs=EPSG:26915&spatialIndex=no&subsetIndex=no&watchFile=no', 'TARGET_OUT_GRID' : 'TEMPORARY_OUTPUT', 'TARGET_TEMPLATE' : None, 'TARGET_USER_FITS' : 0, 'TARGET_USER_SIZE' : 0.005, 'TARGET_USER_XMIN TARGET_USER_XMAX TARGET_USER_YMIN TARGET_USER_YMAX' : '-93.137700000,-93.137400000,45.008800000,45.009000000 [EPSG:26915]', 'WEIGHT' : 0 }

grid_gridding "Natural Neighbour" -TARGET_DEFINITION 0 -SHAPES "C:/Users/EVG/AppData/Local/Temp/processing_cqmvma/d04b7f8a78374cd7b19d97efc90135ce/SHAPES.shp" -FIELD "altitude (m)" -METHOD 1 -TARGET_USER_XMIN -93.1377 -TARGET_USER_XMAX -93.1374 -TARGET_USER_YMIN 45.0088 -TARGET_USER_YMAX 45.009 -TARGET_USER_SIZE 0.005 -TARGET_USER_FITS 0 -WEIGHT 0.0 -TARGET_OUT_GRID "C:/Users/EVG/AppData/Local/Temp/processing_cqmvma/3c02fc10349c4e929b80d64fd7598377/TARGET_OUT_GRID.sdat"

C:UsersEVGDocuments>set SAGA=C:/PROGRA~1/QGIS3~1.16/appssaga-ltr

C:UsersEVGDocuments>set SAGA_MLB=C:/PROGRA~1/QGIS3~1.16/appssaga-ltrmodules

C:UsersEVGDocuments>PATH=C:PROGRA~1QGIS3~1.16appsqgisbin;C:PROGRA~1QGIS3~1.16appsPython37;C:PROGRA~1QGIS3~1.16appsPython37Scripts;C:PROGRA~1QGIS3~1.16appsqt5bin;C:PROGRA~1QGIS3~1.16bin;C:Windowssystem32;C:Windows;C:Windowssystem32WBem;C:PROGRA~1QGIS3~1.16appsPython37libsite-packagespywin32_system32;C:PROGRA~1QGIS3~1.16appsPython37libsite-packagesnumpy.libs;C:/PROGRA~1/QGIS3~1.16/appssaga-ltr;C:/PROGRA~1/QGIS3~1.16/appssaga-ltrmodules

C:UsersEVGDocuments>call saga_cmd grid_gridding "Natural Neighbour" -TARGET_DEFINITION 0 -SHAPES "C:/Users/EVG/AppData/Local/Temp/processing_cqmvma/d04b7f8a78374cd7b19d97efc90135ce/SHAPES.shp" -FIELD "altitude (m)" -METHOD 1 -TARGET_USER_XMIN -93.1377 -TARGET_USER_XMAX -93.1374 -TARGET_USER_YMIN 45.0088 -TARGET_USER_YMAX 45.009 -TARGET_USER_SIZE 0.005 -TARGET_USER_FITS 0 -WEIGHT 0.0 -TARGET_OUT_GRID "C:/Users/EVG/AppData/Local/Temp/processing_cqmvma/3c02fc10349c4e929b80d64fd7598377/TARGET_OUT_GRID.sdat"
____________________________

##### ## ##### ##
### ### ## ###
### # ## ## #### # ##
### ##### ## # #####
##### # ## ##### # ##
____________________________

SAGA Version: 2.3.2 (64 bit)

____________________________
library path: C:PROGRA~1QGIS3~1.16appssaga-ltrmodules
library name: grid_gridding
library : Gridding
tool : Natural Neighbour
author : O.Conrad (c) 2008
processors : 8 [8]
____________________________

Load shapes: C:/Users/EVG/AppData/Local/Temp/processing_cqmvma/d04b7f8a78374cd7b19d97efc90135ce/SHAPES.shp...




Parameters


Points: SHAPES
Attribute: type
Target Grid System: user defined
Cellsize: 0.005000
Left: -93.137400
Right: -93.137400
Bottom: 45.009000
Top: 45.009000
Fit: nodes
Method: Sibson
Minimum Weight: 0.000000



Save grid: C:/Users/EVG/AppData/Local/Temp/processing_cqmvma/3c02fc10349c4e929b80d64fd7598377/TARGET_OUT_GRID.sdat...



C:UsersEVGDocuments>exit
Execution completed in 6.89 seconds
Results:
{'TARGET_OUT_GRID': 'C:/Users/EVG/AppData/Local/Temp/processing_cqmvma/3c02fc10349c4e929b80d64fd7598377/TARGET_OUT_GRID.sdat'}

Loading resulting layers
The following layers were not correctly generated.
• C:/Users/EVG/AppData/Local/Temp/processing_cqmvma/296776b11bdb4764bdd1bbc3ab300d77/TARGET_OUT_GRID.sdat
You can check the 'Log Messages Panel' in QGIS main window to find more information about the execution of the algorithm.

One Answer

The output extend you defined is quite weird:

'TARGET_USER_XMIN TARGET_USER_XMAX TARGET_USER_YMIN TARGET_USER_YMAX' : '-93.137700000,-93.137400000,45.008800000,45.009000000 [EPSG:26915]'

You use EPSG:26915 (NAD83 / UTM zone 15N), but your coordinates seem to be in EPSG:4326. The coordinates -93.137700000,45.008800000 in EPSG:26915 are located somewhere in the Pacific, several hundred km west of the Galapagos islands. If I'm not mistaken, you are interested in an area around Minneapolis, however. This corresponds to your coordinates if they are interpreted in EPSG:4326. -93.137700000,45.008800000, this place here - this is the southwestern corner (x_min, y_min) of your extent, right? However, if I have a look at the north eastern corner of your extent (interpreting the coordinate values in EPSG:4326), it is just about 20 meters away. So something is wrong anyway with these coordinates.

The coordinates -93.137700000,45.008800000 in EPSG:26915 correspond to -0.0008 degrees N, -97.4883 degrees W in EPSG:4326, this point here out in the ocean. You have defined this as the corner in the SW of your extent (and similarly the other coordinates).

Do never interpret coordinates from one CRS with another one: they have to be re-calculated (transformed) between different CRS. You should re-define your extend in the CRS your layer is. See details here to understand why your coordinates are somewhere out in the Pacific: https://gis.stackexchange.com/a/383437/88814

Answered by Babel on February 13, 2021

Add your own answers!

Ask a Question

Get help from others!

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP