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