names = ['Bias', 'MSE', '$r^2$'] with open('spec_atten_stats.tex', 'w') as fp: for lam_text in wavelengths: for pol,abbr in [('Horizontal','H'), ('Vertical', 'V'), ('Differential', None)]: for ind,name in enumerate(names): fp.write('''\\begin{{frame}} \\frametitle{{Summary: {name} at {lam}-band for {pol} Polarization}} \\begin{{center}} \\begin{{tabular}}{{| c | c | c | c | c |}} \\hline Experiment & Linear & ZPHI & SC & MSC \\\\ \\hline \\hline\n'''.format(name=name, lam=lam_text, pol=pol)) for exp in exps: fp.write(' %s & %.4f & %.4f & %.4f & %.4f \\\\\n'%(exp, stats[lam_text, exp, abbr, 'Linear'][ind], stats[lam_text, exp, abbr, 'ZPHI'][ind], stats[lam_text, exp, abbr, 'SC'][ind], stats[lam_text, exp, abbr, 'MSC'][ind])) fp.write(r''' \hline \end{tabular} \end{center} \end{frame} ''') if __name__ == '__main__': data_cache = load_model_experiments('ref_runs/') make_figs(data_cache)
# Generate PPI plots of Z and Attenuation for (calc, zphi, sc) import matplotlib.pyplot as plt from matplotlib import rcParams from disser.experiments import (script_args, load_model_experiments, process_all_atten) from disser.atten import attenAlgs from disser.datatypes import Reflectivity, Attenuation, PlotInfoContext, PhiDP from disser.plots.defaults import multipanel_cbar_each rcParams['figure.subplot.top'] = 0.95 rcParams['figure.subplot.bottom'] = 0.07 args = script_args().parse_args() data_cache = load_model_experiments('ref_runs', glob='Cband*') data = data_cache[('C','Control')] data.fields.default_keys['source'] = 'average' attenAlgs.run(['ZPHI', 'SC'], data, var='H') moms = [(Reflectivity, 'average'), (PhiDP, 'average'), (Attenuation, 'calc'), (Attenuation, 'ZPHI'), (Attenuation, 'SC')] #moms = [(Attenuation, 'SC'), (Attenuation, 'ZPHI')] data.fields.default_keys['pol'] = 'H' with PlotInfoContext(wavelength=data.wavelength): for dt,src in moms: fig = plt.figure(figsize=(8, 6)) mom = data.fields.grab(dt, source=src) multipanel_cbar_each(fig, (1, 1), [mom], data) if args.save: fig.savefig('{.name}_{:s}.png'.format(dt, src))