[GRASS-user] Wrong projection info after importing Geotif

Markus Neteler neteler at osgeo.org
Thu Oct 23 08:22:52 EDT 2008


Agustin,

On Wed, Oct 22, 2008 at 8:19 PM, Agustin Lobo <aloboaleu at gmail.com> wrote:
> Markus,
>
> You are right that the file I'm using is not quite the
> one of the link that I included in my message. I actually
> passed that map to R, and just redifined:
>>b <-  npp_geotiff
>>b at data$band1[b at data$band1  <0 ]<- 0
>
> and then wrote b as npp.tif (I used QGIS manageR)
>
> But I do not
> think that this is a matter of concern. The point is that
> the geographic info displayed by gdal for npp.tif is correct
> (the "90d 2'30.00"S" problem is not in this file)

Can you please post this output, or, better, the npp.tif file?
It's hard to test without.

> while the region set by grass from it with
> g.region rast=npp.tif
> is not.
>
> I guess that the r.in.gdal I have is older (and wrong) than yours,

I don't think so (I am not aware of relevant changes).

> mine is the one in QGIS 1.0preview:
> GRASS > g.version
> GRASS 6.3.0 (2008)
>
> Thus I don't think there is any way to solve this, except
> installing a 6.4 version, import ant then copy to the LOCATIOn and mapset I
> would use with QGIS. Am I wrong?

Before doing that, let me try your file.

Markus


More information about the grass-user mailing list