Methodology of image analysis for study of the vertisols moisture content

The main problem to study vertical drainage from the moisture distribution, on a vertisol profile, is searching for suitable methods using these procedures. Our aim was to design a digital image processing methodology and its analysis to characterize the moisture content distribution of a vertisol profile. In this research, twelve soil pits were excavated on a ba re Mazic Pellic Vertisols ix of them in May 13/2011 and the rest in May 19 /2011 after a moderate rainfall event. Digital RGB images were taken from each vertisol pit using a Kodak? camera selecting a size of 1600x945 pixels. Each soil image was processed to homogenized brightness and then a spatial filter with several window sizes was applied to select the optimum one. The RGB image obtained were divided in each matrix color selecting the best thresholds for each one, maximum and minimum, to be applied and get a digital binary pattern. This one was analyzed by estimating two fractal scaling exponents box counting dimension D BC) and interface fractal dimension (D) In addition, three pre-fractal scaling coefficients were determinate at maximum resolution: total number of boxes intercepting the foreground pattern (A), fractal lacunarity (?1) and Shannon entropy S1). For all the images processed the spatial filter 9x9 was the optimum based on entropy, cluster and histogram criteria. Thresholds for each color were selected based on bimodal histograms.

Diagnostic ultrasound activation of contrast agent gas bodies induces capillary rupture in mice

Interaction of diagnostic ultrasound with gas bodies produces a useful contrast effect in medical images, but the same interaction also represents a mechanism for bioeffects. Anesthetized hairless mice were scanned by using a 2.5-MHz transducer (610-ns pulses with 3.6-kHz repetition frequency and 61-Hz frame rate) after injection of Optison and Evans blue dye. Petechial hemorrhages (PHs) in intestine and abdominal muscle were counted 15 min after exposure to characterize capillary rupture, and Evans blue extravasation was evaluated in samples of muscle tissue. For 5 ml⋅kg-1 contrast agent and exposure to 10 alternating 10-s on and off periods, PH counts in muscle were approximately proportional to the square of peak negative pressure amplitude and were statistically significant above 0.64 MPa. PH counts in intestine and Evans blue extravasation into muscle tissue were significant above 1.0 MPa. The PH effect in muscle was proportional to contrast dose and was statistically significant for the lowest dose of 0.05 ml⋅kg-1. The effects decreased nearly to sham levels if the exposure was delayed 5 min. The PH effect in abdominal muscle was significant and statistically indistinguishable for uninterrupted 100-s exposure, 10-s exposure, 100 scans repeated at 1 Hz, and even for a single scan. The results confirms a previous report of PH induction by diagnostic ultrasound with contrast agent in mammalian skeletal muscle [Skyba, D. M., Price, R. J., Linka, A. Z., Skalak, T. C. & Kaul, S. (1998) Circulation 98, 290–293].

Appetite of an epiphyte: Quantitative monitoring of bacterial sugar consumption in the phyllosphere

We report here the construction, characterization, and application of a bacterial bioreporter for fructose and sucrose that was designed to monitor the availability of these sugars to microbial colonizers of the phyllosphere. Plasmid pPfruB-gfp[AAV] carries the Escherichia coli fruB promoter upstream from the gfp[AAV] allele that codes for an unstable variant of green fluorescent protein (GFP). In Erwinia herbicola, this plasmid brings about the accumulation of GFP fluorescence in response to both fructose and sucrose. Cells of E. herbicola (pPfruB-gfp[AAV]) were sprayed onto bean plants, recovered from leaves at various time intervals after inoculation, and analyzed individually for GFP content by quantitative analysis of digital microscope images. We observed a positive correlation between single-cell GFP accumulation and ribosomal content as determined by fluorescence in situ hybridization, indicating that foliar growth of E. herbicola occurred at the expense of fructose and/or sucrose. One hour after inoculation, nearly all bioreporter cells appeared to be actively engaged in fructose consumption. This fraction dropped to approximately 11% after 7 h and to ≈1% a day after inoculation. This pattern suggests a highly heterogeneous availability of fructose to individual E. herbicola cells as they colonize the phyllosphere. We estimated that individual cells were exposed to local initial fructose abundances ranging from less than 0.15 pg fructose to more than 4.6 pg.

An O(N log N) algorithm for shape modeling.

We present a shape-recovery technique in two dimensions and three dimensions with specific applications in modeling anatomical shapes from medical images. This algorithm models extremely corrugated structures like the brain, is topologically adaptable, and runs in O(N log N) time, where N is the total number of points in the domain. Our technique is based on a level set shape-recovery scheme recently introduced by the authors and the fast marching method for computing solutions to static Hamilton-Jacobi equations.

