IV. 2MASS Data Processing


4. Point Source Detection and Photometry

a. Frame Source Detection and Aperture Photometry

Point sources were initially detected and photometered in each instrumentally-corrected Read_1 and Read_2-Read_1 frame. These frame positions and magnitudes were used to reconstruct the relative and absolute positions of the frames on the sky. The positions and magnitudes from the Read_2-Read_1 frames were not carried further in the 2MASS pipeline, being re-estimated subsequent to frame position reconstruction using the combined frames. For those sources saturated in the Read_2-Read_1 frames, the positions and aperture magnitudes estimated from the individual Read_1 frames were combined and carried forward in the pipeline and used in the catalog.

Point sources were detected in the frames by identifying local maxima. Positions were measured for these detections using a maximum-likelihood estimator, and magnitudes were estimated using aperture photometry within a 4´´ (two camera pixel) radius aperture.

The values for all the pixels were histogrammed, and the median sky value and the noise were estimated. The noise was estimated as one-half the distance between the 15.87% and the 84.13% quantiles of the histogram.

All the pixels in the frame were examined to identify local maxima. For each local maximum exceeding the local background by more than Tpk (where Tpk=4), the following computations were performed:

The centroid of the 3×3 pixel patch centered on the local maximum was first evaluated:
 
 

(Eq. IV.4.a.1)
(Eq. IV.4.a.2)

where:

(Eq. IV.4.a.3)

and (xc,yc) are relative to the center of the local maximum pixel, and dij are the data values less the local background. If Fs<Tsnr (where Tsnr=10), or any of the pixels in the 3×3 patch were masked off, the local maximum was not processed further. If no saturated pixels were encountered, the position was estimated as described in the following.

The local centroid was used as the initial source position, which was then refined using the following maximum likelihood position estimation technique, with the cost function on each axis computed by:
 

(Eq. IV.4.a.4)
(Eq. IV.4.a.5)

where:

P = the normalized point spread function (PSF) and Ft is the template amplitude,
 

(Eq. IV.4.a.6)

Cx and Cy were followed to their zero crossings, and then interpolated to solve for (xc,yc).

The fraction of flux in the peak pixel (fpeak) was computed by:
 

(Eq. IV.4.a.7)

If the following criteria were satisfied, the fpeak value is entered into a histogram used for estimating the seeing and interpolating an appropriate PSF:

The aperture photometry module was called to compute the aperture magnitude at the estimated position.  The aperture photometry on a frame was performed by summing pixels entirely within the aperture and interpolating pixels partially within the aperture. The sky background for each object was computed in an annulus with an inner radius of 24.0´´ and an outer radius of 30.0´´. Pixels in the sky annulus were entirely included or excluded based on the distance of their centers from the source. The sky value was estimated by first excluding saturated, masked, or unreasonably low pixels. A -trimmed median of the surviving sky pixels is then used as the sky estimate. Aperture photometry and positions from the maximum six possible overlapping frames were later combined using an unweighted average for the Read_1 data.

The 3×3 pixel patch centered on the peak pixel along with the coordinates of the peak pixel and the aperture magnitude and associated error were written to a temporary file so that the detection parameters could be re-estimated in a second pass using a PSF matched to the seeing.

Then:

The aperture photometry measured in this step for non-saturated Read_2-Read_1 sources from individual frames was not carried forward as part of the final source characteristics.

i. Bright Source (READ_1) Aperture Photometry

The photometric dynamic range of 2MASS was extended by the use of the 51 ms Read_1 exposures. Sources start to saturate on the 1.3 s Read_2-Read_1 exposures at magnitude levels of approximately 9.5, 9.0 and 8.5 at J, H, and Ks, respectively. For objects that were found to have one or more saturated pixels within the measurement aperture on the Read_2-Read_1 frames, the "default magnitude" quoted in the Point Source Catalog records was taken from the aperture photometry from the Read_1 frames. This is indicated by a value of "1" in the "rd_flg" parameter in the point source records, for the appropriate bands. Positions measured from the Read_1 frames were used in the final source position estimation only if Read_2-Read_1 profile fit results were not available.
 

ii. Faint Source Detection

The fainter, and thus majority of sources found by 2MASS are detected from the Atlas Images. Each Atlas Image is convolved with a zero-sum 4´´ FWHM Gaussian over a 13-pixel sub-array. The resulting zero-sum filtered image is thresholded, and for each maximum over threshold, a detection is identified and a rough position estimate is computed from the corrected centroid. This detection process is in all essential aspects identical to that used in the FIND subprocess of the DAOPHOT II program developed by Peter Stetson (P.B. Stetson 1991). The detections list is sent to the software module that computes the running estimate of the seeing during a scan, and to the photometry routines that compute the refined estimates of flux and position. The detection threshold used is 3.0 times the estimated point source noise level for the Atlas Image. The noise level is estimated as (1/0.7) times the difference between the 32.22% and the 50% quantiles of the zero-sum filtered image histogram. This noise estimator is sensitive to confusion noise and the detection threshold (in units of flux) thus increases in areas of high source density such as the galactic plane.

