Processing Log

In order to get anything useful out of our images, we have to correct for sytematic errors or "noise". Take for example, this raw image of our galaxy:

Raw Image

This is a raw image (ngc2903c1-50-300-001), with certain noise points enhanced (the red line on the left was added to make the dark overscan region visible on this page). Processing the images can get rid of most of this noise. Image is displayed north down, east left due to ImageJ's tendency to flip fits files.

The same basic steps are applied to all images, including calibration frames. The following is a discription of some of the common processes:

Overscan Subtract

This is the first step applied to all images, including the calibration masters. On the left edge of the raw image is a band of pixels about 6 wide and as long as the image called the overscan region. After the image has been exposed and the CCD chip is being read out, it adds these pixels which contain information on the bias offset at that time. This step includes measuring this region, taking the mean value, and subtracting it from the entire image. Then, starting at x = 8 pixels, cropping the image to 520x512 pixels, thereby removing the overscan region

Registration

This is the process of aligning images so that each area on the image, say a star, is at the same pixel location on each image. There are two plugins for ImageJ (created by Phillipe Thévenaz, available here)that perform this function: TurboReg, which aligns separate images, and StackReg, which aligns all the images in a stack at once.

Averaging

This process simply takes a bunch of images and creates a new image. Each pixel on the new image has the mean value of each corresponding pixel on the images that created it. This has a tendancy of eliminating random noise.

We can get rid of noise by taking calibration frames which we average and apply to the images. The types of frames we take are:

Bias

These are 10 0s exposures with the shutter on the camera closed. Basically it's just a readout of the noise generated by the voltage in the chip. An Overscan Subtract is performed on the frames and then they are averaged to produce the Master Bias.

Darks

These correct for noise caused by chip temperature and "hot pixels". A CCD chip works by collecting electrons released by the photoelectric effect (photons hitting the chip), but electrons can also be released by heat, which adds noise. Hot pixels are "bad" pixels; they read too high on each image. They are exposures with the shutter closed for the same length of time as the images to be calibrated. 5 of these are taken, an Overscan subtract is Master Bias is subtracted from them before averaging to create the Master Dark frame. While there are occasional random extremely hot pixels (ie on the order of hundreds or thousands in a region where the pixel counts are on the order of 1 or tens) these are not caused by the chip and so are not corrected for as they are truly random events.

Flats

These correct for uneven lighting and artifacts (shadows due to dust and stuff on optical surfaces) caused by the telescope-filter-camera optics and position of the guide camera. There are three different methods of taking them, we used two. Dark Sky Flats (feb6/7): take 5-10 images of the sky around the object, each for the same amount of time as the object data, shifting the field each time so that any stars are on different pixels. Twilight Flats (mar15/16): at sunset/sunrise take images of the sky opposite the setting/rising sun at ~ 25° altitude, and expose until the maximum counts hit 20 000. These require their own darks. To process both, an overscan subract is performed on each before subtracting the Master Dark and Master Bias frames. Then the image is normalized (ie every pixel divided) by the mean, and then combined by median to create the Master Flat frame. Combining by median is the same process as averaging, except each pixel of the new image is the median value of the corresponding pixels in the component images.

For a better look at the final images and some of the composite frames click here.

February 6/7 Data

  1. Bias frames converted from 16- bit to 32-bit and Master Bias created as described above
  2. Dark frames converted from 16-bit to 32-bit. Master Dark created as described above
  3. Flat frames converted from 16-bit to 32-bit and Master Dark and Master Bias subtracted from flat frames after overscan subtract. Background extractor plugin used to normalize frames by mean and combine by median into Master Flat frame
  4. Object images imported as a stack and converted from 16-bit to 32-bit
  5. Overscan Subtract performed on each slice
  6. Master Dark and Master Bias frames subtracted from each
  7. Each image divided by Master Flat
  8. Images registered to # 1 using StackReg
  9. Images Averaged

Comments: Noise added by Master Flat division due to the guide camera being moved between our images and flat fields

March 15/16 Clear Data

  1. Bias frames converted from 16- bit to 32-bit and Master Bias created as described above
  2. Dark frames converted from 16-bit to 32-bit. 3s Master Dark and 60s Master Dark created as described above
  3. Flat frames converted from 16-bit to 32-bit and 3s Master Dark and Master Bias subtracted from flat frames after overscan subtract. Background extractor plugin used to normalize frames by mean and combine by median into Master Flat frame
  4. Object images imported as a stack and converted from 16-bit to 32-bit
  5. Overscan Subtract performed on each slice
  6. 60s Master Dark and Master Bias frames subtracted from each
  7. Each image divided by Master Flat
  8. Images registered to # 1 using StackReg
  9. Removed images 2, 3, 7, 8, 11, 13, 14, 16, 17, 18, 22, 23, 24, 25, 26, 27, 28, and 31 from the stack due to very fuzzy stars and random extremely hot pixels in the region of the galaxy.
  10. Remaining images averaged

Comments: Unfortunately, we failed to notice the presence of fringes (ripple caused by light reflecting off the back of the CCD chip and interfering with incoming light) in the images, and so were unable to correct for them.

March 15/16 Blue Data

  1. Used same Bias and Dark frames as previous
  2. Blue filter flat frames converted from 16-bit to 32-bit. 60s Master Dark and Master Bias subtracted from flat frames after overscan subtract. Blue Master Flat created by previous method
  3. Blue images imported as a stack and converted from 16- to 32-bit
  4. Overscan Subtract performed on each slice
  5. 60s Master Dark and Master Bias subtracted from each
  6. Each image divided by Blue Master Flat
  7. Stack of images duplicated
  8. Difficulty in StackReg aligning the frames, applied a median filter to duplicate stack and then used StackReg. Resulting frames Averaged to create first image.
  9. Second stack split into individual images and then manually aligned to the first image using TurboReg.
  10. Averaging the image produced cold pixels (black areas) all over the image, frames combined by median instead.
First Image (Median filter) Second Image (Med Combined )

 

 

RawBasic ImageManipulation
BackIndex

_________________________________________

Banner image is a slice of HST image by NASA/ESA