def read_area_source(area_source_file, discretisation=200.):
"""Calls OpenQuake parsers to read area source model
from source_mode.xml type file, convert to point sources
and write to a new nrml source model file.
:params area_source_file:
nrml format file of the area source
:params discretisation:
Grid size (km) for the area source discretisation,
which defines the distance between resulting point
sources.
"""
converter = SourceConverter(50,
10,
width_of_mfd_bin=0.1,
area_source_discretization=discretisation)
parser = SourceModelParser(converter)
print[method
for method in dir(parser)] # if callable(getattr(parser, method))]
try:
sources = parser.parse_sources(area_source_file)
except AttributeError: # Handle version 2.1 and above
sources = []
groups = parser.parse_src_groups(area_source_file)
for group in groups:
for source in group:
sources.append(source)
for source in sources:
print source.polygon
return sources
def check_round_trip(self, fname):
parser = SourceModelParser(SourceConverter(50., 1., 10, 0.1, 10.))
groups = parser.parse_src_groups(fname)
fd, name = tempfile.mkstemp(suffix='.xml')
with os.fdopen(fd, 'wb'):
write_source_model(name, groups, 'Test Source Model')
if open(name).read() != open(fname).read():
raise Exception('Different files: %s %s' % (name, fname))
os.remove(name)
def area2pt_source(area_source_file, discretisation=200.):
"""Calls OpenQuake parsers to read area source model
from source_mode.xml type file, convert to point sources
and write to a new nrml source model file.
:params area_source_file:
nrml format file of the area source
:params discretisation:
Grid size (km) for the area source discretisation,
which defines the distance between resulting point
sources.
"""
converter = SourceConverter(50,
10,
width_of_mfd_bin=0.1,
area_source_discretization=discretisation)
parser = SourceModelParser(converter)
print[method
for method in dir(parser)] # if callable(getattr(parser, method))]
try:
sources = parser.parse_sources(area_source_file)
except AttributeError: # Handle version 2.1 and above
sources = []
groups = parser.parse_src_groups(area_source_file)
for group in groups:
for source in group:
sources.append(source)
name = 'test_point_model'
new_pt_sources = {}
for source in sources:
pt_sources = area_to_point_sources(source)
for pt in pt_sources:
pt.source_id = pt.source_id.replace(':', '')
pt.name = pt.name.replace(':', '_')
try:
new_pt_sources[pt.tectonic_region_type].append(pt)
except KeyError:
new_pt_sources[pt.tectonic_region_type] = [pt]
# print [method for method in dir(pt) if callable(getattr(pt, method))]
# print [attribute for attribute in dir(pt)]
nrml_pt_file = area_source_file[:-4] + '_pts.xml'
source_group_list = []
id = 0
for trt, sources in new_pt_sources.iteritems():
source_group = SourceGroup(trt, sources=sources, id=id)
id += 1
source_group_list.append(source_group)
write_source_model(nrml_pt_file, source_group_list, name='Leonard2008')
def read_simplefault_source(simplefault_source_file):
"""Read nrml source model into simpmle fault objects
"""
converter = SourceConverter(50, 2, width_of_mfd_bin=0.1,
area_source_discretization=200.)
parser = SourceModelParser(converter)
try:
sources = parser.parse_sources(simplefault_source_file)
except AttributeError: # Handle version 2.1 and above
sources = []
groups = parser.parse_src_groups(simplefault_source_file)
for group in groups:
for source in group:
sources.append(source)
for fault in sources:
print fault.mfd.max_mag
return sources
def read_pt_source(pt_source_file):
"""Read nrml source model into pt source objects
"""
converter = SourceConverter(50,
10,
width_of_mfd_bin=0.1,
area_source_discretization=200.)
parser = SourceModelParser(converter)
try:
sources = parser.parse_sources(pt_source_file)
except AttributeError: # Handle version 2.1 and above
sources = []
groups = parser.parse_src_groups(pt_source_file)
for group in groups:
for source in group:
sources.append(source)
return sources
