[QGIS-trac] Re: [Quantum GIS] #1079: on the fly projections do not work

[Quantum GIS]

#1079: on the fly projections do not work
--------------------------------------------------------------+-------------
        Reporter:  pcav                                       |         Owner:  mhugent      
            Type:  defect                                     |        Status:  reopened     
        Priority:  critical: causes crash or data corruption  |     Milestone:  Version 1.0.0
       Component:  Projection Support                         |       Version:  HEAD         
      Resolution:                                             |      Keywords:  srs.db datum 
Platform_version:                                             |      Platform:  All          
        Must_fix:  Yes                                        |   Status_info:  1            
--------------------------------------------------------------+-------------
Changes (by neteler):

  * status:  closed => reopened
  * resolution:  fixed =>

Comment:

 I confirm that QGIS does not understand data in a GRASS location with
 towgs parameter:

 {{{
 Warning: QgsCoordinateReferenceSystem::getRecord failed :  select * from
 tbl_srs where
 parameters='+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0
 +ellps=intl
 +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +units=m +no_defs'

 GRASS 6.4.svn (pat): > g.proj -w
 PROJCS["Transverse Mercator",
     GEOGCS["international",
         DATUM["Monte_Mario",
             SPHEROID["International_1924",6378388,297],
             TOWGS84[-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",9],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",1500000],
     PARAMETER["false_northing",0],
     UNIT["Meter",1]]

 g.proj -wf > /tmp/gb1.prj
 testepsg /tmp/gb1.prj
 Validate Succeeds.
 WKT[/tmp/gb1.prj] =
 PROJCS["Transverse Mercator",
     GEOGCS["international",
         DATUM["Monte_Mario",
             SPHEROID["International_1924",6378388,297],
             TOWGS84[-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",9],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",1500000],
     PARAMETER["false_northing",0],
     UNIT["Meter",1]]

 Simplified WKT[/tmp/gb1.prj] =
 PROJCS["Transverse Mercator",
     GEOGCS["international",
         DATUM["Monte_Mario",
             SPHEROID["International_1924",6378388,297],
             TOWGS84[-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",9],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",1500000],
     PARAMETER["false_northing",0],
     UNIT["Meter",1]]

 Old Style WKT[/tmp/gb1.prj] = PROJCS["Transverse
 Mercator",GEOGCS["international",DATUM["Monte_Mario",
 SPHEROID["International_1924",6378388,297]],PRIMEM["Greenwich",0],
 UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],
 PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",9],
 PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",1500000],
 PARAMETER["false_northing",0],UNIT["Meter",1]]

 ESRI'ified WKT[/tmp/gb1.prj] =
 PROJCS["Transverse Mercator",
     GEOGCS["international",
         DATUM["D_Monte_Mario",
             SPHEROID["International_1924",6378388,297]],
         PRIMEM["Greenwich",0],
         UNIT["Degree",0.017453292519943295]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",9],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",1500000],
     PARAMETER["false_northing",0],
     UNIT["Meter",1]]
 PROJ.4 rendering of [/tmp/gb1.prj] = +proj=tmerc +lat_0=0 +lon_0=9
 +k=0.9996 +x_0=1500000 +y_0=0 +ellps=intl
 +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +units=m +no_defs
 }}}

 Since GDAL considers this projection '''valid''', it must  be a QGIS
 problem. I take
 liberty to reopen the ticket.

 Markus

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