...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, TargetGrid)

Sorry for spamming. 


On Mon, 12 Feb 2018 08:03:41 +0000
Stefano Sofia <stefano.so...@regione.marche.it> wrote:

> 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 snowfall values on mountain regions, and in this 
> case altitude is an important factor, the use of a Digital Elevation Model 
> might make the difference.
> Looking in the web I found some important material about accessing elevation 
> data in R with the "elevatr" package (by J.W.Hollister), I am reading it.
> My concern would then be to interpolate snowfall values on a Digital 
> Elevation Model. Did somebody already use R for these kinds of calculations? 
> Could somebody share with me some useful hints?
> Thank you for your attention
> Stefano Sofia
>          (oo)
> --oOO--( )--OOo----------------
> Stefano Sofia PhD
> Area Meteorologica e  Area nivologica - Centro 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----------------
> ________________________________
> AVVISO IMPORTANTE: Questo messaggio di posta elettronica può contenere 
> informazioni confidenziali, pertanto è destinato solo a persone autorizzate 
> alla ricezione. I messaggi di posta elettronica per i client di Regione 
> Marche possono contenere informazioni confidenziali e con privilegi legali. 
> Se non si è il destinatario specificato, non leggere, copiare, inoltrare o 
> archiviare questo messaggio. Se si è ricevuto questo messaggio per errore, 
> inoltrarlo al mittente ed eliminarlo completamente dal sistema del proprio 
> computer. Ai sensi dell’art. 6 della DGR n. 1394/2008 si segnala che, in caso 
> di necessità ed urgenza, la risposta al presente messaggio di posta 
> elettronica può essere visionata da persone estranee al destinatario.
> IMPORTANT NOTICE: This e-mail message is intended to be received only by 
> persons entitled to receive the confidential information it may contain. 
> E-mail messages to clients of Regione Marche may contain information that is 
> confidential and legally privileged. Please do not read, copy, forward, or 
> store this message unless you are an intended recipient of it. If you have 
> received this message in error, please forward it to the sender and delete it 
> completely from your computer system.
> --
> This message was scanned by Libra ESVA and is believed to be clean.
>       [[alternative HTML version deleted]]
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo


Dr. Daniel Knitter
CRC1266 -- Scales of Transformation
Project A2 "Integrative Modeling of Socio-Environmental Dynamics"

Christian-Albrechts-Universität zu Kiel
Department of Geography
Physical Geography -- Landscape Ecology and Geoinformation
Ludewig-Meyn-Str. 14, room 109
D - 24118 Kiel, Germany

Phone: +49 431 880 2941
Email: knit...@geographie.uni-kiel.de

Attachment: pgpRJnJh9g8zo.pgp
Description: OpenPGP digital signature

R-sig-Geo mailing list

Reply via email to