[R-sig-Geo] Extract all point in a quadrats by spatstat package

Marcelino de la Cruz Rot marcelino.delacruz at urjc.es
Tue Nov 28 17:49:18 CET 2017


Maybe

cut(X, as.tess(qX))

or

split.ppp(X, as.tess(qX))

?
Cheers,

Marcelino



El 28/11/2017 a las 16:37, Mohammad Tanvir Ahamed via R-sig-Geo escribió:
> Hi,
>
> With the following code i can divides window into quadrats and counts the numbers of points in each quadrat.
>
> library(spatstat)
> X <- runifpoint(50)
> quadratcount(X)
> quadratcount(X, 4, 5)
> quadratcount(X, xbreaks=c(0, 0.3, 1), ybreaks=c(0, 0.4, 0.8, 1))
> qX <-  quadratcount(X, 4, 5)
> plot(X)
> plot(qX, add=TRUE)
>
> But I want to mark each  quadrats  and select/ extract only those points by selecting quadrats id .
> Can anyone plese help me regarding this issue ?
> Thanks in advance .
>
>
>
> Regards.............
> Tanvir Ahamed
> Stockholm, Sweden     |  mashranga at yahoo.com
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo


-- 
Marcelino de la Cruz Rot
Depto. de Biología y Geología
Física y Química Inorgánica
Universidad Rey Juan Carlos
Móstoles España



More information about the R-sig-Geo mailing list