Introductie

Geospatiale data visualiseren in Python

Mary van Valkenburg

Data Science Program Manager, Nashville Software School

Locatie

Cholera-uitbraak in Londen in 1854
600+ doden

miasmawolk

Snows stipkaart

stipkaart van choleradoden in 1864

https://myweb.rollins.edu/jsiry/London-Cholera-John_Snow-copy.jpg

Wat je in deze cursus leert

Geografische punten plotten als scatterplots
Geometrieën plotten met geopandas
Een GeoDataFrame maken uit een DataFrame
Gegevenssets ruimtelijk joinen
Een straatkaart aan je plots toevoegen
Wanneer en hoe je een choropleth maakt

Lengte- en breedtegraad

spreidingsplot van lengte vader en zoon

aarde met lengte- en breedtegraad-rasterlijnen

plt.scatter(schools.Longitude, 
            schools.Latitude, 
            c = 'darkgreen', 
            marker = 'p')
plt.show()

eenvoudige scatterplot van schoollocaties

plt.scatter(schools.Longitude, schools.Latitude, 
            c = 'darkgreen', marker = 'p')
plt.xlabel('Longitude')
plt.ylabel('Latitude')
plt.title('Nashville Public Schools')
plt.grid()
plt.show()

scatterplot van schoollocaties met labels, titel, rasterlijnen

Lengte- en breedtegraad extraheren

bus_stops.head()

Stop ID    StopName       Location
4431      MCC5_11    (36.16659, -86.781996)
588       CHA7AWN    (36.165, -86.78406)
590       CHA8AWN    (36.164393, -86.785451)
541       CXONGULC   (36.162249, -86.790464)
5231      7AVUNISM   (36.163822, -86.783791)

Lengte- en breedtegraad extraheren

bus_stops['lat'] = [loc[0] for loc in bus_stops.Location]
bus_stops['lng'] = [loc[1] for loc in bus_stops.Location]
bus_stops.head()

Stop ID    StopName  Location                  lat         lng
4431      MCC5_11   (36.16659, -86.781996)    36.16659    -86.781996    
588       CHA7AWN  (36.165, -86.78406)        36.165      -86.78406
590       CHA8AWN  (36.164393, -86.785451)    36.164393   -86.785451
541       CXONGULC  (36.162249, -86.790464)   36.162249   -86.790464
5231      7AVUNISM  (36.163822, -86.783791)   36.163822   -86.783791

lng en lat extraheren met reguliere expressies

bus_stops2.head()

Stop ID        Location            
4431       MCC - BAY 11\nNashville, TN\n(36.16659, -86.78199)
588        CHARLOTTE AVE\nNashville, TN\n(36.165, -86.78406)
590        CHARLOTTE AV\nNashville, TN\n(36.164393, -86.785451)
541        CHARLOTTE\nNashville, TN\n(36.162249, -86.790464)
5231       Nashville, TN\n(36.163822, -86.783791)

lng en lat extraheren met reguliere expressies

lat_lng_pattern = re.compile(r'\((.*),\s*(.*)\)', flags=re.MULTILINE)
def extract_lat_lng(address):
    try:
        lat_lng_match = lat_lng_pattern.search(address)
        lat = float(lat_lng_match.group(1))
        lng = float(lat_lng_match.group(2))
        return (lat, lng)
    except:
        return (np.NaN, np.NaN)

lat_lngs = [extract_lat_lng(location)for location in \
           bus_stops2.loc[:, 'Location']]
bus_stops2['lat'] = [lat for lat, lng in lat_lngs]
bus_stops2['lng'] = [lng for lat, lng in lat_lngs]

Open data van Nashville

schermafbeelding van nashville.data.gov

reuzenkippenkop

Laten we oefenen!

Geospatiale data visualiseren in Python