## Not run:
# # 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"))
# ## End(Not run)
Run the code above in your browser using DataLab