[R-sig-Geo] more doubts with spTransform()

2009-07-14 Thread Agustin Lobo
Inverse transform now: anyellaWGS84 - readOGR(dsn=., layer=parcellesAnyellaWGS84) OGR data source with driver: ESRI Shapefile Source: ., layer: parcellesAnyellaWGS84 with 6 rows and 7 columns Feature type: wkbPoint with 2 dimensions proj4string(anyellaWGS84) [1] +proj=utm +zone=31

Re: [R-sig-Geo] more doubts with spTransform()

2009-07-14 Thread Roger Bivand
On Tue, 14 Jul 2009, Agustin Lobo wrote: Inverse transform now: anyellaWGS84 - readOGR(dsn=., layer=parcellesAnyellaWGS84) OGR data source with driver: ESRI Shapefile Source: ., layer: parcellesAnyellaWGS84 with 6 rows and 7 columns Feature type: wkbPoint with 2 dimensions

Re: [R-sig-Geo] more doubts with spTransform()

2009-07-14 Thread Tomislav Hengl
Hi Agustin, You should really work with some referent point for which you know both local and geographic (WGS84) coordinates. Here is an illustration: http://spatial-analyst.net/wiki/index.php?title=MGI_/_Balkans_coordinate_systems#Validation_of_CRS_p arameters Note that I use Google Earth

[R-sig-Geo] getting subfolders in a directory on FTP (MODIS)

2009-07-14 Thread alexandre villers
Good morning, I know this is not exactly a question for R SIG GEO, but as it may interest beginners like me who work with MODIS data... (and I'm hoping I'm not the only beginner on the list !) I adapted the script of Tomislav Hengl (again thank you for your wiki page!

Re: [R-sig-Geo] kml generation question

2009-07-14 Thread Matt Oliver
Turns out the best idea was Barry's..Simply co-opt some existing kml code and make it do what you want with your image http://www.ceoe.udel.edu/cms/moliver/EDA_20090619.170.1758.n18.kml On Mon, Jul 13, 2009 at 12:12 PM, Roger Bivand roger.biv...@nhh.no wrote: On Mon, 13 Jul 2009, Barry

Re: [R-sig-Geo] Convert SpatialGridDataFrame to gridTopology object

2009-07-14 Thread Roger Bivand
On Tue, 14 Jul 2009, Paulo E. Cardoso wrote: How can I convert one to another? library(sp) data(meuse.grid) coordinates(meuse.grid) - c(x, y) gridded(meuse.grid) - TRUE fullgrid(meuse.grid) - TRUE class(meuse.grid) grd - slot(meuse.grid, grid) class(grd) Paulo E. Cardoso

Re: [R-sig-Geo] getting subfolders in a directory on FTP (MODIS)

2009-07-14 Thread Barry Rowlingson
On Tue, Jul 14, 2009 at 2:12 PM, Tomislav Henglhe...@spatial-analyst.net wrote: To tell you honestly, I could not figure out how to fetch names of sub-folders from an ftp (maybe it is simply not possible?), Apparently the FTP spec doesn't specify how to return a list of directories, which,

Re: [R-sig-Geo] getting subfolders in a directory on FTP (MODIS)

2009-07-14 Thread Tomislav Hengl
Hey Barry, Your code works fine (at least with the NASA ftps). Thanks! I knew it - everything is possible in R :) T. Hengl http://spatial-analyst.net -Original Message- From: b.rowling...@googlemail.com [mailto:b.rowling...@googlemail.com] On Behalf Of Barry Rowlingson Sent:

Re: [R-sig-Geo] getting subfolders in a directory on FTP (MODIS)

2009-07-14 Thread alexandre villers
Dear Barry, Thank you very much for this ! This will really help ! Best regards Alex Alexandre Villers PhD. Student Team Biodiversity CEBC CNRS 79360 Beauvoir sur Niort __ Information from ESET Mail Security, version of virus signature database 4242 (20090714) __ The message

[R-sig-Geo] split SpatialPolygon after unionSpatialPolygons()

2009-07-14 Thread pecardoso
I've used unionSpatialPolygons to dissolve polygons and get a grid as a spatialpolygon. A simple example: grelha-readGDAL(paste(C:/...,teste.bmp,sep=)) grd - slot(grelha, grid) pol - as.SpatialPolygons.GridTopology(grd) pol - unionSpatialPolygons(pol, SGDF2PCT(grelha)$idx) When the original