Numerical methods for solar tomography in the STEREO era

In this thesis, we propose several advances in the numerical and computational algorithms that are used to determine tomographic estimates of physical parameters in the solar corona. We focus on methods for both global dynamic estimation of the coronal electron density and estimation of local transient phenomena, such as coronal mass ejections, from empirical observations acquired by instruments onboard the STEREO spacecraft. We present a first look at tomographic reconstructions of the solar corona from multiple points-of-view, which motivates the developments in this thesis. In particular, we propose a method for linear equality constrained state estimation that leads toward more physical global dynamic solar tomography estimates. We also present a formulation of the local static estimation problem, i.e., the tomographic estimation of local events and structures like coronal mass ejections, that couples the tomographic imaging problem to a phase field based level set method. This formulation will render feasible the 3D tomography of coronal mass ejections from limited observations. Finally, we develop a scalable algorithm for ray tracing dense meshes, which allows efficient computation of many of the tomographic projection matrices needed for the applications in this thesis.

Dynamic Imaging in Electrical Impedance Tomography of the Human Chest With Online Transition Matrix Identification

One of the electrical impedance tomography objectives is to estimate the electrical resistivity distribution in a domain based only on electrical potential measurements at its boundary generated by an imposed electrical current distribution into the boundary. One of the methods used in dynamic estimation is the Kalman filter. In biomedical applications, the random walk model is frequently used as evolution model and, under this conditions, poor tracking ability of the extended Kalman filter (EKF) is achieved. An analytically developed evolution model is not feasible at this moment. The paper investigates the identification of the evolution model in parallel to the EKF and updating the evolution model with certain periodicity. The evolution model transition matrix is identified using the history of the estimated resistivity distribution obtained by a sensitivity matrix based algorithm and a Newton-Raphson algorithm. To numerically identify the linear evolution model, the Ibrahim time-domain method is used. The investigation is performed by numerical simulations of a domain with time-varying resistivity and by experimental data collected from the boundary of a human chest during normal breathing. The obtained dynamic resistivity values lie within the expected values for the tissues of a human chest. The EKF results suggest that the tracking ability is significantly improved with this approach.

Capillaries within compartments: Microvascular interpretation of dynamic positron emission tomography data

Measurement of exchange of substances between blood and tissue has been a long-lasting challenge to physiologists, and considerable theoretical and experimental accomplishments were achieved before the development of the positron emission tomography (PET). Today, when modeling data from modern PET scanners, little use is made of earlier microvascular research in the compartmental models, which have become the standard model by which the vast majority of dynamic PET data are analysed. However, modern PET scanners provide data with a sufficient temporal resolution and good counting statistics to allow estimation of parameters in models with more physiological realism. We explore the standard compartmental model and find that incorporation of blood flow leads to paradoxes, such as kinetic rate constants being time-dependent, and tracers being cleared from a capillary faster than they can be supplied by blood flow. The inability of the standard model to incorporate blood flow consequently raises a need for models that include more physiology, and we develop microvascular models which remove the inconsistencies. The microvascular models can be regarded as a revision of the input function. Whereas the standard model uses the organ inlet concentration as the concentration throughout the vascular compartment, we consider models that make use of spatial averaging of the concentrations in the capillary volume, which is what the PET scanner actually registers. The microvascular models are developed for both single- and multi-capillary systems and include effects of non-exchanging vessels. They are suitable for analysing dynamic PET data from any capillary bed using either intravascular or diffusible tracers, in terms of physiological parameters which include regional blood flow. (C) 2003 Elsevier Ltd. All rights reserved.

Differentiation of Metabolically Distinct Areas within Head and Neck Region using Dynamic 18F-FDG Positron Emission Tomography Imaging

