search for: quadratcount

Hello, I am trying to plot a quadratcount object over a ppp object in the spatstat package. I ultimately want to get something like this http://bg9.imslab.co.jp/Rhelp/R-2.4.0/src/library/spatstat/man/images/big_quadratcount_001.png http://bg9.imslab.co.jp/Rhelp/R-2.4.0/src/library/spatstat/man/images/big_quadratcount_001.png See http:/...

Hi, I know there are older threads discussing the quadratcount function in spatstat. Unfortunately, I could not find a solution to my problem there. I'm analyzing a point pattern in an irregular polygonal window. Both the window (an entire country) and the points are projected using WGS84. When I do quadratcount with only one quadrat for the entire cou...

Hi all, I am using quadratcount in spatstat to divide a window containing a point pattern into a grid of quadrats containing the intensity of points in each quadrat. My data is in UTM co-ordinates. My window is defined as follows: >p15<-ppp(x,y,window=owin(c(341710,342100),c(3126465,3126780)),marks=NUL L, checks=TRUE) G...

Hi, Using the quadrats and quadratcount functions in spatstat, when I go to plot either of these, I get the quadrats coloured by their identity, i.e., using a color ramp applied to the sequence of quadrats. This only happens when the quadrats are applied to an owin which is polygonal, i.e., when I have an irregularly shaped study ar...

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 ?...

I have photographs of plots that look like so: http://r.789695.n4.nabble.com/file/n4631960/Untitled.jpg I need to divide it up so each circle has an equal area surrounding it. So into 20 equal segments, each of which contains a circle. Quadratcount is not sufficient because if I divide it up into 36 equal quadrats, some quadrats do not contain one of the circles. I'm not even sure how to do it mathematically, let alone using R. Can anyone help? -- View this message in context: http://r.789695.n4.nabble.com/Quadrat-counting-with-spat...

...pixel", "pixels") c <- im(b, xcol=seq(ncol(b)), yrow=seq(nrow(b)), unitname="pixels") d <- tess(image = c) plot(d) #I can then plot my ppp pattern on top of my tessellation using: plot(d) plot(data02, add=TRUE) #So far so good, but when I run for example: qb <- quadratcount(data02, tess = d) #I get the following error: Warning message: In quadratcount.ppp(data02, tess = d) : Tessellation does not contain all the points of X #When I investigate further the following is returned for each object: > d Tessellation Tessellation is determined by a factor-valued im...

...ow","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 I have checked and the extent of the tessellated elevation image is different from the original eleva...

...ow","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 I have checked and the extent of the tessellated elevation image is different from the original eleva...

...nty so it had a defined border which was somewhat rectangular. Could it be because the shapefile I am using for NYC has too many borders? Is it the coords code that is incorrect? Any advice or suggestions would be greatly appreciated. In case it is helpful, the code that follows is: qtype_1<-quadratcount(type_1,10,10) plot(type_1,pch="+",main=Title") plot(type_1,add=TRUE,col="red",ltw=2,cex=1.2) polygon(coords,lwd=2) points(type_1,pch=".") Thanks, Abby

...have a huge amount of points over a 2D space. I divide my space in a grid with a given resolution,say, 100m. The main loop that I am not sure if it's mandatory or not (any alternative is welcomed) is to go through EACH cell/pixel that contains at least 2 points (right now I am using the method quadratcount within the package spatstat). Inside this loop, thus for each one of this non empty cells, I have to find and keep only a maximum of 10 Male-Female pairs that are within 3 meters from each other. The 3-meter buffer can be done using the "disc" function within spatstat. To select points f...

...rantastic) updates this week (Crantastic)   51. Re: Sapply (Rolf Turner)   52. Re: Sapply (hadley wickham)   53. Aggregating irregular time series (R_help Help)   54. Re: SVM coefficients (Steve Lianoglou)   55. Re: [R-SIG-Finance] Aggregating irregular time series       (Gabor Grothendieck)   56.  Quadratcount image plot colors (phil smart)   57. RE xcel - two problems - Running RExcel & Working Directory       error (Raoul)   58. Re : Re:  Best R text editors? (Wolfgang RAFFELSBERGER)   59. clarificatin on validate.ols method='cross' (Dylan Beaudette)   60. Re: RE xcel - two problems - Runni...