[R-sig-Geo] gstat for interpolation

[Ashton Shortridge]

Dear Minti,

with only 8 stations there is not enough data to construct reliable covariance 
(variogram) models. Eight stations simply isn't very many for reliably 
modeling spatial structure. You would be much better off just using inverse 
distance weighting (IDW), which gstat can also do, or even nearest-neighbor 
interpolation. 

Second, it would be best to project the latitudes and longitudes to a planar 
coordinate system appropriate for your study region. I think gstat may be able 
to deal with spherical distances, but I haven't any experience with that. The 
rgdal package is one of several that can do this in R; it can also be done in 
GIS or with other geospatial tools.

Yours,

Ashton


On 09/07/12, Mintewab Bezabih, wrote:
> Dear the R-sig-geo users,
> 
> 
> I have been trying to interpolate rainfall and temperature data for 8
> stations and that needs to be interpolated at a farm level so that I
> have farm level temperature and rainfall observations.
> 
> I am trying to use R for that. First off is kriging in gstat something
> that I can use for that purpose? I have latitude and longitude
> information. Can somebody give me some tips on how I can go about
> that?
> 
> many thanks in advance
> minti
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo


-----
Ashton Shortridge
Associate Professor			ashton at msu.edu
Dept of Geography			http://www.msu.edu/~ashton
235 Geography Building		ph (517) 432-3561
Michigan State University		fx (517) 432-1671