Multiresolution-based image fusion with additive wavelet decomposition

The standard data fusion methods may not be satisfactory to merge a high-resolution panchromatic image and a low-resolution multispectral image because they can distort the spectral characteristics of the multispectral data. The authors developed a technique, based on multiresolution wavelet decomposition, for the merging and data fusion of such images. The method presented consists of adding the wavelet coefficients of the high-resolution image to the multispectral (low-resolution) data. They have studied several possibilities concluding that the method which produces the best results consists in adding the high order coefficients of the wavelet transform of the panchromatic image to the intensity component (defined as L=(R+G+B)/3) of the multispectral image. The method is, thus, an improvement on standard intensity-hue-saturation (IHS or LHS) mergers. They used the ¿a trous¿ algorithm which allows the use of a dyadic wavelet to merge nondyadic data in a simple and efficient scheme. They used the method to merge SPOT and LANDSATTM images. The technique presented is clearly better than the IHS and LHS mergers in preserving both spectral and spatial information.

Image restoration in astronomy: a Bayesian perspective

When preparing an article on image restoration in astronomy, it is obvious that some topics have to be dropped to keep the work at reasonable length. We have decided to concentrate on image and noise models and on the algorithms to find the restoration. Topics like parameter estimation and stopping rules are also commented on. We start by describing the Bayesian paradigm and then proceed to study the noise and blur models used by the astronomical community. Then the prior models used to restore astronomical images are examined. We describe the algorithms used to find the restoration for the most common combinations of degradation and image models. Then we comment on important issues such as acceleration of algorithms, stopping rules, and parameter estimation. We also comment on the huge amount of information available to, and made available by, the astronomical community.

Introduction of sensor spectral response into image fusion methods. Application to wavelet-based methods

Usual image fusion methods inject features from a high spatial resolution panchromatic sensor into every low spatial resolution multispectral band trying to preserve spectral signatures and improve spatial resolution to that of the panchromatic sensor. The objective is to obtain the image that would be observed by a sensor with the same spectral response (i.e., spectral sensitivity and quantum efficiency) as the multispectral sensors and the spatial resolution of the panchromatic sensor. But in these methods, features from electromagnetic spectrum regions not covered by multispectral sensors are injected into them, and physical spectral responses of the sensors are not considered during this process. This produces some undesirable effects, such as resolution overinjection images and slightly modified spectral signatures in some features. The authors present a technique which takes into account the physical electromagnetic spectrum responses of sensors during the fusion process, which produces images closer to the image obtained by the ideal sensor than those obtained by usual wavelet-based image fusion methods. This technique is used to define a new wavelet-based image fusion method.

Robust volume-targeted balanced steady-state free-precession coronary magnetic resonance angiography in a breathhold at 3.0 Tesla: a reproducibility study.

BACKGROUND: Transient balanced steady-state free-precession (bSSFP) has shown substantial promise for noninvasive assessment of coronary arteries but its utilization at 3.0 T and above has been hampered by susceptibility to field inhomogeneities that degrade image quality. The purpose of this work was to refine, implement, and test a robust, practical single-breathhold bSSFP coronary MRA sequence at 3.0 T and to test the reproducibility of the technique. METHODS: A 3D, volume-targeted, high-resolution bSSFP sequence was implemented. Localized image-based shimming was performed to minimize inhomogeneities of both the static magnetic field and the radio frequency excitation field. Fifteen healthy volunteers and three patients with coronary artery disease underwent examination with the bSSFP sequence (scan time = 20.5 ± 2.0 seconds), and acquisitions were repeated in nine subjects. The images were quantitatively analyzed using a semi-automated software tool, and the repeatability and reproducibility of measurements were determined using regression analysis and intra-class correlation coefficient (ICC), in a blinded manner. RESULTS: The 3D bSSFP sequence provided uniform, high-quality depiction of coronary arteries (n = 20). The average visible vessel length of 100.5 ± 6.3 mm and sharpness of 55 ± 2% compared favorably with earlier reported navigator-gated bSSFP and gradient echo sequences at 3.0 T. Length measurements demonstrated a highly statistically significant degree of inter-observer (r = 0.994, ICC = 0.993), intra-observer (r = 0.894, ICC = 0.896), and inter-scan concordance (r = 0.980, ICC = 0.974). Furthermore, ICC values demonstrated excellent intra-observer, inter-observer, and inter-scan agreement for vessel diameter measurements (ICC = 0.987, 0.976, and 0.961, respectively), and vessel sharpness values (ICC = 0.989, 0.938, and 0.904, respectively). CONCLUSIONS: The 3D bSSFP acquisition, using a state-of-the-art MR scanner equipped with recently available technologies such as multi-transmit, 32-channel cardiac coil, and localized B0 and B1+ shimming, allows accelerated and reproducible multi-segment assessment of the major coronary arteries at 3.0 T in a single breathhold. This rapid sequence may be especially useful for functional imaging of the coronaries where the acquisition time is limited by the stress duration and in cases where low navigator-gating efficiency prohibits acquisition of a free breathing scan in a reasonable time period.

Combined functional and viability cardiac MR imaging in a single breathhold

