[GRASS-user] Wrong projection info after importing Geotif

[Agustin Lobo]

Markus, I already did use -d, I copy the output below. But
the problem that worries me most is not that the default region
be wrong, but that r.in.gdal gets confused and interfered by
the default or active region while I think that the only geographic info 
that should
matter for r.in.gdal is the one in the geotiff file,
which is correct.

$ g.region -d
$ g.region -p
projection: 3 (Latitude-Longitude)
zone:       0
datum:      wgs84
ellipsoid:  wgs84
north:      90N
south:      90S
west:       180E
east:       180E
nsres:      0:14:25.153538
ewres:      0:21:36
rows:       749
cols:       1000
cells:      749000
$ r.in.gdal input=/media/mifat32/MASTER_ICTA2007_2008/2008b/npp.tif 
output=npp

GRASS_INFO_MESSAGE(18845,1): Projection of input dataset and current 
location appear to match

$ g.region rast=npp
$ g.region -p

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

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

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

Thanks,
Agus

Markus Neteler wrote:
> Agus,
> 
> On Mon, Oct 27, 2008 at 10:17 PM, Agustin Lobo <Agustin.Lobo at ija.csic.es> wrote:
>> 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"
> 
> Yes (@devs: suggest the -d flag?)
> 
>> First problem is that I cannot fix this region as explained in
>> http://lists.osgeo.org/pipermail/grass-user/2008-October/047198.html
> 
> If you looks closely, there was suggested:
> 
> g.region -d
> 
> Did you really run that (g.region set to default region)?
> GRASS doesn't automatically recover from a illegal
> region setting.
> 
> Then I suppose that it will work.
> 
> Markus
> 
>> 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
>>
> 
> 
> 

-- 
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