Skip to content
Prev 4889 / 29559 Next

e00 to SpatialPolygons: how to set the correct projection

On 26.01.2009, at 12:06, Paul Hiemstra wrote:

        
On 26.01.2009, at 12:06, Roger Bivand wrote:
Thanks a lot for both hints! I got it ;)

I did this:

 > e00proj <- CRS(proj4string(readOGR(e00File, "ARC")))
 > e00proj
CRS arguments:
  +proj=lcc +lat_1=25 +lat_2=47 +lat_0=10 +lon_0=110 +x_0=0 +y_0=0  
+ellps=clrk66 +units=m +no_defs

 > NMG <- pal2SpatialPolygons(arcs, pal, IDs=IDs, proj4string = e00proj)
 > NMG <- spTransform(NMG, CRS("+proj=longlat +ellps=WGS84  
+datum=WGS84 +no_defs +towgs84=0,0,0"))


--Hans

**********************************************************
Hans-Joerg Bibiko
Max Planck Institute for Evolutionary Anthropology
Department of Linguistics
Deutscher Platz 6     phone:   +49 (0) 341 3550 341
D-04103 Leipzig       fax:     +49 (0) 341 3550 333
Germany               e-mail:  bibiko at eva.mpg.de