def extract_hmaps(dstore, what):
"""
Extracts hazard maps. Use it as /extract/hmaps?imt=PGA
"""
info = get_info(dstore)
if what == '': # npz exports for QGIS
sitecol = dstore['sitecol']
mesh = get_mesh(sitecol, complete=False)
dic = _get_dict(dstore, 'hmaps-stats',
{imt: info['poes']
for imt in info['imtls']}, info['stats'])
yield from hazard_items(dic,
mesh,
investigation_time=info['investigation_time'])
return
params = parse(what, info)
if 'imt' in params:
[imt] = params['imt']
m = info['imt'][imt]
s = slice(m, m + 1)
else:
s = ALL
if params['rlzs']:
dset = dstore['hmaps-rlzs']
for k in params['k']:
yield 'rlz-%03d' % k, hdf5.extract(dset, ALL, k, s, ALL)[:, 0]
else:
dset = dstore['hmaps-stats']
stats = list(info['stats'])
for k in params['k']:
yield stats[k], hdf5.extract(dset, ALL, k, s, ALL)[:, 0]
yield from params.items()
def extract_hmaps(dstore, what):
"""
Extracts hazard maps. Use it as /extract/hmaps?imt=PGA
"""
info = get_info(dstore)
if what == '': # npz exports for QGIS
sitecol = dstore['sitecol']
mesh = get_mesh(sitecol, complete=False)
dic = _get_dict(dstore, 'hmaps-stats',
{imt: info['poes'] for imt in info['imtls']},
info['stats'])
yield from hazard_items(
dic, mesh, investigation_time=info['investigation_time'])
return
params = parse(what, info)
if 'imt' in params:
[imt] = params['imt']
m = info['imt'][imt]
s = slice(m, m + 1)
else:
s = ALL
if params['rlzs']:
dset = dstore['hmaps-rlzs']
for k in params['k']:
yield 'rlz-%03d' % k, hdf5.extract(dset, ALL, k, s, ALL)[:, 0]
else:
dset = dstore['hmaps-stats']
stats = list(info['stats'])
for k in params['k']:
yield stats[k], hdf5.extract(dset, ALL, k, s, ALL)[:, 0]
yield from params.items()
def extract_hcurves(dstore, what):
"""
Extracts hazard curves. Use it as /extract/hcurves?kind=mean or
/extract/hcurves?kind=rlz-0, /extract/hcurves?kind=stats,
/extract/hcurves?kind=rlzs etc
"""
info = get_info(dstore)
if what == '': # npz exports for QGIS
sitecol = dstore['sitecol']
mesh = get_mesh(sitecol, complete=False)
dic = _get_dict(dstore, 'hcurves-stats', info['imtls'], info['stats'])
yield from hazard_items(dic,
mesh,
investigation_time=info['investigation_time'])
return
params = parse(what, info)
if 'imt' in params:
[imt] = params['imt']
slc = info['imtls'](imt)
else:
slc = ALL
sids = params.get('site_id', ALL)
if params['rlzs']:
dset = dstore['hcurves-rlzs']
for k in params['k']:
yield 'rlz-%03d' % k, hdf5.extract(dset, sids, k, slc)[:, 0]
else:
dset = dstore['hcurves-stats']
stats = list(info['stats'])
for k in params['k']:
yield stats[k], hdf5.extract(dset, sids, k, slc)[:, 0]
yield from params.items()
def extract_hcurves(dstore, what):
"""
Extracts hazard curves. Use it as /extract/hcurves?kind=mean or
/extract/hcurves?kind=rlz-0, /extract/hcurves?kind=stats,
/extract/hcurves?kind=rlzs etc
"""
info = get_info(dstore)
if what == '': # npz exports for QGIS
sitecol = dstore['sitecol']
mesh = get_mesh(sitecol, complete=False)
dic = _get_dict(dstore, 'hcurves-stats', info['imtls'], info['stats'])
yield from hazard_items(
dic, mesh, investigation_time=info['investigation_time'])
return
params = parse(what, info)
if 'imt' in params:
[imt] = params['imt']
slc = info['imtls'](imt)
else:
slc = ALL
sids = params.get('site_id', ALL)
if params['rlzs']:
dset = dstore['hcurves-rlzs']
for k in params['k']:
yield 'rlz-%03d' % k, hdf5.extract(dset, sids, k, slc)[:, 0]
else:
dset = dstore['hcurves-stats']
stats = list(info['stats'])
for k in params['k']:
yield stats[k], hdf5.extract(dset, sids, k, slc)[:, 0]
yield from params.items()
def extract_hmaps(dstore, what):
"""
Extracts hazard maps. Use it as /extract/hmaps?imt=PGA
"""
oq = dstore['oqparam']
stats = oq.hazard_stats()
num_rlzs = len(dstore['weights'])
if what == '': # npz exports for QGIS
sitecol = dstore['sitecol']
mesh = get_mesh(sitecol, complete=False)
dic = _get_dict(dstore, 'hmaps-stats',
{imt: oq.poes
for imt in oq.imtls}, stats)
yield from hazard_items(dic,
mesh,
investigation_time=oq.investigation_time)
return
params = parse(what, stats, num_rlzs)
if 'imt' in params:
[imt] = params['imt']
m = list(oq.imtls).index(imt)
s = slice(m, m + 1)
else:
s = ALL
for k, i in params['kind'].items():
if k.startswith('rlz-'):
yield k, hdf5.extract(dstore['hmaps-rlzs'], ALL, i, s, ALL)[:, 0]
else:
yield k, hdf5.extract(dstore['hmaps-stats'], ALL, i, s, ALL)[:, 0]
yield from params.items()
def extract_hcurves(dstore, what):
"""
Extracts hazard curves. Use it as /extract/hcurves?kind=mean or
/extract/hcurves?kind=rlz-0, /extract/hcurves?kind=stats,
/extract/hcurves?kind=rlzs etc
"""
oq = dstore['oqparam']
num_rlzs = len(dstore['weights'])
stats = oq.hazard_stats()
if what == '': # npz exports for QGIS
sitecol = dstore['sitecol']
mesh = get_mesh(sitecol, complete=False)
dic = _get_dict(dstore, 'hcurves-stats', oq.imtls, stats)
yield from hazard_items(dic,
mesh,
investigation_time=oq.investigation_time)
return
params = parse(what, stats, num_rlzs)
if 'imt' in params:
[imt] = params['imt']
slc = oq.imtls(imt)
else:
slc = ALL
sids = params.get('site_id', ALL)
for k, i in params['kind'].items():
if k.startswith('rlz-'):
yield k, hdf5.extract(dstore['hcurves-rlzs'], sids, i, slc)[:, 0]
else:
yield k, hdf5.extract(dstore['hcurves-stats'], sids, i, slc)[:, 0]
yield from params.items()
def extract_uhs(dstore, what):
"""
Extracts uniform hazard spectra. Use it as /extract/uhs?kind=mean or
/extract/uhs?kind=rlz-0, etc
"""
info = get_info(dstore)
if what == '': # npz exports for QGIS
sitecol = dstore['sitecol']
mesh = get_mesh(sitecol, complete=False)
dic = {}
for stat, s in info['stats'].items():
hmap = dstore['hmaps-stats'][:, s]
dic[stat] = calc.make_uhs(hmap, info)
yield from hazard_items(dic,
mesh,
investigation_time=info['investigation_time'])
return
params = parse(what, info)
periods = []
for m, imt in enumerate(info['imtls']):
if imt == 'PGA' or imt.startswith('SA'):
periods.append(m)
if 'site_id' in params:
sids = params['site_id']
else:
sids = ALL
if params['rlzs']:
dset = dstore['hmaps-rlzs']
for k in params['k']:
yield ('rlz-%03d' % k, hdf5.extract(dset, sids, k, periods,
ALL)[:, 0])
else:
dset = dstore['hmaps-stats']
stats = list(info['stats'])
for k in params['k']:
yield stats[k], hdf5.extract(dset, sids, k, periods, ALL)[:, 0]
yield from params.items()
def extract_uhs(dstore, what):
"""
Extracts uniform hazard spectra. Use it as /extract/uhs?kind=mean or
/extract/uhs?kind=rlz-0, etc
"""
info = get_info(dstore)
if what == '': # npz exports for QGIS
sitecol = dstore['sitecol']
mesh = get_mesh(sitecol, complete=False)
dic = {}
for stat, s in info['stats'].items():
hmap = dstore['hmaps-stats'][:, s]
dic[stat] = calc.make_uhs(hmap, info)
yield from hazard_items(
dic, mesh, investigation_time=info['investigation_time'])
return
params = parse(what, info)
periods = []
for m, imt in enumerate(info['imtls']):
if imt == 'PGA' or imt.startswith('SA'):
periods.append(m)
if 'site_id' in params:
sids = params['site_id']
else:
sids = ALL
if params['rlzs']:
dset = dstore['hmaps-rlzs']
for k in params['k']:
yield ('rlz-%03d' % k,
hdf5.extract(dset, sids, k, periods, ALL)[:, 0])
else:
dset = dstore['hmaps-stats']
stats = list(info['stats'])
for k in params['k']:
yield stats[k], hdf5.extract(dset, sids, k, periods, ALL)[:, 0]
yield from params.items()
def extract_uhs(dstore, what):
"""
Extracts uniform hazard spectra. Use it as /extract/uhs?kind=mean or
/extract/uhs?kind=rlz-0, etc
"""
oq = dstore['oqparam']
num_rlzs = len(dstore['weights'])
stats = oq.hazard_stats()
if what == '': # npz exports for QGIS
sitecol = dstore['sitecol']
mesh = get_mesh(sitecol, complete=False)
dic = {}
for s, stat in enumerate(stats):
hmap = dstore['hmaps-stats'][:, s]
dic[stat] = calc.make_uhs(hmap, oq)
yield from hazard_items(dic,
mesh,
investigation_time=oq.investigation_time)
return
params = parse(what, stats, num_rlzs)
periods = []
for m, imt in enumerate(oq.imtls):
if imt == 'PGA' or imt.startswith('SA'):
periods.append(m)
if 'site_id' in params:
sids = params['site_id']
else:
sids = ALL
for k, i in params['kind'].items():
if k.startswith('rlz-'):
yield k, hdf5.extract(dstore['hmaps-rlzs'], sids, i, periods,
ALL)[:, 0]
else:
yield k, hdf5.extract(dstore['hmaps-stats'], sids, i, periods,
ALL)[:, 0]
yield from params.items()