Lossless Astronomical Image Compression and the Effects of Noise

We compare a variety of lossless image compression methods on a large sample of astronomical images and show how the compression ratios and speeds of the algorithms are affected by the amount of noise in the images. In the ideal case where the image pixel values have a random Gaussian distribution, the equivalent number of uncompressible noise bits per pixel is given by Nbits =log2(sigma * sqrt(12)) and the lossless compression ratio is given by R = BITPIX / Nbits + K where BITPIX is the bit length of the pixel values and K is a measure of the efficiency of the compression algorithm. We perform image compression tests on a large sample of integer astronomical CCD images using the GZIP compression program and using a newer FITS tiled-image compression method that currently supports 4 compression algorithms: Rice, Hcompress, PLIO, and GZIP. Overall, the Rice compression algorithm strikes the best balance of compression and computational efficiency; it is 2--3 times faster and produces about 1.4 times greater compression than GZIP. The Rice algorithm produces 75%--90% (depending on the amount of noise in the image) as much compression as an ideal algorithm with K = 0. The image compression and uncompression utility programs used in this study (called fpack and funpack) are publicly available from the HEASARC web site. A simple command-line interface may be used to compress or uncompress any FITS image file.Comment: 20 pages, 9 figures, to be published in PAS

Optimal Compression of Floating-point Astronomical Images Without Significant Loss of Information

We describe a compression method for floating-point astronomical images that gives compression ratios of 6 -- 10 while still preserving the scientifically important information in the image. The pixel values are first preprocessed by quantizing them into scaled integer intensity levels, which removes some of the uncompressible noise in the image. The integers are then losslessly compressed using the fast and efficient Rice algorithm and stored in a portable FITS format file. Quantizing an image more coarsely gives greater image compression, but it also increases the noise and degrades the precision of the photometric and astrometric measurements in the quantized image. Dithering the pixel values during the quantization process can greatly improve the precision of measurements in the images. This is especially important if the analysis algorithm relies on the mode or the median which would be similarly quantized if the pixel values are not dithered. We perform a series of experiments on both synthetic and real astronomical CCD images to quantitatively demonstrate that the magnitudes and positions of stars in the quantized images can be measured with the predicted amount of precision. In order to encourage wider use of these image compression methods, we have made available a pair of general-purpose image compression programs, called fpack and funpack, which can be used to compress any FITS format image.Comment: Accepted PAS

The Future is Hera! Analyzing Astronomical Over the Internet

Hera is the data processing facility provided by the High Energy Astrophysics Science Archive Research Center (HEASARC) at the NASA Goddard Space Flight Center for analyzing astronomical data. Hera provides all the pre-installed software packages, local disk space, and computing resources need to do general processing of FITS format data files residing on the users local computer, and to do research using the publicly available data from the High ENergy Astrophysics Division. Qualified students, educators and researchers may freely use the Hera services over the internet of research and educational purposes

DDR2 controls breast tumor stiffness and metastasis by regulating integrin mediated mechanotransduction in CAFs

Biomechanical changes in the tumor microenvironment influence tumor progression and metastases. Collagen content and fiber organization within the tumor stroma are major contributors to biomechanical changes (e., tumor stiffness) and correlated with tumor aggressiveness and outcome. What signals and in what cells control collagen organization within the tumors, and how, is not fully understood. We show in mouse breast tumors that the action of the collagen receptor DDR2 in CAFs controls tumor stiffness by reorganizing collagen fibers specifically at the tumor-stromal boundary. These changes were associated with lung metastases. The action of DDR2 in mouse and human CAFs, and tumors in vivo, was found to influence mechanotransduction by controlling full collagen-binding integrin activation via Rap1-mediated Talin1 and Kindlin2 recruitment. The action of DDR2 in tumor CAFs is thus critical for remodeling collagen fibers at the tumor-stromal boundary to generate a physically permissive tumor microenvironment for tumor cell invasion and metastases

Feasibility and performances of compressed-sensing and sparse map-making with Herschel/PACS data

The Herschel Space Observatory of ESA was launched in May 2009 and is in operation since. From its distant orbit around L2 it needs to transmit a huge quantity of information through a very limited bandwidth. This is especially true for the PACS imaging camera which needs to compress its data far more than what can be achieved with lossless compression. This is currently solved by including lossy averaging and rounding steps on board. Recently, a new theory called compressed-sensing emerged from the statistics community. This theory makes use of the sparsity of natural (or astrophysical) images to optimize the acquisition scheme of the data needed to estimate those images. Thus, it can lead to high compression factors. A previous article by Bobin et al. (2008) showed how the new theory could be applied to simulated Herschel/PACS data to solve the compression requirement of the instrument. In this article, we show that compressed-sensing theory can indeed be successfully applied to actual Herschel/PACS data and give significant improvements over the standard pipeline. In order to fully use the redundancy present in the data, we perform full sky map estimation and decompression at the same time, which cannot be done in most other compression methods. We also demonstrate that the various artifacts affecting the data (pink noise, glitches, whose behavior is a priori not well compatible with compressed-sensing) can be handled as well in this new framework. Finally, we make a comparison between the methods from the compressed-sensing scheme and data acquired with the standard compression scheme. We discuss improvements that can be made on ground for the creation of sky maps from the data.Comment: 11 pages, 6 figures, 5 tables, peer-reviewed articl

