[gdal-dev] Reprojecting GeoTIFFs with Overviews?

[Jason Greenlaw - NOAA Affiliate]

Using gdalwarp I'm attempting to reproject a 4-band (RGBA) GeoTIFF (in
EPSG:4269) containing 4 overviews to web mercator (EPSG:3857).  While the
reprojection seems to work fine for the image itself, it appears the
overviews are not being reprojected/carried over to the new reprojected
GeoTIFF.

I've tried specifying -co COPY_SRC_OVERVIEWS=YES to gdalwarp, but it still
results in the overviews being lost.

I am using GDAL 1.9.2.

Is this the expected behavior?  Is there perhaps some way to ask gdalwarp
to also reproject the overviews?  Or do I have to reproject each overview
individually, then merge them into a single GeoTIFF with gdal_translate?

Any help/information would be greatly appreciated.

Here is the gdalwarp command I'm using (esri-3857.prf contains custom
EPSG:3857 WKT to be compatible with ESRI software):

---------------------------------------------------------------------
gdalwarp -s_srs EPSG:4269 -t_srs esri-3857.prf -r bilinear -co COMPRESS=LZW
-co COPY_SRC_OVERVIEWS=YES -overwrite NGOFS_WL_4269.tif NGOFS_WL_3857.tif
---------------------------------------------------------------------

Below are the output of gdalinfo for both the source and reprojected
GeoTIFFs.

---------------------------------------------------------------------
$ gdalinfo NGOFS_WL_4269.tif
Driver: GTiff/GeoTIFF
Files: NGOFS_WL_4269.tif
Size is 3000, 1255
Coordinate System is `'
Origin = (-97.401951999999994,30.701952500000001)
Pixel Size = (0.003900000000000,-0.003903000000000)
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( -97.4019520,  30.7019525)
Lower Left  ( -97.4019520,  25.8036875)
Upper Right ( -85.7019520,  30.7019525)
Lower Right ( -85.7019520,  25.8036875)
Center      ( -91.5519520,  28.2528200)
Band 1 Block=3000x1 Type=Byte, ColorInterp=Red
  Overviews: 1500x627, 750x313, 375x156, 187x78
  Mask Flags: PER_DATASET ALPHA
  Overviews of mask band: 1500x627, 750x313, 375x156, 187x78
Band 2 Block=3000x1 Type=Byte, ColorInterp=Green
  Overviews: 1500x627, 750x313, 375x156, 187x78
  Mask Flags: PER_DATASET ALPHA
  Overviews of mask band: 1500x627, 750x313, 375x156, 187x78
Band 3 Block=3000x1 Type=Byte, ColorInterp=Blue
  Overviews: 1500x627, 750x313, 375x156, 187x78
  Mask Flags: PER_DATASET ALPHA
  Overviews of mask band: 1500x627, 750x313, 375x156, 187x78
Band 4 Block=3000x1 Type=Byte, ColorInterp=Alpha
  Overviews: 1500x627, 750x313, 375x156, 187x78
---------------------------------------------------------------------
---------------------------------------------------------------------
$ gdalinfo NGOFS_WL_3857.tif
Driver: GTiff/GeoTIFF
Files: NGOFS_WL_3857.tif
Size is 2940, 1391
Coordinate System is:
PROJCS["WGS 84 / Pseudo-Mercator",
    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"]],
    PROJECTION["Mercator_1SP"],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0
+lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext
 +no_defs"],
    AUTHORITY["EPSG","3857"]]
Origin = (-10842735.698910875245929,3572289.254035749472678)
Pixel Size = (443.061998399919901,-443.061998399919901)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (-10842735.699, 3572289.254) ( 97d24' 7.03"W, 30d31'59.96"N)
Lower Left  (-10842735.699, 2955990.014) ( 97d24' 7.03"W, 25d39' 5.65"N)
Upper Right (-9540133.424, 3572289.254) ( 85d42' 1.72"W, 30d31'59.96"N)
Lower Right (-9540133.424, 2955990.014) ( 85d42' 1.72"W, 25d39' 5.65"N)
Center      (-10191434.561, 3264139.634) ( 91d33' 4.37"W, 28d 7'12.80"N)
Band 1 Block=2940x1 Type=Byte, ColorInterp=Red
  Mask Flags: PER_DATASET ALPHA
Band 2 Block=2940x1 Type=Byte, ColorInterp=Green
  Mask Flags: PER_DATASET ALPHA
Band 3 Block=2940x1 Type=Byte, ColorInterp=Blue
  Mask Flags: PER_DATASET ALPHA
Band 4 Block=2940x1 Type=Byte, ColorInterp=Alpha
---------------------------------------------------------------------

--
Jason Greenlaw
Contractor, ERT, Inc.
NOAA/NOS/OCS/CSDL
http://nowcoast.noaa.gov
Jason.Greenlaw at noaa.gov
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20130110/f0ddd8e2/attachment.html>