Restoration of Hyperspectral Astronomical Data with Spectrally Varying Blur

This book is a collection of 19 articles which reflect the courses given at the Collège de France/Summer school “Reconstruction d'images − Applications astrophysiques“ held in Nice and Fréjus, France, from June 18 to 22, 2012. The articles presented in this volume address emerging concepts and methods that are useful in the complex process of improving our knowledge of the celestial objects, including Earth

Image reconstruction in optical interferometry: Benchmarking the regularization

With the advent of infrared long-baseline interferometers with more than two telescopes, both the size and the completeness of interferometric data sets have significantly increased, allowing images based on models with no a priori assumptions to be reconstructed. Our main objective is to analyze the multiple parameters of the image reconstruction process with particular attention to the regularization term and the study of their behavior in different situations. The secondary goal is to derive practical rules for the users. Using the Multi-aperture image Reconstruction Algorithm (MiRA), we performed multiple systematic tests, analyzing 11 regularization terms commonly used. The tests are made on different astrophysical objects, different (u,v) plane coverages and several signal-to-noise ratios to determine the minimal configuration needed to reconstruct an image. We establish a methodology and we introduce the mean-square errors (MSE) to discuss the results. From the ~24000 simulations performed for the benchmarking of image reconstruction with MiRA, we are able to classify the different regularizations in the context of the observations. We find typical values of the regularization weight. A minimal (u,v) coverage is required to reconstruct an acceptable image, whereas no limits are found for the studied values of the signal-to-noise ratio. We also show that super-resolution can be achieved with increasing performance with the (u,v) coverage filling. Using image reconstruction with a sufficient (u,v) coverage is shown to be reliable. The choice of the main parameters of the reconstruction is tightly constrained. We recommend that efforts to develop interferometric infrastructures should first concentrate on the number of telescopes to combine, and secondly on improving the accuracy and sensitivity of the arrays.Comment: 15 pages, 16 figures; accepted in A&

Milli-arcsecond images of the Herbig Ae star HD 163296

The very close environments of young stars are the hosts of fundamental physical processes, such as planet formation, star-disk interactions, mass accretion, and ejection. The complex morphological structure of these environments has been confirmed by the now quite rich data sets obtained for a few objects by near-infrared long-baseline interferometry. We gathered numerous interferometric measurements for the young star HD163296 with various interferometers (VLTI, IOTA, KeckI and CHARA), allowing for the first time an image independent of any a priori model to be reconstructed. Using the Multi-aperture image Reconstruction Algorithm (MiRA), we reconstruct images of HD 163296 in the H and K bands. We compare these images with reconstructed images obtained from simulated data using a physical model of the environment of HD 163296. We obtain model-independent $H$ and $K$-band images of the surroundings of HD 163296. The images present several significant features that we can relate to an inclined asymmetric flared disk around HD 163296 with the strongest intensity at about 4-5 mas. Because of the incomplete spatial frequency coverage, we cannot state whether each of them individually is peculiar in any way. For the first time, milli-arcsecond images of the environment of a young star are produced. These images confirm that the morphology of the close environment of young stars is more complex than the simple models used in the literature so far.Comment: 11 pages, 10 figures, accepted A&A pape

A scene model of exosolar systems for use in planetary detection and characterisation simulations

Instrumental projects that will improve the direct optical finding and characterisation of exoplanets have advanced sufficiently to trigger organized investigation and development of corresponding signal processing algorithms. The first step is the availability of field-of-view (FOV) models. These can then be submitted to various instrumental models, which in turn produce simulated data, enabling the testing of processing algorithms. We aim to set the specifications of a physical model for typical FOVs of these instruments. The dynamic in resolution and flux between the various sources present in such a FOV imposes a multiscale, independent layer approach. From review of current literature and through extrapolations from currently available data and models, we derive the features of each source-type in the field of view likely to pass the instrumental filter at exo-Earth level. Stellar limb darkening is shown to cause bias in leakage calibration if unaccounted for. Occurrence of perturbing background stars or galaxies in the typical FOV is unlikely. We extract galactic interstellar medium background emissions for current target lists. Galactic background can be considered uniform over the FOV, and it should show no significant drift with parallax. Our model specifications have been embedded into a Java simulator, soon to be made open-source. We have also designed an associated FITS input/output format standard that we present here.Comment: 9 pages (+5 of appendices), 7 figures, accepted for publication in Astronomy & Astrophysic

Super-resolution in map-making based on a physical instrument model and regularized inversion. Application to SPIRE/Herschel

