[gdal-dev] Reprojecting GeoTIFFs with Overviews?

[Norman Vine]

Jason

you need to make overviews for the reprojected image

see 
http://www.gdal.org/gdaladdo.html



On Jan 10, 2013, at 12:20 PM, Jason Greenlaw - NOAA Affiliate <jason.greenlaw at noaa.gov> wrote:

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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20130110/fecba09a/attachment-0001.html>