[GRASS-user] Error while using r.in.gdal for USGS landsat7 images

Hamish hamish_b at yahoo.com
Thu Jul 8 07:41:40 EDT 2010


António wrote:
> I'm using r.in.gdal to import Landsat7 images downloaded
> from USGS website. These landsat images are distributed in
> UTM-WGS84 for the specific zone. in this case I'm working
> with images from Brazil (WGS84 UTM ZONE 24S). So first I
> created  an Location based on EPSG:32724.
> Here goes the gdalinfo command for the image:
> /Driver: GTiff/GeoTIFF
> Files: L5216065_06520070405_B10.TIF
> Size is 8171, 6991
> Coordinate System is:
> PROJCS["WGS 84 / UTM zone 24N",
>    GEOGCS["WGS 84",
>        DATUM["WGS_1984",
>            SPHEROID["WGS
> 84",6378137,298.257223563,
>            
>    AUTHORITY["EPSG","7030"]],
>        
>    AUTHORITY["EPSG","6326"]],
>        PRIMEM["Greenwich",0],
>    
>    UNIT["degree",0.0174532925199433],
>        AUTHORITY["EPSG","4326"]],
>    PROJECTION["Transverse_Mercator"],
>    PARAMETER["latitude_of_origin",0],
>    PARAMETER["central_meridian",-39],
>    PARAMETER["scale_factor",0.9996],
>    PARAMETER["false_easting",500000],
>    PARAMETER["false_northing",0],
>    UNIT["metre",1,
>        AUTHORITY["EPSG","9001"]],
>    AUTHORITY["EPSG","32624"]]
> Origin = (443999.999999999941792,-694500.000000000000000)
> Pixel Size = (30.000000000000000,-30.000000000000000)/
> 
> As can be seen it's North instead of South. Nevertheless,
> r.in.gdal imports it and it even says: Projection of input
> dataset and current location appear to match
> But when I try to display I get an error window stating:
> Execution failed: 'd.rast -o map=import at PERMANENT --q'
> Details:
> Error: Unable to open raster map <import at PERMANENT>
> And nothing else is stated


what does "g.list rast" say? how about r.info?


Hamish



      


More information about the grass-user mailing list