...sorry, I missed one step: you need to extract the elevation values to your
points via extract function from the raster package:
YourPoints$YourExtractedElevationValues <- extract(x = DEM, y = YourPoints)
autoKrige(YourMeasurements ~ YourExtractedElevationValues,
YourMeasurementLocations, Tar
Hi Stefano,
assuming you have a DEM (downloaded e.g. from https://earthexplorer.usgs.gov/
or using the getData function from the raster package) as SpatialGridDataFrame
loaded via readGDAL from the rgdal package you can use
autoKrige(YourMeasurements ~ DEM, YourMeasurementLocations, TargetGrid
Dear list users,
I have to produce rainfall maps. I know how to create a (bi-dimensional) grid
and interpolate rainfall values (from automatic rain gauges) on that grid using
Kriging:
myinterpolation <- autoKrige(myrainfall_cumulate ~ 1, rain_data_UTM, mygrid)
In reality I am dealing with snowf