[R-sig-Geo] read prj file to CRS

[Roger Bivand]

On Thu, 17 May 2007, Patrick Giraudoux wrote:

> Dear all,
> 
> Definitely, I am not that clever with projection systems: I would like 
> to read a *.prj file (part of a set of ESRI shapefiles) and get it in a 
> CRS compatible format (handled by the  PROJ.4 projection system) and 
> possibly added to sp objects where needed. What would would be the 
> simplest way to do that? I am messing with rgdal and sp without result...
> 

Yes, good idea, given that the reverse is provided in showWKT() in rgdal. 
As things are, you'll need to read it as part of a shapefile - readOGR() 
will read it when given a shapefile to read, or readGDAL() with an Arc 
Ascii grid, driver AAIGrid. I'll look at putting a simple *.prj reader 
into rgdal, but for now giving an arbitrary shapefile or Ascii grid the 
same name should do it.

Hope this helps,

Roger

> Cheers,
> 
> Patrick
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> 

-- 
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: Roger.Bivand at nhh.no