Positron Emission Tomography (PET) using ¹⁸F-FDG is playing a vital role in the diagnosis and treatment planning of cancer. However, the most widely used radiotracer, ¹⁸F-FDG, is not specific for tumours and can also accumulate in inflammatory lesions as well as normal physiologically active tissues making diagnosis and treatment planning complicated for the physicians. Malignant, inflammatory and normal tissues are known to have different pathways for glucose metabolism which could possibly be evident from different characteristics of the time activity curves from a dynamic PET acquisition protocol. Therefore, we aimed to develop new image analysis methods, for PET scans of the head and neck region, which could differentiate between inflammation, tumour and normal tissues using this functional information within these radiotracer uptake areas. We developed different dynamic features from the time activity curves of voxels in these areas and compared them with the widely used static parameter, SUV, using Gaussian Mixture Model algorithm as well as K-means algorithm in order to assess their effectiveness in discriminating metabolically different areas. Moreover, we also correlated dynamic features with other clinical metrics obtained independently of PET imaging. The results show that some of the developed features can prove to be useful in differentiating tumour tissues from inflammatory regions and some dynamic features also provide positive correlations with clinical metrics. If these proposed methods are further explored then they can prove to be useful in reducing false positive tumour detections and developing real world applications for tumour diagnosis and contouring.

Cryo-electron tomography of Kaposi's sarcoma-associated herpesvirus capsids reveals dynamic scaffolding structures essential to capsid assembly and maturation.

Kaposi's sarcoma-associated herpesvirus (KSHV) is a recently discovered DNA tumor virus that belongs to the gamma-herpesvirus subfamily. Though numerous studies on KSHV and other herpesviruses, in general, have revealed much about their multilayered organization and capsid structure, the herpesvirus capsid assembly and maturation pathway remains poorly understood. Structural variability or irregularity of the capsid internal scaffolding core and the lack of adequate tools to study such structures have presented major hurdles to earlier investigations employing more traditional cryo-electron microscopy (cryoEM) single particle reconstruction. In this study, we used cryo-electron tomography (cryoET) to obtain 3D reconstructions of individual KSHV capsids, allowing direct visualization of the capsid internal structures and systematic comparison of the scaffolding cores for the first time. We show that B-capsids are not a structurally homogenous group; rather, they represent an ensemble of "B-capsid-like" particles whose inner scaffolding is highly variable, possibly representing different intermediates existing during the KSHV capsid assembly and maturation. This information, taken together with previous observations, has allowed us to propose a detailed pathway of herpesvirus capsid assembly and maturation.

Evaluation of liver parenchyma and perfusion using dynamic contrast-enhanced computed tomography and contrast-enhanced ultrasonography in captive green iguanas (Iguana iguana) under general anesthesia.

BACKGROUND Contrast-enhanced diagnostic imaging techniques are considered useful in veterinary and human medicine to evaluate liver perfusion and focal hepatic lesions. Although hepatic diseases are a common occurrence in reptile medicine, there is no reference to the use of contrast-enhanced ultrasound (CEUS) and contrast-enhanced computed tomography (CECT) to evaluate the liver in lizards. Therefore, the aim of this study was to evaluate the pattern of change in echogenicity and attenuation of the liver in green iguanas (Iguana iguana) after administration of specific contrast media. RESULTS An increase in liver echogenicity and density was evident during CEUS and CECT, respectively. In CEUS, the mean ± SD (median; range) peak enhancement was 19.9% ± 7.5 (18.3; 11.7-34.6). Time to peak enhancement was 134.0 ± 125.1 (68.4; 59.6-364.5) seconds. During CECT, first visualization of the contrast medium was at 3.6 ± 0.5 (4; 3-4) seconds in the aorta, 10.7 ± 2.2 (10.5; 7-14) seconds in the hepatic arteries, and 15 ± 4.5 (14.5; 10-24) seconds in the liver parenchyma. Time to peak was 14.1 ± 3.4 (13; 11-21) and 31 ± 9.6 (29; 23-45) seconds in the aorta and the liver parenchyma, respectively. CONCLUSION CEUS and dynamic CECT are practical means to determine liver hemodynamics in green iguanas. Distribution of contrast medium in iguana differed from mammals. Specific reference ranges of hepatic perfusion for diagnostic evaluation of the liver in iguanas are necessary since the use of mammalian references may lead the clinician to formulate incorrect diagnostic suspicions.

Towards clinical computed ultrasound tomography in echo-mode: Dynamic range artefact reduction

