US20010053202A1 - Densitometry adapter for compact x-ray fluoroscopy machine - Google Patents
Densitometry adapter for compact x-ray fluoroscopy machine Download PDFInfo
- Publication number
- US20010053202A1 US20010053202A1 US09/745,704 US74570400A US2001053202A1 US 20010053202 A1 US20010053202 A1 US 20010053202A1 US 74570400 A US74570400 A US 74570400A US 2001053202 A1 US2001053202 A1 US 2001053202A1
- Authority
- US
- United States
- Prior art keywords
- positioner
- image
- ray
- recited
- detector
- 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.)
- Granted
Links
- 238000002594 fluoroscopy Methods 0.000 title claims abstract description 23
- 238000000326 densiometry Methods 0.000 title abstract description 8
- NJPPVKZQTLUDBO-UHFFFAOYSA-N novaluron Chemical compound C1=C(Cl)C(OC(F)(F)C(OC(F)(F)F)F)=CC=C1NC(=O)NC(=O)C1=C(F)C=CC=C1F NJPPVKZQTLUDBO-UHFFFAOYSA-N 0.000 claims abstract description 34
- 238000005259 measurement Methods 0.000 claims abstract description 28
- 238000012937 correction Methods 0.000 claims abstract description 20
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims description 112
- 230000005855 radiation Effects 0.000 claims description 26
- 230000037182 bone density Effects 0.000 claims description 13
- 230000009977 dual effect Effects 0.000 claims description 13
- 239000000463 material Substances 0.000 claims description 12
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 8
- 244000309466 calf Species 0.000 claims description 8
- 230000009467 reduction Effects 0.000 claims description 6
- 210000001519 tissue Anatomy 0.000 claims description 6
- 230000002238 attenuated effect Effects 0.000 claims description 3
- 238000001739 density measurement Methods 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 abstract description 24
- 230000008569 process Effects 0.000 description 75
- 238000012545 processing Methods 0.000 description 11
- 230000008859 change Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 10
- 210000004872 soft tissue Anatomy 0.000 description 10
- 230000000694 effects Effects 0.000 description 9
- 230000000875 corresponding effect Effects 0.000 description 8
- 238000012360 testing method Methods 0.000 description 8
- 210000000245 forearm Anatomy 0.000 description 7
- 230000009466 transformation Effects 0.000 description 7
- 235000004443 Ricinus communis Nutrition 0.000 description 6
- 240000000136 Scabiosa atropurpurea Species 0.000 description 6
- 235000012771 pancakes Nutrition 0.000 description 6
- 230000002829 reductive effect Effects 0.000 description 6
- 240000000528 Ricinus communis Species 0.000 description 5
- 230000002596 correlated effect Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000012935 Averaging Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000013507 mapping Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000001276 controlling effect Effects 0.000 description 3
- 210000003414 extremity Anatomy 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000004313 glare Effects 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000005266 casting Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000001816 cooling Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 210000002414 leg Anatomy 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000036961 partial effect Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- 210000000689 upper leg Anatomy 0.000 description 2
- 208000010392 Bone Fractures Diseases 0.000 description 1
- 229910000760 Hardened steel Inorganic materials 0.000 description 1
- 208000001132 Osteoporosis Diseases 0.000 description 1
- 229920005372 Plexiglas® Polymers 0.000 description 1
- 239000006096 absorbing agent Substances 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 229910021417 amorphous silicon Inorganic materials 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000007687 exposure technique Methods 0.000 description 1
- 239000010408 film Substances 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000005058 metal casting Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007491 morphometric analysis Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000000399 orthopedic effect Effects 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 239000004926 polymethyl methacrylate Substances 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000011410 subtraction method Methods 0.000 description 1
- 239000013077 target material Substances 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 238000004846 x-ray emission Methods 0.000 description 1
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/42—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4225—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using image intensifiers
-
- 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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4405—Constructional features of apparatus for radiation diagnosis the apparatus being movable or portable, e.g. handheld or mounted on a trolley
-
- 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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
- A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
- A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
-
- 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/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- 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
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
-
- 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/54—Control of apparatus or devices for radiation diagnosis
- A61B6/548—Remote control of the apparatus or devices
-
- 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/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G06T5/80—
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
- H04N23/80—Camera processing pipelines; Components thereof
- H04N23/81—Camera processing pipelines; Components thereof for suppressing or minimising disturbance in the image signal generation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N3/00—Scanning details of television systems; Combination thereof with generation of supply voltages
- H04N3/10—Scanning details of television systems; Combination thereof with generation of supply voltages by means not exclusively optical-mechanical
- H04N3/16—Scanning details of television systems; Combination thereof with generation of supply voltages by means not exclusively optical-mechanical by deflecting electron beam in cathode-ray tube, e.g. scanning corrections
- H04N3/22—Circuits for controlling dimensions, shape or centering of picture on screen
- H04N3/23—Distortion correction, e.g. for pincushion distortion correction, S-correction
- H04N3/233—Distortion correction, e.g. for pincushion distortion correction, S-correction using active elements
- H04N3/2335—Distortion correction, e.g. for pincushion distortion correction, S-correction using active elements with calculating means
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/30—Transforming light or analogous information into electric information
- H04N5/32—Transforming X-rays
- H04N5/321—Transforming X-rays with video transmission of fluoroscopic images
- H04N5/325—Image enhancement, e.g. by subtraction techniques using polyenergetic X-rays
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05G—X-RAY TECHNIQUE
- H05G1/00—X-ray apparatus involving X-ray tubes; Circuits therefor
- H05G1/08—Electrical details
- H05G1/26—Measuring, controlling or protecting
- H05G1/30—Controlling
- H05G1/38—Exposure time
- H05G1/42—Exposure time using arrangements for switching when a predetermined dose of radiation has been applied, e.g. in which the switching instant is determined by measuring the electrical energy supplied to the tube
- H05G1/44—Exposure time using arrangements for switching when a predetermined dose of radiation has been applied, e.g. in which the switching instant is determined by measuring the electrical energy supplied to the tube in which the switching instant is determined by measuring the amount of radiation directly
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05G—X-RAY TECHNIQUE
- H05G1/00—X-ray apparatus involving X-ray tubes; Circuits therefor
- H05G1/08—Electrical details
- H05G1/60—Circuit arrangements for obtaining a series of X-ray photographs or for X-ray cinematography
-
- 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
-
- 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/40—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4035—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis the source being combined with a filter or grating
-
- 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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4488—Means for cooling
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/419—Imaging computed tomograph
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/612—Specific applications or type of materials biological material
Definitions
- the invention relates generally to x-ray equipment and in particular to an adapter for bone density measurements as may be used with compact x-ray fluoroscopy equipment used for orthopedic and similar procedures.
- the present invention provides an adapter that may convert a compact fluoroscopy machine or other mobile x-ray source into a precision quantitative densitometer suitable for measuring bone mass or density such as may be helpful in the treatment and detection of osteoporosis.
- the invention provides a stand to be used with a fluoroscopy machine, the stand having a cradle for accurately locating the x-ray source and/or detector with respect to either the patient's forearm or foot.
- Special software is loaded to the computer of the fluoroscopy machine to operate the fluoroscopy machine in a quantitative dual energy mode and to adapt the fluoroscopy data to densiometric data.
- the invention includes a provision for a separate digital dual energy detector and if necessary an associated processing computer.
- FIG. 3 is a top view of the C-arm system of FIG. 1 with the articulated arm in partial phantom showing the four-bar linkage of the arm for extending the C-arn toward and away from the cart;
- FIG. 9 is a flow chart of a method of the present invention providing differently weighted noise reduction to different areas of the image based on motion in the areas of the image;
- FIG. 16 is a first embodiment of the second block of the flow chart of FIG. 14;
- FIG. 17 is a second embodiment of the second block of the flow chart of FIG. 14;
- FIG. 20 is a flow chart of the steps performed by the computer in correcting and transforming the image of FIGS. 11 and 19;
- FIG. 22 is a flow chart of the steps of calculating and removing scatter using the occluder of FIG. 21;
- FIG. 32 is a flow chart of software executed by a processing computer associated with the x-ray detector for providing quantitative densiometric data
- the articulated arm 19 may be slid horizontally toward the front of the cart 12 to provide a second positioning axis 25 of the x-ray machine 10 .
- a first pulley 28 is rotatively fixed in a vertical plane, attached to the portion of the slide 24 that may move with respect to the cart 12 , and is pivotally attached to a rigid arm 30 extending toward the front of the cart 12 .
- the other end of the rigid arm 30 supporting a second pulley 32 is also mounted to swivel with respect to arm 30 .
- motion of the collar 52 may be had in a vertical manner by means of the parallelogram linkage formed by pulleys 28 and 32 of the articulated arm 19 as shown in FIG. 2.
- Forward and backward motion away from and toward the cart 12 may be had by the second four bar linkage formed from pulleys 36 and 38 .
- a slight pivoting of the C-arm 56 about a vertical axis slightly to the rear of the collar 52 and concentric with the axis of pulley 38 may be had by means of the rotation between collar 52 and 50 of FIG. 4.
- Greater rotation of the C-arm about the vertical axis passing through pivot point 15 may be had by rotation of the cart about one of its stationary castors 14 .
- considerable flexibility in positioning the C-arm may be had.
- the x-ray tube 68 is connected to an x-ray tube power supply 74 which separately controls the current and voltage to the x-ray tube 68 based on signals received from the computer 22 as will be described.
- the control signals to the x-ray tube power supply 74 are encoded on a fiber optic within the cable 60 to be noise immune.
- Low voltage conductors are also contained within cable 60 to provide power to the x-ray tube power supply 74 from a low voltage power supply 76 positioned on the cart 12 .
- the CCD camera 84 provides a complete set of radiation data for an entire image 86 (a frame) periodically once every “frame interval” so that real-time image of a patient placed within the x-ray beam 80 may be obtained.
- Typical frame rates are in the order of thirty frames per second or thirty complete readouts of the CCD detector area to the computer 22 each second.
- a typical image 86 there will be some stationary object 90 such as bone and some moving object 92 such as a surgical instrument such as a catheter.
- some moving object 92 such as a surgical instrument such as a catheter.
- the bone 90 remains in the same place relative to the edge of the image 86 and 86 ′, however the surgical instrument 92 has moved. Accordingly, some pixels 88 ′ show no appreciable change between images 86 and 86 ′, whereas some other pixels 88 ′′ show a significant change between images 86 and images 86 ′.
- the computer 22 executes a stored program to compare current pixels of the image 86 ′ to the last pixels obtained from image 86 as indicated by process block 94 .
- This comparison is on a pixel by pixel basis with only corresponding pixels in the images 86 and 86 ′ compared.
- the difference between the values of the pixels 88 reflecting a difference in the amount of x-ray flux received at the CCD camera 84 , is mapped to a weight between zero and one, with greater difference between pixels 88 in these two images corresponding to larger values of this weight w. This mapping to the weighting is shown at process block 96 .
- the present invention develops an average image combining the values of the pixels acquired in each frame 86 , 86 ′ in which those pixels in the current image 86 ′ which exhibit very little change between images 86 and 86 ′ contribute equally to the average image, but those pixels in the current image 86 ′ that exhibit a great degree of change between images 86 and 86 ′ are given a substantially greater weight in the average image.
- a compromise is reached between using historical data to reduce noise and using current data so that the image accurately reflects changes.
- the value of each pixel displayed in the image is computed as follows.
- the present inventors have adopted a real-time digital re-mapping of radiation data from the CCD camera 84 to the image 86 to correct for any pincushion-type distortion.
- This remapping requires the imaging of the rectilinear grid 100 and an interpolation of the position of the radiation data received from the CCD camera 84 to new locations on the image 86 ′′ according to that test image.
- this remapping may be done on a real-time basis with good accuracy.
- these parameters may be input to the computer 22 and used at a transformation of received image 86 into image data 164 as indicated by process block 168 .
- new parameters of the polynomials may be entered by means of hand-held remote control wand 63 shown in FIG. 1.
- the entire area of the histogram 122 is assumed to represent the imaged object only (no background pixels). Such a situation represents an image of raw radiation only or a high dose image of a thin object with possible clipping.
- a reduced exposure rate will result as will be understood from the following description and the peak classification process, to now be described, is skipped.
- this right most peak 124 is compared against three empirically derived parameters indicated in the following Table 1: TABLE 1 Minimum Slope Range (MSR) Minimum necessary pixel range for which the slope of the peak must be monitonically increasing. Histogram Noise Level (HNL) Minimum height of the maximum value of the peak. Maximum Raw Radiation Width Maximum width of the detected peak (MRRW) with respect to the width of the entire histogram.
- MSR Minimum Slope Range
- HNL Histogram Noise Level
- MRRW Maximum width of the detected peak
- MSR and HNL are intended to prevent noise peaks and peaks caused by bad imaging elements in the CCD camera 84 or quantization of the video signal in the A to D conversion from being interpreted as background pixels.
- Peaks 124 with a suitable stretch of monotonically increasing slope 131 (shown in FIG. 13) according to the MSR value and that surpass the histogram noise level HNL 133 are evaluated against the MRRW parameter.
- This third evaluation compares the width 135 of the histogram 122 against the width of the entire histogram 122 .
- the MRRW value is intended to detect situations where the imaged object completely fills the imaging field and hence there are no unattenuated x-ray beams or background pixels being detected.
- a valid peak 124 will normally have a width 135 more than 33% of the total width of the histogram 122 .
- the current exposure rate is next compared to a reference exposure rate, in the preferred embodiment being 1.0 mR per frame, however this value may be refined after further clinical testing. If at process block 148 , the current exposure rate is within a “half fine-tune range” of the reference exposure rate, then the program proceeds to process block 150 , a fine tuning process block, and the amperage provided to the x-ray tube are adjusted in accordance to the disparity between the amperage and reference exposure rate. That is, if the current exposure is greater than the reference exposure rate, the amperage to the x-ray tube is reduced. The new value of amperage is compared against a predetermined range of amperage values (maximum beam current and minimum beam current values) so that the amperage value may never vary outside of this range.
- a predetermined range of amperage values maximum beam current and minimum beam current values
- mA is the amperage provided to the x-ray tube
- kVp is the voltage provided to the x-ray tube
- this value of ‘n’ is checked against threshold values intended to detect whether an erroneous value of n has been produced as a result of ‘clipping’ in the radiation data used to calculate exposure.
- clipping occurs when an increased dose of an element of the CCD camera produces no increase in the camera's output.
- n is tested to see if it is less than zero. This value of n is realized when the exposure rate of the second frame changes in the opposite direction of the tube voltage. This suggests a clipped histogram and therefore the program branches back to process block 154 to obtain a new second frame. This condition may also arrive from object motion between the first and second frame.
- the image 208 is then used to derive a scatter map.
- the amount of scatter at a given point will be a function of how many x-ray photons are received at points adjacent to the given point. For example, comparing the image 208 to a theoretical scatterless image 228 generally in an attenuated region 230 of the image 208 (e.g., under the spine 200 ), scatter will increase the apparent value in the image 208 as a result of radiation from nearby low attenuation regions scattering into the high attenuation region 230 . Conversely the apparent value at a low attenuation region 232 will be decreased because of the scatter into the high attenuation region.
- the mobile fluoroscopy machine 310 includes a mobile cart 312 supporting a computer 314 and monitor and keyboard 317 for receiving and processing digital x-ray image data.
- the cart 312 supports on one side an articulating arm assembly 316 terminating in a rotatable C-arm 318 .
- the C-arm supports, at the ends of the C, an image intensifier 320 and an x-ray source 322 opposed along an axis 324 so that the x-ray source 322 projects a cone-beam of x-ray radiation toward the image intensifier 320 along axis 324 .
- the pedestal 326 also provides on its upper surface a channel 330 extending across the axis 324 between the image intensifier 320 and the x-ray source 322 when the latter are positioned on the pedestal 326 .
- the channel 330 may receive a limb positioner 332 such as may be adapted to support a patient's foot or arm across the axis 324 for densiometric measurement.
- the pedestal 326 may be weighted so as to provide a stable surface for support of the x-ray source 322 and image intensifier 320 and to provide adequate support for the patient's limb.
- the height of the pedestal 326 is selected to be suitable for either arm or foot imaging.
- data is collected for three distinct images 384 with: 1) no x-ray exposure, 2) high energy x-ray exposure, and 3) low energy x-ray exposure.
- Each of the exposures is preserved as a separate image file in the memory of the computer 314 .
- the first exposure is used for correction routines to be described; the latter two exposures are used to deduce bone density according to methods well known in the art in which variations in high energy and low energy absorption are used to deduce the Compton scattering and atomic number of the material lying between the x-ray source 322 and the image intensifier 320 . As is understood in the art, these two measurements allow the amount of bone as opposed to soft tissue located in that image region to be accurately measured.
- the images 384 are log corrected reflecting the fact that attenuation is exponentially related to thickness.
- the images are now related to thickness, a dimension which will be important in the ultimate bone-density determination.
- the identified speckle of process block 418 is corrected by eliminating these identified pixels from subsequent calculation or by replacing them with a local average value.
Abstract
An accessory for fluoroscopy equipment is provided to support the x-ray tube and detector on a pedestal with respect to a patient limb for quantitative bone densitometry measurement. Software loaded into the associated digital imaging fluoroscopy equipment provides necessary correction of the images for the quantitative accuracy needed for bone densitometry. Alternatively, a specialized detector or extremely low form factor image intensifier may be inserted in the pedestal to be used in lieu of the fluoroscopy equipment detector. A similar software correction is performed on an associated computer when a separate detector must be used.
Description
- This application is based on U.S. provisional application No. 60/080,164 filed Mar. 31, 1998 and is a continuation in part of U.S. application Ser. No. 08/814,800 filed Mar. 10, 1997 and is a continuation in part of U.S. application Ser. No. 09/006,358 filed Jan. 13, 1998 which is a continuation-in -part of PCT Application Ser. No. 97/02770 designating the United States filed Feb. 21, 1997 claiming the benefit of provisional application No. 60/011,993 filed Feb. 21, 1996. This provisional application is incorporated by reference herein.
- The invention relates generally to x-ray equipment and in particular to an adapter for bone density measurements as may be used with compact x-ray fluoroscopy equipment used for orthopedic and similar procedures.
- Portable x-ray fluoroscopy machines provide an x-ray source held in opposition to an electronic image detector, typically on a C-arm, so that x-rays from the x-ray source are received by the image detector. The C-arm may slide through a collar so as to allow it to be rotated to different angles about the patient. Further, the collar may be supported by a pivoting arm providing additional freedom in the positioning of the C-arm.
- When the C-arm is correctly positioned, the x-ray source is activated and x-rays pass through the patient to be received by the image detector which provides electronic signals to a video monitor. For larger mobile C-arm systems, the video monitor is typically held on a separate cart or may be suspended from the ceiling on a fixed bracket to be connected to the mobile unit when the mobile unit is in place.
- With improvements in electronic hardware and in particular the development of compact image intensifiers and CCD video cameras, it has become possible to build extremely compact mobile C-arm systems. Such systems may make use of increasingly powerful desktop computer technology for image processing and other tasks and may use compact digital printers for producing images.
- The present invention provides an adapter that may convert a compact fluoroscopy machine or other mobile x-ray source into a precision quantitative densitometer suitable for measuring bone mass or density such as may be helpful in the treatment and detection of osteoporosis.
- The invention provides a stand to be used with a fluoroscopy machine, the stand having a cradle for accurately locating the x-ray source and/or detector with respect to either the patient's forearm or foot. Special software is loaded to the computer of the fluoroscopy machine to operate the fluoroscopy machine in a quantitative dual energy mode and to adapt the fluoroscopy data to densiometric data. For fluoroscopy equipment not providing for digital imaging or that may not be easily operated in a dual energy mode, the invention includes a provision for a separate digital dual energy detector and if necessary an associated processing computer.
- Other objects, advantages, and features of the present invention will become apparent from the following specification when taken in conjunction with the accompanying drawings.
- FIG. 1 is a perspective view of the fluoroscopy machine suitable for use with the present invention showing a C-arm supporting an image intensifier/video camera and x-ray tube in opposition for rotation in a vertical plane, the C-arm held along a mid-line of a cart by an articulated arm attached to the side of the cart;
- FIG. 2 is a side view in elevation of the cart of FIG. 1 showing a slide attaching the articulated arm to the side of the cart and showing a four-bar linkage motion of the arm for elevation of the C-arm;
- FIG. 3 is a top view of the C-arm system of FIG. 1 with the articulated arm in partial phantom showing the four-bar linkage of the arm for extending the C-arn toward and away from the cart;
- FIG. 4 is a detail fragmentary view of an outer pivot of the articulated arm attached to the C-arm such as allows limited pivoting of a plane of rotation of the C-arm about a vertical axis;
- FIG. 5 is a detail view of the C-arn of FIG. 1 and the attached x-ray tube assembly showing the electrical cabling providing power to an x-ray tube power supply fitting into a groove in the C-arm and showing an abutment of the anode of the x-ray tube against the metal casting of the C-arm for heat sinking purposes;
- FIG. 6 is a schematic block diagram of the fluoroscopy machine of FIG. 1 showing the path of control of a remote x-ray tube power supply by a microprocessor and the receipt of data from the image intensifier/video camera by the microprocessor for image processing;
- FIGS. 7 and 8 are simplified images such as may be obtained by the system of FIG. 1 showing portions of the image having moving elements and portions having stationary elements;
- FIG. 9 is a flow chart of a method of the present invention providing differently weighted noise reduction to different areas of the image based on motion in the areas of the image;
- FIG. 10 is a figure similar to that of FIG. 7 showing an image of a rectilinear grid as affected by pincushion distortion in the image intensifier and video camera optics such as may provide a confusing image of a surgical tool being manipulated in real-time;
- FIG. 11 is a figure similar to FIG. 10 showing equal areas of the image that encompass different areas of the imaged object, such variation as may affect quantitative bone density readings;
- FIG. 12 is a plot of raw image data from the image intensifier/video camera as is translated into pixel brightness in the images of FIGS. 7, 8,10, and 11 by the microprocessor of FIG. 6 according to a non-linear mapping process such as provides noise equilibrium in the images and maximum dynamic range for clinical data;
- FIG. 13 is a histogram plotting values of data from the image intensifier/video camera versus the frequency of occurrence of data values showing an isolated Gaussian distribution at the right most side representing unattenuated x-ray values;
- FIG. 14 is a flowchart describing the steps taken by the programmed microprocessor of FIG. 6 to identify background pixels and remove them from a calculation of exposure rate used for controlling the remote x-ray tube power supply of FIG. 6;
- FIG. 15 is a detailed block diagram of the first block of the flow chart of FIG. 14;
- FIG. 16 is a first embodiment of the second block of the flow chart of FIG. 14;
- FIG. 17 is a second embodiment of the second block of the flow chart of FIG. 14;
- FIG. 18 is a detailed flow chart of the third block of the flow chart of FIG. 14;
- FIG. 19 is a schematic representation of a distorted image of FIG. 11 and a schematic representation of a corresponding undistorted image showing the variables used in the mathematical transformation of the distorted image to correct for rotation and distortion;
- FIG. 20 is a flow chart of the steps performed by the computer in correcting and transforming the image of FIGS. 11 and 19;
- FIG. 21 is a perspective view of an occluder placed in an x-ray beam prior to an imaged object and used for calculating scatter;
- FIG. 22 is a flow chart of the steps of calculating and removing scatter using the occluder of FIG. 21;
- FIG. 23 is a cross-sectional view through the occluder of an imaged object of FIG. 21 along line23-23, aligned with a graph depicting attenuation of x-rays as a function distance along the line of cross-section as well as theoretical attenuation without scatter and scatter components;
- FIG. 24 is a graphical representation of an adjustment of calculated scatter from the image of FIG. 23 based on normalizing points established by the occluder of FIG. 21;
- FIG. 25 is a perspective view similar to that of FIG. 1 showing a C-arm system similar to that of FIG. 1 in position on the densitometry cradle of the present invention to provide a beam of x-rays along a horizontal axis across the top of the cradle;
- FIG. 26 is a fragmentary exploded view of the x-ray detector and x-ray source of the C-arm system of FIG. 25 removed from the cradle and showing a removable foot positioner also removed from the cradle;
- FIG. 27 is an enlarged, fragmentary view of FIG. 26 showing the alternative fitting of a forearm positioner or the foot positioner within a channel of the cradle along the path of x-rays between the x-ray source and x-ray detector;
- FIG. 28 is a top plan view of the cradle of FIG. 25 but with the C-arm removed, showing the forearm positioner in use with a patient's arm;
- FIG. 29 is a side elevational view of the cradle of FIG. 28;
- FIG. 30 is a top plan view of the cradle of FIG. 28 but with the foot positioner in use with a patient's leg and showing the use of an auxiliary pancake image intensifier for use with fluoroscopy or other x-ray equipment not having suitable dual energy or digital imaging capabilities;
- FIG. 31 is a cross sectional view of the pancake image intensifier of FIG. 30 with the protective shrouding removed and taken along a plane including the x-ray beam axis showing the compact, high distortion configuration and an attached processing computer;
- FIG. 32 is a flow chart of software executed by a processing computer associated with the x-ray detector for providing quantitative densiometric data;
- FIG. 33 is a block diagram of the steps of scatter correction of the present invention; and
- FIG. 34 is a block diagram illustrating the removal of line correlated noise per the present invention.
- Referring now to FIG. 1, an
x-ray machine 10 per the present invention includes a generally box-shaped cart 12 havingcastors 14 extending downward from its four lower corners. Thecastors 14 have wheels rotating about a generally horizontal axis, and swiveling about a generally vertical axis passing along the edges of thecart 12.Castors 14, as are understood in the art, may be locked against swiveling and/or against rotation. - With one
castor 14 locked and the others free to rotate and swivel, apivot point 15 for thecart 12 is established with respect to the floor such as may be used as afirst positioning axis 11 for thex-ray machine 10. - Positioned on the top of the
cart 12 is aturntable 16 holding avideo monitor 18 and attachedkeyboard 20 for swiveling about a vertical axis for convenience of the user. Thevideo monitor 18 and thekeyboard 20 may swivel separately so that one operator may view thevideo monitor 18 while a second operates thekeyboard 20. - The
video monitor 18 and thekeyboard 20 allow for control of acomputer 22 contained in a shelf on thecart 12 open from the front of thecart 12. Thecomputer 22 may include a general microprocessor-type processor 23 and aspecialized image processor 27 for particular functions as will be described. Thecomputer 22 further includes a number of interface boards allowing it to provide control signals to various components of thex-ray machine 10 as will be described and to receive x-ray image data. In addition, thecomputer 22 receives signals from afoot switch 61 that is used to activate the x-ray system for a brief exposure. Control of thecomputer 22 may also be accomplished through aremote control wand 63 of a type known in the art. - Referring now also to FIG. 2, attached to the right side of the
cart 12 is ahorizontal slide 24 positioned to provide anattachment point 26 for an articulatedarm 19 supporting a substantially circular C-arm 56, which in turn holds anx-ray tube 68 and animage intensifier 82 and camera 84, in opposition, and facing each other as will be described below. The function of the x-ray tube, the image intensifier and the camera are well known in the prior art in the use of mobile C-arm type x-ray devices used for image display and are described in U.S. Pat. No. 4,797,907 hereby incorporated by reference as part of the prior art. The C-arm may be mass balanced, that is to say its weight may be distributed to reduce its tendency to rotate throughcollar 54 so that minimal frictional pressure may be used to prevent it from moving. - The articulated
arm 19 may be slid horizontally toward the front of thecart 12 to provide asecond positioning axis 25 of thex-ray machine 10. Afirst pulley 28 is rotatively fixed in a vertical plane, attached to the portion of theslide 24 that may move with respect to thecart 12, and is pivotally attached to arigid arm 30 extending toward the front of thecart 12. The other end of therigid arm 30 supporting asecond pulley 32 is also mounted to swivel with respect toarm 30. Abelt 34 wraps around a portion of the circumference of each ofpulleys pulleys arm 30 about thecenter point 26 ofpulley 28 causes rotation ofpulley 32 so that it maintains a fixed rotational orientation with respect to thecart 12 aspulley 32 and hence C-arm 56 is moved up and down along athird axis 37. The linkage, so created, is a variation of the “four bar linkage” well known in the art. - Helical tension springs (not shown for clarity) balance the
pulley 32 in rotative equilibrium aboutpoint 26 against the weight of the articulatedarm 19, C-arm 56, and other devices attached to thearm 19. - Attached to
pulley 32 is athird pulley 36 extending in a generally horizontal plane perpendicular to the plane ofpulley 32. Thethird pulley 32 is attached pivotally to a secondrigid arm 40 which at its other end holds anotherpulley 38 positioned approximately at themidline 41 of thecart 12. Themidline 41 symmetrically divides the left and right sides of thecart 12. - Portions of the circumference of
pulleys belt 44 so as to form a second four barlinkage allowing pulley 38 to move toward and away from thecart 12, along afourth positioning axis 45, withpulley 38 and C-arm 56 maintaining their rotational orientation with respect tocart 12. - Referring now to FIG. 4,
pulley 38 includes acenter shaft member 50 having a coaxialouter collar 52 to whichbelt 44 is attached. Astop 55 attached to theshaft 50 limits the motion of thecollar 52 in rotation with respect to theshaft 50 to approximately 26 degrees. Frictional forces betweenshaft 50 andcollar 52cause shaft 50 to maintain its rotational orientation with respect tocollar 52 and hence with respect topulley 36 until sufficient force is exerted onshaft 50 to displace it with respect tocollar 52. Thus pressure on the C-arm 56 can provide some pivoting motion of the C-arm about the axis of the pulley along the fifthpositional axis 55. - Referring now to FIGS. 1, 3 and4, attached to the
shaft 50 is a C-arm collar 52 supporting the arcuate C-arm 56 curving through an approximately 180 degree arc in a vertical plane substantially aligned with themidline 41 of thecart 12 as has been mentioned. Theshaft 50 may connect tocollar 52 so that the latter may swivel in about a horizontal axis bisecting the circle of the C-arm 56. This axis may be aligned with the center of mass of the C-arm 56 so that there is not a self-righting tendency of the C-arm or the axis may be placed above the axis of the C-arm so as to provide for a beneficial self righting action. This motion is orthogonal to that provided by motion ofshaft 50 and may augment that provided by thecastors 14. Techniques of balancing the C-arm in its various rotational modes, when this is desired, is taught by U.S. Pat. No. 5,038,371 to Janssen issued Aug. 6th, 1991 and hereby incorporated by reference as exemplifying the known prior art understood to all those of ordinary skill in the art. - As described above, motion of the
collar 52 may be had in a vertical manner by means of the parallelogram linkage formed bypulleys arm 19 as shown in FIG. 2. Forward and backward motion away from and toward thecart 12 may be had by the second four bar linkage formed frompulleys arm 56 about a vertical axis slightly to the rear of thecollar 52 and concentric with the axis ofpulley 38 may be had by means of the rotation betweencollar pivot point 15 may be had by rotation of the cart about one of itsstationary castors 14. Thus, considerable flexibility in positioning the C-arm may be had. - Referring now to FIG. 5, the C-
arm 56 is an aluminum casting having formed along its outer circumference achannel 58 into which acable 60 may be run as will be described. C-arm 56 has a generally rectangular cross-section taken along a line of radius of the C-arm arc. Each corner of that rectangular cross-section holds a hardenedsteel wire 62 to provide a contact point forcorner bearings 64 within thecollar 52. Thecorner bearings 64 support the C-arm 56 but allow movement of the C-arm 56 along its arc through thecollar 52. - A
cable guide pulley 66 positioned over thechannel 58 and having a concave circumference feeds thecable 60 into thechannel 58 as the C-arm moves preventing tangling of thecable 60 or its exposure at the upper edge of the C-arm 56 when the C-arm 56 is rotated. The excess length ofcable 60 loops out beneath thecollar 52. - Referring now to FIGS. 5 and 6, the C-arm supports at one end a generally
cylindrical x-ray tube 68 having acathode 70 emitting a stream of electrons against a fixedanode 72. The conversion efficiencies of x-ray tubes are such that theanode 72 can become quite hot and typically requires cooling. In the present invention, theanode 72 is positioned to be bolted against the aluminum casting of the C-arm 56 thereby dissipating its heat into a large conductive metal structure of the C-arm 56. - The
x-ray tube 68 is connected to an x-raytube power supply 74 which separately controls the current and voltage to thex-ray tube 68 based on signals received from thecomputer 22 as will be described. The control signals to the x-raytube power supply 74 are encoded on a fiber optic within thecable 60 to be noise immune. Low voltage conductors are also contained withincable 60 to provide power to the x-raytube power supply 74 from a lowvoltage power supply 76 positioned on thecart 12. - During operation, an
x-ray beam 80 emitted from thex-ray tube 68 passes through a patient (not shown) and is received by animage intensifier 82 and recorded by a charge couple device (“CCD”) camera 84 such as is well known in the art. The camera provides digital radiation values to thecomputer 22 inversely proportional to the x-ray absorption of the imaged object for processing as will be described below. Each radiation value is dependent on the intensity of x-ray radiation received at a specific point on the imaging surface of theimage intensifier 82. - Referring now to FIGS. 6 and 7, the data collected by the CCD camera84 may be used to provide an
image 86 displayed onvideo monitor 18. As will be described in more detail below, the CCD camera receiving a light image from theimage intensifier 82 at a variety of points, provides data to the computer which maps the data from the CCD camera 84 to apixel 88 in theimage 86. For convenience, the data from the CCD camera 84 will also be termed radiation data reflecting the fact that there is not necessarily a one-to-one correspondence between data detected by the CCD camera 84 andpixels 88 displayed on thevideo monitor 18. - The CCD camera84 provides a complete set of radiation data for an entire image 86 (a frame) periodically once every “frame interval” so that real-time image of a patient placed within the
x-ray beam 80 may be obtained. Typical frame rates are in the order of thirty frames per second or thirty complete readouts of the CCD detector area to thecomputer 22 each second. - Each frame of data is stored in the memory of the
computer 22 and held until after complete storage of the next frame of data. The memory of thecomputer 22 also holds an average frame of data which represents an historical averaging of frames of data as will now be described and which is normally used to generate the image on thevideo monitor 18. - In a
typical image 86, there will be somestationary object 90 such as bone and some movingobject 92 such as a surgical instrument such as a catheter. In asecond image 86′ taken one frame after theimage 86, thebone 90 remains in the same place relative to the edge of theimage surgical instrument 92 has moved. Accordingly, somepixels 88′ show no appreciable change betweenimages other pixels 88″ show a significant change betweenimages 86 andimages 86′. - Referring now to FIG. 9, as data arrives at the
computer 22, thecomputer 22 executes a stored program to compare current pixels of theimage 86′ to the last pixels obtained fromimage 86 as indicated byprocess block 94. This comparison is on a pixel by pixel basis with only corresponding pixels in theimages pixels 88, reflecting a difference in the amount of x-ray flux received at the CCD camera 84, is mapped to a weight between zero and one, with greater difference betweenpixels 88 in these two images corresponding to larger values of this weight w. This mapping to the weighting is shown atprocess block 96. - Thus
pixels 88″, whose value changes almost by the entire range of pixel values betweenimages pixels 88′ which have no change betweenimages pixels 88 being neither unchanged nor radically changed will receive a value somewhere between zero and one. - Generally, because the amount of x-ray fluence in the
beam 80 is maintained at a low level to reduce the dose to the patient, theimages images pixel 88 over a number of frames of acquisition effectively increasing the amount of x-ray contributing to the image of that pixel. - Nevertheless, this averaging process tends to obscure motion such as exhibited by the
surgical instrument 92. Accordingly, the present invention develops an average image combining the values of the pixels acquired in eachframe current image 86′ which exhibit very little change betweenimages current image 86′ that exhibit a great degree of change betweenimages - P i=(1−w)P i−1 +wP i,t (1)
- where Pi−1 is a pixel in the previous average image, w is the weighting factor described above and Pi,t is the current data obtained from the CCD camera 84. This effective merger of the new data and the old data keyed to the change in the data is shown at
process block 98. - Referring now to FIG. 10, an
image 86″ of arectilinear grid 100 positioned in thex-ray beam 80 will appear to have a barrel or pincushion shape caused by distortion of theimage intensifier 82 and the optics of the CCD camera 84. During a real-time use of theimage 86″ by a physician, this distortion may cause confusion by the physician controlling atool 102. For example,tool 102 may be a straight wire shown by the dotted line, but may display animage 86′ as a curved wire whose curvature changes depending on the position of thetool 102 within theimage 86. This distortion thus may provide an obstacle to a physician attempting to accurately place thetool 102 with respect to an object within theimage 86′. - Referring now to FIG. 1, the distortion of
image 86″ also means that two equal area regions of interest 105 (equal in area with respect to the image) do not encompass equal areas of thex-ray beam 80 received by theimage intensifier 82. Accordingly, if the data from the CCD camera 84 is used for quantitative purposes, for example to deduce bone density, this distortion will cause an erroneous variation in bone density unrelated to the object being measured. - Accordingly, the present inventors have adopted a real-time digital re-mapping of radiation data from the CCD camera84 to the
image 86 to correct for any pincushion-type distortion. This remapping requires the imaging of therectilinear grid 100 and an interpolation of the position of the radiation data received from the CCD camera 84 to new locations on theimage 86″ according to that test image. By using digital processing techniques in adedicated image processor 27, this remapping may be done on a real-time basis with good accuracy. - Referring to FIG. 19, there are two types of distortion, isotropic and anisotropic. Isotropic distortion is rotationally symmetric (e.g. like barrel and pin cushion distortion). Anisotropic distortion is not rotationally symmetric. Both types of distortion and rotation are so-called third order aberrations which can be written in the form:
- Dx=r 2(Du−dv) (2)
- Dy=r 2(Dv+du) (3)
- where Dx and Dy are pixel shifts due to distortion; r is the distance from the correct position to the optical axis and D and d are distortion coefficients which are constant and u and v are correct pixel positions.
- Referring also FIG. 2, received
image 86 may exhibit pin cushion distortion evident if animage 86 of therectilinear grid 100 is made. The distortion is caused by the pixel shifts described above. - Equations 1 and 2 may be rewritten as third order two-dimensional polynomials, the case for equation (1) following:
- x=(a x +e x v+i x v 2 +m x v 3)+(b x +f x v+j x v 2 +n x v 3)u+(c x +g x v+k x v 2 +o x v 3)u 2+(d x +h x v+l x v 2 +p x v 3)u 3 (4)
- In these polynomials, ax and ay govern the x and y translation of the image, ex and by take care of scaling the output image, while ey and bx enable the output image to rotate. The remaining higher order terms generate perspective, sheer and higher order distortion transformations as will be understood to those of ordinary skill in the art. Thirty-two parameters are required for the two, third order polynomials. These parameters may be extracted by a program executed by the computer in an off-line (non-imaging) mode after imaging the known
grid 100 and comparing the distorted image of thegrid 100 to the knowngrid 100 to deduce the degrees of distortion. - Referring now to FIG. 19 in a first step of the correction process, the
grid 100 is imaged as indicated byprocess block 160 to determine the exact type of distortion present and to obtain values for the coefficients a through p of the above referenced polynomial equations. - At
process block 166, these parameters may be input to thecomputer 22 and used at a transformation of receivedimage 86 intoimage data 164 as indicated byprocess block 168. For rotation of theimage 164, new parameters of the polynomials may be entered by means of hand-heldremote control wand 63 shown in FIG. 1. - The transformation process generally requires a determination of the pixel shift for each
radiation pixel 163 of theinput image 86 which in turn requires an evaluation of the polynomials whose coefficients have been input. A number of techniques are known to evaluate such polynomials including a forward differencing technique or other techniques known in the art. These transformations provide values of u and v for animage pixel 170 corresponding to aparticular radiation pixel 163. - After the transformation of
process block 168, the u, v locations of the radiation pixels will not necessarily be centered at a pixel location defined by the hardware of the video monitor 18 which usuallyspaces pixels 170 at equal distances along a Cartesian axis. Accordingly, the transformed pixels must be interpolated to actual pixel locations as indicated byprocess block 172. - A number of interpolation techniques are well known including bilateral and closest neighbor interpolation, however in the preferred embodiment, a high resolution cubic spline function is used. A given value of an interpolated pixel170 (Pint) is deduced from a 4×4 block of transform pixels (Pij) in which it is centered as follows:
- Pint =f(n−2)X l +f(n−1)X 2 +f(n)X 3 +f(n+1)X 4 (5)
- where:
- X i =f(m−2)P i,1 +f(m−1)P i,2 +f(m)P i,3 +f(m+1)P i,4 (6)
- where:
- f(x)=(a+2)x 3+−(a+3)x 2+1 for x ε[0,1];
- f(x)=ax 3+−5ax 2+8ax−4a for x ε[1,2]; (7)
- f(x) is symmetrical about zero. In the preferred embodiment a=−0.5
- and where m and n are fractions indicating the displacement of the neighboring pixels Pi,j with respect to Pint in the x and y directions, respectively.
- At process block180, the transformed and interpolated image is displayed.
- Referring now to FIG. 12, the radiation data from the CCD camera84 are mapped to the brightness of the pixels of the
image 86 according to a second transformation. In the preferred embodiment, this mapping between CCD radiation data and image pixel brightness follows a nonlinear curve 103 based on the hyperbolic tangent and being asymptotically increasing to the maximum CCD pixel value. This curve is selected from a number of possibilities so that equally wide bands ofimage pixel brightness - The noise in the
image 86 is further reduced by controlling the fluence of thex-ray beam 80 as a function of the density of tissue of the patient within thebeam 80. This density is deduced from theimage 86 itself produced by the CCD camera 84. In response to the image data, a control signal is sent via the fiber optic strand within thecable 60 to the x-raytube power supply 74 positioned adjacent to the x-ray tube 68 (shown in FIG. 5). By positioning the x-raytube power supply 74 near thex-ray tube 68, extremely rapid changes in the power supplied to thex-ray tube 68 may be obtained. Distributed capacitances along high tension cables connecting thex-ray tube 68 to a stationary x-ray tube power supply are thus avoided in favor oflow voltage cable 60, and the shielding and inflexibility problems with such high tension cables are also avoided. - Referring now to FIGS. 13 and 14, a determination of the proper control signal to send to the x-ray
tube power supply 74 begins by analyzing theimage data 86 as shown inprocess block 120. The goal is to provide for proper exposure of an arbitrary object placed within thex-ray beam 80 even if it does not fill the field of view of the CCD camera 84. For this reason, it is necessary to eliminate consideration of the data from the CCD camera 84 that form pixels in the image corresponding to x-rays that bypass the imaged object and are unattenuated (“background pixels”). These background pixels may be arbitrarily distributed in theimage 86 and therefore, this identification process identifies these pixels based on their value. To do this, thecomputer 22 collects the values of the pixels from the CCD camera 84 in ahistogram 122 where the pixels are binned according to their values to create a multiple peaked plot. The horizontal axis of thehistogram 122 may for example be from 0 to 255 representing 8 bits of gray scale radiation data and the vertical axis may be a number of pixels having a particular value. - If there is a histogram value at horizontal value 255, and the maximum gray scale exposure recorded, the entire area of the
histogram 122 is assumed to represent the imaged object only (no background pixels). Such a situation represents an image of raw radiation only or a high dose image of a thin object with possible clipping. In assuming that thewhole histogram 122 may be used to calculate technique without removal of background pixels, a reduced exposure rate will result as will be understood from the following description and the peak classification process, to now be described, is skipped. - Otherwise, if there are no pixels with the maximum value of 225, the present invention identifies one
peak 124 in thehistogram 122 as background pixels indicated byprocess block 120 in FIG. 14. In identifying thispeak 124, thecomputer 22 examines thehistogram 122 from the brightest pixels (rightmost) to the darkest pixels (leftmost) assuming that the brightest pixels are more likely to be the unattenuated background pixels. Theprocess block 120 uses several predetermined user settings as will be described below to correctly identify thepeak 124. - Once the
peak 124 has been identified, the pixels associated with that peak are removed perprocess block 126 by thresholding or subtraction. In the thresholding process, pixels above athreshold value 138 below thepeak 124 are considered to be background pixels and are omitted from an exposure rate calculation. In the subtraction method, thepeak 124 itself is used as a template to identify pixels which will be removed. - At
process block 128, an exposure rate is calculated based on the values of the pixels in the remaining histogram data and atprocess block 130, an amperage and voltage value are transmitted via thecable 60 to the x-ray tube power supply and used to change the power to the x-ray tube. Generally, if the exposure rate is above a predetermined value, the amperage and voltage are adjusted to cut the x-ray emission from the x-ray tube, whereas if the exposure rate is below the predetermined value, the amperage and voltage are adjusted to boost the exposure rate to the predetermined value. - Referring now to FIGS. 13, 14 and15, the process of identifying background pixels will be explained in more detail.
Process block 120 includes as a first step, an identification of a rightmost peak 124 in the histogram 122 (shown in FIG. 13) as indicated bysubprocess block 132. - At succeeding
subprocess block 134, this rightmost peak 124 is compared against three empirically derived parameters indicated in the following Table 1:TABLE 1 Minimum Slope Range (MSR) Minimum necessary pixel range for which the slope of the peak must be monitonically increasing. Histogram Noise Level (HNL) Minimum height of the maximum value of the peak. Maximum Raw Radiation Width Maximum width of the detected peak (MRRW) with respect to the width of the entire histogram. - Specifically at
subprocess block 134, each identifiedpeak 124 is tested against the three parameters indicated in Table 1. In the description in Table 1, “width” refers to the horizontal axis of thehistogram 122 and hence a range of pixel values, whereas “height” refers to a frequency of occurrence for pixels within that range, i.e., the vertical axis of thehistogram 122. - These first two tests, MSR and HNL, are intended to prevent noise peaks and peaks caused by bad imaging elements in the CCD camera84 or quantization of the video signal in the A to D conversion from being interpreted as background pixels.
-
Peaks 124 with a suitable stretch of monotonically increasing slope 131 (shown in FIG. 13) according to the MSR value and that surpass the histogramnoise level HNL 133 are evaluated against the MRRW parameter. This third evaluation compares thewidth 135 of thehistogram 122 against the width of theentire histogram 122. The MRRW value is intended to detect situations where the imaged object completely fills the imaging field and hence there are no unattenuated x-ray beams or background pixels being detected. Avalid peak 124 will normally have awidth 135 more than 33% of the total width of thehistogram 122. - At
decision block 136 if the peak 124 passes the above tests, the program proceeds to process block 126 as indicated in FIG. 14. Otherwise, the program branches back to process block 132 and the next peak to the left is examined against the tests ofprocess block 134 until a passing peak is found or no peak is found. If no peak is found, it is assumed that there are no background pixels and a raw exposure value is calculated from all pixels as described above. - Assuming that a peak124 passes the tests of Table 1, then at process block 126 background pixels identified by the
peak 124 selected at process block 120 are eliminated. - In a first method of eliminating background pixels indicated at FIG. 16, a
magnitude threshold 138 within thehistogram 122 is identified. Pixels having values above this threshold will be ignored for the purpose of selecting an exposure technique. Thethreshold 138 is established by identifying thecenter 140 of the peak 124 (its maximum value) and subtracting from the value of the center a value s being the distance between the start of the peak 124 as one moves leftward and the maximum 140. The area under thehistogram 122 for values lower than thethreshold 138 is computed to deduce a raw exposure value which will be used as described below. - In a second embodiment, the shape of the
histogram peak 124 from the start of the peak as one moves leftward to its maximum 140 is reflected about a vertical line passing through the maximum 140 and subtracted from thehistogram peak 124 to the left of the vertical line. This approach assumes that thepeak 124 of the background pixels is symmetrical and thus this method better accommodates some overlap between the object pixels and the background pixels in thehistogram 122. Again, the remaining pixels of thehistogram 122 are summed (by integration of the area under thehistogram 122 minus the area of the peak 124 as generated by the reflection) to provide a raw exposure value. - Referring now to FIG. 18, the raw exposure value is transformed by the known transfer characteristics of the CCD camera (relating actual x-ray dose to pixel value) to produce a calculated current exposure rate as indicated at
process block 144. - Referring to process block146, the current exposure rate is next compared to a reference exposure rate, in the preferred embodiment being 1.0 mR per frame, however this value may be refined after further clinical testing. If at
process block 148, the current exposure rate is within a “half fine-tune range” of the reference exposure rate, then the program proceeds to process block 150, a fine tuning process block, and the amperage provided to the x-ray tube are adjusted in accordance to the disparity between the amperage and reference exposure rate. That is, if the current exposure is greater than the reference exposure rate, the amperage to the x-ray tube is reduced. The new value of amperage is compared against a predetermined range of amperage values (maximum beam current and minimum beam current values) so that the amperage value may never vary outside of this range. - If at
decision block 148, the current exposure rate is outside of the half fine tune range established atdecision block 148, a more substantial adjustment process is undertaken. Generally, the exposure provided by an x-ray system will follow the following equation: - X≈s·mA·kVp n. (8)
- where:
- s is seconds,
- mA is the amperage provided to the x-ray tube,
- kVp is the voltage provided to the x-ray tube, and
- n is a power factor dependent on the geometry of the machine and the particular kind of object being imaged.
- Generally, the value of n will not be known in advance. Accordingly in the more substantial correction process, n is deduced by obtaining two different exposures for equal predetermined intervals with different kVp values so that the value of n may be deduced.
- At
decision block 152, it is determined whether a first or second reference exposure is to be obtained. If the first reference exposure was just obtained, the program proceeds to process block 154 and a new value of kVp is determined for a second exposure. In this case, the first exposure used will be that which was employed to produce thehistogram 122 as previously described. - If the comparison of process block148 indicated that the exposure rate was too high, a lower kVp value is selected; and conversely, if the exposure at process block 148 indicated the exposure was too low, an increased value of kVp is provided. The new kVp value for the second exposure must be within a predetermined range of kVp values established by the user. Mathematically, the kVp value selected may be described as:
- kVp 2 =kVp 1 +a(dkVp) (9)
- where a is a step factor and
- dkVp is a minimum practical change in tube voltage.
- Two preferred means of selecting may be used: one providing linear and one providing logarithmic scaling. Such scaling techniques are well understood to those of ordinary skill in the art.
- If at
decision block 152, a second frame has already been taken with the new voltage value, then the program proceeds to process block 156 and the value of n in equation (9) is calculated. If the value of amperage is held constant between the first and second frame, the value of n may be determined according to the following equation: - where X1 and X2 are the measured exposure rates at the first and second frames, respectively and
- kVp1 and kVp2 are the two x-ray tube voltages during the first and second frames.
- At
process block 158, this value of ‘n’ is checked against threshold values intended to detect whether an erroneous value of n has been produced as a result of ‘clipping’ in the radiation data used to calculate exposure. As is understood in the art, clipping occurs when an increased dose of an element of the CCD camera produces no increase in the camera's output. - At
decision block 158, if the value of n calculated atprocess block 156 is greater than or equal to one, it is assumed to be valid and the program proceeds to process block 160 where kVp and mA are adjusted by setting mA equal to a maximum reference value and calculating kVp according to the following equation: - where kVPnew and mAnew are the settings for the next frame to be shot.
-
- If the value of n in
decision block 158 is less than one, then atprocess block 162, n is tested to see if it is less than zero. This value of n is realized when the exposure rate of the second frame changes in the opposite direction of the tube voltage. This suggests a clipped histogram and therefore the program branches back to process block 154 to obtain a new second frame. This condition may also arrive from object motion between the first and second frame. - On the other hand, if at
decision block 162, n is not less than zero (e.g. n is between zero and 1), the program proceeds to process block 166. Here it is assumed that because the sensitivity of the exposure rate on change in kVp is low, there may be some partial clipping. New values of kVp and mA are then computed and used with the previous second frame values to calculate a new n as follows. Generally, if kVp and mA are high, they are both lowered and if kVp and mA are low, they are both raised. - Referring now to FIG. 1, the image produced by the present invention may be used for quantitative analysis including, for example, that of making a bone density measurement. It is known to make bone density analyses from x-ray images through the use of dual energy techniques in which the voltage across the x-ray tube is changed or a filter is periodically placed within the x-ray beam to change the spectrum of the x-ray energy between two images. The two images may be mathematically processed to yield information about different basis materials within the image object (e.g. bone and soft tissue). For these quantitative measurements, it is desirable to eliminate the effect of scatter.
- Referring now to FIG. 23 in imaging a patient's
spine 200, for example,x-rays 202 are directed from anx-ray source 201 through thepatient 199 to pass throughsoft tissue 204 surrounding aspine 200. Certain of thex-rays 202 are blocked by thespine 200 and others pass through thespine 200 to be recorded at theimage intensifier 206. Anattenuation image 208 measured by an image intensifier measures those x-rays passing through the patient 109. - A
portion 210 of the attenuation image directly beneath thespine 200 records not only thosex-rays 202 passing through thespine 200 and thesoft tissue 204 above and below it, but also scatteredx-rays 212 directed, for example, throughsoft tissue 204 to the side of thespine 200 but then scattered by the soft tissue to proceed at an angle to theportion 210 of theattenuation image 208 beneath thespine 200. Because thescattered x-rays 212 do not carry information about the attenuation of thespine 200, they are desirably removed from theimage 208 prior to its use in quantitative measurement. - For this purpose, the present invention uses an
occluder 214 being an x-ray transparent plate such as may be constructed of Plexiglas and incorporating into its body, a plurality of x-ray blocking lead pins 216. Preferably these pins are placed so as to projectimages 218 onto theimage 208 received by theimage intensifier 206 in positions outside animage 220 of thespine 200. Generally therefore, thepins 216 are placed at the periphery of theoccluder 214. Thepins 216 are sized so as to substantially block all direct x-rays from passing through them but so that theirimages 218 include a significant portion ofscattered x-rays 212. - Referring now to FIG. 22 at a first step of a scatter reduction operation with the
occluder 214 of FIG. 21, an image is acquired of the imaged object, for example, thespine 200 and its surrounding soft tissue 204 (not shown in FIG. 21) including theimages 218 of thepins 216. This acquisition is indicated by process block 221 of FIG. 22. - The
pins 216 are held in predetermined locations with respect to theimage 208 so that theirimages 218 may be readily and automatically identified. Preferably thepins 216 are placed at the interstices of a Cartesian grid, however, other regular patterns may be chosen. Theimage 208 may be corrected for pincushion type distortion, as described above, so that the locations of thepins 216 may be readily located in the image based on their known positions in theoccluder 214. - At each
pin image 218, avalue 222 indicating the magnitude of the received x-rays, shown in FIG. 23, may be ascertained. Thisvalue 222 measures the scatter received in the vicinity ofimage 218 caused generally by the effect of thesoft tissue 204 and possible secondary scatter effects in theimage intensifier 206.Values 222 are recorded, as indicated byprocess block 224, for eachpin image 218. From these values, a set of normalizing points are established. - The
image 208 is then used to derive a scatter map. Referring to FIG. 23, generally the amount of scatter at a given point will be a function of how many x-ray photons are received at points adjacent to the given point. For example, comparing theimage 208 to atheoretical scatterless image 228 generally in anattenuated region 230 of the image 208 (e.g., under the spine 200), scatter will increase the apparent value in theimage 208 as a result of radiation from nearby low attenuation regions scattering into thehigh attenuation region 230. Conversely the apparent value at alow attenuation region 232 will be decreased because of the scatter into the high attenuation region. - A map of the scattered radiation may thus be modeled by “blurring” the
image 208. This blurring can be accomplished by a low pass filtering of theimage 208, i.e., convolving theimage 208 with a convolution kernel having rectangular dimensions corresponding to the desired low pass frequency cut off. The effect is an averaging of theimage 208 producingscatter map 234. - The image used to produce the
scatter map 234 is anattenuation image 208 obtained from thepatient 199 without theoccluder 214 in place, or may be animage 208 including theimages 218 of thepins 216 but with thelatter images 218 removed based on knowledge of their location. This removal ofimages 218 may substitute values of theimage 208 atpoints 239 on either side of theimages 218. The process of driving the scatter map from the image is indicated by process block 235 of FIG. 24. - Next as indicated by
process block 237, thescatter map 234 is fit to the normalizingpoints 222 previously determined atprocess block 224. - Referring to FIG. 24, the
scatter map 234 is thus normalized so that theportions 238 of thescatter map 236 located near the places where theimages 218 would fall are givenvalues 222 as determined atprocess block 224. This involves a simple shifting up or down of thescatter map 236 and may employ a “least square” fit to shift thescatter map 236 tomultiple values 222 obtained from eachpin 216. As adjusted, thescatter map 236 is then subtracted from theimage 208 to eliminate or reduce the scatter in that image as indicated byprocess block 239. - The effect of subtracting a low pass filtered or blurred image properly normalized to actual scatter is to sharpen up the
image 208 but also to preserve its quantitative accuracy. Thus the present invention differs from prior art scatter reduction techniques in that it both addresses the variation in scatter across the image caused by attenuation of x-rays by the imaged object but also incorporates accurate measurements of scatter in certain portions of the image. - Referring now to FIG. 25, a
mobile fluoroscopy machine 310 suitable for use with the present invention is similar to that which has been described above with respect to FIG. 1 with exceptions that will be apparent from context. - The
mobile fluoroscopy machine 310 includes amobile cart 312 supporting acomputer 314 and monitor andkeyboard 317 for receiving and processing digital x-ray image data. Thecart 312 supports on one side an articulatingarm assembly 316 terminating in a rotatable C-arm 318. The C-arm supports, at the ends of the C, animage intensifier 320 and anx-ray source 322 opposed along anaxis 324 so that thex-ray source 322 projects a cone-beam of x-ray radiation toward theimage intensifier 320 alongaxis 324. - The articulating
arm assembly 316 is connected to the C-arm 318 through one or morepivotal links 327 so that theaxis 324 may be positioned to be horizontal approximately two feet above the floor to rest upon or be supported against the upper end of a supportingpedestal 326 or may be attached to thecart 312. Referring also to FIG. 26, thepedestal 326 includes a hemicylindricallyconcave cradle 328 at its upper surface to receive a lower portion of thecylindrical image intensifier 320 when the C-arm is so positioned to rest against thepedestal 326. - The
pedestal 326 also provides on its upper surface achannel 330 extending across theaxis 324 between theimage intensifier 320 and thex-ray source 322 when the latter are positioned on thepedestal 326. Thechannel 330 may receive alimb positioner 332 such as may be adapted to support a patient's foot or arm across theaxis 324 for densiometric measurement. Thepedestal 326 may be weighted so as to provide a stable surface for support of thex-ray source 322 andimage intensifier 320 and to provide adequate support for the patient's limb. The height of thepedestal 326 is selected to be suitable for either arm or foot imaging. - Referring now to FIG. 27 and30, the
channel 330, extending substantially perpendicularly toaxis 324 and has ahorizontal bottom surface 333 pierced by two vertically extendingguide holes 334 which may be used to receive and position corresponding pins 337 on one of twolimb positioners 332. Afoot positioner 336, as shown in FIG. 26, provides a paddedcalf support plate 338 fitting adjacent to thebottom surface 333 and an upwardly extendingsole support 340 forming an obtuse angle with respect to thecalf support plate 338. Acushion 342 on thecalf support plate 338 may be adjusted so as to allow the patient's leg to extend upward somewhat from vertical for comfort. Gussets 344 span the angle between thesole plate 340 andcalf support plate 338 to fix them in relative position but includeapertures 346 to allow for the free passage of x-rays through a portion alongaxis 324 where the os calcis of the heal will be located. - When positioned within the
channel 330, thefoot positioner 336 is also supported by upwardly extendingchannel sidewalls 348 which serve further to provide an alignment surface for the imaging face of theimage intensifier 320 or other detector array and on the other side, an alignment surface for an emitting face of thex-ray source 322.Channel sidewalls 348 are generally radio translucent so as to permit the passage of x-rays therethrough, but may include: calibration materials such as are well known in the art for calibrating dual energy devices, antiscatter grids also well known in the art, or occluders for evaluating scatter as have been described above or in the parent applications hereby incorporated by reference. - Referring to FIG. 27 and28, for forearm imaging, the
foot positioner 336 is removed and apalm support 352 is inserted by means ofpin 337 in one of theholes 334 so as to locate a user's arm resting against thebottom surface 333 with the user's palm against thepalm support 352 such that the bones of the forearm are placed along theaxis 324 for imaging. - Referring to FIG. 28 and29, the
hemispherical support cradle 328 may include three radially inwardly extendingribs 354 attached by means of screws or the like to be replaceable. Two of theribs 354 are positioned in a horizontal plane to substantially bisect theimage intensifier 320 when it is placed within thecradle 328. Thethird rib 354 is positioned at the bottom of thecradle 328 and is opposed by arotating locking collar 358 which may be used to further secure theimage intensifier 320 within thecradle 328. The front edge of the image intensifier is abutted against the upright face of the dividing barrier so as to precisely locate it alongaxis 324. The inner edges of theseribs 354 define aninner radius 356 of lesser diameter than thecradle 328 that by proper design of theribs 354 may be adjusted to conform to the outer surface of aparticular image intensifier 320. - Referring now again to FIG. 30, some fluoroscopy equipment will not permit digital imaging or the necessary dual energy control needed for densitometry. Accordingly, an
independent detector array 360 may be placed within thecradle 328 in lieu of theimage intensifier 320. Thisdetector array 360 may be a pair of stimulable phosphor plates as are understood in the art with intermediate filtering so as to provide dual energy readings with a polychromatic x-ray source. In this way a switching of voltage on thex-ray source 322, as described above, can be avoided. Alternatively, thedetector array 360 may be a large area solid state detector or scanning detector assembly such as are understood in the art including those constructed of amorphous silicon and thin film transistor technology or those employing active pixel technology in which C-MOS integrated circuit fabrication techniques are employed. These detectors may be used with a switchedx-ray source 322 to provide dual energy imaging or may be used in a stacked configuration with intermediate filtering so as to provide separate energy measurements, or may be used in a side-by-side configuration with interleaved detector elements filtered so as to be selectively sensitive to different energies. - In the preferred embodiment, and as shown in FIG. 31, the
independent detector array 360 is a “pancake”image intensifier 361, suitably small so as to fit within the space between aconventional image intensifier 320 and thex-ray source 322. Referring to both FIG. 30 and 31, thepancake image intensifier 361 includes avacuum bowl 362 having a planarfront surface 364 for receiving x-rays 366 (normally through the channel sidewalls 348 of the stand 326). - According to conventional design, the
x-rays 366 pass through thefront surface 364 of thevacuum bowl 362 to strike atarget material 368 to ejectelectrons 370 into the volume of thebowl 362. Focusingelectrodes 372 direct the electrons to aphosphor 374 where an image is formed to be received byimaging array 375 such as a CCD array or camera. The image area of thephosphor 374 is much smaller than thefront surface 364 so as to reduce the image size to one compatible with the camera. In the present invention the distance B between thetarget 368 and the imaging array 375 (including any optical path through one or more focusing lenses) is less than or equal to the radial dimension. A of thefront surface 364 gives thepancake image intensifier 361 an extremely short form factor suitable for practice with the present invention. - Hitherto, such form factors were avoided because they are known to result in severe distortion of the image formed on the
phosphor 374. This distortion is accommodated in the present invention by means of digital image processing incomputer 314 which receives digitized pixel data from scanningelectronics 378 connected to theimaging array 375 and corrects it according to the correction process described above with respect to the pin cushion correction. Accordingly, the addition of digital signal processing allows for production ofpancake image intensifier 361 in which the separation of the imaging optics from the front of the image intensifier is much reduced. - The above adapter may be modified to use in femur imaging. In this case the
pedestal 326 may be eliminated in favor of a positioner (not shown) attached to theimage intensifier 320 orx-ray source 322 directly. In the former case, the positioner may provide for a fixed air gap between the patient and theimage intensifier 320 to reduce received scatter. So as to allow free manipulation of the C-arm 318, the positioner may be a lightweight plastic radiolucent material and may optionally include a calibration system such as a flip in phantom for calibration of the dual energy readings and occluders for scatter correction as has been described above. Collimation and/or a separate solid state dual energy image detectors may also be held by the positioner whose outer surface may guide the positioning of theC arm 318 to the necessary orientation which need not be horizontal but may be vertical for fore arm measurements or the like. For femur measurements, the patient may stand and the C-arm 318 manipulated appropriately as guided by the positioner. - Referring now to FIG. 31 and32,
computer 314 includes aprocessor 380 andmemory 382, the latter of which receives raw image data in the form of pixels having spatial locations and brightnessvalues forming images 384.Memory 382 also includes aprocessing program 386 providing a general interface and control of the operation of thefluoroscopy machine 310 and a processing ofimages 384 so as to provide a quantitative measure of bone isolated from soft tissue. - The
processing program 386 can be simply loaded into thecomputer 314 for thefluoroscopy machine 310 when thepedestal 326 is to be employed with afluoroscopy machine 310 providing digital imaging and x-ray voltage control. If anindependent detector array 360 is required, theprogram 386 may be executed on acomputer 314 associated with thatindependent detector array 360. - At a first step in the
program 386, indicated byprocess block 388, the operator of thefluoroscopy machine 310, having indicated a desire to perform densitometry and having positioned the C-arm in thepedestal 326, enters patient data that will be used to identify theimage 384 to be collected. - At succeeding
process block 390, data is collected for threedistinct images 384 with: 1) no x-ray exposure, 2) high energy x-ray exposure, and 3) low energy x-ray exposure. Each of the exposures is preserved as a separate image file in the memory of thecomputer 314. The first exposure is used for correction routines to be described; the latter two exposures are used to deduce bone density according to methods well known in the art in which variations in high energy and low energy absorption are used to deduce the Compton scattering and atomic number of the material lying between thex-ray source 322 and theimage intensifier 320. As is understood in the art, these two measurements allow the amount of bone as opposed to soft tissue located in that image region to be accurately measured. The data is acquired directly from theindependent detector array 360 or in the event that stimulable plates are used, a reader may be attached to thecomputer 314 so as to acquire the necessary pixel data of animage 384. In the same way a conventional photographic film/filter plate arrangement may be used. - At
next process block 393, each of these images is corrected for non-linearity of the detector such as may be determined empirically at an earlier time by testing the detector according to methods well known in the art, and which is a function of the detector and the technology used by the detector. Generally the testing exposes the detector to different fluences of x-rays and measures the output of the detector and the correction is intended to ensure that, for example, a doubling of fluence results in a doubling of detector output after correction. The correction is generally simply a scaling of each of the images by a factor that is a function of the pixel value for each pixel and possibly the location of the pixel. - At
process block 395, noise related to the particular line of the detector is removed. Referring to FIG. 34, theimaging array 375 provides a matrix ofdetector elements 392 arranged in rows and columns. Normally, either rows or columns are ganged together to be read out by dedicated read outelectronics 391 spanning a particular row or column. The read out electronics introduces noise which is imposed upon eachdetector element 392 of that row and which is thus line correlated, that is, more highly correlated withother detector elements 392 of the line thandetector elements 392 of different lines. To eliminate line correlated noise, onedetector element 392 in each line is blocked by alead mask 394 so as to be shielded from x-rays. Apixel value 396 from this blockeddetector element 392 will provide a value that varies according to the line noise thus a line correlatednoise value 398 may be deduced and subtracted from the pixel values 400 of theother detector elements 392 in the line. - Referring again to FIG. 32, at a succeeding
process block 402, veiling glare is removed and the field is flattened. This former correction attempts to eliminate blurring of the image such as may be caused by scatter or similar effects within theimaging array 375. Glare refers generally to a reading that would be obtained underdetector elements 392 that were wholly shadowed by an occluding absorber on the surface of thearray 375. The glare is a function of the detector technology and is reduced by a deconvolution process based on an empirically derived deconvolution kernel according to a number of techniques well known in the art. - Also at process block402 the field is flattened which is to say the gain variation of the
detector elements 392 are normalized according to an empirically derived normalization map determined at the factory by exposing the detector to a uniform x-ray elimination and noting variation and intensities reported in the pixel values 400. At this time, dark currents from thedetector elements 392 may also be eliminated as determined from the no-exposure x-ray image taken atprocess block 390. - Referring now to process block404, a dynamic scatter correction may be employed as has been previously described with respect to FIGS. 21-23. Alternatively referring also to FIG. 33, a dynamic scatter correction may be employed in which the data of the
image 384 is analyzed so as to create ahistogram 406 ofpixel values 400 for the entire image. The histogram may be divided into regions 408 (five equal regions in the preferred embodiment) corresponding roughly x-ray paths through: 1) air-only, 2) thin tissue, 3) thick tissue only, 4) thin bone and, and 5) thick-bone. Each of these materials will exhibit a different scattering and hence a different empirically derivedscatter kernel 410 may be assigned to eachregion 408 with generally the lower density regions having narrower kernels commensurate with less scatter. - The selected
scatter kernel 410 may be scaled by the pixel value of theimage 384 on a pixel by pixel basis and that kernel, so scaled, applied to adeconvolver 412 used to deconvolve theimage 384 to produce adeconvolved image 414. A number of techniques of deconvolution are well known in the art using a fixed scatter kernel and these same techniques may be used with thevariable scatter kernel 410 described here. During deconvolution thekernel 410 will be sequentially applied to a set of adjacent pixel values determined by the width of the kernel. The center pixel value at any step of the deconvolution will be used to scale the kernel and to identify the region of the histogram for the purpose of selecting thekernel 410. In an alternative embodiment, the kernel may be fixed and simply scaled by the value of the centermost pixel during de-convolution. - Referring again to FIG. 32, at
process block 416, theimages 384 are log corrected reflecting the fact that attenuation is exponentially related to thickness. The images are now related to thickness, a dimension which will be important in the ultimate bone-density determination. - At following
process block 418, speckle may be identified forcertain x-ray detectors 360 that are subject to extremely high readout values caused by noise which is possibly related to direct x-ray irradiation of the detector element. Speckle is identified by a simple thresholding process. - At
next process block 420, path length correction may be performed based on the geometry of the particular C-arm such as may vary path length and magnification across the image as is well understood in the art. - Similarly at succeeding
process block 422, beam hardening, the well known effect of a spectral shift in a polyenergetic x-ray beam as it passes through different thicknesses of material, and a Heel effect correction may be made, the Heel effect correction referring to a variation in the spectrum of an x-ray beam as a function of its angle in the cone of x-ray beams. Both of the corrections are known in the art, but must be employed in the present invention in order to provide suitable quantitative accuracy for densitometry. - At
process block 424, the identified speckle ofprocess block 418 is corrected by eliminating these identified pixels from subsequent calculation or by replacing them with a local average value. - The entire image may then be averaged or low-passed filtered at process block426 so as to further reduce noise and to eliminate unneeded resolution.
- The images are then processed according to well understood techniques to produce a bone mineral density value at
process block 428. This bone mineral density value indicates the amount of bone material at each pixel of the image largely independent of surrounding soft tissue. The pixel image may be analyzed in a number of methods but most simply, as indicated byprocess block 430, by defining either automatically or manually a desired region of interest within the image and making a measurement of total bone density within that region. Automated techniques may look for a local maximum or minimum of bone density or may use image recognition type techniques to locate reproducibly a particular region of the forearm or os calcis. Morphometric analysis may be applied to the image to detect bone fracture and other techniques such as texture analysis may be performed according to methods well known in the art. The results of the analyses and images so processed may be displayed by thecomputer 314. - It is thus envisioned that the present invention is subject to many modifications which will become apparent to those of ordinary skill in the art. Accordingly, it is intended that the present invention not be limited to the particular embodiment illustrated herein, but embraces all such modified forms thereof as come within the scope of the following claim.
Claims (27)
1. A positioner for a fluoroscopy machine of a type having a radiation source and detector separated along a beam axis and mounted on an articulated arm, the adapter comprising:
a weighted base supporting the positioner upon a floor;
a pedestal extending upward from the base to provide receiving supports for at least one of the radiation source and detector so that when the at least one of the radiation source and detector is received by the receiving supports, the beam axis is in predetermined orientation with respect to the top of the pedestal; and
a limb positioner attached to the pedestal between the receiving supports so that a patient's limb held by the limb positioner is intersected by the beam axis.
2. The positioner as recited in wherein the limb positioner is removable and is sized to receive the patent's limb, wherein the limb is an arm, a hand, a leg or a foot.
claim 1
3. The positioner as recited in wherein the limb positioner is a vertically extending palm support.
claim 1
4. The positioner as recited in wherein the palm support is attached to a mounting pin sized and positioned to fit within a guide bore in the pedestal.
claim 3
5. The positioner as recited in wherein the palm support is removable.
claim 4
6. The positioner as recited in wherein the limb positioner is a foot cradle having:
claim 1
a calf support plate adjacent to the top of the pedestal having at least one mounting pin depending downward sized and positioned to fit within a guide bore in the pedestal;
a sole support plate joined to the calf support plate to define an obtuse angle; and
side gussets spanning the calf support plate and the sole support plate, the gussets having apertures positioned to be at the beam axis when the foot cradle is mounted to the pedestal.
7. The positioner as recited in wherein the foot cradle includes padding material at top surfaces of the calf support plates.
claim 6
8. The positioner as recited in wherein the foot cradle is removable.
claim 6
9. The positioner as recited in including a calibration material attached to the pedestal within the beam axis.
claim 1
10. The positioner as recited in including an anti-scatter grid attached to the pedestal within the beam axis.
claim 1
11. The positioner as recited in including an occlude attached to the pedestal within the beam axis.
claim 1
12. The positioner as recited in wherein the receiving supports are disposed axially along a radius of a cavity, the receiving supports extending into the cavity to define a radius of lesser diameter than the cavity according to the size of the radiation detector.
claim 1
13. The positioner as recited in further comprising an index guide adjacent to the limb positioner for properly aligning the radiation detector with the limb positioner on the pedestal.
claim 1
14. The positioner as recited in , further comprising:
claim 1
an independent detector array located along the beam axis between the limb positioner and the radiation detector when the radiation detector is disposed in the receiving supports, the detector array producing attenuated dual energy signals; and
a processor receiving the attenuated dual energy signals and calculating bone density measurements.
15. The positioner as recited in wherein the detector array is a set of stimulable phosphor plates.
claim 14
16. The positioner as recited in wherein the detector array is solid state detector.
claim 14
17. The positioner as recited in wherein the detector array is a stacked linear array scanning detector.
claim 14
18. A method of scatter correction for x-ray measurements comprising the steps of:
(a) obtaining an x-ray measurement providing a plurality of x-ray attenuation values over an area;
(b) preparing a histogram of the x-ray attenuation values, the histogram indicating the frequency of occurrence of attenuation values over a range of attenuations;
(c) dividing the attenuation range into at least two regions corresponding to material types;
(d) associating a different deconvolution kernel with each region; and
(e) deconvolving the x-ray measurement with the different deconvolution kernels depending on the regions of the deconvolved attenuation values.
19. The method of wherein the deconvolution kernels define a width of deconvolution measured as a number adjacent x-ray measurements within the deconvolution and wherein the different deconvolution kernels have different widths.
claim 18
20. The method of wherein the regions corresponding to greater attenuations, and hence to imaged materials having greater attenuation, are associated with deconvolution kernels having greater widths.
claim 19
21. The method of wherein the x-ray measurement is of a human and wherein the attenuation range is divided into regions corresponding to materials selected from the group consisting of: air-only, thin tissue, thick tissue only, thin bone and, and thick-only.
claim 18
22. The method of wherein the x-ray measurement is of a human and wherein the deconvolved x-ray measurement is further analyzed to provide a quantitative measurement of an imaged object.
claim 18
23. The method of wherein the quantitative measurement is an indication of bone density.
claim 22
24. The method of wherein during step (e) the deconvolution kernels are further weighted according to the intensity of the deconvolved attenuation values.
claim 18
25. A method of scatter correction for x-ray measurements comprising the steps of:
(a) obtaining an x-ray measurement providing a plurality of x-ray attenuation measurements over an area;
(b) preparing at least one deconvolution kernel for scatter reduction in the x-ray measurements;
(c) during the application of the deconvolution kernel to the x-ray measurements, applying a scale factor to the deconvolution kernel based on the value of the x-ray measurements being deconvolved.
26. The method of wherein the x-ray measurement is of a human and wherein the deconvolved x-ray measurement is further analyzed to provide a quantitative measurement of an imaged object.
claim 25
27. The method of wherein the quantitative measurement is an indication of bone density.
claim 26
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/745,704 US6315445B1 (en) | 1996-02-21 | 2000-12-21 | Densitometry adapter for compact x-ray fluoroscopy machine |
Applications Claiming Priority (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US1199396P | 1996-02-21 | 1996-02-21 | |
PCT/US1997/002770 WO1997031335A1 (en) | 1996-02-21 | 1997-02-21 | X-ray imaging system |
US09/006,358 US6007243A (en) | 1996-02-21 | 1998-01-13 | Combined mobile x-ray imaging system and monitor cart |
US8016498P | 1998-03-13 | 1998-03-13 | |
US09/281,518 US6215846B1 (en) | 1996-02-21 | 1999-03-30 | Densitometry adapter for compact x-ray fluoroscopy machine |
US09/745,704 US6315445B1 (en) | 1996-02-21 | 2000-12-21 | Densitometry adapter for compact x-ray fluoroscopy machine |
Related Parent Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US09/006,358 Division US6007243A (en) | 1996-02-21 | 1998-01-13 | Combined mobile x-ray imaging system and monitor cart |
US09/281,518 Division US6215846B1 (en) | 1996-02-21 | 1999-03-30 | Densitometry adapter for compact x-ray fluoroscopy machine |
Publications (2)
Publication Number | Publication Date |
---|---|
US6315445B1 US6315445B1 (en) | 2001-11-13 |
US20010053202A1 true US20010053202A1 (en) | 2001-12-20 |
Family
ID=27617734
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US09/745,704 Expired - Fee Related US6315445B1 (en) | 1996-02-21 | 2000-12-21 | Densitometry adapter for compact x-ray fluoroscopy machine |
Country Status (1)
Country | Link |
---|---|
US (1) | US6315445B1 (en) |
Cited By (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080279328A1 (en) * | 2005-11-18 | 2008-11-13 | Koninklijke Philips Electronics N.V. | Systems and Methods Using X-Ray Tube Spectra For Computed Tomography Applications |
US20100046705A1 (en) * | 2007-03-30 | 2010-02-25 | General Electric Company | Image acquisition and processing chain for dual-energy radiography using a portable flat panel detector |
US20100158341A1 (en) * | 2008-12-18 | 2010-06-24 | Siemens Medical Solutions Usa, Inc. | System for Automatically Generating a Mask for Digital Subtraction Angiography |
US8660226B2 (en) * | 2011-01-19 | 2014-02-25 | General Electric Company | Systems and methods for multichannel noise reduction |
CN103767723A (en) * | 2012-10-25 | 2014-05-07 | 南京普爱射线影像设备有限公司 | Cone beam CT three-dimensional digital imaging method based on C-shaped arm |
US20150036910A1 (en) * | 2010-03-23 | 2015-02-05 | Hologic, Inc. | Measuring intramuscular fat |
US20160140708A1 (en) * | 2013-07-23 | 2016-05-19 | Fujifilm Corporation | Radiation-image processing device and method |
US20170026592A1 (en) * | 2015-07-24 | 2017-01-26 | Lytro, Inc. | Automatic lens flare detection and correction for light-field images |
US20170098315A1 (en) * | 2015-10-06 | 2017-04-06 | Rigaku Corporation | Analyzer, analysis method and analysis program of bone mineral density |
WO2017192597A1 (en) * | 2016-05-03 | 2017-11-09 | Rapiscan Systems, Inc. | Radiation signal processing system |
US10275892B2 (en) | 2016-06-09 | 2019-04-30 | Google Llc | Multi-view scene segmentation and propagation |
US10275898B1 (en) | 2015-04-15 | 2019-04-30 | Google Llc | Wedge-based light-field video capture |
US10298834B2 (en) | 2006-12-01 | 2019-05-21 | Google Llc | Video refocusing |
US10334151B2 (en) | 2013-04-22 | 2019-06-25 | Google Llc | Phase detection autofocus using subaperture images |
US10341632B2 (en) | 2015-04-15 | 2019-07-02 | Google Llc. | Spatial random access enabled video system with a three-dimensional viewing volume |
US10354399B2 (en) | 2017-05-25 | 2019-07-16 | Google Llc | Multi-view back-projection to a light-field |
US10412373B2 (en) | 2015-04-15 | 2019-09-10 | Google Llc | Image capture for virtual reality displays |
US10419737B2 (en) | 2015-04-15 | 2019-09-17 | Google Llc | Data structures and delivery methods for expediting virtual reality playback |
US10440407B2 (en) | 2017-05-09 | 2019-10-08 | Google Llc | Adaptive control for immersive experience delivery |
US10444931B2 (en) | 2017-05-09 | 2019-10-15 | Google Llc | Vantage generation and interactive playback |
US10469873B2 (en) | 2015-04-15 | 2019-11-05 | Google Llc | Encoding and decoding virtual reality video |
US10470705B2 (en) | 2004-10-04 | 2019-11-12 | Hologic, Inc. | Estimating visceral fat by dual-energy x-ray absorptiometry |
US10474227B2 (en) | 2017-05-09 | 2019-11-12 | Google Llc | Generation of virtual reality with 6 degrees of freedom from limited viewer data |
US10540818B2 (en) | 2015-04-15 | 2020-01-21 | Google Llc | Stereo image generation and interactive playback |
US10545215B2 (en) | 2017-09-13 | 2020-01-28 | Google Llc | 4D camera tracking and optical stabilization |
US10546424B2 (en) | 2015-04-15 | 2020-01-28 | Google Llc | Layered content delivery for virtual and augmented reality experiences |
US10552947B2 (en) | 2012-06-26 | 2020-02-04 | Google Llc | Depth-based image blurring |
US10565734B2 (en) | 2015-04-15 | 2020-02-18 | Google Llc | Video capture, processing, calibration, computational fiber artifact removal, and light-field pipeline |
US10567464B2 (en) | 2015-04-15 | 2020-02-18 | Google Llc | Video compression with adaptive view-dependent lighting removal |
US10594945B2 (en) | 2017-04-03 | 2020-03-17 | Google Llc | Generating dolly zoom effect using light field image data |
US10646159B2 (en) | 2004-10-04 | 2020-05-12 | Hologic, Inc. | Visceral fat measurement |
US10679361B2 (en) | 2016-12-05 | 2020-06-09 | Google Llc | Multi-view rotoscope contour propagation |
US10965862B2 (en) | 2018-01-18 | 2021-03-30 | Google Llc | Multi-camera navigation interface |
US11328446B2 (en) | 2015-04-15 | 2022-05-10 | Google Llc | Combining light-field data with active depth data for depth map generation |
US11701079B2 (en) | 2013-11-25 | 2023-07-18 | Hologic, Inc. | Bone densitometer |
US11717244B2 (en) | 2015-02-26 | 2023-08-08 | Hologic, Inc. | Methods for physiological state determination in body scans |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100419344B1 (en) * | 2001-09-10 | 2004-02-19 | 주식회사 오스테오시스 | Apparatus for measuring density of bone |
US6754306B2 (en) * | 2001-11-21 | 2004-06-22 | The United States Of America As Represented By The Secretary Of The Army | Portable medical digital radiography assembly |
US6648825B1 (en) * | 2002-05-23 | 2003-11-18 | Koninklijke Philips Electronics N.V. | Diagnostic ultrasound system cart with swiveling control panel |
US7215806B1 (en) * | 2003-05-01 | 2007-05-08 | Analogic Corporation | Adaptive iterative bone correction in computerized tomography |
US6968042B2 (en) * | 2003-09-12 | 2005-11-22 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for target angle heel effect compensation |
US7600294B2 (en) * | 2006-09-19 | 2009-10-13 | Dell Products L.P. | Chassis foot and optional castor assembly |
FI125531B (en) * | 2010-04-29 | 2015-11-13 | Planmed Oy | Medical x-ray equipment |
FI123379B (en) * | 2010-04-29 | 2013-03-15 | Planmed Oy | Medical computed tomography equipment |
CN103747734B (en) * | 2011-07-19 | 2016-03-09 | 株式会社日立医疗器械 | The control method of X-ray imaging apparatus and X-ray generator |
US9044186B2 (en) | 2012-06-25 | 2015-06-02 | George W. Ma | Portable dual-energy radiographic X-ray perihpheral bone density and imaging systems and methods |
EP2749666B1 (en) * | 2012-12-31 | 2019-10-09 | Aktiebolaget SKF | Thermo-mechanical process for martensitic bearing steels and steel bearing component |
US9275769B2 (en) * | 2013-03-14 | 2016-03-01 | Pcc Structurals, Inc. | Marking template for radiography |
JP6283875B2 (en) * | 2013-09-05 | 2018-02-28 | キヤノンメディカルシステムズ株式会社 | Medical image processing apparatus, X-ray diagnostic apparatus, and X-ray computed tomography apparatus |
KR101370496B1 (en) * | 2013-09-26 | 2014-03-06 | 한국건설기술연구원 | Estimation method for volume fractions of constituent materials in a unit voxel of x-ray ct image |
CN104970814A (en) * | 2014-04-04 | 2015-10-14 | 株式会社东芝 | C-shaped arm device, C-shaped arm retraining device and X-ray diagnosis device |
US10448910B2 (en) * | 2016-02-03 | 2019-10-22 | Globus Medical, Inc. | Portable medical imaging system |
US10827971B2 (en) * | 2017-12-20 | 2020-11-10 | Howmedica Osteonics Corp. | Virtual ligament balancing |
CN116485686B (en) * | 2023-06-19 | 2023-08-29 | 青岛国源中创电气自动化工程有限公司 | Sewage treatment image enhancement method by activated sludge process |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3639764A (en) * | 1969-04-02 | 1972-02-01 | Packard Instrument Co Inc | Bone-density measuring instrument |
US4829549A (en) | 1985-06-19 | 1989-05-09 | Vogel John M | Densitometer for scanning os calcis for predicting osteoporosis |
US5042489A (en) | 1988-05-11 | 1991-08-27 | Lunar Corporation | Ultrasonic densitometer device and method |
FR2633793A1 (en) | 1988-07-01 | 1990-01-05 | Gen Electric Cgr | CORRECTION METHOD FOR DISTORTION OF RADIOLOGICAL IMAGES |
US5136743A (en) * | 1991-07-22 | 1992-08-11 | The United States Of America As Represented By The Secretary Of The Air Force | Distal radioulnar joint stress platform |
JPH06153077A (en) | 1992-11-10 | 1994-05-31 | Hitachi Medical Corp | X-ray television camera device |
BE1007766A3 (en) | 1993-11-10 | 1995-10-17 | Philips Electronics Nv | Method and apparatus for computer tomography. |
US5479471A (en) * | 1994-04-22 | 1995-12-26 | Buckland; Peter E. | Foot holding and posing apparatus for X-ray examination |
US5878108A (en) | 1995-11-30 | 1999-03-02 | Hitachi Medical Corporation | Method for generating X-ray image and apparatus therefor |
US6215846B1 (en) * | 1996-02-21 | 2001-04-10 | Lunar Corporation | Densitometry adapter for compact x-ray fluoroscopy machine |
US5949846A (en) * | 1997-02-03 | 1999-09-07 | Hologic, Inc. | Bone densitometry using x-ray imaging systems |
FR2759800B1 (en) | 1997-02-17 | 1999-03-26 | Commissariat Energie Atomique | METHOD FOR CORRECTING THE DIFFUSED FLOW IN DIGITAL RADIOGRAPHY IMAGES |
US5748704A (en) * | 1997-03-10 | 1998-05-05 | Lunar Corporation | Peripheral bone densitometer |
US5903008A (en) | 1997-07-02 | 1999-05-11 | General Electric Company | Scatter correction methods and systems in single photon emission computed tomography |
US6061419A (en) | 1998-08-25 | 2000-05-09 | General Electric Company | Methods and apparatus for noise compensation in an imaging system |
-
2000
- 2000-12-21 US US09/745,704 patent/US6315445B1/en not_active Expired - Fee Related
Cited By (48)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10646159B2 (en) | 2004-10-04 | 2020-05-12 | Hologic, Inc. | Visceral fat measurement |
US10470705B2 (en) | 2004-10-04 | 2019-11-12 | Hologic, Inc. | Estimating visceral fat by dual-energy x-ray absorptiometry |
US20080279328A1 (en) * | 2005-11-18 | 2008-11-13 | Koninklijke Philips Electronics N.V. | Systems and Methods Using X-Ray Tube Spectra For Computed Tomography Applications |
US10298834B2 (en) | 2006-12-01 | 2019-05-21 | Google Llc | Video refocusing |
US20100046705A1 (en) * | 2007-03-30 | 2010-02-25 | General Electric Company | Image acquisition and processing chain for dual-energy radiography using a portable flat panel detector |
US7997799B2 (en) | 2007-03-30 | 2011-08-16 | General Electric Company | Image acquisition and processing chain for dual-energy radiography using a portable flat panel detector |
US20100158341A1 (en) * | 2008-12-18 | 2010-06-24 | Siemens Medical Solutions Usa, Inc. | System for Automatically Generating a Mask for Digital Subtraction Angiography |
US10515451B2 (en) * | 2010-03-23 | 2019-12-24 | Hologic, Inc. | Measuring intramuscular fat |
US11625824B2 (en) | 2010-03-23 | 2023-04-11 | Hologic, Inc. | Measuring intramuscular fat |
US9865050B2 (en) * | 2010-03-23 | 2018-01-09 | Hologic, Inc. | Measuring intramuscular fat |
US20180189948A1 (en) * | 2010-03-23 | 2018-07-05 | Hologic, Inc. | Measuring intramuscular fat |
US20150036910A1 (en) * | 2010-03-23 | 2015-02-05 | Hologic, Inc. | Measuring intramuscular fat |
US8660226B2 (en) * | 2011-01-19 | 2014-02-25 | General Electric Company | Systems and methods for multichannel noise reduction |
US10552947B2 (en) | 2012-06-26 | 2020-02-04 | Google Llc | Depth-based image blurring |
CN103767723A (en) * | 2012-10-25 | 2014-05-07 | 南京普爱射线影像设备有限公司 | Cone beam CT three-dimensional digital imaging method based on C-shaped arm |
US10334151B2 (en) | 2013-04-22 | 2019-06-25 | Google Llc | Phase detection autofocus using subaperture images |
US20160140708A1 (en) * | 2013-07-23 | 2016-05-19 | Fujifilm Corporation | Radiation-image processing device and method |
US10475180B2 (en) * | 2013-07-23 | 2019-11-12 | Fujifilm Corporation | Radiation-image processing device and method |
US11701079B2 (en) | 2013-11-25 | 2023-07-18 | Hologic, Inc. | Bone densitometer |
US11717244B2 (en) | 2015-02-26 | 2023-08-08 | Hologic, Inc. | Methods for physiological state determination in body scans |
US10412373B2 (en) | 2015-04-15 | 2019-09-10 | Google Llc | Image capture for virtual reality displays |
US11328446B2 (en) | 2015-04-15 | 2022-05-10 | Google Llc | Combining light-field data with active depth data for depth map generation |
US10341632B2 (en) | 2015-04-15 | 2019-07-02 | Google Llc. | Spatial random access enabled video system with a three-dimensional viewing volume |
US10275898B1 (en) | 2015-04-15 | 2019-04-30 | Google Llc | Wedge-based light-field video capture |
US10567464B2 (en) | 2015-04-15 | 2020-02-18 | Google Llc | Video compression with adaptive view-dependent lighting removal |
US10419737B2 (en) | 2015-04-15 | 2019-09-17 | Google Llc | Data structures and delivery methods for expediting virtual reality playback |
US10565734B2 (en) | 2015-04-15 | 2020-02-18 | Google Llc | Video capture, processing, calibration, computational fiber artifact removal, and light-field pipeline |
US10546424B2 (en) | 2015-04-15 | 2020-01-28 | Google Llc | Layered content delivery for virtual and augmented reality experiences |
US10469873B2 (en) | 2015-04-15 | 2019-11-05 | Google Llc | Encoding and decoding virtual reality video |
US10540818B2 (en) | 2015-04-15 | 2020-01-21 | Google Llc | Stereo image generation and interactive playback |
US9979909B2 (en) * | 2015-07-24 | 2018-05-22 | Lytro, Inc. | Automatic lens flare detection and correction for light-field images |
US10205896B2 (en) * | 2015-07-24 | 2019-02-12 | Google Llc | Automatic lens flare detection and correction for light-field images |
US20170026592A1 (en) * | 2015-07-24 | 2017-01-26 | Lytro, Inc. | Automatic lens flare detection and correction for light-field images |
US20170098315A1 (en) * | 2015-10-06 | 2017-04-06 | Rigaku Corporation | Analyzer, analysis method and analysis program of bone mineral density |
US9940737B2 (en) * | 2015-10-06 | 2018-04-10 | Rigaku Corporation | Analyzer, analysis method and analysis program of bone mineral density |
GB2565026A (en) * | 2016-05-03 | 2019-01-30 | Rapiscan Systems Inc | Radiation signal processing system |
CN109073572A (en) * | 2016-05-03 | 2018-12-21 | 拉皮斯坎系统股份有限公司 | radiation signal processing system |
GB2565026B (en) * | 2016-05-03 | 2021-08-18 | Rapiscan Systems Inc | Radiation signal processing system |
WO2017192597A1 (en) * | 2016-05-03 | 2017-11-09 | Rapiscan Systems, Inc. | Radiation signal processing system |
US10275892B2 (en) | 2016-06-09 | 2019-04-30 | Google Llc | Multi-view scene segmentation and propagation |
US10679361B2 (en) | 2016-12-05 | 2020-06-09 | Google Llc | Multi-view rotoscope contour propagation |
US10594945B2 (en) | 2017-04-03 | 2020-03-17 | Google Llc | Generating dolly zoom effect using light field image data |
US10474227B2 (en) | 2017-05-09 | 2019-11-12 | Google Llc | Generation of virtual reality with 6 degrees of freedom from limited viewer data |
US10444931B2 (en) | 2017-05-09 | 2019-10-15 | Google Llc | Vantage generation and interactive playback |
US10440407B2 (en) | 2017-05-09 | 2019-10-08 | Google Llc | Adaptive control for immersive experience delivery |
US10354399B2 (en) | 2017-05-25 | 2019-07-16 | Google Llc | Multi-view back-projection to a light-field |
US10545215B2 (en) | 2017-09-13 | 2020-01-28 | Google Llc | 4D camera tracking and optical stabilization |
US10965862B2 (en) | 2018-01-18 | 2021-03-30 | Google Llc | Multi-camera navigation interface |
Also Published As
Publication number | Publication date |
---|---|
US6315445B1 (en) | 2001-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6315445B1 (en) | Densitometry adapter for compact x-ray fluoroscopy machine | |
US6215846B1 (en) | Densitometry adapter for compact x-ray fluoroscopy machine | |
US6298109B1 (en) | X-ray imaging system | |
US6282261B1 (en) | Multi-mode x-ray image intensifier system | |
AU682480B2 (en) | Compact C-arm tomographic bone scanning system | |
US7697661B2 (en) | Method for obtaining a tomosynthesis image | |
US6990171B2 (en) | System and method of determining a user-defined region-of-interest of an imaging subject for x-ray flux management control | |
US5617462A (en) | Automatic X-ray exposure control system and method of use | |
EP1016375B1 (en) | Imaging system for generating high quality images | |
US7068750B2 (en) | System and method of x-ray flux management control | |
US7409034B2 (en) | Radiation tomography apparatus | |
US7289596B2 (en) | Radiation tomography apparatus and scan condition setting device | |
US7313217B2 (en) | System and method of collecting imaging subject positioning information for x-ray flux control | |
US8494118B2 (en) | Radiographic apparatus | |
JPH08509881A (en) | Method and apparatus for segmented radiography scanning | |
US7636422B2 (en) | X-ray CT apparatus and X-ray tube current determining method | |
AU8338891A (en) | Method for improving the dynamic range of an imaging system | |
JP4159188B2 (en) | Tube current adjusting method and apparatus, and X-ray CT apparatus | |
JP2003070780A (en) | Method to obtain data measured by computer tomograph and computer tomograph | |
EP1149559A1 (en) | Radiological imaging device, method and program of control associated with the device | |
US6870898B1 (en) | Computed tomography apparatus with automatic parameter modification to prevent impermissible operating states | |
EP0988830A2 (en) | Methods and apparatus for indirect high voltage verification in an X-ray imaging system | |
JPH1057361A (en) | Radiograph | |
JP4648355B2 (en) | Tube current adjusting method and apparatus, and X-ray CT apparatus | |
Xu et al. | Arbitrary shape region-of-interest fluoroscopy system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
CC | Certificate of correction | ||
FPAY | Fee payment |
Year of fee payment: 4 |
|
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: 8 |
|
REMI | Maintenance fee reminder mailed | ||
LAPS | Lapse for failure to pay maintenance fees | ||
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20131113 |