QGIS version: 3.34.2-Prizren QGIS code revision: 7d199797fc Qt version: 5.15.3 Python version: 3.9.5 GDAL version: 3.8.2 GEOS version: 3.12.1-CAPI-1.18.1 PROJ version: Rel. 9.3.1, December 1st, 2023 PDAL version: 2.6.0 (git-version: 3fced5) Algorithm started at: 2023-12-27T23:38:15 Algorithm 'Stream and Catchment Delineation' starting… Input parameters: { 'INPUT_CRS' : QgsCoordinateReferenceSystem('EPSG:3857'), 'INPUT_DEM' : 1, 'INPUT_EXTENT' : '1416721.999200000,1438659.174700000,5190004.766000000,5200973.353800000 [EPSG:3857]', 'INPUT_KEY' : '9cfea4184505b3ed4186*****', 'INPUT_OUTLET' : '1427598.155018,5195571.221617 [EPSG:3857]', 'INPUT_THRESHOLD' : 6, 'INPUT_TOLERANCE' : 250, 'OUTPUT_CATCHMENT' : 'ogr:dbname=\'C:/Users/viking/Documents/pcrastertry/5t/vik.gpkg\' table="Catchment" (geom)', 'OUTPUT_DEM' : 'C:/Users/viking/Documents/pcrastertry/5t/dem.map', 'OUTPUT_FLOWDIRECTION' : 'C:/Users/viking/Documents/pcrastertry/5t/flowdirection.map', 'OUTPUT_FOLDER' : 'C:\\Users\\viking\\Documents\\pcrastertry\\5t', 'OUTPUT_STREAMS' : 'ogr:dbname=\'C:/Users/viking/Documents/pcrastertry/5t/vik.gpkg\' table="Streams" (geom)' } Coordinates written to C:\Users\viking\Documents\pcrastertry\5t\coordinates.csv Downloading DEM... 8 KB downloaded 120 KB downloaded Successfully downloaded https://portal.opentopography.org/API/globaldem?demtype=SRTMGL1&south=42.19347168356631&north=42.2664303914047&west=12.726630252297245&east=12.923695252717865&outputFormat=GTiff&API_Key=9cfea4184505b3ed4*** Results: {'OUTPUT': 'C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/985f646564b84774a0c36824f2dc8709/OUTPUT.nc'} Results: {'OUTPUT': 'OUTPUT_e6d902c7_e8d6_4d18_8929_e96f20ede107'} Results: {'OUTPUT': 'C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/985f646564b84774a0c36824f2dc8709/OUTPUT.nc'} Reprojecting DEM... GDAL command: gdalwarp -overwrite -t_srs EPSG:3857 -dstnodata -9999.0 -tr 30.0 30.0 -r near -of netCDF C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/985f646564b84774a0c36824f2dc8709/OUTPUT.nc C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/3074e234236e407695c6d0c679874e1d/OUTPUT.nc GDAL command output: Creating output file that is 731P x 366L. Processing C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/985f646564b84774a0c36824f2dc8709/OUTPUT.nc [1/1] : 0Using internal nodata values (e.g. -32768) for image C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/985f646564b84774a0c36824f2dc8709/OUTPUT.nc. ...10...20...30...40...50...60...70...80...90...100 - done. Process completed successfully Results: {'OUTPUT': 'C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/3074e234236e407695c6d0c679874e1d/OUTPUT.nc'} Converting DEM to PCRaster format... Results: {'OUTPUT': 'C:/Users/viking/Documents/pcrastertry/5t/dem.map'} Calculating flow direction (this can take some time!)... Calculating Strahler orders... Deriving stream network raster... Snapping outlet to stream... Delineating catchment raster... Thinning the stream network raster... g.proj -c wkt="C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/5d9e04915ed34ade8ba2c43a14acc170/crs.prj" r.in.gdal input="C:\Users\viking\Documents\pcrastertry\5t\river.map" band=1 output="rast_658c6820078642" --overwrite -o g.region n=5200988.00331236 s=5190008.00331236 e=1438640.5917915162 w=1416710.5917915162 res=30.0 r.thin input=rast_658c6820078642 iterations=200 output=output744a5fc4802041deab786908da6d66c4 --overwrite g.region raster=output744a5fc4802041deab786908da6d66c4 r.out.gdal -t -m input="output744a5fc4802041deab786908da6d66c4" output="C:\Users\viking\AppData\Local\Temp\processing_yNrCwH\f83232ca79324394ab4e400283dbfe1e\output.nc" format="netCDF" --overwrite WARNING: Concurrent mapset locking is not supported on Windows C:\Users\viking\Documents>chcp 1252 1>NUL C:\Users\viking\Documents>g.proj -c wkt="C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/5d9e04915ed34ade8ba2c43a14acc170/crs.prj" WARNING: WKT validation grammar errors: WARNING: whitespace only string ERROR: WKT not recognized: C:\Users\viking\Documents>r.in.gdal input="C:\Users\viking\Documents\pcrastertry\5t\river.map" band=1 output="rast_658c6820078642" --overwrite -o No projection information available Over-riding projection check Importing raster map ... 0..3..6..9..12..15..18..21..24..27..30..33..36..39..42..45..48..51..54..57..60..63..66..69..72..75..78..81..84..87..90..93..96..99..100 C:\Users\viking\Documents>g.region n=5200988.00331236 s=5190008.00331236 e=1438640.5917915162 w=1416710.5917915162 res=30.0 C:\Users\viking\Documents>r.thin input=rast_658c6820078642 iterations=200 output=output744a5fc4802041deab786908da6d66c4 --overwrite Raster map - 366 rows X 731 columns Bounding box: l = 5, r = 681, t = 3, b = 367 Pass number 1 Deleted 274 pixels Pass number 2 Deleted 0 pixels Thinning completed successfully. Output map 366 rows X 731 columns Window 366 rows X 731 columns C:\Users\viking\Documents>g.region raster=output744a5fc4802041deab786908da6d66c4 C:\Users\viking\Documents>r.out.gdal -t -m input="output744a5fc4802041deab786908da6d66c4" output="C:\Users\viking\AppData\Local\Temp\processing_yNrCwH\f83232ca79324394ab4e400283dbfe1e\output.nc" format="netCDF" --overwrite Checking GDAL data type and nodata value... 2..5..8..11..14..17..20..23..26..29..32..35..38..41..44..47..50..53..56..59..62..65..68..71..74..77..80..83..86..89..92..95..98..100 Using GDAL data type Input raster map contains cells with NULL-value (no-data). The value 255 will be used to represent no-data values in the input map. You can specify a nodata value with the nodata option. WARNING: Unable to set projection Exporting raster data to netCDF format... 2..5..8..11..14..17..20..23..26..29..32..35..38..41..44..47..50..53..56..59..62..65..68..71..74..77..80..83..86..89..92..95..98..100 r.out.gdal complete. File created. C:\Users\viking\Documents>exit Results: {'output': 'C:\\Users\\viking\\AppData\\Local\\Temp\\processing_yNrCwH\\f83232ca79324394ab4e400283dbfe1e\\output.nc'} Vectorizing the stream network... g.proj -c wkt="C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/e37aaae119a2460681afa068611a9892/crs.prj" r.in.gdal input="C:\Users\viking\AppData\Local\Temp\processing_yNrCwH\f83232ca79324394ab4e400283dbfe1e\output.nc" band=1 output="rast_658c68217d0a53" --overwrite -o g.region n=5200988.00331236 s=5190008.00331236 e=1438640.59179152 w=1416710.59179152 res=30.0 r.to.vect input=rast_658c68217d0a53 type="line" column="value" -v output=output4a4be05f907149bdbb247f12656a58b6 --overwrite v.out.ogr type="auto" input="output4a4be05f907149bdbb247f12656a58b6" output="C:\Users\viking\AppData\Local\Temp\processing_yNrCwH\0bc0ec727dd044e3a1b11c43def6a085\output.gpkg" format="GPKG" --overwrite WARNING: Concurrent mapset locking is not supported on Windows C:\Users\viking\Documents>chcp 1252 1>NUL C:\Users\viking\Documents>g.proj -c wkt="C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/e37aaae119a2460681afa068611a9892/crs.prj" WARNING: WKT validation grammar errors: WARNING: whitespace only string ERROR: WKT not recognized: C:\Users\viking\Documents>r.in.gdal input="C:\Users\viking\AppData\Local\Temp\processing_yNrCwH\f83232ca79324394ab4e400283dbfe1e\output.nc" band=1 output="rast_658c68217d0a53" --overwrite -o No projection information available Over-riding projection check Importing raster map ... 0..3..6..9..12..15..18..21..24..27..30..33..36..39..42..45..48..51..54..57..60..63..66..69..72..75..78..81..84..87..90..93..96..99..100 C:\Users\viking\Documents>g.region n=5200988.00331236 s=5190008.00331236 e=1438640.59179152 w=1416710.59179152 res=30.0 C:\Users\viking\Documents>r.to.vect input=rast_658c68217d0a53 type="line" column="value" -v output=output4a4be05f907149bdbb247f12656a58b6 --overwrite Extracting lines... 2..5..8..11..14..17..20..23..26..29..32..35..38..41..44..47..50..53..56..59..62..65..68..71..74..77..80..83..86..89..92..95..98..100 100 Updating attributes... 0..33..66..100 Building topology for vector map ... Registering primitives... r.to.vect complete. C:\Users\viking\Documents>v.out.ogr type="auto" input="output4a4be05f907149bdbb247f12656a58b6" output="C:\Users\viking\AppData\Local\Temp\processing_yNrCwH\0bc0ec727dd044e3a1b11c43def6a085\output.gpkg" format="GPKG" --overwrite Exporting 133 features... 5..11..17..23..29..35..41..47..53..59..65..71..77..83..89..95..100 v.out.ogr complete. 133 features (Line String type) written to (GPKG format). C:\Users\viking\Documents>exit Results: {'output': 'C:\\Users\\viking\\AppData\\Local\\Temp\\processing_yNrCwH\\0bc0ec727dd044e3a1b11c43def6a085\\output.gpkg'} Results: {'OUTPUT': 'Assigned_CRS_e30993f9_81bc_4999_95dc_44188462f84f'} makeValid failed for feature 14 makeValid failed for feature 16 makeValid failed for feature 20 makeValid failed for feature 22 makeValid failed for feature 24 makeValid failed for feature 30 makeValid failed for feature 32 makeValid failed for feature 36 makeValid failed for feature 41 makeValid failed for feature 43 makeValid failed for feature 45 makeValid failed for feature 48 makeValid failed for feature 51 makeValid failed for feature 57 makeValid failed for feature 61 makeValid failed for feature 63 makeValid failed for feature 66 makeValid failed for feature 77 makeValid failed for feature 85 makeValid failed for feature 88 makeValid failed for feature 95 makeValid failed for feature 110 makeValid failed for feature 116 makeValid failed for feature 122 Results: {'OUTPUT': 'C:/Users/viking/Documents/pcrastertry/5t/vik.gpkg|layername=Streams'} Polygonizing the catchment boundary... GDAL command: gdal_polygonize.bat C:\Users\viking\Documents\pcrastertry\5t\catchment.map -b 1 -f "Interlis 2" C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/322fca63ab5546be823900f460cdb91a/OUTPUT.xtf OUTPUT DN GDAL command output: Creating output C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/322fca63ab5546be823900f460cdb91a/OUTPUT.xtf of format Interlis 2. Traceback (most recent call last): File "C:\PROGRA~1\QGIS33~2.0\apps\Python39\lib\site-packages\osgeo_utils\auxiliary\gdal_argparse.py", line 226, in main self.doit(**kwargs) File "C:\PROGRA~1\QGIS33~2.0\apps\Python39\lib\site-packages\osgeo_utils\gdal_polygonize.py", line 351, in doit return gdal_polygonize(**kwargs) File "C:\PROGRA~1\QGIS33~2.0\apps\Python39\lib\site-packages\osgeo_utils\gdal_polygonize.py", line 162, in gdal_polygonize dst_layer = dst_ds.CreateLayer( AttributeError: 'NoneType' object has no attribute 'CreateLayer' Process returned error code 1 Results: {'OUTPUT': 'C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/322fca63ab5546be823900f460cdb91a/OUTPUT.xtf'} Unable to execute algorithm Could not load source layer for INPUT: C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/322fca63ab5546be823900f460cdb91a/OUTPUT.xtf not found Unable to execute algorithm Could not load source layer for INPUT: C:/Users/viking/AppData/Local/Temp/processing_yNrCwH/322fca63ab5546be823900f460cdb91a/OUTPUT.xtf not found Execution failed after 20.51 seconds Loading resulting layers Algorithm 'Stream and Catchment Delineation' finished