[R-sig-Geo] Spatial data interpolation on R

[Thiago Veloso]

  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
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>




      
-------------- 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>