import folium
import google.datalab.bigquery as bq

bq_data = bq.Query(“””
select
 status.hour,
 stations.station_id,
 stations.latitude,
 stations.longitude,
 status.bikes_available
FROM `bigquery-public-data.san_francisco.bikeshare_stations` stations
JOIN (
select
 station_id,
 extract(hour from time) as hour,
 extract(dayofweek from time) as dow,
 sum(if(bikes_available > 2, 1, 0)) as bikes_available
 from `bigquery-public-data.san_francisco.bikeshare_status`
 group by station_id, hour
 having hour=18
) status on stations.station_id = status.station_id
“””)
df = bq_data.execute(output_options=bq.QueryOutput.dataframe()).result()

max_amount = float(df['bikes_available'].max())
hmap = folium.Map(location=[37.333988, -121.894902], zoom_start=9)
hm_wide = HeatMap(zip(df.latitude.values, df.longitude.values, df.bikes_available.values),
                  min_opacity=0.2,
                  max_val = max_amount,
                  radius=17, blur=15,
                  max_zoom=1,
                 )
hmap.add_child(hm_wide)