[gdal-dev] Reprojecting GeoTIFFs with Overviews?
Jason Greenlaw - NOAA Affiliate
jason.greenlaw at noaa.gov
Thu Jan 10 10:03:33 PST 2013
Norman,
Unfortunately that's not an option - In this case, my overviews are
generated separately as individual .ovr files, e.g.:
NGOFS_WL.tif
NGOFS_WL.tif.ovr
NGOFS_WL.tif.ovr.ovr
NGOFS_WL.tif.ovr.ovr.ovr
NGOFS_WL.tif.ovr.ovr.ovr.ovr
I'm then stitching them together into a single image with gdal_translate
-co COPY_SRC_OVERVIEWS=YES before attempting to reproject the result.
Here's an example of what it looks like:
http://jgreenlaw.org/gdal/ngofs-wl.png
In each overview, the "streaklet" arrows actually need to be rendered at
the proper resolution; I can't simply interpolate the original image.
Perhaps I'm out of luck and need to reproject each of the .ovr files and
the .tif files separately before combining them?
Jason
--
Jason Greenlaw
Contractor, ERT, Inc.
NOAA/NOS/OCS/CSDL
http://nowcoast.noaa.gov
Jason.Greenlaw at noaa.gov
(603) 862-3560
On Thu, Jan 10, 2013 at 12:37 PM, Norman Vine <nhv at cape.com> wrote:
> 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/dcb9eede/attachment.html>
More information about the gdal-dev
mailing list