Computed ultrasound tomography in echo-mode (CUTE) allows imaging the speed of sound inside tissue using hand-held pulse-echo ultrasound. This technique is based on measuring the changing local phase of beamformed echoes when changing the transmit beam steering angle. Phantom results have shown a spatial resolution and contrast that could qualify CUTE as a promising novel diagnostic modality in combination with B-mode ultrasound. Unfortunately, the large intensity range of several tens of dB that is encountered in clinical images poses difficulties to echo phase tracking and results in severe artefacts. In this paper we propose a modification to the original technique by which more robust echo tracking can be achieved, and we demonstrate in phantom experiments that dynamic range artefacts are largely eliminated. Dynamic range artefact reduction also allowed for the first time a clinical implementation of CUTE with sufficient contrast to reproducibly distinguish the different speed of sound in different tissue layers of the abdominal wall and the neck.

When static meets dynamic: Comparing cone-beam computed tomography and acoustic reflection for upper airway analysis

INTRODUCTION: Upper airway measurement can be important for the diagnosis of breathing disorders. Acoustic reflection (AR) is an accepted tool for studying the airway. Our objective was to investigate the differences between cone-beam computed tomography (CBCT) and AR in calculating airway volumes and areas. METHODS: Subjects with prescribed CBCT images as part of their records were also asked to have AR performed. A total of 59 subjects (mean age, 15 ± 3.8 years) had their upper airway (5 areas) measured from CBCT images, acoustic rhinometry, and acoustic pharyngometry. Volumes and minimal cross-sectional areas were extracted and compared with software. RESULTS: Intraclass correlation on 20 randomly selected subjects, remeasured 2 weeks apart, showed high reliability (r >0.77). Means of total nasal volume were significantly different between the 2 methods (P = 0.035), but anterior nasal volume and minimal cross-sectional area showed no differences (P = 0.532 and P = 0.066, respectively). Pharyngeal volume showed significant differences (P = 0.01) with high correlation (r = 0.755), whereas pharyngeal minimal cross-sectional area showed no differences (P = 0.109). The pharyngeal volume difference may not be considered clinically significant, since it is 758 mm3 for measurements showing means of 11,000 ± 4000 mm3. CONCLUSIONS: CBCT is an accurate method for measuring anterior nasal volume, nasal minimal cross-sectional area, pharyngeal volume, and pharyngeal minimal cross-sectional area.

Mechanical ventilation for the treatment of severe excessive dynamic airway collapse

[Excerpt] We read with interest the case report by Ismael et al1 describing a patient with Sjo¨gren’s syndrome and cystic lung disease who could not be weaned from a ventilator due to severe central excessive dynamic airway collapse (EDAC) of the lower part of the trachea and proximal bronchi. EDAC corresponds to the expiratory bulging of the tracheobronchial wall without known airway structural abnormalities, leading to a decrease of at least 50% in internal diameter.2 It is a rare and underdiagnosed entity, commonly confused with other respiratory diseases such as asthma and COPD. Although noninvasive procedures such as cervicothoracic computed tomography scan on inspiration and expiration may suggest the disorder, the accepted standard method for diagnosis is bronchoscopy.3-7 (...).

Octopus Standard Automated Perimetry, Pulsar Perimetry, Moorfields Motion Displacement Test and Heidelberg Retinal Tomography in Glaucoma Detection - A Comparison of Diagnostic Accuracy

Purpose: To investigate the accuracy of 4 clinical instruments in the detection of glaucomatous damage. Methods: 102 eyes of 55 test subjects (Age mean = 66.5yrs, range = [39; 89]) underwent Heidelberg Retinal Tomography (HRTIII), (disc area<2.43); and standard automated perimetry (SAP) using Octopus (Dynamic); Pulsar (TOP); and Moorfields Motion Displacement Test (MDT) (ESTA strategy). Eyes were separated into three groups 1) Healthy (H): IOP<21mmHg and healthy discs (clinical examination), 39 subjects, 78 eyes; 2) Glaucoma suspect (GS): Suspicious discs (clinical examination), 12 subjects, 15 eyes; 3) Glaucoma (G): progressive structural or functional loss, 14 subjects, 20 eyes. Clinical diagnostic precision was examined using the cut-off associated with the p<5% normative limit of MD (Octopus/Pulsar), PTD (MDT) and MRA (HRT) analysis. The sensitivity, specificity and accuracy were calculated for each instrument. Results: See table Conclusions: Despite the advantage of defining glaucoma suspects using clinical optic disc examination, the HRT did not yield significantly higher accuracy than functional measures. HRT, MDT and Octopus SAP yielded higher accuracy than Pulsar perimetry, although results did not reach statistical significance. Further studies are required to investigate the structure-function correlations between these instruments.

