Re: [R-sig-Geo] Prevent minor discrepancy in pixel resolution using Spatstat

2010-02-17 Thread Roger Bivand
On Wed, 17 Feb 2010, Tyler Dean Rudolph wrote: When I specify eps=25 (m resolution) in a call to distmap... library(spatstat) temp <- data.frame(x=-670049.2, y=752814.6, X2=-669961.9, Y2=752648.2) subwin<-owin(xrange=c(-100, 100) + sort(c(temp$x, temp$X2)), yrange=c(-100, 100) + sort

Re: [R-sig-Geo] Rpy with ArcMap or Spatial Data

2010-02-17 Thread Roger Bivand
On Wed, 17 Feb 2010, Alexander Brenning wrote: Hi, there is an ArcGIS toolbox written by Mark Janikas from ESRI that links ArcGIS/ArcMap/ArcCatalog with R through Python. This is more like a proof-of-principle, it allows to access a point clustering method from R from within ArcGIS. Don't as

Re: [R-sig-Geo] Rpy with ArcMap or Spatial Data

2010-02-17 Thread Alexander Brenning
Hi, there is an ArcGIS toolbox written by Mark Janikas from ESRI that links ArcGIS/ArcMap/ArcCatalog with R through Python. This is more like a proof-of-principle, it allows to access a point clustering method from R from within ArcGIS. Don't ask me for the details, but at least one of the tw

[R-sig-Geo] spatial regression for large data sets

2010-02-17 Thread Seth J Myers
Hi everyone, I'm trying to model a binary response using logistic regression for a large data set with spatial autocorrelation issues. The mixed models in SAS and R that can include spatially correlated errors cannot handle the large NxN matrix needed for their methods. Past around 700 meters

[R-sig-Geo] Prevent minor discrepancy in pixel resolution using Spatstat

2010-02-17 Thread Tyler Dean Rudolph
When I specify eps=25 (m resolution) in a call to distmap... > library(spatstat) > temp <- data.frame(x=-670049.2, y=752814.6, X2=-669961.9, Y2=752648.2) >subwin<-owin(xrange=c(-100, 100) + sort(c(temp$x, temp$X2)), yrange=c(-100, 100) + sort(c(temp$y, temp$Y2))) >trajpsp<-psp(te

Re: [R-sig-Geo] Merge polygons

2010-02-17 Thread Robert J. Hijmans
Hi Tyler, Your example would have been more easily reproducible if it had started with something like: temp = list(x=-300358, y=748992.6, X2=-300450.8, Y2=748764.8) library(spatstat) ... library(maptools) union = unionSpatialPolygons(polybuffs, ID=rep(1, times=length(polybu...@polygons)) ) plo

[R-sig-Geo] Merge polygons

2010-02-17 Thread Tyler Dean Rudolph
I've been attempting for many moons to develop a function which will allow me to place a buffer of adjustable radius around an animal's movement trajectory with which I may extract underlying pixel values from a raster of environmental data. I would like to be able to extract the habitat values us

Re: [R-sig-Geo] Calculate location given observation distance and angle from transect line

2010-02-17 Thread Robert J. Hijmans
Corrie: I think you can use the function 'destPoint' in 'geosphere', It computes (for longitude/latitude coordinate data) the coordinates of a destination point given a start point, initial bearing, and distance. Robert On Wed, Feb 17, 2010 at 8:00 AM, Corrie Curtice wrote: > Hello, > > I've tu

Re: [R-sig-Geo] compute buffer from point shapefile to have shapefile

2010-02-17 Thread Roger Bivand
On Wed, 17 Feb 2010, Tyler Dean Rudolph wrote: Here is another way to compute buffered polygons around points: ## One record (point) at a time (buffer of 25m radius) temp xy 1 -300358 748992.6 library(spatstat) discbuff<-disc(radius=25, centre=c(temp$x, temp$y)) Note that fr

[R-sig-Geo] Calculate location given observation distance and angle from transect line

2010-02-17 Thread Corrie Curtice
Hello, I've turned up little in searching for a package that already has this function. Dspat is close but I don't think it does quite what I want. I have observation data gathered from a ship transect doing a zig-zag pattern. I have the current location of the ship, the ship's bearing, the ang

Re: [R-sig-Geo] compute buffer from point shapefile to have shapefile

2010-02-17 Thread Tyler Dean Rudolph
Here is another way to compute buffered polygons around points: ## One record (point) at a time (buffer of 25m radius) > temp xy 1 -300358 748992.6 library(spatstat) discbuff<-disc(radius=25, centre=c(temp$x, temp$y)) library(sp) polybuff<-SpatialPolygons(list(Polygons(list(Polyg