[R-sig-Geo] SpatStat: creating window from quadrats

[marcelino.delacruz at upm.es]

Hi Kevin,

My answer with an example:

 W <- square(10) # the big window
 Z <- quadrats(W, 4, 5) # the quadrats

 # point pattern of boarddamage in the big window:
 boardd <- rpoispp(0.1, win=W)

# A splited point pattern. Each element is a quadrat with at least one
point
 boardquadrats <- split(boardd, Z, drop=T)

Cheers, Marcelino





Con fecha 3/8/2011, "Kévin Morelle" <Kevin.Morelle at ulg.ac.be> escribió:

>Dear R-sig-geo list members,
>
>I'm using the SpatStat R packages for analysing patterns of wild boar 
>damages, considered as my point pattern dataset.
>This dataset is quite large, covering almost whole Belgium. I would like 
>to make my analysis by quadrats of 10x10km. I thus created quadrats 
>covering my points dataset using quadrats function.
>There come my question. I would like to retain as my window, only the 
>quadrats that contain damage points and remove the one that are empty. I 
>have the feeling this question is trivial but I couldn't find any way to 
>do that (apart from creating a grid shape in GIS, but would prefer not...).
>Thanks for your help,
>Regards,
>
>Kevin
>
>-- 
>Morelle Kevin
>Université de Liège - Gembloux Agro-Bio Tech
>Unité de Gestion des Ressources forestières et des Milieux naturels
>--------
>University of Liège - Gembloux Agro-Bio Tech
>Unit of Forest and Nature Management
>--------
>2, Passage des Déportés
>B-5030  GEMBLOUX
>BELGIQUE
>
>Tél : + 32 (0)81 62.22.14
>Fax : + 32 (0)81 62.23.01
>E-mail : kevin.morelle at ulg.ac.be
>Web : http://www.fsagx.ac.be/gf/
>
>_______________________________________________
>R-sig-Geo mailing list
>R-sig-Geo at r-project.org
>https://stat.ethz.ch/mailman/listinfo/r-sig-geo