[Gdal-dev] GeoTIFF from gdalwarp doesn't load properly on ArcMap

[Diego Berge]

Hello,

    After reprojecting a GTiff using gdalwarp, when trying to load the 
resulting file in ArcMap, it reports the file's projection information 
as "This file has no associated projection information".  It does 
*appear* (I still need to confirm, but the coordinates look sensible) to 
load the file in the correct place, though.

This is the input file data:
------------------------------------
$ gdalinfo etm165_47.tif
Driver: GTiff/GeoTIFF
Size is 15108, 14175
Coordinate System is:
PROJCS["WGS 84 / UTM zone 39N",
     GEOGCS["WGS 84",
         DATUM["WGS_1984",
             SPHEROID["WGS 84",6378137,298.2572235629972,
                 AUTHORITY["EPSG","7030"]],
             AUTHORITY["EPSG","6326"]],
         PRIMEM["Greenwich",0],
         UNIT["degree (supplier to define 
representation)",0.01745329251994328],
         AUTHORITY["EPSG","4326"]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",51.00000000000005],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",500000],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]],
     AUTHORITY["EPSG","32639"]]
Origin = (-135105.000000,2190330.000000)
Pixel Size = (15.00000000,-15.00000000)
Metadata:
   TIFFTAG_IMAGEDESCRIPTION=Blue raster from e165_47
Corner Coordinates:
Upper Left  ( -135105.000, 2190330.000) ( 44d56'54.52"E, 19d42'23.00"N)
Lower Left  ( -135105.000, 1977705.000) ( 45d 0'59.27"E, 17d47'43.24"N)
Upper Right (   91515.000, 2190330.000) ( 47d 6'10.94"E, 19d45'58.89"N)
Lower Right (   91515.000, 1977705.000) ( 47d 8'49.23"E, 17d50'56.76"N)
Center      (  -21795.000, 2084017.500) ( 46d 3'13.07"E, 18d46'56.70"N)
Band 1 Block=15108x1 Type=Byte, ColorInterp=Red
Band 2 Block=15108x1 Type=Byte, ColorInterp=Green
Band 3 Block=15108x1 Type=Byte, ColorInterp=Blue
------------------------------------

This is the transformation command:

$ gdalwarp -t_srs "+proj=lcc +lat_1=33N +lat_2=17N +lat_0=25d05'22.236N 
+lon_0=48E +x_0=0 +y_0=0 +a=6378388 +b=6356912 
+towgs84=-143.0,-236.0,7.0" etm165_47.tif lcc165_47.tif

I've also tried the roughly equivalent:

$ gdalwarp -t_srs EPSG:2318 etm165_47.tif lcc165_47_epsg.tif

This is the result in the first case (lcc165_47.tif):
------------------------------------
$ gdalinfo lcc165_47.tif
Driver: GTiff/GeoTIFF
Size is 15299, 14376
Coordinate System is:
PROJCS["unnamed",
     GEOGCS["unnamed",
         DATUM["unknown",
             SPHEROID["unnamed",6378388,297.0007450176942]],
         PRIMEM["Greenwich",0],
         UNIT[,0.0174532925199433]],
     PROJECTION["Lambert_Conformal_Conic_2SP"],
     PARAMETER["standard_parallel_1",33],
     PARAMETER["standard_parallel_2",17],
     PARAMETER["latitude_of_origin",25.08951000004958],
     PARAMETER["central_meridian",48],
     PARAMETER["false_easting",0],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]]]
Origin = (-318119.632160,-584327.455283)
Pixel Size = (14.89885658,-14.89885658)
Corner Coordinates:
Upper Left  ( -318119.632, -584327.455) ( 44d56'54.19"E, 19d44'10.26"N)
Lower Left  ( -318119.632, -798513.417) ( 44d59'38.71"E, 17d47'40.12"N)
Upper Right (  -90182.025, -584327.455) ( 47d 8'5.21"E, 19d45'58.29"N)
Lower Right (  -90182.025, -798513.417) ( 47d 8'51.87"E, 17d49'26.16"N)
Center      ( -204150.829, -691420.436) ( 46d 3'22.44"E, 18d47'0.65"N)
Band 1 Block=256x256 Type=Byte, ColorInterp=Red
   Overviews: 7650x7188, 3825x3594, 1913x1797, 957x899, 479x450, 240x225
Band 2 Block=256x256 Type=Byte, ColorInterp=Green
   Overviews: 7650x7188, 3825x3594, 1913x1797, 957x899, 479x450, 240x225
Band 3 Block=256x256 Type=Byte, ColorInterp=Blue
   Overviews: 7650x7188, 3825x3594, 1913x1797, 957x899, 479x450, 240x225
------------------------------------

And this is the result for the second case (lcc165_47_epsg.tif):
------------------------------------
$ gdalinfo lcc165_47_epsg.tif
Driver: GTiff/GeoTIFF
Size is 15299, 14375
Coordinate System is:
PROJCS["Ain el Abd / Aramco Lambert",
     GEOGCS["Ain el Abd",
         DATUM["unknown",
             SPHEROID["unretrievable - using WGS84",6378137,298.257223563,
                 AUTHORITY["EPSG","0"]],
             AUTHORITY["EPSG","6204"]],
         PRIMEM["Greenwich",0],
         UNIT["degree (supplier to define 
representation)",0.01745329251994328],
         AUTHORITY["EPSG","4204"]],
     PROJECTION["Lambert_Conformal_Conic_2SP"],
     PARAMETER["standard_parallel_1",17.00000000000002],
     PARAMETER["standard_parallel_2",33.00000000000004],
     PARAMETER["latitude_of_origin",25.08951000000002],
     PARAMETER["central_meridian",48.00000000000005],
     PARAMETER["false_easting",0],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]],
     AUTHORITY["EPSG","2318"]]
Origin = (-318184.383793,-584287.299047)
Pixel Size = (14.89931340,-14.89931340)
Corner Coordinates:
Upper Left  ( -318184.384, -584287.299) ( 44d56'51.48"E, 19d44'11.19"N)
Lower Left  ( -318184.384, -798464.929) ( 44d59'36.04"E, 17d47'41.21"N)
Upper Right (  -90239.788, -584287.299) ( 47d 8'3.08"E, 19d45'59.26"N)
Lower Right (  -90239.788, -798464.929) ( 47d 8'49.77"E, 17d49'27.29"N)
Center      ( -204212.086, -691376.114) ( 46d 3'20.04"E, 18d47'1.68"N)
Band 1 Block=15299x1 Type=Byte, ColorInterp=Red
Band 2 Block=15299x1 Type=Byte, ColorInterp=Green
Band 3 Block=15299x1 Type=Byte, ColorInterp=Blue
------------------------------------

The input file loads without any problems, btw.  If needed, I could put 
a reduced version of the files (originals are about 600Mb) on an ftp site.

Any ideas?

Regards,
Diego Berge.