The combination of cardiac viability and functional information enhances the identification of different heart tissues in the setting of ischemic heart disease. A method has recently been proposed for obtaining black-blood delayed-enhancement (DE) viability images using the stimulated-echo acquisition mode (STEAM) MRI pulse sequence in a single short breathhold. The method was validated against conventional inversion-recovery (IR) DE images for identifying regions of myocardial infarction (MI). The method was based on the acquisition of three consecutive images of the same anatomical slice. One image has T(1)-weighted contrast in which infarction appears bright. The two other images are used to construct an anatomical image of the heart, which is combined with the first image to produce a black-blood viability image. However, using appropriate modulation and demodulation frequencies, the latter two images bear useful information about myocardial deformation that results in a cardiac strain-encoding (SENC) functional image. In this work, a method is proposed for obtaining three consecutive SENC images in a single acquisition that can be combined to produce a composite image of the heart, which shows both functional and viability information. The proposed technique reduces scan time by one-half, compared with separate acquisitions of functional and viability images, and alleviates misregistration problems caused by separate breathholds.

Aortic vessel wall magnetic resonance imaging at 3.0 Tesla: a reproducibility study of respiratory navigator gated free-breathing 3D black blood magnetic resonance imaging.

The purpose of this study was to evaluate a free-breathing three-dimensional (3D) dual inversion-recovery (DIR) segmented k-space gradient-echo (turbo field echo [TFE]) imaging sequence at 3T for the quantification of aortic vessel wall dimensions. The effect of respiratory motion suppression on image quality was tested. Furthermore, the reproducibility of the aortic vessel wall measurements was investigated. Seven healthy subjects underwent 3D DIR TFE imaging of the aortic vessel wall with and without respiratory navigator. Subsequently, this sequence with respiratory navigator was performed twice in 10 healthy subjects to test its reproducibility. The signal-to-noise (SNR), contrast-to-noise ratio (CNR), vessel wall sharpness, and vessel wall volume (VWV) were assessed. Data were compared using the paired t-test, and the reproducibility of VWV measurements was evaluated using intraclass correlation coefficients (ICCs). SNR, CNR, and vessel wall sharpness were superior in scans performed with respiratory navigator compared to scans performed without. The ICCs concerning intraobserver, interobserver, and interscan reproducibility were excellent (0.99, 0.94, and 0.95, respectively). In conclusion, respiratory motion suppression substantially improves image quality of 3D DIR TFE imaging of the aortic vessel wall at 3T. Furthermore, this optimized technique with respiratory motion suppression enables assessment of aortic vessel wall dimensions with high reproducibility.

Tandem repeat sequence variation as causative cis-eQTLs for protein-coding gene expression variation: the case of CSTB.

Association studies have revealed expression quantitative trait loci (eQTLs) for a large number of genes. However, the causative variants that regulate gene expression levels are generally unknown. We hypothesized that copy-number variation of sequence repeats contribute to the expression variation of some genes. Our laboratory has previously identified that the rare expansion of a repeat c.-174CGGGGCGGGGCG in the promoter region of the CSTB gene causes a silencing of the gene, resulting in progressive myoclonus epilepsy. Here, we genotyped the repeat length and quantified CSTB expression by quantitative real-time polymerase chain reaction in 173 lymphoblastoid cell lines (LCLs) and fibroblast samples from the GenCord collection. The majority of alleles contain either two or three copies of this repeat. Independent analysis revealed that the c.-174CGGGGCGGGGCG repeat length is strongly associated with CSTB expression (P = 3.14 × 10(-11)) in LCLs only. Examination of both genotyped and imputed single-nucleotide polymorphisms (SNPs) within 2 Mb of CSTB revealed that the dodecamer repeat represents the strongest cis-eQTL for CSTB in LCLs. We conclude that the common two or three copy variation is likely the causative cis-eQTL for CSTB expression variation. More broadly, we propose that polymorphic tandem repeats may represent the causative variation of a fraction of cis-eQTLs in the genome.

The tomato genome sequence providies insights into fleshy fruit evolution

Tomato (Solanum lycopersicum) is a major crop plant and a model system for fruit development. Solanum is one of the largest angiosperm genera1 and includes annual and perennial plants from diverse habitats. Here we present a high-quality genome sequence of domesticated tomato, a draft sequence of its closest wild relative, Solanum pimpinellifolium2, and compare them to each other and to the potato genome (Solanum tuberosum). The two tomato genomes show only 0.6% nucleotide divergence and signs of recent admixture, but show more than 8% divergence from potato, with nine large and several smaller inversions. In contrast to Arabidopsis, but similar to soybean, tomato and potato small RNAs map predominantly to gene-rich chromosomal regions, including gene promoters. The Solanum lineage has experienced two consecutive genome triplications: one that is ancient and shared with rosids, and a more recent one. These triplications set the stage for the neofunctionalization of genes controlling fruit characteristics, such as colour and fleshiness.

Semi-supervised remote sensing image classification with cluster kernels

