[QGIS-trac] [Quantum GIS] #3672: region error when creating a GRASS location?

[Quantum GIS]

#3672: region error when creating a GRASS location?
------------------------------------------------+---------------------------
   Reporter:  lutra                             |              Owner:               
       Type:  bug                               |             Status:  new          
   Priority:  major: does not work as expected  |          Milestone:  Version 1.7.0
  Component:  GRASS                             |            Version:  Trunk        
   Keywords:                                    |   Platform_version:               
   Platform:  All                               |           Must_fix:  Yes          
Status_info:  0                                 |  
------------------------------------------------+---------------------------
 Hi, I have seen this on both Linux and Windows.

 If I open a "old" mapset (don't remember exactly when it was created),
 that was created with EPSG 40000, '''g.region -p''' will return as
 expected


 {{{
 projection: 99 (Transverse Mercator)
 zone:       0
 datum:      rome40
 ellipsoid:  international
 north:      4834913.78536118
 south:      4803364.45415656
 west:       1650370.47056664
 east:       1689398.62729587
 nsres:      24.99947005
 ewres:      25.00202225
 rows:       1262
 cols:       1561
 cells:      1969982
 }}}


 if I create new location with EPSG 40000 then g.region -p will return


 {{{
 projection: 99 (Transverse Mercator)
 zone:       0
 datum:      towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68
 ellipsoid:  international
 north:      5336010
 south:      4004460
 west:       1240010
 east:       2430820
 nsres:      1191.01073345
 ewres:      1192.002002
 rows:       1118
 cols:       999
 cells:      1116882
 }}}

 that I don't know if it is correct.


 Now if I create a location in EPSG 3003 (that should be 40000==3003+towgs)
 the result of g.region -p is


 {{{
 projection: 99 (Transverse Mercator)
 zone:       0
 datum:      ** unknown (default: WGS84) **
 ellipsoid:  international
 north:      5336090
 south:      4004520
 west:       1239980
 east:       2430820
 nsres:      1191.02862254
 ewres:      1190.84
 rows:       1118
 cols:       1000
 cells:      1118000
 }}}


 that seems wrong.


 The same result happens also with other projection, for example epsg 3763
 (or all other CRS system I use normally)


 {{{

 projection: 99 (Transverse Mercator)
 zone:       0
 datum:      ** unknown (default: WGS84) **
 ellipsoid:  grs80
 north:      -497.84044931
 south:      -292994.94849781
 west:       -98340.22283752
 east:       105157.76515519
 nsres:      4.99995056
 ewres:      4.99995056
 rows:       58500
 cols:       40700
 cells:      2380950000
 }}}

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