The list of detections from each image is passed to a subroutine which estimates the seeing shape parameter which is then passed along with the list of detections to the profile fit photometry subprocess described in the next section of this documentation.
 

iii.  Very Bright Source (Saturated READ_1) Detection and Photometry

Stars brighter than approximately magnitude=4 to 5 saturated even the 51 ms Read_1 exposures. Special techniques have been utilized for photometry and astrometry on these objects, extending the dynamic range another ~8 magnitudes to the brightest objects in the 2MASS sky. The corresponding "rd_flg" value for these measurements is "3". The photometry and astrometry for these objects is of reduced accuracy compared to the unsaturated objects. 2MASS photometric and astrometric accuracy for these objects were required to be sufficient to allow reliable identification of associated bright star artifacts, but accuracies achieved may be sufficient for additional purposes. It should be noted that many of the brightest near infrared stars are variables and may also have high proper motions.

Photometry on these stars has been done using a 1d radial profile fit to the unsaturated wings of the point spread function, out to a snr limited radius (or 60 arcseconds for the brightest objects). Positions in the individual data frames were estimated by computing the flux weighted centroid of the pixels over a snr threshold  The single frame positions for each band were combined during position reconstruction in a manner analogous to that used for unsaturated Read_1 measurements.

The brightest star in the 2MASS Catalog is alpha Orionis (M2 Ia), with J= -2.989, H=-4.007, Ks= -4.378 mag.

Saturated Read_1 Frame Position Estimation

The positions for stars saturated in the Read_1 frames were estimated using a flux weighted centroid. When a saturated pixel was encountered in an Read_1 frame during Read_1 source extraction, connected pixels over a threshold (500dn) were examined and the flux weighted centroid over those pixels was computed. A second centroid was also computed for that subset of pixels also over a second higher threshold (20000dn). If more than 25 pixels were counted over the second threshold, the second centroid was used for the source position, if less, the first centroid was used. The second centroid avoided  problems due to asymmetries in the outer halo of very bright stars, and allowed measurements to be made close to frame edges for those stars. These 2 centroids were found to be essentially identical for sources with 15 to 50 pixels over the second threshold. This technique gives good astrometry for sources from those just over the Read_1 saturation threshold to the brightest sources in the 2MASS Catalog.

Combination of the frame positions and subsequent astrometric processing is discussed in the position reconstruction section of this documentation.

Saturated Read_1 Photometry

Photometry for a star saturated in Read_1 was done using a 1d radial profile fit to the combined data from all 6 frames covering that star in a scan. When saturated pixels were encountered in an Read_1 frame during Read_1 source extraction, the frame data was written to a temporary file for later examination. After all the Read_1 sources were identified in the single data frames, all the data for each saturated star in the scan was registered and converted to radial profiles. The radial profile data was fit to analytic models for the camera and band being processed. The radial profile models are in the following form:
 

(Eq. IV.4.a.8)

The logarithm of that portion of the radial profile data that was unsaturated and of sufficient signal to noise ratio was fit to the logarithm of the model profile using least-square minimization of the errors. The square root of the variance of the residuals from this fit was converted to a magnitude and is reported as the magnitude sigma. It should be noted that this is derived from the fit residuals and is not a formal error.

The bright source model profiles do not account for profile variations due to seeing. A corrective model for estimated magnitude dependence on seeing has been developed but was not applied during 2MASS processing. Later analysis found some deficiencies in the corrective model for seeing. Seeing corrections can be as large as several tenths of a magnitude. Discussion and analysis of these issues may be found here.

A detailed description of these algorithms and the model parameters used may be found here.

Although there are issues with the seeing corrections for the saturated Read_1 magnitudes, accuracies of 0.1-0.2 magnitudes have been demonstrated throughout most of the range of 8 magnitudes, as well as astrometry on the order of 170 mas.

The confusion radius is larger for bright stars than for the rest of the catalog. Although many saturated doubles are separated down to 7 arcsec, a few very bright saturated doubles with larger separations were not resolved. The worst example known is alpha Centauri AB with components 17.7 arcsec apart which were not separated in the 2MASS catalog. A small number of other bright stars are known to be missing measurements in one or two bands or missing entirely from the 2MASS catalog due to processing problems.

[Last Updated: 2003 Feb 27; by R. Cutri & E. Kopan]


Previous page.Next page.
Return to Explanatory Supplement TOC Page.