bioNMF: a versatile tool for non-negative matrix factorization in biology

Medical imaging has become an absolutely essential diagnostic tool for clinical practices; at present, pathologies can be detected with an earliness never before known. Its use has not only been relegated to the field of radiology but also, increasingly, to computer-based imaging processes prior to surgery. Motion analysis, in particular, plays an important role in analyzing activities or behaviors of live objects in medicine. This short paper presents several low-cost hardware implementation approaches for the new generation of tablets and/or smartphones for estimating motion compensation and segmentation in medical images. These systems have been optimized for breast cancer diagnosis using magnetic resonance imaging technology with several advantages over traditional X-ray mammography, for example, obtaining patient information during a short period. This paper also addresses the challenge of offering a medical tool that runs on widespread portable devices, both on tablets and/or smartphones to aid in patient diagnostics.

Simulação computacional de sistemas de registro de imagens mamográficas

Neste projeto foi desenvolvido um método computacional para verificação da melhor combinação tela intensificadora - filme para exames mamográficos através do estudo de suas características sensitométricas. O software, desenvolvido em ambiente Delphi para windows, apresenta na tela do microcomputador a imagem a ser obtida para cada tipo de combinação tela intensificadora - filme, utilizando imagens de \"Phantoms\" e de mamas reais. Em razão da ampla quantidade de fatores que influenciam a imagem mamográfica final, tais como magnificação, característica dos filmes e telas intensificadoras e condições da processadora, o método proposto pode proporcionar uma ampla avaliação da qualidade dos sistemas de imagem mamográfica de uma forma simples, rápida e automática, através de procedimentos de simulação computacional. A simulação investigou a influência que um determinado sistema de registro exerce sobre a qualidade da imagem, possibilitando conhecer previamente a imagem final a ser obtida com diferentes equipamentos e sistemas de registro. Dentre os sistemas investigados, três filmes (Kodak Min R 2000, Fuji UM MA-HC e Fuji ADM) e duas telas intensificadoras (Kodak Min R 2000 e Fuji AD Mammo Fine), aquele que apresentou melhores resultados, com melhor qualidade de imagens e menor exposição à paciente foi o de tela Min R 2000 com filme Min R 2000 da Kodak.

Ice rafted debris distribution in sediment core CRP-1

Lonestone abundances in CRP-1 were investigated using three methods: core examination at Cape Roberts Camp, analysis of digital core images and follow-up core examination. For all images of split-core, we determined size and depth of every detectable lonestone larger than 3 mm. Lonestone abundance decreases exponentially with clast size. Although no significant depth-dependent variations in lonestone size distribution were detected, a strong 0.5-0.7 m abundance periodicity, of unknown origin, is evident within diamicts. Lonestone volume percentage was estimated from size distribution: most size classes contribute approximately the same volume to the total. Sizes >16 mm have rare enough lonestones that their counts are nonrepresentative when based on short intervals of split core. This problem does not affect total counts significantly, but the volume analysis needs to be confined to <= 6 mm lonestones to avoid instability induced by rare and nonrepresentative larger lonestones. If lonestone abundance can be used as an indicator of glacial proximity, then our CRP-1 lonestone abundance logs confirm the overall character of previously inferred variations in relative distance to the ice margin. Large-scale changes in lonestone abundance also reflect the CRP-1 sequence stratigraphy, with individual sequences generally characterised by basal lonestone-rich diamict overlain by lonestone-poor sands and muds. The relationship between glacial proximity and lonestone abundance within diamicts and within sand-mud intervals is, however, less certain. For example, two or three gradual lonestone increases may indicate regressions during glacial advances, in contrast to the more common CRP-l pattern of dominantly transgressive sequences.

Sea ice habitat characteristics and number of narwhal (Monodon monoceros) sightings in Disko Bay, West Greenland

