� 13 febbraio 2018 17.21
A: Stefano Sofia
Cc: Daniel Knitter; r-sig-geo@r-project.org
Oggetto: Re: [R-sig-Geo] Interpolating snowfall values on a Digital Elevation
Model
You can't use a lat/long coordinate system when kriging because the concept of
distance is ambiguous. Convert all your data
tro Funzionale
> Servizio Protezione Civile - Regione Marche
> Via del Colle Ameno 5
> 60126 Torrette di Ancona, Ancona
> Uff: 071 806 7743
> E-mail: stefano.so...@regione.marche.it
> ---Oo-oO
>
> Da: R-sig-Geo
niel Knitter
[knit...@geographie.uni-kiel.de]
Inviato: lunedì 12 febbraio 2018 9.16
A: r-sig-geo@r-project.org
Oggetto: Re: [R-sig-Geo] Interpolating snowfall values on a Digital Elevation
Model
...sorry, I missed one step: you need to extract the elevation values to your
points via ex
...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,
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,
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