Ecosyste.ms: Awesome

An open API service indexing awesome lists of open source software.

Awesome Lists | Featured Topics | Projects

https://github.com/roblabs/gdal-geopdf

Processing GeoPDFs using GDAL
https://github.com/roblabs/gdal-geopdf

forest-service gdal geospat-pdf geotiff processing-geopdfs quadrangle-maps topo usgs

Last synced: 5 days ago
JSON representation

Processing GeoPDFs using GDAL

Host: GitHub
URL: https://github.com/roblabs/gdal-geopdf
Owner: roblabs
License: mit
Created: 2017-10-14T03:30:09.000Z (about 7 years ago)
Default Branch: master
Last Pushed: 2021-10-21T18:24:55.000Z (about 3 years ago)
Last Synced: 2024-12-22T20:06:40.378Z (11 days ago)
Topics: forest-service, gdal, geospat-pdf, geotiff, processing-geopdfs, quadrangle-maps, topo, usgs
Size: 1.43 MB
Stars: 11
Watchers: 3
Forks: 1
Open Issues: 4
Metadata Files:
- Readme: readme.md
- License: LICENSE

Awesome Lists containing this project

README

# GeoPDFs

This project shows how to process GeoPDFs from the [USGS](https://viewer.nationalmap.gov/basic/) and [USFS](https://data.fs.usda.gov/geodata/rastergateway/states-regions/states.php).

## US Topo from US Geologic Survey

### US Topo of Wizard Island at Crater Lake National Park

![](assets/readme-wizard-island.png)

### Process Examples

#### `gdalinfo`

Inspecting the response from `gdalinfo` with the option `-mdd` yields some interesting metadata.

```bash
# get Geospatial PDF of Crater Lake West
# Metadata from nationalmap.gov
# USGS US Topo 7.5-minute map for Crater Lake West, OR 2017
# Published Date: 2017-03-23
# Format: GeoPDF (27.87 MB), Extent: 7.5 x 7.5 minute

# Make a temp directory and save as a temp file name
mkdir -p tmp
wget -O tmp/test.pdf https://prd-tnm.s3.amazonaws.com/StagedProducts/Maps/USTopo/PDF/OR/OR_Crater_Lake_West_20170323_TM_geo.pdf

gdalinfo -mdd LAYERS tmp/test.pdf

# output
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.Geographic_Names
LAYER_05_NAME=Map_Frame.Structures
LAYER_06_NAME=Map_Frame.Transportation

LAYER_12_NAME=Map_Frame.PLSS
LAYER_13_NAME=Map_Frame.Wetlands
LAYER_14_NAME=Map_Frame.Hydrography
LAYER_15_NAME=Map_Frame.Terrain
LAYER_16_NAME=Map_Frame.Terrain.Contours
LAYER_17_NAME=Map_Frame.Terrain.Shaded_Relief
LAYER_18_NAME=Map_Frame.Woodland
LAYER_19_NAME=Map_Frame.Boundaries

LAYER_29_NAME=Images
LAYER_30_NAME=Images.Orthoimage
LAYER_31_NAME=Barcode

```

##### PDF Dimensions

* PDF dimensions

```
Driver: PDF/Geospatial PDF
Files: tmp/test.pdf
Size is 3412, 4350
```

* Opening in a usual PDF viewer shows the image dimensions

```
22.75 × 29 inches
```

* Doing the math yields 150 pixels per inch.

#### `gdal_translate`

* Convert to GeoTIFF,

``` bash
DPI=150
gdal_translate tmp/test.pdf tmp/test.$DPI.tif \
-co COMPRESS=LZW \
--config GDAL_PDF_DPI $DPI

# Extract
gdalwarp tmp/test.$DPI.tif tmp/test.$DPI.wizard-island.tif \
-t_srs EPSG:4326 -dstalpha \
-co COMPRESS=LZW \
-te -122.17833 42.92361 -122.13799 42.95766

gdal_translate tmp/test.pdf tmp/test.$DPI.Shaded_Relief.tif \
-co COMPRESS=LZW \
--config GDAL_PDF_LAYERS "Map_Collar.Map_Elements,Map_Frame.Terrain.Shaded_Relief" \
--config GDAL_PDF_DPI $DPI

gdalwarp tmp/test.$DPI.Shaded_Relief.tif tmp/test.$DPI.Shaded_Relief.wizard-island.tif \
-co COMPRESS=LZW \
-t_srs EPSG:4326 -dstalpha \
-te -122.17833 42.92361 -122.13799 42.95766

```

### Links

* Larry Moore, USGS, 2016 — [Converting US Topo GeoPDF Layers to GeoTIFF](https://nationalmap.gov/ustopo/documents/ustopo2gtif_current.pdf)
* Andrew Burnes, FOSS4G North America, 2016 — [Using GDAL
to Translate
US Topo GeoPDFs](https://2016.foss4g-na.org/sites/default/files/slides/using-gdal-to-translate-us-topo-geopdf.pdf)
* [roblabs/gdal](https://hub.docker.com/r/roblabs/gdal/)

### Background
> The term “US Topo” refers specifically to quadrangle topographic maps published in 2009 and later. [source][1]

Excerpt from October 2017

> In 2017, the US Topo map production system was redesigned and modernized to provide a system that facilitates long term goals for more efficient production and continued product improvements. The new system produces maps in a format that uses a different georeferencing mechanism compliant with ISO 32000. The new products can continue to be viewed and printed with Adobe Reader or any comparable PDF viewing software. [source][2]

-----

## FSTopo from US Forest Service

### Wizard Island at Crater Lake National Park

![](assets/readme-fstopo-wizard-island.png)

### Process Examples

#### `gdalinfo`
```bash
# get GeoPDF of Crater Lake West
wget https://data.fs.usda.gov/geodata/rastergateway/data/42122/fstopo/425212207_Crater_Lake_West_FSTopo.pdf

gdalinfo -mdd LAYERS 425212207_Crater_Lake_West_FSTopo.pdf

```

Comments attempting to document the significant layers. Some sublayers removed; run `gdalinfo` for the full details.

```bash
Metadata (LAYERS):
LAYER_00_NAME=Other_5 # Declination, scale and legend
LAYER_01_NAME=Quadrangle_Extent

LAYER_06_NAME=Adjacent_Quadrangle_Diagram

LAYER_12_NAME=Other_2 # Quad label
LAYER_13_NAME=Quadrangle
LAYER_14_NAME=Quadrangle.Other
LAYER_15_NAME=Quadrangle.Labels

LAYER_17_NAME=Quadrangle.Neatline/Mask

LAYER_21_NAME=Quadrangle.Text
LAYER_22_NAME=Quadrangle.Road_Shields
LAYER_23_NAME=Quadrangle.Buildings
LAYER_24_NAME=Quadrangle.Culture_Features
LAYER_25_NAME=Quadrangle.Drainage_Features
LAYER_26_NAME=Quadrangle.Spot_Elevations
LAYER_27_NAME=Quadrangle.Geodetic_Control
LAYER_28_NAME=Quadrangle.Recreation/Facility_Features
LAYER_29_NAME=Quadrangle.Linear_Transportation_Features
LAYER_30_NAME=Quadrangle.Linear_Culture_Features
LAYER_31_NAME=Quadrangle.PLSS
LAYER_32_NAME=Quadrangle.Linear_Drainage_Features
LAYER_33_NAME=Quadrangle.Contour_Labels
LAYER_34_NAME=Quadrangle.Contours

LAYER_38_NAME=Quadrangle.UTM_Grid

LAYER_42_NAME=Quadrangle.Large_Buildings
LAYER_43_NAME=Quadrangle.Drainage_Polygons
LAYER_44_NAME=Quadrangle.WoodlandUSGS_P
```

#### `gdal_translate`

``` bash
DPI=300

gdal_translate 425212207_Crater_Lake_West_FSTopo.pdf 425212207_Crater_Lake_West_FSTopo.$DPI.tif \
-co COMPRESS=LZW \
--config GDAL_PDF_DPI $DPI

gdalwarp 425212207_Crater_Lake_West_FSTopo.$DPI.tif 425212207_Crater_Lake_West_FSTopo.$DPI.wizard-island.tif \
-co COMPRESS=LZW \
-t_srs EPSG:4326 -dstalpha \
-te -122.17833 42.92361 -122.13799 42.95766
```

### Links

* [FSTopo](https://data.fs.usda.gov/geodata/vector/index.php)
* [roblabs/gdal](https://hub.docker.com/r/roblabs/gdal/)

### Background

> FSTopo is the Forest Service Primary Base Map Series (1:24,000 scale for the lower 48 and Puerto Rico, 1:63,360 for Alaska) quadrangle maps. FSTopo products cover the US Forest Service lands.

[1]: https://nationalmap.gov/ustopo/index.html
[2]: https://nationalmap.gov/ustopo/about.html