Multislice tomographic imaging and analysis of human breast-equivalent phantoms and biological tissues

A laser transillumination tomographic system, consisting of electrical, optical, mechanical, and software components, to obtain multislice images of tissue-equivalent breast phantoms and biological tissues, is developed. The tissue-equivalent phantoms are prepared from paraffin wax mixed with wax color pigments by matching their surface backscattered profiles as measured by multiprobe laser reflectometer, with that of respective tissues. The optical parameters of these phantoms are determined by matching their reflectance profiles with that as obtained by Monte Carlo simulation of optical scattering. For multislice tomographic analysis conical breast phantoms of height 80.0 mm and 80.0 mm base diameter with inclusions of different optical properties and dimensions are developed. The resolution of the inclusions in the tomograms depends on their sizes and optical parameters. The minimum size of the inclusion as detected by this procedure in a slice of diameter 50.0 mm is 3.0 mm. The structural variation as observed in the tomograms of phantoms of combination of biological tissues indicates its possible applications in detecting the abnormalities developing in human healthy soft tissues.     

A viewpoint determination system for stenosis diagnosis andquantification in coronary angiographic image acquisition

This paper describes the usefulness of computer assistance in the acquisition of “good” images for stenosis diagnosis and quantification in coronary angiography. The system recommends the optimal viewpoints from which stenotic lesions can be observed clearly based on images obtained from initial viewpoints. First, the viewpoint dependency of the apparent severity of a stenotic lesion is experimentally analyzed using software phantoms in order to show the seriousness of the problem. The implementation of the viewpoint determination system is then described. The system provides good user-interactive tools for the semi-automated estimation of the orientation and diameter of stenotic segments and the three-dimensional (3-D) reconstruction of vessel structures. Using these tools, viewpoints that will not give rise to foreshortening and vessel overlap can be efficiently determined. Experiments using real coronary angiograms show the system to be capable of the reliable diagnosis and quantification of stenosis     

激光照射血液疗法中光在血液中的传输特性

本文采用和Monte Carlo方法研究了激光血管内照射,在不同的激光入射角、不同的激光发散角和不同的血管直径、以及入射点偏离血管中心轴线情况下,光在血管液中的分布或吸收状况。研究说明,激光血管内照射应确保入射光纤与被照射血管相平行,从光纤输出的激光要有较小的发散角,被照射血管的直径应大于6mm,使激光能量得以充分吸收,并可减小入射点偏离中心轴线对光能量吸收的影响。     

Multiple Sensor Optical Thermometry System for Application in Clinical Hyperthermia

The thermometry system described is based upon the temperature dependence of the band edge absorption of infrared light in GaAs crystal. The design of the thermometry was completed, and the system was subjected to an extensive evaluation, including testing with tissue phantoms and microwave applicators. The system has up to 12 temperature sensors which are packaged in two basic probe coinfigurations: a single-sensor probe with a length of 1.2 m and a diameter of 0.6 mm; and a four-sensor linear array probe with a length of 1.2 m, diameter of 1.1 mm, and spacing of 1.5 cm between adjacent sensors. Results of thermometry evaluation are presented, including data on automatic calibration, temperature accuracy and stability, and EMI protection.     

Three-dimensional finite-element analyses for radio-frequencyhepatic tumor ablation

Radio-frequency (RF) hepatic ablation, offers an alternative method for the treatment of hepatic malignancies. We employed finite-element method (FEM) analysis to determine tissue temperature distribution during RF hepatic ablation. We constructed three-dimensional (3-D) thermal-electrical FEM models consisting of a four-tine RF probe, hepatic tissue, and a large blood vessel (10-mm diameter) located at different locations. We simulated our FEM analyses under temperature-controlled (90 degrees C) 8-min ablation. We also present a preliminary result from a simplified two-dimensional (2-D) FEM model that includes a bifurcated blood vessel. Lesion shapes created by the four-tine RF probe were mushroom-like, and were limited by the blood vessel. When the distance of the blood vessel was 5 mm from the nearest distal electrode 1) in the 3-D model, the maximum tissue temperature (hot spot) appeared next to electrodes A. The location of the hot spot was adjacent to another electrode 2) on the opposite side when the blood vessel was 1 mm from electrode A. The temperature distribution in the 2-D model was highly nonuniform due to the presence of the bifurcated blood vessel. Underdosed areas might be present next to the blood vessel from which the tumor can regenerate.     

Reconstruction and quantification of the carotid artery bifurcation from 3-D ultrasound images