A semisupervised support vector machine is presented for the classification of remote sensing images. The method exploits the wealth of unlabeled samples for regularizing the training kernel representation locally by means of cluster kernels. The method learns a suitable kernel directly from the image and thus avoids assuming a priori signal relations by using a predefined kernel structure. Good results are obtained in image classification examples when few labeled samples are available. The method scales almost linearly with the number of unlabeled samples and provides out-of-sample predictions.

Mesozoic sequence of Fuerteventura (Canary Islands): Witness of Early Jurassic sea-floor spreading in the central Atlantic

The Fuerteventura Jurassic sedimentary succession consists of oceanic and elastic deposits, the latter derived from the southwestern Moroccan continental margin. Normal mid-oceanic-ridge basalt (N-MORB) flows and breccias are found at the base of the sequence and witness sea-floor spreading events in the central Atlantic. These basalts were extruded in a postrift environment (post-late Pliensbachian), We propose a Toarcian age for the Atlantic oceanic floor in this region, on the basis of the presence higher up in the sequence of the Bositra buchi filament microfacies (Aalenian-Bajocian) and of elastic deposits reflecting tectono-eustatic events (e.g,, late Toarcian to mid-Callovian erosion of the rift shoulder). The S-l sea-floor oceanic magnetic anomaly west of Fuerteventura is therefore at least Toarcian in age. The remaining sequence records Atlantic-Tethyan basinal facies (e.g., Callovian-Oxfordian red clays, Aptian-Albian black shales) alternating with elastic deposits (e.g., Kimmeridgian-Berriasian periplatform calciturbidites and a Lower Cretaceous deep-sea fan system). The Fuerteventura N-MORB outcrops represent the only Early Jurassic oceanic basement described so far in the central Atlantic. They are covered by a 1600 m, nearly continuous sedimentary sequence which extends to Upper Cretaceous facies.

Functional analysis of the 5' flanking sequence of a vaccinia virus late gene.

A series of mutations, including 5' and 3' deletions, as well as insertions were introduced into the 5' flanking nucleotide sequence of a vaccinia virus late gene. This DNA has been shown previously to contain all the necessary elements for correct regulation of the gene most probably transcribed by the viral RNA polymerase. To facilitate the assays, the mutated DNA was fused to the chloramphenicol acetyltransferase gene and inserted into the genome of live vaccinia virus. The effects of the mutations on expression of the chimeric gene were studied by both enzyme assays and nuclease S1 analysis. The results showed that 5' deletions up to about 15 bp from the putative initiation site of transcription still yielded high levels of gene expression. All mutations, however, that deleted the authentic late mRNA start site, abolished promoter activity.

Transfer insensitive labeling technique (TILT): application to multislice functional perfusion imaging.

Cerebral blood flow can be studied in a multislice mode with a recently proposed perfusion sequence using inversion of water spins as an endogenous tracer without magnetization transfer artifacts. The magnetization transfer insensitive labeling technique (TILT) has been used for mapping blood flow changes at a microvascular level under motor activation in a multislice mode. In TILT, perfusion mapping is achieved by subtraction of a perfusion-sensitized image from a control image. Perfusion weighting is accomplished by proximal blood labeling using two 90 degrees radiofrequency excitation pulses. For control preparation the labeling pulses are modified such that they have no net effect on blood water magnetization. The percentage of blood flow change, as well as its spatial extent, has been studied in single and multislice modes with varying delays between labeling and imaging. The average perfusion signal change due to activation was 36.9 +/- 9.1% in the single-slice experiments and 38.1 +/- 7.9% in the multislice experiments. The volume of activated brain areas amounted to 1.51 +/- 0.95 cm3 in the contralateral primary motor (M1) area, 0.90 +/- 0.72 cc in the ipsilateral M1 area, 1.27 +/- 0.39 cm3 in the contralateral and 1.42 +/- 0.75 cm3 in the ipsilateral premotor areas, and 0.71 +/- 0.19 cm3 in the supplementary motor area.

One hundred base pairs of 5' flanking sequence of a vaccinia virus late gene are sufficient to temporally regulate late transcription.

A vaccinia virus late gene coding for a major structural polypeptide of 11 kDa was sequenced. Although the 5' flanking gene region is very A+T rich, it shows little homology either to the corresponding region of vaccinia early genes or to consensus sequences characteristic of most eukaryotic genes. Three DNA fragments (100, 200, and 500 base pairs, respectively), derived from the flanking region and including the late gene mRNA start site, were inserted into the coding sequence of the vaccinia virus thymidine kinase (TK) early gene by homologous in vivo recombination. Recombinants were selected on the basis of their TK- phenotype. Cells were infected with the recombinant viruses and RNA was isolated at 1-hr intervals. Transcripts initiating either from the TK early promoter, or from the late gene promoter at its authentic position, or from the translocated late gene promoters within the early gene were detected by nuclease S1 mapping. Early after infection, only transcripts from the TK early promoter were detected. Later in infection, however, transcripts were also initiated from the translocated late promoters. This RNA appeared at the same time and in similar quantities as the RNA from the late promoter at its authentic position. No quantitative differences in promoter efficiency between the 100-, 200-, and 500-base-pair insertions were observed. We conclude that all necessary signals for correct regulation of late-gene expression reside within only 100 base pairs of 5' flanking sequence.