In drastega/specklestar: Reduction of Speckle Data from BTA 6-m Telescope

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

Workflow for reduction of speckle interferometric binary stars

Reduction of speckle images of binary and multiple stars using algorithm described in 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)

knitr::include_graphics("Power_spectrum.png")

drastega/specklestar documentation built on Nov. 20, 2023, 8:21 p.m.