[gdal-dev] problem with a geotiff

[Thomas Juntunen]

Hi Philippe,

It works for me:
>>> from osgeo import gdal
>>> gdal.VersionInfo()
'2020100'
>>> ds = gdal.Open(r"E:\16MAR08101212-M1BS-055675519010_01_P001.TIF")
>>> ds
<osgeo.gdal.Dataset; proxy of <Swig Object of type 'GDALDatasetShadow *' at
0x00000000010A6CC0> >
>>> ds.GetProjectionRef()
'GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS
84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]]'

>From gdalinfo:
E:\>gdalinfo E:\16MAR08101212-M1BS-055675519010_01_P001.TIF
Driver: GTiff/GeoTIFF
Files: E:\16MAR08101212-M1BS-055675519010_01_P001.TIF
Size is 6616, 6903
Coordinate System is `'
GCP Projection =
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
GCP[  0]: Id=1, Info=
          (0.5,0.5) -> (13.1567807900763,32.9252312619272,0)
GCP[  1]: Id=2, Info=
          (6615.5,0.5) -> (13.2984103237501,32.9244584648933,0)
GCP[  2]: Id=3, Info=
          (6615.5,6902.5) -> (13.2983785314726,32.7993021397123,0)
GCP[  3]: Id=4, Info=
          (0.5,6902.5) -> (13.1568844718844,32.8000571640568,0)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_COPYRIGHT=(C) COPYRIGHT 2016 DigitalGlobe, Inc., Longmont CO USA
80503
  TIFFTAG_DATETIME=2016:09:02 06:26:36
  TIFFTAG_IMAGEDESCRIPTION={
  bandList =
  [
    2;
    3;
    4;
    5;
  ]
}
  TIFFTAG_MAXSAMPLEVALUE=2047
  TIFFTAG_MINSAMPLEVALUE=1
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 6903.0)
Upper Right ( 6616.0,    0.0)
Lower Right ( 6616.0, 6903.0)
Center      ( 3308.0, 3451.5)
Band 1 Block=6616x79 Type=UInt16, ColorInterp=Red
Band 2 Block=6616x79 Type=UInt16, ColorInterp=Green
Band 3 Block=6616x79 Type=UInt16, ColorInterp=Blue
Band 4 Block=6616x79 Type=UInt16, ColorInterp=Undefined

On Tue, Sep 25, 2018 at 11:59 AM, <mauge.philippe at free.fr> wrote:

> Please, I'd like to submit you an isssue about a geotiff file.
> Go to the web site http://www.digitalglobe.com/resources/product-samples/
> tripoli-lybya
> Choose to download the second file: System-Ready (1B), 50 cm
> Unzip it and try to load the file 16MAR08101212-P1BS-055675519010_01_P001.TIF
> in GDAL
> I can't load it in GDAL because no projection is returned by the fonction
> GetProjectionRef()
> It can be loaded only if you create manually a .prj and a .tfw file.
> The file is OK because it is correctly loaded in Global Mapper.
> Usually I can load geotiff files in GDAL, but this file won't.
> Can you try to load it and may be tell me what's wrong with it ?
>
> Best regards
> Philippe
>
>
> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient> Garanti
> sans virus. www.avast.com
> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient>
> <#m_117792916060688249_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>



-- 
*The right coordinate system can turn an impossible problem into two really
hard problems. -- Charlie Pellerin*

Thomas Juntunen
GIS Specialist
Polar Geospatial Center
R280 Learning and Environmental Science
1954 Buford Ave
University of Minnesota
St. Paul, MN 55108
612-626-0505
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20180925/ad28f4bb/attachment.html>