Bug report #14238

Processing v.surf.rst / v.surf.rst.line not working

Added by Pedro Venâncio about 4 years ago. Updated about 4 years ago.

Status:Closed
Priority:Normal
Assignee:Giovanni Manghi
Category:Processing/GRASS
Affected QGIS version:master Regression?:No
Operating System: Easy fix?:No
Pull Request or Patch supplied:No Resolution:
Crashes QGIS or corrupts data:No Copied to github as #:22235

Description

Processing GRASS modules v.surf.rst / v.surf.rst.line don't work, both with GRASS6 or GRASS7.

The log message:

Algorithm v.surf.rst.line - Spatial approximation and topographic analysis using regularized spline with tension. starting...
g.proj -c proj4="+proj=tmerc +lat_0=39.66825833333333 +lon_0=-8.133108333333334 +k=1 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs" 
v.in.ogr min_area=0.0001 snap=-1 dsn="/home/pedro/SIG/Testes/CNFinal" layer="CNfinal_sample" output=tmp1454683157134 --overwrite -o
g.region n=134611.252663 s=133245.488458 e=73164.2735396 w=71504.4906524 res=20
v.surf.rst input="tmp1454683157134" zcolumn="ELEVATION" tension="40" segmax="40" npmin="300" dmin="0.5" dmax="2.5" zmult="1" theta="0" scalex="0" elev=elevdcc5a11accca4be3b5aded058aec2f1b slope=slopedcc5a11accca4be3b5aded058aec2f1b aspect=aspectdcc5a11accca4be3b5aded058aec2f1b pcurv=pcurvdcc5a11accca4be3b5aded058aec2f1b tcurv=tcurvdcc5a11accca4be3b5aded058aec2f1b mcurv=mcurvdcc5a11accca4be3b5aded058aec2f1b --overwrite
g.region rast=elevdcc5a11accca4be3b5aded058aec2f1b
r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=elevdcc5a11accca4be3b5aded058aec2f1b output="/home/pedro/SIG/Testes/CNFinal/test.tif" 
g.region rast=slopedcc5a11accca4be3b5aded058aec2f1b
r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=slopedcc5a11accca4be3b5aded058aec2f1b output="/tmp/processingbe93564e8ce2457b9646e9d175ad31cd/c2444c08a20649938abd323659f35f1c/slope.tif" 
g.region rast=aspectdcc5a11accca4be3b5aded058aec2f1b
r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=aspectdcc5a11accca4be3b5aded058aec2f1b output="/tmp/processingbe93564e8ce2457b9646e9d175ad31cd/7a4cb0683d784e2f927dcb262230bedb/aspect.tif" 
g.region rast=pcurvdcc5a11accca4be3b5aded058aec2f1b
r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=pcurvdcc5a11accca4be3b5aded058aec2f1b output="/tmp/processingbe93564e8ce2457b9646e9d175ad31cd/188e99fb373b49eab4b1e48368ede32f/pcurv.tif" 
g.region rast=tcurvdcc5a11accca4be3b5aded058aec2f1b
r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=tcurvdcc5a11accca4be3b5aded058aec2f1b output="/tmp/processingbe93564e8ce2457b9646e9d175ad31cd/75eb02c2fb01464e9c6cf26ee690c9d9/tcurv.tif" 
g.region rast=mcurvdcc5a11accca4be3b5aded058aec2f1b
r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=mcurvdcc5a11accca4be3b5aded058aec2f1b output="/tmp/processingbe93564e8ce2457b9646e9d175ad31cd/21dd6a7be100411684de1cfc6078e06c/mcurv.tif" 
GRASS commands output:
Cleaning up temporary files ... 
Starting GRASS ... 
Executing '/home/pedro/.qgis2//processing/grass_batch_job.sh' ... 
Welcome to GRASS 6.4.3 (2013) 

GRASS_INFO_MESSAGE(4024,1): Default region was updated to the new projection, but if you have multiple mapsets `g.region -d` should be run in each to update the region from the default 
GRASS_INFO_END(4024,1) 

GRASS_INFO_MESSAGE(4024,2): Projection information updated 
GRASS_INFO_END(4024,2) 

GRASS_INFO_MESSAGE(4025,1): Over-riding projection check 
GRASS_INFO_END(4025,1) 

GRASS_INFO_MESSAGE(4025,2): Layer: CNfinal_sample 
GRASS_INFO_END(4025,2) 