Increase of [(18)F]FLT tumor uptake in vivo mediated by FdUrd: toward improving cell proliferation positron emission tomography.

PURPOSE: 3'-deoxy-3'-[(18)F]fluorothymidine ([(18)F]FLT), a cell proliferation positron emission tomography (PET) tracer, has been shown in numerous tumors to be more specific than 2-deoxy-2-[(18)F]fluoro-D-glucose ([(18)F]FDG) but less sensitive. We studied the capacity of a nontoxic concentration of 5-fluoro-2'-deoxyuridine (FdUrd), a thymidine synthesis inhibitor, to increase uptake of [(18)F]FLT in tumor xenografts. METHODS: The duration of the FdUrd effect in vivo on tumor cell cycling and thymidine analogue uptake was studied by varying FdUrd pretreatment timing and holding constant the timing of subsequent flow cytometry and 5-[(125)I]iodo-2'-deoxyuridine biodistribution measurements. In [(18)F]FLT studies, FdUrd pretreatment was generally performed 1 h before radiotracer injection. [(18)F]FLT biodistributions were measured 1 to 3 h after radiotracer injection of mice grafted with five different human tumors and pretreated or not with FdUrd and compared with [(18)F]FDG tumor uptake. Using microPET, the dynamic distribution of [(18)F]FLT was followed for 1.5 h in FdUrd pretreated mice. High-field T2-weighted magnetic resonance imaging (MRI) and histology were used comparatively in assessing tumor viability and proliferation. RESULTS: FdUrd induced an immediate increase in tumor uptake of 5-[(125)I]iodo-2'-deoxyuridine, that vanished after 6 h, as also confirmed by flow cytometry. Biodistribution measurements showed that FdUrd pretreatment increased [(18)F]FLT uptake in all tumors by factors of 3.2 to 7.8 compared with controls, while [(18)F]FDG tumor uptake was about fourfold and sixfold lower in breast cancers and lymphoma. Dynamic PET in FdUrd pretreated mice showed that [(18)F]FLT uptake in all tumors increased steadily up to 1.5 h. MRI showed a well-vascularized homogenous lymphoma with high [(18)F]FLT uptake, while in breast cancer, a central necrosis shown by MRI was inactive in PET, consistent with the histomorphological analysis. CONCLUSION: We showed a reliable and significant uptake increase of [(18)F]FLT in different tumor xenografts after low-dose FdUrd pretreatment. These results show promise for a clinical application of FdUrd aimed at increasing the sensitivity of [(18)F]FLT PET.

Three-dimensional ordered-subset expectation maximization iterative protocol for evaluation of left ventricular volumes and function by quantitative gated SPECT: a dynamic phantom study.

The purposes of this study were to characterize the performance of a 3-dimensional (3D) ordered-subset expectation maximization (OSEM) algorithm in the quantification of left ventricular (LV) function with (99m)Tc-labeled agent gated SPECT (G-SPECT), the QGS program, and a beating-heart phantom and to optimize the reconstruction parameters for clinical applications. METHODS: A G-SPECT image of a dynamic heart phantom simulating the beating left ventricle was acquired. The exact volumes of the phantom were known and were as follows: end-diastolic volume (EDV) of 112 mL, end-systolic volume (ESV) of 37 mL, and stroke volume (SV) of 75 mL; these volumes produced an LV ejection fraction (LVEF) of 67%. Tomographic reconstructions were obtained after 10-20 iterations (I) with 4, 8, and 16 subsets (S) at full width at half maximum (FWHM) gaussian postprocessing filter cutoff values of 8-15 mm. The QGS program was used for quantitative measurements. RESULTS: Measured values ranged from 72 to 92 mL for EDV, from 18 to 32 mL for ESV, and from 54 to 63 mL for SV, and the calculated LVEF ranged from 65% to 76%. Overall, the combination of 10 I, 8 S, and a cutoff filter value of 10 mm produced the most accurate results. The plot of the measures with respect to the expectation maximization-equivalent iterations (I x S product) revealed a bell-shaped curve for the LV volumes and a reverse distribution for the LVEF, with the best results in the intermediate range. In particular, FWHM cutoff values exceeding 10 mm affected the estimation of the LV volumes. CONCLUSION: The QGS program is able to correctly calculate the LVEF when used in association with an optimized 3D OSEM algorithm (8 S, 10 I, and FWHM of 10 mm) but underestimates the LV volumes. However, various combinations of technical parameters, including a limited range of I and S (80-160 expectation maximization-equivalent iterations) and low cutoff values (< or =10 mm) for the gaussian postprocessing filter, produced results with similar accuracies and without clinically relevant differences in the LV volumes and the estimated LVEF.

