US4956859A - Source filter for radiographic imaging - Google Patents
Source filter for radiographic imaging Download PDFInfo
- Publication number
- US4956859A US4956859A US07/321,989 US32198989A US4956859A US 4956859 A US4956859 A US 4956859A US 32198989 A US32198989 A US 32198989A US 4956859 A US4956859 A US 4956859A
- Authority
- US
- United States
- Prior art keywords
- photons
- source
- radiation
- filter element
- filter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/50—Clinical applications
- A61B6/505—Clinical applications involving diagnosis of bone
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/04—Positioning of patients; Tiltable beds or the like
- A61B6/0407—Supports, e.g. tables or beds, for the body or parts of the body
- A61B6/0421—Supports, e.g. tables or beds, for the body or parts of the body with immobilising means
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21K—TECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
- G21K1/00—Arrangements for handling particles or ionising radiation, e.g. focusing or moderating
- G21K1/10—Scattering devices; Absorbing devices; Ionising radiation filters
Definitions
- the technical field of this invention is radiology and, in particular, bone absorptiometry by radiographic measurements.
- osteoporosis The depletion of bone mineral content, typically referred to as osteoporosis, is a common consequence of a variety of diseases and natural aging process.
- bone minerals can be lost as the result of drugs, stress, dietary deficiencies, pregnancy or lactation.
- skeletal bone mass drops below the level necessary to provide mechanical support, the depletion of bone mineral content becomes an important cause of morbidity, particularly in elderly patients.
- BMC bone mineral content
- a system for analyzing biological structures by photon absorptiometry including a radiation source emitting photons and at least one source filter which operates to reduce the emission of high energy photons, thereby creating a sharp drop-off at the high end of the detected energy spectrum.
- the filter radiation is received by a detector which determines the spatial intensity of the radiation over the image area and preferably has a characteristic response which is relatively insensitive to low energy photons, such that the source filter and detector cooperate to measure the intensity of photons within a narrowed band of energy levels.
- the radiation source is a broad band X-ray source emitting X-rays at energy levels up to about ⁇ KeV and the source filter creates a sharp drop-off of photons having energy levels above about 50 Kev.
- the detector can be a multiple wire proportional chamber (MWPC) device filled with an ionizable gas, or alternatively, a scintillation screen or a photosimulatable luminescence measuring device.
- MWPC device which is insensitive to most photons having energies less than 35 KeV.
- the system can also include one or more additional source filter elements to further narrow the energy band of the photons.
- one or more metal filter elements ranging from 0.01 to 1.0 millimeters in thickness, are employed to sharpen the energy band and provide for highest transmissivity in the range of about 42 to about 48 KeV.
- a composite K-edge filter using 2 or more, preferably 3, pure metals stacked on one another is particularly useful in modifying the raw spectrum of a radiation source to achieve greater spatial resolution in the detector.
- metals with atomic numbers from 45 to 52 usually two metals
- at least one metal whose atomic number ranges from 57 to 65 give an optimum spatial resolution in Xe gas detectors.
- cadmium (Cd), silver (Ag), and gadolinium (Gd) can be particularly useful in combination.
- cadmium (Cd), tin (Sn) and samarium (Sm) form another useful combination filter.
- the result of using the composite filter is twofold: first, by blocking low energy photons, it reduces the radiation dose otherwise absorbed by the patient and, second, by narrowing the energy band of the transmitted photons, it dramatically diminished the projected electron range in the detector, itself, thereby improving spatial resolution.
- Bone mineral content measurements typically involve measurement of the intensity of photons through a tissue equivalent bolus which is the same thickness as patient's limb and which contains a material having an absorption characteristic which is the same as tissue free of fat.
- the fat content can be a significant (e.g., five to fifteen percent) absorptiometry factor in heavy individuals; and when it varies, it will lead to errors in measurement of BMC changes from year to year.
- the composite filter removes a large percentage of the photons in the range 20 to 35 Kev where the difference between the absorption coefficients for fat and tissue is largest, the observed differences in photon intensity through fat-containing tissue will be reduced substantially for the filtered x-ray spectra.
- FIG. 1 is an isometric view of a radiographic imaging system according to the present invention
- FIG. 2 is a schematic block diagram of a radiographic imaging system according to the invention.
- FIG. 3 is a more detailed schematic diagram of the radiation source and detection of the system of FIG. 2;
- FIG. 4 is a more detailed schematic illustration of the beam collimator of FIG. 3;
- FIG. 5 is an exploded schematic illustration of the beam filter of FIG. 3;
- FIG. 6 is a schematic side view of a multi-wire detector module useful in the system of FIG. 2;
- FIG. 7 is a more detailed schematic illustration of the cathode and anode grids of the detector of FIG. 6;
- FIG. 8 is a graph showing a typical raw energy spectrum of an X-ray source, such as shown in FIG. 3;
- FIG. 9 is a graph showing a typical energy transmission spectrum of a beam filter, such as shown in FIG. 3;
- FIG. 10 is a graph showing a typical filtered beam spectrum after passage through a beam filter, such as shown in FIG. 3;
- FIG. 11 is a graph showing typical detection efficiencies versus beam energy for a detector, such as shown in FIG. 3;
- FIG. 12 is a graph showing the overall energy response of a radiographic imaging system according to the invention.
- FIG. 13 is a more detailed schematic top view of a limb positioning and calibration device of the system of FIG. 2;
- FIG. 14 is a schematic side view of the positioning and calibration device of FIG. 13 in a closed position
- FIG. 15 is a schematic side view of the Positioning and calibration device of FIG. 13 in a closed position
- FIG. 16 is a schematic block diagram of the data read-out circuitry of FIG. 6;
- FIG. 17 is a more detailed diagram of a gated event digitizer as shown in FIG. 16;
- FIG. 18 is a graphic illustration of the priority encoding operation performed by the circuitry of FIG. 16 to obtain an event coordinate;
- FIG. 19 is a schematic block diagram of data acquisition circuitry useful in the system of FIG. 2;
- FIG. 20A is a timing diagram for a programmable logic device as shown in FIG. 19 and further illustrated schematically in FIG. 20B;
- FIG. 21 is a schematic illustration of a data flow path within the circuitry of FIG. 19;
- FIG. 22 is a schematic illustration of a patient X-ray image (showing non-imaged bones in phantom) in accordance with the present invention.
- FIG. 23 is a more detailed view of an X-ray image illustrating a method of identifying bone reference points according to the present invention.
- FIG. 24 is another detailed view of an X-ray image illustrating a method of defining bone analysis regions according to the present invention.
- FIG. 1 is an isometric view of a radiographic imaging system 10 for analyzing biological structures according to the present invention consisting of a radiation source 12, which transmits a beam of radiation to detector module 14. Disposed between the radiation source 12 and the detector module 14 is the limb of a patient 16, which is secured in place by a limb positioning apparatus 18 and registration means 15A, 15B which prevents the limb from pivoting about the positioning apparatus 18 in operation.
- the detector module 14 and the limb positioning apparatus 18 and the registration means 15A, 15B conveniently can be incorporated into a housing 20 so that the patient can sit while the procedure is being conducted.
- FIG. 2 is a more detailed schematic block diagram of the radiographic imaging system as shown, including radiation source 12, detector module 14, positioning apparatus 18, data acquisition module 22, data processing module 24, data display module 26 and data storage module 28.
- Radiation emitted by the radiation source 12 passes through the patient's limb (e.g., the wrist) 16 and is detected by detector module 14 and converted therein to digital data.
- the data is then transferred to the acquisition module 22 where a radiographic image of the patient's limb is built up over time.
- the images formed in the acquisition module can be analyzed in the data processing module 24 (e.g., to detect osteoporosis or other changes in bone mineral content), displayed by the display module 26 and/or stored for subsequent analysis or display in data storage module 28.
- FIGS. 3-7 are more detailed schematic illustrations of the radiation source 12 and detector module 14.
- radiation source 12 includes a radiation head assembly 30 which houses a radiation emitter 32, such as a radioisotope or, preferably, an X-ray tube operable to generate photon radiation having maximum energy levels, for example, on the order of about 45 to about 70 KeV and, preferably, about 55 KeV.
- a radiation emitter 32 such as a radioisotope or, preferably, an X-ray tube operable to generate photon radiation having maximum energy levels, for example, on the order of about 45 to about 70 KeV and, preferably, about 55 KeV.
- X-ray tubes are available commercially from a variety of sources.
- One such commercially available tube is the Model OIX-15 dental X-ray tube (Eureka Company, Chicago, Ill.).
- a detector module 14 which preferably is an area detector capable of determining the spatial intensity of radiation emitted by the radiation source 12 and transmitted through the patient's limb 16.
- One such area detector is a multi-wire, proportional chamber detector, described in more detail below in connection with FIGS. 6 and 7.
- Other area detectors which can be employed in the present invention include scintillation screen devices, employing gadolinium sulfate or lanthanum sulfate crystals and associated electronic imaging (i.e., CCD) elements, or photosimulatable luminescence (PSL) measuring devices. Detectors of these types are commercially available, for example, from Kodak (Rochester, N.Y.) or Fuji (Tokyo, Japan).
- the radiation head assembly 30 shown in FIG. 3 also includes a threaded barrel 33 into which can be inserted a collimator assembly 34 and a filter assembly 38.
- a collimator assembly 34 is shown consisting of a machined metal (i.e., lead) body 37 having a central hole 35 and peripheral threads 36 so that it can be screwed into the threaded barrel 33 of the radiation head assembly.
- FIG. 5 is an exploded schematic illustration of the beam filter assembly 38, including a filter housing 40 and cover plate 42 which can be connected by bolts 44 on the housing through holes 46 on the cover plate via nuts or other suitable attachment means. Disposed within the filter assembly is one or more filter elements. In the illustration, three such filter elements 48, 50 and 52 are shown. In one preferred embodiment of the invention, three filter elements of gadolinium, silver and tin, each having a thickness on the order of about 0.08 millimeters can be employed.
- At least one filter element which operates to reduce the transmission of high energy photons in conjunction with a detector module having a characteristic response which is relatively insensitive to low energy photons, such that the source filter and the detector module cooperate to measure the intensity of photons within a narrowed band of energy levels.
- the detector is chosen such that it has an energy response which increases efficiency rapidly above a defined energy level, known as the "K-edge", for a particular detector material.
- K-edge a defined energy level
- Xenon gas-based detectors have a K-edge at about 35 KeV.
- Barium-based PSL imaging plates have a similar K-edge at about 37 KeV.
- Gadolinium scintillation screens have a K-edge at about 50 KeV.
- the filter is chosen, such that its K-edge is higher than that of the detector, such that the resultant system acts as a bandpass filter with a lower limit defined by the K-edge of the detector and an upper limit defined by the K-edge of the filter.
- At least one filter element has an elemental composition consisting of one or more of the lanthanide metals, such as lanthanum, cerium, praseodymium, neodymium, promethium, samarium, europium, gadolinium, terbium, dysprosium, holmium, erbium, thulium, ytterbium, lutetium or alloys thereof.
- This filter element can be employed alone or in combination with one or additional filter elements having an elemental composition of silver, cadmium, indium, tin or alloys thereof.
- one particularly useful combination of filter elements is the ternary combination of gadolinium, silver and tin.
- Each filter element can have a thickness ranging from about 0.01 to about 1.0 millimeters, preferably ranging from about 0.05 to about 0.5 millimeters.
- the collimator 34 is employed to eliminate scattering effects that may be caused during the filtering processes and produces a beam of radiation of appropriate size and shape for passage through the patient limb undergoing radiographic imaging to the detector module.
- a multiwire proportional chamber (MWPC) device 80 is illustrated as part of detector module 14.
- the chamber 80 is filled with a gas 82, which typically comprises xenon (Xe) and one or more buffer or additive gases, and contains an X-cathode grid 84, a Y-cathode grid 86 and an anode grid 88.
- the MWPC device 80 can be incorporated into a housing 78 having a window or otherwise radiation-transparent upper surface 79.
- Voltage source 90 produces a positive charge on the anode grid 88.
- the cathode grids 84, 86 are connected to data readout circuitry 92 to sense ionization events which occur when the photon radiation ionizes the xenon gas molecules within the chamber 80 and when the electrons, resulting from such ionization, drift to the positive anodes where charge multiplication takes place.
- FIG. 7 the operation of the MWPC device 80 is further illustrated schematically.
- five X-cathode and five Y-cathode lines 84, 86 are disposed in orthogonal relationship (each of which is actually composed of a series of parallel wires which are ganged together) forming a 5 ⁇ 5 coordinate matrix.
- the X-grid 84 and the Y-grid 86 are separated from each other by an anode grid 88.
- anode grid 88 When a photon of sufficient energy strikes a xenon gas molecule, the gas molecule is ionized, producing an electron-ion pair. The electron drifts towards the anode where the intense-field around the anode results in charge multiplication.
- the resulting avalanche 94 of electrons is sensed as a positive charge on the nearest wires of both the X-cathode and the Y-cathode grids, as well as a negative charge on the anode grid.
- the coincidence of the X-cathode signals, X i , the Y-cathode signals, Y i , and the anode signal, a permits the detector to accurately pinpoint the location of the radiation-induced ionization event.
- the grids are larger and finer than shown in FIG. 7 to provide higher spatial resolution.
- one embodiment for wrist imaging can employ 160 parallel Y-cathode wires and 80 parallel X-cathode wires, which are ganged as pairs to provide a 40 ⁇ 80 location grid.
- the spacing between adjacent wires can be approximately 1 millimeter.
- the cathode grids can be constructed on radiation-transparent (e.g., polyimide) printed circuit boards, for example, having copper strips 4-5 millimeters wide separated from each other by gaps of about 0.1 to about 0.5 millimeters in width.
- the anode grid 88 can be similarly constructed of a plurality of thin parallel (e.g., 20 micrometer gold plated tungsten) wires.
- the anode wires can be all wired together and, in operation, are all at the same voltage determined by the voltage source 90.
- the voltage will depend upon the pressure of the Xe gas within the chamber, as well as the spacing between the anode wires, and the anode-to-cathode spacing. For a typical pressure of 3 atmospheres and a spacing between anode wires of 2 millimeters, a positive voltage on the order of about 2,500 to about 5,000 volts relative to the cathode wires can be employed.
- FIG. 8-12 illustrate how the filter elements and detector cooperate to narrow the energy band of detected radiation in the present invention.
- the preferred energy band ranges from about 35 KeV to about 50 KeV.
- the graphs illustrate X-ray energies from about 27 KeV to 55 KeV; photons having energies below 27 KeV are typically absorbed by the patient's body and do not reach the detector.
- Xe gas is used as the radiation detecting element.
- Photons of energy below 34 Kev can ionize the Xe atoms and produce L-shell electrons, having energies equal to the energy of the photon minus 5.0 Kev (L-edge of Xe); likewise, photons with energy greater than 34.5 Kev can ionize the Xe atoms and produce K-shell electrons, the kinetic energy of which is the difference between the photon energy and 34.5 Kev, the K-edge for Xe.
- an electron with a given kinetic energy will be able to travel a finite range in the gas before it collides and ionizes other Xe atoms very close to the anode wire.
- the range of the electron can have a very deleterious effect on the spatial resolution of the image obtained by detecting the radiation transmitted though an object (i.e., bone) via the MWPC detector. For example, if a photon interacted at a point in the gas, created a free electron, and the electron traveled to another point where it caused an avalanche at a chamber wire, the initial photon interaction site in the gas could be one or two wire spacings away from it initial interaction site. The result is that information as to where the initial interaction occurred then becomes very imprecise and, hence, lost of spatial resolution in the image of the object occurs.
- a composite K-edge filter using 2 or more, preferably 3, pure metals stacked on one another is particularly useful in modifying the raw spectrum of a radiation source to achieve greater spatial resolution in the detector.
- metals with atomic numbers from 45 to 52 usually two metals
- at least one metal whose atomic number ranges from 57 to 65 give an optimum spatial resolution in Xe gas detectors.
- cadmium (Cd), silver (Ag), and gadolinium (Gd) can be particularly useful in combination.
- cadmium (Cd), tin (Sn) and samarium (Sm) form another useful combination filter.
- the transmission efficiency of a Cd-Ag-Gd filter combination is illustrated in FIG. 9. As can be seen, the filter combination provides for highest transmissivity in the range of about 42 to about 48 KeV.
- the results of using this Cd-Ag-Gd filter on the raw spectrum of the X-ray tube is shown in FIG. 10 where a distinctive narrowing of the energy band can be seen. Similar results in terms of limiting the upper range of photons can be achieved for higher X-ray source voltages than 55 Kv by modifying the thickness of the filter elements.
- the characteristic response of the detector can also be used to further narrow the energy band.
- FIG. 11 the efficiency of an MWPC device is plotted, showing that the Xe gas chamber is most effective in detecting photons having energies above 35 KeV and below 50 KeV.
- the overall response as illustrated in FIG. 12, is a well-defined energy band with sharp edges.
- the ultimate result of using the composite filter is threefold: first, by blocking low energy photons, it reduces the radiation dose otherwise absorbed by the patient and, second, by narrowing the energy band of the photons from the source, it improves the precision of measurement of the X-ray absorption coefficient for the tissue and bone and, thrid, dramatically diminishes the projected electron range in the detector, itself, thereby improving spatial resolution.
- Table 1 The particular advantages of the composite filter can be seen in Table 1 below, where the X-ray tube was controlled so as to impose an upper limit on the photon energy spectrum of 55 KeV and various filter combinations were compared under the same conditions (i.e., same total metal thickness regardless of composition):
- Another result of the composite filter lies in its ability to reduce the error typically inherent in bone mineral content (BMC) measurements due to the presence of fat in the soft tissues of the wrist and forearm.
- BMC bone mineral content
- to make bone mineral content measurements one typically measures the intensity of photons through a tissue equivalent bolus which is the same thickness as patient's limb and which contains a material having an absorption characteristic which is the same as tissue free of fat.
- the fat content can be significant (e.g., five to fifteen percent) in heavy individuals; and when it varies, it will lead to errors in measurement from year to year.
- the composite filter removes a large percentage of the photons in the range 20 to 35 Kev where the difference between the absorption coefficients for fat and tissue is largest, the observed differences in photon intensity through fat-containing tissue will be reduced substantially for the filtered X-ray spectra.
- FIGS. 13-15 illustrate a limb positioning and calibration apparatus 18 according to the present invention.
- a top view of the apparatus 18 is shown, including end blocks 60A, 60B, a plurality of top posts 66, a plurality of tracks 78, bottom plate 74 (which is composed of a tissue equivalent material), spring-loaded biasing elements 68A, 68B and two side blocks 70, 72, which engage the sides of the patient's limb in operation.
- the side blocks are composed of two different materials; for example, side block 70 can be composed of a bone equivalent material and side block 72 can be composed of a tissue equivalent material. Such materials are well known in the state-of-the-art of bone densitometry.
- FIGS. 14 and 15 show further features of the limb positioning and calibration apparatus 18 and illustrate its operation.
- the apparatus is shown in its opened position affixed to the housing 20 and disposed below the radiation source 12.
- the apparatus further includes a top cover 62 (at least a portion of which is composed of a tissue equivalent material 64) which is attached to the end blocks 60A, 60B via a plurality of post clamps 76.
- the apparatus 18 is shown in its closed position.
- the limb 16 is positioned in the apparatus 18 such that side blocks 70 and 72 are pressed by the biasing elements 68A, 68B against the sides of the limb 16 and preferably conform to its shape.
- the blocks 70, 72 can have a minimum gap, no tension, of about 3 centimeters when in the closed position.
- the blocks 70, 72 are separated and the limb (e.g. a wrist) is placed in the gap and then the blocks are allowed to close against the limb due to force of the biasing elements 68A, 68B.
- the blocks can be formed from a flexible plastic shell and filled with fluid, both the plastic and fluid being chosen to provide a calibration equivalence (e.g. bone equivalence in the case of block 70 and tissue equivalence in the case block 72).
- the patient's elbow is positioned so that the wrist and forearm are in approximately a straight line.
- the liquid-filled flexible material contained in blocks 70 and 72 will then form a tight junction and eliminate the gaps near the edge of the wrist that otherwise would occur.
- the bottom part of the wrist will lie on bottom plate 74, which is also preferably a bed of tissue-equivalent material encased in a plastic covering.
- a top plate 62 which includes an upper tissue equivalence material is lowered onto the wrist 16.
- the top and bottom plates 64, 74 of the apparatus 18 are made flush with the top and bottom of the wrist to give the same thickness as the wrist for the side blocks 70, 72.
- top plate 62 This is accomplished by forcing the top plate 62 down to compress the flexible side blocks on each side of the wrist.
- the vertical force necessary to compress the side block via the coverplate can be that applied by the operator or the coverplate itself can also be spring-loaded (not shown).
- a thin layer of flexible filled material 64 may also be incorporated in the bottom surface of the top plate 62 so as to ensure a precise fit. Then the coverplate can be secured by simply turning the screws 76 that lock the coverplate in place via the posts 66.
- the side blocks 70, 72 preferably incorporate two separate equivalence materials, one for tissue and one for bone, to permit calibration of the system. Such calibration is important to verify results from one machine to another or from measurement session to another, perhaps a year or more later.
- the system can be calibrated to negate possible changes in the energy spectrum of the radiation source or in the detector efficiency.
- calibration can be achieved by taking the natural logarithm of the intensity of radiation transmitted through the two equivalence materials. The ratio of these values can then be used to establish a measurement constant as follows:
- I b is the intensity of radiation transmitted through a first (e.g., bone equivalence) material
- I t is the intensity of radiation transmitted through a second (e.g., tissue equivalence) material
- U b is the bone absorption coefficient
- U t is the tissue absorption coefficient
- k is the measurement constant at which the system desirable operates.
- the calibration apparatus of the present invention thus permits the operator to measure the ratio of U b /U t and then dynamically adjust the system prior to patient imaging until a desired measurement constant is achieved. Such adjustments will typically involve varying the intensity of radiation source (e.g. by adjusting the voltage applied to the X-ray tube) or modifying the sensitivity of the detector module. Alternatively, any deviation from the desired measurement constant can instead be noted and used subsequently by the data acquisition or data processing modules to apply a correction factor to the data values. (It should also be apparent that the calibration system of the present invention is useful not only in single energy photo systems but also in so-called "dual photon" systems which detect radiation at two different radiation levels to image bones or conduct bone mineral content analyses.)
- the calibration system of the present invention also permits the user to accurately determine the thickness of the limb without resort to physical measurements (e.g., with a caliber). By comparing the logarithmic intensity values of radiation transmitted through the two calibration materials, the thickness of the materials can be calculated. Since the clamping means 18 will constraint the materials to have the same thickness as the limb, the limb thickness can also be inferred.
- this thickness T can be calculated according to the following equation: ##EQU1## where ln (I b ) is the natural logarithm of the measurement radiation intensity transmitted through a first calibration material (e.g., a bone equivalence medium) and ln (I t ) is the natural logarithm of the measurement radiation intensity transmitted through a second calibration material (e.g., a tissue equivalent material and U b and U t are the radiation absorption coefficients of the first and second calibration materials, respectively.
- a first calibration material e.g., a bone equivalence medium
- ln (I t ) is the natural logarithm of the measurement radiation intensity transmitted through a second calibration material
- U b and U t are the radiation absorption coefficients of the first and second calibration materials, respectively.
- FIG. 16 is a more detailed block diagram of the data readout circuitry 100 shown schematically in FIG. 6.
- the data readout circuitry 100 includes gated event digitizers 102A and 102B, highest priority encoders 104A and 104B, lowest priority encoders 106A and 106B and adders 108A and 108B.
- Each of the X-cathode signal lines are connected to gated event digitizer 102A which produces a synchronized set of digital signals representative of the X-cathode wires which have "sensed" a radiation event.
- the highest priority encoder 104A selects the highest X-cathode wire for which a positive digital signal has been generated and produces a digital number corresponding to the highest wire.
- the lowest priority encoder 106A selects the lowest X-cathode wire for which a positive digital signal has been generated and produces a similar digital number corresponding to the lowest wire. By summing the outputs of encoders 104A and 106A, and dividing by two, the origin of the radiation event on the X-coordinate axis can be determined.
- a similar operation is performed by digitizer 102B, encoders 104B and 106B, and adder 108B to arrive at the origin of the radiation event on the Y-coordinate axis.
- digitizer 102B encoders 104B and 106B
- adder 108B to arrive at the origin of the radiation event on the Y-coordinate axis.
- digitizer 102A includes a plurality of amplitude discriminator circuits 120A-120E. (Although five discriminator circuits are shown for purposes of illustration, the total number of discriminators in practice will equal the total number of X-cathode lines.) The lines which feed the discriminator circuits will be impressed with analog electrical signals of varying amplitudes depending upon their distance from the radiation event. The discriminator circuits are set by an adjustable threshold voltage to produce a digital output signal when the associated cathode line "senses", an electrical signal above the threshold.
- the digital outputs of the discriminator circuits 120A-E are registered by latch 124, which in turn, is triggered by a digital enablement signal generated by a similar discriminator circuit 122 connected to the anode grid.
- latch 124 which in turn, is triggered by a digital enablement signal generated by a similar discriminator circuit 122 connected to the anode grid.
- the coincidence of the anode and cathode signals above their respective thresholds yields a set of gated outputs to the encoders 104A and 106A of FIG. 16, thereby permitting a determination of the X-axis location.
- the "center of cluster" principle for determining radiation event locations in the priority encoders is shown in more detail in FIG. 18. Pulses of varying amplitudes are sensed by the N x-cathode lines. To determine the center of the signal cluster (X) and, hence, the origin on the X-axis of the radiation event, the priority encoders first determine (X1) and (X2), the locations of the end points of the pulse distribution. The center (X) can be calculated as follows: ##EQU2##
- FIG. 19 is a more detailed block diagram of the data acquisition module 22 shown schematically in FIG. 2.
- the data acquisition module 22 includes programmable logic device 130, address multiplexers 132A, 132B, incrementer 134, data multiplexers 136A, 136B, memory and multiplexer controllers 138A, 138B, dynamic random access memory banks 140A, 140B, logic controller 142 and a computer interface controller 144.
- the data acquisition module is formed from two separate memory subassemblies, each consisting, for example, of a standard 512 ⁇ 512 ⁇ 16 bit memory operated by a common programmed logic device 130 and logic controller 142.
- the logic controller 142 directs new events to the other memory bank and proceeds to transfer the contents of the first memory bank through the computer interface controller 144 to a data processing module, a data display module or a data storage module.
- the memory banks 140A, 140B each can each be constructed as a single 512 ⁇ 512 pixel image memory, or they can be partitioned to represent multiple smaller images, such as for example, four 256 ⁇ 256 pixel images or sixteen 128 ⁇ 128 images. Any of these arrangements permits the system to take multiple frames of data sequentially in time; each memory bank (or partition) thus acts as a histogram creating an image based upon radiation events recorded over time. Since the module 22 permits one memory bank to be transferred as the other is filled, it is possible to process or store data continuously.
- the interface controller 144 is preferably a small computer systems interface (SCSI) which is compatible with a wide variety of computer hardware, thereby allowing data transfer continuously, in either frame or list mode, to an external computer memory and/or an external magnetic or optical memory hard disk. Moreover, since the SCSI standard supports multiple SCSI devices on a single bus, it is possible to have several data acquisition systems operating in parallel from the same or multiple sources.
- SCSI small computer systems interface
- the logic controller 142, incrementer 134 and programmed logic device 130 of FIG. 19 cooperate to acquire and transfer data from the detector module, which can be a multiwire proportional chamber device, as described above, or another type of imaging device, such as a gamma camera or other photon imager.
- the logic controller 142 can be constructed from a general purpose microprocessor, such as the HD64180 microprocessor manufactured by Hitachi (San Jose, Calif.) and/or other commercially available parts.
- the microprocessor is preferably configured to define a direct memory data access channel for fast transfer of data between the memory banks 140A, 140B and the computer interface controller 144, and serves to refresh the dynamic random access memories periodically as well as initialize the programmed logic device 130, and multiplexers 132 and 136 and activate the interface controller 144.
- the microprocessor is programmed to perform tasks in response to requests by either the computer interface controller or other events that signal the end of a data acquisition cycle.
- the computer interface controller 144 will alert the logic controller 142 of a pending host computer request. This request can be, for example, an order to commence acquiring data from a wire chamber detector.
- the logic controller 142 will then decide which of the memory banks is to be used for this data acquisition cycle and set up the memory controllers 138 for that purpose.
- the logic controller 142 will also signal the programmed logic device 130 to begin acquiring data as described below.
- the logic controller 142 can initiate an internal timer to stop the data acquisition process after a predefined time has elapsed. Alternatively, it can program an event counter to stop data acquisition after a predefined number of events have been recorded.
- a host computer may also request from the logic controller 142 the initiation of data transfer from a previously acquired frame.
- the logic controller 142 will then direct the memory controllers 138 and multiplexers 136 and 132 to connect the corresponding memory bank 132A or 132B to the data path of the logic interface controller 144.
- the logic interface controller 142 will then set up the parameters of a direct memory access controller to sequentially read the locations of the memory bank and simultaneously write them to the computer interface controller 144.
- the incrementer 134 can be implemented by a parallel load counter, such as the model 74F779 counter manufactured by National Semiconductor, Fairchild Division (Mountainview, Calif.).
- the programmed logic device 130 can also be implemented with a commercially available, programmable logic device such as the model 20R4 device manufactured by Advanced Micro Devices (Sunnyvale, Calif.).
- FIG. 19 The various multiplexers and memory elements also shown in FIG. 19 are conventional in design and can be implemented with a wide variety of commercially available components. It should be clear that various changes, additions and subtractions to the circuitry of FIG. 19 can be made by those skilled in the art without departing from the spirit or scope of the invention.
- the memory and mutiplexer control functions of elements 138A, 138B can be supported by separate bit slice processor elements or can be implemented by firmware as part of the logic controller 142 or the programmed logic device 130.
- the functions of the components illustrated in FIG. 19 can be combined into larger integrated circuits or divided into smaller processing elements.
- the programmed logic device 130 is implemented as a state machine with eight states.
- the device 130 permits very fast acquisition of histogramming information by addressing a location in one of the memory banks, which corresponds to the particular coordinate location in the detector module where the radiation event occurred, retrieving the contents of the memory location and presenting it to the incrementer 134 which adds a "1" to those contents, and then rewrites the new number to the addressed memory contents.
- programmed logic device 130 can be further understood by reference to the timing diagrams of FIG. 20A and further illustrated by the schematic diagram of FIG. 20B. As shown, the state machine has the following sequence:
- the state machine In the idle mode, the state machine cycles between state #0 and state #1. If the detector module signals a new event by asserting the INFO READY line during state #0, the state machine switches to the active mode.
- State #1 is a timing safeguard to insure that the INFO READY signal has been received.
- state #2 the horizontal (or x) information is requested from the detector module, which is latched as row information in the DRAM-based memory bank 140.
- the vertical (or y) information is requested from the detector module, which is latched as column information in the memory bank 140.
- the present contents at the location X-Y in the memory bank 140 are retrieved.
- the output buffers are enabled, and the incrementer is ordered to load the contents (READ MEMORY signal).
- State #4 is a timing delay to ensure that the memory contents have been loaded into the incrementer.
- State #6 is again a timing delay to insure that the incrementer has completed it task.
- the programmed logic device can be driven at any clock rate compatible with the other components of the system. For example, when a 10 MegaHz clock is employed, the state machine can cycle through its eight states as shown in FIG. 20A in four clock cycles, or approximately 50 nanoseconds per state. A faster clock can also be employed to further reduce the data acquisition time.
- the input signal labeled REFADV indicates to the state machine that a DRAM-refresh cycle is about to take place.
- the state machine refrains from initiating the 8 state sequence while this signal is asserted.
- the DATA READY signal used by the wire chamber system to indicate availability of data, may be withdrawn by state #3, and must be withdrawn by state #7.
- the general purpose logic controller 142 While reading the histogramming memory 140, the general purpose logic controller 142, and not the programmed logic device 130, preferably, generates the necessary control signals to retrieve the information from the memory 140 and presents it to the processor interface 142. When data transfer from a memory element to the computer interface 144 is taking place, no data acquisition can be performed to that histogram memory and vice versa.
- the data acquisition module 22 provides for simultaneous buildup of a histogram, and transfer of data to the interface 144 of a previously-acquired histogram, by employing two identical histogramming memory banks 140A, 140B. This prevents loss of data acquisition information while reading previously acquired information.
- a plurality of two-position multiplexers or digital switches are incorporated for this purpose. They provide a connection of each of the histogramming memory banks 140A, 140B to the data acquisition or data transfer address and data lines.
- the multiplexers, for one of the histogramming memory banks, are always in the opposite position, as the ones for the remaining histogram memory.
- the one exception is during the DRAM refresh cycle (generated by the logic controller 142), in which both address multiplexers 132A, 132B are connected to the logic controller 142.
- FIG. 21 illustrates the data flow path during a typical portion of the acquisition cycle.
- the acquisition module 22 is shown in a state where memory bank 140A is acquiring data, while memory bank 140B is transferring data through the processor interface 142.
- the data processing module 24, a data display module 26 and a data storage module 28 can all be assembled from commercially available components.
- a wide variety of general purpose microcomputers, minicomputers or mainframe computers can be employed to process the data in the data processing module 24.
- a wide variety of computer monitors or higher resolution video display equipment can be employed in the display module 26.
- the image presented to the operator can be composed of pixels having gray-tone levels or color values which are representative of the detector responses at given areas of the patient's limb being imaged.
- the data storage module 28 can also be of conventional design and rely on any one of a number of known storage media, such as magnetic disks, diskettes, tape cassettes, optical disks, or solid-state random access memory elements.
- the data processing module 24 and the display module 26 can be used cooperatively to perform measurements of bone density on a particular bone structure or bone structures of the patient's limb being imaged.
- FIG. 22 an overall schematic view of the presentation of a patient's limb 16, in this instance, a wrist, to the system is shown.
- the bone structure of the patient's limb 16 in shown in phantom, except for that portion 150 which is actually viewed by the detector.
- the two principal bone structures in the wrist image are the radius 152 and the ulna 154.
- the data processing module of the present invention employs the styloid tips of the radius and the ulna bones as reference points, and establishes measurement areas in select regions of those bones that contain high percentages of trabecular bone, rather than cortical bone. These select regions, in the case of the radius and the ulna, have been determined to lie at specific distances from the respective styloid tips of these bones.
- measurement areas are selected that contain 60% or more of trabecular bone.
- the percentage of trabecular bone present in the radius and the ulna beyond 30 millimeters from the styloid tip of each bone drops below 10% and the percentage of cortical bone for both the radius and the ulna is approximately 95% at points beyond 40 millimeters from the styloid tip.
- the region associated with this 60% limit of trabecular bone usually begins at 2 millimeters and 4 millimeters, respectively, from the styloid tips of the radius and the ulna, respectively.
- the highest trabecular content by volume is typically found about 9 millimeters to 20 millimeters from the styloid tip in the radius and between about 4 millimeters and 12 millimeters from the styloid tip in the ulna.
- bone mineral content measurement regions are selected within each of these areas of high trabecular bone content.
- the operator can perform a bone mineral content measurement by employing the data processing module 24 and data display module 26 together.
- the system is initialized and an image is displayed to the clinician via a conventional computer monitor.
- the resolution of the display is not critical, so long as the bone edges are identifiable by the operator.
- the operator is prompted to use a cursor control device, such as a conventional mouse or joy stick, to place the styloid tip region markers around the respective styloid tips of the radius and the ulna.
- a cursor control device such as a conventional mouse or joy stick
- These markers may include, for example, hollow rectangular region markers, similar to styloid tip markers 166 and 168 shown in FIG. 23, which depicts a more detailed radiological image of a human wrist.
- the data processing module executing the method steps of this invention reads the position of the styloid tip region markers to generate initializing data for the location of the styloid tips.
- edges surrounding the styloid tips within each marker region can be determined by techniques known in the art.
- Bone edge finding can be implemented in accordance with any conventional edge finding algorithms, such as those disclosed in J. F. Canny, Tech. Rep. 720. Artificial Intelligence Laboratory, MIT (1983); "A Computational Approach To Edge Detection,” Vol. PAMI-8 IEEE Trans. Pattern Anal. Machine Intell. (1986) or F. Bergholm, "Edge Focusing,” Vol. PAMI-9 IEEE Trans. Pattern Anal. Machine Intell. (1987).
- Preferred edge finding algorithms include those which reduce noise contributions in the images and optimize the location of single edges rather than multiple spurious edges.
- the edge finding step is followed by a positional search executed on the edges of the styloid tip to identify one or more pixels that define the extremity of the styloid tip on each bone. Because the precise point-by-point shape of the respective radial and ulnal styloid tips can differ from individual to individual, it can be useful to employ more than one form of positional search or a curve fitting routine to determine the extremity of each styloid tip.
- the search or curved fitting routine can include a technique which finds the minimum-gradient tangent between points defining the edge, or one which locates the intersection of curves fitted to the styloid tip edges.
- the invention can utilize known techniques for initializing edge finding on the outer edge of the radius or ulna and progressing up and around the styloid tips to determine the edges of the bone extremities.
- utilization of operator input reduces processing time.
- the step of responding to operator input increases accuracy over a fully automated system, because it is likely that bone edges other than those of the radius and ulna will be present within the styloid tip region of interest.
- the phenomenon is illustrated in FIG. 23, where the edges of the scaphoid bone 156 near the radius, and the lunate bone 158 near the ulna, are proximate to the styloid tips 162 and 170.
- the styloid tip locations have been determined, one or more pixels identifying the extremity of the tip can be highlighted on the display for inspection by the operator, to confirm proper identification of the styloid tips.
- the data processing module executes a further edge finding step to determine the inner and outer longitudinal edges of the radius and the ulna.
- the term "longitudinal direction” is meant to define the direction corresponding to the long axis of the bone. For purposes of illustration, this longitudinal axis corresponds to the Y-coordinate axis of the detector.
- the initial starting position can be, for example, 3 centimeters from the styloid tip extremity of the radius. This initial position is preferable to ensure that the edge finding occurs in a region having an interosseous gap between the radius and the ulna.
- the interosseous separation between the ulna and the radius can range from 0.5-2.5 centimeters, depending upon the size and shape of the bones.
- edge finding process by initiating the edge finding process at a selected displacement, for example, approximately 3 centimeters, from the extremity of the styloid tip of the radius, ensures that the edge finding will occur in a region having a gap between the radius and the ulna.
- Longitudinal edge finding thus progresses along the bone from the initial position for a selected distance in the direction of a styloid tip and from the initial position for a selected distance in the opposite direction, to identify the edges of the bones on either side of the initial starting position the selected distance can be, for example, 2 centimeters or less from this initial position.
- each bone can be computed at selected Y-axis positions.
- a horizontal line referred to as the "region base line,” is then established and displayed at a Y-axis position on each bone having a selected displacement from the extremity of the respective styloid tip.
- a bone mineral content measurement can then be made on each bone. The measurement regions can be highlighted on display for the inspection of the operator.
- the measurement regions 170 and 172 can be, for example, rectangular in form having a selected height and a selected width based on the bone widths calculated previously.
- the measurement regions can extend from bone edge to bone edge, as does measurement region 172 on the ulna 154 in FIG. 24 or, the measurement regions can be reduced by a selected amount, to eliminate the outer and inner edges of the imaged bone from the BMC measurement region.
- Area 170 on the radius 152 in FIG. 24 is an example of a narrowed measurement region.
- narrowing the measurement region provides increased sensitivity in BMC measurements by reducing inaccuracies in measurement and computation which can result from highly irregular bone shapes.
- the BMC values can be derived in a known manner, by summing pixel intensity values in each measurement region. Bone density values can also be obtained by dividing the BMC value by the size of the measurement regions. The area covered by the measurement regions is proportional to the number of pixels in the measurement region. Those skilled in the art will appreciate that BMC values are proportional to the sum of the values of the data points in the measurement region, which in turn are represented by the pixel intensity values in the measurement region. In particular, the digital number associated with each pixel is typically proportional to the magnitude of the radiation received by the detector in the region corresponding to that pixel.
- Processing can be terminated when the BMC and density values have been obtained. Alternatively, the process can be re-initiated for processing of additional radiographic images.
- the present invention also permits a number of additional analyses to be performed from the stored data. For example, in successive bone mineral measurements of the same patient, the currently generated image data can be compared with previously stored images. A comparison can be made between the current orientation of the radius and ulna and the orientation of the radius and ulna represented in the stored data.
- a correction factor is established for the X and Y coordinate factors, to correlate previously generated and currently generated images and to correct small errors in the positioning of the patient's limb.
- Two-dimensional images such as the bone joint images processed in accordance with the present invention, have defined orientations within an X-Y orthogonal coordinate axis system. Differences in the orientation of the ulna and the radius in successive images can be detected with reference to this X-Y coordinate system.
- One technique is to evaluate the gradient at various points along the outer edges of each bone, to identify the flattest portion of the outer edge of each bone. Curve-fitted portions of the currently generated and previously generated bone joint images can than be compared to determine whether small translations and/or rotations of the imaged bones have occurred from one measurement to the next.
- rotational transformation should be corrected from image to image.
- Corrections for rotational transformations can be made by computing the angles between the flat portion of the outer bone edges for each of the first and second images. Based on the computed angles, corrections for both the X- and Y-axis positional pixel values can be provided by known computational techniques. Evaluation of these angles, and the positional values of the extremities of the styloid tips, permits images to be superimposed, or otherwise placed into point-to-point registration, to ensure that the vertical and horizontal bone distances are consistent from image to image for a given individual.
Abstract
Description
TABLE 1 ______________________________________ FILTER EFFECTS ON ELECTRON RANGE IN XENON Filter 1Filter 2Filter 3 No Filter (Ag) (Ag + Cd) (Ag + Cd + Gd) ______________________________________ electron 0.92 mm 0.56 mm 0.47 mm 0.36 mm range ______________________________________
ln (I.sub.b /I.sub.t)=U.sub.b /U.sub.t =k (1)
Claims (22)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US07/321,989 US4956859A (en) | 1989-03-10 | 1989-03-10 | Source filter for radiographic imaging |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US07/321,989 US4956859A (en) | 1989-03-10 | 1989-03-10 | Source filter for radiographic imaging |
Publications (1)
Publication Number | Publication Date |
---|---|
US4956859A true US4956859A (en) | 1990-09-11 |
Family
ID=23252917
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US07/321,989 Expired - Lifetime US4956859A (en) | 1989-03-10 | 1989-03-10 | Source filter for radiographic imaging |
Country Status (1)
Country | Link |
---|---|
US (1) | US4956859A (en) |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5090040A (en) * | 1989-03-10 | 1992-02-18 | Expert Image Systems, Inc. | Data acquisition system for radiographic imaging |
EP0499975A1 (en) * | 1991-02-20 | 1992-08-26 | Matsushita Electric Industrial Co., Ltd. | X-ray apparatus employing a K-edge filter |
EP0595733A2 (en) * | 1992-10-22 | 1994-05-04 | Eastman Kodak Company | Imaging combination for detecting soft tissue anomalies |
US5487661A (en) * | 1993-10-08 | 1996-01-30 | Dentsply International, Inc. | Portable dental camera and system |
JP2001190531A (en) * | 1999-11-01 | 2001-07-17 | General Electric Co <Ge> | Imaging system including radiation filter for x-ray imaging |
US20020044676A1 (en) * | 2000-09-07 | 2002-04-18 | Guo-Qing Wei | Method and system for automatic computed radiography (CR) image composition by white band detection and consistency rechecking |
US6405068B1 (en) | 1997-02-04 | 2002-06-11 | Osteometer Meditech A/S | Diagnostic imaging in relation to arthritic conditions |
US20030015208A1 (en) * | 2001-05-25 | 2003-01-23 | Philipp Lang | Methods to diagnose treat and prevent bone loss |
US20030112921A1 (en) * | 2000-10-11 | 2003-06-19 | Philipp Lang | Methods and devices for analysis of x-ray images |
US20040136500A1 (en) * | 2002-07-26 | 2004-07-15 | Jeol Ltd. | X-ray analyzer |
US20050259793A1 (en) * | 2004-05-19 | 2005-11-24 | Yeo In H | Medical phantom, holder and method of use thereof |
US20050286678A1 (en) * | 2004-06-28 | 2005-12-29 | Pediatric Imaging Technology, Llc | Ionizing radiation imaging system and method with decreased radiation dose |
US20070047794A1 (en) * | 2000-10-11 | 2007-03-01 | Philipp Lang | Methods and devices for analysis of x-ray images |
US20070274443A1 (en) * | 2006-05-25 | 2007-11-29 | Schlumberger Technology Corporation | Apparatus and method for fluid density determination |
US20080031412A1 (en) * | 2003-09-19 | 2008-02-07 | Imaging Therapeutics, Inc. | Method for Bone Structure Prognosis and Simulated Bone Remodeling |
US20090015119A1 (en) * | 2007-07-11 | 2009-01-15 | Maytag Corp. | Outer support body for a drawer-type dishwasher |
US20090207970A1 (en) * | 2000-08-29 | 2009-08-20 | Imaging Therapeutics, Inc. | Methods and Devices for Quantitative Analysis of X-Ray Images |
US20100098212A1 (en) * | 2000-08-29 | 2010-04-22 | Imaging Therapeutics, Inc. | Methods and Devices for Quantitative Analysis of X-Ray Images |
US7840247B2 (en) | 2002-09-16 | 2010-11-23 | Imatx, Inc. | Methods of predicting musculoskeletal disease |
US7995822B2 (en) | 2003-03-25 | 2011-08-09 | Imatx, Inc. | Methods for the compensation of imaging technique in the processing of radiographic images |
US8000441B2 (en) | 2000-08-29 | 2011-08-16 | Imatx, Inc. | Calibration devices and methods of use thereof |
US8073521B2 (en) | 2003-09-19 | 2011-12-06 | Imatx, Inc. | Method for bone structure prognosis and simulated bone remodeling |
US8600124B2 (en) | 2004-09-16 | 2013-12-03 | Imatx, Inc. | System and method of predicting future fractures |
US8639009B2 (en) | 2000-10-11 | 2014-01-28 | Imatx, Inc. | Methods and devices for evaluating and treating a bone condition based on x-ray image analysis |
JP2014042769A (en) * | 2012-08-29 | 2014-03-13 | Fujifilm Corp | Bone mineral quantity analysis method and bone mineral quantity analysis system, and recording medium |
US8939917B2 (en) | 2009-02-13 | 2015-01-27 | Imatx, Inc. | Methods and devices for quantitative analysis of bone and cartilage |
US8965075B2 (en) | 2002-09-16 | 2015-02-24 | Imatx, Inc. | System and method for predicting future fractures |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3715588A (en) * | 1970-10-26 | 1973-02-06 | Norland Corp | Bone mineral analyzer |
US3772521A (en) * | 1971-08-30 | 1973-11-13 | Univ California | Radiation camera and delay line readout |
-
1989
- 1989-03-10 US US07/321,989 patent/US4956859A/en not_active Expired - Lifetime
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3715588A (en) * | 1970-10-26 | 1973-02-06 | Norland Corp | Bone mineral analyzer |
US3772521A (en) * | 1971-08-30 | 1973-11-13 | Univ California | Radiation camera and delay line readout |
Non-Patent Citations (4)
Title |
---|
"Balanced Filters for Silver KαX-Rays", Berman et al., Review of Scientific Instruments, vol. 41, #6, 6/1970, pp. 870-873. |
"Definition and Measurement of Means of Filtered X Radiation", Grudskii et al., Raduim Institute Leningrad, Aug. 1974. |
Balanced Filters for Silver K X Rays , Berman et al., Review of Scientific Instruments, vol. 41, 6, 6/1970, pp. 870 873. * |
Definition and Measurement of Means of Filtered X Radiation , Grudskii et al., Raduim Institute Leningrad, Aug. 1974. * |
Cited By (59)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5090040A (en) * | 1989-03-10 | 1992-02-18 | Expert Image Systems, Inc. | Data acquisition system for radiographic imaging |
EP0499975A1 (en) * | 1991-02-20 | 1992-08-26 | Matsushita Electric Industrial Co., Ltd. | X-ray apparatus employing a K-edge filter |
US5285489A (en) * | 1991-02-20 | 1994-02-08 | Matsushita Electric Industrial Co., Ltd. | K-edge filter and an x-ray apparatus employing such a filter |
US5365567A (en) * | 1991-02-20 | 1994-11-15 | Matsushita Electric Industrial Co., Ltd. | K-edge filter and X-ray apparatus employing the same |
EP0595733A2 (en) * | 1992-10-22 | 1994-05-04 | Eastman Kodak Company | Imaging combination for detecting soft tissue anomalies |
EP0595733A3 (en) * | 1992-10-22 | 1994-10-26 | Eastman Kodak Co | Imaging combination for detecting soft tissue anomalies. |
US6132211A (en) * | 1993-10-08 | 2000-10-17 | Dentsply International Inc. | Portable dental camera, system and method |
US5836762A (en) * | 1993-10-08 | 1998-11-17 | Dentsply International Inc. | Portable dental camera system and method |
US5487661A (en) * | 1993-10-08 | 1996-01-30 | Dentsply International, Inc. | Portable dental camera and system |
US6405068B1 (en) | 1997-02-04 | 2002-06-11 | Osteometer Meditech A/S | Diagnostic imaging in relation to arthritic conditions |
JP2001190531A (en) * | 1999-11-01 | 2001-07-17 | General Electric Co <Ge> | Imaging system including radiation filter for x-ray imaging |
US6418193B1 (en) * | 1999-11-01 | 2002-07-09 | General Electric Company | Imaging system including radiation filter for x-ray imaging |
US8649481B2 (en) | 2000-08-29 | 2014-02-11 | Imatx, Inc. | Methods and devices for quantitative analysis of X-ray images |
US8588365B2 (en) | 2000-08-29 | 2013-11-19 | Imatx, Inc. | Calibration devices and methods of use thereof |
US8068580B2 (en) | 2000-08-29 | 2011-11-29 | Imatx, Inc. | Methods and devices for quantitative analysis of x-ray images |
US8031836B2 (en) | 2000-08-29 | 2011-10-04 | Imatx, Inc. | Methods and devices for quantitative analysis of x-ray images |
US8000441B2 (en) | 2000-08-29 | 2011-08-16 | Imatx, Inc. | Calibration devices and methods of use thereof |
US20100098212A1 (en) * | 2000-08-29 | 2010-04-22 | Imaging Therapeutics, Inc. | Methods and Devices for Quantitative Analysis of X-Ray Images |
US20090207970A1 (en) * | 2000-08-29 | 2009-08-20 | Imaging Therapeutics, Inc. | Methods and Devices for Quantitative Analysis of X-Ray Images |
US20020044676A1 (en) * | 2000-09-07 | 2002-04-18 | Guo-Qing Wei | Method and system for automatic computed radiography (CR) image composition by white band detection and consistency rechecking |
US6757418B2 (en) * | 2000-09-07 | 2004-06-29 | Siemens Corporate Research, Inc. | Method and system for automatic computed radiography (CR) image composition by white band detection and consistency rechecking |
US9767551B2 (en) | 2000-10-11 | 2017-09-19 | Imatx, Inc. | Methods and devices for analysis of x-ray images |
US9275469B2 (en) | 2000-10-11 | 2016-03-01 | Imatx, Inc. | Methods and devices for evaluating and treating a bone condition on x-ray image analysis |
US7660453B2 (en) | 2000-10-11 | 2010-02-09 | Imaging Therapeutics, Inc. | Methods and devices for analysis of x-ray images |
US20100197639A1 (en) * | 2000-10-11 | 2010-08-05 | Imaging Therapeutics, Inc. | Methods and Devices for Analysis of X-Ray Images |
US8625874B2 (en) | 2000-10-11 | 2014-01-07 | Imatx, Inc. | Methods and devices for analysis of x-ray images |
US20030112921A1 (en) * | 2000-10-11 | 2003-06-19 | Philipp Lang | Methods and devices for analysis of x-ray images |
US20070047794A1 (en) * | 2000-10-11 | 2007-03-01 | Philipp Lang | Methods and devices for analysis of x-ray images |
US8913818B2 (en) * | 2000-10-11 | 2014-12-16 | Imatx, Inc. | Methods and devices for evaluating and treating a bone condition based on X-ray image analysis |
US8639009B2 (en) | 2000-10-11 | 2014-01-28 | Imatx, Inc. | Methods and devices for evaluating and treating a bone condition based on x-ray image analysis |
US8000766B2 (en) | 2001-05-25 | 2011-08-16 | Imatx, Inc. | Methods to diagnose treat and prevent bone loss |
US9267955B2 (en) | 2001-05-25 | 2016-02-23 | Imatx, Inc. | Methods to diagnose treat and prevent bone loss |
US20030015208A1 (en) * | 2001-05-25 | 2003-01-23 | Philipp Lang | Methods to diagnose treat and prevent bone loss |
US20040136500A1 (en) * | 2002-07-26 | 2004-07-15 | Jeol Ltd. | X-ray analyzer |
US6968043B2 (en) * | 2002-07-26 | 2005-11-22 | Jeol Ltd. | X-ray analyzer |
US9460506B2 (en) | 2002-09-16 | 2016-10-04 | Imatx, Inc. | System and method for predicting future fractures |
US8965075B2 (en) | 2002-09-16 | 2015-02-24 | Imatx, Inc. | System and method for predicting future fractures |
US7840247B2 (en) | 2002-09-16 | 2010-11-23 | Imatx, Inc. | Methods of predicting musculoskeletal disease |
US8818484B2 (en) | 2002-09-16 | 2014-08-26 | Imatx, Inc. | Methods of predicting musculoskeletal disease |
US9155501B2 (en) | 2003-03-25 | 2015-10-13 | Imatx, Inc. | Methods for the compensation of imaging technique in the processing of radiographic images |
US7995822B2 (en) | 2003-03-25 | 2011-08-09 | Imatx, Inc. | Methods for the compensation of imaging technique in the processing of radiographic images |
US8781191B2 (en) | 2003-03-25 | 2014-07-15 | Imatx, Inc. | Methods for the compensation of imaging technique in the processing of radiographic images |
US8290564B2 (en) | 2003-09-19 | 2012-10-16 | Imatx, Inc. | Method for bone structure prognosis and simulated bone remodeling |
US8073521B2 (en) | 2003-09-19 | 2011-12-06 | Imatx, Inc. | Method for bone structure prognosis and simulated bone remodeling |
US20080031412A1 (en) * | 2003-09-19 | 2008-02-07 | Imaging Therapeutics, Inc. | Method for Bone Structure Prognosis and Simulated Bone Remodeling |
US20050259793A1 (en) * | 2004-05-19 | 2005-11-24 | Yeo In H | Medical phantom, holder and method of use thereof |
WO2006004638A2 (en) * | 2004-06-28 | 2006-01-12 | Pediatric Imaging Technnology, Llc | Ionizing radiation imaging system and method with decreased radiation dose |
WO2006004638A3 (en) * | 2004-06-28 | 2006-08-10 | Pediatric Imaging Technnology | Ionizing radiation imaging system and method with decreased radiation dose |
US20050286678A1 (en) * | 2004-06-28 | 2005-12-29 | Pediatric Imaging Technology, Llc | Ionizing radiation imaging system and method with decreased radiation dose |
US7203271B2 (en) * | 2004-06-28 | 2007-04-10 | Pediatric Imaging Technology, Llc | Ionizing radiation imaging system and method with decreased radiation dose |
US20070177714A1 (en) * | 2004-06-28 | 2007-08-02 | Pediatric Imaging Technology, Llc | Ionizing Radiation Imaging System and Method With Decreased Radiation Dose |
US8965087B2 (en) | 2004-09-16 | 2015-02-24 | Imatx, Inc. | System and method of predicting future fractures |
US8600124B2 (en) | 2004-09-16 | 2013-12-03 | Imatx, Inc. | System and method of predicting future fractures |
WO2007138402A1 (en) * | 2006-05-25 | 2007-12-06 | Schlumberger Technology B.V. | Apparatus and method for fluid density determination |
US20070274443A1 (en) * | 2006-05-25 | 2007-11-29 | Schlumberger Technology Corporation | Apparatus and method for fluid density determination |
US7507952B2 (en) | 2006-05-25 | 2009-03-24 | Schlumberger Technology Corporation | Apparatus and method for fluid density determination |
US20090015119A1 (en) * | 2007-07-11 | 2009-01-15 | Maytag Corp. | Outer support body for a drawer-type dishwasher |
US8939917B2 (en) | 2009-02-13 | 2015-01-27 | Imatx, Inc. | Methods and devices for quantitative analysis of bone and cartilage |
JP2014042769A (en) * | 2012-08-29 | 2014-03-13 | Fujifilm Corp | Bone mineral quantity analysis method and bone mineral quantity analysis system, and recording medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5090040A (en) | Data acquisition system for radiographic imaging | |
US4956859A (en) | Source filter for radiographic imaging | |
US5005195A (en) | Digital readout system for radiographic imaging | |
US5005196A (en) | Limb positioning and calibration apparatus for radiographic imaging | |
Phelps et al. | ECAT: a new computerized tomographic imaging system for positron-emitting radiopharmaceuticals | |
Demirkaya et al. | Image processing with MATLAB: applications in medicine and biology | |
Evans et al. | Performance evaluation of the PC-2048: a new 15-slice encoded-crystal PET scanner for neurological studies | |
US7620444B2 (en) | Systems and methods for improving usability of images for medical applications | |
US5103823A (en) | Motion correction in cardiology inspection systems | |
US3965353A (en) | Cross-sectional X-ray emission imaging system | |
EP0106267A2 (en) | Apparatus for the diagnosis of body structures into which a gamma-emitting radioactive isotope has been introduced | |
US3714429A (en) | Tomographic radioisotopic imaging with a scintillation camera | |
JPH0614914A (en) | Apparatus for correcting motion of patient | |
US7409240B1 (en) | System and method for imaging myocardial infarction | |
Hoffman et al. | A new tomograph for quantitative positron emission computed tomography of the brain | |
Baily | Video techniques for x‐ray imaging and data extraction from roentgenographic and fluoroscopic presentations | |
US6140649A (en) | Imaging attenuation correction employing simultaneous transmission/emission scanning | |
Koeppe et al. | Instrumentation for positron emission tomography: tomographs and data processing and display systems | |
Ben Younes et al. | A fully automated contour detection algorithm the preliminary step for scatter and attenuation compensation in SPECT | |
Chang et al. | New methods of examining gamma camera collimators | |
Thompson et al. | Imaging performance of a dynamic positron emission tomograph: Positome IIIp | |
Tsujimoto et al. | Fluctuation of quantitative values on acquisition time and the reconstruction conditions in 99mTc-SPECT | |
Williams et al. | Design of the neuro-ECAT®: A high-resolution, high efficiency positron tomograph for imaging the adult head or infant torso | |
Baru et al. | Digital x-ray imaging installation for medical diagnostics | |
Phelps et al. | Design considerations in positron computed tomography (PCT) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: LAHIVE, JOHN A., LIEPMANN, W. HUGO; COCKFIELD, JAM Free format text: SECURITY INTEREST;ASSIGNOR:EXPERT IMAGE SYSTEMS, INC.;REEL/FRAME:005091/0702 Effective date: 19890310 |
|
AS | Assignment |
Owner name: EXPERT IMAGE SYSTEMS, INC., MASSACHUSETTS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST.;ASSIGNORS:LANZA, RICHARD C.;VOTANO, JOSEPH R.;REEL/FRAME:005308/0011 Effective date: 19890310 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
AS | Assignment |
Owner name: LAHIVE & COCKFIELD, MASSACHUSETTS Free format text: SECURITY INTEREST;ASSIGNOR:EXPERT IMAGE SYSTEMS, INC.;REEL/FRAME:005483/0718 Effective date: 19890721 |
|
AS | Assignment |
Owner name: ENGELLENNER, THOMAS J., MASSACHUSETTS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST.;ASSIGNOR:EXPERT IMAGE SYSTEMS, INC., A CORP. OF MA;REEL/FRAME:006258/0027 Effective date: 19910923 |
|
FEPP | Fee payment procedure |
Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
FEPP | Fee payment procedure |
Free format text: PAT HOLDER CLAIMS SMALL ENTITY STATUS - SMALL BUSINESS (ORIGINAL EVENT CODE: SM02); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
AS | Assignment |
Owner name: HOLOGIC, INC., MASSACHUSETTS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:ENGELLENNER, THOMAS J.;REEL/FRAME:007815/0113 Effective date: 19951028 |
|
FPAY | Fee payment |
Year of fee payment: 8 |
|
REMI | Maintenance fee reminder mailed | ||
SULP | Surcharge for late payment | ||
FEPP | Fee payment procedure |
Free format text: PAT HOLDER NO LONGER CLAIMS SMALL ENTITY STATUS, ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: STOL); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
FPAY | Fee payment |
Year of fee payment: 12 |
|
AS | Assignment |
Owner name: GOLDMAN SACHS CREDIT PARTNERS L.P., CALIFORNIA Free format text: PATENT SECURITY AGREEMENT;ASSIGNOR:HOLOGIC, INC.;REEL/FRAME:020018/0818 Effective date: 20071022 Owner name: GOLDMAN SACHS CREDIT PARTNERS L.P.,CALIFORNIA Free format text: PATENT SECURITY AGREEMENT;ASSIGNOR:HOLOGIC, INC.;REEL/FRAME:020018/0818 Effective date: 20071022 |
|
AS | Assignment |
Owner name: GOLDMAN SACHS CREDIT PARTNERS L.P., AS COLLATERAL Free format text: PATENT SECURITY AGREEMENT;ASSIGNOR:HOLOGIC, INC.;REEL/FRAME:021301/0796 Effective date: 20080717 |
|
AS | Assignment |
Owner name: CYTYC SURGICAL PRODUCTS II LIMITED PARTNERSHIP, MA Free format text: TERMINATION OF PATENT SECURITY AGREEMENTS AND RELEASE OF SECURITY INTERESTS;ASSIGNOR:GOLDMAN SACHS CREDIT PARTNERS, L.P., AS COLLATERAL AGENT;REEL/FRAME:024944/0315 Effective date: 20100819 Owner name: CYTYC CORPORATION, MASSACHUSETTS Free format text: TERMINATION OF PATENT SECURITY AGREEMENTS AND RELEASE OF SECURITY INTERESTS;ASSIGNOR:GOLDMAN SACHS CREDIT PARTNERS, L.P., AS COLLATERAL AGENT;REEL/FRAME:024944/0315 Effective date: 20100819 Owner name: BIOLUCENT, LLC, CALIFORNIA Free format text: TERMINATION OF PATENT SECURITY AGREEMENTS AND RELEASE OF SECURITY INTERESTS;ASSIGNOR:GOLDMAN SACHS CREDIT PARTNERS, L.P., AS COLLATERAL AGENT;REEL/FRAME:024944/0315 Effective date: 20100819 Owner name: THIRD WAVE TECHNOLOGIES, INC., WISCONSIN Free format text: TERMINATION OF PATENT SECURITY AGREEMENTS AND RELEASE OF SECURITY INTERESTS;ASSIGNOR:GOLDMAN SACHS CREDIT PARTNERS, L.P., AS COLLATERAL AGENT;REEL/FRAME:024944/0315 Effective date: 20100819 Owner name: SUROS SURGICAL SYSTEMS, INC., INDIANA Free format text: TERMINATION OF PATENT SECURITY AGREEMENTS AND RELEASE OF SECURITY INTERESTS;ASSIGNOR:GOLDMAN SACHS CREDIT PARTNERS, L.P., AS COLLATERAL AGENT;REEL/FRAME:024944/0315 Effective date: 20100819 Owner name: CYTYC SURGICAL PRODUCTS III, INC., MASSACHUSETTS Free format text: TERMINATION OF PATENT SECURITY AGREEMENTS AND RELEASE OF SECURITY INTERESTS;ASSIGNOR:GOLDMAN SACHS CREDIT PARTNERS, L.P., AS COLLATERAL AGENT;REEL/FRAME:024944/0315 Effective date: 20100819 Owner name: CYTYC PRENATAL PRODUCTS CORP., MASSACHUSETTS Free format text: TERMINATION OF PATENT SECURITY AGREEMENTS AND RELEASE OF SECURITY INTERESTS;ASSIGNOR:GOLDMAN SACHS CREDIT PARTNERS, L.P., AS COLLATERAL AGENT;REEL/FRAME:024944/0315 Effective date: 20100819 Owner name: DIRECT RADIOGRAPHY CORP., DELAWARE Free format text: TERMINATION OF PATENT SECURITY AGREEMENTS AND RELEASE OF SECURITY INTERESTS;ASSIGNOR:GOLDMAN SACHS CREDIT PARTNERS, L.P., AS COLLATERAL AGENT;REEL/FRAME:024944/0315 Effective date: 20100819 Owner name: CYTYC SURGICAL PRODUCTS LIMITED PARTNERSHIP, MASSA Free format text: TERMINATION OF PATENT SECURITY AGREEMENTS AND RELEASE OF SECURITY INTERESTS;ASSIGNOR:GOLDMAN SACHS CREDIT PARTNERS, L.P., AS COLLATERAL AGENT;REEL/FRAME:024944/0315 Effective date: 20100819 Owner name: R2 TECHNOLOGY, INC., CALIFORNIA Free format text: TERMINATION OF PATENT SECURITY AGREEMENTS AND RELEASE OF SECURITY INTERESTS;ASSIGNOR:GOLDMAN SACHS CREDIT PARTNERS, L.P., AS COLLATERAL AGENT;REEL/FRAME:024944/0315 Effective date: 20100819 Owner name: HOLOGIC, INC., MASSACHUSETTS Free format text: TERMINATION OF PATENT SECURITY AGREEMENTS AND RELEASE OF SECURITY INTERESTS;ASSIGNOR:GOLDMAN SACHS CREDIT PARTNERS, L.P., AS COLLATERAL AGENT;REEL/FRAME:024944/0315 Effective date: 20100819 |