def create_png(url, outfile):
quakes=get_eartquake_data(url)
print(quakes[0].__dict__)
#basemap
mpl.rcParams["figure.figsize"]="20,16"
m = Basemap(projection="kav7", lon_0=-90,resolution="l", area_thresh=1000.0)
m.drawcoastlines()
m.drawcountries()
m.mapboundary(fill_color="0.3")
m.drawparallels(np.arange(-90.,99.,30.))
junk = m.drawmeridians(np.arange(-180.,180.,60))
start_day = quakes[-1].timestamp[:10]
end_day = quakes[0].timestamp[:10]
quakes.sort(key=lambda x:x.magnitude,reverse=True)
for q in quakes:
x,y =m(q.lon,q.lat)
mcolor,msize = get_marker(q.magnitude)
m.plot(x,y,mcolor,msize)
plt.title("Earthquakes {0} to {1}".format(start_day, end_day))
plt.savefig(outfile)
def Map(lon,
lat,
dpi=300,
projection='cyl',
resolution='l',
continent=True,
geoCoordinates=True,
mapboundary=False):
m = Basemap(projection=projection,
llcrnrlon=lon.min(),
llcrnrlat=lat.min(),
urcrnrlon=lon.max(),
urcrnrlat=lat.max(),
resolution=resolution)
m.lat = lat
m.lon = lon
m.x, m.y = m(lon, lat)
m.dpi = dpi
m.cs = None
m.continent = continent
m.geoCoordinates = geoCoordinates
m.colorbar = False
m.mapboundary = mapboundary
# create figure
plt.close()
figureSize = 3, 3
m.fig = plt.figure(num=1, figsize=figureSize, dpi=dpi)
return m
def Map(lon, lat, dpi=300, projection='cyl', resolution='l',
continent=True, geoCoordinates=True, mapboundary = False):
m = Basemap(projection=projection, llcrnrlon=lon.min(),
llcrnrlat=lat.min(), urcrnrlon=lon.max(),
urcrnrlat=lat.max(), resolution=resolution)
m.lat = lat
m.lon = lon
m.x, m.y = m(lon, lat)
m.dpi = dpi
m.cs = None
m.continent = continent
m.geoCoordinates = geoCoordinates
m.colorbar = False
m.mapboundary = mapboundary
# create figure
plt.close()
figureSize = 3,3
m.fig = plt.figure(num=1, figsize=figureSize, dpi=dpi)
return m