Increase of [(18)F]FLT Tumor Uptake In Vivo Mediated by FdUrd: Toward Improving Cell Proliferation Positron Emission Tomography.

PURPOSE: 3'-deoxy-3'-[(18)F]fluorothymidine ([(18)F]FLT), a cell proliferation positron emission tomography (PET) tracer, has been shown in numerous tumors to be more specific than 2-deoxy-2-[(18)F]fluoro-D: -glucose ([(18)F]FDG) but less sensitive. We studied the capacity of a nontoxic concentration of 5-fluoro-2'-deoxyuridine (FdUrd), a thymidine synthesis inhibitor, to increase uptake of [(18)F]FLT in tumor xenografts. METHODS: The duration of the FdUrd effect in vivo on tumor cell cycling and thymidine analogue uptake was studied by varying FdUrd pretreatment timing and holding constant the timing of subsequent flow cytometry and 5-[(125)I]iodo-2'-deoxyuridine biodistribution measurements. In [(18)F]FLT studies, FdUrd pretreatment was generally performed 1 h before radiotracer injection. [(18)F]FLT biodistributions were measured 1 to 3 h after radiotracer injection of mice grafted with five different human tumors and pretreated or not with FdUrd and compared with [(18)F]FDG tumor uptake. Using microPET, the dynamic distribution of [(18)F]FLT was followed for 1.5 h in FdUrd pretreated mice. High-field T2-weighted magnetic resonance imaging (MRI) and histology were used comparatively in assessing tumor viability and proliferation. RESULTS: FdUrd induced an immediate increase in tumor uptake of 5-[(125)I]iodo-2'-deoxyuridine, that vanished after 6 h, as also confirmed by flow cytometry. Biodistribution measurements showed that FdUrd pretreatment increased [(18)F]FLT uptake in all tumors by factors of 3.2 to 7.8 compared with controls, while [(18)F]FDG tumor uptake was about fourfold and sixfold lower in breast cancers and lymphoma. Dynamic PET in FdUrd pretreated mice showed that [(18)F]FLT uptake in all tumors increased steadily up to 1.5 h. MRI showed a well-vascularized homogenous lymphoma with high [(18)F]FLT uptake, while in breast cancer, a central necrosis shown by MRI was inactive in PET, consistent with the histomorphological analysis. CONCLUSION: We showed a reliable and significant uptake increase of [(18)F]FLT in different tumor xenografts after low-dose FdUrd pretreatment. These results show promise for a clinical application of FdUrd aimed at increasing the sensitivity of [(18)F]FLT PET.

Dynamic model of the left ventricle for use in simulation of myocardial perfusion SPECT and gated SPECT

Simulation is a useful tool in cardiac SPECT to assess quantification algorithms. However, simple equation-based models are limited in their ability to simulate realistic heart motion and perfusion. We present a numerical dynamic model of the left ventricle, which allows us to simulate normal and anomalous cardiac cycles, as well as perfusion defects. Bicubic splines were fitted to a number of control points to represent endocardial and epicardial surfaces of the left ventricle. A transformation from each point on the surface to a template of activity was made to represent the myocardial perfusion. Geometry-based and patient-based simulations were performed to illustrate this model. Geometry-based simulations modeled ~1! a normal patient, ~2! a well-perfused patient with abnormal regional function, ~3! an ischaemic patient with abnormal regional function, and ~4! a patient study including tracer kinetics. Patient-based simulation consisted of a left ventricle including a realistic shape and motion obtained from a magnetic resonance study. We conclude that this model has the potential to study the influence of several physical parameters and the left ventricle contraction in myocardial perfusion SPECT and gated-SPECT studies.