Workflow for reduction of speckle interferometric binary stars
For reduction of speckle images of binary and multiple stars we use algorithm described in paper Pluzhnik (2005).
Calculate statistics for object file
obj_filename <- file.choose() # Choose file with binary star
obj_stat <- speckle_stat(obj_filename)Calculate statistics for dark file
dark_filename <- file.choose() # Choose file with dark frames
dark_stat <- speckle_stat(dark_filename)Calculate statistics for flat file
flat_filename <- file.choose() # Choose file with flat field frames
flat_stat <- speckle_stat(flat_filename)Calculate average dark frame
zero_matrix <- matrix(0, 512, 512)
mid_dark <- middle_frame(filename = dark_filename, subtrahend = zero_matrix)Calculate average flat field frame with average dark subtraction
mid_flat <- middle_frame(flat_filename, subtrahend = mid_dark)Calculate power spectrum with subtraction of neighboring frame
ps_diff <- speckle_ps_diff(obj_filename)Calculate acf
acf <- speckle_acf(ps_diff)Calculate power spectrum
ps <- speckle_ps(obj_filename, dark = mid_dark, flat = mid_flat)
Fig. 1. Power spectrum of star ADS 15182 in 550/20 nm band.