We investigate super-resolution methods for image reconstruction from data provided by a family of scanning instruments like the Herschel observatory. To do this, we constructed a model of the instrument that faithfully reflects the physical reality, accurately taking the acquisition process into account to explain the data in a reliable manner. The inversion, ie the image reconstruction process, is based on a linear approach resulting from a quadratic regularized criterion and numerical optimization tools. The application concerns the reconstruction of maps for the SPIRE instrument of the Herschel observatory. The numerical evaluation uses simulated and real data to compare the standard tool (coaddition) and the proposed method. The inversion approach is capable to restore spatial frequencies over a bandwidth four times that possible with coaddition and thus to correctly show details invisible on standard maps. The approach is also applied to real data with significant improvement in spatial resolution.Comment: Astronomy & Astrophysic

A low optical depth region in the inner disk of the HerbigAe star HR5999

Circumstellar disks surrounding young stars are known to be the birthplaces of planets, and the innermost astronomical unit is of particular interest. We present new long-baseline spectro-interferometric observations of the HerbigAe star, HR5999, obtained in the H and K bands with the AMBER instrument at the VLTI, and aim to produce near-infrared images at the sub-AU spatial scale. We spatially resolve the circumstellar material and reconstruct images using the MiRA algorithm. In addition, we interpret the interferometric observations using models that assume that the near-infrared excess is dominated by the emission of a circumstellar disk. We compare the images reconstructed from the VLTI measurements to images obtained using simulated model data. The K-band image reveals three main elements: a ring-like feature located at ~0.65 AU, a low surface brightness region inside, and a central spot. At the maximum angular resolution of our observations (1.3 mas), the ring is resolved while the central spot is only marginally resolved, preventing us from revealing the exact morphology of the circumstellar environment. We suggest that the ring traces silicate condensation, i.e., an opacity change, in a circumstellar disk around HR 5999. We build a model that includes a ring at the silicate sublimation radius and an inner disk of low surface brightness responsible for a large amount of the near-infrared continuum emission. The model successfully fits the SED, visibilities, and closure phases, and provides evidence of a low surface brightness region inside the silicate sublimation radius. This study provides additional evidence that in HerbigAe stars, there is material in a low surface brightness region, probably a low optical depth region, located inside the silicate sublimation radius and of unknown nature.Comment: 11 pages, 10 figure

First astronomical unit scale image of the GW Ori triple. Direct detection of a new stellar companion

Young and close multiple systems are unique laboratories to probe the initial dynamical interactions between forming stellar systems and their dust and gas environment. Their study is a key building block to understanding the high frequency of main-sequence multiple systems. However, the number of detected spectroscopic young multiple systems that allow dynamical studies is limited. GW Orionis is one such system. It is one of the brightest young T Tauri stars and is surrounded by a massive disk. Our goal is to probe the GW Orionis multiplicity at angular scales at which we can spatially resolve the orbit. We used the IOTA/IONIC3 interferometer to probe the environment of GW Orionis with an astronomical unit resolution in 2003, 2004, and 2005. By measuring squared visibilities and closure phases with a good UV coverage we carry out the first image reconstruction of GW Ori from infrared long-baseline interferometry. We obtain the first infrared image of a T Tauri multiple system with astronomical unit resolution. We show that GW Orionis is a triple system, resolve for the first time the previously known inner pair (separation $\rho\sim$1.4 AU) and reveal a new more distant component (GW Ori C) with a projected separation of $\sim$8 AU with direct evidence of motion. Furthermore, the nearly equal (2:1) H-band flux ratio of the inner components suggests that either GW Ori B is undergoing a preferential accretion event that increases its disk luminosity or that the estimate of the masses has to be revisited in favour of a more equal mass-ratio system that is seen at lower inclination. Accretion disk models of GW Ori will need to be completely reconsidered because of this outer companion C and the unexpected brightness of companion B.Comment: 5 pages, 9 figures, accepted Astronomy and Astrophysics Letters. 201

3-D deconvolution of hyper-spectral astronomical data

In this paper we present a general forward fitting method for multichannel image restoration based on regularized chi2. We introduce separable regularizations that account for the dynamic of the model and take advantage of the continuities present in the data, leaving only two hyper-parameters to tune. We illustrate a practical implementation of this method in the context of host galaxy subtraction for the Nearby SuperNova factory. We show that the image restoration obtained fulfills the stringent requirements on bias and photometricity needed by this program. The reconstruction yields sub-percent integrated residuals in all the synthetic filters considered both on real and simulated data. Even though our implementation is tied to the SNfactory data, the method translates to any hyper-spectral data. As such, it is of direct relevance to several new generation instruments like MUSE. Also, this technique could be applied to multi-band astronomical imaging for which image reconstruction is important, for example to increase image resolution for weak lensing surveys.Comment: 14 pages, 12 figures, 3 tables. Accepted for publication in MNRA

