[R-sig-Geo] sampling a raster using sample points

2011-01-14 Thread Stratford, Jeffrey
Hi everyone, Does anyone have any suggestions on sampling a raster (NLCD land cover map, 30 m res) using circular buffers centered on 100+ sample points? The NLCD map has classifications such as high density urban, medium density urban, crop, etc. The goal is to produce a data set with the perce

[R-sig-Geo] difficulty converting DDD to UTM (PBSmapping/convUL)

2011-01-23 Thread Stratford, Jeffrey
Hi everyone, I am attempting to convert points in DDD to UTM (zone 18) using PBSmapping and I'm getting an error message. The code and output: library(PBSmapping) dd_sosp <- read.csv("g:\\sosp\\2010\\sosp2010coord.csv", header=T) attr(dd_sosp, "zone") <- 18 str(dd_sosp) 'data.frame': 109 obs

[R-sig-Geo] difficulty converting DDD to UTM): closer

2011-01-23 Thread Stratford, Jeffrey
Hi everyone, I've gotten closer to a solution but there's still something wonky going on. Here's the new code I'm using dd_sosp <- read.csv("g:\\sosp\\2010\\sosp2010coord.csv", header=T) names(dd_sosp) <- c("X","Y") attr(dd_sosp, "projection") <- c("LL") sosp.utm <- convUL(dd_sosp)

[R-sig-Geo] creating an adjacency matrix in spdep

2011-01-30 Thread Stratford, Jeffrey
Hi everyone, Is it possible to create an adjacency matrix required for WinBUGS/GeoBUGS in spdep from a csv file containing UTM coordinates? I've gone through the spdep documentation and it looks like the equivalent there is col.gal.nb? Is that right? I see where that object is used but n

Re: [R-sig-Geo] creating an adjacency matrix in spdep

2011-02-01 Thread Stratford, Jeffrey
help! I'm going through your spatial book and it's really good - lots of tips! -Original Message- From: Roger Bivand [mailto:roger.biv...@nhh.no] Sent: Monday, January 31, 2011 3:30 AM To: Stratford, Jeffrey Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] creating an adjacency ma

[R-sig-Geo] Error in anyDuplicated.default(coords) message with trinb

2011-07-07 Thread Stratford, Jeffrey
Hi everyone, I'm attempting to do some spatial analyses with a relatively simple data set and I'm running into some issues. Here's my code: sosp09 <- read.csv("f:\\sosp\\2009\\sosp09.means.csv", header=T) sp.point <- cbind(sosp09$x, sosp09$y) colnames(sp.point) <- c("LONG","LAT") # Pro