Noise Variance Estimation In Signal Processing

We present a new method of estimating noise variance. The method is applicable for 1D and 2D signal processing. The essence of this method is estimation of the scatter of normally distributed data with high level of outliers. The method is applicable to data with the majority of the data points having no signal present. The method is based on the shortest half sample method. The mean of the shortest half sample (shorth) and the location of the least median of squares are among the most robust measures of the location of the mode. The length of the shortest half sample has been used as the measurement of the data scatter of uncontaminated data. We show that computing the length of several sub samples of varying sizes provides the necessary information to estimate both the scatter and the number of uncontaminated data points in a sample. We derive the system of equations to solve for the data scatter and the number of uncontaminated data points for the Gaussian distribution. The data scatter is the measure of the noise variance. The method can be extended to other distributions

Fast Direct Plane-to-Plane Coordinate Transformations

A set of highly efficient coordinate transformations has been derived to convert between Cartesian coordinates in two different projection planes. The gain in execution speed is as great as a factor of 20 compared to the two‐step process of going from Cartesian coordinates in one plane to the celestial coordinates, and then to Cartesian coordinates in the second plane. The new coordinate transformations are successfully used by the Spitzer Space Telescope mosaicker MOPEX

Mosaicking with MOPEX

We present MOPEX - a software package for image mosaicking and point source extraction. MOPEX has been developed for the pitzer Space Telescope. This paper concentrates on the mosaicking aspects of the package. MOPEX features the use of several interpolation techniques, coaddition schemes, and robust and flexible outlier detection based on spatial and temporal filtering. A number of original algorithms have been designed and implemented in MOPEX. Among them is direct plane-to-plane coordinate transformation, which allows at least an order of magnitude speed up in performing coordinate transformation by bypassing the sky coordinates. The dual outlier detection makes possible outlier detection in the areas of even minimal redundancy. Image segmentation based on adaptive thresholding is used for object detection, which is part of outlier detection. Efficient use of computer memory allows mosaicking of data sets of very deep coverage of thousands of images per pointing, as well as areas of sky covering many square degrees. Although designed for Spitzer data, MOPEX does not require any Spitzer-specific fits header keywords to run, and can be applied to other data, that have standard header information on the image geometry and pointing. The package is available for distribution at http://ssc.spitzer.caltech.edu/postbcd/

Efficient Mosaicking of Spitzer Space Telescope Images

A parallel version of the MOPEX software, which generates mosaics of infrared astronomical images acquired by the Spitzer Space Telescope, extends the capabilities of the prior serial version. In the parallel version, both the input image space and the output mosaic space are divided among the available parallel processors. This is the only software that performs the point-source detection and the rejection of spurious imaging effects of cosmic rays required by Spitzer scientists. This software includes components that implement outlier-detection algorithms that can be fine-tuned for a particular set of image data by use of a number of adjustable parameters. This software has been used to construct a mosaic of the Spitzer Infrared Array Camera Shallow Survey, which comprises more than 17,000 exposures in four wavelength bands from 3.6 to 8 m and spans a solid angle of about 9 square degrees. When this software was executed on 32 nodes of the 1,024-processor Cosmos cluster computer at NASA s Jet Propulsion Laboratory, a speedup of 8.3 was achieved over the serial version of MOPEX. The performance is expected to improve dramatically once a true parallel file system is installed on Cosmos

Point Source Extraction with MOPEX

MOPEX (MOsaicking and Point source EXtraction) is a package developed at the Spitzer Science Center for astronomical image processing. We report on the point source extraction capabilities of MOPEX. Point source extraction is implemented as a two step process: point source detection and profile fitting. Non-linear matched filtering of input images can be performed optionally to increase the signal-to-noise ratio and improve detection of faint point sources. Point Response Function (PRF) fitting of point sources produces the final point source list which includes the fluxes and improved positions of the point sources, along with other parameters characterizing the fit. Passive and active deblending allows for successful fitting of confused point sources. Aperture photometry can also be computed for every extracted point source for an unlimited number of aperture sizes. PRF is estimated directly from the input images. Implementation of efficient methods of background and noise estimation, and modified Simplex algorithm contribute to the computational efficiency of MOPEX. The package is implemented as a loosely connected set of perl scripts, where each script runs a number of modules written in C/C++. Input parameter setting is done through namelists, ASCII configuration files. We present applications of point source extraction to the mosaic images taken at 24 and 70 micron with the Multiband Imaging Photometer (MIPS) as part of the Spitzer extragalactic First Look Survey and to a Digital Sky Survey image. Completeness and reliability of point source extraction is computed using simulated data.Comment: 20 pages, 13 Postscript figures, accepted for publication in PAS

The Spitzer Survey of the Small Magellanic Cloud: Discovery of Embedded Protostars in the HII Region NGC 346

