how about this...
power<- seq(1.5, 4.0, by=0.1)
for (i in 1:length(power)){
rsaga.grid.calculus(in.grids = c("DCM_1.sgrd"),
out.grid = paste("DCM_1_power",i,".sgrd",sep=""),
formula = paste("a^",power[i],sep=""))
}
Alessandro wrote:
Hi Thanx Alex and Mike
SORRY I am
Hi Thanx Alex and Mike
SORRY I am a brand-new
I tried this formula
for (pow in seq(1.5, 4.0, by = 0.1))
rsaga.grid.calculus(in.grids = c("DCM_1.sgrd"),
out.grid = "DCM_1_power15.sgrd",
formula = paste("a^", pow, sep=""))
but the code create only one raster (DCM_1_power
