CLIMADA中重现期输出与影响超越曲线不匹配问题咨询
CLIMADA中重现期影响与超越频率曲线不匹配问题排查求助
我在CLIMADA中针对单一暴露线计算重现期影响,已得到对应输出,但针对同一数据集生成的**超越频率曲线(Exceedance frequency curve)**与重现期数据完全不匹配。尝试调整重现期输入范围后,问题仍未解决。
以下是用于输出重现期的代码:
def local_exceedance_imp(self, return_periods=(25, 50, 100, 250)): """Compute exceedance impact map for given return periods. Requires attribute imp_mat. Parameters ---------- return_periods : Any, optional return periods to consider Dafault is (25, 50, 100, 250) Returns ------- np.array """ LOGGER.info('Computing exceedance impact map for return periods: %s', return_periods) if self.imp_mat.size == 0: raise ValueError('Attribute imp_mat is empty. Recalculate Impact' 'instance with parameter save_mat=True') num_cen = self.imp_mat.shape[1] imp_stats = np.zeros((len(return_periods), num_cen)) cen_step = CONFIG.max_matrix_size.int() // self.imp_mat.shape[0] if not cen_step: raise ValueError('Increase max_matrix_size configuration parameter to > ' f'{self.imp_mat.shape[0]}') # separte in chunks chk = -1 for chk in range(int(num_cen / cen_step)): self._loc_return_imp(np.array(return_periods), self.imp_mat[:, chk * cen_step:(chk + 1) * cen_step].toarray(), imp_stats[:, chk * cen_step:(chk + 1) * cen_step]) self._loc_return_imp(np.array(return_periods), self.imp_mat[:, (chk + 1) * cen_step:].toarray(), imp_stats[:, (chk + 1) * cen_step:]) return imp_stats
请求协助排查该不匹配问题的原因。
内容的提问来源于stack exchange,提问作者CLIMAGUI
相关产品推荐
相关产品推荐

