if (FALSE) {
# IDW interpolation of attribute "z" from the point shapefile
# 'points.shp' to a grid with the same extent and resolution
# as the (pre-existing) geology grid:
rsaga.inverse.distance("points", "dem", field = "z", maxdist = 1000,
target = rsaga.target(target="target.grid",
target.grid = "geology"))
}
Run the code above in your browser using DataLab