I downloaded this scene from earthexplorer, and everything is fine:
$ gdalinfo LC80080472013174LGN00_B1.TIF
Driver: GTiff/GeoTIFF
Files: LC80080472013174LGN00_B1.TIF
./LC80080472013174LGN00_MTL.txt
Size is 7571, 7411
Coordinate System is:
PROJCS["WGS 84 / UTM zone 19N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-69],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
AUTHORITY["EPSG","32619"]]
Origin = (145185.000000000000000,2189715.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Metadata:
AREA_OR_POINT=Point
METADATATYPE=ODL
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 145185.000, 2189715.000) ( 72d23' 7.97"W, 19d46'16.43"N)
Lower Left ( 145185.000, 1967385.000) ( 72d20'44.51"W, 17d45'55.15"N)
Upper Right ( 372315.000, 2189715.000) ( 70d13' 8.50"W, 19d47'56.92"N)
Lower Right ( 372315.000, 1967385.000) ( 70d12'16.74"W, 17d47'24.74"N)
Center ( 258750.000, 2078550.000) ( 71d17'19.50"W, 18d47' 4.68"N)
Band 1 Block=7571x1 Type=UInt16, ColorInterp=Gray
$ gdalinfo --version
GDAL 2.1.0dev, released 2015/99/99
What is your gdal version? Maybe the Landsat geotiff was modified from
some software?
Best regards,
Dmitry
17.01.2016 20:15, Nevzat Guler пишет:
Thank you Dmitry,
What is the best way to extract the information we got using gdalinfo.exe by
using GDAL classes, within a C++ code? As I mentioned, I am using
GetGeoTransform and GetProjectionRef() but the results miss many of the
information (specifically the map info elements as I explained in the previous
email).
On another note, I actually get a different output from gdalinfo.exe using my
file (please see the printed result below). I got the geotiff file from my
colleague. When I open the file with ENVI, I get all the information, including
GEOGCS information.
Can it be that my file is just corrupted and that's why I don't get all the
information from GetGeoTransform and GetProjectionRef()? But then how would
ENVI is able to extract this information if it were not in the file?
Thanks again, - Nevzat
C:\Cubes\Landsat\LC80080472013174LGN00>gdalinfo LC80080472013174LGN00_B1.TIF
Driver: GTiff/GeoTIFF
Files: LC80080472013174LGN00_B1.TIF
.\LC80080472013174LGN00_MTL.txt
Size is 7571, 7411
Coordinate System is:
PROJCS["unnamed",
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-69],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (145185.000000000000000,2189715.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Metadata:
AREA_OR_POINT=Point
METADATATYPE=ODL
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 145185.000, 2189715.000)
Lower Left ( 145185.000, 1967385.000)
Upper Right ( 372315.000, 2189715.000)
Lower Right ( 372315.000, 1967385.000)
Center ( 258750.000, 2078550.000)
Band 1 Block=7571x1 Type=UInt16, ColorInterp=Gray
-----Original Message-----
From: gdal-dev [mailto:[email protected]] On Behalf Of Dmitry
Baryshnikov
Sent: Sunday, January 17, 2016 6:57 AM
To: [email protected]
Subject: Re: [gdal-dev] map info and projections from GeoTIFF file to create
ENVI header
Hi Nevzat,
The landsat geotiff already have CRS embeded in file
$ gdalinfo LC80150362013097LGN03_B1.TIF
Driver: GTiff/GeoTIFF
Files: LC80150362013097LGN03_B1.TIF
./LC80150362013097LGN03_MTL.txt
Size is 7281, 7261
Coordinate System is:
PROJCS["WGS 84 / UTM zone 17N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
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"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
AUTHORITY["EPSG","32617"]]
Origin = (706185.000000000000000,3943515.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Metadata:
AREA_OR_POINT=Point
METADATATYPE=ODL
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 706185.000, 3943515.000) ( 78d43'24.74"W, 35d36'50.15"N) Lower Left ( 706185.000,
3725685.000) ( 78d46'35.65"W, 33d39' 3.82"N) Upper Right ( 924615.000, 3943515.000) ( 76d19' 2.05"W,
35d32'39.74"N) Lower Right ( 924615.000, 3725685.000) ( 76d25'33.46"W, 33d35'10.99"N)
Center ( 815400.000, 3834600.000) ( 77d33'38.95"W, 34d36'17.35"N)
Band 1 Block=7281x1 Type=UInt16, ColorInterp=Gray
And this is right CRS.
Where do you get you geotiff?
Best regards,
Dmitry
17.01.2016 03:41, Nevzat Guler пишет:
Dear experts,
I am reading LANDSAT GeoTIFF files and combining bands to create and
ENVI cube.
I am using GDAL to read the GeoTIFF file, and my own custom cube
writer to write the ENVI style cube. Transferring the map information
proved to be a bit convoluted. I seek your help to figure out this.
My first problem is to access map info from GeoTIFF file. The map info
includes:
map info = {UTM, 1.000, 1.000, 145185.000, 2189715.000,
3.0000000000e+001, 3.0000000000e+001, 19, North, WGS-84, units=Meters}
I am able to get the location (x,y), and pixel size (x,y) from :
double adfGeoTransform[6];
srcDS->GetGeoTransform(adfGeoTransform);
But, I could not find any way to get the other variables in the
string: UTM, pixel tie points, Unit and DATUM. Is there any way to get them?
Finally, for the Coordinate System String, GDAL gives me:
PROJCS["unnamed",PROJECTION["Transverse_Mercator"],PARAMETER["latitude
_of_or
igin
",0],PARAMETER["central_meridian",-69],PARAMETER["scale_factor",0.9996
],PARA
METE
R["false_easting",500000],PARAMETER["false_northing",0],UNIT["Meter",1
]]
But I should actually get
PROJCS["UTM_Zone_19N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS
84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,
0,0],A
UTHO
RITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UN
IT["de
gree
",0.0174532925199433,AUTHORITY["EPSG","9108"]],AUTHORITY["EPSG","4326"
]],PRO
JECT
ION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER
["cent
ral_
meridian",-69],PARAMETER["scale_factor",0.9996],PARAMETER["false_easti
ng",50
0000
],PARAMETER["false_northing",0],UNIT["Meter",1]]
GDAL gives me PROJCS, "unnamed". And no GEOGCS information is given.
I am able to get the coordinate system properties by using the
following commands but I first need to set few things to use exportToWTK.
OGRSpatialReference oSRS;
char *pszWKT = NULL;
oSRS.SetProjCS("UTM_Zone_19N");
oSRS.SetWellKnownGeogCS("WGS84");
oSRS.SetUTM(19, TRUE);
oSRS.SetWellKnownGeogCS("WGS84");
oSRS.exportToWkt(&pszWKT);
printf("%s\n", pszWKT);
However, this requires me to set SetProjCS, SetWellKnownGeogCS and
SetUTM first. I don't know these before I read the file. Is there any
way to get them from the file so that I don't need prior knowledge to
create the header. When I pull the GeoTIFF file into ENVI, all these
information is already there. If this list is only or development
purposes, please let me know if there is any other email list that can
give me guidance on this quest.
I appreciate your help and suggestions.
Best regards,
- Nevzat
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
