Robustness of speckle imaging techniques applied to horizontal imaging scenarios

Atmospheric turbulence near the ground severely limits the quality of imagery acquired over long horizontal paths. In defense, surveillance, and border security applications, there is interest in deploying man-portable, embedded systems incorporating image reconstruction methods to compensate turbulence effects. While many image reconstruction methods have been proposed, their suitability for use in man-portable embedded systems is uncertain. To be effective, these systems must operate over significant variations in turbulence conditions while subject to other variations due to operation by novice users. Systems that meet these requirements and are otherwise designed to be immune to the factors that cause variation in performance are considered robust. In addition robustness in design, the portable nature of these systems implies a preference for systems with a minimum level of computational complexity. Speckle imaging methods have recently been proposed as being well suited for use in man-portable horizontal imagers. In this work, the robustness of speckle imaging methods is established by identifying a subset of design parameters that provide immunity to the expected variations in operating conditions while minimizing the computation time necessary for image recovery. Design parameters are selected by parametric evaluation of system performance as factors external to the system are varied. The precise control necessary for such an evaluation is made possible using image sets of turbulence degraded imagery developed using a novel technique for simulating anisoplanatic image formation over long horizontal paths. System performance is statistically evaluated over multiple reconstruction using the Mean Squared Error (MSE) to evaluate reconstruction quality. In addition to more general design parameters, the relative performance the bispectrum and the Knox-Thompson phase recovery methods is also compared. As an outcome of this work it can be concluded that speckle-imaging techniques are robust to the variation in turbulence conditions and user controlled parameters expected when operating during the day over long horizontal paths. Speckle imaging systems that incorporate 15 or more image frames and 4 estimates of the object phase per reconstruction provide up to 45% reduction in MSE and 68% reduction in the deviation. In addition, Knox-Thompson phase recover method is shown to produce images in half the time required by the bispectrum. The quality of images reconstructed using Knox-Thompson and bispectrum methods are also found to be nearly identical. Finally, it is shown that certain blind image quality metrics can be used in place of the MSE to evaluate quality in field scenarios. Using blind metrics rather depending on user estimates allows for reconstruction quality that differs from the minimum MSE by as little as 1%, significantly reducing the deviation in performance due to user action.

Comprehensive review and application of particle image velocimetry

For a fluid dynamics experimental flow measurement technique, particle image velocimetry (PIV) provides significant advantages over other measurement techniques in its field. In contrast to temperature and pressure based probe measurements or other laser diagnostic techniques including laser Doppler velocimetry (LDV) and phase Doppler particle analysis (PDPA), PIV is unique due to its whole field measurement capability, non-intrusive nature, and ability to collect a vast amount of experimental data in a short time frame providing both quantitative and qualitative insight. These properties make PIV a desirable measurement technique for studies encompassing a broad range of fluid dynamics applications. However, as an optical measurement technique, PIV also requires a substantial technical understanding and application experience to acquire consistent, reliable results. Both a technical understanding of particle image velocimetry and practical application experience are gained by applying a planar PIV system at Michigan Technological University’s Combustion Science Exploration Laboratory (CSEL) and Alternative Fuels Combustion Laboratory (AFCL). Here a PIV system was applied to non-reacting and reacting gaseous environments to make two component planar PIV as well as three component stereographic PIV flow field velocity measurements in conjunction with chemiluminescence imaging in the case of reacting flows. This thesis outlines near surface flow field characteristics in a tumble strip lined channel, three component velocity profiles of non-reacting and reacting swirled flow in a swirl stabilized lean condition premixed/prevaporized-fuel model gas turbine combustor operating on methane at 5-7 kW, and two component planar PIV measurements characterizing the AFCL’s 1.1 liter closed combustion chamber under dual fan driven turbulent mixing flow.

Ultraviolet digital imaging of volcanic plumes : implementation and application to magmatic processes at basaltic volcanoes

