Try the chull function to create a border.
On Thu, Oct 1, 2009 at 12:39 AM, Freddy López <freddy.vat...@gmail.com> wrote: > Hello folks, > > I'm working with geostatistics data and I'm using mainly geoR library. > Now, in some examples, is useful to use the 'borders' to delimitate > our kriging process, e.g., in the krige.conv() function. To steal a > Ribeiro & Diggle's example > (http://leg.ufpr.br/geoR/tutorials/aula050819.R), we have some lines > such as > > >>Ksat.in <- krige.conv(Ksat, loc=grid0, borders=Ksat$borders,krige= > krige.control(obj=Ksat.ef)) >>image(Ksat.in, col=gray(seq(1, 0.1, l=21))) > > > in which is nice to have Ksat$borders because we can have our > prediction inside of these borders (and it looks better). > > Now, suppose we have an object of the class geodata but it has not got > borders. > > We can visualisate I'm saying with a little modification of latter example > > >> Ksat.in <- krige.conv(coords=Ksat$coords,data=Ksat$data, > loc=grid0,krige=krige.control(obj=Ksat.ef)) >> image(Ksat.in, col=gray(seq(1, 0.1, l=21))) > > > which looks as an rectangle with no reference to something. > > That's my problem: I have an object of the class geodata (indeed, is a > worked shapefile, because I'm treating with counties) and I have done > all good up to the kriging stage because I don't know how to extract > or create borders. I mean, I can get a rectangular kriging, but it is > not the aim. > > Again, excuse my clumsy english. > > Thanks. > > > -- > 'qué importa la tristeza si hubo en el tiempo > alguien que se dijo feliz,' > > JL Borges > > [[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 > > -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40" S 49° 16' 22" O _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo