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()) # %%