Spatial data interpolation on R
??Dear all,
??Thank you very much for those who read and tried to help. However, I'm really a begginer on R and especially geostatistics.?I'm trying to follow "gstat" package's guide, but I'm facing some trouble when trying to reproduce the commands created for the meuse dataset with my own dataset.
??I wonder if anybody could just help me in this task. In order to better illustrate my problem, I am attaching a figure which shows my spatial domain. The greater contour is Rio Grande do Sul, Brazilian southernmost state. The highlighted region in its interior represents the largest soybean producer area (my area of study). The finer grid are satellite precipitation values and the red dots are the
station.
?I'm also attaching the two files that contains the gauge and cmorph (satellite) data I mentioned (both files in ascii format). As you can note, I've already organized them in the following structure:
LON ? ? ? ? ? ? ? ? ? ? ? ? ? ? LAT ? ? ? ? ? ? ? ? ? ? ? ? ? ? precipitation record
??Satellite data spatial ranges from longitude -55.6951 to -50.0929 and latitude -29.6907 to -27.1437, with increments ("grid size") of 8km (0.072?, as you can note in the file). Station data varies from longitude -55.2670 to -50.0660 and latitude -29.2680 to -27.1920, disposed randomly.
??Thus, all I need is to interpolate the station values into the same grid of cmorph (satellite) in order to perform a proper comparison between these two sources. Quite challenging problem, huh?
??Could you help me to do this using gstat on R?? Please let me know if I have skipped any
detail...
??Best wishes,
??Thiago
--- On Wed, 2/6/10, jgarcia at ija.csic.es <jgarcia at ija.csic.es> wrote:
From: jgarcia at ija.csic.es <jgarcia at ija.csic.es> Subject: Re: [R-sig-Geo] Spatial data interpolation on R To: "Thiago Veloso" <thi_veloso at yahoo.com.br> Cc: r-sig-geo at stat.math.ethz.ch Date: Wednesday, 2 June, 2010, 15:09 You should conduct a block kriging from the point sparse data to the regular grid (the domain of the satellite images). Try, e.g., gstat Javier ///
Dear R colleagues! ?I?d like to start my participation in this list by describing my current problem: inside my area of study I need to compare precipitation data from two different
sources: both station (total of 86) and a grid (at 8km) of
satellite estimates. ??My specific objective is to interpolate the station data into a regular grid in the same resolution of the satellite estimates, preferentially having control of the spatial domain (lat/lon coordinates). As far as I know this is the correct way of making such comparison. ??Could anybody please point directions to perform this task using R? I?m such a beginner that I don?t even know if ? there?s a package designed to create regular grids from "random" data (interpolating by kriging or other technique). Usage examples would be welcomed as well! Thanks in advance, ?Thiago. ??? [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo
mailing list
R-sig-Geo at stat.math.ethz.ch
-------------- next part -------------- An HTML attachment was scrubbed... URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20100604/34e23d0b/attachment.html> -------------- next part -------------- A non-text attachment was scrubbed... Name: example.png Type: image/png Size: 29202 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20100604/34e23d0b/attachment.png> -------------- next part -------------- A non-text attachment was scrubbed... Name: gauge_09.csv Type: text/csv Size: 2858 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20100604/34e23d0b/attachment.bin> -------------- next part -------------- A non-text attachment was scrubbed... Name: cmorph_09.csv Type: text/csv Size: 41088 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20100604/34e23d0b/attachment-0001.bin>