[GRASS-user] Wrong projection info after importing Geotif

[Agustin Lobo]

Markus,

I made the mapset using the QGIS plugin Grass /New mapset.
I'm not sure if I made something wrong or it was a problem of
the plugin(I'll clarify this), but the  fact is that I ended up with
this (wrong) region:
~$ g.region -p

GRASS_INFO_ERROR(9537,1): region for current mapset is invalid

GRASS_INFO_ERROR(9537,1): line 9: <e-w resol:  0>

GRASS_INFO_ERROR(9537,1): run "g.region"

First problem is that I cannot fix this region as explained in
http://lists.osgeo.org/pipermail/grass-user/2008-October/047198.html

Second problem is that I run r.in.gdal as I mentioned to import the file 
that I've put
in  http://aloboaleu.googlepages.com/forMarkus
and get a wrong region for the raster as well:
r.in.gdal input=/media/mifat32/MASTER_ICTA2007_2008/2008b/npp.tif output=npp
~$ g.region rast=npp
~$ g.region -p

GRASS_INFO_ERROR(19559,1): region for current mapset is invalid

GRASS_INFO_ERROR(19559,1): line 9: <e-w resol:  0>

Nevertheless, according to gdalinfo, the geotif file is correct
(see my last message for the gdalinfo output).

  Probably, r.in.gdal would have worked fine if the default region had 
been correct, but this is no the point (I think). r.in.gdal should set 
the region for the raster according to the geotiff settings, which are 
correct.

So there is something wrong with r.in.gdal?

Thanks for your patience!

Agus

Markus Neteler wrote:
> Agus,
> 
> On Sat, Oct 25, 2008 at 10:48 AM, Agustin Lobo <Agustin.Lobo at ija.csic.es> wrote:
>> 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?
> 
> since I still don't know which steps you really did (apparently
> you created the location from the GeoTIFF file?)
> I don't know how to help. Maybe I just missed that.
> 
>> 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
> 
> OK, but in QGIS, what did you do to set up the default region?
> 
>> and, by some reason, was wrong and cannot modify it)
> 
> To reproduce it, I would need the steps in a detailed way...
> (say, due to work overload I am most motivated with the
> full procedure :-).
> 
> 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