Re: [R-sig-Geo] Mapping exponentials from GWR

2011-08-15 Thread Roger Bivand
On Mon, 15 Aug 2011, Michelle Convey wrote: Hi there, I am working on a project with the spgwr package.  I am able to map the coefficients from my GWR analysis and am able to get summary data on the exponentials on the same data, but am not sure how to map them.  I'm probably missing somethi

Re: [R-sig-Geo] anisotropy modeling [SEC=UNCLASSIFIED]

2011-08-15 Thread Jin.Li
Hi All, Following Edzer's examples, I have tried to test the effects of the specification of alpha in variogram on the results of fit.variogram as below: > v = variogram(zinc~1,meuse) > v0 = variogram(zinc~1,meuse, alpha = 0, tol.hor = 45) > v90 = variogram(zinc~1,meuse, alpha = 90, tol.

Re: [R-sig-Geo] geo meetup at UseR! 2011

2011-08-15 Thread Pierre Roudier
Hi Barry, I'll be happy to join :) I'm sure there's some more spatial useRs out there in Warwick at the moment, Cheers, Pierre 2011/8/14 Barry Rowlingson : >  If any R-sig-geo readers are going to UseR! at Warwick U next week, > there's a plan to meet up on Tuesday evening for a chat about all

Re: [R-sig-Geo] Changing the NA value

2011-08-15 Thread Brian Oney
Hi Robert, thanks for the full report. The suggestion to write a raster file with a certain NA values was all I needed. Thinking about it now, other users would probably also like to be able to determine the NA-value when writing files, so I wouldn't expect a common need. I will try to stick to

[R-sig-Geo] Segmentation of SpatialLines at Intersection

2011-08-15 Thread Cam Plouffe
Hi everybody, I have a SpatialLines object that I created by using the gBoundary() function on a SpatialPolygons object. I've attached a plot ("SpatialLines Segmentation.png") of the SpatialLines object to this email. What I would like to do, is segment these SpatialLines wherever there is an in

Re: [R-sig-Geo] Changing the NA value

2011-08-15 Thread Robert J. Hijmans
Brian, If you simply want to change the NA flag (i.e. the value used used in a _file_ to represent NA), you should use the NAvalue function NAvalue(rast) <- -9995 But note that this is only to be used when there is an error in the file metadata, or other reason why values that should be NA are n

Re: [R-sig-Geo] Error when projecting polygons into Azimuthal Equidistant Projection

2011-08-15 Thread Roger Bivand
On Mon, 15 Aug 2011, bart wrote: Hi All First a short introduction to what i'm trying to do, I'm trying to find the length of a line going through a point intersecting with the coast lines. My thought was to find the right polygon and then re-project it to an Azimuthal Equidistant Projection

[R-sig-Geo] gstat: Specify user defined variogram or specify covariance matrices

2011-08-15 Thread Jonathan Kennel
Hello, As a learning experience I have made a kriger in R and have been exploring different variogram models. I really like the speed and simplicity of entering data to the gstat package which solves much faster than my toy kriger. I have taken a brief look at the gstat C source code where the v

[R-sig-Geo] Error when projecting polygons into Azimuthal Equidistant Projection

2011-08-15 Thread bart
Hi All First a short introduction to what i'm trying to do, I'm trying to find the length of a line going through a point intersecting with the coast lines. My thought was to find the right polygon and then re-project it to an Azimuthal Equidistant Projection on which an intersection with the

[R-sig-Geo] Mapping exponentials from GWR

2011-08-15 Thread Michelle Convey
Hi there, I am working on a project with the spgwr package.  I am able to map the coefficients from my GWR analysis and am able to get summary data on the exponentials on the same data, but am not sure how to map them.  I'm probably missing something quite obvious as I am fairly new to this, bu

Re: [R-sig-Geo] Reading a WMS layer as a sp object

2011-08-15 Thread Manuel Schneider
Following Michael's advice I'm trying to execute the first example on http://www.gdal.org/frmt_wms.html and run gdal_translate <- "C:/Progs/FWTools2.4.7/bin/gdal_translate.exe" download.file("http://www.gdal.org/frmt_wms_onearth_global_mosaic.xml";, destfile="C:/temp/test.xml") wms_xml <- "C:/temp

Re: [R-sig-Geo] A one to many merge involving a spatial data frame object

2011-08-15 Thread Paul Hiemstra
On 08/09/2011 08:23 PM, Dan Putler wrote: > All, > > I'm working with the 2010 TIGER/Line shapefile road data, and I'm > dealing with alternative road names for the same road segment. There > are two relevant files in the TIGER collection, the first is the road > edge shapefile, which has a unique