[mapserver-users] HOW-TO convert USGS topo to GTiff

[Mike Flannigan]

Nice write up.  Thanks for sharing that.
I didn't know about turning off layers.

This is the GDAL command I like to use.

gdal_translate -of GTiff bm.pdf bm.tif --config GDAL_PDF_DPI 300 -co 
"COMPRESS=JPEG" -co "JPEG_QUALITY=85"

It appears to provide a good trade-off between file
size and image quality.



Mike


On 3/27/2014 2:00 PM, mapserver-users-request at lists.osgeo.org wrote:
> I have had several people send me follow-up emails about converting the USGS eTopos (you can download these from store.usgs.gov), so I thought it might make sense to just post this to the list as a HOW-TO.  Apologies for cross-post.
>
> The challenge is to take the GeoPDF that you download from USGS and make it into a series of TIFFs you can mosaic together.  This requires at least GDAL 1.8 and for full functionality you want 1.10.
>
> First, the simplest syntax would be
>
> gdal_translate topo.pdf topo.tif
>
> That would give you a GeoTiff with all the layers turned on, rendered at 150 dpi.  Easy enough.
>
> However, depending on what you want, you may want to alter that.  For example the GeoPDF renders differently as you zoom in; it has scale-dependent rendering.  So you may want to change the dpi to a higher number.  I did all mine at 300 dpi which seemed to fit a 24k rendering in ArcMap and looks nice on the screen.
>
> gdal_translate --config GDAL_PDF_DPI 300 topo.pdf topo.tif
>
> That still has all the layers turned on though.  You may want some of them off.  For example, you may want to turn off the NAIP data and the UTM grid (like I did).  If so you can specify certain LAYERS to turn off, first you have to figure out what they are with this (this requires GDAL 1.10):
>
> gdal_info -mdd LAYERS topo.pdf
>
> will return the metadata and a list of layers something like this:
> <snip>
> Metadata (LAYERS):
>    LAYER_00_NAME=Map_Collar
>    LAYER_01_NAME=Map_Collar.Map_Elements
>    LAYER_02_NAME=Map_Frame
>    LAYER_03_NAME=Map_Frame.Projection_and_Grids
>    LAYER_04_NAME=Map_Frame.Projection_and_Grids.Projection_Coordinate_Values
>    LAYER_05_NAME=Map_Frame.Projection_and_Grids.Geographic_and_Grid_Ticks
>    LAYER_06_NAME=Map_Frame.Projection_and_Grids.Projection_Line_Mask
>    LAYER_07_NAME=Map_Frame.Projection_and_Grids.Grid_Lines
> <snip>
>    LAYER_28_NAME=Map_Frame.Land_Cover.Woodland
>    LAYER_29_NAME=Images
>    LAYER_30_NAME=Images.Orthoimage
> <snip>
>
> So I'm looking to get rid of layer 7 and layer 30.  The syntax for that uses the layer name like this:
>
> gdal_translate --config GDAL_PDF_DPI 300 --config GDAL_PDF_LAYERS_OFF Map_Frame.Projection_and_Grids.Grid_Lines,Images.Orthoimage topo.pdf topo.tif
>
> This only works for the new 2011/12/13 maps.  The historic maps are scanned images, so they don't have any real layers.  They are just one raster image in a PDF.  So you would just do
> gdal_translate --config GDAL_PDF_DPI 300 topo.pdf topo.tif.
>
> What you get in both cases is an output GeoTIFF of the topo map.  The map will have all the collar info, so you will then have to crop the images, remove the collars, and mosaic the slivers to get a seamless mosaic, but that is just done the same way as any other imagery collection.  That process is to simply:
>                  Crop the image to its minimum size with gdalwarp
>                  Burn some null value to the remaining collar slivers with gdal_rasterize -i -burn <etc.>
>                  Mosaic the neighboring tile data into the sliver with gdalwarp
>
> One other cool thing about the GeoPDF with its scale-dependent rendering is that you could create topos at different dpis (scales) and thus make your overviews look like what a 100K or 250K topo would look like (well with the new ones anyway, not the historic ones).
>
> Don't forget to RTF(riendly!)M http://www.gdal.org/frmt_pdf.html
>
>