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/