[Qgis-developer] SEXTANTE - using calculator tools from SAGA and
GRASS
Pedro Venâncio
pedrongvenancio at yahoo.com
Wed May 2 13:22:40 PDT 2012
Hi,
I'm trying to use the calculator tools both from SAGA and GRASS in SEXTANTE, and both are giving me errors.
For example, making the difference between two rasters in SAGA, with Grid Difference module, all runs ok.
Calculating the difference between the same rasters, but with Grid Calculator module, gives me an error. Looking at the log files, I noticed that the only difference is in the name of the input files, that in the case of Grid Difference uses the .sgrd temporary files, and in the case of Grid Calculator uses the name of the input files, which I put in the formula:
Grid Calculator:
SAGA execution commands
libio_gdal 0 -GRIDS "/home/pedro/sextante/tempdata/1335987049.2352.sgrd" -FILES "/home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif"
libgrid_tools "Resampling" -INPUT "/home/pedro/sextante/tempdata/1335987049.2352.sgrd" -TARGET 0 -SCALE_UP_METHOD 4 -SCALE_DOWN_METHOD 4 -USER_XMIN 67004.8257193 -USER_XMAX 113984.825719 -USER_YMIN 93577.2639045 -USER_YMAX 146317.263904 -USER_SIZE 30.0 -USER_GRID "/home/pedro/sextante/tempdata/1335987049.2353.sgrd"
libio_gdal 0 -GRIDS "/home/pedro/sextante/tempdata/1335987049.2354.sgrd" -FILES "/home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2.tif"
libgrid_tools "Resampling" -INPUT "/home/pedro/sextante/tempdata/1335987049.2354.sgrd" -TARGET 0 -SCALE_UP_METHOD 4 -SCALE_DOWN_METHOD 4 -USER_XMIN 67004.8257193 -USER_XMAX 113984.825719 -USER_YMIN 93577.2639045 -USER_YMAX 146317.263904 -USER_SIZE 30.0 -USER_GRID "/home/pedro/sextante/tempdata/1335987049.2355.sgrd"
libgrid_calculus "Grid Calculator" -GRIDS "/home/pedro/sextante/tempdata/1335987049.2353.sgrd;/home/pedro/sextante/tempdata/1335987049.2355.sgrd" -FORMULA "mdt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif-mdt_aster_pt_D73_GeoTIFF_clip2" -RESULT "/home/pedro/sextante/tempdata/1335987049.250.tif.sgrd"
libio_gdal 1 -GRIDS "/home/pedro/sextante/tempdata/1335987049.250.tif.sgrd" -FORMAT 1 -TYPE 0 -FILE "/home/pedro/sextante/tempdata/1335987049.250.tif"
Grid Difference:
SAGA execution commands
libio_gdal 0 -GRIDS "/home/pedro/sextante/tempdata/1335987103.7557.sgrd" -FILES "/home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif"
libgrid_tools "Resampling" -INPUT "/home/pedro/sextante/tempdata/1335987103.7557.sgrd" -TARGET 0 -SCALE_UP_METHOD 4 -SCALE_DOWN_METHOD 4 -USER_XMIN 67004.8257193 -USER_XMAX 113984.825719 -USER_YMIN 93577.2639045 -USER_YMAX 146317.263904 -USER_SIZE 30.0 -USER_GRID "/home/pedro/sextante/tempdata/1335987103.7558.sgrd"
libio_gdal 0 -GRIDS "/home/pedro/sextante/tempdata/1335987103.7559.sgrd" -FILES "/home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2.tif"
libgrid_tools "Resampling" -INPUT "/home/pedro/sextante/tempdata/1335987103.7559.sgrd" -TARGET 0 -SCALE_UP_METHOD 4 -SCALE_DOWN_METHOD 4 -USER_XMIN 67004.8257193 -USER_XMAX 113984.825719 -USER_YMIN 93577.2639045 -USER_YMAX 146317.263904 -USER_SIZE 30.0 -USER_GRID "/home/pedro/sextante/tempdata/1335987103.7560.sgrd"
libgrid_calculus "Grid Difference" -A "/home/pedro/sextante/tempdata/1335987103.7558.sgrd" -B "/home/pedro/sextante/tempdata/1335987103.7560.sgrd" -C "/home/pedro/sextante/tempdata/1335987103.7555.tif.sgrd"
libio_gdal 1 -GRIDS "/home/pedro/sextante/tempdata/1335987103.7555.tif.sgrd" -FORMAT 1 -TYPE 0 -FILE "/home/pedro/sextante/tempdata/1335987103.7555.tif"
Grid Calculator error:
SAGA execution console output
_____________________________________________
##### ## ##### ##
### ### ## ###
### # ## ## #### # ##
### ##### ## # #####
##### # ## ##### # ##
_____________________________________________
library path:/usr/lib/saga
library name:libio_gdal
module name :GDAL: Import Raster
author :(c) 2007 by O.Conrad (A.Ringeler)
_____________________________________________
go...
Parameters
Grids: No objects
Files: "/home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif"
loading: /home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif
Driver: /home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif/
Cells: x 1566, y 1758
Bands: 1
Transformation x' = 67004.825719 + x * 30.000000 + y * 0.000000
Transformation y' = 146317.263904 + x * 0.000000 + y * -30.000000
/
Save grid: /home/pedro/sextante/tempdata/1335987049.2352.sgrd...
ready
okay
_____________________________________________
##### ## ##### ##
### ### ## ###
### # ## ## #### # ##
### ##### ## # #####
##### # ## ##### # ##
_____________________________________________
library path:/usr/lib/saga
library name:libgrid_tools
module name :Resampling
author :O.Conrad (c) 2003
_____________________________________________
go...
Load grid: /home/pedro/sextante/tempdata/1335987049.2352.sgrd...
ready
okay
Parameters
Grid system: 30.000000; 1566x 1758y; 67019.825719x 93592.263904y
Grid: 1335987049.2352
Additional Grids: No objects
Additional Grids: No objects
Preserve Data Type: no
Target Grid: user defined
ready
Save grid: /home/pedro/sextante/tempdata/1335987049.2353.sgrd...
ready
okay
_____________________________________________
##### ## ##### ##
### ### ## ###
### # ## ## #### # ##
### ##### ## # #####
##### # ## ##### # ##
_____________________________________________
library path:/usr/lib/saga
library name:libio_gdal
module name :GDAL: Import Raster
author :(c) 2007 by O.Conrad (A.Ringeler)
_____________________________________________
go...
Parameters
Grids: No objects
Files: "/home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2.tif"
loading: /home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2.tif
Driver: /home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2.tif/
Cells: x 1566, y 1758
Bands: 1
Transformation x' = 67004.825719 + x * 30.000000 + y * 0.000000
Transformation y' = 146317.263904 + x * 0.000000 + y * -30.000000
/
Save grid: /home/pedro/sextante/tempdata/1335987049.2354.sgrd...
ready
okay
_____________________________________________
##### ## ##### ##
### ### ## ###
### # ## ## #### # ##
### ##### ## # #####
##### # ## ##### # ##
_____________________________________________
library path:/usr/lib/saga
library name:libgrid_tools
module name :Resampling
author :O.Conrad (c) 2003
_____________________________________________
go...
Load grid: /home/pedro/sextante/tempdata/1335987049.2354.sgrd...
ready
okay
Parameters
Grid system: 30; 1566x 1758y; 67019.825719x 93592.263904y
Grid: 1335987049.2354
Additional Grids: No objects
Additional Grids: No objects
Preserve Data Type: no
Target Grid: user defined
ready
Save grid: /home/pedro/sextante/tempdata/1335987049.2355.sgrd...
ready
okay
error: Error in formula at position 0
[m] dt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif-mdt_aster_pt_D73_GeoTIFF_clip2
function not found
error: executing module [Grid Calculator]
_____________________________________________
##### ## ##### ##
### ### ## ###
### # ## ## #### # ##
### ##### ## # #####
##### # ## ##### # ##
_____________________________________________
library path:/usr/lib/saga
library name:libgrid_calculus
module name :Grid Calculator
author :Copyrights (c) 2003 by Andre Ringeler
_____________________________________________
go...
Load grid: /home/pedro/sextante/tempdata/1335987049.2353.sgrd...
ready
okay
Load grid: /home/pedro/sextante/tempdata/1335987049.2355.sgrd...
ready
okay
Parameters
Grid system: 30; 1567x 1759y; 67004.825719x 93577.263904y
Grids: 2 objects (1335987049.2353, 1335987049.2355))
Grids from different Systems: No objects
Result: Result
Formula: mdt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif-mdt_aster_pt_D73_GeoTIFF_clip2
Name: Calculation
Take Formula: no
Use NoData: no
/
-
error: Grid file could not be opened.
error: input file [/home/pedro/sextante/tempdata/1335987049.250.tif.sgrd]
error: empty input list [GRIDS]
error: executing module [GDAL: Export Raster]
_____________________________________________
##### ## ##### ##
### ### ## ###
### # ## ## #### # ##
### ##### ## # #####
##### # ## ##### # ##
_____________________________________________
library path:/usr/lib/saga
library name:libio_gdal
module name :GDAL: Export Raster
author :O.Conrad (c) 2007
_____________________________________________
go...
Load grid: /home/pedro/sextante/tempdata/1335987049.250.tif.sgrd...
failed
So, what's the syntax I should use in Formula?
The GRASS r.mapcalculator module
GRASS execution commands
g.region n=146317.263904 s=93577.2639045 e=113984.825719 w=67004.8257193 res=30.0
r.in.gdal input="/home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2.tif" band=1 out=tmp13359846747445 --overwrite -o
r.in.gdal input="/home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif" band=1 out=tmp13359846747446 --overwrite -o
r.in.gdal input="/home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif" band=1 out=tmp13359846747447 --overwrite -o
r.in.gdal input="/home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif" band=1 out=tmp13359846747448 --overwrite -o
r.in.gdal input="/home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif" band=1 out=tmp13359846747449 --overwrite -o
r.in.gdal input="/home/pedro/SIG/ASTER_GDEM/mdt_aster_pt_D73_GeoTIFF/mdt_aster_pt_D73_GeoTIFF_clip2_fillnulls.tif" band=1 out=tmp13359846747450 --overwrite -o
r.mapcalculator amap=tmp13359846747445 bmap=tmp13359846747450 cmap=tmp13359846747450 dmap=tmp13359846747450 emap=tmp13359846747450 fmap=tmp13359846747450 formula=amap-bmap outfile=outfile --overwrite
g.region rast=outfile
r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=outfile output="/home/pedro/sextante/tempdata/1335984674.7343.tif"
gives me this error:
GRASS execution console output
Cleaning up temporary files ...
Starting GRASS ...
Executing '/home/pedro/sextante/grass_batch_job.sh' ...
Welcome to GRASS 6.4.1 (2011)
ATENÃÃO: Over-riding projection check
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.in.gdal completo. Raster map <tmp13359846747445> created.
ATENÃÃO: Over-riding projection check
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.in.gdal completo. Raster map <tmp13359846747446> created.
ATENÃÃO: Over-riding projection check
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.in.gdal completo. Raster map <tmp13359846747447> created.
ATENÃÃO: Over-riding projection check
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.in.gdal completo. Raster map <tmp13359846747448> created.
ATENÃÃO: Over-riding projection check
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.in.gdal completo. Raster map <tmp13359846747449> created.
ATENÃÃO: Over-riding projection check
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.in.gdal completo. Raster map <tmp13359846747450> created.
r.mapcalc "outfile" = "(amap-bmap)"
Invalid map <amap>
Invalid map <bmap>
Parse error
ERRO:Calculating outfile. Try expert mode.
ERRO:Raster map <outfile> not found
ERRO:Raster map or group <outfile> not found
Closing monitors ...
Cleaning up temporary files ...
Batch job '/home/pedro/sextante/grass_batch_job.sh' (defined in GRASS_BATCH_JOB variable) was executed.
Goodbye from GRASS GIS
In the module interface it's necessary to choose files for all fields (amap, bmap, cmap, dmap, emap, fmap), but in this case I just need the first and the second, because this calculation only involve two raster files. Could this be the source of the error?
Thank you very much!
Pedro Venâncio
More information about the Qgis-developer
mailing list