[GRASS-dev] [GRASS GIS] #821: g.region vect=3Dmap breaks region

[GRASS GIS]

#821: g.region vect=3Dmap breaks region
----------------------+-----------------------------------------------------
 Reporter:  hamish    |       Owner:  grass-dev at lists.osgeo.org
     Type:  defect    |      Status:  new                      
 Priority:  major     |   Milestone:  6.4.0                    
Component:  default   |     Version:  svn-develbranch6         
 Keywords:  g.region  |    Platform:  Linux                    
      Cpu:  x86-32    |  
----------------------+-----------------------------------------------------
 Hi,

 In a lat/lon location if you import a lon,lat,z file:
 {{{
 G65> awk '{print $1 "|" $2 "|" $3}' fort.100 | \
        v.in.ascii out=fort_100 x=1 y=2 z=3 -t -z
 }}}

 then set the region to match the vector map's bounds:
 {{{
 G65> g.region v=fort_100
 }}}

 the region will be broken:
 {{{
 G65> g.region -p
 ERROR: region for current mapset is invalid
        line 18: <t-b resol: 1:00:00.0196>
        run "g.region"
 }}}

 removing this from the end of $MAPSET/WIND gets you operational again:
 {{{
 top:        4560.63
 bottom:     -7.39548
 cols3:      8
 rows3:      9
 depths:     4568
 e-w resol3: 1:02:32.898867
 n-s resol3: 1:03:00.119988
 t-b resol:  1:00:00.0196
 }}}


 t-b resolution (elevation) obviously shouldn't be in degrees..
 ?


 haven't had the connectivity to test 6.4 or 7 yet.


 thanks,
 Hamish

-- 
Ticket URL: <https://trac.osgeo.org/grass/ticket/821>
GRASS GIS <http://grass.osgeo.org>