[R-sig-Geo] [R] Spatstat: tessellation problems.

[Adrian.Baddeley at csiro.au]

Hello all

This is a bug (related to changes in the handling of factors in R)
 
It will be fixed in the next version of spatstat (1.20-4) which will be released shortly

regards
Adrian Baddeley

Neba Funwi-Gabga <fusigabsmail at gmail.com> writes:

> I am trying to do a quadrat count defined by covariate properties in
> spatstat. I have read my elevation raster into R (from ascii) and converted
> to class "im" for use in spatstat.
> I have also created a tessellation of the elevation image into 4 categories
> as follows:

>> >elev.cut<-cut(elevation.im, breaks=elev.quantiles, 
>>   labels=c("Very low","Low","High","Very High"))
>> >elev.tess<-tess(image=elev.cut)

> I would like to do a quadrat count to know how many points lie in each of
> the elevation categories I have created, and as far as I know, this function
> should give that:

>>  >elev.quad.count<-quadratcount(points, tess=elev.tess)

> But when I do this, I get a warning message saying:

> Warning message:
> In quadratcount.ppp(points, tess = aspect.tess) :
>   Tessellation does not contain all the points of X