Add interstates to pyart radar maps

Python
1

@mgrover1 maybe we need to have this built into pyart since the roads feature provided by cartopy is way to many roads. But here is my hack at getting interstates (e.g., i-90 etc) into my radar maps:

!wget https://www2.census.gov/geo/tiger/TIGER2016/PRIMARYROADS/tl_2016_us_primaryroads.zip
!unzip tl_2016_us_primaryroads.zip

import cartopy.io.shapereader as shpreader
from cartopy.feature import ShapelyFeature
import cartopy.crs as ccrs
import cartopy.feature as cfeature

# Read shape file
reader = shpreader.Reader("./tl_2016_us_primaryroads.shp")
names = []
geoms = []
for rec in reader.records():
    if (rec.attributes['FULLNAME'][0]=='I'):
        names.append(rec)
        geoms.append(rec.geometry)

#make interstate feature
shape_feature = ShapelyFeature(geoms, ccrs.PlateCarree(), edgecolor='darkblue', lw=0.5,facecolor='none')

Then you call it in your plot block like this:

ax.add_feature(shape_feature)

Thoughts?

1 Like
2

Thanks for sharing this! I think adding this into one of the toolkits would be great!! An INTERSTATES feature?