Magmatic volatiles play a crucial role in volcanism, from magma production at depth to generation of seismic phenomena to control of eruption style. Accordingly, many models of volcano dynamics rely heavily on behavior of such volatiles. Yet measurements of emission rates of volcanic gases have historically been limited, which has restricted model verification to processes on the order of days or longer. UV cameras are a recent advancement in the field of remote sensing of volcanic SO2 emissions. They offer enhanced temporal and spatial resolution over previous measurement techniques, but need development before they can be widely adopted and achieve the promise of integration with other geophysical datasets. Large datasets require a means by which to quickly and efficiently use imagery to calculate emission rates. We present a suite of programs designed to semi-automatically determine emission rates of SO2 from series of UV images. Extraction of high temporal resolution SO2 emission rates via this software facilitates comparison of gas data to geophysical data for the purposes of evaluating models of volcanic activity and has already proven useful at several volcanoes. Integrated UV camera and seismic measurements recorded in January 2009 at Fuego volcano, Guatemala, provide new insight into the system’s shallow conduit processes. High temporal resolution SO2 data reveal patterns of SO2 emission rate relative to explosions and seismic tremor that indicate tremor and degassing share a common source process. Progressive decreases in emission rate appear to represent inhibition of gas loss from magma as a result of rheological stiffening in the upper conduit. Measurements of emission rate from two closely-spaced vents, made possible by the high spatial resolution of the camera, help constrain this model. UV camera measurements at Kilauea volcano, Hawaii, in May of 2010 captured two occurrences of lava filling and draining within the summit vent. Accompanying high lava stands were diminished SO2 emission rates, decreased seismic and infrasonic tremor, minor deflation, and slowed lava lake surface velocity. Incorporation of UV camera data into the multi-parameter dataset gives credence to the likelihood of shallow gas accumulation as the cause of such events.

INVESTIGATING THE ROLE OF THE CONTACT LINE IN HETEROGENEOUS NUCLEATION WITH HIGH SPEED IMAGING

While nucleation of solids in supercooled liquids is ubiquitous [15, 65, 66], surface crystallization, the tendency for freezing to begin preferentially at the liquid-gas interface, has remained puzzling [74, 18, 68, 69, 51, 64, 72, 16]. Here we employ high-speed imaging of supercooled water drops to study the phenomenon of heterogeneous surface crystallization. Our geometry avoids the "point-like contact" of prior experiments by providing a simple, symmetric contact line (triple line defined by the substrate-liquid-air interface) for a drop resting on a homogeneous silicon substrate. We examine three possible mechanisms that might explain these laboratory observations: (i) Line Tension at the triple line, (ii) Thermal Gradients within the droplets and (iii) Surface Texture. In our first study we record nearly perfect spatial uniformity in the immersed (liquid-substrate) region and, thereby, no preference for nucleation at the triple line. In our second study, no influence of thermal gradients on the preference for freezing at the triple line was observed. Motivated by the conjectured importance of line tension (τ) [1, 66] for heterogeneous nucleation, we also searched for evidence of a transition to surface crystallization at length scales on the order of δ ∼ τ/σ, where σ is the surface tension [14]; poorly constrained τ [49] leads to δ ranging from microns to nanometers. We demonstrate that nano-scale texture causes a shift in the nucleation to the three-phase contact line, while micro-scale texture does not. The possibility of a critical length scale has implications for the effectiveness of nucleation catalysts, including formation of ice in atmospheric clouds [7].

Velocity recovery cycles of human muscle action potentials and their sensitivity to ischemia

This study was undertaken to test whether recovery cycle measurements can provide useful information about the membrane potential of human muscle fibers. Multifiber responses to direct muscle stimulation through needle electrodes were recorded from the brachioradialis of healthy volunteers, and the latency changes measured as conditioning stimuli were applied at interstimulus intervals of 2-1000 ms. In all subjects, the relative refractory period (RRP), which lasted 3.27 +/- 0.45 ms (mean +/- SD, n = 12), was followed by a phase of supernormality, in which the velocity increased by 9.3 +/- 3.4% at 6.1 +/- 1.3 ms, and recovered over 1 s. A broad hump of additional supernormality was seen at around 100 ms. Extra conditioning stimuli had little effect on the early supernormality but increased the later component. The two phases of supernormality resembled early and late afterpotentials, attributable respectively to the passive decay of membrane charge and potassium accumulation in the t-tubules. Five minutes of ischemia progressively prolonged the RRP and reduced supernormality, confirming that these parameters are sensitive to membrane depolarization. Velocity recovery cycles may provide useful information about altered muscle membrane potential and t-tubule function in muscle disease. Muscle Nerve, 2008.