Three-dimensional (3-D) ultrasound is a relatively new technique, which is well suited to imaging superficial blood vessels, and potentially provides a useful, noninvasive method for generating anatomically realistic 3-D models of the peripheral vasculature. Such models are essential for accurate simulation of blood flow using computational fluid dynamics (CFD), but may also be used to quantify atherosclerotic plaque more comprehensively than routine clinical methods. In this paper, we present a spline-based method for reconstructing the normal and diseased carotid artery bifurcation from images acquired using a freehand 3-D ultrasound system. The vessel wall (intima-media interface) and lumen surfaces are represented by a geometric model defined using smoothing splines. Using this coupled wall-lumen model, we demonstrate how plaque may be analyzed automatically to provide a comprehensive set of quantitative measures of size and shape, including established clinical measures, such as degree of (diameter) stenosis. The geometric accuracy of 3-D ultrasound reconstruction is assessed using pulsatile phantoms of the carotid bifurcation, and we conclude by demonstrating the in vivo application of the algorithms outlined to 3-D ultrasound scans from a series of patient carotid arteries.     

Confocal microwave imaging for breast cancer detection: delay-multiply-and-sum image reconstruction algorithm

A new image reconstruction algorithm, termed as delay-multiply-and-sum (DMAS), for breast cancer detection using an ultra-wideband confocal microwave imaging technique is proposed. In DMAS algorithm, the backscattered signals received from numerical breast phantoms simulated using the finite-difference time-domain method are time shifted, multiplied in pair, and the products are summed to form a synthetic focal point. The effectiveness of the DMAS algorithm is shown by applying it to backscattered signals received from a variety of numerical breast phantoms. The reconstructed images illustrate improvement in identification of embedded malignant tumors over the delay-and-sum algorithm. Successful detection and localization of tumors as small as 2 mm in diameter are also demonstrated.     

Radar-Based Breast Cancer Detection Using a Hemispherical Antenna Array—Experimental Results

In this contribution, an ultrawideband (UWB) microwave system for breast cancer detection is presented. The system is based on a novel hemispherical real-aperture antenna array, which is employed in a multi-static radar-based detection system. The array consists of 16 UWB aperture-coupled stacked-patch antennas located on a section of a hemisphere. The radar system is designed to be used with realistic three-dimensional (3D) breast phantoms, which have been developed, as well as with real breast cancer patients during initial clinical trials. Images are formed using two different beamforming algorithms and the performance of these algorithms is firstly compared through numerical simulation. Experimental results for the same beamforming techniques are then presented, demonstrating the successful detection of 4 and 6 mm diameter spherical tumors in the curved breast phantom.      

Medical imaging with a microwave tomographic scanner

A microwave tomographic scanner for biomedical applications is presented. It consists of a 64-element circular array with a useful diameter of 20 cm. Electronically scanning the transmitting and receiving antennas allows multiview measurements with no mechanical movement. Imaging parameters-a spatial resolution of 7 mm and a contrast resolution of 1% for a measurement time of 3 s-are appropriate for medical use. Measurements on tissue-simulating phantoms and volunteers, together with numerical simulations, are presented to assess the system for absolute imaging of tissue distribution and for differential imaging of physiological, pathological, and induced changes in tissues     

Blood Velocity Measurement by Fiber Optic Laser Doppler Anemometry

Laser Doppler anemometry has been used to measure blood velocities by using a fine fiber optic probe both to deliver laser light into the blood and to receive reflected light from red blood cells at the tip of the probe. In vitro, with the probe aligned in the direction of flow, nonlinear calibration is possible to velocities of 15 cm/s in a 3 mm diameter vessel. When the probe is aligned against the flow, measurement of velocity is linear to at least 1 m/s. In vivo recordings of coronary and aortic blood velocity are produced by performing a fast Fourier transform (FET) in real time.     

Adaptive approach to accurate analysis of small-diameter vessels incineangiograms

In coronary vessels smaller than 1 mm in diameter, it is difficult to accurately identify lumen borders using existing border detection techniques. Computer-detected diameters of small coronary vessels are often severely overestimated due to the influence of the imaging system point spread function and the use of an edge operator designed for a broad range of diameter vessel sizes. Computer-detected diameters may be corrected if a calibration curve for the X-ray system is available. Unfortunately, the performance of this postprocessing diameter correction approach is severely limited by the presence of image noise. The authors report here a new approach that uses a two-stage adaption of edge operator parameters to optimally match the edge operator to the local lumen diameter. In the first stage, approximate lumen diameters are detected using a single edge operator in a half-resolution image. Depending on the approximate lumen size, one of three edge operators is selected for the second full-resolution stage in which left and right coronary borders are simultaneously identified. The method was tested in a set of 72 segments of nine angiographic phantom vessels with diameters ranging from 0.46 to 4.14 mm and in 82 clinical coronary angiograms. Performance of the adaptive simultaneous border detection method was compared to that of a conventional border detection method and to that of a postprocessing diameter correction border detection method. Adaptive border detection yielded significantly improved accuracy in small phantom vessels and across all vessel sizes in comparison to the conventional and postprocessing diameter correction methods (p<0.001 in all cases). Adaptive simultaneous coronary border detection provides both accurate and robust quantitative analysis of coronary vessels of all sizes     

