[R-sig-Geo] Fill in empty links with nearest neighbors

2011-07-14 Thread Kitty Lee
Hi. I have a shapefile. I guess because of cartography problem, the border of some of the units that are supposed to be physically connected are 'snapped' properly. As a result, when I use poly2nb, there are empty links. I know I can set zero.policy=T in nb2listw. But I'm trying to use skater()

Re: [R-sig-Geo] Long time listener, first time caller

2011-07-14 Thread Lyndon Estes
Hi Robert, On my Mac the plot displays fine with no lines. R version 2.13.1 (2011-07-08) Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) [13] raster_1.8-39 sp_0.9-83 Cheers, Lyndon On Thu, Jul 14, 2011 at 8:54 PM, Michael Sumner wrote: > I see the same on Windows 7 64-bit, with 32-bit

Re: [R-sig-Geo] classify a raster by unique value

2011-07-14 Thread Robert Hijmans
> I've a raster stack with 3 layers. > > I need to classify it assigning one class to each unique combination of > values. Perhaps the below works. It may be very slow if unq has many rows. unq <- unique(as.matrix(data_stack)) fun <- function(x, ...) { r <- rep(NA, nrow(x))

Re: [R-sig-Geo] Long time listener, first time caller

2011-07-14 Thread Michael Sumner
I see the same on Windows 7 64-bit, with 32-bit R and 64bit R. There's no need for the transpose/flip btw: image(volcano) but raster graphics don't product it, so it's the underlying image C code (which is rect under the hood) image(volcano, useRaster = TRUE) plot(0, xlim = c(0, 1), ylim = c(0

Re: [R-sig-Geo] Long time listener, first time caller

2011-07-14 Thread Adam Sparks
On my Linux system with 2.13.1 the R graphical display shows white lines. If I generate/export a PNG using the png() function there are no white lines. > sessionInfo() R version 2.13.1 (2011-07-08) Platform: x86_64-unknown-linux-gnu (64-bit) On Fri, Jul 15, 2011 at 8:40 AM, Robert Hijmans wro

Re: [R-sig-Geo] Long time listener, first time caller

2011-07-14 Thread Robert Hijmans
> Lines appear to be a GUI issue, went back to 2.13.0 and plot of raster looks fine. When plotting a Raster* object, raster calls the 'image' function. spplot uses 'levelplot'. Both are affected by these white lines in R 2.13.1 on Windows, I think (anyone seeing this on Mac or Linux? or not seein

Re: [R-sig-Geo] plotting geographic coordinates on a map

2011-07-14 Thread Robert Hijmans
> I simply would like to plot these on a map of the US. I've tried all sorts of approaches, including >> map('usa', project='albers', par=c(39,34)) >> points(x=data$LONGITUDE, y=data$LATITUTUDE, col='red') You are plotting longitude/latitude data on an map with Albers projection. Won't give you

Re: [R-sig-Geo] plot3D raster coordinates

2011-07-14 Thread Robert Hijmans
> (plot3D is in fact in the rasterVis package now) > > There's an adjustment factor additional to the zfac parameter. It is > impossible to override it. The solution I've found is to comment these > lines > and run the function under another name (these lines are present 2 times). > >#

Re: [R-sig-Geo] Long time listener, first time caller

2011-07-14 Thread dman
Ah-ha, NA's are the cells outside my mask. Lines appear to be a GUI issue, went back to 2.13.0 and plot of raster looks fine. D -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Long-time-listener-first-time-caller-tp6585078p6585234.html Sent from the R-sig-geo mailing

Re: [R-sig-Geo] Long time listener, first time caller

2011-07-14 Thread dman
@SM - I am in fact using 2.13.1 and will look into the image() call. The code should be attached, here it is again http://r-sig-geo.2731867.n2.nabble.com/file/n6585138/dist_1km_code.jpg -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Long-time-listener-first-time-caller

Re: [R-sig-Geo] Long time listener, first time caller

2011-07-14 Thread steven mosher
Hi, please post code to replicate the problem. also include sessionInfo(). how are you plotting the raster, are you using 2.13.1? ( there was a bug report that looks somewhat like yours on that today using the image() call. More details and we can help. Thanks Robert Hijmans University of Cal

[R-sig-Geo] Long time listener, first time caller

2011-07-14 Thread David Kennedy
Hello there, I'm using "raster" to great effect in a little something I'm building, but having a time trying to figure out this current issue. My intent is to build a "distance raster" from a center point out to 1km with cell size of 30x30. I had a few runs yesterday where this worked, but toda

Re: [R-sig-Geo] predict() LDL error

2011-07-14 Thread Edzer Pebesma
The error message is an indication that you're predicting (kriging) with a singular covariance matrix. The most frequent reasons for this are duplicate observations [1] or perfect correlation in the coregionalization model [2]. For [1], see ?remove.duplicates, for [2] see argument correct.diagonal

[R-sig-Geo] spTransform to convert projected coords into geog coord

2011-07-14 Thread Robin W Hunnewell
Re: this previous post, sorry my fault... fn applies to Spatial* objects, not "ppp" objects - my apologies.Date: 14 July 2011 2:50:42 PMTo: r-sig-geo@r-project.orgSubject: [R-sig-Geo] spTransform to convert projected coords into geog coordHello, I get an error when calling function 'spTransform' --

[R-sig-Geo] spTransform to convert projected coords into geog coord

2011-07-14 Thread Robin W Hunnewell
Hello, I get an error when calling function 'spTransform' -- from package {rgdal}. I've used spTransform recently w/ no problem, to convert coordinates from one CRS to coords of another. Has this function been deprecated, or is there some other method?  My code below. I'll lean in for a dope slap i

Re: [R-sig-Geo] Download MODIS images

2011-07-14 Thread ari
Hi, If I want to use ocean data (for example- Sea Surface Temperature, Chlorophyll-a Concententration etc.) or in others words download MOD 28 and MOD 21 data using the same script for a paticular block, which ftp server should I use? Thanks, A -- View this message in context: http://r-sig-geo.

[R-sig-Geo] predict() LDL error

2011-07-14 Thread Matevž Pavlič
Hi all, i am trying to use categorical kriging for prediction of soil classes in an area. I am pretty new to using cat. kriging in R soI probably am making some beginers mistakes. Attached is a reproducable sample with sample data. When i try to use predict() i geta n error- I am pretty

Re: [R-sig-Geo] plotting geographic coordinates on a map

2011-07-14 Thread Ashton Shortridge
On 2011-07-14, atlasrrg, wrote: > Hello, > > I am pretty sure that I am making a novice mistake here as I've only just > begun exploring spatial data in R. I have a data.csv file that contains > geographic coordinates as follows: > > SEQUENCE,longitude,latitude,altitude > AF196835_Pc_40.77_73.96_

[R-sig-Geo] plotting geographic coordinates on a map

2011-07-14 Thread atlasrrg
Hello, I am pretty sure that I am making a novice mistake here as I've only just begun exploring spatial data in R. I have a data.csv file that contains geographic coordinates as follows: SEQUENCE,longitude,latitude,altitude AF196835_Pc_40.77_73.96_1999.71,-73.96,40.77,18 AF202541_Hs_42.76_75.82