similar to: re ference

Hello list, I need to keep track of objects that are related to particular observations. In this case, I need to keep track of polygons that are associated with observations. What I would ideally have is one column of a dataframe hold a "polygonal" object (from the spatstat package). My question: I seem to have managed to do it, but as I haven't read that dataframes are supposed

Hi, How can I, from a single "driver" file, source other files in such a way that I can access their functions with parameters defined in my "driver" file? I wish to do this to avoid creating a single, self-contained but clunky piece of code. I have searched and found functions such as file(), pipe(), open(), found manuals on creating my own packages, etc. but I become

I have a list of polygons generated by the contourLines() command (each object of the list is a list in itself with two objects: a vector of x values, and a vector of y values for each vertex). I wish to convert that list into a gpc.poly object of multiple contours. How do I do this? gpclib apparently has no method of coercing lists into the gpc.poly object type. As well, when I have a

I have used polygon() to mark the confidence limits of a survival curve. In another project, I have used the col parameter to fill my boxplots. The poly() description refers to filling but actually produces hatching (i.e. lines ). boxplot() does truly fill the boxes with a colour or shades of grey (e.g. col="red"). My novices perception of R graphics is: If you can hatch it, you

I am a new user of the R spatstat package and am having problems creating a polygonal observation window with owin(). Code follows: library("maps") library ("sp")` library("spatstat") mass.map <- map("state", "massachusetts:main", fill=T) # This returns a data frame includding x and y components that form a polygon of massachusetts mainland`

Hi, I have the following: for(j in (y.raw+1):(rownum-1)){ valsum<-tstfframed[min.x,j]+tstfframed[min.x,j+1] if(valsum == 1){ cat("valsum loop") int.num<-int.num+1 } } but I get the error message: "Error in if (valsum == 1) { : argument is of length zero". I checked whether to see if all the stated variables have associated values (and are not NAs) and

Hi there, Is there any package which can calculate the overlapping area of two polygons or output the coordinates of the overlapping of two polygons? Thanks. Haiyong

I was just thinking last night, how about this: For each frame, detect edges and divide the frame into a bunch of non-overlapping polygons. Extract textures from the first frame, then save changes in polygon vertex coordinates only, unless the textures change too much in which case new textures (which could be wavelet compressed) would be stored. Also, if the movements of the vertices would

Hola He hecho el típico mapa con gadm y spplot: library(sp) ES <- local(get(load(url("http://biogeo.ucdavis.edu/data/gadm2/R/ESP_adm2.RData ")))) acerco las canarias: A <- which(ES en data$NAME_1 == "Islas Canarias") L <- length(ES en polygons[[A]]@Polygons) L <- length(ES en polygons[[A]]@Polygons) for (i in 1:L){ ES en polygons[[A]]@Polygons[[i]]@coords

Hola, Yo lo resuelvo usando grid.rect. Tienes una posible solución aquí: https://github.com/oscarperpinan/spacetime-vis/blob/master/choropleth.R#L216 (que es una versión actualizada de este artículo https://procomun.wordpress.com/2012/02/18/maps_with_r_1/) Saludos. Oscar. ----------------------------------------------------------------- Oscar Perpiñán Lamigueiro Dpto. Ing. Eléctrica,

Background: I have an algorithm which produces a large number of small polygons (of the spatial kind) which I would like to use within R using objects from sp. I can't predict the exact number of polygons a-priori, the polygons will be grouped into regions, and each region will be filled sequentially, so an appropriate C++ 'framework' (for the point of illustration) might be: typedef

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 list, I have a question regarding the included polygon in the khat function of the splancs library defining the area where points appear. I have not only one simple polygon included, my map includes several islands. I read my map which was a shape file by readShapeSpatial of the maptools library. Then I split up the SpatialPolygonsDataFrame into the included polygons by coor <- NULL

Hello list, I'm just getting started with using R - I have been trying over the past day or so to work out a method for generating voronoi polygons for PostGIS using SQL. I was able to put together a procedure which works relatively well, but is somewhat inefficient. Someone on the PostGIS list pointed me to the deldir() function in R, for which I can export a text file with x/y coordinates

I have an XML file that has within it the coordinates of some polygons that I would like to extract and use in R. The polygons are nested rather deeply. For example, I found by trial and error that I can extract the coordinates of one of them using functions from the XML package: doc <- xmlInternalTreeParse('doc.kml') docroot <- xmlRoot(doc) pgon <-

Dear all, I'm trying to replicate the panel function behaviour that Deepayan and others discuss here: https://stat.ethz.ch/pipermail/r-help/2007-April/130779.html and Oscar Lamigueiro refers to here and provides sample code for: http://www.r-bloggers.com/confidence-bands-with-lattice-and-r/ where a polygon is generated around lines on panels to display the confidence interval bounds for

Scotland is the second feature in the UK data, so get it and split this one MULTIPOLYGON feature into individual POLYGONS scot = st_cast(the_uk$geometry[2],"POLYGON") # which is the largest polygon? which.max(st_area(scot)) [1] 1 # the first one. ok... plot(scot[[1]]) # mainland # add the rest of the islands for context, in grey, maybe to show they're outside our study area:

Dear list I have two datasets, one with a large number of point coordinates, the other with spatial polygons. I want to measure the shortest distance from each point to the border of all polygons, and then store the results as a column in a data frame. (Or a vector that can be turned into a column in a data frame). The points and polygons represent different types of events, happening at

Dear R-users, I’m currently trying to extract the value of a slot (area) but can’t find out how to do that. >str(overlperc) List of 1 $ :Formal class 'SpatialPolygons' [package "sp"] with 4 slots .. ..@ polygons :List of 1 .. .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots .. .. .. .. ..@ Polygons :List of 1 .. .. .. .. .. ..$

Hello dear members of R-help and R-core mailing list, I am not sure if this request is a "ticket" that should be filled somewhere outside the mailing list. If so, I apologize for not doing and would like to know where I should have filled it. And to the subject matter: I would like to use a command like this: plot(c(1,8), 1:2, type="n") polygon(1:7, c(2,1,2,NA,2,1,2),