def test_stats_models():
conf_file = os.path.join(upone, 'defaultconfigfiles', 'models',
'zhu_2015.ini')
conf = ConfigObj(conf_file)
data_path = os.path.join(datadir, 'loma_prieta', 'model_inputs')
# Check slopefile trimming
conf['zhu_2015']['slopefile'] = 'global_gted_maxslope_30c.flt'
conf = correct_config_filepaths(data_path, conf)
# Run with divfactor of 1
conf['zhu_2015']['divfactor'] = '1.'
shakefile = os.path.join(datadir, 'loma_prieta', 'grid.xml')
lm = LM.LogisticModel(shakefile, conf, saveinputs=True)
#maplayers1 = lm.calculate()
conf_file = os.path.join(upone, 'defaultconfigfiles', 'models',
'zhu_2017_coastal.ini')
conf = ConfigObj(conf_file)
data_path = os.path.join(datadir, 'loma_prieta', 'model_inputs')
conf['zhu_2017_coastal']['slopefile'] = 'global_gted_maxslope_30c.flt'
conf = correct_config_filepaths(data_path, conf)
# Run with divfactor of 1
conf['zhu_2017_coastal']['divfactor'] = '1.'
shakefile = os.path.join(datadir, 'loma_prieta', 'grid.xml')
lm = LM.LogisticModel(shakefile, conf, saveinputs=True)
maplayers2 = lm.calculate()
# Change shakemap name so that it doesn't stomp on the other
maplayers2['model']['description']['shakemap'] = '19891018000415_ver2'
#model_list = [maplayers1, maplayers2]
#test_dict1 = assess_models.concatenateModels(model_list)
#test_dict2 = assess_models.concatenateModels(model_list, astitle='model')
# I think this test is freezing on travis
# tmp = assess_models.modelSummary(test_dict2, showplots=False,
# summary_figure=False,
# individual_plots=False)
# np.testing.assert_allclose(tmp[0][0], 0.025677016713957716)
# np.testing.assert_allclose(tmp[1][0], 0.00098462898029272805)
hagg = stats.computeHagg(maplayers2['model']['grid'])
np.testing.assert_allclose(hagg['hagg_0.00g'], 65.85571, atol=0.001)
parea = stats.computeParea(maplayers2['model']['grid'], probthresh=0.2)
np.testing.assert_allclose(parea, 77.03092, atol=0.001)
stats2 = stats.computeStats(maplayers2['model']['grid'],
probthresh=0.2,
shakefile=shakefile,
shakethreshtype='pga',
shakethresh=20.,
statprobthresh=0.0)
np.testing.assert_allclose(stats2['Max'], 0.41278, atol=0.001)
np.testing.assert_allclose(stats2['Median'], 0.00032317817, rtol=0.001)
np.testing.assert_allclose(stats2['Std'], 0.04855, atol=0.0001)
np.testing.assert_allclose(stats2['hagg_0.20g'], 55.47086, atol=0.001)
def test_stats_models():
conf_file = os.path.join(upone, 'defaultconfigfiles', 'models',
'zhu_2015.ini')
conf = ConfigObj(conf_file)
data_path = os.path.join(datadir, 'loma_prieta', 'model_inputs')
# Check slopefile trimming
conf['zhu_2015']['slopefile'] = 'global_gted_maxslope_30c.flt'
conf = correct_config_filepaths(data_path, conf)
# Run with divfactor of 1
conf['zhu_2015']['divfactor'] = '1.'
shakefile = os.path.join(datadir, 'loma_prieta', 'grid.xml')
lm = LM.LogisticModel(shakefile, conf, saveinputs=True)
#maplayers1 = lm.calculate()
conf_file = os.path.join(upone, 'defaultconfigfiles', 'models',
'zhu_2017_coastal.ini')
conf = ConfigObj(conf_file)
data_path = os.path.join(datadir, 'loma_prieta', 'model_inputs')
conf['zhu_2017_coastal']['slopefile'] = 'global_gted_maxslope_30c.flt'
conf = correct_config_filepaths(data_path, conf)
# Run with divfactor of 1
conf['zhu_2017_coastal']['divfactor'] = '1.'
shakefile = os.path.join(datadir, 'loma_prieta', 'grid.xml')
lm = LM.LogisticModel(shakefile, conf, saveinputs=True)
maplayers2 = lm.calculate()
# Change shakemap name so that it doesn't stomp on the other
maplayers2['model']['description']['shakemap'] = '19891018000415_ver2'
hagg = stats.computeHagg(maplayers2['model']['grid'])
np.testing.assert_allclose(hagg['hagg_0.00g'], 5.12723506006, atol=0.001)
stats2 = stats.computeStats(maplayers2['model']['grid'],
shakefile=shakefile,
shakethreshtype='pga',
shakethresh=20.,
probthresh=0.0)
np.testing.assert_allclose(stats2['Max'], 0.4105819792026343, atol=0.001)
np.testing.assert_allclose(stats2['Median'],
0.34855563636356035,
rtol=0.01)
np.testing.assert_allclose(stats2['Std'], 0.04494649549605273, atol=0.001)
np.testing.assert_allclose(stats2['hagg_0.20g'],
2.845917180389227,
atol=0.001)