class Analysis(object):
def __init__(self):
self.name = None
self.catchment = Catchment("River Town", "River Burn")
self.db_session = db.Session()
self.gauged_catchments = CatchmentCollections(self.db_session)
self.qmed = None
def finish(self):
self.db_session.close()
def run_qmed_analysis(self):
self.qmed_analysis = QmedAnalysis(self.catchment, self.gauged_catchments)
self.results = self.qmed_analysis.results_log
self.results['qmed_all_methods'] = self.qmed_analysis.qmed_all_methods()
def run_growthcurve(self):
results = {}
analysis = GrowthCurveAnalysis(self.catchment, self.gauged_catchments, results_log=results)
gc = analysis.growth_curve()
aeps = [0.5, 0.2, 0.1, 0.05, 0.03333, 0.02, 0.01333, 0.01, 0.005, 0.002, 0.001]
growth_factors = gc(aeps)
flows = growth_factors * self.qmed
results['aeps'] = aeps
results['growth_factors'] = growth_factors
results['flows'] = flows
self.results['gc'] = results
class Analysis(object):
def __init__(self):
self.name = None
self.catchment = Catchment("River Town", "River Burn")
self.db_session = db.Session()
self.gauged_catchments = CatchmentCollections(self.db_session)
self.qmed = None
def finish(self):
self.db_session.close()
def run_qmed_analysis(self):
self.qmed_analysis = QmedAnalysis(self.catchment,
self.gauged_catchments)
self.results = self.qmed_analysis.results_log
self.results['qmed_all_methods'] = self.qmed_analysis.qmed_all_methods(
)
def run_growthcurve(self):
results = {}
analysis = GrowthCurveAnalysis(self.catchment,
self.gauged_catchments,
results_log=results)
gc = analysis.growth_curve()
aeps = [
0.5, 0.2, 0.1, 0.05, 0.03333, 0.02, 0.01333, 0.01, 0.005, 0.002,
0.001
]
growth_factors = gc(aeps)
flows = growth_factors * self.qmed
results['aeps'] = aeps
results['growth_factors'] = growth_factors
results['flows'] = flows
self.results['gc'] = results
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())