GRASS_INFO_MESSAGE(4025,3): Counting polygons for 43 features... 
GRASS_INFO_END(4025,3) 

GRASS_INFO_MESSAGE(4025,4): Importanto feições do mapa 43... 
GRASS_INFO_END(4025,4) 

GRASS_INFO_MESSAGE(4025,5): ----------------------------------------------------- 
GRASS_INFO_END(4025,5) 

GRASS_INFO_MESSAGE(4025,6): Building topology for vector map ... 
GRASS_INFO_END(4025,6) 

GRASS_INFO_MESSAGE(4025,7): Registering primitives... 
GRASS_INFO_END(4025,7) 

GRASS_INFO_MESSAGE(4025,8): 43 primitives registered 
GRASS_INFO_END(4025,8) 

GRASS_INFO_MESSAGE(4025,9): 3800 vertices registered 
GRASS_INFO_END(4025,9) 

GRASS_INFO_MESSAGE(4025,10): Number of nodes: 83 
GRASS_INFO_END(4025,10) 

GRASS_INFO_MESSAGE(4025,11): Number of primitives: 43 
GRASS_INFO_END(4025,11) 

GRASS_INFO_MESSAGE(4025,12): Number of points: 0 
GRASS_INFO_END(4025,12) 

GRASS_INFO_MESSAGE(4025,13): Number of lines: 43 
GRASS_INFO_END(4025,13) 

GRASS_INFO_MESSAGE(4025,14): Number of boundaries: 0 
GRASS_INFO_END(4025,14) 

GRASS_INFO_MESSAGE(4025,15): Number of centroids: 0 
GRASS_INFO_END(4025,15) 

GRASS_INFO_MESSAGE(4025,16): Number of areas: - 
GRASS_INFO_END(4025,16) 

GRASS_INFO_MESSAGE(4025,17): Number of isles: - 
GRASS_INFO_END(4025,17) 

GRASS_INFO_MESSAGE(4025,18): Building topology for vector map ... 
GRASS_INFO_END(4025,18) 

GRASS_INFO_MESSAGE(4025,19): Registering primitives... 
GRASS_INFO_END(4025,19) 

GRASS_INFO_MESSAGE(4025,20): 43 primitives registered 
GRASS_INFO_END(4025,20) 

GRASS_INFO_MESSAGE(4025,21): 3800 vertices registered 
GRASS_INFO_END(4025,21) 

GRASS_INFO_MESSAGE(4025,22): Building areas... 
GRASS_INFO_END(4025,22) 

GRASS_INFO_MESSAGE(4025,23): 0 areas built 
GRASS_INFO_END(4025,23) 

GRASS_INFO_MESSAGE(4025,24): 0 isles built 
GRASS_INFO_END(4025,24) 

GRASS_INFO_MESSAGE(4025,25): Attaching islands... 
GRASS_INFO_END(4025,25) 

GRASS_INFO_MESSAGE(4025,26): Attaching centroids... 
GRASS_INFO_END(4025,26) 

GRASS_INFO_MESSAGE(4025,27): Number of nodes: 83 
GRASS_INFO_END(4025,27) 

GRASS_INFO_MESSAGE(4025,28): Number of primitives: 43 
GRASS_INFO_END(4025,28) 

GRASS_INFO_MESSAGE(4025,29): Number of points: 0 
GRASS_INFO_END(4025,29) 

GRASS_INFO_MESSAGE(4025,30): Number of lines: 43 
GRASS_INFO_END(4025,30) 

GRASS_INFO_MESSAGE(4025,31): Number of boundaries: 0 
GRASS_INFO_END(4025,31) 

GRASS_INFO_MESSAGE(4025,32): Number of centroids: 0 
GRASS_INFO_END(4025,32) 

GRASS_INFO_MESSAGE(4025,33): Number of areas: 0 
GRASS_INFO_END(4025,33) 

GRASS_INFO_MESSAGE(4025,34): Number of isles: 0 
GRASS_INFO_END(4025,34) 

GRASS_INFO_MESSAGE(4028,1): Porcentagem completada: 
GRASS_INFO_END(4028,1) 

GRASS_INFO_MESSAGE(4028,2): Reading features from vector map ... 
GRASS_INFO_END(4028,2) 
Points are too concentrated -- please increase DMIN 

GRASS_INFO_ERROR(4030,1): Raster map not found 
GRASS_INFO_END(4030,1) 

