Displaying 1 result from an estimated 1 matches for "acmi50".
Did you mean:
acmi2
2004 Mar 04
0
Extracting Krig results
...Frankly, all I really need is to query the Krig on a specific XY pair
and return a result. I can carry the rest. I am using tps from the
fields package.
Here's the code I use to create the Krig object:
=========== Snip! ===============================
dCurr <- sqlQuery(ds,"select x,y,acmi50 from public.cc_output" )
fitCurr <- Tps(dCurr[1:2],dCurr$acmi50, scale.type="unscaled")
surface(fitCurr,type="C", extrap=TRUE, main="Median ACMI, 1961-1990",
xlab="", ylab="", levels=evens)
===============================================
U...