I understand now. I am working with LiDAR data and I have a problem of
memory. If I subset (ex: 1% of total) i can run Tps function.

is It possibile to set a radius search as in IDW?

gianni


On Sun, Oct 3, 2010 at 6:10 PM, Edzer Pebesma <edzer.pebesma@uni-muenster.de
> wrote:

> Sorry this went to the wrong thread.
>
> I don't know. I can't reproduce your example, and a simple example on
> the meuse data set:
>
> library(sp)
> library(fields)
> data(meuse)
> Tps(data.frame(meuse$x,meuse$y), meuse$zinc)
>
> seemed to work. You could provide a reproducable example, or you could
> contact the author(s) of the fields package, I believe they're not
> listening to this list.
> --
> Edzer
>
>
> On 10/03/2010 05:51 PM, gianni lavaredo wrote:
> > Dear Researchers,
> >
> >
> > I have a problem with Tps function. I tried several days and I am worked
> > saturday night and I am working sunday, but I have always the same Error.
> I
> > didn't find a solution in the web sorry.
> >
> > I have a "SpatialPointsDataFrame" with X, Y and H value
> >
> >> class( input.train)
> > [1] "SpatialPointsDataFrame"
> > attr(,"package")
> > [1] "sp"
> >
> > # extract X and Y
> > x <- input.train@coords[,1]
> > y <- input.train@coords[,2]
> > xy <- data.frame(x,y)
> >
> > # extract the value
> > v <- input.train@data[,c("H")]
> >
> > Tps(xy,v)
> >
> > Errore in rep(0, n1 * n2) : argomento 'times' non valido
> > Errore in out$W2 %d*% do.call(out$cov.function.name, c(out$args, list(x1
> =
> > out$knots,  :
> >   error in evaluating the argument 'y' in selecting a method for function
> > '%d*%'
> > Errore in t(out$W2 %d*% do.call(out$cov.function.name, c(out$args,
> list(x1 =
> > out$knots,  :
> >   error in evaluating the argument 'x' in selecting a method for function
> > 't'
> >
> >
> > Sorry and thanks for Help
> > Gianni
> >
> >       [[alternative HTML version deleted]]
> >
> > _______________________________________________
> > R-sig-Geo mailing list
> > R-sig-Geo@stat.math.ethz.ch
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
> --
> Edzer Pebesma
> Institute for Geoinformatics (ifgi), University of Münster
> Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
> 8333081, Fax: +49 251 8339763  http://ifgi.uni-muenster.de
> http://www.52north.org/geostatistics      e.pebesma@wwu.de
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

	[[alternative HTML version deleted]]

