Compression-based Image Registration

Image registration is an important component of image analysis used to align two or more images. In this paper, we present a new framework for image registration based on compression. The basic idea underlying our approach is the conjecture that two images are correctly registered when we can maximally compress one image given the information in the other. The contribution of this paper is twofold. First, we show that the image registration process can be dealt with from the perspective of a compression problem. Second, we demonstrate that the similarity metric, introduced by Li et al., performs well in image registration. Two different versions of the similarity metric have been used: the Kolmogorov version, computed using standard real-world compressors, and the Shannon version, calculated from an estimation of the entropy rate of the images

Cytotoxic effects of FGF2-saporin on bovine epithelial lens cells in vitro.

PURPOSE: To test the ability of two preparations of FGF2-saporin, either FGF2 chemically conjugated to saporin (FGF2-SAP) or genetically engineered FGF2-saporin (rFGF2-SAP) to inhibit the growth of bovine epithelial lens (BEL) cells in vitro when in solution and when immobilized on heparin surface-modified (HSM) polymethylmethacrylate (PMMA) intraocular lenses (IOLs). METHOD: Bovine epithelial lens cells were incubated with various concentrations FGF2-saporin for as long as 4 days. The number of surviving cells was determined by counting the number of nuclei. Because FGF2 binds to heparin, FGF2-saporin was incubated with HSM PMMA IOLs; excess toxin was washed off, and the BEL cells were grown on the FGF2-saporin-treated IOLs (HSM and non-HSM) for 4 days. Cell density was determined by image analysis. RESULTS: Both FGF2-SAP and rFGF2-SAP were highly cytotoxic (nM range), with rFGF2-SAP 10 times less active than FGF2-SAP. FGF2-saporin bound to the surface of HSM IOLs and eluted by 2M NaCl retained its activity. Toxin bound to HSM IOLs killed more than 90% of the BEL cells placed on the IOL surface within 4 days. The ability of FGF2-saporin to prevent the growth of cells on the IOL surface was strictly dependent on the presence of heparin on the IOL. CONCLUSIONS: FGF2-saporin is bound to HSM PMMA IOLs and prevents the growth of epithelial cells on the surface of the lens.

Comparative Study of Chemoembolization, Loadable Beads: In vitro Drug Release and Physical Properties of DC Bead and Hepasphere, Loaded with Doxorubicin and Irinotecan

Purpose To characterize in vitro the loadability, physical properties, and release of irinotecan and doxorubicin from two commercially available embolization microspheres. Materials and Methods DC Bead (500-700 μm) and Hepasphere (400-600 μm) microspheres were loaded with either doxorubicin or irinotecan solutions. Drug amount was quantified with spectrophotometry, bead elasticity was measured under compression, and bead size and loading homogeneity were assessed with microscopy image analysis. Drug release was measured over 1-week periods in saline by using a pharmacopeia flow-through method. Results Almost complete drug loading was obtained for both microsphere types and drugs. Doxorubicin-loaded DC Beads maintained their spherical shape throughout the release. In contrast, Hepaspheres showed less homogeneous doxorubicin loading and, after release, some fractured microspheres. Incomplete doxorubicin release was observed in saline over 1 week (27% ± 2 for DC beads and 18% ± 7 for Hepaspheres; P = .013). About 75% of this amount was released within 2.2 hours for both beads. For irinotecan, complete release was obtained for both types of beads, in a sustained manner over 2-3 hours for DC Beads, and in a significantly faster manner as a 7-minute burst for Hepaspheres. Conclusions The two drug-eluting microspheres could be efficiently loaded with both drugs. Incomplete doxorubicin release was attributed to strong drug-bead ionic interactions. Weaker interactions were observed with irinotecan, which led to faster drug release.

Automatic classification of MR scans in Alzheimer's disease.

