National Forests layer

Properties:

  • forestname: Name of forest from source USFS GIS dataset
  • gis_acres: # of acres of forest from source USFS GIS dataset
  • geometry: Polygon/MultiPolygon geometry of forest from source USFS GIS dataset. Only polygons that intersect the trail of interest are included, but some national forests, like Inyo National Forest, are quite wide-ranging, so geometries are kept in Nevada I think because it has the same name.
  • length: length of trail in national forest in meters
  • wiki_image: url to Wikipedia image. I try to select the best image but it can be difficult sometimes.
  • wiki_url: url to wikipedia page
  • wiki_summary: summary of wikipedia page. Usually this is the first paragraph.
  • official_url: url to National Forest homepage
# Make temp directory
mkdir -p tmp
# Generate national forest polygons
python code/main.py export national-forests \
    `# trail code, i.e. 'pct'` \
    -t pct > tmp/nationalforests.geojson
# Generate point labels from those polygons
python code/main.py geom polylabel \
    `# include only the forestname attribute` \
    -y forestname \
    `# only keep labels for polygons that are >=30% of MultiPolygon area` \
    --rank-filter 0.2 \
    tmp/nationalforests.geojson > tmp/nationalforests_label.geojson

Run tippecanoe on the GeoJSON to create vector tiles 

rm -rf tmp/nationalforests_tiles
tippecanoe \
    `# Guess appropriate max zoom` \
    -zg \
    `# Export tiles to directory` \
    -e tmp/nationalforests_tiles \
    `# Input geojson` \
    -L'{"file":"tmp/nationalforests.geojson", "layer":"nationalforests"}' \
    -L'{"file":"tmp/nationalforests_label.geojson", "layer":"nationalforests_label"}'

Convert the exported metadata.json to a JSON file conforming to the Tile JSON spec 

python code/main.py util metadata-json-to-tile-json \
    `# Set tileset name` \
    --name 'National Forests' \
    `# Set attribution string` \
    --attribution '<a href="https://www.nps.gov/" target="_blank">© USFS</a>' \
    `# tile url paths` \
    --url 'https://tiles.nst.guide/pct/nationalforest/{z}/{x}/{y}.pbf' \
    `# Output file path` \
    -o tmp/nationalforests.json \
    `# input JSON file` \
    tmp/nationalforests_tiles/metadata.json
# Remove the unneeded `metadata.json`
rm tmp/nationalforests_tiles/metadata.json

Remove existing vector tiles 

aws s3 rm \
    --recursive \
    s3://tiles.nst.guide/pct/nationalforest/

Add new vector tiles 

aws s3 cp \
    tmp/nationalforests_tiles s3://tiles.nst.guide/pct/nationalforest/ \
    --recursive \
    --content-type application/x-protobuf \
    --content-encoding gzip \
    `# Set to public read access` \
    --acl public-read \
    `# two hour cache; one day swr` \
    --cache-control "public, max-age=7200, stale-while-revalidate=86400"
aws s3 cp \
    tmp/nationalforests.json s3://tiles.nst.guide/pct/nationalforest/tile.json \
    `# Set to public read access` \
    --acl public-read \
    --content-type application/json