[GRASS-user] Error in Projection from Landsat images

[Micha Silver]

Hi

Kim Besson wrote:

> Greetings
> This topic was being discussed in converting Vector to raster and them 
> I realized that it might be better to change it to here:
> I did
> gdalinfo to my landsat data is this:
> Driver: GTiff/GeoTIFF
> Files: L.tif
> Size is 7662, 7109
> Coordinate System is:
> LOCAL_CS[" Geocoding information not available Projection
 From this it looks like gdal was not able to identify the coordinate 
system.
There are lots of ETRS_1989 projections. Do you have any other way to 
know what the original CRS is?
(BTW, original Landsat imagery are always in some UTM  zone coordinates. 
Is that the case here? or is this some derived and reprojected TIFF?)
Once you know the correct CRS, you'll be able to create a matching 
Location/mapset, import the vector layer, then convert to raster, etc...

> Name = ETRS_1989_Transverse_Mercator Units = meters Geo",
>     UNIT["unknown",1]]
> Origin = (-52163.169232331842068,23020.214927441673353)
> Pixel Size = (30.000000000000000,-30.000000000000000)
> Metadata:
>   TIFFTAG_SOFTWARE=IMAGINE TIFF Support
> Copyright 1991 - 1999 by ERDAS, Inc. All Rights Reserved
> @(#)$RCSfile: etif.c $ $Revision: 1.10.1.9.1.9.2.11 $ $Date: 
> 2004/09/15 18:42:01EDT $
>   TIFFTAG_XRESOLUTION=1
>   TIFFTAG_YRESOLUTION=1
>   TIFFTAG_RESOLUTIONUNIT=1 (unitless)
>   AREA_OR_POINT=Area
> Image Structure Metadata:
>   INTERLEAVE=PIXEL
> Corner Coordinates:
> Upper Left  (  -52163.169,   23020.215) 
> Lower Left  (  -52163.169, -190249.785) 
> Upper Right (  177696.831,   23020.215) 
> Lower Right (  177696.831, -190249.785) 
> Center      (   62766.831,  -83614.785) 
> Band 1 Block=64x64 Type=Byte, ColorInterp=Red
> Band 2 Block=64x64 Type=Byte, ColorInterp=Green
> Band 3 Block=64x64 Type=Byte, ColorInterp=Blue
> Band 4 Block=64x64 Type=Byte, ColorInterp=Undefined
> Band 5 Block=64x64 Type=Byte, ColorInterp=Undefined
> Band 6 Block=64x64 Type=Byte, ColorInterp=Undefined
>
> It has geographic information right? and it's supposed to be 
> recognized by GRASS.
> So I decided to create a new Location.
> 1- Define name of the location (landsat02)
> 2- Choose method for creating a new location: Read Projection and 
> datums from a georreferenced file==> I select the file
> I choose no in defining region
> 3- Enter into location
> 4- do g.region -p and I get:
> projection: 0 (x,y)
> zone:       0
> north:      23020.21492744
> south:      -190249.78507256
> west:       -52163.16923233
> east:       177696.83076767
> nsres:      30
> ewres:      30
> rows:       7109
> cols:       7662
> cells:      54469158
>
> Everything else seems ok except the Projection that is (x,y). This 
> means that it's not recognizing the coordinate system or am I doing 
> something wrong?
>
> Thank you
> Kim
>
>
>
> This mail was received via Mail-SeCure System.
> ------------------------------------------------------------------------
>
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>
> This mail was received via Mail-SeCure System.
>
>
>   


-- 
Micha Silver
Arava Development Co. +972-52-3665918
http://www.surfaces.co.il