The Disk and Dark Halo Mass of the Barred Galaxy NGC 4123. I. Observations

The non-circular streaming motions in barred galaxies are sensitive to the mass of the bar and can be used to lift the degeneracy between disk and dark matter halo encountered when fitting axisymmetric rotation curves of disk galaxies. In this paper, we present photometric and kinematic observations of NGC 4123, a barred galaxy of modest size (V_rot = 130 km/sec, L = 0.7 L_*), which reveal strong non-circular motions. The bar has straight dust lanes and an inner Lindblad resonance. The disk of NGC 4123 has no sign of truncation out to 10 scale lengths, and star-forming regions are found well outside R_25. A Fabry-Perot H-alpha velocity field shows velocity jumps of >100 km/sec at the location of the dust lanes within the bar, indicating shocks in the gas flow. VLA observations yield the velocity field of the H I disk. Axisymmetric mass models yield good fits to the rotation curve outside the bar regionfor disk I-band M/L of 2.25 or less, and dark halos with either isothermal or power-law profiles can fit the data well. In a companion paper, we model the full 2-D velocity field, including non-circular motions, to determine the stellar M/L and the mass of the dark halo.Comment: accepted by ApJ, 16 pages, 9 figures (1 color), uses emulateapj.sty, onecolfloat.st

Evidence That the P\u3csub\u3ei\u3c/sub\u3e Release Event Is the Rate-Limiting Step in the Nitrogenase Catalytic Cycle

Nitrogenase reduction of dinitrogen (N2) to ammonia (NH3) involves a sequence of events that occur upon the transient association of the reduced Fe protein containing two ATP molecules with the MoFe protein that includes electron transfer, ATP hydrolysis, Pi release, and dissociation of the oxidized, ADP-containing Fe protein from the reduced MoFe protein. Numerous kinetic studies using the nonphysiological electron donor dithionite have suggested that the rate-limiting step in this reaction cycle is the dissociation of the Fe protein from the MoFe protein. Here, we have established the rate constants for each of the key steps in the catalytic cycle using the physiological reductant flavodoxin protein in its hydroquinone state. The findings indicate that with this reductant, the rate-limiting step in the reaction cycle is not protein–protein dissociation or reduction of the oxidized Fe protein, but rather events associated with the Pi release step. Further, it is demonstrated that (i) Fe protein transfers only one electron to MoFe protein in each Fe protein cycle coupled with hydrolysis of two ATP molecules, (ii) the oxidized Fe protein is not reduced when bound to MoFe protein, and (iii) the Fe protein interacts with flavodoxin using the same binding interface that is used with the MoFe protein. These findings allow a revision of the rate-limiting step in the nitrogenase Fe protein cycle

The Progenitors of Type Ia Supernovae: Are They Supersoft Sources?

In a canonical model, the progenitors of Type Ia supernovae (SNe Ia) are accreting, nuclear-burning white dwarfs (NBWDs), which explode when the white dwarf reaches the Chandrasekhar mass, M_C. Such massive NBWDs are hot (kT ~100 eV), luminous (L ~ 10^{38} erg/s), and are potentially observable as luminous supersoft X-ray sources (SSSs). During the past several years, surveys for soft X-ray sources in external galaxies have been conducted. This paper shows that the results falsify the hypothesis that a large fraction of progenitors are NBWDs which are presently observable as SSSs. The data also place limits on sub-M_C models. While Type Ia supernova progenitors may pass through one or more phases of SSS activity, these phases are far shorter than the time needed to accrete most of the matter that brings them close to M_C.Comment: submitted to ApJ 18 November 2009; 17 pages, 2 figure

The Progenitors of Type Ia Supernovae: II. Are they Double-Degenerate Binaries? The Symbiotic Channel

In order for a white dwarf (WD) to achieve the Chandrasekhar mass, M_C, and explode as a Type Ia supernova (SNIa), it must interact with another star, either accreting matter from or merging with it. The failure to identify the types of binaries which produce SNeIa is the "progenitor problem". Its solution is required if we are to utilize the full potential of SNeIa to elucidate basic cosmological and physical principles. In single-degenerate models, a WD accretes and burns matter at high rates. Nuclear-burning WDs (NBWDs) with mass close to M_C are hot and luminous, potentially detectable as supersoft x-ray sources (SSSs). In previous work we showed that > 90-99% of the required number of progenitors do not appear as SSSs during most of the crucial phase of mass increase. The obvious implication is that double-degenerate (DD) binaries form the main class of progenitors. We show in this paper, however, that many binaries that later become DDs must pass through a long-lived NBWD phase during which they are potentially detectable as SSSs. The paucity of SSSs is therefore not a strong argument in favor of DD models. Those NBWDs that are the progenitors of DD binaries are likely to appear as symbiotic binaries for intervals > 10^6 years. In fact, symbiotic pre-DDs should be common, whether or not the WDs eventually produce SNeIa. The key to solving the progenitor problem lies in understanding the appearance of NBWDs. Most do not appear as SSSs most of the time. We therefore consider the evolution of NBWDs to address the question of what their appearance may be and how we can hope to detect them.Comment: 24 pages; 5 figures; submitted to Ap