[QGIS-trac] Re: [Quantum GIS] #1537: Units unknown in GRASS mapset creation for EPSG 26745

[Quantum GIS]

#1537: Units unknown in GRASS mapset creation for EPSG 26745
-----------------------------------------------------+----------------------
        Reporter:  jctull                            |         Owner:  nobody                                        
            Type:  bug                               |        Status:  new                                           
        Priority:  major: does not work as expected  |     Milestone:  Version 1.0.2                                 
       Component:  GRASS                             |       Version:  1.0.0                                         
      Resolution:                                    |      Keywords:  grass mapset creation epsg:26745 units unknown
Platform_version:                                    |      Platform:  All                                           
        Must_fix:  No                                |   Status_info:  0                                             
-----------------------------------------------------+----------------------
Comment (by wonder):

 Indeed it looks to be OGR problem as Markus suggested. From my tests - the
 proj4 string for this projection (epsg 26745) is:

 {{{
 +proj=lcc +lat_1=35.46666666666667 +lat_2=34.03333333333333 +lat_0=33.5
 +lon_0=-118 +x_0=609601.2192024384 +y_0=0 +ellps=clrk66 +datum=NAD27
 +to_meter=0.3048006096012192 +no_defs
 }}}

 When you import it to OGR and export as WKT, you'll get:

 {{{
 PROJCS["unnamed",
   GEOGCS["NAD27",
     DATUM["North_American_Datum_1927",
       SPHEROID["Clarke
 1866",6378206.4,294.978698213898,AUTHORITY["EPSG","7008"]],
       TOWGS84[-3,142,183,0,0,0,0],
       AUTHORITY["EPSG","6267"]],
     PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],
     UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9108"]],
     AUTHORITY["EPSG","4267"]],
     PROJECTION["Lambert_Conformal_Conic_2SP"],
     PARAMETER["standard_parallel_1",35.46666666666667],
     PARAMETER["standard_parallel_2",34.03333333333333],
     PARAMETER["latitude_of_origin",33.5],
     PARAMETER["central_meridian",-118],
     PARAMETER["false_easting",2000000],
     PARAMETER["false_northing",0],
     UNIT["unknown",0.3048006096012192]]
 }}}

 You can see that the units are set to unknown and not feet.

 Probably this could be fixed in OGR... Frank?

-- 
Ticket URL: <http://trac.osgeo.org/qgis/ticket/1537#comment:12>
Quantum GIS <http://qgis.org>
Quantum GIS is an Open Source GIS viewer/editor supporting OGR, PostGIS, and GRASS formats