def obsolete_plot_hazard_map(indir, datafile, out_dir, output_file, save_file=None,
title=None, np_posn=None, s_posn=None, cb_label=None,
cb_steps=None, invert=False, colourmap=None,
annotate=None):
"""Plot a hazard map from gridded data.
indir input directory
data name of the data file to plot (in indir)
out_dir general output directory
output_file name of map output file to create in 'out_dir' directory
save_file name of map data output file to create in 'out_dir' directory
title title to put on the graph
np_posn place to put a north pointer symbol at
s_posn place to put a length scale at
cb_label label text to put on the colour bar
(if not defined, no colour bar)
cb_steps sequence of required values in colourbar
colourmap the base colourmap to use
annotate user annotations required
"""
if cb_steps is None:
cb_steps = [0.1,0.15,0.20,0.30,0.40,0.50]
data = cld.calc_load_data(datafile, invert=invert)
pgxc.plot_gmt_xyz_contour(data, output_file,
title=title, np_posn=np_posn,
s_posn=s_posn, cb_label=cb_label,
cb_steps=cb_steps, annotate=annotate)
def do_example(datafile, outfile, invert=False):
# ignore_value = 0.025
#
annotations = []
# if ignore_value:
# annotations.append(('text', (143.6,-37.10),
# 'Values <= %.2f are ignored' % ignore_value))
data = cld.calc_load_data(datafile, invert=invert)
# TODO : add call to calc_ignore
pgxc.plot_gmt_xyz_contour(data, outfile,
title='Bedrock Hazard, RP=?, RSA=?',
np_posn='NE', s_posn='SE',
cb_label='Acceleration (g)',
cb_steps=[0.1,0.15,0.20,0.30,0.40,0.50],
annotate=annotations, linewidth=0.0)
def do_example(datafile, outfile, invert=False):
# ignore_value = 0.025
#
annotations = []
# if ignore_value:
# annotations.append(('text', (143.6,-37.10),
# 'Values <= %.2f are ignored' % ignore_value))
data = cld.calc_load_data(datafile, invert=invert)
# TODO : add call to calc_ignore
pgxc.plot_gmt_xyz_contour(data,
outfile,
title='Bedrock Hazard, RP=?, RSA=?',
np_posn='NE',
s_posn='SE',
cb_label='Acceleration (g)',
cb_steps=[0.1, 0.15, 0.20, 0.30, 0.40, 0.50],
annotate=annotations,
linewidth=0.0)