Haim Azhari

Super-resolution in PET imaging

This paper demonstrates a super-resolution method for improving the resolution in clinical positron emission tomography (PET) scanners. Super-resolution images were obtained by combining four data sets with spatial shifts between consecutive acquisitions and applying an iterative algorithm. Super-resolution attenuation corrected PET scans of a phantom were obtained using the two-dimensional and three-dimensional (3-D) acquisition modes of a clinical PET/computed tomography (CT) scanner (Discovery LS, GEMS). In a patient study, following a standard /sup 18/F-FDG PET/CT scan, a super-resolution scan around one small lesion was performed using axial shifts without increasing the patient radiation exposure. In the phantom study, smaller features (3 mm) could be resolved axially with the super-resolution method than without (6 mm). The super-resolution images had better resolution than the original images and provided higher contrast ratios in coronal images and in 3-D acquisition transaxial images. The coronal super-resolution images had superior resolution and contrast ratios compared to images reconstructed by merely interleaving the data to the proper axial location. In the patient study, super-resolution reconstructions displayed a more localized /sup 18/F-FDG uptake. A new approach for improving the resolution of PET images using a super-resolution method has been developed and experimentally confirmed, employing a clinical scanner. The improvement in axial resolution requires no changes in hardware.

Correction of Heart Motion Due to Respiration in Clinical Myocardial Perfusion SPECT Scans Using Respiratory Gating

Several studies have described nonuniform blurring of myocardial perfusion imaging (MPI) due to respiration. This article describes a technique for correcting the respiration effect and assesses its effectiveness in clinical studies. Methods: Simulated phantoms, physical phantoms, and patient scans were used in this study. A heart phantom, which oscillated back and forth, was used to simulate respiration. The motion was measured on a γ-camera supporting list-mode functionality synchronized with an external respiratory strap or resistor sensor. Eight clinical scans were performed using a 1-d 99mTc-sestamibi protocol while recording the respiratory signal. The list-mode capability along with the strap or sensor signals was used to generate respiratory bin projection sets. A segmentation process was used to detect the shift between the respiratory bins. This shift was further projected to the acquired projection images for correction of the respiratory motion. The process was applied to the phantom and patient studies, and the rate of success of the correction was assessed using the conventional bulls eye maps. Results: The algorithm provided a good correction for the phantom studies. The shift after the correction, measured by a fitted ellipsoid, was <1 mm in the axial direction. The average motion due to respiration in the clinical studies was 9.1 mm in the axial direction. The average shift between the respiratory phases was reduced to 0.5 mm after correction. The maximal change in the bulls eye map for the clinical scans after the correction was 6%, with a mean of 3.75%. The postcorrection clinical summed perfusion images were more uniform, consistent, and, for some patients, clinically significant when compared with the images before correction for respiration. Conclusion: Myocardial motion generated by respiration during MPI SPECT affects perfusion image quality and accuracy. Motion introduced by respiration can be corrected using the proposed method. The degree of correction depends on the patient respiratory pattern and can be of clinical significance in certain cases.

Reconstruction in diffraction ultrasound tomography using nonuniform FFT