There is a paucity of information on abundance, densities, and habitat selection of narwhals Monodon monoceros in the offshore pack ice of Baffin Bay, West Greenland, despite the critical importance of winter foraging regions and considerable sea ice declines in the past decades. We conducted a double-platform visual aerial survey over a narwhal wintering ground to obtain pack ice densities and develop the first fully corrected abundance estimate using point conditional mark-recapture distance sampling. Continuous video recording and digital images taken along the trackline allowed for in situ quantification of winter narwhal habitat and for the estimation of fine-scale narwhal habitat selection and habitat-specific sighting probabilities. Abundance at the surface was estimated at 3484 (coefficient of variation [CV] = 0.46) including whales missed by observers. The fully corrected abundance of narwhals was 18 044 (CV = 0.46), or approximately one-quarter of the entire Baffin Bay population. The narwhal wintering ground surveyed (~9500 km**2) had 2.4 to 3.2% open water based on estimates from satellite imagery (NASA Moderate Resolution Imaging Spectroradiometer) and 1565 digital photographic images collected on the trackline. Thus, the ~18 000 narwhals had access to 233 km**2 of open water, resulting in an average density of ~77 narwhals/km**2 open water. Narwhal sighting probability near habitats with <10% or 10 to 50% open water was significantly higher than sighting probability in habitats with >50% open water, suggesting narwhals select optimal foraging areas in dense pack ice regardless of open water availability. This study provides the first quantitative ecological data on densities and habitat selection of narwhals in pack ice foraging regions that are rapidly being altered with climate change.

Capillary break-up rheometry of low-viscosity elastic fluids

We investigate the dynamics of the capillary thinning and break-up process for low viscosity elastic fluids such as dilute polymer solutions. Standard measurements of the evolution of the midpoint diameter of the necking fluid filament are augmented by high speed digital video images of the break up dynamics. We show that the successful operation of a capillary thinning device is governed by three important time scales (which characterize the relative importance of inertial, viscous and elastic processes), and also by two important length scales (which specify the initial sample size and the total stretch imposed on the sample). By optimizing the ranges of these geometric parameters, we are able to measure characteristic time scales for tensile stress growth as small as 1 millisecond for a number of model dilute and semi-dilute solutions of polyethylene oxide (PEO) in water and glycerol. If the final aspect ratio of the sample is too small, or the total axial stretch is too great, measurements are limited, respectively, by inertial oscillations of the liquid bridge or by the development of the well-known beads-on-a-string morphology which disrupt the formation of a uniform necking filament. By considering the magnitudes of the natural time scales associated with viscous flow, elastic stress growth and inertial oscillations it is possible to construct an operability diagram characterizing successful operation of a capillary break-up extensional rheometer. For Newtonian fluids, viscosities greater than approximately 70 mPas are required; however for dilute solutions of high molecular weight polymer, the minimum Viscosity is substantially lower due to the additional elastic stresses arising from molecular extension. For PEO of molecular weight 2.10(6) g/mol, it is possible to measure relaxation times of order 1 ms in dilute polymer solutions with zero-shear-rate viscosities on the order of 2-10 mPas.

Automatic initialisation of 3D deformable models for cartilage segmentation

Deformable models are a highly accurate and flexible approach to segmenting structures in medical images. The primary drawback of deformable models is that they are sensitive to initialisation, with accurate and robust results often requiring initialisation close to the true object in the image. Automatically obtaining a good initialisation is problematic for many structures in the body. The cartilages of the knee are a thin elastic material that cover the ends of the bone, absorbing shock and allowing smooth movement. The degeneration of these cartilages characterize the progression of osteoarthritis. The state of the art in the segmentation of the cartilage are 2D semi-automated algorithms. These algorithms require significant time and supervison by a clinical expert, so the development of an automatic segmentation algorithm for the cartilages is an important clinical goal. In this paper we present an approach towards this goal that allows us to automatically providing a good initialisation for deformable models of the patella cartilage, by utilising the strong spatial relationship of the cartilage to the underlying bone.

Objective analysis of toric intraocular lens rotation and centration

PURPOSE: To assess the repeatability of an objective image analysis technique to determine intraocular lens (IOL) rotation and centration. SETTING: Six ophthalmology clinics across Europe. METHODS: One-hundred seven patients implanted with Akreos AO aspheric IOLs with orientation marks were imaged. Image quality was rated by a masked observer. The axis of rotation was determined from a line bisecting the IOL orientation marks. This was normalized for rotation of the eye between visits using the axis bisecting 2 consistent conjunctival vessels or iris features. The center of ovals overlaid to circumscribe the IOL optic edge and the pupil or limbus were compared to determine IOL centration. Intrasession repeatability was assessed in 40 eyes and the variability of repeated analysis examined. RESULTS: Intrasession rotational stability of the IOL was ±0.79 degrees (SD) and centration was ±0.10 mm horizontally and ±0.10 mm vertically. Repeated analysis variability of the same image was ±0.70 degrees for rotation and ±0.20 mm horizontally and ±0.31 mm vertically for centration. Eye rotation (absolute) between visits was 2.23 ± 1.84 degrees (10%>5 degrees rotation) using one set of consistent conjunctival vessels or iris features and 2.03 ± 1.66 degrees (7%>5 degrees rotation) using the average of 2 sets (P =.13). Poorer image quality resulted in larger apparent absolute IOL rotation (r =-0.45,P<.001). CONCLUSIONS: Objective analysis of digital retroillumination images allows sensitive assessment of IOL rotation and centration stability. Eye rotation between images can lead to significant errors if not taken into account. Image quality is important to analysis accuracy.

