Displaying 1 result from an estimated 1 matches for "border_limite".
2002 Apr 16
0
still have problem with krige and border option
Dear all,
I would like to estimate rainfall with the krige function. First, I produce a
polygon of my region of interest (where poly_test.txt is a x,y suite of
points defining a polygon, obtained with the grass v.out.ascii command) :
user>border_limite<-read.table("/home/lionel/rwork/poly_test.txt",header=FALSE)
user>polygone<-list(x=border_limite[,1],y=border_limite[,2])
So now I have a list of x,y coordinates to limit my work area.
If I call the function krige with extrap=FALSE like this
user>grid2.krige<-krige(grid2....