histoDemo.py

import opals
from opals import Histo, Import
 
imp=Import.Import(inFile="fullwave.fwf").run()
his=Histo.Histo()
his.commons.screenLogLevel = opals.Types.LogLevel.error
his.inFile = ["fullwave.odm"]
his.probabilities = [0.01,0.05,0.10,0.90,0.95,0.99]
his.sampleRange=[280,320]
his.binWidth=5
his.plotFile="off"  # do not create a svg file
his.run()
stats = his.histogram[0]
print("Selected height statistics of dataset: {}".format(his.inFile[0]))
print("Min   : {:0.2f}".format(stats.getMin()))
print("Median: {:0.2f}".format(stats.getMedian()))
print("Max   : {:0.2f}".format(stats.getMax()))
print("q0.95 : {:0.2f}".format(stats.getQuantiles()[4][1]))
print("\nHistogram: lower height bound [m] - abs. #samples : rel. #samples [%]")
bAbs, bRel = stats.getAbsBins(),stats.getBins()
for i in range(len(bAbs)):
    print("  {:0.1f} - {:d}: {:0.2f}".format(bAbs[i][0], bAbs[i][1], bRel[i][1]*100.))