To be diagnostically useful, structural MRI must reliably distinguish Alzheimer's disease (AD) from normal aging in individual scans. Recent advances in statistical learning theory have led to the application of support vector machines to MRI for detection of a variety of disease states. The aims of this study were to assess how successfully support vector machines assigned individual diagnoses and to determine whether data-sets combined from multiple scanners and different centres could be used to obtain effective classification of scans. We used linear support vector machines to classify the grey matter segment of T1-weighted MR scans from pathologically proven AD patients and cognitively normal elderly individuals obtained from two centres with different scanning equipment. Because the clinical diagnosis of mild AD is difficult we also tested the ability of support vector machines to differentiate control scans from patients without post-mortem confirmation. Finally we sought to use these methods to differentiate scans between patients suffering from AD from those with frontotemporal lobar degeneration. Up to 96% of pathologically verified AD patients were correctly classified using whole brain images. Data from different centres were successfully combined achieving comparable results from the separate analyses. Importantly, data from one centre could be used to train a support vector machine to accurately differentiate AD and normal ageing scans obtained from another centre with different subjects and different scanner equipment. Patients with mild, clinically probable AD and age/sex matched controls were correctly separated in 89% of cases which is compatible with published diagnosis rates in the best clinical centres. This method correctly assigned 89% of patients with post-mortem confirmed diagnosis of either AD or frontotemporal lobar degeneration to their respective group. Our study leads to three conclusions: Firstly, support vector machines successfully separate patients with AD from healthy aging subjects. Secondly, they perform well in the differential diagnosis of two different forms of dementia. Thirdly, the method is robust and can be generalized across different centres. This suggests an important role for computer based diagnostic image analysis for clinical practice.

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.

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.

A rapid and effective vignetting correction for quantitative microscopy

Images acquired using optical microscopes are inherently subject to vignetting effects due to imperfect illumination and image acquisition. However, such vignetting effects hamper accurate extraction of quantitative information from biological images, leading to less effective image segmentation and increased noise in the measurements. Here, we describe a rapid and effective method for vignetting correction, which generates an estimate for a correction function from the background fluorescence without the need to acquire additional calibration images. We validate the usefulness of this algorithm using artificially distorted images as a gold standard for assessing the accuracy of the applied correction and then demonstrate that this correction method enables the reliable detection of biologically relevant variation in cell populations. A simple user interface called FlattifY was developed and integrated into the image analysis platform YeastQuant to facilitate easy application of vignetting correction to a wide range of images.

Method to assess and optimise dependability of complex macro-systems: application to a railway signalling system

Achieving a high degree of dependability in complex macro-systems is challenging. Because of the large number of components and numerous independent teams involved, an overview of the global system performance is usually lacking to support both design and operation adequately. A functional failure mode, effects and criticality analysis (FMECA) approach is proposed to address the dependability optimisation of large and complex systems. The basic inductive model FMECA has been enriched to include considerations such as operational procedures, alarm systems. environmental and human factors, as well as operation in degraded mode. Its implementation on a commercial software tool allows an active linking between the functional layers of the system and facilitates data processing and retrieval, which enables to contribute actively to the system optimisation. The proposed methodology has been applied to optimise dependability in a railway signalling system. Signalling systems are typical example of large complex systems made of multiple hierarchical layers. The proposed approach appears appropriate to assess the global risk- and availability-level of the system as well as to identify its vulnerabilities. This enriched-FMECA approach enables to overcome some of the limitations and pitfalls previously reported with classical FMECA approaches.

Soil infiltration based on bp neural network and grey relational analysis

Soil infiltration is a key link of the natural water cycle process. Studies on soil permeability are conducive for water resources assessment and estimation, runoff regulation and management, soil erosion modeling, nonpoint and point source pollution of farmland, among other aspects. The unequal influence of rainfall duration, rainfall intensity, antecedent soil moisture, vegetation cover, vegetation type, and slope gradient on soil cumulative infiltration was studied under simulated rainfall and different underlying surfaces. We established a six factor-model of soil cumulative infiltration by the improved back propagation (BP)-based artificial neural network algorithm with a momentum term and self-adjusting learning rate. Compared to the multiple nonlinear regression method, the stability and accuracy of the improved BP algorithm was better. Based on the improved BP model, the sensitive index of these six factors on soil cumulative infiltration was investigated. Secondly, the grey relational analysis method was used to individually study grey correlations among these six factors and soil cumulative infiltration. The results of the two methods were very similar. Rainfall duration was the most influential factor, followed by vegetation cover, vegetation type, rainfall intensity and antecedent soil moisture. The effect of slope gradient on soil cumulative infiltration was not significant.

A post-processing pipeline to reconstruct the developing fetal brain using low-resolution MRI

