[R-sig-Geo] Creating borders.-

[Henrique Dallazuanna]

Try the chull function to create a border.


On Thu, Oct 1, 2009 at 12:39 AM, Freddy López <freddy.vate01 at 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 at 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