Two-Phase Flow In Microchannels: Morphology And Interface Phenomena

The existence and morphology, as well as the dynamics of micro-scale gas-liquid interfaces is investigated numerically and experimentally. These studies can be used to assess liquid management issues in microsystems such as PEMFC gas flow channels, and are meant to open new research perspectives in two-phase flow, particularly in film deposition on non-wetting surfaces. For example the critical plug volume data can be used to deliver desired length plugs, or to determine the plug formation frequency. The dynamics of gas-liquid interfaces, of interest for applications involving small passages (e.g. heat exchangers, phase separators and filtration systems), was investigated using high-speed microscopy - a method that also proved useful for the study of film deposition processes. The existence limit for a liquid plug forming in a mixed wetting channel is determined by numerical simulations using Surface Evolver. The plug model simulate actual conditions in the gas flow channels of PEM fuel cells, the wetting of the gas diffusion layer (GDL) side of the channel being different from the wetting of the bipolar plate walls. The minimum plug volume, denoted as critical volume is computed for a series of GDL and bipolar plate wetting properties. Critical volume data is meant to assist in the water management of PEMFC, when corroborated with experimental data. The effect of cross section geometry is assessed by computing the critical volume in square and trapezoidal channels. Droplet simulations show that water can be passively removed from the GDL surface towards the bipolar plate if we take advantage on differing wetting properties between the two surfaces, to possibly avoid the gas transport blockage through the GDL. High speed microscopy was employed in two-phase and film deposition experiments with water in round and square capillary tubes. Periodic interface destabilization was observed and the existence of compression waves in the gas phase is discussed by taking into consideration a naturally occurring convergent-divergent nozzle formed by the flowing liquid phase. The effect of channel geometry and wetting properties was investigated through two-phase water-air flow in square and round microchannels, having three static contact angles of 20, 80 and 105 degrees. Four different flow regimes are observed for a fixed flow rate, this being thought to be caused by the wetting behavior of liquid flowing in the corners as well as the liquid film stability. Film deposition experiments in wetting and non-wetting round microchannels show that a thicker film is deposited for wetting conditions departing from the ideal 0 degrees contact angle. A film thickness dependence with the contact angle theta as well as the Capillary number, in the form h_R ~ Ca^(2/3)/ cos(theta) is inferred from scaling arguments, for contact angles smaller than 36 degrees. Non-wetting film deposition experiments reveal that a film significantly thicker than the wetting Bretherton film is deposited. A hydraulic jump occurs if critical conditions are met, as given by a proposed nondimensional parameter similar to the Froude number. Film thickness correlations are also found by matching the measured and the proposed velocity derived in the shock theory. The surface wetting as well as the presence of the shock cause morphological changes in the Taylor bubble flow.

Parametric estimation of pulse arrival time: a robust approach to pulse wave velocity

Pulse wave velocity (PWV) is a surrogate of arterial stiffness and represents a non-invasive marker of cardiovascular risk. The non-invasive measurement of PWV requires tracking the arrival time of pressure pulses recorded in vivo, commonly referred to as pulse arrival time (PAT). In the state of the art, PAT is estimated by identifying a characteristic point of the pressure pulse waveform. This paper demonstrates that for ambulatory scenarios, where signal-to-noise ratios are below 10 dB, the performance in terms of repeatability of PAT measurements through characteristic points identification degrades drastically. Hence, we introduce a novel family of PAT estimators based on the parametric modeling of the anacrotic phase of a pressure pulse. In particular, we propose a parametric PAT estimator (TANH) that depicts high correlation with the Complior(R) characteristic point D1 (CC = 0.99), increases noise robustness and reduces by a five-fold factor the number of heartbeats required to obtain reliable PAT measurements.

