def _run_qmed_analysis(self):
results = {}
analysis = QmedAnalysis(self.catchment, self.gauged_catchments, results_log=results)
self.qmed = analysis.qmed()
results['qmed'] = self.qmed
self.results['qmed'] = results
def test_qmed_two_donors(self):
analysis = QmedAnalysis(self.catchment, CatchmentCollections(self.db_session), year=2000)
donors = analysis.find_donor_catchments()[0:2] # 17001, 10001
result = analysis.qmed(method='descriptors', donor_catchments=donors)
# exp(ln(0.61732109) + 0.34379622 * 0.55963062 + 0.00102012 * 0.02991561) =
# exp(ln(0.61732109) + 0.192429411) = 0.748311028
self.assertAlmostEqual(result, 0.748311028, places=5)
def _run_qmed_analysis(self):
results = {}
analysis = QmedAnalysis(self.catchment,
self.gauged_catchments,
results_log=results)
self.qmed = analysis.qmed()
results['qmed'] = self.qmed
self.results['qmed'] = results
def analyse_catchment(catchment, gauged_catchments):
result = {}
result['id'] = catchment.id
analysis = QmedAnalysis(catchment)
result['qmed_amax'] = analysis.qmed(method='amax_records')
result['qmed_descr'] = analysis.qmed(method='descriptors')
result['qmed_descr_1999'] = analysis.qmed(method='descriptors_1999')
analysis.gauged_catchments = gauged_catchments
donors = analysis.find_donor_catchments()
analysis.idw_power = 2
result['qmed_descr_idw'] = analysis.qmed(method='descriptors', donor_catchments=donors)
analysis.idw_power = 3
result['qmed_descr_idw3'] = analysis.qmed(method='descriptors', donor_catchments=donors)
analysis.donor_weighting = 'equal'
result['qmed_descr_first'] = analysis.qmed(method='descriptors', donor_catchments=donors[0:1])
return result
def analyse_catchment(catchment, gauged_catchments):
result = {}
result['id'] = catchment.id
analysis = QmedAnalysis(catchment)
result['qmed_amax'] = analysis.qmed(method='amax_records')
result['qmed_descr'] = analysis.qmed(method='descriptors')
result['qmed_descr_1999'] = analysis.qmed(method='descriptors_1999')
analysis.gauged_catchments = gauged_catchments
donors = analysis.find_donor_catchments()
analysis.idw_power = 2
result['qmed_descr_idw'] = analysis.qmed(method='descriptors',
donor_catchments=donors)
analysis.idw_power = 3
result['qmed_descr_idw3'] = analysis.qmed(method='descriptors',
donor_catchments=donors)
analysis.donor_weighting = 'equal'
result['qmed_descr_first'] = analysis.qmed(method='descriptors',
donor_catchments=donors[0:1])
return result
from floodestimation.loaders import load_catchment
from floodestimation import db
from floodestimation.collections import CatchmentCollections
from floodestimation.analysis import QmedAnalysis
db_session = db.Session()
dee_catchment = load_catchment('nith_cds.cd3')
gauged_catchments = CatchmentCollections(db_session)
qmed_analysis = QmedAnalysis(dee_catchment, gauged_catchments)
print(qmed_analysis.qmed())
print(qmed_analysis.methods)
print(qmed_analysis.qmed_all_methods())
print(qmed_analysis.urban_adj_factor())
print(qmed_analysis.find_donor_catchments(5, 200.0))
qmed_analysis.idw_power = 1.5
print(qmed_analysis.idw_power)
donors = qmed_analysis.find_donor_catchments(5, 200.0)
for donor in donors:
Q = QmedAnalysis(donors[0], gauged_catchments)
print(donor, qmed_analysis._error_correlation(donor), Q.qmed_all_methods())
db_session.close()
from floodestimation.loaders import load_catchment
from floodestimation import db
from floodestimation.collections import CatchmentCollections
from floodestimation.analysis import QmedAnalysis
db_session = db.Session()
dee_catchment = load_catchment('nith_cds.cd3')
gauged_catchments = CatchmentCollections(db_session)
qmed_analysis = QmedAnalysis(dee_catchment, gauged_catchments)
print(qmed_analysis.qmed())
print(qmed_analysis.methods)
print(qmed_analysis.qmed_all_methods())
print(qmed_analysis.urban_adj_factor())
print(qmed_analysis.find_donor_catchments(5, 200.0))
qmed_analysis.idw_power = 1.5
print(qmed_analysis.idw_power)
donors = qmed_analysis.find_donor_catchments(5, 200.0)
for donor in donors:
Q = QmedAnalysis(donors[0], gauged_catchments)
print(donor,qmed_analysis._error_correlation(donor),Q.qmed_all_methods())
