| 2013-11-19 15:21:25 |
Damiano Monelli |
bug |
|
|
added bug |
| 2013-12-12 11:07:54 |
Damiano Monelli |
description |
When deriving a hazard map or when doing a disaggregation calculation, the ground motion value associated to a probability of exceedance value, is computed by performing a liner interpolation between the two closest probability of exceedance values (see https://github.com/gem/oq-engine/blob/master/openquake/engine/calculators/hazard/classical/post_processing.py#L77, https://github.com/gem/oq-engine/blob/master/openquake/engine/calculators/hazard/disaggregation/core.py#L137.
For more precise results the interpolation should be done in log space, that is instead of:
numpy.interp(poes, curve[::-1], imls)
we should use:
numpy.interp(poes, numpy.log10(curve[::-1]), numpy.log10(imls)) |
When deriving a hazard map or when doing a disaggregation calculation, the ground motion value associated to a probability of exceedance value, is computed by performing a liner interpolation between the two closest probability of exceedance values (see https://github.com/gem/oq-engine/blob/master/openquake/engine/calculators/hazard/classical/post_processing.py#L77, https://github.com/gem/oq-engine/blob/master/openquake/engine/calculators/hazard/disaggregation/core.py#L137.
For more precise results the interpolation should be done in log space, that is instead of:
numpy.interp(poes, curve[::-1], imls)
we should use:
numpy.exp(numpy.interp(numpy.log10(poes), numpy.log10(curve[::-1]),
numpy.log10(imls))) |
|
| 2013-12-12 11:16:26 |
Damiano Monelli |
description |
When deriving a hazard map or when doing a disaggregation calculation, the ground motion value associated to a probability of exceedance value, is computed by performing a liner interpolation between the two closest probability of exceedance values (see https://github.com/gem/oq-engine/blob/master/openquake/engine/calculators/hazard/classical/post_processing.py#L77, https://github.com/gem/oq-engine/blob/master/openquake/engine/calculators/hazard/disaggregation/core.py#L137.
For more precise results the interpolation should be done in log space, that is instead of:
numpy.interp(poes, curve[::-1], imls)
we should use:
numpy.exp(numpy.interp(numpy.log10(poes), numpy.log10(curve[::-1]),
numpy.log10(imls))) |
When deriving a hazard map or when doing a disaggregation calculation, the ground motion value associated to a probability of exceedance value, is computed by performing a liner interpolation between the two closest probability of exceedance values (see https://github.com/gem/oq-engine/blob/master/openquake/engine/calculators/hazard/classical/post_processing.py#L77, https://github.com/gem/oq-engine/blob/master/openquake/engine/calculators/hazard/disaggregation/core.py#L137.
For more precise results the interpolation should be done in log space, that is instead of:
numpy.interp(poes, curve[::-1], imls)
we should use:
10 ** (numpy.interp(numpy.log10(poes), numpy.log10(curve[::-1]),
numpy.log10(imls))) |
|
| 2013-12-12 16:13:35 |
Michele Simionato |
oq-engine: assignee |
|
Michele Simionato (michele-simionato) |
|
| 2013-12-12 16:13:37 |
Michele Simionato |
oq-engine: milestone |
|
1.0.1 |
|
| 2013-12-12 16:13:41 |
Michele Simionato |
oq-engine: importance |
Undecided |
High |
|
| 2013-12-12 16:13:45 |
Michele Simionato |
oq-engine: status |
New |
Fix Committed |
|
| 2014-12-15 11:55:55 |
Daniele Viganò |
oq-engine: status |
Fix Committed |
Fix Released |
|
