similar to: tessellation from biological data in spatstat

Hello R Users, 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. Now I have point data of class "ppp" which window is the same extent as the elevation image. I can conveniently plot the image and the points on the image as follows:

Hello R Users, 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. Now I have point data of class "ppp" which window is the same extent as the elevation image. I can conveniently plot the image and the points on the image as follows:

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 area. There doesn't seem to be any obvious way to

This is enough to enable tessellation support on nvc0. It seems to work a lot better on my GF108 than GK208. I suspect that there's some sort of scheduling shenanigans that need to be adjusted for kepler+. Or perhaps some shader header things. Even with the GF108, I still get occasional blue triangles in Heaven, but I get a *ton* of them on the GK208 -- seemingly the same issue, but it's

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

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 country it holds all my 154 points:

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) Giving me a distance of 390m in the 'x' direction and 315m in

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

There's Renka's STRIPACK, and TRIPACK, respectively, ACM TOMS Algorithms 772 and 751, and there's the R package "deldir" which does the Delaunay for a plane, but does anyone have or know of the tessellation in R for a sphere? Also, is there a standard indexing scheme for Delaunay facets, and perhaps of edges in such facets? I'd expect that to be a publication reference,

Hello, I've been debugging a few different tessellation shader issues with nouveau, but let's start small. I see this issue on my GK208 with high frequency, and I *think* I've seen it once or twice on my GF108, but it's exceedingly rare, if it does happen. I don't have a GK10x to test on, unfortunately, but I assume it'll have the same issue as the GK208. The issue is

Dear R users, I wish to run spatial point pattern analysis (e.g. pair correlation function, mark correlation function) for which I need to create an observation window (window=owin) from which the spatial analysis is generated. The command I used to create this observation window as follows: X1<- ppp(x, y, window=owin(c(80.58,144.96),c(101.06,165.13)),

One additional observation that I just made is that on GK208, the blob apparently doesn't use the result of S2R Rx, SR_INVOCATION_ID wholesale in TCS. It either passes it through a I2I.S32.S32 Rx, |Rx| (i.e. absolute value), or even more paradoxically, shl 2; shr 2; which removes the top *2* bits, rather than just the top 1. However I see no such behaviour on GF108. I'm going to test out

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

Okay, been working with tripack, seems the most mature package for this. Got it to work well with their test data set - data(tritest). When i tried random numbers to explore further, i am getting some results that don't reconcile. example run this: library(tripack) y <- runif(100) x <- runif(100) vm <- voronoi.mosaic(x,y) plot(vm) par(new=T) plot(x,y,col='blue') when

On 02/06/17 01:17, Lluis.Hurtado at uv.es wrote: > Dear all, > > I am currently working with the spatstat package, using windows and pixel images. > > First: > > My aim is to transform a shapefile (see attached) into a pixel image. > > My idea is to start transforming the shapefile into a Spatial Polygon file: > > x <-

Dear All; Is it possible to create a list of lists (I am sure it is) along these lines; I have a dataframe data02 that holds a lot of information, and the first column is ³date² I have a list of dates in; data03<-c(date1,.....,daten) And would like to create a list; data04 <- subset(data02, date == data03[1,]) Ie. data04 holds the data from data02 that matches a date in data03 How do

Hi, i want to parse a XML-File. I made some Tutorial but with my special Format it don't work. An Example of my format: <?xml version="1.0" encoding="ISO-8859-1"?> <mzML xmlns="http://psi.hupo.org/ms/mzml" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://psi.hupo.org/ms/mzml

hi! sometimes we've on some client-servers hanging uninterruptible processes ("ps aux" stat is on "D" ) and on one the CPU wait I/O grows within some minutes to 100%. you are not able to kill such processes - also "kill -9" doesnt work - when you connect via "strace" to such an process, you wont see anything and you cannot detach it again. there

Backtrace: ... 0x00007ffff57fc392 in __assert_fail () from /lib64/libc.so.6 0x00007ffff0cf5bec in nvc0_shader_stage (pipe=<optimized out>) at ./nvc0/nvc0_context.h:204 nvc0_set_constant_buffer (pipe=0x631080, shader=<optimized out>, index=<optimized out>, cb=0x0) at nvc0/nvc0_state.c:771 0x00007ffff0a2cf63 in st_destroy_context (st=0x68a9f0) at state_tracker/st_context.c:382

I've recently setup a distributed/replicated cluster and have had an issue with seeing the directories on the cluster. Also, a df -h only shows data from one of the three bricks. The strange behavior doesn't end there. If I log into the 'primary' server as root, then do an ls on the client, the directories appear. However, df -h is still incorrect. I'm not sure exactly