Mapping and GIS Related Scripts:
--------------------------------

airlineroutes.py:
    Plot US airline routes on a background map.
    Mainly useful as an illustration of plotting from geopandas.

demraytrace.py:
    Raytrace a Digital Elevation Model file using povray,
    creating a 360-degree panorama.
    demproj.py is similar but pointless, making lots of smaller
    povray projections that you can stitch together.

firedata.py:
    Track progress of a wildfire. Download MODIS and VIIRS hotspot data
    plus perimeter data, and massage it into something pytopo can show.

geoipcompare.py:
    Check several different geoip services and log the varying
    answers they give.

geojsonmap.py:
    How to plot a geojson file of polygons exported from QGIS.

grassviewshed.py:
    Use GRASS r.viewshed for a viewshed analysis on a DEM input file.

greatcircle.py:
    Calculate bearings for a great circle trip between two coordinates
    on a spherical planet.

plotrivers.py:
    Follow the tutorial at
    https://towardsdatascience.com/creating-beautiful-river-maps-with-python-37c9b5f5b74c
    updated to work with current data files, and showing North America.
    Try to fix some matplotlib problems that show up.
    *Very slow.*

polidistmap.py:
    Show political districts -- or any other geojson data -- on a map
    using folium to create a leaflet map.
    Also has examples for using gdal to translate shapefile to geojson.
    (But see also districtmaps in my webapps repo for a more up-to-date
    and maintained solution, which I use on the LWVNM
    https://www.lwvnm.org/districtmaps/ page.)

shapemerge.py:
    Starting from a shapefile of lots of small shapes, merge (dissolve)
    them into a smaller number of shapes. For instance, coalesce
    voting precincts into voting districts each containing several
    precincts. Uses fiona and shapely.

splitgeojson.py:
    Given a geojson that has two tags of interest, split by one of
    those tags. For instance, the NM SOS county commission district
    GIS has all counties in one file, with tags DIST_NUM and COUNTY,
    and it's useful to split by COUNTY into a file per county.

viewshed2view.py:
    Take a viewshed (like one generated by GRASS r.viewshed) and
    turn it into a 360 image showing what's visible.

waymaker.py:
    Take a file and of descriptions, multi-line and separated by blank lines,
    and turn it into a collection of GPX waypoints suitable for import into
    any mapping program (e.g. Osmand or PyTopo).