GRASS_INFO_ERROR(4031,1): Raster map or group not found 
GRASS_INFO_END(4031,1) 

GRASS_INFO_ERROR(4032,1): Raster map not found 
GRASS_INFO_END(4032,1) 

GRASS_INFO_ERROR(4033,1): Raster map or group not found 
GRASS_INFO_END(4033,1) 

GRASS_INFO_ERROR(4034,1): Raster map not found 
GRASS_INFO_END(4034,1) 

GRASS_INFO_ERROR(4035,1): Raster map or group not found 
GRASS_INFO_END(4035,1) 

GRASS_INFO_ERROR(4036,1): Raster map not found 
GRASS_INFO_END(4036,1) 

GRASS_INFO_ERROR(4037,1): Raster map or group not found 
GRASS_INFO_END(4037,1) 

GRASS_INFO_ERROR(4038,1): Raster map not found 
GRASS_INFO_END(4038,1) 

GRASS_INFO_ERROR(4039,1): Raster map or group not found 
GRASS_INFO_END(4039,1) 

GRASS_INFO_ERROR(4040,1): Raster map not found 
GRASS_INFO_END(4040,1) 

GRASS_INFO_ERROR(4041,1): Raster map or group not found 
GRASS_INFO_END(4041,1) 
Closing monitors ... 
Cleaning up temporary files ... 
Batch job '/home/pedro/.qgis2//processing/grass_batch_job.sh' (defined in GRASS_BATCH_JOB variable) was executed. 
Goodbye from GRASS GIS 
Cleaning up temporary files ... 
Starting GRASS ... 
Executing '/home/pedro/.qgis2//processing/grass_batch_job.sh' ... 
Welcome to GRASS 6.4.3 (2013) 

GRASS_INFO_ERROR(4096,1): Raster map not found 
GRASS_INFO_END(4096,1) 

GRASS_INFO_ERROR(4097,1): Raster map or group not found 
GRASS_INFO_END(4097,1) 

GRASS_INFO_ERROR(4098,1): Raster map not found 
GRASS_INFO_END(4098,1) 

GRASS_INFO_ERROR(4099,1): Raster map or group not found 
GRASS_INFO_END(4099,1) 

GRASS_INFO_ERROR(4100,1): Raster map not found 
GRASS_INFO_END(4100,1) 

GRASS_INFO_ERROR(4101,1): Raster map or group not found 
GRASS_INFO_END(4101,1) 

GRASS_INFO_ERROR(4102,1): Raster map not found 
GRASS_INFO_END(4102,1) 

GRASS_INFO_ERROR(4103,1): Raster map or group not found 
GRASS_INFO_END(4103,1) 

GRASS_INFO_ERROR(4104,1): Raster map not found 
GRASS_INFO_END(4104,1) 

GRASS_INFO_ERROR(4105,1): Raster map or group not found 
GRASS_INFO_END(4105,1) 

GRASS_INFO_ERROR(4106,1): Raster map not found 
GRASS_INFO_END(4106,1) 

GRASS_INFO_ERROR(4107,1): Raster map or group not found 
GRASS_INFO_END(4107,1) 
Closing monitors ... 
Cleaning up temporary files ... 
Batch job '/home/pedro/.qgis2//processing/grass_batch_job.sh' (defined in GRASS_BATCH_JOB variable) was executed. 
Goodbye from GRASS GIS 
Converting outputs
Loading resulting layers

The following layers were not correctly generated.
Interpolated RST
You can check the log messages to find more information about the execution of the algorithm

Looking at the log message, it says: Points are too concentrated -- please increase DMIN

I've tried many values for dmin, but I always got the same error.

Sample data attached.

sample_contour_lines.zip (49.6 KB) Pedro Venâncio, 2016-02-05 06:47 AM

History

#1 Updated by Giovanni Manghi about 4 years ago

  • Target version deleted (Version 2.14)
  • Assignee changed from Victor Olaya to Giovanni Manghi
  • Status changed from Open to Feedback

Hi Pedro, just tested with your dataset.... for me for example dmin=50 works fine (not tested other values).

#2 Updated by Pedro Venâncio about 4 years ago

  • Status changed from Feedback to Closed

Hi Giovanni,

You are right, I've tested with several dmin values and it works with values from 17.

Testing the complete dataset this value needs to be increased further.

So I think that this is not actually a bug.

Sorry for the noise! Closing...

Also available in: Atom PDF