earthquakes.csv

glue_world_map.py

from glue import custom_viewer
from glue.core.subset import RoiSubsetState

from mpl_toolkits.basemap import Basemap

GREEN = (229/255., 245/255., 224/255.)
BLUE = (222/255., 235/255., 247/255.)

world_map = custom_viewer('World Map',
                          longitude='att(longitude)',
                          latitude='att(latitude)',
                          magnitude='att(mag)')

@world_map.make_selector
def make_selector(roi, longitude=None, latitude=None, **data):

    state = RoiSubsetState()
    state.roi = roi
    state.xatt = longitude.id
    state.yatt = latitude.id

    return state


@world_map.plot_data
def show_data(ax, longitude=None, latitude=None, magnitude=None, style=None, **kwargs):
    x, y = ax.m(longitude.values, latitude.values)
    return ax.m.scatter(x, y, c=style.color, s=magnitude.values**3, alpha=0.5, edgecolor='none')


@world_map.plot_subset
def show_subset(ax, longitude=None, latitude=None, magnitude=None, style=None, **kwargs):
    x, y = ax.m(longitude.values, latitude.values)
    return ax.m.scatter(x, y, c=style.color, s=magnitude.values**3, alpha=1.0, edgecolor='none')


@world_map.setup
def show_world(ax):

    m = Basemap(projection='mill',llcrnrlat=-80,urcrnrlat=80,\
                llcrnrlon=-180,urcrnrlon=180,lat_ts=20,resolution='c', ax=ax)
    m.drawcoastlines()
    m.fillcontinents(color=GREEN,lake_color=BLUE)
    m.drawmapboundary(fill_color=BLUE)

    ax.m = m

glue_world_map_2.py

from glue import custom_viewer
from glue.core.subset import RoiSubsetState

GREEN = (229 / 255., 245 / 255., 224 / 255.)
BLUE = (222 / 255., 235 / 255., 247 / 255.)

world_map = custom_viewer('World Map',
                          longitude='att(longitude)',
                          latitude='att(latitude)',
                          magnitude='att(mag)')


@world_map.make_selector
def make_selector(state, roi):

    sub = RoiSubsetState()
    sub.roi = roi
    sub.xatt = state.longitude.id
    sub.yatt = state.latitude.id

    return sub


@world_map.plot_data
def show_data(state):
    m = state.m
    longitude = state.longitude
    latitude = state.latitude
    magnitude = state.magnitude
    style = state.style

    x, y = m(longitude.values, latitude.values)
    return m.scatter(x, y, c=style.color, s=magnitude.values ** 3, alpha=0.5, edgecolor='none')


@world_map.plot_subset
def show_subset(state):
    m = state.m
    longitude = state.longitude
    latitude = state.latitude
    magnitude = state.magnitude
    style = state.style

    x, y = m(longitude.values, latitude.values)
    return m.scatter(x, y, c=style.color, s=magnitude.values ** 3, alpha=1.0, edgecolor='none')


@world_map.setup
def show_world(state, ax):

    from mpl_toolkits.basemap import Basemap
    m = Basemap(projection='mill', llcrnrlat=-80, urcrnrlat=80,
                llcrnrlon=-180, urcrnrlon=180, lat_ts=20, resolution='c', ax=ax)
    m.drawcoastlines()
    m.fillcontinents(color=GREEN, lake_color=BLUE)
    m.drawmapboundary(fill_color=BLUE)

    state.m = m


class WorldMap(CustomViewer):

    longitude = 'att(longitude)'
    latitude = 'att(latitude)'
    magnitude = 'att(mag)'

    def setup(self):

        ax = self.axes

        from mpl_toolkits.basemap import Basemap
        m = Basemap(projection='mill', llcrnrlat=-80, urcrnrlat=80,
                    llcrnrlon=-180, urcrnrlon=180, lat_ts=20, resolution='c', ax=ax)
        m.drawcoastlines()
        m.fillcontinents(color=GREEN, lake_color=BLUE)
        m.drawmapboundary(fill_color=BLUE)
        self.m = m

    def plot_subset(self):
        m = self.m
        longitude = self.longitude
        latitude = self.latitude
        magnitude = self.magnitude
        style = self.style

        x, y = m(longitude.values, latitude.values)
        return m.scatter(x, y, c=style.color, s=magnitude.values ** 3, alpha=1.0, edgecolor='none')

    def plot_data(self):
        m = self.m
        longitude = self.longitude
        latitude = self.latitude
        magnitude = self.magnitude
        style = self.style

        x, y = m(longitude.values, latitude.values)
        return m.scatter(x, y, c=style.color, s=magnitude.values ** 3, alpha=0.5, edgecolor='none')

    def make_selector(self, roi):

        state = RoiSubsetState()
        state.roi = roi
        state.xatt = self.longitude.id
        state.yatt = self.latitude.id

        return state

glue_world_map_3.py

from glue import custom_viewer
from glue.core.subset import RoiSubsetState

GREEN = (229 / 255., 245 / 255., 224 / 255.)
BLUE = (222 / 255., 235 / 255., 247 / 255.)

world_map = custom_viewer('World Map',
                          longitude='att(longitude)',
                          latitude='att(latitude)',
                          magnitude='att(mag)')


@world_map.make_selector
def make_selector(roi, longitude, latitude):

    state = RoiSubsetState()
    state.roi = roi
    state.xatt = longitude.id
    state.yatt = latitude.id

    return state


@world_map.plot_data
def show_data(longitude, latitude, magnitude, style, state):
    x, y = state.m(longitude.values, latitude.values)
    return state.m.scatter(x, y, c=style.color, s=magnitude.values ** 3, alpha=0.5, edgecolor='none')


@world_map.plot_subset
def show_subset(longitude, latitude, magnitude, style, state):
    x, y = state.m(longitude.values, latitude.values)
    return state.m.scatter(x, y, c=style.color, s=magnitude.values ** 3, alpha=1.0, edgecolor='none')


@world_map.setup
def show_world(axes, state):
    from mpl_toolkits.basemap import Basemap

    m = Basemap(projection='mill', llcrnrlat=-80, urcrnrlat=80,
                llcrnrlon=-180, urcrnrlon=180, lat_ts=20, resolution='c', ax=axes)
    m.drawcoastlines()
    m.fillcontinents(color=GREEN, lake_color=BLUE)
    m.drawmapboundary(fill_color=BLUE)

    state.m = m


class WorldMap(CustomViewer):

    longitude = 'att(longitude)'
    latitude = 'att(latitude)'
    magnitude = 'att(mag)'

    def show_world(self, axes):
        from mpl_toolkits.basemap import Basemap

        m = Basemap(projection='mill', llcrnrlat=-80, urcrnrlat=80,
                    llcrnrlon=-180, urcrnrlon=180, lat_ts=20, resolution='c', ax=axes)
        m.drawcoastlines()
        m.fillcontinents(color=GREEN, lake_color=BLUE)
        m.drawmapboundary(fill_color=BLUE)

        self.m = m

    def plot_subset(self, longitude, latitude, magnitude, style):
        x, y = self.m(longitude.values, latitude.values)
        return self.m.scatter(x, y, c=style.color, s=magnitude.values ** 3, alpha=1.0, edgecolor='none')

    def plot_data(self, longitude, latitude, magnitude, style):
        x, y = self.m(longitude.values, latitude.values)
        return self.m.scatter(x, y, c=style.color, s=magnitude.values ** 3, alpha=0.5, edgecolor='none')

    def make_selector(self, roi, longitude, latitude):

        state = RoiSubsetState()
        state.roi = roi
        state.xatt = longitude.id
        state.yatt = latitude.id

        return state