Image quality assessment in digital mammography: part II. NPWE as a validated alternative for contrast detail analysis

Assessment of image quality for digital x-ray mammography systems used in European screening programs relies mainly on contrast-detail CDMAM phantom scoring and requires the acquisition and analysis of many images in order to reduce variability in threshold detectability. Part II of this study proposes an alternative method based on the detectability index (d') calculated for a non-prewhitened model observer with an eye filter (NPWE). The detectability index was calculated from the normalized noise power spectrum and image contrast, both measured from an image of a 5 cm poly(methyl methacrylate) phantom containing a 0.2 mm thick aluminium square, and the pre-sampling modulation transfer function. This was performed as a function of air kerma at the detector for 11 different digital mammography systems. These calculated d' values were compared against threshold gold thickness (T) results measured with the CDMAM test object and against derived theoretical relationships. A simple relationship was found between T and d', as a function of detector air kerma; a linear relationship was found between d' and contrast-to-noise ratio. The values of threshold thickness used to specify acceptable performance in the European Guidelines for 0.10 and 0.25 mm diameter discs were equivalent to threshold calculated detectability indices of 1.05 and 6.30, respectively. The NPWE method is a validated alternative to CDMAM scoring for use in the image quality specification, quality control and optimization of digital x-ray systems for screening mammography.     

Amorphous selenium flat panel detectors for digital mammography: validation of a NPWE model observer with CDMAM observer performance experiments

Model observers have been developed which incorporate a specific imaging task, system performance, and human observer characteristics and can potentially overcome some of the limitations in using detective quantum efficiency for optimization and comparison of detectors. In this paper, a modified nonprewhitening matched filter (NPWE) model observer was developed and validated to predict object detectability for an amorphous selenium (a-Se) direct flat-panel imager (FPI) where aliasing is severe. A preclinical a-Se digital mammography FPI with 85 microm pixel size was used in this investigation. Its physical imaging properties including modulation transfer function (MTF), noise power spectrum, and DQE were fully characterized. An observer performance study was conducted by imaging the CDMAM 3.4 contrast-detail phantom designed specifically for digital mammography and presenting these images to a panel of seven observers. X-ray attenuation and scatter due to the phantom were determined experimentally for use in development of the model observer. The observer study results were analyzed via threshold averaging and signal detection theory (SDT) based techniques to produce contrast-detail curves where threshold contrast is plotted as a function of disk diameter. Validity of the model was established using SDT analysis of the experimental data. The effect of aliasing on the detectability of small diameter disks was determined using the NPWE model observer. The signal spectrum was calculated using the presampling MTF of the detector with and without including the aliased terms. Our results indicate that the NPWE model based on Fourier domain parameters provides reasonable prediction of object detectability for the signal-known-exactly task in uniform image noise for a-Se direct FPI.     

Comparison of low contrast detectability between a digital amorphous silicon and a screen-film based imaging system for thoracic radiography.

Low contrast threshold detectability is investigated theoretically and experimentally for an amorphous silicon (a-Si) x-ray detector designed for digital radiography and for a standard thoracic screen-film combination. A theoretical signal-to-noise ratio is described with a human observer signal detection model. It is calculated using the detective quantum efficiency (DQE) and the modulation transfer function of the imaging system, as well as a spatial response function for the human visual system. Using a four-alternative forced choice observer perception study, the threshold contrasts of disk shaped objects ranging in size from 0.3 to 4.0 mm are determined. Significantly better contrast detectability is obtained from the digital detector, which is attributed to its higher DQE. On average, the disk shaped objects can be detected at 45% less contrast than required for screen-film. With no free parameters the experimental data agree well with those predicted by the observer model. Based upon the data, the model predicts that x-ray exposure for the a-Si detector only needs to be 30% of the exposure required to perform equally well on the contrast-detail detectability task using screen-film.     

New low-contrast resolution phantoms for computed tomography

Computed tomography (CT) has been established as a major imaging modality in diagnostic radiology. Accordingly, acceptance testing and quality control of CT scanners is of great importance. While most procedures and phantoms for testing are widely accepted, there is still discussion and uncertainty about low-contrast (LC) sensitivity. In our opinion this unsatisfactory situation is caused at least in part by the lack of suitable phantoms for LC resolution measurements. We investigated the commonly used phantoms for LC detectability, the Catphan, and for LC resolution, the ATS phantom. While the Catphan showed stable object contrasts, the ATS phantom's measured contrast exhibited a strong dependence on temperature and x-ray quality. Based on newly developed polyurethane resin materials, we designed and tested a LC resolution phantom with several different contrast steps. The object contrasts showed no dependence on temperature and beam quality. The new LC resolution phantom proved to be very suitable for measuring a scanner's low-contrast sensitivity in the image plane, one of the most important image quality parameters. To assess LC resolution in three dimensions we designed an additional phantom with rows of spherical objects. A first prototype was evaluated in a multicenter study. The setup proved to be very helpful to quantify the in-plane and axial LC sensitivity of spiral CT scan modes.     

Automated analysis of phantom images for the evaluation of long-term reproducibility in digital mammography

The performance of an automatic software package was evaluated with phantom images acquired by a full-field digital mammography unit. After the validation, the software was used, together with a Leeds TORMAS test object, to model the image acquisition process. Process modelling results were used to evaluate the sensitivity of the method in detecting changes of exposure parameters from routine image quality measurements in digital mammography, which is the ultimate purpose of long-term reproducibility tests. Image quality indices measured by the software included the mean pixel value and standard deviation of circular details and surrounding background, contrast-to-noise ratio and relative contrast; detail counts were also collected. The validation procedure demonstrated that the software localizes the phantom details correctly and the difference between automatic and manual measurements was within few grey levels. Quantitative analysis showed sufficient sensitivity to relate fluctuations in exposure parameters (kV(p) or mAs) to variations in image quality indices. In comparison, detail counts were found less sensitive in detecting image quality changes, even when limitations due to observer subjectivity were overcome by automatic analysis. In conclusion, long-term reproducibility tests provided by the Leeds TORMAS phantom with quantitative analysis of multiple IQ indices have been demonstrated to be effective in predicting causes of deviation from standard operating conditions and can be used to monitor stability in full-field digital mammography.     

Low-dose CT reconstruction via edge-preserving total variation regularization

High radiation dose in computed tomography (CT) scans increases the lifetime risk of cancer and has become a major clinical concern. Recently, iterative reconstruction algorithms with total variation (TV) regularization have been developed to reconstruct CT images from highly undersampled data acquired at low mAs levels in order to reduce the imaging dose. Nonetheless, the low-contrast structures tend to be smoothed out by the TV regularization, posing a great challenge for the TV method. To solve this problem, in this work we develop an iterative CT reconstruction algorithm with edge-preserving TV (EPTV) regularization to reconstruct CT images from highly undersampled data obtained at low mAs levels. The CT image is reconstructed by minimizing energy consisting of an EPTV norm and a data fidelity term posed by the x-ray projections. The EPTV term is proposed to preferentially perform smoothing only on the non-edge part of the image in order to better preserve the edges, which is realized by introducing a penalty weight to the original TV norm. During the reconstruction process, the pixels at the edges would be gradually identified and given low penalty weight. Our iterative algorithm is implemented on graphics processing unit to improve its speed. We test our reconstruction algorithm on a digital NURBS-based cardiac-troso phantom, a physical chest phantom and a Catphan phantom. Reconstruction results from a conventional filtered backprojection (FBP) algorithm and a TV regularization method without edge-preserving penalty are also presented for comparison purposes. The experimental results illustrate that both the TV-based algorithm and our EPTV algorithm outperform the conventional FBP algorithm in suppressing the streaking artifacts and image noise under a low-dose context. Our edge-preserving algorithm is superior to the TV-based algorithm in that it can preserve more information of low-contrast structures and therefore maintain acceptable spatial resolution.     

Dose efficiency and low-contrast detectability of an amorphous silicon x-ray detector for digital radiography

The effect of dose reduction on low-contrast detectability is investigated theoretically and experimentally for a production grade amorphous silicon (a-Si) x-ray detector and compared with a standard thoracic screen-film combination. A non-prewhitening matched filter observer model modified to include a spatial response function and internal noise for the human visual system (HVS) is used to calculate a signal-to-noise ratio (SNR) related to object detectability. Other inputs to the SNR calculation are the detective quantum efficiency (DQE) and the modulation transfer function (MTF) of the imaging system. Besides threshold detectability, the model predicts the equivalent perception dose ratio (EPDR), which is the fraction of the screen film exposure for which the digital detector provides equal detectability. Images of a contrast-detail phantom are obtained with the digital detector at dose levels corresponding to 27%, 41%, 63% and 100% of the dose used for screen-film. The images are used in a four-alternative forced choice (4-AFC) observer perception study in order to measure threshold detectability. A statistically significant improvement in contrast detectability is measured with the digital detector at 100% and 63% of the screen-film dose. There is no statistical difference between screen-film and digital at 41% of the dose. On average, the experimental EPDR is 44%, which agrees well with the model prediction of 40%.     

Effect of Dose Reduction on the Ability of Digital Mammography to Detect Simulated Microcalcifications

The purpose of this article was to report the relationship between radiation dose and the ability of sentence digital mammography to detect microcalcifications. All images were acquired by computed radiography and an anthropomorphic breast phantom. The tube voltage and anode/filter combination used were 28 kVp and Mo/Mo. Simulated microcalcifications with an approximate diameter of 250–350 μm were positioned on the phantom. Groups of six microcalcifications were arranged in one of two patterns, a line cluster 1 cm long or a hexagonal cluster 4 mm wide. One of the six microcalcifications was removed to create a negative control. Each cluster was placed on 25 different points. Four levels of milliampere-second (mAs) values were applied: 100%, 50%, 25%, and 12.5%. Five staff radiologists participated in an observer performance test. All observers used a workstation with a 3-megapixel monochrome LCD monitor. The areas under the receiver-operating characteristics curves (AUC) were used to compare diagnostic performance among the four doses. The overall AUC scores were 0.97 with 100% mAs, 0.93 (n.s.) with 50%, 0.90 (p < 0.05) with 25%, and 0.81 (p < 0.01) with 12.5% mAs. Among the negative series, the percentage of images on which observers were able to identify the removed microcalcification point decreased from 88.8% with 100% mAs to 83.6% (n.s.) with 50%, 74.8% (p < 0.001) with 25%, and 67.2% (p < 0.001) with 12.5% mAs. A certain level of dose reduction in digital mammography may be an option.Key words: Digital mammography, computed radiography, observer performance, radiation dose, ROC-based analysis, phantoms, imaging     

Automatic polyp detection and measurement with computed tomographic colonography: A phantom study

Purpose

The purpose of this study is to assess the performance of computer-aided detection (CAD) software in detecting and measuring polyps for CT Colonography, based on an in vitro phantom study.

Material and methods

A colon phantom was constructed with a PVC pipe of 3.8 cm diameter. Nine simulated polyps of various sizes (3.2mm-25.4mm) were affixed inside the phantom that was placed in a water bath. The phantom was scanned on a 64-slice CT scanner with tube voltage of 120 kV and current of 205 mAs. Two separate scans were performed, with different slice thickness and reconstruction interval. The first scan (thin) had a slice thickness of 1mm and reconstruction interval 0.5mm. The second scan (thick) had a slice thickness of 2mm and reconstruction interval of 1mm. Images from both scans were processed using CT Colonography software that automatically segments the colon phantom and applies CAD that automatically highlights and provides the size (maximum and minimum diameters, volume) of each polyp. Two readers independently measured each polyp (two orthogonal diameters) using both 2D and 3D views. Readers’ manual measurements (diameters) and automatic measurements from CAD (diameters and volume) were compared to actual polyp sizes as measured by mechanical calipers.

Results

All polyps except the smallest (3.2mm) were detected by CAD. CAD achieved 100% sensitivity in detecting polyps ≥6mm. Mean errors in CAD automated volume measurements for thin and thick slice scans were 8.7% and 6.8%, respectively. Almost all CAD and manual readers’ 3D measurements overestimated the size of polyps to variable extent. Both over- and underestimation of polyp sizes were observed in the readers’ manual 2D measurements. Overall, Reader 1 (expert) had smaller mean error than Reader 2 (non-expert).

Conclusion

CAD provided accurate size measurements for all polyps, and results were comparable to the two readers'' manual measurements     

CHO数字观察器在PET图像滤波方法评估中的应用

在临床应用中需要限制扫描时间和药物剂量,这往往会使正电子发射断层扫描(PET)的图像的分辨率变低,噪声变多。为提供可供临床诊断的图像,去噪是一个必须的手段,而在重建后增加一个滤波器是目前最常用的去噪方法。因此对不同滤波器滤波效果的比较是PET图像重建中的重要环节,其中最关键的是滤波参数的选取。目前采用的信噪比(SNR)以及恢复系数(RC)等评估方法可以用来非定量地选取参数,研究者们只能凭经验选取最优参数。而通道化霍特林观察器(CHO)作为一个比较通用的数字观察器,已被用于与PET图像质量相关的各种参数的选择,如重建算法参数、系统设计参数、临床协议参数等,然而其在评估不同滤波方法对图像重建质量的影响中的应用研究还比较少。通过比较CHO计算得到的ROC(receiver operating characteristic)曲线下面积(area under the ROC curve,AUC),选择两种常用的滤波器(即高斯滤波器和非局部均值(Non-Local Mean, NLM)滤波器)的最优参数,并评估它们在PET中的滤波效果。结果表明,对于13 mm球体,σ为1.1~1.4的高斯滤波器和f为0.5~0.9的NLM滤波器可以达到最大的检测能力值,而对于10 mm球体,σ为1.4~2.0的高斯滤波器和f为0.5~0.9的NLM滤波器可以达到最大的检测能力值。虽然两个滤波器所对应的AUC值都能高达0.9,但是NLM滤波器的AUC值高于高斯滤波器。通过IEC图像和病人图像也能发现,NLM滤波后的PET图像中的亮点比高斯滤波的更加清晰,噪声更少。该结论和传统滤波器评估方法得到的结论一致,这说明在PET的病灶检测任务中,CHO能够准确地比较这两种滤波器的性能。     

Experimental investigation of the dose and image quality characteristics of a digital mammography imaging system

Our purpose in this study was to investigate the image quality and absorbed dose characteristics of a digital mammography imaging system with a CsI scintillator, and to identify an optimal x-ray tube voltage for imaging simulated masses in an average size breast with 50% glandularity. Images were taken of an ACR accreditation phantom using a LORAD digital mammography system with a Mo target and a Mo filter. In one experiment, exposures were performed at 80 mAs with x-ray tube voltages varying between 24 and 34 kVp. In a second experiment, the x-ray tube voltage was kept constant at 28 kVp and the technique factor was varied between 5 and 500 mAs. The average glandular dose at each x-ray tube voltage was determined from measurements of entrance skin exposure and x-ray beam half-value layer. Image contrast was measured as the fractional digital signal intensity difference for the image of a 4 mm thick acrylic disk. Image noise was obtained from the standard deviation in a uniformly exposed region of interest expressed as a fraction of the background intensity. The measured digital signal intensity was proportional to the mAs and to the kVp5.8. Image contrast was independent of mAs, and dropped by 21% when the x-ray tube voltage increased from 24 to 34 kVp. At a constant x-ray tube voltage, image noise was shown to be approximately proportional to (mAs)(-05), which permits the image contrast to noise ratio (CNR) to be modified by changing the mAs. At 80 mAs, increasing the x-ray tube voltage from 24 to 34 kVp increased the CNR by 78%, and increased the average glandular dose by 285%. At a constant lesion CNR, the lowest average glandular dose value occurred at 27.3 kVp. Increasing or decreasing the x-ray tube voltage by 2.3 kVp from the optimum kVp increased the average glandular dose values by 5%. These results show that imaging simulated masses in a 4.2 cm compressed breast at approximately 27 kVp with a Mo/Mo target/filter results in the lowest average glandular dose.     

A search for improved technique factors in paediatric fluoroscopy

A Monte Carlo computational model of a fluoroscopic imaging chain was used for deriving optimal technique factors for paediatric fluoroscopy. The optimal technique was defined as the one that minimizes the absorbed dose (or dose rate) in the patient with a constraint of constant image quality. Image quality was assessed for the task of detecting a detail in the image of a patient-simulating phantom, and was expressed in terms of the ideal observer's signal-to-noise ratio (SNR) for static images and in terms of the accumulating rate of the square of SNR for dynamic imaging. The entrance air kerma (or air kerma rate) and the mean absorbed dose (or dose rate) in the phantom quantified radiation detriment. The calculations were made for homogeneous phantoms simulating newborn, 3-, 10- and 15-year-old patients, barium and iodine contrast material details, several x-ray spectra, and for imaging with or without an antiscatter grid. The image receptor was modelled as a CsI x-ray image intensifier (XRII). For the task of detecting low- or moderate-contrast iodine details, the optimal spectrum can be obtained by using an x-ray tube potential near 50 kV and filtering the x-ray beam heavily. The optimal tube potential is near 60 kV for low- or moderate-contrast barium details, and 80-100 kV for high-contrast details. The low-potential spectra above require a high tube load, but this should be acceptable in paediatric fluoroscopy. A reasonable choice of filtration is the use of an additional 0.25 mm Cu, or a suitable K-edge filter. No increase in the optimal tube potential was found as phantom thickness increased. With the constraint of constant low-contrast detail detectability, the mean absorbed doses obtained with the above spectra are approximately 50% lower than those obtained with the reference conditions of 70 kV and 2.7 mm Al filter. For the smallest patient and x-ray field size, not using a grid was slightly more dose-efficient than using a grid, but when the patient size and field size were increased a fibre interspaced grid resulted in lower doses than imaging without a grid. For a 15-year-old patient the mean absorbed doses were up to 40% lower with this grid than without the grid.     

呼吸时相和模体位置对4D-CT图像均匀性的影响

目的:通过分析Catphan模体CTP486模块的四维CT（4D-CT）图像,研究呼吸时相和模体位置对4D-CT图像均匀性的影响。方法:使用西门子Sensation Open CT模拟机和瓦里安RPM系统,获取Catphan 504模体CTP486模块的4D-CT图像。对3种模体位置情形进行研究。情形A:模体悬空放置;情形T:模体下有一个碳纤维CT平板床;情形B+T:模体下有一个碳纤维固定底板和一个CT平板床。每种情形重复3次测量。使用DoseLab放疗验证软件分析4D-CT图像,得到图像均匀性值（U）。对每一套4D-CT重建10个呼吸时相（0%, 10%,[ ?], 90%）的CT图像序列。对U值进行统计分析和比较研究。结果:4D-CT所有CT序列的U值均小于5 HU。每种情形得到30个U值数据,情形A、情形T和情形B+T 3种情形U值分别为（1.44±0.79）、（1.91±0.91）和（1.99±0.77） HU。统计U值对应的4个边缘感兴趣区（ROI）出现的次数和比例,情形T中ROI 22（9点钟方向）出现13次（占比43.33%）,情形B+T中ROI 23（12点钟方向）出现13次（占比43.33%）。结论:本研究中4D-CT的图像均匀性满足使用要求,不同呼吸时相4D-CT的图像均匀性值不同,模体摆放位置对4D-CT的图像均匀性有一定影响,碳纤维CT平板床和固定底板的存在使4D-CT 的图像均匀性变差。     

DoseLab软件检测CT图像噪声的程序改进及应用分析

目的:对DoseLab软件进行程序改进,增加检测CT图像噪声的功能,对改进的程序进行测试分析。方法:首先,通过使用圆的内接多边形顶点位置计算公式,得到圆内接正三十二边形顶点坐标值。然后,在DoseLab软件Catphan 504模体CTP486模块的图像分析程序中,添加一个正三十二边形的感兴趣区（ROI）,用于检测CT图像噪声。选取2018年每月由西门子CT模拟机日常质量检测（DQC）程序得到的水模体两个层面（S3和S4）的CT图像,对DoseLab改进程序进行测试。对DoseLab改进程序和DQC程序得到的CT图像噪声数据,进行统计分析和比较研究。结果:根据公式计算得到了半径4 cm圆的内接正三十二边形的32个顶点的坐标值,该多边形ROI的面积为49.94 cm2。计算DoseLab改进程序和DQC程序得到的CT图像噪声的差异（ΔN）。在120 kV情形,S3和S4层的ΔN值分别为（0.06±0.07） HU和（0.03±0.09） HU;在140 kV情形,S3和S4层的ΔN值分别为（0.10±0.09） HU和（0.08±0.09） HU。结论:通过添加正三十二边形ROI得到的DoseLab改进程序,可以自动分析水模体和Catphan模体,得到CT图像噪声数据。     

Grid removal and impact on population dose in full-field digital mammography

The study purpose was to determine the impact of anti-scatter grid removal on patient dose, in full field digital mammography. Dose saving, phantom based, was evaluated with the constraint that images acquired with and without grid would provide the same contrast-to-noise ratio (CNR). The digital equipment employed a flat panel detector with cesium iodide for x-ray to light conversion, 100 microm pixel size; the x-ray source was a dual-track tube with selectable filtration. Poly(methyl-emathocrylate) (PMMA) layers in the range 20-70 mm were used to simulate the absorption of different breast thickness, while two Al foils, 0.1 and 0.2 mm thick were used to provide a certain CNR. Images with grid were acquired with the same beam quality as selected in full automatic exposure mode and the mAs levels as close as possible, and the CNR measured for each thickness between 20 and 70 mm. Phantom images without grid were acquired in manual exposure mode, by selecting the same anode/filter combination and kVp as the image with grid at the same thickness, but varying mAs from 10 to 200. For each thickness, an image without aluminum was acquired for each mAs value, in order to obtain a flat image to be used to subtract the scatter nonuniformity from the phantom images. After scatter subtraction, the CNR was measured on images without grid. The mAs value that should be set to acquire a phantom image without grid so that it has the same CNR as the corresponding grid image was calculated. Therefore, mAs reduction percentage was determined versus phantom thickness. Results showed that dose saving was lower than 30% for PMMA equivalent breast thinner than 40 mm, decreased below 10% for intermediate thickness (45-50 mm), but there was no dose gain for thickness beyond 60 mm. By applying the mAs reduction factors to a clinical population derived from a data base of 4622 breasts, dose benefit was quantified in terms of population dose. On the average, the overall dose reduction was about 8%. It was considered small, not sufficient to justify a clinical implementation, and the anti-scatter grid was maintained.     

Receiver operating characteristic and location analysis of simulated near-infrared tomography images

Receiver operating characteristic (ROC) analysis was performed on simulated near-infrared tomography images, using both human observer and contrast-to-noise ratio (CNR) computational assessment, for application in breast cancer imaging. In the analysis, a nonparametric approach was applied for estimating the ROC curves. Human observer detection of objects had superior capability to localize the presence of heterogeneities when the objects were small with high contrast, with a minimum detectable threshold of CNR near 3.0 to 3.3 in the images. Human observers were able to detect heterogeneities in the images below a size limit of 4 mm, yet could not accurately find the location of these objects when they were below 10 mm diameter. For large objects, the lower limit of a detectable contrast limit was near 10% increase relative to the background. The results also indicate that iterations of the nonlinear reconstruction algorithm beyond 4 did not significantly improve the human detection ability, and degraded the overall localization ability for the objects in the image, predominantly by increasing the noise in the background. Interobserver variance performance in detecting objects in these images was low, suggesting that because of the low spatial resolution, detection tasks with NIR tomography is likely consistent between human observers.     

Stereomammography: evaluation of depth perception using a virtual 3D cursor

We are evaluating the usefulness of stereomammography in improving breast cancer diagnosis. One area that we are investigating is whether the improved depth perception associated with stereomammography might be significantly enhanced with the use of a virtual 3D cursor. A study was performed to evaluate the accuracy of absolute depth measurements made in stereomammograms with such a cursor. A biopsy unit was used to produce digital stereo images of a phantom containing 50 low contrast fibrils (0.5 mm diam monofilaments) at depths ranging from 1 to 11 mm, with a minimum spacing of 2 mm. Half of the fibrils were oriented perpendicular (vertical) and half parallel (horizontal) to the stereo shift direction. The depth and orientation of each fibril were randomized, and the horizontal and vertical fibrils crossed, simulating overlapping structures in a breast image. Left and right eye images were generated by shifting the x-ray tube from +2.5 degrees to -2.5 degrees relative to the image receptor. Three observers viewed these images on a computer display with stereo glasses and adjusted the position of a cross-shaped virtual cursor to best match the perceived location of each fibril. The x, y, and z positions of the cursor were indicated on the display. The z (depth) coordinate was separately calibrated using known positions of fibrils in the phantom. The observers analyzed images of two configurations of the phantom. Thus, each observer made 50 vertical filament depth measurements and 50 horizontal filament depth measurements. These measurements were compared with the true depths. The correlation coefficients between the measured and true depths of the vertically oriented fibrils for the three observers were 0.99, 0.97, and 0.89 with standard errors of the estimates of 0.39 mm, 0.83 mm, and 1.33 mm, respectively. Corresponding values for the horizontally oriented fibrils were 0.91, 0.28, and 0.08, and 1.87 mm, 4.19 mm, and 3.13 mm. All observers could estimate the absolute depths of vertically oriented objects fairly accurately in digital stereomammograms; however, only one observer was able to accurately estimate the depths of horizontally oriented objects. This may relate to different aptitudes for stereoscopic visualization. The orientations of most objects in actual mammograms are combinations of horizontal and vertical. Further studies are planned to evaluate absolute depth measurements of fibrils oriented at various intermediate angles and of objects of different shapes. The effects of the shape and contrast of the virtual cursor and the stereo shift angle on the accuracy of the depth measurements will also be investigated.     

Diffuse optical tomographic reconstruction using multifrequency data

We investigated the use of multifrequency diffuse optical tomography (MF-DOT) data for the reconstruction of the optical parameters. The experiments were performed in a 63 mm diameter cylindrical phantom containing a 15 mm diameter cylindrical object. Modulation frequencies ranging from 110 MHz to 280 MHz were used in the phantom experiments changing the contrast in absorption of the object with respect to the phantom while keeping the scattering value the same. The diffusion equation was solved using the finite element method. The sensitivity information from each frequency was combined to form a single Jacobian. The inverse problem was solved iteratively by minimizing the difference between the measurements and forward problem using single and multiple modulation frequency data. A multiparameter Tikhonov scheme was used for regularization. The phantom results show that the peak absorption coefficient in a region of interest was obtained with an error less then 5% using two-frequency reconstruction for absorption contrast values up to 2.2 times higher than background and 10% for the absorption contrast values larger than 2.2. The use of two-frequency data is sufficient to improve the quantitative accuracy compared with the single frequency reconstruction with appropriate selection of these frequencies.     

Effects of magnification and zooming on depth perception in digital stereomammography: an observer performance study

We are evaluating the application of stereoscopic imaging to digital mammography. In the current study, we investigated the effects of magnification and zooming on depth perception. A modular phantom was designed which contained six layers of 1-mm-thick Lexan plates, each spaced 1 mm apart. Eight to nine small, thin nylon fibrils were pasted on each plate in horizontal or vertical orientations such that they formed 25 crossing fibril pairs in a projected image. The depth separation between each fibril pair ranged from 2 to 10 mm. A change in the order of the Lexan plates changed the depth separation of the two fibrils in a pair. Stereoscopic image pairs of the phantom were acquired with a GE full-field digital mammography system. Three different phantom configurations were imaged. All images were obtained using a Rh target/Rh filter spectrum at 30 kVp tube potential and a +/- 3 stereo shift angle. Images were acquired in both contact and 1.8X magnification geometry and an exposure range of 4 to 63 mAs was employed. The images were displayed on a Barco monitor driven by a Metheus stereo graphics board and viewed with LCD stereo glasses. Five observers participated in the study. Each observer visually judged whether the vertical fibril was in front of or behind the horizontal fibril in each fibril pair. It was found that the accuracy of depth discrimination increased with increasing fibril depth separation and x-ray exposure. The accuracy was not improved by electronic display zooming of the contact stereo images by 2X. Under conditions of high noise (low mAs) and small depth separation between the fibrils, the observers' depth discrimination ability was significantly better in stereo images acquired with geometric magnification than in images acquired with a contact technique and displayed with or without zooming. Under our experimental conditions, a 2 mm depth discrimination was achieved with over 60% accuracy on contact images with and without zooming, and with over 90% accuracy on magnification images. This study indicates that stereoscopic imaging, especially with magnification, may be useful for visualizing the spatial distribution of microcalcifications in a cluster and for differentiating overlapping tissues from masses on mammograms.     

Effect of dose metrics and radiation risk models when optimizing CT x-ray tube voltage

We investigated the effect of different CT dose metrics, as well as the implications of various radiation risk models, on the optimization of x-ray tube voltage (kV) in CT. Soft tissue attenuation characteristics and noise levels, obtained from CT scans of a Rando phantom, were used to compute contrast-to-noise ratios (CNR) at x-ray tube voltages between 80 and 140 kV. Four CT dose metrics were evaluated: (a) CTDI(air), (b) weighted CTDI(w), (c) organ dose (D(organ)), and (d) effective dose (E). All doses were obtained using the ImPACT CT Dosimetry software package. Soft tissue CNR was adjusted by the modification of the mAs by assuming that CNR(2) was proportional to mAs. Optimization criteria were: (a) maintaining a constant CNR at each kV and identifying the value that minimizes patient dose; and (b) maintaining a constant dose at each kV and identifying the value that maximizes CNR. We also investigated the implication for optimization strategies assuming that radiation risk is proportional to E(n), with n varying between 0 and 2. Optimizing with respect to phantom measurements (i.e., CTDI(air) and CTDI(w)) could generate results that differed quantitatively and qualitatively from those obtained using patient doses (i.e., D(organ) and E). For head CT scans, 140 kV offered the lowest patient doses as well as the highest CNR, whereas in abdominal scans 80 kV was optimal. Use of an optimal kV for CT imaging over current practice of using 120 kV might reduce patient doses by 10-15%, or improve CNR by 5-10%. Assuming that the risk was proportional to E(n) made no difference to the optimal kV for positive values of n up to 2. We conclude that (a) CT optimization with respect to kV should generally be performed with respect to the patient effective dose, (b) neither CTDI(air) nor the body CTDI(w) are appropriate for use in CT optimization, (c) the range of current radiation risk models should not affect the optimal kV value in CT imaging.