1. Field of the Invention
The disclosure in the present application is in the field of 3-D volumetric medical imaging, and, in particular, is directed to an improved device and method to determine tissue properties in the body of a subject.
2. Description of the Related Art
Many diagnostic 3-D imaging devices are used in modem medicine for many types of exams, which are performed by subjective viewing of the images on either film or more commonly on electronic displays. This subjective viewing uses assumed quantitative image pixels, which define boundaries of tissues, organs or foreign masses, and discrimination of tissue types based on grey scale or color images, representative of the underlying tissue properties. Identification of diagnostic details is fundamentally dependent upon the detection of image detail edges. Diagnostic interpretation or measurements from these images assume an appropriate relationship of the image to the tissue property. All of these imaging devices however have significant limitations for quantitative measurements because the relationship is not sufficiently accurate or defined. The embodiments disclosed in the present application are directed at improving these limitations while reducing costs, time and effort in the clinical setting.
Measurement, display, and analysis of tissue properties from medical images have many diagnostic benefits in the living subject. These include such measures as density, mass, volume, edges, etc.; image display reference values, translucency, texture; tissue water or hydrogen content, electron density, blood vessel iodine contrast media concentration and density or blood flow, tissue iron content, fat/muscle ratio, air/tissue ratio, and the like. Modern imaging devices enable the potential to identify and quantify, either automatically or manually, most organs and tissues of the body. The imaging devices to which the improvements disclosed herein are directed provide a set of volumetric images acquired from multiple projection angles, which allow reconstruction of images in various planes. Such devices include dual-energy and single-energy CT scanners, rotational C-arms, MRI devices, x-ray tomosynthesis and 3D DXA (dual-energy X-absorptometry) in which the x-ray source is moved to provide various projections, such as a Hologic 3-D DXA device in development. Several tissue measurements of interest include bone density, lung nodule density, cardiac, aortic and vascular calcifications, vascular soft plaque, fat measurements, muscle mass, lung volume and density (emphysema and the like), liver iron content, perfusion and blood flow, organ volume, density and mass, contrast angiography and the like. Tissue corrections for radiotherapy dose calculations are based on the electron density. Accurate and conveniently available measures of true tissue properties and their change with disease conditions or therapy allow diagnostic analysis of images and new diagnostic criteria not currently possible. The methods disclosed in this application provide improved accuracy and greater ease in all of these potential applications.
Radiologists routinely make subjective and even quantitative measurements of foreign masses, tissues or organs by manually placing cursors to define the 2-D extent of the target. If the window and/or level (brightness and contrast) are changed in the display, the apparent size of the target changes because the boundary is not discrete and is moved in or out of the display range. Thus, the measured object size is frequently inaccurate and will vary from operator to operator and from scanner to scanner depending on the display conditions and scanner properties. In addition, the process to set and adjust the window and level requires operator time and is currently very inefficient. Electronic image data are frequently erased, and only the films retained for the medical records.
Prior art methods have allowed CT scanners to be used as quantitative instruments for bone density measurements in quantitative computerized tomography (QCT) by the use of calibration phantoms scanned simultaneously with the patient (simultaneous calibration). Such phantoms have greatly aided the standardization accuracy and reproducibility of bone density measurements. Representative prior art methods have been disclosed in U.S. Pat. No. 4,233,507 to Volz), U.S. Pat. No. 4,985,906 to Arnold, U.S. Pat. No. 5,335,260 to Arnold, U.S. Pat. No. 4,782,502 to Schulz, U.S. Pat. No. 4,651,335 to Kalender, U.S. Pat. No. 4,870,666 to Lonn, and others. These prior phantom designs included, for example, samples within pillows and positioned beside the subject; a flexible phantom positioned on top of the subject; a flexible phantom positioned under the subject with means to force contact to patient; samples within a slot between two couch pads; and in a removable, rigid structure or table top section in the original Volz patent.
The Volz patent used samples with 5 cm cross-sectional diameters and subsequent improvement patents used similar or larger sample areas because of the need to have sufficient numbers of pixels in each CT slice region of interest (ROI) for statistical measurements. The base material and calibration samples were large and dense and increased the dose and image artifacts. U.S. Pat. No. 4,985,906 to Arnold discloses improvements to prior phantoms by using tissue equivalent base material with reduced size and mass to minimize beam hardening effects.
In U.S. Pat. No. 6,990,222 to Arnold discloses a method for hybrid calibration using simultaneous phantom calibration along with internal tissue references of the individual patient. U.S. Pat. No. 6,990,222 is incorporated by reference herein in its entirety. More recently, fast multi-detector MDCT scanners have been used for coronary and aortic calcium analysis with or without simultaneous phantom calibration.
The prior calibration methods were successful because CT numbers, (Hounsfield Units, HU), are only estimates of the x-ray attenuation coefficients of tissue relative to water as the manufacturer's calibration reference material. CT numbers fail to be truly quantitative for several reasons. For example, the tissue attenuation coefficients are photon energy dependent, and the x-ray beam energy spectra are not measured or known for individual patients. Further, many beam energy spectra exist in each CT slice (e.g., a unique spectrum for each path length through the patient and seen at a particular detector element, and a unique spectrum for each view through the patient). The beam spectrum changes with the thickness and composition of tissues in the path length. The quantities of fat, soft tissue, air, and bone vary with each projection. X-ray tube filtration to shape the beam intensity also changes the beam spectrum resulting in variations in tissue attenuation based on locations within the field of view. Scattered radiation is also variable and dependent on some of these parameters. Manufacturers' calibrations (for example, CT number and beam hardening corrections in current practice and scatter correction) are based on idealized phantoms scanned independently from the patient, which are often circular in shape and composed of water, plastics, or other homogeneous, synthetic materials. These differ significantly from the shape and varied composition of real patients. Image pixel intensities vary from image to image, and are dependent on table height, position in the beam, scanner drift, tube changes, manufacture reconstruction software, body region thickness and volume, field of view, and sometimes even the time of day as the imagers warm up.
Measurement errors also result from organ motion artifacts, tissue heterogeneity (fat, muscle, blood mixtures in sub-voxel volumes), gantry wobble and vibrations of the x-ray source, escape of K x-rays in the detectors, reconstruction algorithm, electronic noise, local beam hardening, scattered radiation and structure artifacts (usually from the presence of larger bones, air volumes and table structures). As a result of the measurement errors, tissue densities vary with imaging device, image acquisition parameters, gating and motion, etc., some of which have been reported. In addition, measurements will vary with patient size and body type which is poorly recognized. Because of the finite detector element size, the x-ray tube focal spot size and geometry, source movement during scanning and image back projection reconstruction imperfections, the finite number of views and source and/or table movement, the reconstruction of objects is not exact resulting in edge blurring by the point spread function (PSF) and loss of accurate tissue density representations at edges of objects. As a result, bone density and calcium scores vary significantly with different devices, over time and between different institutions, patient body composition and imaging techniques. The above-referenced measurements are also dependent on this listing of sources of errors.
Prior methods for simultaneous CT calibrations in bone density measurements have used calcium phantoms which are removed after each exam. This requires the operator to place the phantom on the CT table before and to remove the phantom after each quantitative study and then place the phantom in storage. Further, the operator is required to position the phantom in a surrounding foam pad for patient comfort and then to position the phantom and pad and patient in the CT field-of-view (FOV). It is important to reproduce the position of the phantom in the FOV and to avoid its movement during the scans. Misalignments cause variable results and failure to reproduce positioning on follow-up exams months or years later can reduce measurement precision. Unlike these prior methods, the device and methods disclosed in the present application overcomes all of these limitations. Further, the operator does not need to be concerned with phantom positioning or repositioning and thus saves time during each patient procedure.
With prior methods, any study without the phantom cannot then or later be analyzed. With the device and method disclosed herein, any prior study from months to years can be analyzed at any later date since the references may always be present in all images and all studies. Increased radiation doses, costs, times and efforts of repeating scans for quantitative analysis are avoided when existing scans made for other exams can be reprocessed. Since many medical conditions such as osteoporosis, coronary or vascular disease, emphysema, body composition, etc., are chronic and slowly changing conditions, later measurements are likely to remain highly relevant even if analyzed months or years later.
Prior art methods employ large cross-sectional area phantoms with higher density base materials. The large sample sizes were required because software methods operated on usually 2 to 4 reconstructed thick slices by manual or later automated methods. (See, for example, U.S. Pat. No. 4,922,915 to Arnold.) These features are now known to cause beam hardening and scatter artifacts and cause streaks from the reference phantom itself. The larger/dense phantom increases radiation dose to the patient and/or reduces image quality by adding noise and scatter. Because of the automated software methods disclosed herein, much smaller, cross-sectional reference samples can be used with great advantage.
Prior methods and phantoms create patient discomfort from the hard material and sharp edges as well as large size. Patient discomfort often leads to patient motion during the scan creating additional errors. The device and method disclosed herein overcome patient discomfort. Automated software methods to locate and measure reproducibly the very small diameter samples were not available in prior methods. Further, finding and measuring these small samples which are now not rigid and which extend over long distances required advanced methods disclosed herein. With the advent of MDCT scanners with several hundred CT images per exam, calibration on all images in background automatically was required.
Prior methods used phantoms of reduced length, which were significantly shorter than the length of the patient's torso. The required large sample sizes and removal of the phantom between studies resulted in short phantoms. The phantom had to be repositioned between multiple scans, such as hip and thoracic vertebral BMD scans. In addition, the operator must be aware of the position of the phantom and/or scout scan the patient to verify phantom positioning under the anatomical region to be scanned. This requires added time and radiation dose to the patient. The disclosed devices and method overcome these limitations.
Prior methods have been limited primarily to bone density measurements using calcium equivalent phantom samples. Additional reference and calibration materials and new calibrations and image corrections are possible with the currently disclosure method and system.
Prior methods to correct CT scanners for gantry, tube and table motion and methods to measure or correct for image reconstruction, detector variations, system motion and wobble, have used phantoms with pins, spheres and small diameter wires contained within idealized phantoms. These phantoms and test objects have been scanned independent of the patients. These methods, while providing the ability to correct for scanner imperfections in idealized conditions, are not necessarily representative of the imaging device when the patients are present and are not representative of the device at the time of specific clinical exams. X-ray imagers and CT scanners are known to drift, change sensitivity with use and environmental conditions, patient size and weight, (CT table movement, backlash and bending), and scanner and scanning parameters. For all of these and additional reasons, current test phantoms and methods for image corrections have limitations because they have been applied without the patient present. The disclosed device and methods by including similar phantom targets for image corrections simultaneously with each individual patient exam and at the relevant time of the exam overcomes many of these limitations.
An aspect in accordance with the invention disclosed herein is a calibration and reference couch pad, patient support device or tabletop that contains a variety of reference sample materials and test targets. The samples and targets are permanently enclosed within the pad, device or tabletop at fixed locations and are scanned simultaneously with patient scans. The samples and targets are small in cross-sectional size and area which provides patient comfort minimizes patient movement and avoids image artifacts. The samples and targets are further manufactured with high reproducibility and homogeneity and of sufficient length or repeating placements to always be present in the image for all torso exams without operator effort or placement or alignment verification. Any scan on any body part and at any time can be post-analyzed for a variety of tissue and image quality measurements. Patient comfort and radiation dose are improved. Operator time and effort is reduced. Costs to patients or facilities for repeat or additional scans at a later date are also reduced.
One aspect of the disclosure is advanced software methods to identify and measure the very small size calibration samples and to do so over extended sample lengths. The samples being of a smaller cross-sectional area and extended length tend to move and flex with patient weight. The automated software methods were required to operate in background without operate assistance and to reproducible identify, center measurement regions and apply the corrections and calibrations in acceptable times on large image sets.
Additional reference samples and test targets facilitate new and innovative analysis, corrections and calibrations for a variety of clinical exams or scanner imperfections. Calibration of in vivo iodine contrast distribution and blood flow can be achieved with simultaneous calibration for the first time. Dual-energy CT or 3-D DXA can be corrected for imperfect energy subtraction and associated algorithms for such purposes as CT angiography, bone removal or tissue type separations. Bone density and vascular calcifications can be calibrated for any region of the body and without patient or phantom repositioning and without operator concern for alignments during the scans.
Image reconstruction imperfections can be analyzed and corrected which are representative of the actual patient and scan of interest. Image quality, motion errors, image blur from the system represented by the device PS, table, gantry and X-ray source motion and backlash can be determined and corrected and under the specific conditions of the scans to improve image quality and quantitative measurements.
Other aspects of the disclosed embodiments are improved calibration and reference devices and methods which allow improvements in current 3-D x-ray exams as well as development of new quantitative exams. In addition, improved references and methods for corrections of device imperfections are disclosed which are device, technique, patient and time specific.
Prior to the disclosure herein, uses of simultaneous calibration phantoms or methods for measurements of soft tissue density, non-calcified vascular plaque (so-called soft plaque), CT angiography, lung density, dimensional measurements or correction of image imperfections were not known. It is desirable to provide improved calibration methods for all the tissues of the body.
Without limiting the scope of the invention, the foregoing aspects and other aspects in accordance with embodiments disclosed herein are described herein in connection with the accompanying drawings, images, and figures in which:
An important determination of the present disclosure is that the energy X-ray attenuation response (slope of the regression calibration equation) can be measured independently from the intercept of the regression calibration equation (i.e., the bias component). An additional determination of this discloser is that the slope is relatively independent of the sample size and surrounding material. The hybrid calibration method in combination with the disclosed device and methods has been shown to be superior to the use of external phantoms only and particularly so when large and dense phantoms are used.
Another important determination of the present disclosure is that Z-axis integration of pixels in the small cross-sectional samples can be used preferably to phantoms with large samples and single slices. Practical use of the presently disclosed device incorporates innovative software methods which automatically detect and measure the small samples, which also may move somewhat from scan to scan. Prior methods measure bone density on single reconstructed CT slices, typically using 2 to 5 slices for each session. The slices were 8-10 mm thick and were positioned or reconstructed through 2 to 5 individual vertebrae. Locating and positioning regions-of-interest (ROI) measurement areas within the larger phantom samples required much less complex software methods. The ROIs could be manually adjusted in practical times, and allowed much slower run times or real time operation with operator present.
Innovative and advanced software methods are used to detect and measure the small samples of the disclosed system without operator assistance and running in background mode. The software methods automatically locate and determine measurement regions within each sample and on every CT slice. CT studies containing hundreds of axial images and even greater than 1,000 images can readily be analyzed without operator interaction. The automated software (“The Phantom Finder”) operates in background mode and can archive region placements on hundreds of images in a few seconds (e.g., on the order of 5 seconds). The auto Phantom Finder is particularly advantageous for small calibration and test samples, for high resolution (small voxel, thin CT slices and/or smaller pixels) images and volumetric scans of large body regions or whole body scans. The software allows variable Z-axis integration, which can be adapted to different scan parameters, different patient size and body mass and resulting beam hardening, scatter and image noise effects.
A representative embodiment of the improved full-torso phantom is shown in
As illustrated in
In certain preferred embodiments, the calibration pad 10 further includes a wire 40 (e.g., a small diameter aluminum wire) that extends for the full length of the calibration pad. In certain preferred embodiments, the calibration pad further includes a metallic or increased density material in the form of a ribbon (or sheet) 42 that extends for the full length of the pad. In certain preferred embodiments, the calibration pad further includes a plurality of metallic pins 44 (one shown in
In an alternative embodiment shown in cross section in
Using the calibration pad 10 of
Calibration measurements are expressed in a regression equation of the form y=Io+S×HUs. The slope, S, of the regression curves is related to the effective beam energy used to make the image. For each different KVps, filtration of the primary beam, and different patient sizes and compositions, the slope varies. The intercept, Io, of the regression equation is a measure of the bias and is determined form the patient's own blood/muscle tissues in the hybrid calibration method. The slope, S, will vary with location in a given patient due to differing amounts and densities of tissue being present in the beam.
Voxels within the boundaries of the reference samples are automatically recorded and are then used to calculate a histogram. The histogram may be smoothed or a Gaussian curve in accordance with known methods. The ROI readings from several neighboring slices may be integrated for improved statistics. The ends of the histogram may be clipped to remove partial volume or blurred voxel readings from the margins of samples.
The calculated mode is next combined with the phantom regression equation to create the final hybrid calibration equation. The calibration equation is an expression of the form:
Td=(HUb−HU′b)+S×HU
In the foregoing expression, Td is the tissue density of the target tissue or organ to be analyzed. The tissue density Td may be expressed as corrected HU values or Td may be represented in tissue density units such as grams per cubic centimeter (g/cc). The parameter HUb is the previously determined or known CT number of representative blood/muscle, or other tissues. The parameter HU′b is the measured reference tissue density, which in this example may be the mode of the histogram of the heart and blood, the mode of fat or a calculated soft tissue representative (STR) calculated from the fat and muscle/blood of the patient. The parameter S is the slope of the regression equation measured from the calibration phantom samples. The parameter HU is the scanner CT number or may be a CT number measure of the target tissue to be analyzed. Note that the target tissue may be coronary calcifications, an organ, an unknown mass, or the edge of a target region or organ.
When the detail is imaged in, for example, a CT scanner (e.g., the CT scanner 14 of
The calcium reference samples 34 and 36 in the calibration pad 10 have a composition closely equivalent to calcifications within the patient since both are largely composed of calcium hydroxyapatite (CaHA). Any change in effective x-ray beam energy from scanner drifts, tube changes, different CT scanners, filtration and highly important, differences in patient compositions are largely corrected for by the calibration process which occurs on every CT image. Likewise the iodine reference sample 38 is closely representative of iodine contrast media injected into the blood vessels of the patient and blood perfusion into tissues.
Dual energy CT (DECT) used in angiography studies and with or without vascular calcifications suffers from imperfect subtractions and imperfect image reconstructions as discussed above. Computation methods to improve the separation of two targets or tissues are needed and are dependent on the specific attenuation of specific patients. The availability of both calcium and iodine references in DECT images with known properties allows corrections to be made for improvements. Dual energy imaging in other devices such as 3-D DXA, tomosynthesis, rotational x-ray systems, and the like, can be improved by like processing with known references present in the images.
Scattered radiation in newer MDCT scanners with 16 to 256 detector rows is a major image degradation factor. By design, the MDCT scanners image a much larger volume of tissue in each rotation resulting in greatly increased scatter. Scatter removal methods including grids, and post processing methods are imperfect resulting in loss of image detail and quantitative accuracy of the post-processed image. The available of reference samples and test objects in each image along with the patients provides known references at locations which can be determined. Reconstructed images can be forward-projected for corrections and calibrations. The processing uses the known references as the gold standard for that image by processing until the desired image properties are achieved. The processed images can then be back-projected again following the corrections. Alternatively, image processing can be performed directly on each CT view angle to achieve corrections based on a set ideal of the known reference. Current scatter correction methods are arbitrary and based on the CT scanner and exposing parameters. Corrections made based on the specific patients anatomy and body composition can improve scatter correction methods.
In a first step 130 of the flow chart 120 of
In a step 136, a lowest standard deviation (stdev) is determined, and a new start location is determined and centered at the center pixel of the ROI with the lowest stdev. The stdev is determined for each of 4 new adjacent ROIs centered at the new start location. As illustrated in a step 140, the steps 132, 134 and 136 are repeated until no ROI position is found with a lower stdev than in the preceding iteration.
The center of the ROI determined in the step 140 is saved as an object center location (Lo) used for additional image processing. In a step 142, a new search ROI is located at Lo. The new search ROI has a radius equal to 4 times the object radius. In a step 144, edge detection is performed in the new search ROI. Then, in a step 146, the Hough Transform is applied with a diameter equal to the known object diameter. Following a first preprocessing step with an edge detection algorithm with constraints optimized to the expected shape and densities of the samples, candidate edges are applied to the Hough Transform to extract the desired feature of sample x-y Center-of-Area (CoA). Deviations from circular shaped edges are automatically excluded in preference for those with best fits. From the selected edge boundaries and satisfying sanity checks, the x, y, and z locations of the selected outlines are determined. A circular ROI of area smaller than the known sample is then centered at the CoA in a step 150. Preferably, the ROI has a radius of approximately 70% of the known object radius (e.g., 0.707 times the known object radius). A goodness-of-fit criteria is calculated based on the standard deviation of the pixel reading within the positioned ROI which serves to verify centering of the ROI at or near the center of the small samples. It is important that pixel readings at and near the edges of the samples are not employed in the final calibrations. The pixel values in the final ROI are recorded in a step 152. As illustrated by a step 154, the preceding steps of the method are repeated for each of the objects which may be present in each image. For example, in the illustrated embodiment of the calibration pad 10, the images of the 4 round samples are located and measured. As illustrated by a step 156, the method is repeated for each of N images which contain the reference objects. The software that implements the method thus records n ROI values per image times the number of images in a study.
As illustrated by a step 160, the pixel values may be summed in a certain number of the final ROIs of adjacent images to reduce noise effects due to the very small samples. Then, in a step 162, a histogram of the pixel values of the sum of n adjacent ROIs of each object is then calculated. After completing the foregoing steps, an additional step of fitting a Gaussian curve to the histogram distribution of pixel readings may be implemented. The two tails of the histogram are preferably clipped to a level of approximately 10%. The final sample ROI reading is the average of the pixel intensities within the so positioned and processed ROI. The mean and stdev of the pixel values are calculated from the sum of each of the ROIs of each object in a step 164. These values are the sample readings used for calibration or correction. The values are determined automatically in background mode without operator interaction. In a step 166, the values are provided as the return values from the method, and the method then stops.
One skilled in art will appreciate that the foregoing embodiments are illustrative of the present invention. The present invention can be advantageously incorporated into alternative embodiments and a variety of imaging devices which produce axial images or volumetric scans while remaining within the spirit and scope of the present invention, as defined by the appended claims.
The present application claims the benefit of priority under 35 USC §119(e) to U.S. Provisional Application No. 61/118,438, filed on Nov. 27, 2008.
Number | Name | Date | Kind |
---|---|---|---|
4233507 | Volz | Nov 1980 | A |
4649561 | Arnold | Mar 1987 | A |
4663772 | Mattson et al. | May 1987 | A |
4724110 | Arnold | Feb 1988 | A |
4782502 | Schulz | Nov 1988 | A |
4870666 | Lonn | Sep 1989 | A |
4922915 | Arnold et al. | May 1990 | A |
4985906 | Arnold | Jan 1991 | A |
5034969 | Ozaki | Jul 1991 | A |
5068788 | Goodenough et al. | Nov 1991 | A |
5222021 | Feldman et al. | Jun 1993 | A |
5235628 | Kalender | Aug 1993 | A |
5335260 | Arnold | Aug 1994 | A |
5442674 | Picard et al. | Aug 1995 | A |
5521955 | Gohno et al. | May 1996 | A |
5577089 | Mazess | Nov 1996 | A |
5696805 | Gaborski et al. | Dec 1997 | A |
5712892 | Weil et al. | Jan 1998 | A |
5757877 | Wilting | May 1998 | A |
5774519 | Lindstrom et al. | Jun 1998 | A |
5782762 | Vining | Jul 1998 | A |
5891030 | Johnson et al. | Apr 1999 | A |
6026142 | Gueziec et al. | Feb 2000 | A |
6226350 | Hsieh | May 2001 | B1 |
6233304 | Hu et al. | May 2001 | B1 |
6243437 | Hu et al. | Jun 2001 | B1 |
6278761 | Kim et al. | Aug 2001 | B1 |
6302582 | Nord et al. | Oct 2001 | B1 |
6320931 | Arnold | Nov 2001 | B1 |
6421552 | Hsieh | Jul 2002 | B1 |
6639965 | Hsieh et al. | Oct 2003 | B1 |
6674834 | Acharya et al. | Jan 2004 | B1 |
6697451 | Acharya et al. | Feb 2004 | B2 |
6789943 | Zapalac | Sep 2004 | B2 |
6990222 | Arnold | Jan 2006 | B2 |
7127096 | Kaufman et al. | Oct 2006 | B2 |
7251306 | Sauer et al. | Jul 2007 | B2 |
7409035 | Kaufman et al. | Aug 2008 | B2 |
7471765 | Jaffray et al. | Dec 2008 | B2 |
7583778 | Mori | Sep 2009 | B2 |
20030048867 | Acharya et al. | Mar 2003 | A1 |
20030095693 | Kaufman et al. | May 2003 | A1 |
Number | Date | Country | |
---|---|---|---|
61118438 | Nov 2008 | US |