Matched filter estimation of serial blood vessel diameters from video images

A method for making a contiguous series of blood vessel diameter estimates from digitized images is proposed. It makes use of a vessel intensity profile model based on the vessel geometry and the physics of the imaging process, providing estimates of far greater accuracy than previously obtained. A variety of techniques are used to reduce the computational demand. The method includes the generation of measurement estimation error, which is important in determining total vessel patency as well as providing a basic measure of diameter estimate accuracy.     

Three-dimensional trajectory assessment of an IVUS transducer from single-plane cineangiograms: a phantom study

The recovery of the three-dimensional (3-D) path of the transducer used during an intravascular ultrasound (IVUS) examination is of primary importance to assess the exact 3-D shape of the vessel under study. Traditionally, the reconstruction is done by simply stacking the images during the pullback, or more recently using biplane angiography to recover the vessel curvature. In this paper, we explain, how single-plane angiography can be used with two projection models, to perform this task. Two types of projection geometry are analyzed: weak-perspective and full-perspective. In weak-perspective projection geometry, the catheter path can be reconstructed without prior transducer depth information. With full-perspective projection geometry, precise depth location of reference points are needed in order to minimize the error of the recovered transducer angle of incidence. The transducer angulation reconstruction is based on the foreshortening effect as seen from the X-ray images. By comparing the measured to the true transducer length, we are able to get its incidence angle. The transducer trajectory is reconstructed by stitching together the different estimated angulations obtained from each image in a cineangiogram sequence. The method is described and validated on two helical vessel phantoms, giving on average a reconstructed path that is less than 2 mm distant from the true path when using full-perspective projection.     

Quantitative coronary angiography with deformable spline models

Although current edge-following schemes can be very efficient in determining coronary boundaries, they may fail when the feature to be followed is disconnected (and the scheme is unable to bridge the discontinuity) or branch points exist where the best path to follow is indeterminate. Here, the authors present new deformable spline algorithms for determining vessel boundaries, and enhancing their centerline features. A bank of even and odd S-Gabor filter pairs of different orientations are convolved with vascular images in order to create an external snake energy field. Each fitter pair will give maximum response to the segment of vessel having the same orientation as the filters. The resulting responses across filters of different orientations are combined to create an external energy field for snake optimization. Vessels are represented by B-Spline snakes, and are optimized on filter outputs with dynamic programming. The points of minimal constriction and the percent-diameter stenosis are determined from a computed vessel centerline. The system has been statistically validated using fixed stenosis and flexible-tube phantoms. It has also been validated on 20 coronary lesions with two independent operators, and has been tested for interoperator and intraoperator variability and reproducibility. The system has been found to be specially robust in complex images involving vessel branchings and incomplete contrast filling     

A mean blood velocity statistic for the Doppler signal from anarrow ultrasound beam

An easily calculable statistic proportional to the instantaneous spatial mean blood velocity through a vessel cross section is derived from Doppler power spectral estimates for the case where the Doppler beam is assumed to be of negligible thickness compared to the vessel diameter. This is an alternative statistic to that derived where uniform insonation is assumed, an assumption thought to be poorer in many real cases. The main requirement is that the velocity profile is monotonic increasing from the vessel wall to the vessel center. Errors in each statistic are compared for a variety of true beam dimensions, using a variety of velocity profiles, and the new statistic is shown to incur less error for Gaussian beam response profiles with a standard deviation less than 0.4 of the vessel radius, or for rectangular response profiles with a width less than 0.65 of the vessel diameter. If an estimate can be made of the true beam dimensions and vessel diameter, a weighted sum of the two statistics can give a more accurate estimate of mean velocity. The case of a beam displaced from the center of the vessel is also considered, and errors are found to be less than 4% for a displacement of 20% of the vessel radius     

3-D segmentation algorithm of small lung nodules in spiral CT images.

