Rational-operator-based depth-from-defocus approach to scene reconstruction

This paper presents a rational-operator-based approach to depth from defocus (DfD) for the reconstruction of three-dimensional scenes from two-dimensional images, which enables fast DfD computation that is independent of scene textures. Two variants of the approach, one using the Gaussian rational operators (ROs) that are based on the Gaussian point spread function (PSF) and the second based on the generalized Gaussian PSF, are considered. A novel DfD correction method is also presented to further improve the performance of the approach. Experimental results are considered for real scenes and show that both approaches outperform existing RO-based methods

GalPak3D: A Bayesian parametric tool for extracting morpho-kinematics of galaxies from 3D data

We present a method to constrain galaxy parameters directly from three-dimensional data cubes. The algorithm compares directly the data with a parametric model mapped in $x,y,\lambda$ coordinates. It uses the spectral lines-spread function (LSF) and the spatial point-spread function (PSF) to generate a three-dimensional kernel whose characteristics are instrument specific or user generated. The algorithm returns the intrinsic modeled properties along with both an `intrinsic' model data cube and the modeled galaxy convolved with the 3D-kernel. The algorithm uses a Markov Chain Monte Carlo (MCMC) approach with a nontraditional proposal distribution in order to efficiently probe the parameter space. We demonstrate the robustness of the algorithm using 1728 mock galaxies and galaxies generated from hydrodynamical simulations in various seeing conditions from 0.6" to 1.2". We find that the algorithm can recover the morphological parameters (inclination, position angle) to within 10% and the kinematic parameters (maximum rotation velocity) to within 20%, irrespectively of the PSF in seeing (up to 1.2") provided that the maximum signal-to-noise ratio (SNR) is greater than $\sim3$ pixel$^{-1}$ and that the ratio of the galaxy half-light radius to seeing radius is greater than about 1.5. One can use such an algorithm to constrain simultaneously the kinematics and morphological parameters of (nonmerging) galaxies observed in nonoptimal seeing conditions. The algorithm can also be used on adaptive-optics (AO) data or on high-quality, high-SNR data to look for nonaxisymmetric structures in the residuals.Comment: 16 pages, 10 figures, accepted to publication in AJ, revised version after proofs corrections. Algorithm available at http://galpak.irap.omp.e

Data Reduction Pipeline for the CHARIS Integral-Field Spectrograph I: Detector Readout Calibration and Data Cube Extraction

We present the data reduction pipeline for CHARIS, a high-contrast integral-field spectrograph for the Subaru Telescope. The pipeline constructs a ramp from the raw reads using the measured nonlinear pixel response, and reconstructs the data cube using one of three extraction algorithms: aperture photometry, optimal extraction, or $\chi^2$ fitting. We measure and apply both a detector flatfield and a lenslet flatfield and reconstruct the wavelength- and position-dependent lenslet point-spread function (PSF) from images taken with a tunable laser. We use these measured PSFs to implement a $\chi^2$-based extraction of the data cube, with typical residuals of ~5% due to imperfect models of the undersampled lenslet PSFs. The full two-dimensional residual of the $\chi^2$ extraction allows us to model and remove correlated read noise, dramatically improving CHARIS' performance. The $\chi^2$ extraction produces a data cube that has been deconvolved with the line-spread function, and never performs any interpolations of either the data or the individual lenslet spectra. The extracted data cube also includes uncertainties for each spatial and spectral measurement. CHARIS' software is parallelized, written in Python and Cython, and freely available on github with a separate documentation page. Astrometric and spectrophotometric calibrations of the data cubes and PSF subtraction will be treated in a forthcoming paper.Comment: 18 pages, 15 figures, 3 tables, replaced with JATIS accepted version (emulateapj formatted here). Software at https://github.com/PrincetonUniversity/charis-dep and documentation at http://princetonuniversity.github.io/charis-de

Decomposition of AGN host galaxy images

We describe an algorithm to decompose deep images of Active Galactic Nuclei into host galaxy and nuclear components. Currently supported are three galaxy models: A de-Vaucouleurs spheroidal, an exponential disc, and a two-component disc+bulge model. Key features of the method are: (semi-)analytic representation of a possibly spatially variable point-spread function; full two-dimensional convolution of the model galaxy using gradient-controlled adaptive subpixelling; multiple iteration scheme. The code is computationally efficient and versatile for a wide range of applications. The quantitative performance is measured by analysing simulated imaging data. We also present examples of the application of the method to small test samples of nearby Seyfert 1 galaxies and quasars at redshifts z < 0.35.Comment: 12 pages, 15 figures, accepted for publication in MNRA

Probing the LMC age gap at intermediate cluster masses

The LMC has a rich star cluster system spanning a wide range of ages and masses. One striking feature of the LMC cluster system is the existence of an age gap between 3-10 Gyrs. But this feature is not as clearly seen among field stars. Three LMC fields containing relatively poor and sparse clusters whose integrated colours are consistent with those of intermediate age simple stellar populations have been imaged in BVI with the Optical Imager (SOI) at the Southern Telescope for Astrophysical Research (SOAR). A total of 6 clusters, 5 of them with estimated initial masses M < 10^4M_sun, were studied in these fields. Photometry was performed and Colour-Magnitude Diagrams (CMD) were built using standard point spread function fitting methods. The faintest stars measured reach V ~ 23. The CMD was cleaned from field contamination by making use of the three-dimensional colour and magnitude space available in order to select stars in excess relative to the field. A statistical CMD comparison method was developed for this purpose. The subtraction method has proven to be successful, yielding cleaned CMDs consistent with a simple stellar population. The intermediate age candidates were found to be the oldest in our sample, with ages between 1-2 Gyrs. The remaining clusters found in the SOAR/SOI have ages ranging from 100 to 200 Myrs. Our analysis has conclusively shown that none of the relatively low-mass clusters studied by us belongs to the LMC age-gap.Comment: 9 pages, 8 figures. Accepted to MNRA

The Murchison Widefield Array: Design Overview

The Murchison Widefield Array (MWA) is a dipole-based aperture array synthesis telescope designed to operate in the 80-300 MHz frequency range. It is capable of a wide range of science investigations, but is initially focused on three key science projects. These are detection and characterization of 3-dimensional brightness temperature fluctuations in the 21cm line of neutral hydrogen during the Epoch of Reionization (EoR) at redshifts from 6 to 10, solar imaging and remote sensing of the inner heliosphere via propagation effects on signals from distant background sources,and high-sensitivity exploration of the variable radio sky. The array design features 8192 dual-polarization broad-band active dipoles, arranged into 512 tiles comprising 16 dipoles each. The tiles are quasi-randomly distributed over an aperture 1.5km in diameter, with a small number of outliers extending to 3km. All tile-tile baselines are correlated in custom FPGA-based hardware, yielding a Nyquist-sampled instantaneous monochromatic uv coverage and unprecedented point spread function (PSF) quality. The correlated data are calibrated in real time using novel position-dependent self-calibration algorithms. The array is located in the Murchison region of outback Western Australia. This region is characterized by extremely low population density and a superbly radio-quiet environment,allowing full exploitation of the instrumental capabilities.Comment: 9 pages, 5 figures, 1 table. Accepted for publication in Proceedings of the IEE