Rotational and centration stability of an aspheric intraocular lens with a simulated toric design

Purpose: To assess the stability of the Akreos AO intraocular lens (IOL) platform with a simulated toric design using objective image analysis. Setting: Six hospital eye clinics across Europe. Methods: After implantation in 1 eye of patients, IOLs with orientation marks were imaged at 1 to 2 days, 7 to 14 days, 30 to 60 days, and 120 to 180 days. The axis of rotation and IOL centration were objectively assessed using validated image analysis. Results: The study enrolled 107 patients with a mean age of 69.9 years ± 7.7 (SD). The image quality was sufficient for IOL rotation analysis in 91% of eyes. The mean rotation between the first day postoperatively and 120 to 180 days was 1.93 ± 2.33 degrees, with 96% of IOLs rotating fewer than 5 degrees and 99% rotating fewer than 10 degrees. There was no significant rotation between visits and no clear bias in the direction of rotation. In 71% of eyes, the dilation and image quality was sufficient for image analysis of centration. The mean change in centration between 1 day and 120 to 180 days was 0.21 ± 0.11 mm, with all IOLs decentering less than 0.5 mm. There was no significant decentration between visits and no clear bias in the direction of the decentration. Conclusion: Objective analysis of digital retroillumination images taken at different postoperative periods shows the aspheric IOL platform was stable in the eye and is therefore suitable for the application of a toric surface to correct corneal astigmatism.

Tissue Equivalent Phantom Design for Optimization of a Coherent Scatter Imaging System

Scatter in medical imaging is typically cast off as image-related noise that detracts from meaningful diagnosis. It is therefore typically rejected or removed from medical images. However, it has been found that every material, including cancerous tissue, has a unique X-ray coherent scatter signature that can be used to identify the material or tissue. Such scatter-based tissue-identification provides the advantage of locating and identifying particular materials over conventional anatomical imaging through X-ray radiography. A coded aperture X-ray coherent scatter spectral imaging system has been developed in our group to classify different tissue types based on their unique scatter signatures. Previous experiments using our prototype have demonstrated that the depth-resolved coherent scatter spectral imaging system (CACSSI) can discriminate healthy and cancerous tissue present in the path of a non-destructive x-ray beam. A key to the successful optimization of CACSSI as a clinical imaging method is to obtain anatomically accurate phantoms of the human body. This thesis describes the development and fabrication of 3D printed anatomical scatter phantoms of the breast and lung.

The purpose of this work is to accurately model different breast geometries using a tissue equivalent phantom, and to classify these tissues in a coherent x-ray scatter imaging system. Tissue-equivalent anatomical phantoms were designed to assess the capability of the CACSSI system to classify different types of breast tissue (adipose, fibroglandular, malignant). These phantoms were 3D printed based on DICOM data obtained from CT scans of prone breasts. The phantoms were tested through comparison of measured scatter signatures with those of adipose and fibroglandular tissue from literature. Tumors in the phantom were modeled using a variety of biological tissue including actual surgically excised benign and malignant tissue specimens. Lung based phantoms have also been printed for future testing. Our imaging system has been able to define the location and composition of the various materials in the phantom. These phantoms were used to characterize the CACSSI system in terms of beam width and imaging technique. The result of this work showed accurate modeling and characterization of the phantoms through comparison of the tissue-equivalent form factors to those from literature. The physical construction of the phantoms, based on actual patient anatomy, was validated using mammography and computed tomography to visually compare the clinical images to those of actual patient anatomy.

Improving the Accuracy of 3-D Reconstruction in Robotic Vision Applications

Thesis (Ph.D.)--University of Washington, 2016-08

An algoritmh for the estimation of 3D surfaces using competitive learning

This paper addresses the estimation of object boundaries from a set of 3D points. An extension of the constrained clustering algorithm developed by Abrantes and Marques in the context of edge linking is presented. The object surface is approximated using rectangular meshes and simplex nets. Centroid-based forces are used for attracting the model nodes towards the data, using competitive learning methods. It is shown that competitive learning improves the model performance in the presence of concavities and allows to discriminate close surfaces. The proposed model is evaluated using synthetic data and medical images (MRI and ultrasound images).