My best guess is that you need to add your projection to one of the files in /usr/share/gdal or in the projection definition files shipped with rgdal. Please send it upstream if you do.
THK On 9/11/06, Zev Ross <[EMAIL PROTECTED]> wrote: > Hi All, > > I think this is probably a simple one. I'm reading a shapefile, making a > couple of changes and then sending it back out to shapefile. I am making > no changes to the prj information but for some reason it seems that the > projection name is being lost in the process. This is a non-traditional > projection, one used by many California government agencies. It appears > that the actual projection information is being maintained correctly, > but I'd like to keep the name. > > Zev > > Here's the relevant piece of the code: > > myprj<-readOGR("c:/junk/monitors", layer="monitors_all") > showWKT(proj4string(myprj),file="c/junk/monitors_all.prj") > > here is the prj that I'm reading in > > PROJCS["NAD_1983_California_Teale_Albers",GEOGCS["GCS_North_American_1983", > DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]], > PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Albers"], > PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",-4000000.0], > PARAMETER["Central_Meridian",-120.0],PARAMETER["Standard_Parallel_1",34.0], > PARAMETER["Standard_Parallel_2",40.5],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] > > > here is the output prj > > PROJCS["Unknown",GEOGCS["GCS_North_American_1983", > DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137,298.257222101]], > PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Albers"], > PARAMETER["standard_parallel_1",34],PARAMETER["standard_parallel_2",40.5], > PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-120], > PARAMETER["false_easting",0],PARAMETER["false_northing",-4000000],UNIT["Meter",1]] > > > > -- > Zev Ross > ZevRoss Spatial Analysis > 303 Fairmount Ave. > Ithaca, NY 14850 > (607) 277-0004 (phone) > (866) 877-3690 (fax toll-free) > [EMAIL PROTECTED] > www.zevross.com > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Timothy H. Keitt http://www.keittlab.org/ _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo