def do_it():
scc_session = Session(default_config, 'scopesim_cluster')
scc_session.do_it_all()
return scc_session
def do_it():
session_lpc = Session(lowpass_config, image_name_lpc)
session_lpc.correct_scopesim_flux_hack()
#session_lpc.find_stars().select_epsfstars_auto().make_epsf()
session_lpc.do_it_all()
return session_lpc
figa = plots.plot_xy_deviation(result_table_multi_conv)
figb = plots.plot_xy_deviation(result_table_multi_subpixel)
# %%
lowpass_config.oversampling = 2
lowpass_config.cutout_size = 40
lowpass_config.fitshape = 25
#session_std = Session(lowpass_config, 'scopesim_grid_16_perturb2_mag18_24')
#session_known_psf = Session(lowpass_config, 'scopesim_grid_16_perturb2_mag18_24')
session_std = Session(lowpass_config, 'gausscluster_N2000_mag22_subpixel')
session_known_psf = Session(lowpass_config,
'gausscluster_N2000_mag22_subpixel')
session_std.do_it_all()
session_known_psf.epsf = make_anisocado_model()
session_known_psf.do_astrometry()
# %%
figa = plots.plot_xy_deviation(session_std.tables.valid_result_table)
figb = plots.plot_xy_deviation(session_known_psf.tables.valid_result_table)
# %%
fig, axs = plt.subplots(1, 2)
axs[0].imshow(session_known_psf.epsf.data, norm=LogNorm())
axs[1].imshow(session_std.epsf.data, norm=LogNorm())
# %%