Computed tomography (CT) is the most sensitive imaging technique for detecting lung nodules, and is now being evaluated as a screening tool for lung cancer in several large samples studies all over the world. In this report, we describe a semiautomatic method for 3-D segmentation of lung nodules in CT images for subsequent volume assessment. The distinguishing features of our algorithm are the following. 1) The user interaction process. It allows the introduction of the knowledge of the expert in a simple and reproducible manner. 2) The adoption of the geodesic distance in a multithreshold image representation. It allows the definition of a fusion--segregation process based on both gray-level similarity and objects shape. The algorithm was validated on low-dose CT scans of small nodule phantoms (mean diameter 5.3--11 mm) and in vivo lung nodules (mean diameter 5--9.8 mm) detected in the Italung-CT screening program for lung cancer. A further test on small lung nodules of Lung Image Database Consortium (LIDC) first data set was also performed. We observed a RMS error less than 6.6% in phantoms, and the correct outlining of the nodule contour was obtained in 82/95 lung nodules of Italung-CT and in 10/12 lung nodules of LIDC first data set. The achieved results support the use of the proposed algorithm for volume measurements of lung nodules examined with low-dose CT scanning technique.     

A new model-based technique for enhanced small-vessel measurements in X-ray ciné-angiograms

Arterial diameter estimation from X-ray ciné angiograms is important for quantifying coronary artery disease (CAD) and for evaluating therapy. However, diameter measurement in vessel cross sections < or =1.0 mm is associated with large measurement errors. We present a novel diameter estimator which reduces both magnitude and variability of measurement error. We use a parametric nonlinear imaging model for X-ray ciné angiography and estimate unknown model parameters directly from the image data. Our technique allows us to exploit additional diameter information contained within the intensity profile amplitude, a feature which is overlooked by existing methods. This method uses a two-step procedure: the first step estimates the imaging model parameters directly from the angiographic frame and the second step uses these measurements to estimate the diameter of vessels in the same image. In Monte-Carlo simulation over a range of imaging conditions, our approach consistently produced lower estimation error and variability than conventional methods. With actual X-ray images, our estimator is also better than existing methods for the diameters examined (0.4-4.0 mm). These improvements are most significant in the range of narrow vessel widths associated with severe coronary artery disease.     

A system for real-time measurement of the brachial artery diameter in B-mode ultrasound images.

The measurement of the brachial artery diameter is frequently used in clinical studies for evaluating the flow-mediated dilation and, in conjunction with the blood pressure value, for assessing arterial stiffness. This paper presents a system for computing the brachial artery diameter in real-time by analyzing B-mode ultrasound images. The method is based on a robust edge detection algorithm which is used to automatically locate the two walls of the vessel. The measure of the diameter is obtained with subpixel precision and with a temporal resolution of 25 samples/s, so that the small dilations induced by the cardiac cycle can also be retrieved. The algorithm is implemented on a standalone video processing board which acquires the analog video signal from the ultrasound equipment. Results are shown in real-time on a graphical user interface. The system was tested both on synthetic ultrasound images and in clinical studies of flow-mediated dilation. Accuracy, robustness, and intra/inter observer variability of the method were evaluated.     

A System for Real-Time Measurement of the Brachial Artery Diameter in B-Mode Ultrasound Images

The measurement of the brachial artery diameter is frequently used in clinical studies for evaluating the flow-mediated dilation and, in conjunction with the blood pressure value, for assessing arterial stiffness. This paper presents a system for computing the brachial artery diameter in real-time by analyzing B-mode ultrasound images. The method is based on a robust edge detection algorithm which is used to automatically locate the two walls of the vessel. The measure of the diameter is obtained with subpixel precision and with a temporal resolution of 25 samples/s, so that the small dilations induced by the cardiac cycle can also be retrieved. The algorithm is implemented on a standalone video processing board which acquires the analog video signal from the ultrasound equipment. Results are shown in real-time on a graphical user interface. The system was tested both on synthetic ultrasound images and in clinical studies of flow-mediated dilation. Accuracy, robustness, and intra/inter observer variability of the method were evaluated     

Three-dimensional reconstruction of transillumination tomographic images of human breast phantoms by red and infrared lasers

A laser transillumination tomographic system, using red and near infrared lasers, to obtain cross-sectional images of human breast phantoms and human hand is proposed. The scanning assembly is consisting of upward, downward, rotational and pitch movements. The phantoms, made of paraffin wax, agar gel, and milk placed in a glass model, are embedded with abnormalities like blood, water, solid objects, and tissues. By illuminating the phantoms at various heights by either red or infrared laser the projection data are collected. Based on 64 projections the tomogram of each section is constructed. By volume visualization procedure applied to the tomograms the objects of varying composition embedded within the phantoms are detected and their size, shape, and location depth are determined. The cross-sectional image of a human hand obtained by this procedure further shows the possibility of application of this technique for imaging of organs.