SUBSTANCE QUANTIFICATION IN ANIMAL BODIES

Abstract: A system and method is disclosed for utilizing impage pixel value information generated by a digital radiographic system for quantifying the amount of calcium in a particular object (134) of interest within a subject, wherein the subject also con­ tains additional quantities of calcium located at positions (130) such that the additional calcium interferes with the acquisition of data describing the object of interest (134). A cursor (136) configured as a parallelogram is described and pixel values are summed along lines parallel to the sides of the parallelogram. A profile (138) is generated from these summations. A portion of the profile (138) corresponding to pixel lines not intersecting the object of interest (134) is used to estimate the contribution of background and other calcium to the profile as a whole, and this estimated value is subtracted from the total profile. In­ tegration of the remainder of the profile (138) provides a scaler value corresponding to the amount of calcium in the object of interest (134).

Description: This invention relates to an improved method and apparatus especially suited for use in the environment of digital radiography for employing signals normally used for imaging for quantifying the amount of calcium or other known substance in a particular localized object of interest in a body which contains other sites having the same substance at positions such that the other sites interfere with acquisition of data describing the object of interest. Radiography and fluoroscopy are well known imaging techniques. In a conventional radiography system, an x-ray source is actuated to direct a divergent area beam of x-rays through a patient. A cassette containing an x-ray sensitive phosphor screen and light and x-ray sensitive film is positioned in the x-ray path on the side of the patient opposite the source. Radiation passing through the patient's body is attenuated in varying degrees in accordance with the various types of tissue through which the x-rays pass. The attenuated x-rays from the patient emerge in a pattern, and strip the phosphor screen, which in turn exposes the film. The x-ray film is processed to yield a visible image which can be interpreted by a radiologist as defining internal body structure and/or condition of the patient. In conventional fluoroscopy, a continuous or rapidly pulsed area beam of x-rays is directed through the patient's body. An image intensifier tube is positioned in the path of the beam opposite the source with respect to the patient. The image intensifier tube receives the emergent radiation pattern from the patient, and converts it to a small, brightened visible image at an output face. Either a mirror or closed circuit television system views the output face and produces a dynamic real time visual image, such as on a CRT, for interpretation by a radiologist. More recently, digital radiography and fluoroscopy techniques have been developed. In digital radiography, the source directs x-radiation through a patient's body to a detector in the beam path beyond the patient. The detector, by use of appropriate sensor means, responds to incident radiation to produce analog signals representing the sensed radiation image, which signals are converted to digital information and fed to a digital data processing unit. The data processing unit records, and/or processes and enhances the digital data. A display unit responds to the appropriate digital data representing the image to convert the digital information back into analog form and produce a visual display of the patient's internal body structure derived from the acquired image pattern of radiation emergent from the patient's body. The display system can be coupled directly to the digital data processing unit for substantially real time imaging, or can be fed stored digital data from digital storage means such as tapes or discs representing patient images from earlier studies. Digital radiography includes radiographic techniques in which a thin fan beam of x-rays is used. In this technique, often called "scan (or slit) projection radiography" (SPR) a fan beam of x-rays is directed through a patient's body. The fan is scanned across the patient, or the patient is movably interposed between the fan beam x-ray source and an array of individual cellular detector segments which are aligned along an arcuate or linear path. Relative movement is effected between the source-detector arrangement and the patient's body, keeping the detector aligned with the beam, such that a large area of the patient's body is scanned by the fan beam of x-rays. Each of the detector segments produces analog signals indicating characteristics of the received x-rays. These analog signals are digitized and fed to a data processing unit which operates on the data in a predetermined fashion to actuate display apparatus to produce a display image representing the internal structure and/or condition of the patient's body. One of the advantages of digital radiography and fluoroscopy is that the digital image information generated from the emergent radiation pattern incident on the detector can be processed, more easily than analog data, in various ways to enhance certain aspects of the image, to make the image more readily intelligible and to display a wider range of anatomical attenuation differences. An important technique for enhancing a digitally represented image is called "energy subtraction". Energy subtraction exploits energy-related differences in attenuation properties of various types of tissue, such as soft tissue and bone, to derive a material-specific image, mapping substantially only a single material in the body. It is known that different tissues, such as soft tissue (which is mostly water) and bone, exhibit different characteristics in their capabilities to attenuate x-radiation of differing energy levels. It is also known that the capability of soft tissue to attentuate x-radiation is less dependent on the x-ray's energy level than is the capability of bone to attenuate x-rays. Soft tissue shows less change in attenuation capability with respect to energy than does bone. This phenomenon enables performance of energy subtraction. In practicing that technique, pulses of x-rays having alternating higher and lower energy levels are directed through the patient's body. When a lower energy pulse is so generated, the detector and associated digital processing unit cooperate to acquire and store a set of digital data representing the image produced in response to the lower energy pulse. A very short time later, when the higher energy pulse is produced, the detector and digital processing unit again similarly cooperate to acquire and store a set of digital information representing the image produced by the higher energy pulse. In early energy subtraction techniques, the values obtained representing the lower energy image were then simply subtracted from the values representing the higher energy image. Since the attenuation of the lower energy x-rays by the soft tissue is about the same as the attenuation of the higher energy x-rays, subtraction of the lower energy image data from the higher energy image data approximately cancels out the information describing the configuration of the soft tissue. When this information has been so cancelled, substantially all that remains in the image is the representation of bone. In this manner, the contrast and visibility of the bone is substantially enhanced by energy subtraction. Energy subtaction has the advantage of being substantially not subject to motion artifacts resulting from the patient's movement between exposures. The time separating the lower and higher energy image acquisitions is quite short, often less than one sixtieth of a second. Details of energy subtraction techniques in digital radiography and fluoroscopy are set forth in the following technical publications Hall, A.L. et al: "Experimental System for Dual Energy Scanned Projection Radiology". Digital Radiography proc. of the SPIE 314: 155-159, 1981; Summer, F.G. et al: "Abdominal Dual Energy Imaging". Digital Radiography proc. SPIE 314: 172-174, 1981; Blink, N. et al: "Dual Energy Radiography: a Preliminary Study". Digital Radiography proc. SPIE 314: 181-182, 1981; and Lehman, L. A. et al: "Generalized Image Ccmbinations in Dual kVp Digital Radiography", Medical Physics 8: 659-667, 1981. The above Article by Lehman et al describes more recently conceived techniques for modifying the above described simple subtraction technique to enhance the quality of the energy subtracted image. Dual energy subtraction has been accomplished, as noted above, by pulsing an x-ray source in a digital scanning slit device at two kVp's, typically 120 and 80 kVp, and sychronizing the pulses with a rotating filter which hardens the high kVp pulses by filtering out the lower energy x-rays. This results in the patient and x-ray detector sequentially seeing high energy and low energy beams from which the mass per unit area of bone and soft tissue can be solved for. More recently, it has been proposed in energy subtraction to utilize a particular type of dual energy detector assembly which can produce separate signals representing each of lower and higher x-ray energy incident on the dectector. Such a detector assembly enables the practice of energy subtraction without the necessity for switching kVp x-ray output levels, or employing other means for periodically attenuating the x-ray beam, such as rapid interposition and removal of a filter to and from the x-ray path. Such a detector employs a dual layer of phosphor-detector elements, wherein the phosphor material of a first, or front, layer preferentially responds to energy of a relatively lower energy value. A second, or rear, detector layer preferentially responds to x-ray energy in a higher range. Such a detector, and its method of use, is described in patent specification EP-A-1 115125. It is often desirable, in medical imaging, to detect characteristics of lesions such as nodules within the human body. Such nodules often are found in the chest area of the patient, including in the lungs. The nodules are typically composed of both soft tissue and calcium, the calcium occurring in differing amounts among different nodules. The ratio of soft tissue to calcium in a lesion such as a nodule can carry highly significant medical implications. For example, a high degree of calcification of a nodule bears a strong correlation with the likelihood that the nodule is not malignant. Calcification, however, can be an indication of disease other than cancer. Moreover, the progression of calcification often carries significant medical implications as well. Additionally, the relative location of calcification within a nodule can have important significance. One of the problems in attempting to quantify calcification of a lesion or other structure or object of interest within the body is that normal anatomy includes much calcified structure, such as bone. Additionally, even soft tissue in the body contains trace amounts of calcium which shows up in energy subtraction detection as a measurable background value. Accordingly, a calcified lesion or other structure, whose calcium content is sought to be detected, is almost always overlying or overlain by other interfering calcium, making it difficult to distinguish between the calcium of the lesion or structure sought to be investigated and the other calcium of the body. A proposal has been made to utilize computerized tomography (CT) equipment for detecting calcification: See Siegelman, S.S., et al., "CT Of The Solitary Pulmonary Nodule", American Journal of Roentgenology, 1980: 135: 1-13; Tarver, R. D. et al., "Experimental Lung Nodule Model: Numbers, Nodule Size, and Actual Calcium Content"; Journal of Computer Assisted Tomography, 1983; 7 (3): 402-406. It is not believed, however, that such methods actually directly determined the quantification of calcification of nodules. Another proposal involved measuring the difference in optical density on conventional radiographic film exposures between a location within a nodule and a location just outside the nodule, which was said to "allow a quantitative estimate of... [calcification] within the nodule." This method, however, yielded only an estimate, and did not directly measure calcification. Also, it was stated to be inapplicable when another calcified structure, such as a rib, overlay the nodule. See Kruger, A. A. et al "Dual Energy Film Subtraction Technique for Detecting Calcification In Solitary Pulmonary Nodules", Radiology, 1981; 140: 213-219. It is believed that the only known method of directly determining calcification of a nodule, in a practical sense, has been the surgical excision and subsequent chemical assay of the nodule. It is an object of this invention to provide a non-invasive method and system for accurately determining the quantity of calcium or other substance present in a lesion or structure of an animal body, and to distinguish the substance in the lesion or structure of interest from other presence within the body of the same substance. According to one aspect of the present invention, there is provided a method of evaluating the amount of a particular substance within an object of interest in an animal body, wherein the substance in said object of interest is at least partially obscured by an additional quantity of such substance, said method utilizing a digital imaging system having a source and detector means for employing penetrative radiation and means, for generating and storing pixel values including signals corresponding to detected radiation emergent from the animal body, said method being characterized by the steps of: a) actuating the source for causing penetrative radiation to pass through both the object of interest and the obscuring quantity of the substance to emerge in a pattern; b) detecting the emergent pattern; c) generating signals defining a set of pixel values related to the emergent pattern; d) identifying first and second portions of said pattern, whereby a first set of pixel values corresponding to said first pattern portion and a second set of pixel values corresponding to said second pattern portion are identified, said first pattern portion including said object of interest and said additional quantity of said obscuring substance and said second pattern portion not including said object of interest but including a like quantity of said obscuring substance to that included in said first portion; e) summing pixel values of said first set; f) summing pixel values of said second set and thereby producing estimated pixel values corresponding to said additional quantity in said first portion and g) subtracting from the sum produced by summing said pixel values of said first set the sum produced by summing said pixel values of said second set to obtain an indication of said amount. According to a further aspect of the invention there is provided a system for quantifying the amount of a particular substance within an object of interest in an animal body, wherein the substance in said object of interest is partially obscured by an additional quantity of such substance, said system utilizing a digital imaging system (s) including means for causing penetrative radiation to pass through both the object of interest and the obscuring quantity of the substance to emerge in a pattern, means for detecting the emergent pattern, means for generating signals defining a set of pixel values related to the emergent pattern, said system being characterized by means for indentifying first and second portions of said pattern whereby a first set of pixel values corresponds to said first pattern portion and a second set of pixel values corresponds to said second pattern portion wherein said first set of pixel values describes a first pattern portion of the subject including said object of interest and said additional quantity of obscuring substance and said second set of pixel values describes a second pattern portion of said subject not including said object of interest but including a like quantity of said obscuring substance to that included in said first pattern portion, means for summing pixel values of said first set, means for summing pixel values of said second set, and thereby producing estimated pixel values corresponding to said additional quantity in said first pattern portion, and means for subtracting from the resultant produced by summing said pixel values of said first set the sum produced by summing said pixel values of said second set to obtain an indication of said amount. The present invention will be understood in greater detail by reference to the following detailed description, and to the drawings, in which: Figure 1 is a pictorial view illustrating an imaging system in which the present invention is incorporated; FIGURE 2 is a plan view corresponding to the system of FIGURE 1; FIGURE 3 is a more detailed block diagram illustrating a portion of the system of FIGURE 1; FIGURE 4 illustrates a portion of anatomy susceptible of examination by use of the present invention; FIGURES 5 and 6 are detail drawings illustrating the application of one embodiment of the present invention; FIGURES 7-9 are detail and graphical drawings illustrating the application of another embodiment of this invention; and FIGURES 10-14 are detail drawings illustrating the application of additional embodiments of the present invention. Figures 1 and 2 illustrate a slit projection type of digital radiography system S in which the present invention is incorporated. The system S scans an X-ray spread beam approximately one millimeters in thickness about a vertical axis across a patient's chest and detects a pattern of X-rays emergent from the patient's body. Information represented by the detected X-rays is processed and displayed to illustrate a representation of an image of the patient's internal body structure or condition. It should be understood at the outset that, while the present invention is described and capable of performance in the environment of digital radiographic imaging systems, the present invention is not in itself an imaging technique. Rather, and as will become clear from the succeeding description, the invention involves method and apparatus for the use and processing of signals which are normally used in imaging to derive and indicate other information relating to quantification of particular materials in certain specified regions of interest within on animal body. More specifically, the system S includes an X-ray source X affixed to mounting structure M for projecting the spread ham B of X-rays through the body of a patient P, to strike an aligned detector assembly D comprising a plurality of detector elements. The spread beam B is confined by a forward slit K to a substantially vertical plane. The detector assembly D comprises a generally vertically elongated array of individual detector elements V and is aligned with the vertical plane defined by the spread ham B. An aft slit J attached to the detector assembly D serves to further define the spread ham B. The X-ray source X is mounted on the structure M to rotate about a vertical axis, defined in Figure 2 as extending into the plane of the paper. Mechanical linkage L couples the X-ray tube X to the detector array D and slits K and J and causes the detector array to scan behind the patient's body along an arcuate path defined by the arrows A, A' in- order to maintain the detector assembly D aligned with the beam B throughout the scanning rotative motion of the beam. The embodiment of the scanning mechanism is not to be limited to fixed or rigid mechanical linkage connecting the elements to be moved. Servo control and associated power drive apparatus embodiments can also be adapted by those of skill in the art to accomplish the desired scanning. In accordance with another aspect of this embodiment, the X-ray tube X can also be pivoted about its focal spot, to maintain the beam B aligned with the scanning detector. The X-ray source X is controlled by power means to emit the spread beam B as either a continuous X-ray beam or a rapid succession of X-ray pulses. The X-ray tube X and the detector assembly D synchronously scan, about a vertical axis, across the patient from one side of his body to the other. Analog detector outputs from each of the detector elements are periodically sampled. Each sampling produces analog signals representing a portion of image information. Over the course of the scan from one side to the other side, signals ore developed describing a plurality of image lines, which together constitute an area image of the patient's internal body structure. The analog signals produced by the detector assembly are provided to an analog to digital converter C which digitizes the outputs and feeds them to a digital processing and receiving unit DPU. The DPU processes these digitized output signals to construct a digital representation of an image of the patient's internal body structure scanned by the X-ray beam B, on a pixel-by-pixel basis. Digital signals from the DPU are converted to analog form by way of a digital to analog converter DAC, and fed to a display unit T, which, in response, produces an image in visual form corresponding to the image representing signals from the DPU. Digital storage means are provided in conjunction with the DPU in order to digitally store the image representations for future use. In such event, the digitally stored signals are later played through the DPU, converted to analog form, and their corresponding images then displayed. The system S produces at its display an image comprising a rectangular matrix array of individual image pixels. Each pixel is defined by a particular brightness level. The pixels are aligned in respective horizontal rows and vertical columns. The radiation imaging system generally described above is explained in the following publication by Tesic, M.M., et al., "Digital Radiography of the Chest: Design Features and Considerations for a Prototype Unit", Radiology, Vol. 148 No.1, pp. 259-264, July, 1983. The detector assembly of the system S comprises detector structure such as described in the above mentioned patent specification EP-A-1 115125. More specifically, the detector comprises two generally vertically extending rows of detectors, one behind the other with respect of the x-ray source. Each of the front and rear arrays are disposed along a substantially arcuate path defining a portion of a circle oriented in a vertical plane. Having discussed the system S in general terms, reference to the block diagram of Figure 3 describes in more detail the operation of the system. A block 10 represents the system gantry and its associated dual energy detector. The x-ray source is actuated, and the gantry and detector are caused to move relative to the subject to acquire raw data representing the high energy image 12 and the lower energy image 14. This raw data is fed to a data processing unit 16, which is a VAX 11/750 Processing Unit manufactured by Digital Equipment Corporation, and 2 DEC RM80 Disc Drives. The disc drives store the image data to be processed and displayed by the rest of the system. The processing unit services operator requests for specified system functions. Certain computations, to be described in more detail below, are performed here, among which are generation of cursor graphics, execution of reprojection for forward projection operations, and background suppression, discussed in more detail below. In a known manner, the processing unit 16 produces signals representing a material-specific soft tissue image, indicated at block 18, and signals representing a material-specific bone image, represented at block 20. The signals representing the respective soft tissue and bone images are transmitted to an image processor 22. The image processor 22 is a Gould DeAnza IP 8400 Image Processor. The image processor 22 manipulates images and graphics for display. Further information on details and operation of the image processor is available in the following publications, each of which is expressly incorporated by reference: Gould DeAnza P8500 Image Processor Software Manual 800219, 1982; and Gould DeAnza IP8500 Image Processor Hardware Reference Manual 800218-REV C, 1983. An output from the image processor 22 feeds a viewing monitor 24 which produces an analog image corresponding to the processed image data which is manipulated by and output from the image processor 22. The data processor 16 also produces signals representing certain summations of image pixel data which will be described in more detail below, and which are input to a printer 26. Also input to the printer 26 ore background suppression signals. The printer 26 produces tangible representations of summation Profiles, indicated at blocks 28, 30. The printer is a DEC LA120 Console/Printer. The processing unit and image processor receive command instructions from a viewing console unit 32. The viewing console 32 includes a keyboard for typed operator input, and a touch-sensitive plasma panel for information display, control knobs for specifying image window, level and combination angles in known fashion, and a trackball for defining cursor graphics locations. The viewing console also includes three image viewing monitors. The viewing console transmits image window and level signal controls in known fashion to the image processor 22. The viewing console also transmits to the image processor signals defining image cursor positioning. Image cursor positioning signals are also transmitted to the processing unit 16. Calcium calibration data are transmitted and processed in known fashion by the processing unit in a known fashion, as indicated at the block 34. The image processor has, in its commercially available form, inherent cursor creating capabilities. A cursor, for the purposes of this disclosure, means a visible overlay on the image displayed on the analog monitor. The viewing console and image processor identified above also include operator trackball and other controls and circuitry for defining the size, shape and orientation of the cursor. The size, shape and orientation of the cursor are communicated to the processing unit, which includes means for separately processing and handling image representing signals corresponding to the pixels of the image encompassed by the particular cursor selected by the operator. The use of the cursor capability will be discussed in more detail below. One of the functions of which the processing unit is capable is that of separately summing the brightness values for selected image pixels encompassed by a cursor, either in their entirety, or separately along each of a plurality of separate lines of pixels within the cursor. A known technique for performing this funtion is known as "reprojection", and is described in the following Publication, which is hereby expressly incorporated by reference: Joseph, P.M. "An Improved Algorithm for Reprojecting Rays Through Pixel Images", IEEE Transactions On Medical Imaging, Vol. MI-1, No. 3, November, 1982; on alternative technique is described in Newman, W.M., et al., "Principles of Interactive Graphics", 1979, McGraw-Hill, pp. 22-25. Radiography of the type described above is often employed to search for lesions, such as nodules or other objects of interest in the chest region of a patient to be examined. Such nodules sometimes are benign, but sometimes constitute malignant tumors which must be treated or excised. Such nodules are often calcified to a degree which may vary from one nodule to another: The remainder of the nodule which does not consist of calcium has most of the attributes of soft tissue. Additionally, bone contains much calcium which is prominently imaged, and even uncalcified soft tissue contains trace amounts of calcium which con show up in radiographic images as measurable background. Dual energy radiography produces material-specific images of an animal body, and can be set up to map substantially only calcium in the body. Besides making certain pathological conditions easier to recognize, data describing images derived by energy subtraction techniques are also susceptible of processing to provide quantitative information about the water/calcium make-up of structures or objects of interest such a lung nodules. Such quantification determining techniques provide, essentially, a non-invasive chemical assay of the calcium content of a nodule, other lesion, or of normal anatomy. Consider an image of bone derived by the use of digital radiography employing the energy subtraction technique. Each pixel of the image defines a ray path from the x-ray tube focal spot through the patient to the detector surface. The value of each pixel is proportional to the amount of calcium along its ray path. In medical physics, it is often convenient to describe the amount (thickness) of a substance being imaged in terms of mass per unit area, usually given as grams per square centimeter (g/cm²). Pixel values in the bone image are proportional to the mass per unit area of overlying calcium. This technique is described in the following publication which is hereby expressly incorporated by reference: Johns, H.E., et al., "The physics, of Radiology", 4th Edition, 1983, Charles C. Thomas-Publisher, pg. 140. The present invention provides a technique for extracting quantitative data from images of calcium containing nodules, lesions or other body structure as defined by a region of interest. More specifically, the technique is designed to measure the total calcium mass of lung nodules. Consider imaging a lung nodule. In many cases, the image of the nodule will overlie or be overlain by one or more obscuring interfering objects, such as ribs. For the moment, however, consider the simple case where there is no rib superposition. The pixel values in the bone image are proportional to overlying or interfering g/cm² of calcium: W = KP (1) where P is a bone image pixel value, W is the g/cm² of calcium overlying that pixel, and K is a constant which is characteristic of the particular imaging system used. Now let Ap be the area (at the image surface) of a pixel, and let M be the nodule magnification. Then the calcium mass Bj overlying pixel j is given by the relation: B j = (A p /M²) W j = (KA p /M²)P j (2) Assume also, for the moment, that there is no calcium "background" overlying or aligned with the nodule. Then P equals 0 for pixels outside the nodule. Now define a region R which includes the nodule. The total calcium mass B of the nodule is then given by: Equations (2) and (3) indicate how to calculate the calcium mass of a nodule, assuming: (1) the nodule does not overlie any ribs, and (2) the background overlying and surrounding the nodule has no calcium component. The two assumptions required, however, to derive equations (2) and (3) are seldom true in practice. The soft tissue background overlying and surrounding a nodule does have a small calcium component. Moreover, tests have suggested that between 50% and 75% of all nodules overlie at least one rib. "The term "overlie" as used here means that she shadowgraphic image of the nodule or object of interest is at least partially overlapping with the shadowgraphic image of the other, or interfering, structure, such as a rib as referred to here. This invention includes a technique for quantifying the amount of calcium in a nodule, lesion or other body structure of interest, while distinguishing the calcium of the structure of interest sought to be evaluated from the calcium background and from the calcium of any overlying ribs or other interfering or obscuring calcium. This technique is illustrated, for example, in Figures 4-6. Figure 4 illustrates a portion 100 of a rib of a human patient to be examined, and a nodule 102 Partially overlying the rib. As a point of reference, x-rays produced by the system described above are considered as passing through the rib and nodule in a direction perpendicular to the plane of the paper. The detector assembly is not shown. One method of quantifying the calcium in the nodule 102 is illustrated in Figures 5 and 6. Figures 5 and 6 illustrate the rib 100, the nodule 102, and a rectangular box 106. The box 106 is to be considered as a cursor appearing on the display screen of the main imaging system described above at a location corresponding to its location illustrated in Figures 5 and 6, respectively. The cursors size, orientation, and position are controllable by an operator in known fashion. The purpose of the cursor is to define a region of interest around the nodule. The main imaging system, as noted above, is in known fashion capable of separately processing, enhancing and/or handling data from image pixels within the boundaries of the regions defined by the cursors. In Figures 5 and 6, the cursors are rectangular, having their side pairs aligned in horizontal and vertical directions. The cursors can, however, be of any mutually congruent geometry. The system S is capable of establishing multiple cursors simultaneously, and the cursors illustrated in 5 and 6 are established in that manner. One method of evaluating the amount of calcium in the nodule 102 is to compare two sets of data, i.e., the data acquired in the region defined by the cursor 106 at the Figure 5 location, and one set acquired in the region defined by the cursor located as in Figure 6. This will be described in more detail below. The cursor of Figure 5 is established to encompass all of the nodule 102, a portion of the rib, and other background. The cursor of Figure 6 is established to encompass a portion of the rib, and other background, and does not encompass any portion of the nodule 102. The rib portion 108 defined by the cursor positioned as in Figure 6 is substantially the same size, shape and location relative to the rib edge as the region 109 of the rib encompassed by the cursor of Figure 5. Thus, the cursor depicted in Figure 5 defines a region whose image pixels describe the contribution to radiation attenuation brought about by the relatively uniform soft tissue calcium traced background, by the rib, and by the nodule. The region defined by the cursor as positioned in Figure 6 exhibits radiation attenuation effects brought about by only the relatively uniform background, and by the rib. Moveover, the portion of the rib included within the cursor of Figure 6 is substantially congruent, and analogously positioned, as is the rib portion encompassed by the cursor of Figure 5. The data acquired from the pixels defined by the two cursors simultaneously positioned as in Figures 5 and 6, respectively, are separately stored. The pixel values of the pixels encompassed by the cursor of Figure 5 are summed. The pixel values encompassed by the cursor positioned as shown in Figure 6 are separately summed. The sum associated with the cursor of Figure 6 is subtracted from the sum associated with the cursor of Figure 5. Where the data from the regions of the two cursors is derived from information describing the calcium image of the subject, as in energy subtraction, the remainder, a scalar number, is a function of the quantity of calcium in the nodule. The subtraction operation effectively cancels out the contribution of calcium located in the rib, and calcium located in soft tissue in the region. While this embodiment of the invention is primarily described for use in the context of dual energy subtraction, it is contemplated that useful information can be derived by the implementation of this embodiment in the context of single energy exposures. Where only single energy exposure is used, however, the information derived will be a function of the total mass of the nodule, rather than a quantification of its calcium content alone. Another embodiment of a method for assaying the amount of calcium in a nodule is to derive a "local" background from a region defined by a stationary single cursor. Such a technique is illustrated in connection with Figures 7-9. Figure 7 shows a nodule 119 partially overlying a rib 110, and enclosed within a rectangular cursor box 116. As In connection with the previous embodiment, this box is to be thought of as a cursor on a display. The cursor's size, orientation and position are controlled by an operator. The purpose of the cursor is to define a region around the nodule. In this embodiment of the invention, two assumptions are made about the distribution of calcium in the rib: (1) the rib edge 112 is approximately a straight line over the region defined by the cursor and (2) the rib calcium distribution varies substantially only in the direction perpendicular to the rib edge. It is also assumed that the background soft tissue calcium is essentially constant over the region defined by the cursor. The cursor rectangle is oriented with one pair of sides parallel to the rib edge. Within the cursor, pixel sums are generated along lines perpendicular to the rib edge, yielding a summation profile of the cursor region. Such a profile is illustrated in Figure 7, where an elevated portion, or "hump" 118 evidences the calcium in the nodule. The profile also includes "tail" portions 120, 122. Preferably, the image pixel values are generated in the environment of execution of the energy subtraction technique discussed above, such that the resultant summation profile is based on substantially only data representing the content of calcium in the structures encompassed by the cursor. This technique could, however, he used with single energy radiography, but the profile would then depict the attenuation characteristics exhibited by all the materials within the cursor-defined region, rather than on the particular attenuation characteristics of the calcium alone. The hump portion of the profile corresponds to the region in which pixel lines perpendicular to the rib edge intersect the nodule, while the tail portions of the profile correspond to regions within the cursor wherein lines of pixels perpendicular to the rib edge are aligned only with a portion of the rib itself and with background material, and do not intersect the nodule. As an aid in understanding how this technique works, consider Figure 8 which is identical to Figure 7 except that the nodule is absent. Thus, the profile in Figure 8 is derived from only the rib configuration and the uniform calcium background contribution, and is thus approximately uniform along a direction parallel to the rib edge. This uniformity results from the rib edge being substantially a straight line and the rib having substantially uniform calcium distribution along a direction parallel to its edge. In order to obtain a representation which corresponds substantially only to the calcium in the nodule, one simply overlays the profiles of Figures 7 and 8 and subtracts away the contribution to the profile of Figure 7 which results from the rib/background contribution, leaving only the nodule calcium profile. The nodule calcium profile can then be integrated in a known manner to yield a value for total calcium content. It is important to realize that the procedure depicted in Figure 8 is not actually separately performed distinct from the procedure depicted in connection with Figure 7. The Figure 8 procedure is hypothetical only. In practice, the rib profile (the uniform profile shown in Figure 8) can often be estimated with sufficient accuracy in the profile of Figure 7 by simply connecting the tails of the Figure 7 Profile produced by the combination of the rib/background/nodule combination. This procedure is illustrated in Figure 9, where the shaded area corresponds to the integral of the nodule calcium profile. Where the two "tail" portions of the summation profile are approximately horizontal and of equal value, interpolation of the value of the rib/background portion (Figure 8) across the area of the nodule (Figure 7) comprises only a simple step of connecting the two tail portions over the space between them. Where, however, the respective tail portions are of different magnitudes, or are non-uniform along their visible lengths, known interpolation or extrapolation techniques can be employed to obtain a statistically valid estimate or approximation of the rib/background contribution of the Figure 8 profile as it affects the nodule encompassing portion of the total profile as shown in Figure 7. Where S represents the integral of the nodule calcium profile, S is analogous to the pixel station of equation (3), so the total calcium of the nodule is given by: B = (KA p /M²) x S (4) In developing the sums of pixel lines extending perpendicular to the rib edge, it must be kept in mind that, (in the Figure 7 embodiment) unless the rib edge happens to be aligned with the actual pixel rows or columns in the image produced by the system, means must be employed to approximate the pixel line integrals. Orientation of the cursor at a non-parallel angle with respect to the pixel lines or columns necessitates employment of a known technique called "reprojection". Reprojection sums the pixels encompassed within the non-parallel orientated cursor in a way which causes approximations of pixel lines within the cursor. The pixel sums thus derived are really approximations to line integrals, and so the technique amounts to finding a method for estimating an arbitrarily directed line integral. This problem is generally referred to as reprojection or forward projection because it is, in a sense, the reverse operation of back projection, a known technique often used in computerized tomography. Reprojection techniques suitable for incorporation in the practice of this invention are discussed in the above incorporated articles by Joseph and by Newman. Another embodiment of the present invention makes it possible to obtain the advantages of the previously described embodiment without the necessity for employing complex, general reprojection techniques in generating pixel line integrals. In accordance with this embodiment, and with reference to Figure 10, a rib 130 having an edge 132, and a nodule 134, are shown. In the embodiment of Figure 10, however, the cursor is defined as a parallelogram 136. As in the case of Figures 7 and 8, the rib edge 132 is oriented at an angle in the image which is neither vertical nor horizontal, and is thus not parallel to either the columns or rows of pixels. The parallelogram shaped cursor 136 has no right angles. The parallelogram 136 has one of its side pairs oriented substantially parallel to the rib edge 132, while the other side pair is oriented vertically. This means that, if pixel lines parallel to the vertical side pair ore summed across the parallelogram 136, it is not necessary to employ any reprojection techniques, since the summed lines of pixels are parallel to the pixel lines defined in the image detection apparatus. In Figure 10, a profile 138 is illustrated, which is developed in a manner analogous to that described in connection with Figure 7. As in the embodiment described in connection with Figures 7 and 8, the tail portions of the profile 138 can be connected, and the area of the hump above the connected tail portions represents the total calcium content of the nodule 134. The embodiment described in connection with Figure 10 is employed when the slope of the line 132 defined by the rib edge is less than or equal to unity, i.e., when the angle of the line 132 to the horizontal is less than or equal to 45°. Another embodiment is employed when the rib edge defines a line whose slope, with respect to the horizontal,is greater than or equal to unity, i.e., is greater than or equal to 45°. Figure 11 illustrates this embodiment. Figure 11 illustrates a rib 140 having an edge 142 approximately defining a line whose angle to the horizontal is greater than 45°. A nodule 144 is also illustrated. A cursor 146, having the configuration of a parallelogram without right angles is also illustrated. In this embodiment, one of the line pairs defining the parallelogram are substantially parallel to the line 142. The opposite pair of sides are horizontal. In this embodiment, the pixel lines encompassed within the cursor 146 are summed in a horizontal direction. This means that there is in the Figure 11 embodiment also no need to utilize complex reprojection techniques to sum pixel lines, the summation taking place horizontally in the image, along the direction of the pixel rows inherently defined in the display and by the system circuitry. A profile 148 is also shown, representing the result of the summations taking place. The profile 148 is employed in a manner analogous to that discussed in connection with Figures 7 and 8. Thus, utilizing the embodiments of Figures 10 and 11, two of the parallelogram lines defined by the cursor configuration are substantially parallel to the rib edge, while the remaining two lines are parallel to either the pixel rows or pixel columns. In these embodiments, reprojection reduces to the trivial operation of row-wise or column-wise summation. The parallelogram configuration for the cursor affords more positional flexibility than does a rectangle, allowing one to quantify calcium in nodules or other structures that are in difficult locations. The relevant operating principles for constructing the cursor are that: (1) where an interfering or obscuring object defines a generally linear edge, one line pair of the parallelogram should be parallel to that edge, and (2) lines of pixel summations should extend parallel to the other line pair, irrespective of the direction of that other line pair. No other geometrical condition is essential. It is not required that other line pair be parallel to pixel rows or columns, though, as discussed above, it Bay sometimes be desirable to so construct the cursor. Figure 12 shows an instance where a rectangular cursor 150 overlies two ribs 152, 154 even though a nodule 156 over lies only one rib 152. Even if the nodule were not present, the background profile, to which only the rib and background calcium contribute (analogous to that of Figure 8) would not be a straight line. Thus, estimation of that profile from the tails of the profile analogous of that to Figure 7, to which all of the rib, background and nodule contribute, would be inaccurate, making qantification of calcium in the nodule correspondingly inaccurate. Figure 13 indicates how the generalized parallelogram form of cursor affords the flexibility, in a situation such as Figure 12, of structuring the cursor such that it encompasses a portion of only one rib. Calcium quantification accuracy in this instance is maintained. As often happens in practice, Figure 13 shows one side of the parallelogram skirting very closely to the nodule. Because of this, substantially no "tail" portion of the profile will appear corresponding to this side of the nodule. The tail from the other side of the nodule must in such instance be used (e.g. by extrapolation) to estimate the profile analogous of Figure 8 to which only the rib and background contribute, so that it may be subtracted away. This one tail extrapolation method has been found in tests to give good results, although it is considered preferable where available to use the two tail technique in estimating the background contribution. Sometimes a nodule will overlie more than one rib. Quantification even in this instance is still possible, but requires that the nodule he considered as broken-up or segmented into several contiguous regions. Refer to Figure 14. The calcium content of each region of the nodule is calculated as described above, and the results are added to obtain the total calcium content of the nodule. There is a subtlety involved in the segmentation technique of Figure 14. A region 1 includes part of the nodule and part of the upper rib. Quantification of this region poses no problem, being done in a manner analogous to that discussed in connection with Figure 13. However, a region 2, as described in Figure 14, includes a subregion where the upper and lower ribs overlap. The sides (bases) AF and BE of the region 2 parollelogram ore parallel to the edge of the lower rib, but not parallel to the edge of the upper rib. Hence, the portion of the profile contributed by only the background and the rib of region 2 is not flat, and the model upon which the quantification technique is based is not strictly correct. Tests have shown, however, that the segmentation technique depicted in Figure 14 works quite well when estimation of the rib/background contribution is made utilizing two "tails" of each profile. This is believed to be the case because estimation using two tails averages the tails on either side of the nodule, giving an estimated mean for the rib/background profile which is acceptably proximate the actual mean. Generation of cursor graphics will now be discussed. A parallelogram is defined when three of its vertices are defined. An operator can define a parallelogram using the trackball of the viewing console and a point cursor (that is, a small, essentially single pixel cursor). The point cursor is moved around by using the trackball and when it is in a location judged by the operator to be appropriate for a vertex of the parallelogram, the operator in known fashion actuates the viewing console to "write" a dot (vertex) into the display overlay plane. Then the operator can move the point cursor to the next place where a vertex is desired, etc. The operator first defines two such points which define a line segment parallel to a rib edge, as described in more detail below. These points (pixels) are appropriately called P₁ and P₂. Once P₁ and P₂ are defined, the processing unit takes their x-y coordinates and calculates the best (pixel) line connecting these points, in a manner analogous to that disclosed in the Joseph and Newman publications incorporated herein. Techniques for taking the coordinates of two points and generating a best pixel line connecting them are well known, as set forth in the above referenced articles. After the points on this line are calculated, the line is written into the display overlay plane. Next, the operator defines a third point in such a way that the ensuing parallelogram encloses the region of interest appearing in the analog image. Once the third point P₃ is defined, it is a simple matter for the processing unit to calculate the fourth and final point P₄. The calculation of the location of P₄ is accomplished in the following way: let X₁ and Y₁ be the coordinates of P₁, etc. Then, X₄ = (X₃ - X₂) + X₁ Y₄ = (Y₃ - Y₂) + Y₁ Once P₁-P₄ are known, the processing unit uses a known line drawing technique to write the remaining three aides of the parallelogram into the overlay plane. Let us call the base, and the side parallel to the base, (both parallel to the rib edge) both "bases" of the parallelogram, and call the other two sides "sides". Also, call the "length" of the parallelogram the length of (number of pixels in) each of the bases, and the "width" of the Parallelogram the length of (number of pixels in) each of the sides. This nomenclature will be convenient in the discussion below. After the parallelogram is completely defined and displayed, the processing unit then operates to acquire the summation profile in accordance with known technique. Pixel summations are done parallel to the sides of the parallelogram. The technique for profile generation is now discussed. One begins with the pixel line (parallelogram side) defined between P₁ and P₄. One then adds up all the values of the pixels in this line. This represents the first "projection" (point, value, or element) of the profile. Next, one moves to the next pair of points of the respective parallelogram bases, using the line defining technique described above to define a pixel line between these points. One then adds up all the values of the pixels in this new line. This is the second "projection" of the profile. One then continues moving to successive pairs of respective base points until the other side of the parallelogram is reached, i.e., points P₂, P₃. The profile is generated in this way. This is the profile which is then "tail-subtracted" and "integrated", as described in detail in connection with Figures 7 et seq. Where integration of the area under a curve is required, the integration simply amounts to summing up the profile projections after subtracting the tail background. It will be seen, in connection with the above approach to profile generation, wherein pixel lines are approximated by reprojection, that see pixels are used more than once during the summations, and other pixels are not used at all. That is, the total number of pixels used in generating the profile is not necessarily equal to the total number of pixels contained in the parallelogram, due to the line approximation technique which is necessary, and which is discussed in connection with the incorporated publications by Joseph and Newman. Ideally, if one took the profile and summed all the values in the profile without tail background correction, one would get the same answer as summing all the pixels in the parallelogram. Since this is not true in practice, one must apply a correction factor to the result of profile integration. The correction factor is simply the ratio of the actual number of pixels in the parallelogram, including those representing the bases and sides, to the number of pixels used in generating the profile. The number of pixels N used in generating the profile is simply the length times the width of the parallelogram. The actual number of pixels in the parallelogram (A) is calculated from a formula for the area of a parallelogram: A = (X₁ Y₂) + (X₂ Y₃) + (X₃ X₁) - (X₃ Y₂) - (X₂ Y₁) - (X₁ Y₃). Therefore, equation (4) of the disclosure becomes: B = (A/N) (K A p /M²) S (5) Several of the above disclosed equations incorporate a constant K which is derived from, and which is a characteristic of, the particular attributes of the system being used. K is conveniently expressed in units of milligrams per square centimeter (mg/cm²) per pixel count. A step-wedge phantom containing steps of various known thicknesses (mg/cm²) of calcium is imaged, and these thickness are plotted vs. the bone image pixel values that the thicknesses of calcium yield. The plot is a straight line with slope K.

