similar to: extracting coordinates from SpatialPolygonsDataFrame using slot functions

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

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`

Hello, I am using the following code to create ppp files from csv data and map shape files, but I am getting some errors which I have been unable to fix by searching them online: library(spatstat) library(maps) library(maptools) NYC2<-readShapePoly("nybb.shp") # this is a map of the NYC boroughs without waterways and no census tract divisions (but it does include lines separating

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 <- readShapeSpatial("200001441.shp") y <- as(x, "SpatialPolygons") z <- as.owin(y) Given z, I want to identify each

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

Hello All, I'm running Samba 3.0.24 as PDC and file server on Gentoo Linux, AMD64, 2.6.19 kernel, 100MBit/s network, and experience quite slow file transfers from Samba to WinXP SP2 clients: the speed is varying, but is about 1-2Mb/s at best. I spent quite some time investigating the issue, here are the intermediate results: 1) This happens only when transferring _from_ _Samba_ to

Ok guys, as requested, I will add more info so that you understand why a simple vector operation is not possible. It's not easy to explain in few words but let's see. I 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

My sincerest apologies, as this is a very elementary problem, but I have searched through archives and FAQs and on the web, and I am at the end of my own resources. I need to do analysis of a spatial point process occuring with four counties, using spatstat. I've had no trouble importing the shapefiles, but I can't seem to figure out what I need to do in order to use the boundary as an

So I wanted to export a contour3d object into VTK and it worked well enough, except the top of the ball (the contour) was missing. Example below, or gist here: https://gist.github.com/muschellij2/5126544 rm(list=ls()) require(misc3d) f <- function(x, y, z) x^2+y^2+z^2 x <- seq(-2,2,len=20) ball <- contour3d(f,4,x,x,x, draw=FALSE) drawScene.rgl(ball) ### get the triangles allids

Dear Community I am at my wits end and seek advice. My wish is to plot coordinates (x,y in WGS84_UMTS for the ones interested) of sampling points. This I can do by the standard spatstat prodcedure via owin. I then try to add an image, which is a map/satellite photo in the background. Firstly Problem: With the current code, I can not get the edges of the image to match the boarders of the plot

Dear users, I am trying to export polygons from Arcmap into Spatstat to run some simulations using functions available in Spatstat package. One particular area to be exported is formed by a number of polygons defining the external boundaries of the area (as a groups of islands) and a number of polygons inside the previous ones, as ?holes? not to be considered as part of the area. I have

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)),

I am building a mask with area.urb.sp.W <- as(area.urb.sp, "owin") and got the message Erro em owin(poly = opls) : Polygon data contain duplicated vertices, self-intersection and overlaps between polygons how can I solve this problem in R? many thanks IVAN [[alternative HTML version deleted]]

I have a simple problem: I need to load a ERSI shapefile of US states and check whether or not a set of points are within the boundary of these states. I have the shapefile, I have the coordinates but I'm having a great deal of difficulty bringing the two together. The problem is the various GIS packages for R do not play well with each other. sp, shapefiles, maptools, etc all use different

I would like to use the output from the map function in the package maps for use in, say, the spatstat package. I don't quite understand the coordinates for the border of the state: Example: library(maps) iowa<-map("region","iowa) x<-iowa$x y<-iowa$y There are NA's and duplicated coordinates. What would I need to do to use this in the spatstat owin command?

Good day Im working with some coordinates, and want to create a PPP object, I found that error: > Datos=read.table("puntos_texto.txt",dec=".",sep="\t",header=T) > summary(Datos) id y x Min. : 1.0 Min. :1013581 Min. :1177842 1st Qu.: 821.2 1st Qu.:1014442 1st Qu.:1179658 Median

Hello everbody, I am new to this mailing list and hope to find some help. I'm trying to get into the spatstat package and encountered two problems. First a graphical one: There is an example dataset called "finpines" which has several marks (http://www.oga-lab.net/RGM2/func.php?rd_id=spatstat:finpines) When I pass the given code from the website to R data(finpines)

I've written the function below to simulate the mean 1st through nth nearest neighbor distances for a random spatial pattern using the functions nndist() and runifpoint() from spatsat. It works, but runs relatively slowly - would appreciate suggestions on how to speed up this function. Thanks. --Dale library(spatstat) sim.nth.mdist <- function(nth,nsim) { D <- matrix(ncol=nth,

Dear list, I'd like to use the function smooth.map{maps} to smooth out aggregated spatial data. The first argument to this function must be a map object. I have a shape file which I read with readShapeSpatial{maptools} and I got an object of class "SpatialPolygonsDataFrame". Is it possible to create a map object from a SpatialPolygonsDataFrame object? Alternatively, is there any

I have obtained shapefiles for Indian states from here: http://www.maptell.com/index.php?option=com_remository&Itemid=159&func=fileinfo&filecatid=115&parent=category Problem: I want to extract centroid coordinates for each State, but there is some coding problem with the shapefiles that prevents this. #Code: #After extracting the shapefiles from the india_state.zip file, then: