CCD Photometry with Mathematica

Table of Contents

CCDPhotometry`Photometry`

This package contains functions useful in photometry calculations.

IsolateRegion[imobj, {i, j}, r]	Returns a submatrix of the supplied FITS image object that is centered at {i, j} and has width and height, r.
Photometry[imobj, {x, y}, rad]	Sums up all of the pixel values in the specified FITS data object that are centered at {x, y} and that lie within a circular aperture with radius, rad.
Phase[imobj, knownmin, orbper]	Returns the orbital phase of an object based on information in the header of the FITS image object. The HJD of a known minimum, knownmin, and the orbital period in days, orbper, are also required.

Photometry calculations

This loads the package

In[1]:=

This loads the packages for image loading and display.

In[2]:=

In[3]:=

This loads the data and displays it.

In[4]:=

data = ImageLoad[ToFileName["CCDPhotometry", "11080154.fit"]]

Out[4]=

In[5]:=

[Graphics:HTMLFiles/photometry_7.gif]

Here is an isolated subset of the above data.

In[6]:=

Out[6]=

{{6.48855, 6.78287, 4.7379, 4.05342, 3.35701, 6.04707, 5.72238, 3.82413, 9.20134, 5.11404, 6 ... 7587, 6.63357, 7.24095, 4.98828, 3.75972, 5.77005, 2.71831, 0.337188, 7.0739, 9.17154, 10.8554}}

Pixel sums can be performed as follows

In[7]:=

Out[7]=

For sources with known orbital periods, the orbital phase of the data can be calculated.

In[8]:=

Out[8]=

HOME

Created by Mathematica (October 24, 2004)