Claims: A method of evaluating the amount of a particular substance within an object of interest (102, 119, 134, 144, 156) in an animal body (P), wherein the substance in said object of interest (102, 119, 134, 144, 156) is at least partially obscured by an additional quantity (100, 110, 130, 152, 154) of such substance, said method utilizing a digital imaging system (S) having a source (X) and detector (D) means for employing penetrative radiation and means (C, 16) for generating and storing pixel values including signals corresponding to detect radiation emergent from the animal body (P), said method being characterized by the steps of: a) actuating the source (X) for causing penetrative radiation to pass through both the object of interest and the obscuring quantity of the substance to emerge in a pattern; b) detecting the emergent pattern; c) generating signals defining a set of pixel values related to the emergent pattern; d) identifying first and second portions (106, 118, 120, 122) of said pattern, whereby a first set of pixel values corresponding to said first pattern portion and a second set of pixel values corresponding to said second pattern portion are identified, said first pattern portion (106,118) including said object of interest (102) and said additional quantity (109) of said obscuring substance and said second pattern portion (106, 120, 122) not including said object of interest (102) but including a like quantity (108) of said obscuring substance to that included in said first portion; e) summing pixel values of said first set; f) summing pixel values of said second set and thereby producing estimated pixel values corresponding to said additional quantity (109, 110, 130) in said first portion (106, 118); and g) subtracting from the sum produced by summing said pixel values of said first set the sum produced by summing said pixel values of said second set to obtain an indication of said amount. A method as claimed in Claim 1 wherein said first and second pattern portions (106) are mutually substantially congruent. A method as claimed in Claim 1 wherein said pattern portions are defined within a single parallelogram (116, 136, 146, 150) encompassing said object of interest (119, 134, 144, 156) and said additional obscuring quantity (110, 130, 140, 152, 154), wherein said summing steps comprise separately summing pixel values corresponding to respective substantially parallel lines of pixels within said parallelogram and further comprising producing a profile (138, 148) of said summations, said profile defining said first pattern portion (118) corresponding to pixel lines intersecting said object of interest (119, 134, 144, 156) and said second pattern portion (120, 122) corresponding to pixel lines not intersecting said object of interest (119, 134, 144, 156) utilizing said second pattern portion (120, 122) to estimate the contribution to said profile of said obscuring quantity, producing a representation of said summation profile less said estimated value of said contribution of said obscuring quantity and producing a representation of the integral of the remainder of said summation profile. A method according to Claim 3 wherein said summing steps are conducted along lines substantially parallel to said sides of said parallelogram (116, 136, 146, 150). A method according to Claim 3 or Claim 4 wherein said additional obscuring quantity of said substance resides in an interfering object (110, 130, 140) defining an edge (112, 132, 142) which approximates a straight line, said pixel summing steps comprising: summing said pixels along respective lines in a direction substantially perpendicular to said edge (112, 132, 142). A method according to Claim 3 or Claim 4 wherein said additional obscuring quantity of said substance resides in an interfering object (110, 130, 140) defining a substantially straight line (112, 132, 142), and wherein said parallelogram (116, 136, 146) is defined with sides substantially parallel to said edge. A method according to any one of Claims 3 to 6 wherein said system (S) comprises apparatus and circuitry for operating in an energy subtraction mode. A method according to any preceding claim wherein said substance comprises calcium. A method according to Claim 8 wherein a) said object of interest is within a lung, and b) said obscuring substance resides in a rib (100, 110, 130, 140, 152, 154). A system for quantifying the amount of a particular substance within an object of interest (102, 119, 134, 144, 156) in an animal body (P), wherein the substance in said object of interest is partially obscured by an additional quantity (100, 110, 130, 152, 154) of such substance, said system utilizing a digital imaging system (S) including means (X) for causing penetrative radiation to pass through both the object of interest and the obscuring quantity of the substance to emerge in a pattern, means (D) for detecting the emergent pattern, means (C) for generating signals defining a set of pixel values related to the emergent pattern, said system being characterized by means for identifying first and second portions (106,118, 120, 122) of said pattern whereby a first set of pixel values corresponds to said first pattern portion and a second set of pixel values corresponds to said second pattern portion wherein said first set of pixel values describes a first pattern portion (106, 118) of the subject including said object of interest (102) and said additional quantity (109) of obscuring substance and said second set of pixel values describes a second pattern portion (106, 120, 122) of said subject not including said object of interest (102) but including a like quantity (108) of said obscuring substance to that included in said first pattern portion (106, 118), means for summing pixel values of said first set, means for summing pixel values of said second set, and thereby producing estimated pixel values corresponding to said additional quantity (109) in said first pattern portion (106), and means for subtracting from the resultant produced by summing said pixel values of said first set the sum produced by summing said pixel values of said second set to obtain an indication of said amount. A system as claimed in Claim 10 wherein said pattern portions are defined within a single parallelogram (116, 136, 146, 150) encompassing said object of interest (119, 134, 144, 156) and said additional obscuring quantity (110, 130, 140, 152, 154), wherein said summing means include means for separately summing pixel values corresponding to respective substantially parallel lines of said pixels within said parallelogram; and said system further includes means for producing a profile (138, 148) of said summations, said profile defining said first pattern portion (118) corresponding to pixel lines intersecting said object of interest and said second pattern portion (120, 122) corresponding to pixel lines not intersecting the object of interest; means for utilizing said second pattern portion (120, 122) of said profile to estimate the contribution to said additional quantity of said substance located outside said object of interest (119, 134, 144, 156); means for producing a representation of said summation profile less said estimated value of said contribution to the profile of said additional quantity of said substance, and means for producing a representation of the integral of the remainder of said summation profile.

Inventors: BARNES, GARY T., LAURO, KAREN L., SONES, RICHARD A., TESIC, MIKE M.
Applicants: PICKER INTERNATIONAL, INC.
Application year: 1986
Main classification: G06F15
IPC codes: G01N23, A61B6, G06T7