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.000,-30.000)
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
