roblabs / gdal-geopdf Goto Github PK

How to convert the El Capitan GeoTIFF to California Albers

• Get the Forest Service GeoTIFF from El Capitan
• Confirm that the projection is epsg.io/4269
• Convert using GDAL to epsg.io/3310

# Get the Forest Service GeoTIFF for El Capitan
wget https://data.fs.usda.gov/geodata/rastergateway/data/37119/fstopo/373711937_El_Capitan_FSTopo.tif
# epsg.io/4269 
gdalinfo 373711937_El_Capitan_FSTopo.tif | grep AUTHORITY


# Convert using GDAL to epsg.io/3310
gdalwarp -t_srs EPSG:3310 373711937_El_Capitan_FSTopo.tif el-cap.3310.tif

While downloading National Forest System Trails from FSTopo, some of the Pacific Crest Trail is missing in certain National Forests.

• Inyo National Forest near Mount Whitney and Crabtree Ranger Station
  • Shape file data for the PCT is in FSTopo Transportation Line and have FCSUBTYPE of 106, 107 and 525. See #1 for details.
  • No shape file data for National Forest System Trails

PCT Trail data near Mount Whitney and John Muir Trail is only in FSTopo Transportation Line

• S_USA.FSTopo_Transport_LN in red
• S_USA.TrailNFS_Publish in black with no shape file nor trail data

Summary

GeoPDF and Vector Shape file data for 7.5' Palomar Observatory quad show there are differences between the USGS and USFS topo maps.

Filing this issue as a tracking issue to share with others who may be using the US Topo PDF, FSTopo PDF and FSTopo vector data.

See live map at RobLabs.com/usgs

Issues

• Is there any documentation on the FCSUBTYPE? An issue filed by JesseCrocker suggests that FCSUBTYPE = 107 & FCSUBTYPE = 525 are significant: trailbehind/USFSTrailImport#1
• Entirely possible that the next revision of US Topo for Palomar Mountain (and others) will have the update.
• Is it possible that the the USGS vector data is sourced from the FSTopo Transportation Line and are mis-categorized in the FSTopo vector data set?

Observations

• USFS legend shows roads near Barker Valley trailhead as "Unimproved, 4 wheel drive"
• USGS legend shows roads near Barker Valley trailhead as "Local Road", but it does have a legend option for "4WD"
  • USGS does not show the Barker Valley trail
  • USFS Vector file, FSTopo Transportation Line (Date of last refresh: Oct 25, 2017) shows the Barker Valley trail labeled as FCSUBTYPE = 107
  • FSTopo Transportation Line Source

QGIS view of FSTopo Transportation Line, FCSUBTYPE = 107

QGIS view of FSTopo Transportation Line, FCSUBTYPE = 106

Barker Valley trail is just below this highlighted section of the Palomar Divide Road.

FSTopo PALOMAR OBSERVATORY, CA

• This product generated and downloaded on 2/2/2017
• wget https://data.fs.usda.gov/geodata/rastergateway/data/33116/fstopo/331511645_Palomar_Observatory_FSTopo.pdf
  

US Topo PALOMAR OBSERVATORY, CA

• Modification date from PDF proprieties — Mar 28, 2016, 8:40 AM
• wget https://prd-tnm.s3.amazonaws.com/StagedProducts/Maps/USTopo/1/22340/7541782.pdf

Issue

When processing the 2018 Geospatial PDFs from the nationalMap.gov an error occurs that was not seen when processing the same quadrangle from 2015. Posting this issue in case someone else has any research into what the problem may be.

• What — Processing USGS US Topos from 2018 with gdal_translate
• Error — produces ERROR 1: insufficient arguments for Marked Content
• Compare — 2015 USGS US Topo does not produce this same error
• Compute time — the 2015 US Topo also takes less time complete in GDAL compared to 2018's version

Steps to reproduce

2015 GeospatialPDF

# wget
wget 
https://prd-tnm.s3.amazonaws.com/StagedProducts/Maps/USTopo/PDF/CA/CA_Descanso_20150223_TM_geo.pdf

#gdal
gdal_translate CA_Descanso_20150223_TM_geo.pdf CA_Descanso_20150223_TM_geo.tif -co COMPRESS=LZW --config GDAL_PDF_DPI 200

2018 GeospatialPDF

# wget 
wget 
https://prd-tnm.s3.amazonaws.com/StagedProducts/Maps/USTopo/PDF/CA/CA_Descanso_20180904_TM_geo.pdf

#gdal
gdal_translate CA_Descanso_20180904_TM_geo.pdf CA_Descanso_20180904_TM_geo.tif -co COMPRESS=LZW --config GDAL_PDF_DPI 200

# Errors many times; similar to this
ERROR 1: Pos = 47918296, insufficient arguments for Marked Content
ERROR 1: Pos = 47918296, insufficient arguments for Marked Content
ERROR 1: Pos = 47918317, insufficient arguments for Marked Content
ERROR 1: Pos = 47918327, insufficient arguments for Marked Content
ERROR 1: Pos = 47918327, insufficient arguments for Marked Content
ERROR 1: Pos = 47918329, insufficient arguments for Marked Content
ERROR 1: Pos = 47918334, insufficient arguments for Marked Content
ERROR 1: Pos = 47918334, insufficient arguments for Marked Content
ERROR 1: Pos = 47918334, insufficient arguments for Marked Content
ERROR 1: Pos = 47918339, insufficient arguments for Marked Content
ERROR 1: Pos = 47918343, insufficient arguments for Marked Content
ERROR 1: Pos = 47920226, insufficient arguments for Marked Content
ERROR 1: Pos = 47920231, insufficient arguments for Marked Content
ERROR 1: Pos = 47920268, insufficient arguments for Marked Content
ERROR 1: Pos = 47920268, insufficient arguments for Marked Content
ERROR 1: Pos = 47920271, insufficient arguments for Marked Content
ERROR 1: Pos = 47946009, insufficient arguments for Marked Content
ERROR 1: Pos = 47946015, insufficient arguments for Marked Content
ERROR 1: Pos = 47947917, insufficient arguments for Marked Content
ERROR 1: Pos = 47954326, insufficient arguments for Marked Content