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


_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Reply via email to