force_recompute_photometry = False # True # Calculate master dark/flat: if not os.path.exists(master_dark_path) or not os.path.exists( master_flat_path): print('Calculating master flat:') generate_master_flat_and_dark(flat_paths, dark_paths, master_flat_path, master_dark_path) # Do photometry: if not os.path.exists(output_path) or force_recompute_photometry: print('Calculating photometry:') phot_results = photometry(image_paths, master_dark_path, master_flat_path, target_centroid, comparison_flux_threshold, aperture_radii, centroid_stamp_half_width, psf_stddev_init, aperture_annulus_radius, output_path) else: phot_results = PhotometryResults.load(output_path) print('Calculating PCA...') light_curve = PCA_light_curve(phot_results, transit_parameters, plots=True) plt.figure() plt.plot(phot_results.times, light_curve, 'k.') #egress = 2457777.01 #post_egress_std = np.std(light_curve[phot_results.times > egress]) #plt.axvline(egress) plt.xlabel('Time [JD]')
force_recompute_photometry = True # Calculate master dark/flat: if not os.path.exists(master_dark_path) or not os.path.exists( master_flat_path): print('Calculating master flat:') generate_master_flat_and_dark(night_flat_paths, dark_30s_paths, master_flat_path, master_dark_path) # Do photometry: if not os.path.exists(output_path) or force_recompute_photometry: print('Calculating photometry:') phot_results = photometry(image_paths, master_dark_path, master_flat_path, star_positions, aperture_radii, centroid_stamp_half_width, psf_stddev_init, aperture_annulus_radius, output_path, brightest_star_coords_init) else: phot_results = PhotometryResults.load(output_path) # print('Calculating PCA...') target = phot_results.fluxes[:, 0, 0] mean_comparison = phot_results.fluxes[:, 1, 0] light_curve = target / mean_comparison plt.plot(phot_results.times, light_curve / np.median(light_curve), '.') np.savetxt('outputs/20160705.txt', np.vstack([phot_results.times, light_curve]).T)
transit_parameters = params_c output_path = 'outputs/trappist1c_20160619.npz' force_recompute_photometry = False #True # Calculate master dark/flat: if not os.path.exists(master_dark_path): generate_master_dark(dark_paths, master_dark_path) # Do photometry: if not os.path.exists(output_path) or force_recompute_photometry: print('Calculating photometry:') phot_results = photometry(image_paths, master_dark_path, master_flat_path, target_centroid, comparison_flux_threshold, aperture_radii, centroid_stamp_half_width, psf_stddev_init, aperture_annulus_radius, output_path) else: phot_results = PhotometryResults.load(output_path) print('Calculating PCA...') light_curve = PCA_light_curve(phot_results, transit_parameters, validation_duration_fraction=0.5, plots=True, plot_validation=False, validation_time=0.9, outlier_rejection=False) plt.figure() plt.plot(phot_results.times, light_curve, 'k.')