[R-sig-Geo] Best approach to implement Spatial Interpolations

Cesar Martinez Izquierdo cesar.martinez at iver.es
Thu Jun 28 08:36:46 CEST 2007


Edzer J. Pebesma escribió:
> Cesar Martinez Izquierdo wrote:
>> Hello list,
>> we are trying to implement spatial interpolations (IDW and several 
>> types of  Kriging) in our open source GIS application, which is 
>> written in Java.
>>
>> We plan to use JRI to connect with R from Java, and then use the 
>> features of "gstat" R-package to calculate the interpolations.
>>
>> However, we don't want R to draw the results of interpolations, we 
>> just want to get the raw result which we will draw in our GIS client.
>>
>> Now, the questions:
>> - Do you think gstat is the best R package for this task?
>>   
> As the author, I will refrain from commenting on this. It has a large 
> user base, a code that has been open source and stable for a long time, 
> and it has several options not found in other packages, notably block 
> kriging, cokriging, and kriging in a local neighbourhood with an 
> efficient search algorithm (quadtrees) for huge/massive data sets.
>> I've read that it doesn't allow not-projected data. Is there a package 
>> which provides similar algorithms and works with not-projected data?
>>   
> Recently, this functionality has been added, only the quadtree algorithm 
> does not work in this case; I'm still in doubt whether quadtrees are 
> useful on a sphere. Consider this code somewhat experimental, but I'm 
> willing to help debug it if needed. It does, for example not check 
> whether covariance functions provided are positive definite on a sphere, 
> and has little (if any) functions specifically dedicated for this 
> purpose. It does compute great circle distances though.
> 
> Another package that seems to deal with unprojected coordinates that you 
> may want to look into is package fields.
>> - Do you have some documentation regarding this topics (more exactly: 
>> documentation about how to get the results of a kriging, etc with JRI)?
>>
>>   
> Perhaps the JGR user interface is an example?
> 
> I'm very interested in your developments and experiences; please keep 
> me/us updated on this mailing list.
> -- 
> Edzer
>> Any comments and suggestions are welcome.
>>
>> Regards,
>>
>>   

Thanks to everyone for the information.
Following your advices, we have finally decided to use JRI+R+gstat.
You will have some news when there is something to show.
And maybe we drop some more questions here in the meantime...

Regards,

-- 
César Martínez Izquierdo
gvSIG development team
http://www.gvsig.gva.es

IVER T.I.
C/ Salamanca 50
46005 Valencia - Spain
http://www.iver.es




More information about the R-sig-Geo mailing list