We use Spitzer Space Telescope observations from the Spitzer Survey of the Small Magellanic Cloud (S3MC) to study the young stellar content of N66, the largest and brightest HII region in the SMC. In addition to large numbers of normal stars, we detect a significant population of bright, red infrared sources that we identify as likely to be young stellar objects (YSOs). We use spectral energy distribution (SED) fits to classify objects as ordinary (main sequence or red giant) stars, asymptotic giant branch stars, background galaxies, and YSOs. This represents the first large-scale attempt at blind source classification based on Spitzer SEDs in another galaxy. We firmly identify at least 61 YSOs, with another 50 probable YSOs; only one embedded protostar in the SMC was reported in the literature prior to the S3MC. We present color selection criteria that can be used to identify a relatively clean sample of YSOs with IRAC photometry. Our fitted SEDs indicate that the infrared-bright YSOs in N66 have stellar masses ranging from 2 Msun to 17 Msun, and that approximately half of the objects are Stage II protostars, with the remaining YSOs roughly evenly divided between Stage I and Stage III sources. We find evidence for primordial mass segregation in the HII region, with the most massive YSOs being preferentially closer to the center than lower-mass objects. Despite the low metallicity and dust content of the SMC, the observable properties of the YSOs appear consistent with those in the Milky Way. Although the YSOs are heavily concentrated within the optically bright central region of N66, there is ongoing star formation throughout the complex and we place a lower limit on the star formation rate of 3.2 x 10^-3 Msun/yr over the last ~1 Myr.Comment: 13 pages, 5 figures (3 in color), 2 tables. Accepted for publication in Ap

Spitzer 24-micron Time-Series Observations of the Eclipsing M-dwarf Binary GU Bootis

We present a set of {\it Spitzer} 24$\mu$m MIPS time series observations of the M-dwarf eclipsing binary star GU Bo\"otis. Our data cover three secondary eclipses of the system: two consecutive events and an additional eclipse six weeks later. The study's main purpose is the long wavelength (and thus limb darkening-independent) characterization of GU Boo's light curve, allowing for independent verification of the results of previous optical studies. Our results confirm previously obtained system parameters. We further compare GU Boo's measured 24$\mu$m flux density to the value predicted by spectral fitting and find no evidence for circumstellar dust. In addition to GU Boo, we characterize (and show examples of) light curves of other objects in the field of view. Analysis of these light curves serves to characterize the photometric stability and repeatability of {\it Spitzer's} MIPS 24\micron array over short (days) and long (weeks) timescales at flux densities between approximately 300--2,000$\mu$Jy. We find that the light curve root mean square about the median level falls into the 1--4% range for flux densities higher than 1mJy. Finally, we comment on the fluctuations of the 24\micron background on short and long timescales.Comment: ApJ accepted. 10 pages, 12 figure

Normal and unusual transient events in IRAC images

The Spitzer Space Telescope Infrared Array Camera (IRAC) is a four-channel camera that uses two pairs of 256 x 256 pixel InSb and Si:As IBC detectors to provide simultaneous images at 3.6, 4.5, 5.8, and 8 microns. IRAC experiences a flux of cosmic rays that produce transient events in images from each of the arrays, with 5-7 pixels per second being affected in an IRAC integration. The vast majority of these transient events can be adequately characterized so they can be effectively detected and flagged by a pipeline software module. However, because of the nature of the arrays and their arrangement in the camera structure, a small fraction of the cosmic ray hits on IRAC produce transients with unusual morphologies which cannot be characterized in a general way. We present nominal cosmic ray rates observed for IRAC on-orbit and rates observed during a period of elevated solar proton flux following a series of X-class solar flares in late 2003. We also present a guide for observers to help identify unusual transient events in their data. We comment on the physical nature of the production of many o9f these unusual transients and how this mechanism differs from the production of "normal" transient events

Normal and unusual transient events in IRAC images

The Spitzer Space Telescope Infrared Array Camera (IRAC) is a four-channel camera that uses two pairs of 256 x 256 pixel InSb and Si:As IBC detectors to provide simultaneous images at 3.6, 4.5, 5.8, and 8 microns. IRAC experiences a flux of cosmic rays that produce transient events in images from each of the arrays, with 5-7 pixels per second being affected in an IRAC integration. The vast majority of these transient events can be adequately characterized so they can be effectively detected and flagged by a pipeline software module. However, because of the nature of the arrays and their arrangement in the camera structure, a small fraction of the cosmic ray hits on IRAC produce transients with unusual morphologies which cannot be characterized in a general way. We present nominal cosmic ray rates observed for IRAC on-orbit and rates observed during a period of elevated solar proton flux following a series of X-class solar flares in late 2003. We also present a guide for observers to help identify unusual transient events in their data. We comment on the physical nature of the production of many o9f these unusual transients and how this mechanism differs from the production of "normal" transient events

Spitzer IRS 16 micron Observations of the GOODS Fields

We present Spitzer 16 micron imaging of the Great Observatories Origins Deep Survey (GOODS) fields. We survey 150 square arcminutes in each of the two GOODS fields (North and South), to an average 3 sigma depth of 40 and 65 micro-Jy respectively. We detect about 1300 sources in both fields combined. We validate the photometry using the 3-24 micron spectral energy distribution of stars in the fields compared to Spitzer spectroscopic templates. Comparison with ISOCAM and AKARI observations in the same fields show reasonable agreement, though the uncertainties are large. We provide a catalog of photometry, with sources cross correlated with available Spitzer, Chandra, and HST data. Galaxy number counts show good agreement with previous results from ISOCAM and AKARI, with improved uncertainties. We examine the 16 to 24 micron flux ratio and find that for most sources it lies within the expected locus for starbursts and infrared luminous galaxies. A color cut of S_{16}/S_{24}>1.4 selects mostly sources which lie at 1.1<z<1.6, where the 24 micron passband contains both the redshifted 9.7 micron silicate absorption and the minimum between PAH emission peaks. We measure the integrated galaxy light of 16 micron sources, and find a lower limit on the galaxy contribution to the extragalactic background light at this wavelength to be 2.2\pm 0.2$ nW m^{-2} sr^{-1}.Comment: Accepted for Publication in the AJ. 53 preprint pages, including 15 figures and 8 tables. Table 1-4 are truncated in the ms.tex but are included in full in the tar file (and will be available in the online version of the AJ