Re: [R-sig-Geo] I: Interpolating snowfall values on a Digital Elevation Model
_data, grd) > > > --- > > > > The error I get is: > > Error in autoKrige(Cumulata ~ rain_data$ExtractedElevationValues, > rain_data, : > > Either input_data or new_data is in LongLat, please reproject. > >input_data: +init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs > +ellps=WGS84 +towgs84=0,0,0 > >new_data:+init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs > +ellps=WGS84 +towgs84=0,0,0 > > > > but I did pay attention to have the same reference system in rain_data, > drg and the Digital Elevation Model. > > > > Moreover, if I impose the class of the DEM to SpatialGridDataFrame when > I extraxt the elevation points from the DEM I get the following error: > > Error in (function (classes, fdef, mtable) : > > unable to find an inherited method for function ‘extract’ for > signature ‘"SpatialGridDataFrame", "SpatialPoints"’ > > Calls: extract -> > > > > > > Would you please somebody help to show me where is my mistake? > > > > Thank you for all your attention > > Stefano > > > > (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 > > > > Da: R-sig-Geo [r-sig-geo-boun...@r-project.org] per conto di Daniel > 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 extract function from the raster package: > > > > YourPoints$YourExtractedElevationValues <- extract(x = DEM, y = > YourPoints) > > > > autoKrige(YourMeasurements ~ YourExtractedElevationValues, > YourMeasurementLocations, TargetGrid) > > > > > > Sorry for spamming. > > > > Best, > > Daniel > > > > On Mon, 12 Feb 2018 08:03:41 + > > 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 pr
[R-sig-Geo] I: Interpolating snowfall values on a Digital Elevation Model
nts from the DEM I get the following error: > Error in (function (classes, fdef, mtable) : > unable to find an inherited method for function ‘extract’ for signature > ‘"SpatialGridDataFrame", "SpatialPoints"’ > Calls: extract -> > > > Would you please somebody help to show me where is my mistake? > > Thank you for all your attention > Stefano > > (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 > > Da: R-sig-Geo [r-sig-geo-boun...@r-project.org] per conto di Daniel 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 extract function from the raster package: > > YourPoints$YourExtractedElevationValues <- extract(x = DEM, y = YourPoints) > > autoKrige(YourMeasurements ~ YourExtractedElevationValues, > YourMeasurementLocations, TargetGrid) > > > Sorry for spamming. > > Best, > Daniel > > On Mon, 12 Feb 2018 08:03:41 + > 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://urlsand.esvalabs.com/?u=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-geo=52342f8