[GRASS-user] importing geotif in grass, transformation parameters

[Martina Schaefer]

Hi Paul,

here the infos that you asked me for:

>> "WARNING: No projection name! Projection parameters likely to be
>>         meaningless.
> That is of much more concern. Can you try running
> gdalinfo filename.tif
> on the file that has the problem, and posting the results to the list 
> here in the hope that someone can see what is wrong?

here the output
GRASS 6.2.2 (BelcherBathy):~/Work/Devon-Data/bthy > gdalinfo 
00005latlong.tif
Driver: GTiff/GeoTIFF
Size is 15134, 2496
Coordinate System is:
PROJCS["unnamed",
     GEOGCS["unnamed",
         DATUM["unknown",
             SPHEROID["unretrievable - using WGS84",6378137,298.257223563]],
         PRIMEM["Greenwich",0],
         UNIT[,0.0174532925199433]],
     UNIT["unknown",1]]
Origin = (-81.327157999999997,75.767208999999994)
Pixel Size = (0.000049996696181,-0.000049979967949)
Image Structure Metadata:
   COMPRESSION=LZW
Corner Coordinates:
Upper Left  ( -81.3271580,  75.7672090)
Lower Left  ( -81.3271580,  75.6424590)
Upper Right ( -80.5705080,  75.7672090)
Lower Right ( -80.5705080,  75.6424590)
Center      ( -80.9488330,  75.7048340)
Band 1 Block=15134x1 Type=Byte, ColorInterp=Red
Band 2 Block=15134x1 Type=Byte, ColorInterp=Green
Band 3 Block=15134x1 Type=Byte, ColorInterp=Blue
Band 4 Block=15134x1 Type=Byte, ColorInterp=Alpha


in fact, this is the same data as the following one, only transformed to 
LatLong by ArcGIS

GRASS 6.2.2 (BelcherBathy):~/Work/Devon-Data/bthy > gdalinfo 00005_utm.tif
Driver: GTiff/GeoTIFF
Size is 14805, 2508
Coordinate System is:
PROJCS["WGS 84 / UTM zone 17N",
     GEOGCS["WGS 84",
         DATUM["WGS_1984",
             SPHEROID["WGS 84",6378137,298.2572235630016,
                 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",-81],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",500000],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]],
     AUTHORITY["EPSG","32617"]]
Origin = (490944.258245616627391,8409263.051821239292622)
Pixel Size = (1.414845022256486,-5.571261774748564)
Metadata:
   AREA_OR_POINT=Area
   TIFFTAG_SOFTWARE=IMAGINE TIFF Support
Copyright 1991 - 1999 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: etif.c $ $Revision: 1.9.1.3 $ $Date: 2002/07/29 15:39:06EDT $
   TIFFTAG_XRESOLUTION=1
   TIFFTAG_YRESOLUTION=1
   TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
   COMPRESSION=LZW
Corner Coordinates:
Upper Left  (  490944.258, 8409263.052)
Lower Left  (  490944.258, 8395290.327)
Upper Right (  511891.039, 8409263.052)
Lower Right (  511891.039, 8395290.327)
Center      (  501417.649, 8402276.690)
Band 1 Block=14805x1 Type=Byte, ColorInterp=Red
Band 2 Block=14805x1 Type=Byte, ColorInterp=Green
Band 3 Block=14805x1 Type=Byte, ColorInterp=Blue
Band 4 Block=14805x1 Type=Byte, ColorInterp=Undefined



and from the following file I got this error message
"WARNING: No projection name! Projection parameters likely to be
          meaningless.
A datum name nad83 (North_American_Datum_1983) was specified without 
transformation parameters. "
And even if I can open the file, I can't for example project it to 
LongLat, everytime I get error messages about this transformation 
parameter. I suppose the problem is the "unnamed" in the PROJCS, but how 
could I fill it and with what?
This file is orignially a img file that I can open without problems with 
ArcGIS, but grass says
"GDALOpen failed - 4
`oa16787_058.img' not recognised as a supported file format."
so I converted it to tif by img2tif


GRASS 6.2.2 (BelcherBathy):~/Work/Devon-Data/new > gdalinfo oa16787_058.tif
Driver: GTiff/GeoTIFF
Size is 2357, 2264
Coordinate System is:
PROJCS["unnamed",
     GEOGCS["NAD83",
         DATUM["North_American_Datum_1983",
             SPHEROID["unnamed",6378137,298.2572220960423],
             AUTHORITY["EPSG","6269"]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433],
         AUTHORITY["EPSG","4269"]],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]],
     AUTHORITY["EPSG","26817"]]
Origin = (489494.411800000001676,8397247.058800000697374)
Pixel Size = (5.000000000000000,-5.000000000000000)
Metadata:
   AREA_OR_POINT=Area
Corner Coordinates:
Upper Left  (  489494.412, 8397247.059)
Lower Left  (  489494.412, 8385927.059)
Upper Right (  501279.412, 8397247.059)
Lower Right (  501279.412, 8385927.059)
Center      (  495386.912, 8391587.059)
Band 1 Block=64x64 Type=Byte, ColorInterp=Gray


Thanks for all good ideas!

	Martina