Radiologists’ gaze-related parameters combined with image-based features were utilized to classify suspicious mammographic areas ultimately scored as True Positives (TP) and False Positives (FP). Eight breast radiologists read 120 two-view digital mammograms of which 59 had biopsy proven cancer. Eye tracking data was collected and nearby fixations were clustered together. Suspicious areas on mammograms were independently identified based on thresholding an intensity saliency map followed by automatic segmentation and pruning steps. For each radiologist reported area, radiologist’s fixation clusters in the area, as well as neighboring suspicious areas within 2.5° of the center of fixation, were found. A 45-dimensional feature vector containing gaze parameters of the corresponding cluster along with image-based characteristics was constructed. Gaze parameters included total number of fixations in the cluster, dwell time, time to hit the cluster for the first time, maximum number of consecutive fixations, and saccade magnitude of the first fixation in the cluster. Image-based features consisted of intensity, shape, and texture descriptors extracted from the region around the suspicious area, its surrounding tissue, and the entire breast. For each radiologist, a userspecific Support Vector Machine (SVM) model was built to classify the reported areas as TPs or FPs. Leave-one-out cross validation was utilized to avoid over-fitting. A feature selection step was embedded in the SVM training procedure by allowing radial basis function kernels to have 45 scaling factors. The proposed method was compared with the radiologists’ performance using the jackknife alternative free-response receiver operating characteristic (JAFROC). The JAFROC figure of merit increased significantly for six radiologists.
Abstract Although many semi-automated and automated algorithms for breast density assessment have been recently proposed, none of these have been widely accepted. In this study a novel automated algorithm, named iDensity, inspired by the human visual system is proposed for classifying mammograms into four breast density categories corresponding to the Breast Imaging Reporting and Data System (BI-RADS). For each BI-RADS category 80 cases were taken from the normal volumes of the Digital Database for Screening Mammography (DDSM). For each case only the left medio-lateral oblique was utilized. After image calibration using the provided tables of each scanner in the DDSM, the pectoral muscle and background were removed. Images were filtered by a median filter and down sampled. Images were then filtered by a filter bank consisting of Gabor filters in six orientations and 3 scales, as well as a Gaussian filter. Three gray level histogram-based features and three second order statistics features were extracted from each filtered image. Using the extracted features, mammograms were separated initially separated into two groups, low or high density, then in a second stage, the low density group was subdivided into BI-RADS I or II, and the high density group into BI-RADS III or IV. The algorithm achieved a sensitivity of 95% and specificity of 94% in the first stage, sensitivity of 89% and specificity of 95% when classifying BIRADS I and II cases, and a sensitivity of 88% and 91% specificity when classifying BI-RADS III and IV.
Attenuation correction in positron emission tomography brain imaging of freely moving animals can be very challenging since the body of the animal is often within the field of view and introduces a non negligible atten- uating factor that can degrade the quantitative accuracy of the reconstructed images. An attractive approach that avoids the need for a transmission scan involves the generation of the convex hull of the animal’s head based on the reconstructed emission images. However, this approach ignores the potential attenuation introduced by the animal’s body. In this work, we propose a virtual scanner geometry, which moves in synchrony with the animal’s head and discriminates between those events that traverse only the animal’s head (and therefore can be accurately compensated for attenuation) and those that might have also traversed the animal’s body. For each pose a new virtual scanner geometry was defined and therefore a new system matrix was calculated leading to a time-varying system matrix. This new approach was evaluated on phantom data acquired on the microPET Focus 220 scanner using a custom-made rat phantom. Results showed that when the animal’s body is within the FOV and not accounted for during attenuation correction it can lead to bias of up to 10%. On the contrary, at- tenuation correction was more accurate when the virtual scanner was employed leading to improved quantitative estimates (bias <2%), without the need to account for the animal’s body.
List-mode image reconstruction with motion correction is computationally expensive, as it requires projection of hundreds of millions of rays through a 3D array. To decrease reconstruction time it is possible to use symmetric multiprocessing computers or graphics processing units. The former can have high financial costs, while the latter can require refactoring of algorithms. The Xeon Phi is a new co-processor card with a Many Integrated Core architecture that can run 4 multiple-instruction, multiple data threads per core with each thread having a 512-bit single instruction, multiple data vector register. Thus, it is possible to run in the region of 220 threads simultaneously. The aim of this study was to investigate whether the Xeon Phi co-processor card is a viable alternative to an x86 Linux server for accelerating List-mode PET image reconstruction for motion correction. An existing list-mode image reconstruction algorithm with motion correction was ported to run on the Xeon Phi coprocessor with the multi-threading implemented using pthreads. There were no differences between images reconstructed using the Phi co-processor card and images reconstructed using the same algorithm run on a Linux server. However, it was found that the reconstruction runtimes were 3 times greater for the Phi than the server. A new version of the image reconstruction algorithm was developed in C++ using OpenMP for mutli-threading and the Phi runtimes decreased to 1.67 times that of the host Linux server. Data transfer from the host to co-processor card was found to be a rate-limiting step; this needs to be carefully considered in order to maximize runtime speeds. When considering the purchase price of a Linux workstation with Xeon Phi co-processor card and top of the range Linux server, the former is a cost-effective computation resource for list-mode image reconstruction. A multi-Phi workstation could be a viable alternative to cluster computers at a lower cost for medical imaging applications.
Positron emission mammography ("PEM") is a breast imaging modality that typically involves the administration of
relatively high doses of radiotracer. In order to reduce tracer costs and consider PEM for global screening applications, it
would be helpful to reduce the required amount of administered radiotracer so that patient dose would be comparable to
conventional x-ray mammograms.
We performed GATE Monte Carlo investigations of several possible camera configurations. Increasing the detector
thickness from 10 to 30 mm, increasing the camera surface area from 5×20cm2 to 20×20cm2, and applying depth-ofinteraction information to increase the acceptance angle, increased the overall efficiency to radiation emitted from a breast cancer by a factor of 24 as compared to existing commercial systems.
KEYWORDS: Sensors, Collimators, Single photon emission computed tomography, Reconstruction algorithms, Expectation maximization algorithms, Image processing, Monte Carlo methods, Heart, Signal to noise ratio, Spatial resolution
Single photon emission computerized tomographic (SPECT) images often suffer from low resolution and low count
density. To improve spatial resolution of SPECT it is possible to use a pinhole collimator; however, this further reduces
the system sensitivity. A potential solution to this problem is to use coded apertures, which offers increased sensitivity
by using multiple pinholes, at the cost of increased image reconstruction time.
A generic reconstruction algorithm has been developed which allows for arbitrary acquisition geometry via affine
transforms (translation and rotation). The reconstruction process uses a (Siddon) ray projector, the expectation
maximization (EM) algorithm and a 1 to n pinhole position matrix. Iteration times scale as a function of the number of
pinholes in the collimator. Resolution recovery has also been incorporated into the reconstruction algorithm.
The algorithm developed allows for the investigation of optimal imaging settings for small animal imaging. Simulated
acquisitions of an ex-vivo rat heart with 1, 5 and 8 pinholes, over 360 degree acquisition, showing that multi-pinhole
imaging can be successfully applied to rat cardiac imaging. Further refinement of the acquisition parameters, such as
image overlap, collimator pinhole configuration and geometrical imaging configuration, will predict the theoretical
settings for quantitative cardiac multi-pinhole SPECT imaging.
Yi-Hwa Liu, Zakir Sahul, Christopher Weyman, William Ryder, Donald Dione, Lawrence Dobrucki, Choukri Mekkaoui, Matthew Brennan, Xiaoyue Hu, Christi Hawley, Albert Sinusas
We have developed a new single photon emission computerized tomography (SPECT) hotspot quantification
method incorporating extra cardiac activity correction and hotspot normal limit estimation. The method was validated
for estimation accuracy of myocardial tracer focal uptake in a chronic canine model of myocardial infarction (MI). Dogs
(n = 4) at 2 weeks post MI were injected with Tl-201 and a
Tc-99m-labeled hotspot tracer targeted at matrix
metalloproteinases (MMPs). An external point source filled with
Tc-99m was used for a reference of absolute
radioactivity. Dual-isotope (Tc-99m/Tl-201) SPECT images were acquired simultaneously followed by an X-ray CT
acquisition. Dogs were sacrificed after imaging for myocardial gamma well counting. Images were reconstructed with
CT-based attenuation correction (AC) and without AC (NAC) and were quantified using our quantification method.
Normal limits for myocardial hotspot uptake were estimated based on 3 different schemes: maximum entropy,
meansquared-error minimization (MSEM) and global minimization. Absolute myocardial hotspot uptake was quantified from
SPECT images using the normal limits and compared with well-counted radioactivity on a segment-by-segment basis (n = 12 segments/dog). Radioactivity was expressed as % injected dose (%ID). There was an excellent correlation (r = 0.78-0.92) between the estimated activity (%ID) derived using the SPECT quantitative approach and
well-counting, independent of AC. However, SPECT quantification without AC resulted in the significant underestimation of radioactivity. Quantification using SPECT with AC and the MSEM normal limit yielded the best results compared with well-counting. In conclusion, focal myocardial "hotspot" uptake of a targeted radiotracer can be accurately quantified in
vivo using a method that incorporates SPECT imaging with AC, an external reference, background scatter compensation,
and a suitable normal limit. This hybrid SPECT/CT approach allows for the serial non-invasive quantitative evaluation
of molecular targeted tracers in the heart.
Access to the requested content is limited to institutions that have purchased or subscribe to SPIE eBooks.
You are receiving this notice because your organization may not have SPIE eBooks access.*
*Shibboleth/Open Athens users─please
sign in
to access your institution's subscriptions.
To obtain this item, you may purchase the complete book in print or electronic format on
SPIE.org.
INSTITUTIONAL Select your institution to access the SPIE Digital Library.
PERSONAL Sign in with your SPIE account to access your personal subscriptions or to use specific features such as save to my library, sign up for alerts, save searches, etc.