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.green...@noaa.gov
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
