colorbar_style = lbm.ColorbarStyle(title="Rake",
format="%s",
ticks=None,
spacing="proportional")
thematic_color = lbm.ThematicStyleRanges(
[-180, -135, -45, 45, 135, 180],
['yellow', "green", "yellow", "red", "yellow"],
value_key="rake",
colorbar_style=colorbar_style)
#thematic_color = lbm.ThematicStyleGradient([-180, -90, 0, 90, 180], ['yellow', "green", "yellow", "red", "yellow"], labels=["RL", "Nf", "LL", "Tf", "RL"], value_key="rake", colorbar_style=colorbar_style)
thematic_legend_style = None
focmec_style = lbm.FocmecStyle(size=thematic_size,
fill_color=thematic_color,
thematic_legend_style=thematic_legend_style)
layer = lbm.MapLayer(focmec_data, focmec_style)
layers.append(layer)
legend_style = lbm.LegendStyle(location=1)
title_style = lbm.DefaultTitleTextStyle
title_style.weight = "bold"
map = lbm.LayeredBasemap(layers,
title,
projection,
region=region,
title_style=title_style,
graticule_interval=graticule_interval,
resolution=resolution,
legend_style=legend_style)
map.plot(fig_filespec=fig_filespec)
#layers.append(layer)
## ROB earthquake catalog
catalog = seismodb.query_local_eq_catalog(region=region, start_date=1350)
values = {}
values['magnitude'] = catalog.get_magnitudes()
values['depth'] = catalog.get_depths()
values['year'] = [eq.year for eq in catalog]
point_data = lbm.MultiPointData(catalog.get_longitudes(),
catalog.get_latitudes(),
values=values)
legend_label_style = lbm.TextStyle(font_size=10)
thematic_legend_style = lbm.LegendStyle(title="Magnitude",
location=3,
shadow=True,
label_style=legend_label_style,
fancy_box=True,
label_spacing=0.7)
thematic_size = lbm.ThematicStyleGradient([1, 3, 5, 6.3], [1, 4, 10, 15],
labels=["1.0", "3.0", "5.0", "6.3"],
value_key="magnitude")
## Color scale for depth
#colorbar_style = lbm.ColorbarStyle(title="Depth (km)", location="bottom", format="%d", label_size=18, tick_label_size=18)
#thematic_color = lbm.ThematicStyleColormap(value_key="depth")
#thematic_color = lbm.ThematicStyleRanges([0,1,10,25,50], ['red', 'orange', 'yellow', 'green'], value_key="depth", colorbar_style=colorbar_style)
## Color scale for time
## Discontinuous color scale
#cb_title = u"Jaar - Année - Year"
cb_title = "Year"
layers = []
joined_attributes = {}
joined_attributes[attribute] = {
'key': 'ID_ROB',
'values': {rec['id']: rec[attribute]
for rec in db_records}
}
## Commune layer
gis_data = lbm.GisData(gis_filespec, joined_attributes=joined_attributes)
if attribute == "language":
tsi = lbm.ThematicStyleIndividual(["FR", "NL", "DE"], ['r', 'y', 'b'],
value_key=attribute)
thematic_legend_style = lbm.LegendStyle(title=attribute, location=3)
elif attribute == "id_province":
tsi = lbm.ThematicStyleGradient(values=[1, 6, 11],
styles=["r", "g", "b"],
value_key=attribute)
thematic_legend_style = None
colorbar_style = lbm.ColorbarStyle(title=attribute,
ticks=range(1, 11),
spacing="proportional")
tsi.colorbar_style = colorbar_style
elif attribute == "id_main":
values = joined_attributes[attribute]['values'].values()
vmin, vmax = np.min(values), np.max(values)
tsi = lbm.ThematicStyleColormap(value_key=attribute, vmin=vmin, vmax=vmax)
thematic_legend_style = None
colorbar_style = lbm.ColorbarStyle(title=attribute)