Opened 15 years ago

Closed 15 years ago

#821 closed defect (fixed)

g.region vect=3Dmap breaks region

Reported by: hamish Owned by: grass-dev@…
Priority: major Milestone: 6.4.0
Component: Default Version: svn-develbranch6
Keywords: g.region Cc:
CPU: x86-32 Platform: Linux

Description

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

Change History (1)

comment:1 by hamish, 15 years ago

Resolution: fixed
Status: newclosed
  • this bug was also in grass 6.4 & trunk.
  • it only happened with a 3D vector in a lat/lon location
  • now fixed by r40016-8

Hamish

Note: See TracTickets for help on using tickets.