Imaging the dynamical atmosphere of the red supergiant Betelgeuse in the CO first overtone lines with VLTI/AMBER

We present the first 1-D aperture synthesis imaging of the red supergiant Betelgeuse in the individual CO first overtone lines with VLTI/AMBER. The reconstructed 1-D projection images reveal that the star appears differently in the blue wing, line center, and red wing of the individual CO lines. The 1-D projection images in the blue wing and line center show a pronounced, asymmetrically extended component up to ~1.3 stellar radii, while those in the red wing do not show such a component. The observed 1-D projection images in the lines can be reasonably explained by a model in which the CO gas within a region more than half as large as the stellar size is moving slightly outward with 0--5 km s^-1, while the gas in the remaining region is infalling fast with 20--30 km s^-1. A comparison between the CO line AMBER data taken in 2008 and 2009 shows a significant time variation in the dynamics of the CO line-forming region in the photosphere and the outer atmosphere. In contrast to the line data, the reconstructed 1-D projection images in the continuum show only a slight deviation from a uniform disk or limb-darkened disk. We derive a uniform-disk diameter of 42.05 +/- 0.05 mas and a power-law-type limb-darkened disk diameter of 42.49 +/- 0.06 mas and a limb-darkening parameter of (9.7 +/- 0.5) x 10^{-2}. This latter angular diameter leads to an effective temperature of 3690 +/- 54 K for the continuum-forming layer. These diameters confirm that the near-IR size of Betelgeuse was nearly constant over the last 18 years, in marked contrast to the recently reported noticeable decrease in the mid-IR size. The continuum data taken in 2008 and 2009 reveal no or only marginal time variations, much smaller than the maximum variation predicted by the current 3-D convection simulations.Comment: 21 pages, 12 figures, accepted for publication in Astronomy and Astrophysic

Temporal trends of population viral suppression in the context of Universal Test and Treat: the ANRS 12249 TasP trial in rural South Africa

Introduction The universal test‐and‐treat (UTT) strategy aims to maximize population viral suppression (PVS), that is, the proportion of all people living with HIV (PLHIV) on antiretroviral treatment (ART) and virally suppressed, with the goal of reducing HIV transmission at the population level. This article explores the extent to which temporal changes in PVS explain the observed lack of association between universal treatment and cumulative HIV incidence seen in the ANRS 12249 TasP trial conducted in rural South Africa. Methods The TasP cluster‐randomized trial (2012 to 2016) implemented six‐monthly repeat home‐based HIV counselling and testing (RHBCT) and referral of PLHIV to local HIV clinics in 2 × 11 clusters opened sequentially. ART was initiated according to national guidelines in control clusters and regardless of CD4 count in intervention clusters. We measured residency status, HIV status, and HIV care status for each participant on a daily basis. PVS was computed per cluster among all resident PLHIV (≥16, including those not in care) at cluster opening and daily thereafter. We used a mixed linear model to explore time patterns in PVS, adjusting for sociodemographic changes at the cluster level. Results 8563 PLHIV were followed. During the course of the trial, PVS increased significantly in both arms (23.5% to 46.2% in intervention, +22.8, p < 0.001; 26.0% to 44.6% in control, +18.6, p < 0.001). That increase was similar in both arms (p = 0.514). In the final adjusted model, PVS increase was most associated with increased RHBCT and the implementation of local trial clinics (measured by time since cluster opening). Contextual changes (measured by calendar time) also contributed slightly. The effect of universal ART (trial arm) was positive but limited. Conclusions PVS was improved significantly but similarly in both trial arms, explaining partly the null effect observed in terms of cumulative HIV incidence between arms. The PVS gains due to changes in ART‐initiation guidelines alone are relatively small compared to gains obtained by strategies to maximize testing and linkage to care. The achievement of the 90‐90‐90 targets will not be met if the operational and implementational challenges limiting access to care and treatment, often context‐specific, are not properly addressed. Clinical trial number: NCT01509508 (clinicalTrials.gov)/DOH‐27‐0512‐3974 (South African National Clinical Trials Register)