create get_CEJST() function

[ericnost]

Something like this (see Cumulative Impacts notebook):

import zipfile, requests, io

file = 'https://static-data-screeningtool.geoplatform.gov/data-versions/1.0/data/score/downloadable/1.0-shapefile-codebook.zip'

r = requests.get(file)
z = zipfile.ZipFile(io.BytesIO(r.content))
z.extractall("CEJST")
z = zipfile.ZipFile("/content/CEJST/usa.zip")
z.extractall("CEJST")

cejst = geopandas.read_file("/content/CEJST/usa.shp", mask=county)

# interactive widget selection here, to select the CEJST variable to map

map = choropleth(cejst, "DSF_PFS", key_id="GEOID10")

return cejst, map

See EJScreen class in watershed notebook for example of interactive widget selection