The role of the sinuses of Valsalva in aortic root flow dynamics and aortic root surgery: evaluation by magnetic resonance imaging

BACKGROUND AND AIM OF THE STUDY: Combined replacement of the aortic valve and ascending aorta using a composite graft represents the standard treatment for dilated aortic root with concomitant structural damage of the aortic valve, especially when the aortic valve cannot be preserved. Unfortunately, hemodynamic changes associated with prosthetic replacement of the aortic root have not been fully elucidated. The study aim was to compare hemodynamics within the replaced aortic root using either a prosthetic vascular graft with bulges mimicking the sinuses of Valsalva and including a stented pericardial valve, or a straight xenopericardial conduit and a stentless porcine valve. METHODS: Between July 2004 and March 2006, a total of 35 patients (mean age 65.2 years: range: 32-80 years) was enrolled into the present study. Aortic root replacement was performed in nine patients with a Valsalva graft (Gelweave Valsalva; Vascutek, Renfrewshire, UK) including a stented pericardial valve, and in 19 patients with a xenopericardial conduit containing a stentless porcine valve. All patients underwent postoperative magnetic resonance imaging (MRI). A control group of seven patients allowed for comparison with native aortic root hemodynamics. RESULTS: Maximum flow-velocity above the aortic valve as one marker of compliance of the aortic root was slightly higher in patients with a Valsalva graft compared to native aortic roots (1.9 m/s versus 1.3 m/s, p = 0.001), but was significantly lower than in patients with the xenopericardial graft without neo-sinuses (1.3 m/s versus 2.4 m/s, p < 0.001). CONCLUSION: The pre-shaped bulges in the prosthetic Valsalva graft effectively mimic the native sinuses of Valsalva, improve compliance of the aortic root, and result in a more physiologic flow pattern, as demonstrated by postoperative MRI.

Assessment of splanchnic blood flow using magnetic resonance imaging

BACKGROUND AND AIMS: The splanchnic circulation has an important function in the body under both physiological and pathophysiological conditions. Despite its importance, no reliable noninvasive procedures for estimating splanchnic circulation have been established. The aim of this study was to evaluate MRI as a tool for assessing intra-abdominal blood flows of the aorta, portal vein (VPO) and the major intestinal and hepatic vessels. METHODS: In nine healthy volunteers, the proximal aorta (AOP) and distal abdominal aorta (AOD), superior mesenteric artery (SAM), celiac trunk (CTR), hepatic arteries (common and proper hepatic arteries, AHC and AHP, respectively), and VPO were localized on contrast-enhanced magnetic resonance angiography images. Volumetric flow was measured using a two-dimensional cine echocardiogram-gated phase contrast technique. Measurements were taken before and 30 min after continuous intravenous infusion of somatostatin (250 microg/h) and were independently evaluated by two investigators. RESULTS: Blood flow measured by MRI in the VPO, SAM, AOP, AHP, and CTR significantly decreased after drug infusion. Flows in the AOD and AHC showed a tendency to decrease (P>0.05). Interrater agreement on flows in MRI was very good for large vessels (VPO, AOP, and AOD), with a concordance correlation coefficient of 0.94, as well as for smaller vessels such as the CTR, AHC, AHP, and SAM (concordance correlation coefficient =0.78). CONCLUSION: Somatostatin-induced blood flow changes in the splanchnic region were reliably detected by MRI. MRI may be useful for the noninvasive assessment of blood flow changes in the splanchnic region.

3D Ultrasound imaging--a useful non-invasive tool to detect AV fistulas in transplanted kidneys

