search for: owin

...ge, 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 polygon with a single constant value. This is like adding marks to the SpatialPolygons file. Then I want to convert these polygons into an image, such that the value of each pixel corresponds to the value associated to the polygon where the pixel lies. I have...

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` mass.win <-...

...rst: > > 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 polygon with a single constant value. This is like adding marks to the SpatialPolygons file. Then I want to convert these polygons into an image, such that the value of each pixel corresponds to the value associated to the polygon where the pixel lies....

...th the same results. 8) Samba logs looks normal to me, nothing special. Points (3) and (5), and also the article about SO_SNDBUF tuning gave me an idea that there might be smth wrong with ACKs sent from WinXP. I've tried to search for problems with ACK delays in WinXP and found the following articles: http://support.microsoft.com/kb/328890 http://support.microsoft.com/kb/321098 Based on them, I've done the following changes in WinXP registry: - set TcpAckFrequency to 1 to switch off delayed ACKs - set TcpDelAckTicks to 0 to disable delayed ACKs timer These changes did not...

...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 owin. Why I use the Map2poly command, it appears to work, but not when I attempt define the window using the shapefile; it also gives me errors when I try to define a ppp.object. I am using the latest release of Raqua on G4. My process looks like this: #begin code > library(maptools) > Shapefile...

...ge. 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 managed to export one of these areas using owin(poly=list()), including in the list the number of external boundaries and holes, but the number of polygons that constitute each area is large and I was wondering if this could be easily done using owin(mask=my_area), using only one single matrix file "my_area" with TRUE, FALSE values def...

In converting some older code: ## builds an OWIN to generate a set of points on a grid within the OWIN object: bdry <- read.shape( "../../data/boundary.shp" ) outline.verts <- bdry$Shapes[[2]]$verts[1:bdry$Shapes[[2]]$Pstart[2],] outline.order <- nrow( outline.verts ) : 2 bdry.poly <- vector( 2, mode="list" ) bdry....

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 itself (not the entire plot window, solely the coordinate system) The size of the image I have (from GPS wor...

Dear R list, I'm trying to create an object of class "owin" (observation window) in the package spatstat from GIS mapping data. Here's an example of my problem. Everything goes well until the last line of code. I get the error message shown at the bottom: library(spatstat) library(sp) library(maptools) con <- url("http://gadm.org/data/rda/...

Preamble: ========= In spatial point pattern analysis an important component of a point pattern object is the *observation window*. In the spatstat package we have the function owin() to *construct* such a window. It creates an object of class "owin" which can happily exist separately from any point pattern object (object of class "ppp"). However we have at the moment no convenient function to *extract* the observation component of a ppp object, nor ha...

...ge, 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 polygon with a single constant value. This is like adding marks to the SpatialPolygons file. Then I want to convert these polygons into an image, such that the value of each pixel corresponds to the value associated to the polygon where the pixel lies. I have...

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

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)), unitname=c("metres","metres"), marks=dbh) I managed to create the observation window with this com...

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

...pp(Datos$x,Datos$y,c(min(Datos$x)"1177842",max(Datos$x)"1180213"),c(m in(Datos$y)"1013581",max(Datos$y)"1015575")) Error: inesperado string constante en "danta=ppp(Datos$x,Datos$y,c(min(Datos$x)"1177842"" danta=ppp(Datos$x,Datos$y,window=owin(c(min(Datos$x)"1177842",max(Datos$x)"1 180213"),c(min(Datos$y)"1013581",max(Datos$y)"1015575"))) Error: inesperado string constante en "danta=ppp(Datos$x,Datos$y,window=owin(c(min(Datos$x)"1177842"" Thank you so much [[alternative...

Hi! I have no idea if this is maybe an easy task utilizing R since I read there is geographical map data in some package: I have a huge number of geographical points with their coordinates in Germany. Now I want to determine for each point in which "Bundesland" = state it is located. Can anybody tell me if this is doable relatively easy in R and if so give me some hints or links

...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? owin(poly= ?) Thank you. -- Laura ************************************************ Laura Chihara Professor of Mathematics 507-222-4065 (office) Dept of Mathematics 507-222-4312 (fax) Carleton College 1 North College Street Northfield MN 55057

...-list. It is not my role to provide private consulting for you. I am therefore cc-ing the list in this email. (2) It is still not completely clear what you really want; the figure that you attached appears to be a disc with 4 diameters superimposed. So you might be after a single (circular) owin object and a line segment pattern consisting of the 4 diameters. Or you might be after *eight* owin objects, each being one the eight disc-segments into which the diameters divide the disc. I shall assume the latter. To start with define a function, say "wedge": wedge <- functio...

...hich.marks" ist kein Grafikparameter Something like "which.marks" is not a graphic parameter; and the plots for height and diameter show now differences. Furthermore, I want to create a ppp with several marks, but I did not figure out how this works. Trying X <- as.ppp(mydata, owin(c(174, 178), c(29, 33))) just gives the error Error in as.ppp(mydata, owin(c(174, 178), c(29, 33))) : X must be a two-column or three-column data frame The data set looks something like Date X Y Mar1 Mar2 Mar3 1.1. 4 3 50 6 A 2.1. 2 1 40 9...

....$ hole: logi FALSE .. ..$ units :List of 3 .. .. ..$ singular : chr "unit" .. .. ..$ plural : chr "units" .. .. ..$ multiplier: num 1 .. .. ..- attr(*, "class")= chr "units" .. ..- attr(*, "class")= chr "owin" ..$ : logi NA > area.owin(a[1,]$polygons[[1]]) [1] 7850.828 thanks, allie