[GRASS-user] Wrong projection info after importing Geotif

[Agustin Lobo]

Thanks Markus,
but this should not be the correct way for r.in.gdal to work,
should it? I mean that the default is reading the
the geotif with its own region, which should not
be interfered by the current region unless you use
the -o option. Am I wrong?

Anyway, my problem for applying your solution is
described in message 
http://lists.osgeo.org/pipermail/grass-user/2008-October/047198.html

(I set up the default region using qgis and, by some reason, was wrong
and cannot modify it)

Agus

Markus Neteler wrote:
> Agustin,
> 
> I could reproduce the problem...:
> 
> GRASS 6.4.svn (latlong_tbe_climate):~ > r.in.gdal npp.tif out=npp
> ERROR: region for current mapset is invalid
>        line 9: <e-w resol: 0>
>        run "g.region"
> 
> # Now I reset the location which was still messed up from the previous
> # try to import npp_CruP_All_fine_geo.tif:
> 
> GRASS 6.4.svn (latlong_tbe_climate):~ > g.region -dp
> projection: 3 (Latitude-Longitude)
> zone:       0
> datum:      wgs84
> ellipsoid:  wgs84
> north:      50N
> south:      30N
> west:       5E
> east:       20E
> nsres:      1
> ewres:      1
> rows:       20
> cols:       15
> cells:      300
> 
> # Voilà - the import works again:
> 
> GRASS 6.4.svn (latlong_tbe_climate):~ > r.in.gdal npp.tif out=npp
> Projection of input dataset and current location appear to match
>  100%
> <npp> created
> r.in.gdal complete.
> 
> I assume that a 'g.region -d' which cure it also at your end.
> 
> Markus
> 

-- 
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: Agustin.Lobo at ija.csic.es
http://www.ija.csic.es/gt/obster