#821: g.region vect=3Dmap breaks region
----------------------+-----------------------------------------------------
Reporter: hamish | Owner: grass-dev@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>