BACKGROUND: A precise, non-invasive, non-toxic, repeatable, convenient and inexpensive follow-up of renal transplants, especially following biopsies, is in the interest of nephrologists. Formerly, the rate of biopsies leading to AV fistulas had been underestimated. Imaging procedures suited to a detailed judgement of these vascular malformations are to be assessed. METHODS: Three-dimensional (3D) reconstruction techniques of ultrasound flow-directed and non-flow-directed energy mode pictures were compared with a standard procedure, gadolinium-enhanced nuclear magnetic resonance imaging angiography (MRA) using the phase contrast technique. RESULTS: Using B-mode and conventional duplex information, AV fistulas were localized in the upper pole of the kidney transplant of the index patient. The 3D reconstruction provided information about the exact localization and orientation of the fistula in relation to other vascular structures, and the flow along the fistula. The MRA provided localization and orientation information, but less functional information. Flow-directed and non-flow-directed energy mode pictures could be reconstructed to provide 3D information about vascular malformations in transplanted kidneys. CONCLUSION: In transplanted kidneys, 3D-ultrasound angiography may be equally as effective as MRA in localizing and identifying AV malformations. Advantages of the ultrasound method are that it is cheaper, non-toxic, non-invasive, more widely availability and that it even provides more functional information. Future prospective studies will be necessary to evaluate the two techniques further.

Esophagus imaging for catheter ablation of atrial fibrillation: comparison of two methods with showing of esophageal movement

This study aims to evaluate whether visualization and integration of the computed tomography (CT) scan of the left atrium (LA) and the esophagus into the three-dimensional (3D) electroanatomical map the day before ablation is accurate compared with integration of an esophagus tag into the electroanatomic LA map visualizing the anatomic relationship during the radiofrequency ablation or whether esophagus movement prohibits esophagus visualization the day before ablation.

Phase diagram with an enhanced spin-glass region of the mixed Ising–XY magnet LiHo{x}Er{1-x}F{4}

We present the experimental phase diagram of LiHoxEr1-xF4, a dilution series of dipolar-coupled model magnets. The phase diagram was determined using a combination of ac susceptibility and neutron scattering. Three unique phases in addition to the Ising ferromagnet LiHoF4 and the XY antiferromagnet LiErF4 have been identified. Below x = 0.86, an embedded spin-glass phase is observed, where a spin glass exists within the ferromagnetic structure. Below x = 0.57, an Ising spin glass is observed consisting of frozen needlelike clusters. For x ∼ 0.3–0.1, an antiferromagnetically coupled spin glass occurs. A reduction of TC(x) for the ferromagnet is observed which disobeys the mean-field predictions that worked for LiHoxY1-xF4.

Inverse relationship between fractionated electrograms and atrial fibrosis in persistent atrial fibrillation: combined magnetic resonance imaging and high-density mapping

OBJECTIVES This study sought to evaluate the relationship between fibrosis imaged by delayed-enhancement (DE) magnetic resonance imaging (MRI) and atrial electrograms (Egms) in persistent atrial fibrillation (AF). BACKGROUND Atrial fractionated Egms are strongly related to slow anisotropic conduction. Their relationship to atrial fibrosis has not yet been investigated. METHODS Atrial high-resolution MRI of 18 patients with persistent AF (11 long-lasting persistent AF) was registered with mapping geometry (NavX electro-anatomical system (version 8.0, St. Jude Medical, St. Paul, Minnesota)). DE areas were categorized as dense or patchy, depending on their DE content. Left atrial Egms during AF were acquired using a high-density, 20-pole catheter (514 ± 77 sites/map). Fractionation, organization/regularity, local mean cycle length (CL), and voltage were analyzed with regard to DE. RESULTS Patients with long-lasting persistent versus persistent AF had larger left atrial (LA) surface area (134 ± 38 cm(2) vs. 98 ± 9 cm(2), p = 0.02), a higher amount of atrial DE (70 ± 16 cm(2) vs. 49 ± 10 cm(2), p = 0.01), more complex fractionated atrial Egm (CFAE) extent (54 ± 16 cm(2) vs. 28 ± 15 cm(2), p = 0.02), and a shorter baseline AF CL (147 ± 10 ms vs. 182 ± 14 ms, p = 0.01). Continuous CFAE (CFEmean [NavX algorithm that quantifies Egm fractionation] <80 ms) occupied 38 ± 19% of total LA surface area. Dense DE was detected at the left posterior left atrium. In contrast, the right posterior left atrium contained predominantly patchy DE. Most CFAE (48 ± 14%) occurred at non-DE LA sites, followed by 41 ± 12% CFAE at patchy DE and 11 ± 6% at dense DE regions (p = 0.005 and p = 0.008, respectively); 19 ± 6% CFAE sites occurred at border zones of dense DE. Egms were less fractionated, with longer CL and lower voltage at dense DE versus non-DE regions: CFEmean: 97 ms versus 76 ms, p < 0.0001; local CL: 153 ms versus 143 ms, p < 0.0001; mean voltage: 0.63 mV versus 0.86 mV, p < 0.0001. CONCLUSIONS Atrial fibrosis as defined by DE MRI is associated with slower and more organized electrical activity but with lower voltage than healthy atrial areas. Ninety percent of continuous CFAE sites occur at non-DE and patchy DE LA sites. These findings are important when choosing the ablation strategy in persistent AF.