We show an iterative reconstruction framework for diffraction ultrasound tomography. The use of broad-band illumination allows significant reduction of the number of projections compared to straight ray tomography. The proposed algorithm makes use of forward nonuniform fast Fourier transform (NUFFT) for iterative Fourier inversion. Incorporation of total variation regularization allows the reduction of noise and Gibbs phenomena while preserving the edges. The complexity of the NUFFT-based reconstruction is comparable to the frequency-domain interpolation (gridding) algorithm, whereas the reconstruction accuracy (in sense of the L/sup 2/ and the L/sup /spl infin// norm) is better.

Three-dimensional automatic quantitative analysis of intravascular ultrasound images

Intravascular ultrasound (IVUS) has established itself as a useful tool for coronary assessment. The vast amount of data obtained by a single IVUS study renders manual analysis impractical for clinical use. A computerized method is needed to accelerate the process and eliminate user-dependency. In this study, a new algorithm is used to identify the lumen border and the media-adventitia border (the external elastic membrane). Setting an initial surface on the IVUS catheter perimeter and using active contour principles, the surface inflates until virtual force equilibrium defined by the surface geometry and image features is reached. The method extracts these features in three dimensions (3-D). Eight IVUS procedures were performed using an automatic pullback device. Using the ECG signal for synchronization, sets of images covering the entire studied region and corresponding to the same cardiac phase were sampled. Lumen and media-adventitia border contours were traced manually and compared to the automatic results obtained by the suggested method. Linear regression results for vessel area enclosed by the lumen and media-adventitia border indicate high correlation between manual vs. automatic tracings (y = 1.07 x -0.38; r = 0.98; SD = 0.112 mm(2); n = 88). These results indicate that the suggested algorithm may potentially provide a clinical tool for accurate lumen and plaque assessment.

Effects of afterload on regional left ventricular torsion

OBJECTIVE To determine if left ventricular torsion, as measured by magnetic resonance tissue tagging, is afterload dependent in a canine isolated heart model in which neurohumoral responses are absent, and preload is constant. METHODS In ten isolated, blood perfused, ejecting, canine hearts, three afterloads were studied, while keeping preload constant: low afterload, high afterload (stroke volume reduced by approx. 50% of low afterload), and isovolumic loading (infinite afterload). RESULTS There were significant effects of afterload on both torsion (P < 0.05) and circumferential shortening (P < 0.0005). Between low and high afterloads, at the anterior region of the endocardium only, where torsion was maximal, there was a significant reduction in torsion (15.1 +/- 2.2 degrees to 7.8 +/- 1.8 degrees, P < 0.05). Between high afterload and isovolumic loading there was no significant change in torsion (7.8 +/- 1.8 degrees to 6.2 +/- 1.5 degrees, P = NS). Circumferential shortening at the anterior endocardium was significantly reduced both between low and high afterload (-0.19 +/- 0.02 to -0.11 +/- 0.02, P < 0.0005), and also between high afterload and isovolumic loading (-0.11 +/- 0.02 to 0.00 +/- 0.02, P < 0.05). Plots of strains with respect to end-systolic volume demonstrated a reduction in both torsion and shortening with afterload-induced increases in end-systolic volume. Torsion, but not circumferential shortening, persisted at isovolumic loading. CONCLUSIONS Maximal regional torsion of the left ventricle is afterload dependent. The afterload response of torsion appears related to the effects of afterload on end-systolic volume.

Automated Detection and Characterization of Multiple Sclerosis Lesions in Brain MR Images

In the present study an automatic algorithm for detection and contouring of multiple sclerosis (MS) lesions in brain magnetic resonance (MR) images is introduced. This algorithm automatically detects MS lesions in axial proton density, T2-weighted, gadolinium enhanced, and fast fluid attenuated inversion recovery (FLAIR) brain MR images. Automated detection consists of three main stages: (1) detection and contouring of all hyperintense signal regions within the image; (2) partial elimination of false positive segments (defined herein as artifacts) by size, shape index, and anatomical location; (3) the use of an artificial neural paradigm (Back-Propagation) for final removal of artifacts by differentiating them from true MS lesions. The algorithm was applied to 45 images acquired from 14 MS patients. The algorithms sensitivity was 0.87 and the specificity 0.96. In 34 images, 100% of the lesions were detected. The algorithm potentially may serve as a useful preprocessing tool for quantitative MS monitoring via magnetic resonance imaging.

A conical model to describe the nonuniformity of the left ventricular twisting motion

The systolic contraction and fiber shortening in the left ventricle (LV) produces torsional moments in the myocardium, resulting in a gradient of angular displacements about the long axis. This is manifested as a counterclockwise rotation of the apex relative to the base, when viewed from the apex. Recent studies with magnetic resonance imaging (MRI), using noninvasive magnetic tags, have revealed three important properties of the LV twist: (a) The angle of twist (i.e., the angular rotation of a slice relative to the basal slice) is consistently higher at the endocardium as compared to the epicardium; (b) The twist increases towards the apex; and (c) Straight MRI-tagged radial lines at end-diastole (ED) are slightly curved at end-systole (ES), implying a nonlinear transmural variation of the twist. The present study suggests that the geometry of the LV at ES can be represented by a thick-walled hollow cone, and that the transmural twist patterns from ED to ES can be described using the continuum mechanics approach and a small strain analysis of an isotropic cone subjected to external torque. The predicted results are compared with the noninvasive MRI measurements of transmural twist in eight human volunteers. Given the epicardial angles of twist of each slice, the predicted endocardial angles of twist are in good correlation with the experimental findings (r=0.86, slope=1.09, SEE=4.1°). In addition, the model reliably describes the changes in the twist magnitude from apex to base (no significant difference from experimental values,P=0.2), and predicts the curvilinear pattern at ES of the originally straight ED radial lines. Thus, the conical model with uniform properties of the LV, reliably predicts the nonuniformity of the twist patterns, implying that the LV twist is strongly affected by LV geometry.

Regional three-dimensional geometry and function of left ventricles with fibrous aneurysms. A cine-computed tomography study.

BackgroundTo assess the extent and nature of the dysfunction surrounding aneurysms of the left ventricle (LV), we examined the parameters of local and global three-dimensional shape, size, and function of LVs of eight patients with histologically confirmed anterior fibrous aneurysms. Methods and ResultsThree-dimensional reconstructions of each LV were made from 10-12 short-axis fast cine-angiographic computed tomography (cine-CT) slices encompassing the entire heart at end diastole and end systole. Regional three-dimensional wall thickness, thickening, motion, curvature, and stress index were calculated for 84 elements encompassing the entire LV. The aneurysmal border was defined by a sharp decrease in end-diastolic wall thickness and separated the LV into an aneurysmal zone and a normal zone that was further divided into adjacent normal (AN) and remote normal (RN) zones. As expected, thickening was negligible in both the aneurysmal and the border zones. Although both the AN and the RN zones had normal wall thickness (1.05 ± 0.20 and 1.09 ± 0.20 cm, respectively), thickening was depressed in the AN (0.22 ± 0.08 cm) but not the RN (0.44 + 0.19 cm) zones. The size of the dysfunction zone (defined as less than 2 mm thickening) was found to be considerably greater than the anatomic size of the aneurysm (60.9 ± 13.7% versus 33.6 + 7.6% of the left ventricular endocardial area, respectively; p < 0.00l). In addition, the AN zone had a smaller curvature and a higher stress index than the RN zone. ConclusionsLVs with fibrous aneurysms are characterized by a relatively large region of nonfunction that encompasses the thin aneurysmal area and its transitional border zone, a normally functioning remote zone, and an intermediate region of normal wall thickness but with reduced function, which may be attributed to its low curvature and high stress index.

A Method for Characterization of Tissue Elastic Properties Combining Ultrasonic Computed Tomography With Elastography

Objective. The correlation between various diseases and the change in the local mechanical properties of soft tissues has been long known. Over the past 20 years, there have been increasing research efforts to characterize mechanical properties of biological tissues using ultrasonic elastography. However, most of these works were based on characterization of only 1 type of waves (longitudinal or shear). The goal of this work was to devise a comprehensive ultrasound‐based imaging method capable of measuring elastic parameters by combining both backscattered elastography and through‐transmitted ultrasonic computed tomography. Methods. Our suggested technique provides measurements of both longitudinal and shear wave velocities. This enables the noninvasive computation of several tissue elasticity parameters such as Youngs and shear moduli, Poissons ratio, and, more importantly, the bulk modulus, the determination of which requires both wave velocities. Four different phantom types were examined: agar‐gelatin–based phantoms and porcine fat tissue, turkey breast tissue, and bovine liver tissue in vitro specimens. The values of Youngs modulus, the shear modulus, and Poissons ratio were estimated and were consistent with values published in the literature. Results. The average bulk modulus values of the phantoms ± SD were 2.83 ± 0.001, 2.25 ± 0.01, 2.48 ± 0.01, and 2.53 ± 0.02 GPa, respectively. A statistically significant difference (P < .001) in the values of the bulk modulus of the different phantoms was found. Conclusions. The bulk modulus is suitable for differentiation between different tissue types. The obtained results show the feasibility of using a comprehensive ultrasonic imaging technique for noninvasive quantitative tissue characterization.

The reduction of artifacts due to metal hip implants in CT-attenuation corrected PET images from hybrid PET/CT scanners

CT beam hardening artifacts near metal hip implants may erroneously enhance or diminish radiotracer uptake following CT attenuation correction (AC) of PET images. An artifact reduction algorithm (ARA) was developed to reduce metal artifacts in CT-based AC-PET. The algorithm employed a Bayes classifier to identify beam-hardening artifacts, followed by a partial correction of the attenuation map. ARA was implemented on phantom and patient 18F-FDG studies using a clinical PET/CT scanner. In phantom studies ARA successfully removed two artifacts of erroneously elevated uptake near a stainless steel hip prosthesis which were depicted in the standard CT-AC PET. ARA has also identified two targets absent on the scanner PET images. Target-to-background ratios were 1.5–3 times higher for ARA-PET than scanner images. In a patient study, metal artifacts were of lower intensity in ARA-PET as compared to standard images. Potentially, ARA may improve detectability of small lesions located near metal hip implants.