Introduction. Development of the fetal brain surfacewith concomitant gyrification is one of the majormaturational processes of the human brain. Firstdelineated by postmortem studies or by ultrasound, MRIhas recently become a powerful tool for studying in vivothe structural correlates of brain maturation. However,the quantitative measurement of fetal brain developmentis a major challenge because of the movement of the fetusinside the amniotic cavity, the poor spatial resolution,the partial volume effect and the changing appearance ofthe developing brain. Today extensive efforts are made todeal with the âeurooepost-acquisitionâeuro reconstruction ofhigh-resolution 3D fetal volumes based on severalacquisitions with lower resolution (Rousseau, F., 2006;Jiang, S., 2007). We here propose a framework devoted tothe segmentation of the basal ganglia, the gray-whitetissue segmentation, and in turn the 3D corticalreconstruction of the fetal brain. Method. Prenatal MRimaging was performed with a 1-T system (GE MedicalSystems, Milwaukee) using single shot fast spin echo(ssFSE) sequences in fetuses aged from 29 to 32gestational weeks (slice thickness 5.4mm, in planespatial resolution 1.09mm). For each fetus, 6 axialvolumes shifted by 1 mm were acquired (about 1 min pervolume). First, each volume is manually segmented toextract fetal brain from surrounding fetal and maternaltissues. Inhomogeneity intensity correction and linearintensity normalization are then performed. A highspatial resolution image of isotropic voxel size of 1.09mm is created for each fetus as previously published byothers (Rousseau, F., 2006). B-splines are used for thescattered data interpolation (Lee, 1997). Then, basalganglia segmentation is performed on this superreconstructed volume using active contour framework witha Level Set implementation (Bach Cuadra, M., 2010). Oncebasal ganglia are removed from the image, brain tissuesegmentation is performed (Bach Cuadra, M., 2009). Theresulting white matter image is then binarized andfurther given as an input in the Freesurfer software(http://surfer.nmr.mgh.harvard.edu/) to provide accuratethree-dimensional reconstructions of the fetal brain.Results. High-resolution images of the cerebral fetalbrain, as obtained from the low-resolution acquired MRI,are presented for 4 subjects of age ranging from 29 to 32GA. An example is depicted in Figure 1. Accuracy in theautomated basal ganglia segmentation is compared withmanual segmentation using measurement of Dice similarity(DSI), with values above 0.7 considering to be a verygood agreement. In our sample we observed DSI valuesbetween 0.785 and 0.856. We further show the results ofgray-white matter segmentation overlaid on thehigh-resolution gray-scale images. The results arevisually checked for accuracy using the same principlesas commonly accepted in adult neuroimaging. Preliminary3D cortical reconstructions of the fetal brain are shownin Figure 2. Conclusion. We hereby present a completepipeline for the automated extraction of accuratethree-dimensional cortical surface of the fetal brain.These results are preliminary but promising, with theultimate goal to provide âeurooemovieâeuro of the normal gyraldevelopment. In turn, a precise knowledge of the normalfetal brain development will allow the quantification ofsubtle and early but clinically relevant deviations.Moreover, a precise understanding of the gyraldevelopment process may help to build hypotheses tounderstand the pathogenesis of several neurodevelopmentalconditions in which gyrification have been shown to bealtered (e.g. schizophrenia, autismâeuro¦). References.Rousseau, F. (2006), 'Registration-Based Approach forReconstruction of High-Resolution In Utero Fetal MR Brainimages', IEEE Transactions on Medical Imaging, vol. 13,no. 9, pp. 1072-1081. Jiang, S. (2007), 'MRI of MovingSubjects Using Multislice Snapshot Images With VolumeReconstruction (SVR): Application to Fetal, Neonatal, andAdult Brain Studies', IEEE Transactions on MedicalImaging, vol. 26, no. 7, pp. 967-980. Lee, S. (1997),'Scattered data interpolation with multilevel B-splines',IEEE Transactions on Visualization and Computer Graphics,vol. 3, no. 3, pp. 228-244. Bach Cuadra, M. (2010),'Central and Cortical Gray Mater Segmentation of MagneticResonance Images of the Fetal Brain', ISMRM Conference.Bach Cuadra, M. (2009), 'Brain tissue segmentation offetal MR images', MICCAI.

MORPHOLOGICAL AND MICROMORPHOLOGICAL CHANGES IN THE STRUCTURE OF A RHODIC HAPLUDOX AS A RESULT OF AGRICULTURAL MANAGEMENT

In evaluation of soil quality for agricultural use, soil structure is one of the most important properties, which is influenced not only by climate, biological activity, and management practices but also by mechanical and physico-chemical forces acting in the soil. The purpose of this study was to evaluate the influence of conventional agricultural management on the structure and microstructure of a Latossolo Vermelho distroférrico típico (Rhodic Hapludox) in an experimental area planted to maize. Soil morphology was described using the crop profile method by identifying the distinct structural volumes called Morphologically Homogeneous Units (MHUs). For comparison, we also described a profile in an adjacent area without agricultural use and under natural regrowth referred to as Memory. We took undisturbed samples from the main MHUs so as to form thin sections and blocks of soil for micromorphological and micromorphometrical analyses. Results from the application of the crop profile method showed the occurrence of the following structural types: loose (L), fragmented (F) and continuous (C) in both profiles analyzed. In the Memory soil profile, the fragmented structures were classified as Fptμ∆+tf and Fmt∆μ, whose micromorphology shows an enaulic-porphyric (porous) relative distribution with a great deal of biological activity as indicated by the presence of vughs and channels. Lower down, from 0.20 to 0.35 m, there was a continuous soil volume (sub-type C∆μ), with a subangular block microstructure and an enaulic-porphyric relative distribution, though in this case more compact and with aggregate coalescence and less biological activity. The micromorphometrical study of the soil of the Memory Plot showed the predominance of complex pores in NAM (15.03 %), Fmt∆μ (11.72 %), and Fptμ∆+tf (7.73 %), and rounded pores in C∆μ (8.21 %). In the soil under conventional agricultural management, we observed fragmented structures similar to the Memory Plot from 0.02 to 0.20 m, followed by a volume with a compact continuous structure (C∆μ), without visible porosity and with few roots. In the MHUs under conventional management, reduction in the packing pores (40 %) was observed, mainly in the continuous units (C). The microstructure had well-defined blocks, with the occurrence of planar pores and less evidence of biological activity. In conclusion, the morphological and micromorphological analyses of the soil profiles studied offered complementary information regarding soil structural quality, especially concerning the changes in pore types as result of mechanical stress undergone by the soil.

Enhanced oligonucleotide delivery to mouse retinal cells using iontophoresis.

PURPOSE: To study the combination of oligodeoxynucleotides (ODNs) intravitreous injection and saline transpalpebral iontophoresis on the delivery of ODNs to photoreceptors in the newborn rd1/rd1 mice. METHODS: Cathodal or anodal transpalpebral iontophoresis (1.43 mA/cm(2) for 5 min) was applied to eyes of postnatal day 7 (PN7) rd1/rd1 mice immediately before the intravitreous injection of ODNs. The effect of cathodal iontophoresis after ODNs injection was also evaluated. The influence of current intensity (0.5, 1.5, and 2.5 mA) was assayed with cathodal iontophoresis performed prior to ODNs injection. The duration of current-induced facilitation of ODNs delivery to photoreceptors was evaluated for 6 h following iontophoresis. One group of control eyes received cathodal iontophoresis prior to the intravitreous injection of phosphate buffered saline (PBS) or hexachlorofluorescein (Hex). The second control group received ODN or Hex intravitreous injection without iontophoresis. The penetration of fluorescent ODNs in the outer nuclear layer (ONL) was quantified by image analysis of the ONL fluorescence intensity on cryosection microphotographs. Integrity of ODN was assessed using acrylamide gel migration after its extraction from the retina of treated mice. The integrity of retinal structure, 1 and 24 h after iontophoresis, was analyzed using light and electron microscopy. RESULTS: Transpalpebral anodal or cathodal saline iontophoresis enhanced the penetration of ODNs in all retinal layers. Cathodal iontophoresis was more efficient than anodal iontophoresis in enhancing the tissue penetration of the injected ODN. Photoreceptor delivery of ODN was significantly higher when cathodal saline transpalpebral iontophoresis was applied prior than after the injection. The extent of enhanced tissue penetration decreased in parallel to the increased interval between iontophoresis application and the intravitreous injection. Current of 1.5 mA was safe and optimal for the delivery of ODNs to the ONL. One hour after iontophoresis followed by injection, ODN extracted from the retina of treated eyes remained intact. Histology and electron microscopy observations demonstrated that iontophoresis using the optimal parameters did not induce any permanent tissue alterations or structure damage. CONCLUSIONS: Saline transpalpebral iontophoresis facilitates the penetration of injected ODNs in photoreceptors for at least 3 h. This method may be considered for photoreceptor targeted gene therapy.

Photodynamic Diagnosis of Non-muscle-invasive Bladder Cancer with Hexaminolevulinate Cystoscopy: A Meta-analysis of Detection and Recurrence Based on Raw Data.

BACKGROUND: Studies on hexaminolevulinate (HAL) cystoscopy report improved detection of bladder tumours. However, recent meta-analyses report conflicting effects on recurrence. OBJECTIVE: To assess available clinical data for blue light (BL) HAL cystoscopy on the detection of Ta/T1 and carcinoma in situ (CIS) tumours, and on tumour recurrence. DESIGN, SETTING, AND PARTICIPANTS: This meta-analysis reviewed raw data from prospective studies on 1345 patients with known or suspected non-muscle-invasive bladder cancer (NMIBC). INTERVENTION: A single application of HAL cystoscopy was used as an adjunct to white light (WL) cystoscopy. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: We studied the detection of NMIBC (intention to treat [ITT]: n=831; six studies) and recurrence (per protocol: n=634; three studies) up to 1 yr. DerSimonian and Laird's random-effects model was used to obtain pooled relative risks (RRs) and associated 95% confidence intervals (CIs) for outcomes for detection. RESULTS AND LIMITATIONS: BL cystoscopy detected significantly more Ta tumours (14.7%; p<0.001; odds ratio [OR]: 4.898; 95% CI, 1.937-12.390) and CIS lesions (40.8%; p<0.001; OR: 12.372; 95% CI, 6.343-24.133) than WL. There were 24.9% patients with at least one additional Ta/T1 tumour seen with BL (p<0.001), significant also in patients with primary (20.7%; p<0.001) and recurrent cancer (27.7%; p<0.001), and in patients at high risk (27.0%; p<0.001) and intermediate risk (35.7%; p=0.004). In 26.7% of patients, CIS was detected only by BL (p<0.001) and was also significant in patients with primary (28.0%; p<0.001) and recurrent cancer (25.0%; p<0.001). Recurrence rates up to 12 mo were significantly lower overall with BL, 34.5% versus 45.4% (p=0.006; RR: 0.761 [0.627-0.924]), and lower in patients with T1 or CIS (p=0.052; RR: 0.696 [0.482-1.003]), Ta (p=0.040; RR: 0.804 [0.653-0.991]), and in high-risk (p=0.050) and low-risk (p=0.029) subgroups. Some subgroups had too few patients to allow statistically meaningful analysis. Heterogeneity was minimised by the statistical analysis method used. CONCLUSIONS: This meta-analysis confirms that HAL BL cystoscopy significantly improves the detection of bladder tumours leading to a reduction of recurrence at 9-12 mo. The benefit is independent of the level of risk and is evident in patients with Ta, T1, CIS, primary, and recurrent cancer.

Direct three-dimensional myocardial strain tensor quantification and tracking using zHARP.

Images of myocardial strain can be used to diagnose heart disease, plan and monitor treatment, and to learn about cardiac structure and function. Three-dimensional (3D) strain is typically quantified using many magnetic resonance (MR) images obtained in two or three orthogonal planes. Problems with this approach include long scan times, image misregistration, and through-plane motion. This article presents a novel method for calculating cardiac 3D strain using a stack of two or more images acquired in only one orientation. The zHARP pulse sequence encodes in-plane motion using MR tagging and out-of-plane motion using phase encoding, and has been previously shown to be capable of computing 3D displacement within a single image plane. Here, data from two adjacent image planes are combined to yield a 3D strain tensor at each pixel; stacks of zHARP images can be used to derive stacked arrays of 3D strain tensors without imaging multiple orientations and without numerical interpolation. The performance and accuracy of the method is demonstrated in vitro on a phantom and in vivo in four healthy adult human subjects.

Investigation of High Performance Concrete Pavement, MLR-05-01, 2006

The Iowa DOT has been using blended cements in ternary mixes since 1999. Use of these supplementary cementitious materials gives concrete with higher strengths and much lower permeability. Use of these materials has been incorporated for use in High Performance Concrete (HPC) decks to achieve lower permeability and thus long term performance. Since we have been using these materials in paving, it would be informative to determine what concrete pavement properties are enhanced as related to high performance concrete. The air void system was excellent at a spacing factor of 0.0047 in (0.120 mm). AVA spacing factor results are much higher than the hardened air void analysis. Although only 3 samples were tested between the image analysis air content and the RapidAir457, there is pretty good agreement between those test methods. Air void analysis indicates that excessive vibration was not required to place the concrete. Vibration was well within the specification limits with an average of 6683 vpm’s with a standard deviation of 461. Overall ride of the project was very good. The average smoothness for the project was 2.1 in/mile (33.8 mm/km). The International Roughness Index (IRI) was 81 in/mi (1.29 m/km). The compressive strength was 6260 psi (43.2 MPa) at 28 days and 6830 (47.1 MPa) at 56 days. The modulus of rupture by third point loading (MOR-TPL) tested at 28 days was 660 psi (4.55 MPa). The AASHTO T277 rapid chloride permeability results at 28 days using the Virginia cure method correlate fairly well with the 56 and 90 day results with standard curing. The Virginia cure method 28 day results were 2475 coulombs and the standard cure 56 and 90 day test results were 2180 and 2118, respectively.