Integration of merged delayed-enhanced magnetic resonance imaging and multidetector computed tomography for the guidance of ventricular tachycardia ablation: a pilot study

BACKGROUND Delayed enhancement (DE) MRI can assess the fibrotic substrate of scar-related VT. MDCT has the advantage of inframillimetric spatial resolution and better 3D reconstructions. We sought to evaluate the feasibility and usefulness of integrating merged MDCT/MRI data in 3D-mapping systems for structure-function assessment and multimodal guidance of VT mapping and ablation. METHODS Nine patients, including 3 ischemic cardiomyopathy (ICM), 3 nonischemic cardiomyopathy (NICM), 2 myocarditis, and 1 redo procedure for idiopathic VT, underwent MRI and MDCT before VT ablation. Merged MRI/MDCT data were integrated in 3D-mapping systems and registered to high-density endocardial and epicardial maps. Low-voltage areas (<1.5 mV) and local abnormal ventricular activities (LAVA) during sinus rhythm were correlated to DE at MRI, and wall-thinning (WT) at MDCT. RESULTS Endocardium and epicardium were mapped with 391 ± 388 and 1098 ± 734 points per map, respectively. Registration of MDCT allowed visualization of coronary arteries during epicardial mapping/ablation. In the idiopathic patient, integration of MRI data identified previously ablated regions. In ICM patients, both DE at MRI and WT at MDCT matched areas of low voltage (overlap 94 ± 6% and 79 ± 5%, respectively). In NICM patients, wall-thinning areas matched areas of low voltage (overlap 63 ± 21%). In patients with myocarditis, subepicardial DE matched areas of epicardial low voltage (overlap 92 ± 12%). A total number of 266 LAVA sites were found in 7/9 patients. All LAVA sites were associated to structural substrate at imaging (90% inside, 100% within 18 mm). CONCLUSION The integration of merged MDCT and DEMRI data is feasible and allows combining substrate assessment with high-spatial resolution to better define structure-function relationship in scar-related VT.

Human intervertebral disc stiffness correlates better with the Otsu threshold computed from axial T2 map of its posterior annulus fibrosus than with clinical classifications

Degeneration of the intervertebral disc, sometimes associated with low back pain and abnormal spinal motions, represents a major health issue with high costs. A non-invasive degeneration assessment via qualitative or quantitative MRI (magnetic resonance imaging) is possible, yet, no relation between mechanical properties and T2 maps of the intervertebral disc (IVD) has been considered, albeit T2 relaxation time values quantify the degree of degeneration. Therefore, MRI scans and mechanical tests were performed on 14 human lumbar intervertebral segments freed from posterior elements and all soft tissues excluding the IVD. Degeneration was evaluated in each specimen using morphological criteria, qualitative T2 weighted images and quantitative axial T2 map data and stiffness was calculated from the load-deflection curves of in vitro compression, torsion, lateral bending and flexion/extension tests. In addition to mean T2, the OTSU threshold of T2 (TOTSU), a robust and automatic histogram-based method that computes the optimal threshold maximizing the distinction of two classes of values, was calculated for anterior, posterior, left and right regions of each annulus fibrosus (AF). While mean T2 and degeneration schemes were not related to the IVDs' mechanical properties, TOTSU computed in the posterior AF correlated significantly with those classifications as well as with all stiffness values. TOTSU should therefore be included in future degeneration grading schemes.