Re: [R-sig-Geo] Res: converting shapefile/polygon to psp

2008-02-04 Thread Roger Bivand
On Mon, 4 Feb 2008, Milton Cezar Ribeiro wrote: Dear Roger, I am not sure if I really need to convert my polygons to psp format. In fact I have a shapefile with polygons and a pair of points (x0,y0,x1,y1) and what I need is: (1) to checkout if the lines defined by (x0,y0,x1,y1) cross or n

[R-sig-Geo] Res: converting shapefile/polygon to psp

2008-02-04 Thread Milton Cezar Ribeiro
Dear Roger, I am not sure if I really need to convert my polygons to psp format. In fact I have a shapefile with polygons and a pair of points (x0,y0,x1,y1) and what I need is: (1) to checkout if the lines defined by (x0,y0,x1,y1) cross or not any polygon and (2) to identify which polygons

Re: [R-sig-Geo] R to Python - rpy - definition of window pattern, related with 'spatstat library - help needed

2008-02-04 Thread G. Allegri
I think it's better to wrap R scripting and void passing unuseful objects back and forth between Python and R. In your script it doesn't seem that you need direct interaction between Python and R (at least in this part of your script). I usually adopt something like this: r("""CairoPNG(file="/var/

Re: [R-sig-Geo] converting shapefile/polygon to psp

2008-02-04 Thread Roger Bivand
On Mon, 4 Feb 2008, Milton Cezar Ribeiro wrote: > Dear all, > > I have some shapefiles which contents are polygons and I would like to > convert it to psp (spatstat) format, and after run crossdist.psp. I > tryed several times without success. Any idea? Look at the coercion function in maptools

[R-sig-Geo] converting shapefile/polygon to psp

2008-02-04 Thread Milton Cezar Ribeiro
Dear all, I have some shapefiles which contents are polygons and I would like to convert it to psp (spatstat) format, and after run crossdist.psp. I tryed several times without success. Any idea? Miltinho para armazenamento! [[alternative HTML version deleted]] ___

Re: [R-sig-Geo] FYI: Merging GIS and statistics --- RSAGA

2008-02-04 Thread Roger Bivand
On Mon, 4 Feb 2008, Alí Santacruz wrote: Hi, I am trying to use the Grid Calculator in SAGA, through RSAGA but I have problems with the parameter FORMUL, I can´t get the function "a/b" works, the functions "a+b", or "a-b", or "a*b" works, but "a/b" does not, Any help is appreciated, see

Re: [R-sig-Geo] FYI: Merging GIS and statistics --- RSAGA

2008-02-04 Thread Alí Santacruz
Hi, I am trying to use the Grid Calculator in SAGA, through RSAGA but I have problems with the parameter FORMUL, I can´t get the function "a/b" works, the functions "a+b", or "a-b", or "a*b" works, but "a/b" does not, Any help is appreciated, see the code below for the errors reported

Re: [R-sig-Geo] Res: Subsetting shapefile type PolyLine

2008-02-04 Thread Roger Bivand
On Mon, 4 Feb 2008, Rafael Duarte wrote: Obrigado Milton. But: teste2[teste2$ELEVATION == 100] gives me an empty list. And I have data of Elevation = 100. I also tried teste2[teste2$ELEVATION >= 100] teste2[teste2$ELEVATION == "100"] But always an empty list and no plot. Very close, with

Re: [R-sig-Geo] Res: Subsetting shapefile type PolyLine

2008-02-04 Thread Rafael Duarte
Obrigado Milton. But: teste2[teste2$ELEVATION == 100] gives me an empty list. And I have data of Elevation = 100. I also tried teste2[teste2$ELEVATION >= 100] teste2[teste2$ELEVATION == "100"] But always an empty list and no plot. Milton Cezar Ribeiro wrote: > Hi Rafael > i think th

Re: [R-sig-Geo] R to Python - rpy - definition of window pattern, related with 'spatstat library - help needed

2008-02-04 Thread Roger Bivand
On Mon, 4 Feb 2008, volkan kepoglu wrote: > I am new in R and Python. I write a couple of lines. The code does the > followings; Read point shp file, Compute a kernel smoothed intensity > function from a point pattern; density function, and convert pixel > image to spatialgriddataframe and export

[R-sig-Geo] Res: Subsetting shapefile type PolyLine

2008-02-04 Thread Milton Cezar Ribeiro
Hi Rafael i think that it could help you teste2.elev100<-teste2[teste2$ELEVATION=100] plot(teste2,axes=T) plot(teste2.elev100,col=2,add=T) Boa sorte. Miltinho Brazil - Mensagem original De: Rafael Duarte <[EMAIL PROTECTED]> Para: r-sig-geo@stat.math.ethz.ch Enviadas: Segunda-feira, 4

[R-sig-Geo] R to Python - rpy - definition of window pattern, related with 'spatstat library - help needed

2008-02-04 Thread volkan kepoglu
I am new in R and Python. I write a couple of lines. The code does the followings; Read point shp file, Compute a kernel smoothed intensity function from a point pattern; density function, and convert pixel image to spatialgriddataframe and export spatialgriddataframe to tif raster format. the cod

[R-sig-Geo] block kriging of skewed distributions

2008-02-04 Thread Andrew Bevan
Dear All Apologies if the following betrays any naive assumptions but I would be grateful for the advice nonetheless. I am looking to do some block kriging on a pottery counts from an archaeological survey covering the whole surface of an island (ca. 20sq.km, for details see www.ucl.ac.uk/as

[R-sig-Geo] Subsetting shapefile type PolyLine

2008-02-04 Thread Rafael Duarte
Dear all, I am very new to SIG. This is probably a very simple question but I could not find an answer searching the web and the R-SIG emails archive. I have an ESRI shapefile with a coast and bathymetric lines. In the dbf file I have a column “ELEVATION” that gives this information. I import