US20060289765A1 - Method and system for calibrating a computed tomography system - Google Patents
Method and system for calibrating a computed tomography system Download PDFInfo
- Publication number
- US20060289765A1 US20060289765A1 US11/159,706 US15970605A US2006289765A1 US 20060289765 A1 US20060289765 A1 US 20060289765A1 US 15970605 A US15970605 A US 15970605A US 2006289765 A1 US2006289765 A1 US 2006289765A1
- Authority
- US
- United States
- Prior art keywords
- elements
- crosstalk
- accordance
- detector
- determining
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/20—Measuring radiation intensity with scintillation detectors
- G01T1/2018—Scintillation-photodiode combinations
- G01T1/20183—Arrangements for preventing or correcting crosstalk, e.g. optical or electrical arrangements for correcting crosstalk
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
- G01T7/005—Details of radiation-measuring instruments calibration techniques
Definitions
- This invention relates generally to medical imaging systems, and more particularly, to methods and systems for calibrating a Computed Tomography (CT) system.
- CT Computed Tomography
- a CT system generally includes an image acquisition component, a processor and a controller.
- the image acquisition component has an X-ray source and a detector.
- Each detector includes a plurality of elements.
- Each element includes a scintillation crystal and a photosensor.
- Some examples of photosensors include photodiode pixels and phototransistors.
- the photosensors are broadly based on two kinds of technologies, namely, front-illuminated and back-illuminated. The two technologies are mainly distinguished by whether the electrical signal to the photosensors is routed from the front or the rear of the photosensors.
- the X-ray source directs X-ray beams through the object being scanned. These X-ray beams are attenuated by the object and are absorbed by the detector.
- the detector converts the X-ray radiation energy into electric signals. Ideally, the cumulative electric signal collected by the detector is linearly related to X-ray radiation energy that is absorbed by the detector. This criteria is fulfilled when requirements, such as, but not limited to, stability of the detector over time and temperature, sensitivity to focal spot motion and light output over life, are met.
- the electrical signals absorbed by the detector are reconstructed by the controller to form a CT image of the object.
- the calibrations in CT systems are used to remove artifacts from CT images.
- Artifacts such as smudges and center spots in CT images may be caused by crosstalk between the photosensors and/or the scintillation crystals (in a detector).
- the elements of the detector are connected and aligned in two directions, for example, the x and z directions.
- the close proximity of the photosensors causes crosstalk in both these directions.
- the crosstalk in both x and z directions include an average crosstalk and a module-to-module crosstalk.
- the module-to-module (channel-to-channel) crosstalk may be attributed to non-identical relative behavior in detector elements (module).
- the non-identical relative behavior includes factors, such as, fluctuation in the scintillator behavior from one scintillator element of the detector to an adjacent scintillator element and the lateral diffusion of photocarriers in the material of the photosensors.
- the lateral diffusion of photocarriers in the material of the photosensors causes photosensor electrical crosstalk. This results in photocarriers diffusing out of the photosensor collection junction area in which the photocarriers are generated. These photocarriers are then collected by adjacent or neighboring photosensors. This effect is more likely in back-illuminated photosensors than in front-illuminated photosensors, as the thickness of the back-illuminated photosensors increases the diffusion length before the photosensor collection junction.
- Module-to-module z-crosstalk is difficult to correct because of the photosensor alignment requirements of the back-illuminated photosensor.
- the scintillator cell is not aligned (or centered) to the diode active area, the crosstalk from that cell to its neighboring cell will be unbalanced and will cause image artifact when slope anatomy is imaged.
- the amount of crosstalk coming from the neighboring cell to the non-centered cell will be unbalanced when a slope anatomy is imaged.
- the z-crosstalk error is one of the factors that cause smudges and band artifacts in CT images obtained for objects that change in thickness sharply in the z-direction.
- Some known methods address the problem of crosstalk by aligning the scintillator array with the photosensors very accurately. This alignment requires a lot of accuracy in the mechanical dimensions of a scintillator and also needs to meet difficult-to-achieve tolerance requirements of the photosensor-scintillator alignment. Further, some known methods correct only one particular type of crosstalk and are unable to correct the others. For example, boosted reconstruction kernels correct the average crosstalk between the elements of the detector, but do not correct the other types of crosstalks such as differential crosstalk. In addition, known methods do not provide any robust solution for correcting image artifacts in the CT images of objects that have sharp slopes.
- a method and system for calibrating a computed tomography (CT) system includes selectively activating each of a plurality of elements of a detector in the CT system. The method further includes determining for each of the selectively activated elements, a crosstalk effect on elements adjacent to the activated elements, to calibrate the CT system.
- CT computed tomography
- a CT system in yet another exemplary embodiment of the invention, includes an image acquisition component, a controller, and a processor.
- the image acquisition component includes at least one detector having a plurality of elements.
- the controller is configured to selectively activate each of the plurality of elements of the detector.
- the processor is configured to determine for each of the selectively activated elements, a crosstalk effect on adjacent elements of the detector, to calibrate the image acquisition component.
- FIG. 1 is a block diagram of a CT system in accordance with an exemplary embodiment of the invention.
- FIG. 2 is a diagram illustrating a detector present in a CT system, in accordance with FIG. 1 .
- FIG. 3 is a flowchart of a method for calibrating a CT system, in accordance with an exemplary embodiment of the invention.
- FIG. 4 is a flowchart of a method for calibrating a CT system for z-crosstalk errors in accordance with an exemplary embodiment of the invention.
- FIG. 5 is a diagram illustrating a setup of a pre-test bay that is used for the measurement of z-crosstalk errors in accordance with FIG. 3 and FIG. 4 .
- FIG. 6 is a diagram illustrating calculating northward z-crosstalk error values in accordance with various embodiments of the invention.
- FIG. 7 is a diagram illustrating calculating southward z-crosstalk error values for photosensors in accordance with various embodiments of the invention.
- FIG. 8 is a diagram illustrating z-crosstalk error values in northward and southward directions for photosensors in accordance with various embodiments of the invention.
- FIGS. 9 and 10 are tables of collated sample z-crosstalk error values for the North and South directions, respectively, in accordance with FIG. 6 and FIG. 7 , respectively.
- FIG. 11 is a diagram illustrating a first algorithm for calibrating a CT system in accordance with various embodiments of the invention.
- FIG. 12 is a z-crosstalk matrix in accordance with various embodiments of the invention.
- FIG. 13 is a diagram illustrating a second algorithm for calibrating a CT system 300 in accordance with various embodiments of the invention.
- Various embodiments of the invention provide a method and system for calibrating a Computed Tomography (CT) system. Further, various embodiments of the invention provide for correcting crosstalk between the photosensors in the elements of a detector (in the CT system).
- the detector may be, for example, a VCT 64 detector, using multi-slice diodes arrays.
- the crosstalk correction is performed using average and differential crosstalk error values, which are determined prior to the generation of a CT image of an object.
- the absolute and the channel-to-channel crosstalks in the x, z and the diagonal directions of the CT system are corrected by the various embodiments described herein.
- FIG. 1 is a block diagram of a CT system 300 in accordance with an exemplary embodiment of the invention.
- CT system 300 includes an image acquisition component 302 , a controller 304 , a processor 306 , an image reconstruction device 308 , a data storage device 310 , an input device 312 , and a display device 314 .
- Image acquisition component 302 includes a gantry 316 , an X-ray source 318 , a detector 320 , a table 322 , and an object cavity 324 (e.g., a bore).
- Detector 320 includes a plurality of elements 326 . Each element 326 includes at least one scintillation crystal and photosensor.
- Controller 304 includes a table controller 328 , an X-ray controller 330 , a gantry controller 332 , and a data acquisition system (DAS) 334 .
- DAS data acquisition system
- X-ray source 318 and detector 320 are oppositely disposed in gantry 316 such that they are separated by object cavity 324 .
- An object 336 is placed on table 322 , and within object cavity 324 .
- X-ray source 318 and detector 320 are rotatingly disposed relative to gantry 316 and object 336 .
- X-ray source 318 and detector 320 remain stationary.
- X-ray source 318 projects an X-ray beam 338 toward detector 320 .
- X-ray beam 338 is configured to pass through object 336 .
- X-ray beam 338 is collimated to align X-ray beams 338 in an X-Y plane (of a Cartesian coordinate system), referred to as an “imaging plane”.
- X-ray beam 338 passes through object 336 , X-ray beam 338 is attenuated by object 336 .
- An attenuated X-ray beam 340 is absorbed by the scintillator crystals in elements 326 of detector 320 . These scintillation crystals convert the absorbed X-rays into visible light.
- the photosensors in elements 326 convert this light into electric signals. Each photosensor produces an electrical signal proportional to the intensity of attenuated X-ray beam 340 received by that photo sensor.
- Controller 304 is configured to selectively activate each element 326 in detector 320 with X-ray beam 338 .
- Table controller 328 controls the axial and rotational movement, and direction of table 322 .
- X-ray controller 330 provides power and timing signals to X-ray source 318 .
- Gantry controller 332 controls the rotational speed and angular position of X-ray source 318 .
- DAS 334 receives electrical signals produced by elements 326 and digitizes the signals for subsequent processing. The digitized signals are provided to image reconstruction device 308 .
- Image reconstruction device 308 and processor 306 reconstruct a CT image from these signals. The reconstructed CT image is stored in data storage device 310 .
- data storage device 310 may also save any intermediately processed data during the formation of the CT image.
- image reconstruction device 308 determines the effect of crosstalk on the photosensors of elements 326 adjacent to the photosensors of the activated elements 326 in detector 320 .
- Input device 312 is configured to receive manual inputs from a user of CT system 300 .
- the CT image of object 336 is displayed by display device 314 .
- display device 314 is used to determine the calibration values of CT system 300 .
- data storage device 310 may be a magnetic or an optical storage media, such as, but not limited to, a hard disk, a tape drive, a compact disc (CD), and a memory chip.
- computer programs and other instructions may be uploaded into processor 306 through input device 312 .
- input device 312 include, but are not limited to, keys/buttons, audio inputs, and video input devices.
- display device 314 includes one of, but is not limited to, a cathode ray display, a Liquid Crystal Display (LCD), and plasma display.
- LCD Liquid Crystal Display
- FIG. 2 illustrates detector 320 in CT system 300 , in accordance with FIG. 1 .
- detector 320 is configured as a multi-dimensional array of elements 326 .
- Each element in the multi-dimensional array of elements 326 is hereinafter referred to as an element 326 .
- the array is configured as an arch length segment along the x direction and a linear segment along the z direction.
- a collection of elements 326 along the arch length segment of detector 320 is hereinafter referred to as a channel.
- a collection of elements 326 along the linear depth of detector 320 is hereinafter referred to as a row.
- the configuration of detector 320 includes N arrays of elements 326 . Each array has 64 rows and 16 channels. 64 rows allow the imaging of object 336 in 64 slices, in the z direction. The thickness of each row is the slice thickness in the z direction.
- detector 320 may be configured as a combination of multi-dimensional arrays, in accordance with an exemplary embodiment of the invention.
- An example of such a configuration is the combination of a scintillator array and a photosensor array.
- the scintillator array includes a plurality of scintillator crystals arranged in an array.
- the photosensor array includes a plurality of photosensors arranged in an array identical in dimensions to the scintillator array.
- detector 320 includes a scintillating medium across the length of detector 320 .
- Detector 320 further includes a photosensor array that is aligned to the scintillating medium.
- FIG. 3 is a flowchart of a method for calibrating CT system 300 , in accordance with an exemplary embodiment of the invention. Specifically, the flowchart describes a method for measuring the z-crosstalk between various photosensors in elements 326 of detector 320 .
- each element of the detector is selectively activated by X-rays.
- the effect on the photosensors of the elements adjacent to each activated element is measured. For example, the X-ray radiation directed towards a particular element may be spilt onto the photosensors of adjacent elements and is measured accordingly.
- the fraction of signal flowing from the activated pixel to its neighboring pixels is a measure of the crosstalk. If this fraction is known prior to imaging, one can correct the images for this artifact.
- FIG. 4 is a flowchart of a method for calibrating a CT system 300 for z-crosstalk errors, in accordance with an exemplary embodiment of the invention. Specifically, the flowchart describes a method for measuring the z-crosstalk error values for each photosensor in elements 326 of detector 320 , and using these z-crosstalk error values for correcting for z-crosstalk.
- each element 326 of detector 320 is selectively activated by X-rays.
- the northward z-crosstalk error values for each photosensor in detector 320 are calculated.
- the southward z-crosstalks for each photosensor of detector 320 are calculated.
- average northward and southward z-crosstalk error values for each photosensor in detector 320 are calculated.
- the average northward and southward z-crosstalk error values are used for the correction of the z-crosstalk between the photosensors in elements 326 of detector 320 .
- the measurement of the z-crosstalk error values is performed in a pre-test bay before the CT image of object 336 is taken. These z-crosstalk error values are used for the correction of the z-crosstalk and for the calibration of CT system 300 .
- FIG. 5 illustrates a setup of a pre-test bay for the measurement of z-crosstalk error values for the photosensors in a channel in detector 320 in accordance with FIG. 3 and FIG. 4 .
- the setup includes X-ray source 318 , X-ray beam 338 , a slit 702 , and elements 326 .
- slit 702 may be replaced by a pinhole.
- Each element 326 includes a scintillator crystal 704 and a photosensor 706 .
- Detector 320 includes a plurality of photosensors similar to photosensor 706 , which are hereinafter referred to as photosensors 706 .
- the same process is applicable for crosstalk correction in the channel direction (x direction) as well as in the diagonal direction (between x and z direction).
- a channel ‘i’ with R+1 rows of elements 326 are provided.
- Element 326 in channel ‘i’ and row ‘R’ is hereinafter referred to as P(i,R).
- element 326 in channel ‘i’ and row ‘R+1’ is hereinafter referred to as P(i,R+1)
- element 326 in channel ‘i’ and row ‘R ⁇ 1’ is hereinafter referred to as P(i,R ⁇ 1).
- Channel ‘i’ includes, but is not limited to, element P(i,R) between elements P(i,R+1) and P(i,R ⁇ 1).
- the direction from P(i,R) to P(i,R ⁇ 1) is North and the direction from P(i,R) to P(i,R+1) is South. Further, the North and the South directions lie along the z coordinate axis. In an exemplary embodiment of the invention, each row of elements 326 is along the z direction, and each channel of elements 326 is along the x direction.
- Slit 702 is moved along each element 326 in channel ‘i’, in order to measure the crosstalk for each element 326 in that channel ‘i’.
- X-ray beam 338 is focused on, for example, P(i,R) through slit 702 .
- XX-ray beam 338 is absorbed partially by photosensor 706 in P(i,R) and a portion of the radiation is transferred to photosensors 706 in adjacent or neighboring elements P(i,R ⁇ 1) and P(i,R+1), in the North and South directions, respectively and in X-direction as well (P(i ⁇ 1,R) and P(i+1,R)).
- slit 702 has a width equivalent to the pixel size and each element is in the shape of a cube having sides of 0.625 mm according to International Organization for Standardization (ISO) standards.
- ISO International Organization for Standardization
- FIG. 6 illustrates a method for calculating the northward z-crosstalk error values for photosensors 706 present in channel ‘i’ of elements 326 in detector 320 , in accordance with various embodiments of the invention.
- 80% of the signal intensity at P(i,R ⁇ 1) is from adjacent element P(i,R)
- 10% of the signal intensity is from diagonal elements P(i ⁇ 1,R), and P(i+ 1 ,R). Therefore, while calculating the northward z-crosstalk error for P(i,R), only 80% the signal intensity at P(i,R ⁇ 1) is considered.
- the southward z-crosstalk error for P(i,R) is calculated by the following equation: e N ⁇ [ ( 100 * Signal ⁇ ( i , R - 1 ) Signal ⁇ ( i , R ) ) * 0.8 ] ( 1 )
- ‘i’ varies based on the number of elements 326 , for example, from 1 to M, wherein M is the number of channels in the detector and R varies based on the number of element slices, for example, from 1 to 64.
- the northward z-crosstalk on other elements 326 is calculated using equation 1 by varying ‘i’ and ‘R’.
- FIG. 7 illustrates a method for calculating the southward z-crosstalk error values for photosensors 706 in channel ‘i’ of elements 326 in detector 320 , in accordance with various embodiments of the invention.
- X-ray radiations directed on and detected by photosensors 706 of P(i,R), P(i ⁇ 1,R), and P(i+1,R) also may be detected by P(i,R+1). For example, 80% of the contribution may be from adjacent element P(i,R) and 10% contribution each from diagonal elements P(i ⁇ 1,R), and P(i+1,R). Therefore, while calculating the southward z-crosstalk error for P(i,R), only 80% the signal intensity at P(i,R+1) is considered.
- ‘i’ varies from 1 to M, wherein M is the number of detector channels and R varies from 1 to 64.
- the southward z-crosstalk on other elements 326 is calculated using equation 2 by varying ‘i’ and ‘R’.
- FIG. 8 illustrates the z-crosstalk error values in northward and southward directions for photosensors 706 present in channel ‘i’ of elements 326 of detector 320 , in accordance with various embodiments of the invention.
- FIGS. 9 and 10 are tables illustrating collated sample z-crosstalk error values for the northward and southward directions, respectively, in accordance with FIG. 6 and FIG. 7 .
- Each row of elements 326 is referred to as a slice.
- the northward z-crosstalk error value in each element 326 in slice 1 will be zero.
- the southward z-crosstalk error value in each element 326 in slice 64 will be zero.
- the northward z-crosstalk error values over all photosensors 706 are used “as is” or averaged to obtain average northward z-crosstalk error e N and the southward z-crosstalk error values over all photosensors 706 are used “as is” or averaged to obtain average southward z-crosstalk error e S .
- These average z-crosstalk error values are used to correct the z-crosstalks for all photosensors 706 in detector 320 during the generation of the CT image of object 336 .
- the correction of the northward and southward z-crosstalks may be provided according to various embodiments and as described below.
- FIG. 11 illustrates a first algorithm for correcting CT system 300 , in accordance with various embodiments of the invention.
- the algorithm corrects of the z-crosstalk for photosensors 706 in a channel having 64 elements 326 , which are hereinafter referred to as elements 1 - 64 . Further, each element 326 is hereinafter referred to as element R, wherein R ranges from 1 through 64.
- the algorithm may be used to calculate the correction for more than or less than 64 elements 326 in a channel of elements 326 .
- the collection vector intensities D 1 -D 64 do not include crosstalk at elements 1 - 64 . Due to the z-crosstalk in the northward and southward directions between elements 1 - 64 , elements 1 - 64 emit measured intensities S 1 to S 64 (S 1 -S 64 ), respectively.
- the average z-crosstalk error values e N and e S are used for determining collection vector intensities D 1 -D 64 using measured intensities S 1 -S 64 .
- This matrix equation is de-convolved to obtain the collection vector intensities D 1 -D 64 .
- This matrix equation is de-convolved in order to obtain the collection vector intensities values that would have been obtained if there was no crosstalk between photosensors 706 in elements 326 of detector 320 .
- FIG. 12 illustrates a z-crosstalk matrix in accordance with various embodiments of the invention.
- the matrix needs to be solved to determine the collection vector intensities D 1 -D 64 before the crosstalk is determined.
- the collection vector intensities D 1 -D 64 are the crosstalk corrected X-ray signals to be used in the reconstruction of the CT images.
- a second algorithm provides an approximate correction of the z-crosstalk error values between photosensors 706 in elements 326 of detector 320 .
- FIG. 13 illustrates a second algorithm for calibrating CT system 300 , in accordance with various embodiments of the invention. Specifically, the second algorithm describes a method for the correction of z-crosstalk in order to calibrate CT system 300 .
- the correction of elements 1 - 64 is calculated using this algorithm. However, this algorithm may be used to calculate the correction for more than or less than 64 elements in a channel of elements 326 .
- X-ray beams 338 with initial intensities I 1 ⁇ I 64 impinge on and are attenuated by object 336 (shown in FIG. 1 ). After passing through object 336 , attenuated X-ray beams are received by elements 1 - 64 , and have collection vector intensities D 1 -D 64 , respectively.
- the collection vector intensities D 1 -D 64 are the intensities that impinge on elements 1 - 64 , respectively, before the z-crosstalk correction is performed. Due to z-crosstalk in the northward and southward directions between elements 1 - 64 , elements 1 - 64 finally absorb measured intensities S 1 -S 64 , respectively.
- the average z-crosstalk error values e N and e S are used for determining collection vector intensities D 1 -D 64 using measured intensities S 1 -S 64 . Further, for elements R and R+1, ⁇ R(R+1) represents the ratio of difference between collection vector intensities D R and D R+1 , and D R . In various embodiments of the invention, the average z-crosstalk error values e N and e S , and ⁇ R(R+1) are used for correcting the z-crosstalk between elements 1 - 64 .
- D R for element R is a function of S R , ⁇ R(R ⁇ 1) , ⁇ R(R+1) , e N and e S is represented by the following equations:
- D R S R / ( 1 + ⁇ R ⁇ ( R - 1 ) ⁇ e N - ⁇ R ⁇ ( R + 1 ) ⁇ e S ) ) ( 15 )
- ⁇ R ⁇ (R+1) For a flat anatomy where there is no slope, ⁇ R ⁇ (R+1) will be equal to zero and no crosstalk correction will be performed. However, if there is a slope in the anatomy, leading to different X-ray signals in two contiguous rows, then correction will be performed. Because ⁇ R ⁇ (R+1) is based on primary beam signals without crosstalk effects, D R , an approximate value of ⁇ R ⁇ (R+1) can be calculated using the measured intensities S R , which are the X-ray signals that are detected.
- ⁇ R(R+1) ( S R ⁇ S R+1 )/ S R (17)
- the first algorithm and the second algorithm are used for correcting artifacts induced due to z-crosstalk error values between photosensors 706 in detector 320 . Further, these algorithms may be used to correct the x-crosstalk between photosensors 706 in detector 320 .
- Various embodiments of the invention allow for producing accurate CT images without any artifacts or smudges specially, the CT images objects with steep slopes such as liver and the tip of head.
- the z-crosstalk error values e N and e S are used to correct the collection vector intensities D 1 -D 64 . These error values are used to calibrate the CT system while performing image reconstruction.
- One technical effect of the present invention is that accurate manual alignment of the scintillator crystals with the photosensors is not needed in order to obtain accurate images. Also, z-crosstalk error values e N and e S are used to correct the collection vector intensities D 1 -D 64 . These error values are used to calibrate the CT system while performing image reconstruction.
- the various embodiments of the invention and components thereof may be implemented as part of a computer system.
- the computer system may include a computer, an input device, a display unit and an interface, for example, for accessing the Internet.
- the term “computer” may include any processor-based or microprocessor-based system including systems using microcontrollers, reduced instruction set circuits (RISC), application specific integrated circuits (ASICs), logic circuits, and any other circuit or processor capable of executing the functions described herein.
- RISC reduced instruction set circuits
- ASICs application specific integrated circuits
- the above examples are exemplary only, and are thus not intended to limit in any way the definition and/or meaning of the term “computer”.
- the computer system executes a set of instructions that are stored in one or more storage elements, in order to process input data.
- the storage elements may also hold data or other information as desired or needed.
- the storage element may be in the form of an information source or a physical memory element within the processing machine.
- the set of instructions may include various commands that instruct the processing machine to perform specific operations such as the processes of the various embodiments of the invention.
- the set of instructions may be in the form of a software program.
- the software may be in various forms such as system software or application software. Further, the software may be in the form of a collection of separate programs, a program module within a larger program or a portion of a program module.
- the software also may include modular programming in the form of object-oriented programming.
- the processing of input data by the processing machine may be in response to user commands, or in response to results of previous processing, or in response to a request made by another processing machine.
Abstract
Description
- This invention relates generally to medical imaging systems, and more particularly, to methods and systems for calibrating a Computed Tomography (CT) system.
- CT systems may be used for imaging tissues and bones in various body parts. A CT system generally includes an image acquisition component, a processor and a controller. The image acquisition component has an X-ray source and a detector. Each detector includes a plurality of elements. Each element includes a scintillation crystal and a photosensor. Some examples of photosensors include photodiode pixels and phototransistors. The photosensors are broadly based on two kinds of technologies, namely, front-illuminated and back-illuminated. The two technologies are mainly distinguished by whether the electrical signal to the photosensors is routed from the front or the rear of the photosensors.
- The X-ray source directs X-ray beams through the object being scanned. These X-ray beams are attenuated by the object and are absorbed by the detector. The detector converts the X-ray radiation energy into electric signals. Ideally, the cumulative electric signal collected by the detector is linearly related to X-ray radiation energy that is absorbed by the detector. This criteria is fulfilled when requirements, such as, but not limited to, stability of the detector over time and temperature, sensitivity to focal spot motion and light output over life, are met. The electrical signals absorbed by the detector are reconstructed by the controller to form a CT image of the object.
- With advancements in technology in CT systems, there has been a growing need for more precise calibrations and alignments. The calibrations in CT systems are used to remove artifacts from CT images. Artifacts such as smudges and center spots in CT images may be caused by crosstalk between the photosensors and/or the scintillation crystals (in a detector).
- In a detector using back-illuminated photosensors, the elements of the detector are connected and aligned in two directions, for example, the x and z directions. The close proximity of the photosensors causes crosstalk in both these directions. With the higher number of slices in multi-slices detectors and thin slices, the crosstalk becomes an important parameter to optimize and control. The crosstalk in both x and z directions include an average crosstalk and a module-to-module crosstalk. The module-to-module (channel-to-channel) crosstalk may be attributed to non-identical relative behavior in detector elements (module). The non-identical relative behavior includes factors, such as, fluctuation in the scintillator behavior from one scintillator element of the detector to an adjacent scintillator element and the lateral diffusion of photocarriers in the material of the photosensors. The lateral diffusion of photocarriers in the material of the photosensors causes photosensor electrical crosstalk. This results in photocarriers diffusing out of the photosensor collection junction area in which the photocarriers are generated. These photocarriers are then collected by adjacent or neighboring photosensors. This effect is more likely in back-illuminated photosensors than in front-illuminated photosensors, as the thickness of the back-illuminated photosensors increases the diffusion length before the photosensor collection junction.
- Module-to-module z-crosstalk is difficult to correct because of the photosensor alignment requirements of the back-illuminated photosensor. In case the scintillator cell is not aligned (or centered) to the diode active area, the crosstalk from that cell to its neighboring cell will be unbalanced and will cause image artifact when slope anatomy is imaged. At the same time, the amount of crosstalk coming from the neighboring cell to the non-centered cell will be unbalanced when a slope anatomy is imaged.
- Further, the z-crosstalk error is one of the factors that cause smudges and band artifacts in CT images obtained for objects that change in thickness sharply in the z-direction.
- Some known methods address the problem of crosstalk by aligning the scintillator array with the photosensors very accurately. This alignment requires a lot of accuracy in the mechanical dimensions of a scintillator and also needs to meet difficult-to-achieve tolerance requirements of the photosensor-scintillator alignment. Further, some known methods correct only one particular type of crosstalk and are unable to correct the others. For example, boosted reconstruction kernels correct the average crosstalk between the elements of the detector, but do not correct the other types of crosstalks such as differential crosstalk. In addition, known methods do not provide any robust solution for correcting image artifacts in the CT images of objects that have sharp slopes.
- In an exemplary embodiment of the invention, a method and system for calibrating a computed tomography (CT) system is provided. The method includes selectively activating each of a plurality of elements of a detector in the CT system. The method further includes determining for each of the selectively activated elements, a crosstalk effect on elements adjacent to the activated elements, to calibrate the CT system.
- In yet another exemplary embodiment of the invention, a CT system is provided. The CT system includes an image acquisition component, a controller, and a processor. The image acquisition component includes at least one detector having a plurality of elements. The controller is configured to selectively activate each of the plurality of elements of the detector. The processor is configured to determine for each of the selectively activated elements, a crosstalk effect on adjacent elements of the detector, to calibrate the image acquisition component.
-
FIG. 1 is a block diagram of a CT system in accordance with an exemplary embodiment of the invention. -
FIG. 2 is a diagram illustrating a detector present in a CT system, in accordance withFIG. 1 . -
FIG. 3 is a flowchart of a method for calibrating a CT system, in accordance with an exemplary embodiment of the invention. -
FIG. 4 is a flowchart of a method for calibrating a CT system for z-crosstalk errors in accordance with an exemplary embodiment of the invention. -
FIG. 5 is a diagram illustrating a setup of a pre-test bay that is used for the measurement of z-crosstalk errors in accordance withFIG. 3 andFIG. 4 . -
FIG. 6 is a diagram illustrating calculating northward z-crosstalk error values in accordance with various embodiments of the invention. -
FIG. 7 is a diagram illustrating calculating southward z-crosstalk error values for photosensors in accordance with various embodiments of the invention. -
FIG. 8 is a diagram illustrating z-crosstalk error values in northward and southward directions for photosensors in accordance with various embodiments of the invention. -
FIGS. 9 and 10 are tables of collated sample z-crosstalk error values for the North and South directions, respectively, in accordance withFIG. 6 andFIG. 7 , respectively. -
FIG. 11 is a diagram illustrating a first algorithm for calibrating a CT system in accordance with various embodiments of the invention. -
FIG. 12 is a z-crosstalk matrix in accordance with various embodiments of the invention. -
FIG. 13 is a diagram illustrating a second algorithm for calibrating aCT system 300 in accordance with various embodiments of the invention. - Various embodiments of the invention provide a method and system for calibrating a Computed Tomography (CT) system. Further, various embodiments of the invention provide for correcting crosstalk between the photosensors in the elements of a detector (in the CT system). In various exemplary embodiments of the invention, the detector may be, for example, a
VCT 64 detector, using multi-slice diodes arrays. The crosstalk correction is performed using average and differential crosstalk error values, which are determined prior to the generation of a CT image of an object. In various embodiments of the invention, the absolute and the channel-to-channel crosstalks in the x, z and the diagonal directions of the CT system are corrected by the various embodiments described herein. -
FIG. 1 is a block diagram of aCT system 300 in accordance with an exemplary embodiment of the invention.CT system 300 includes animage acquisition component 302, acontroller 304, aprocessor 306, animage reconstruction device 308, adata storage device 310, aninput device 312, and adisplay device 314.Image acquisition component 302 includes agantry 316, anX-ray source 318, adetector 320, a table 322, and an object cavity 324 (e.g., a bore).Detector 320 includes a plurality ofelements 326. Eachelement 326 includes at least one scintillation crystal and photosensor. Some examples of photosensors include, but are not limited to, photodiode pixels and phototransistors.Controller 304 includes atable controller 328, anX-ray controller 330, agantry controller 332, and a data acquisition system (DAS) 334. - In various embodiments,
X-ray source 318 anddetector 320 are oppositely disposed ingantry 316 such that they are separated byobject cavity 324. Anobject 336 is placed on table 322, and withinobject cavity 324. In an exemplary embodiment of the invention,X-ray source 318 anddetector 320 are rotatingly disposed relative to gantry 316 andobject 336. In another exemplary embodiment of the invention,X-ray source 318 anddetector 320 remain stationary. X-raysource 318 projects anX-ray beam 338 towarddetector 320.X-ray beam 338 is configured to pass throughobject 336.X-ray beam 338 is collimated to alignX-ray beams 338 in an X-Y plane (of a Cartesian coordinate system), referred to as an “imaging plane”. WhenX-ray beam 338 passes throughobject 336,X-ray beam 338 is attenuated byobject 336. Anattenuated X-ray beam 340 is absorbed by the scintillator crystals inelements 326 ofdetector 320. These scintillation crystals convert the absorbed X-rays into visible light. The photosensors inelements 326 convert this light into electric signals. Each photosensor produces an electrical signal proportional to the intensity ofattenuated X-ray beam 340 received by that photo sensor. -
Controller 304 is configured to selectively activate eachelement 326 indetector 320 withX-ray beam 338.Table controller 328 controls the axial and rotational movement, and direction of table 322.X-ray controller 330 provides power and timing signals to X-raysource 318.Gantry controller 332 controls the rotational speed and angular position ofX-ray source 318.DAS 334 receives electrical signals produced byelements 326 and digitizes the signals for subsequent processing. The digitized signals are provided to imagereconstruction device 308.Image reconstruction device 308 andprocessor 306 reconstruct a CT image from these signals. The reconstructed CT image is stored indata storage device 310. In one embodiment of the invention,data storage device 310 may also save any intermediately processed data during the formation of the CT image. In addition,image reconstruction device 308 determines the effect of crosstalk on the photosensors ofelements 326 adjacent to the photosensors of the activatedelements 326 indetector 320.Input device 312 is configured to receive manual inputs from a user ofCT system 300. The CT image ofobject 336 is displayed bydisplay device 314. In an exemplary embodiment of the invention,display device 314 is used to determine the calibration values ofCT system 300. - In various embodiments of the invention,
data storage device 310 may be a magnetic or an optical storage media, such as, but not limited to, a hard disk, a tape drive, a compact disc (CD), and a memory chip. In an exemplary embodiment of the invention, computer programs and other instructions may be uploaded intoprocessor 306 throughinput device 312. Examples ofinput device 312 include, but are not limited to, keys/buttons, audio inputs, and video input devices. In various embodiments of the invention,display device 314 includes one of, but is not limited to, a cathode ray display, a Liquid Crystal Display (LCD), and plasma display. -
FIG. 2 illustratesdetector 320 inCT system 300, in accordance withFIG. 1 . In an exemplary embodiment of the invention,detector 320 is configured as a multi-dimensional array ofelements 326. Each element in the multi-dimensional array ofelements 326 is hereinafter referred to as anelement 326. In an embodiment of the invention, the array is configured as an arch length segment along the x direction and a linear segment along the z direction. - A collection of
elements 326 along the arch length segment ofdetector 320 is hereinafter referred to as a channel. A collection ofelements 326 along the linear depth ofdetector 320 is hereinafter referred to as a row. In an exemplary embodiment of the invention, the configuration ofdetector 320 includes N arrays ofelements 326. Each array has 64 rows and 16 channels. 64 rows allow the imaging ofobject 336 in 64 slices, in the z direction. The thickness of each row is the slice thickness in the z direction. - In another embodiment of the invention,
detector 320 may be configured as a combination of multi-dimensional arrays, in accordance with an exemplary embodiment of the invention. An example of such a configuration is the combination of a scintillator array and a photosensor array. The scintillator array includes a plurality of scintillator crystals arranged in an array. The photosensor array includes a plurality of photosensors arranged in an array identical in dimensions to the scintillator array. In alternate embodiments of the invention,detector 320 includes a scintillating medium across the length ofdetector 320.Detector 320 further includes a photosensor array that is aligned to the scintillating medium. -
FIG. 3 is a flowchart of a method for calibratingCT system 300, in accordance with an exemplary embodiment of the invention. Specifically, the flowchart describes a method for measuring the z-crosstalk between various photosensors inelements 326 ofdetector 320. At 502, each element of the detector is selectively activated by X-rays. At 504, the effect on the photosensors of the elements adjacent to each activated element is measured. For example, the X-ray radiation directed towards a particular element may be spilt onto the photosensors of adjacent elements and is measured accordingly. The fraction of signal flowing from the activated pixel to its neighboring pixels is a measure of the crosstalk. If this fraction is known prior to imaging, one can correct the images for this artifact. -
FIG. 4 is a flowchart of a method for calibrating aCT system 300 for z-crosstalk errors, in accordance with an exemplary embodiment of the invention. Specifically, the flowchart describes a method for measuring the z-crosstalk error values for each photosensor inelements 326 ofdetector 320, and using these z-crosstalk error values for correcting for z-crosstalk. At 601, eachelement 326 ofdetector 320 is selectively activated by X-rays. At 602, the northward z-crosstalk error values for each photosensor indetector 320 are calculated. At 604, the southward z-crosstalks for each photosensor ofdetector 320 are calculated. At 606, average northward and southward z-crosstalk error values for each photosensor indetector 320 are calculated. At 608, the average northward and southward z-crosstalk error values are used for the correction of the z-crosstalk between the photosensors inelements 326 ofdetector 320. - In various embodiments of the invention, the measurement of the z-crosstalk error values is performed in a pre-test bay before the CT image of
object 336 is taken. These z-crosstalk error values are used for the correction of the z-crosstalk and for the calibration ofCT system 300. -
FIG. 5 illustrates a setup of a pre-test bay for the measurement of z-crosstalk error values for the photosensors in a channel indetector 320 in accordance withFIG. 3 andFIG. 4 . As shown, the setup includesX-ray source 318,X-ray beam 338, aslit 702, andelements 326. In an embodiment of the invention, slit 702 may be replaced by a pinhole. Eachelement 326 includes ascintillator crystal 704 and aphotosensor 706.Detector 320 includes a plurality of photosensors similar tophotosensor 706, which are hereinafter referred to asphotosensors 706. The same process is applicable for crosstalk correction in the channel direction (x direction) as well as in the diagonal direction (between x and z direction). - A channel ‘i’ with R+1 rows of
elements 326 are provided.Element 326 in channel ‘i’ and row ‘R’ is hereinafter referred to as P(i,R). Similarly,element 326 in channel ‘i’ and row ‘R+1’ is hereinafter referred to as P(i,R+1), andelement 326 in channel ‘i’ and row ‘R−1’ is hereinafter referred to as P(i,R−1). Channel ‘i’ includes, but is not limited to, element P(i,R) between elements P(i,R+1) and P(i,R−1). - In an embodiment of the invention, the direction from P(i,R) to P(i,R−1) is North and the direction from P(i,R) to P(i,R+1) is South. Further, the North and the South directions lie along the z coordinate axis. In an exemplary embodiment of the invention, each row of
elements 326 is along the z direction, and each channel ofelements 326 is along the x direction. -
Slit 702 is moved along eachelement 326 in channel ‘i’, in order to measure the crosstalk for eachelement 326 in that channel ‘i’.X-ray beam 338 is focused on, for example, P(i,R) throughslit 702. However, XX-ray beam 338 is absorbed partially byphotosensor 706 in P(i,R) and a portion of the radiation is transferred tophotosensors 706 in adjacent or neighboring elements P(i,R−1) and P(i,R+1), in the North and South directions, respectively and in X-direction as well (P(i−1,R) and P(i+1,R)). - In an embodiment of the invention, slit 702 has a width equivalent to the pixel size and each element is in the shape of a cube having sides of 0.625 mm according to International Organization for Standardization (ISO) standards.
-
FIG. 6 illustrates a method for calculating the northward z-crosstalk error values forphotosensors 706 present in channel ‘i’ ofelements 326 indetector 320, in accordance with various embodiments of the invention. - X-ray radiations directed on and detected by P(i,R), P(i−1,R), and P(i+1,R), also may be detected by P(i,R−1). In various embodiments of the invention, 80% of the signal intensity at P(i,R−1) is from adjacent element P(i,R), and 10% of the signal intensity is from diagonal elements P(i−1,R), and P(i+1,R). Therefore, while calculating the northward z-crosstalk error for P(i,R), only 80% the signal intensity at P(i,R−1) is considered. Thus, the southward z-crosstalk error for P(i,R) is calculated by the following equation:
In an embodiment of the invention, ‘i’ varies based on the number ofelements 326, for example, from 1 to M, wherein M is the number of channels in the detector and R varies based on the number of element slices, for example, from 1 to 64. The northward z-crosstalk onother elements 326 is calculated usingequation 1 by varying ‘i’ and ‘R’. -
FIG. 7 illustrates a method for calculating the southward z-crosstalk error values forphotosensors 706 in channel ‘i’ ofelements 326 indetector 320, in accordance with various embodiments of the invention. - X-ray radiations directed on and detected by
photosensors 706 of P(i,R), P(i−1,R), and P(i+1,R) also may be detected by P(i,R+1). For example, 80% of the contribution may be from adjacent element P(i,R) and 10% contribution each from diagonal elements P(i−1,R), and P(i+1,R). Therefore, while calculating the southward z-crosstalk error for P(i,R), only 80% the signal intensity at P(i,R+1) is considered. Thus, the southward z-crosstalk error for P(i,R) is calculated by the following equation:
In an embodiment of the invention, ‘i’ varies from 1 to M, wherein M is the number of detector channels and R varies from 1 to 64. The southward z-crosstalk onother elements 326 is calculated usingequation 2 by varying ‘i’ and ‘R’. -
FIG. 8 illustrates the z-crosstalk error values in northward and southward directions forphotosensors 706 present in channel ‘i’ ofelements 326 ofdetector 320, in accordance with various embodiments of the invention. -
FIGS. 9 and 10 are tables illustrating collated sample z-crosstalk error values for the northward and southward directions, respectively, in accordance withFIG. 6 andFIG. 7 . Each row ofelements 326 is referred to as a slice. According toequation 1, the northward z-crosstalk error value in eachelement 326 inslice 1 will be zero. According toequation 2, the southward z-crosstalk error value in eachelement 326 inslice 64 will be zero. - The northward z-crosstalk error values over all
photosensors 706 are used “as is” or averaged to obtain average northward z-crosstalk error eN and the southward z-crosstalk error values over allphotosensors 706 are used “as is” or averaged to obtain average southward z-crosstalk error eS. These average z-crosstalk error values are used to correct the z-crosstalks for allphotosensors 706 indetector 320 during the generation of the CT image ofobject 336. - The correction of the northward and southward z-crosstalks may be provided according to various embodiments and as described below.
-
FIG. 11 illustrates a first algorithm for correctingCT system 300, in accordance with various embodiments of the invention. The algorithm corrects of the z-crosstalk forphotosensors 706 in a channel having 64elements 326, which are hereinafter referred to as elements 1-64. Further, eachelement 326 is hereinafter referred to as element R, wherein R ranges from 1 through 64. In alternate embodiments of the invention, the algorithm may be used to calculate the correction for more than or less than 64elements 326 in a channel ofelements 326. - As shown in
FIG. 11 , X-ray beams 338 with initial intensities I1 to I64 (I1-I64) impinge on and are attenuated by object 336 (shown inFIG. 1 ). After passing throughobject 336, attenuated X-ray beams are received by elements 1-64, respectively, and have collection vector intensities D1 to D64 (D1-D64), respectively. - The collection vector intensities D1-D64 do not include crosstalk at elements 1-64. Due to the z-crosstalk in the northward and southward directions between elements 1-64, elements 1-64 emit measured intensities S1 to S64 (S1-S64), respectively. The average z-crosstalk error values eN and eS are used for determining collection vector intensities D1-D64 using measured intensities S1-S64.
- S1 for
element 1 is a function of D1 and D2 and is represented by the following equation:
S 1=(1−e S)D 1 +e S D 2 (3)
S7 forelement 7 is a function of D6, D7 and D8 and is represented by the following equation:
S 7=(1−e N −e S)D 7 +e N D 6 +e SD8 (4)
Finally, S64 forelement 64 is a function of D63 and D64 and is represented by the following equation:
S 64=(1−e N)D 64 +e N D 63 (5)
Usingequations
S R=(1−e N e S)D R +e N D R−1 +e S D R+1 (6)
Finally, equations relating each S value with each corresponding D value are written in form of the z-crosstalk matrix equation [S]=[ε]*[D].
This matrix equation is de-convolved to obtain the collection vector intensities D1-D64. This matrix equation is de-convolved in order to obtain the collection vector intensities values that would have been obtained if there was no crosstalk betweenphotosensors 706 inelements 326 ofdetector 320. -
FIG. 12 illustrates a z-crosstalk matrix in accordance with various embodiments of the invention. The matrix needs to be solved to determine the collection vector intensities D1-D64 before the crosstalk is determined. The collection vector intensities D1-D64 are the crosstalk corrected X-ray signals to be used in the reconstruction of the CT images. - A second algorithm provides an approximate correction of the z-crosstalk error values between
photosensors 706 inelements 326 ofdetector 320. -
FIG. 13 illustrates a second algorithm for calibratingCT system 300, in accordance with various embodiments of the invention. Specifically, the second algorithm describes a method for the correction of z-crosstalk in order to calibrateCT system 300. - The correction of elements 1-64 is calculated using this algorithm. However, this algorithm may be used to calculate the correction for more than or less than 64 elements in a channel of
elements 326. - As shown, X-ray beams 338 with initial intensities I1−I64 impinge on and are attenuated by object 336 (shown in
FIG. 1 ). After passing throughobject 336, attenuated X-ray beams are received by elements 1-64, and have collection vector intensities D1-D64, respectively. - The collection vector intensities D1-D64 are the intensities that impinge on elements 1-64, respectively, before the z-crosstalk correction is performed. Due to z-crosstalk in the northward and southward directions between elements 1-64, elements 1-64 finally absorb measured intensities S1-S64, respectively. The average z-crosstalk error values eN and eS are used for determining collection vector intensities D1-D64 using measured intensities S1-S64. Further, for elements R and R+1, δR(R+1) represents the ratio of difference between collection vector intensities DR and DR+1, and DR. In various embodiments of the invention, the average z-crosstalk error values eN and eS, and δR(R+1) are used for correcting the z-crosstalk between elements 1-64.
- D1 for
element 1 is a function of S1, δ1−2, eN and eS is represented by the following equations:
Thus, D1 is written as:
D1=S 1/(1+δ1−2(e N −e S)) (9)
Similarly, D63 forelement 63 is a function of S63, δ63−62, δ63−64, eN and eS is represented by the following equations:
wherein in equation (7), eN is 0 and in equation (10) eS is 0. - Therefore, DR for element R is a function of SR, δR(R−1), δR(R+1), eN and eS is represented by the following equations:
δRR+1 values are calculated using the following equation:
δR(R+1)=(D R −D R+1)/D R (16)
In this case, δR−(R+1)=(DR−DR+1)/DR represents the signal gradient from two contiguous rows. For a flat anatomy where there is no slope, δR−(R+1) will be equal to zero and no crosstalk correction will be performed. However, if there is a slope in the anatomy, leading to different X-ray signals in two contiguous rows, then correction will be performed. Because δR−(R+1) is based on primary beam signals without crosstalk effects, DR, an approximate value of δR−(R+1) can be calculated using the measured intensities SR, which are the X-ray signals that are detected. - Therefore, measured intensities S1-S64 are used to calculate the approximate value of δR(R+1), using the following approximate equation:
δR(R+1)=(S R −S R+1)/S R (17)
Thus, the collection vector intensities D1-D64 calculated by this algorithm provide an approximate correction of the crosstalk betweenphotosensors 706 indetector 320. This algorithm is performed before the negative logarithm process in the pre-process flow, just after offset correction. - The first algorithm and the second algorithm are used for correcting artifacts induced due to z-crosstalk error values between
photosensors 706 indetector 320. Further, these algorithms may be used to correct the x-crosstalk betweenphotosensors 706 indetector 320. - Various embodiments of the invention allow for producing accurate CT images without any artifacts or smudges specially, the CT images objects with steep slopes such as liver and the tip of head.
- Various embodiments of the invention do not need very accurate manual alignment of the scintillator crystals with the photosensors in order to obtain accurate images. In various embodiments, the z-crosstalk error values eN and eS are used to correct the collection vector intensities D1-D64. These error values are used to calibrate the CT system while performing image reconstruction.
- One technical effect of the present invention is that accurate manual alignment of the scintillator crystals with the photosensors is not needed in order to obtain accurate images. Also, z-crosstalk error values eN and eS are used to correct the collection vector intensities D1-D64. These error values are used to calibrate the CT system while performing image reconstruction.
- Although the various embodiments are described with respect to medical imaging, it should be understood that the various embodiments described herein are not limited to medical applications, but may be utilized in non-medical applications.
- As used herein, an element or step recited in the singular and preceded with the word “a” or “an” should be understood as not excluding plural elements or steps, unless such exclusion is explicitly recited. Furthermore, references to “an” or “one” “embodiment” of the present invention are not intended to be interpreted as excluding the existence of additional embodiments that also incorporate the recited features. Also as used herein, the phrase “reconstructing an image” is not intended to exclude embodiments of the present invention in which data representing an image is generated, but a viewable image is not. However, many embodiments generate (or are configured to generate) at least one viewable image.
- The various embodiments of the invention and components thereof may be implemented as part of a computer system. The computer system may include a computer, an input device, a display unit and an interface, for example, for accessing the Internet.
- As used herein, the term “computer” may include any processor-based or microprocessor-based system including systems using microcontrollers, reduced instruction set circuits (RISC), application specific integrated circuits (ASICs), logic circuits, and any other circuit or processor capable of executing the functions described herein. The above examples are exemplary only, and are thus not intended to limit in any way the definition and/or meaning of the term “computer”.
- The computer system executes a set of instructions that are stored in one or more storage elements, in order to process input data. The storage elements may also hold data or other information as desired or needed. The storage element may be in the form of an information source or a physical memory element within the processing machine.
- The set of instructions may include various commands that instruct the processing machine to perform specific operations such as the processes of the various embodiments of the invention. The set of instructions may be in the form of a software program. The software may be in various forms such as system software or application software. Further, the software may be in the form of a collection of separate programs, a program module within a larger program or a portion of a program module. The software also may include modular programming in the form of object-oriented programming. The processing of input data by the processing machine may be in response to user commands, or in response to results of previous processing, or in response to a request made by another processing machine.
- While the invention has been described in terms of various specific embodiments, those skilled in the art will recognize that the invention can be practiced with modification within the spirit and scope of the claims.
Claims (21)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/159,706 US7825370B2 (en) | 2005-06-23 | 2005-06-23 | Method and system for calibrating a computed tomography system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/159,706 US7825370B2 (en) | 2005-06-23 | 2005-06-23 | Method and system for calibrating a computed tomography system |
Publications (2)
Publication Number | Publication Date |
---|---|
US20060289765A1 true US20060289765A1 (en) | 2006-12-28 |
US7825370B2 US7825370B2 (en) | 2010-11-02 |
Family
ID=37566227
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/159,706 Active 2025-08-19 US7825370B2 (en) | 2005-06-23 | 2005-06-23 | Method and system for calibrating a computed tomography system |
Country Status (1)
Country | Link |
---|---|
US (1) | US7825370B2 (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110211667A1 (en) * | 2010-02-26 | 2011-09-01 | Abdelaziz Ikhlef | De-populated detector for computed tomography and method of making same |
EP2251680A3 (en) * | 2009-05-13 | 2012-02-08 | Ishida Co., Ltd. | X-ray inspection device |
US20120166128A1 (en) * | 2010-12-28 | 2012-06-28 | Ikhlef Abdelaziz | Method and apparatus for detector calibration |
CN102917642A (en) * | 2010-06-03 | 2013-02-06 | 株式会社日立医疗器械 | X-ray ct device |
US9168008B2 (en) | 2011-11-03 | 2015-10-27 | General Electric Company | Coarse segmented detector architecture and method of making same |
GB2528168A (en) * | 2014-05-19 | 2016-01-13 | Diamond Light Source Ltd | Analysis of signals from pixellated detectors of ionizing radiation |
US20160231441A1 (en) * | 2015-02-05 | 2016-08-11 | Kabushiki Kaisha Toshiba | Radiation detection device, radiation detection method, and computer program product |
US10473796B2 (en) | 2018-02-07 | 2019-11-12 | FMI Medical Systems Co., Ltd. | Scintillating array with alignment features |
US10531848B2 (en) * | 2018-05-16 | 2020-01-14 | FMI Medical Systems Co., Ltd. | Packaging for CT detector |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8988660B2 (en) * | 2011-06-29 | 2015-03-24 | Silicon Laboratories Inc. | Optical detector |
US10433811B2 (en) | 2016-12-23 | 2019-10-08 | General Electric Company | Self-calibrating CT detectors, systems and methods for self-calibration |
CN107271042A (en) * | 2017-04-27 | 2017-10-20 | 中国空空导弹研究院 | A kind of detector crosstalk measuring device, its preparation method and crosstalk measuring method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4897788A (en) * | 1988-04-18 | 1990-01-30 | General Electric Company | Image correction for computed tomography to remove crosstalk artifacts |
US5764720A (en) * | 1996-12-26 | 1998-06-09 | General Electric Company | Methods and apparatus for simplified pre-processing of data in a computed tomography system |
US5828719A (en) * | 1996-12-23 | 1998-10-27 | General Electric Company | Methods and apparatus for modulating data acquisition system gain |
US20070023668A1 (en) * | 2005-07-28 | 2007-02-01 | Sreeram Dhurjaty | Low noise image data capture for digital radiography |
-
2005
- 2005-06-23 US US11/159,706 patent/US7825370B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4897788A (en) * | 1988-04-18 | 1990-01-30 | General Electric Company | Image correction for computed tomography to remove crosstalk artifacts |
US5828719A (en) * | 1996-12-23 | 1998-10-27 | General Electric Company | Methods and apparatus for modulating data acquisition system gain |
US5764720A (en) * | 1996-12-26 | 1998-06-09 | General Electric Company | Methods and apparatus for simplified pre-processing of data in a computed tomography system |
US20070023668A1 (en) * | 2005-07-28 | 2007-02-01 | Sreeram Dhurjaty | Low noise image data capture for digital radiography |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2251680A3 (en) * | 2009-05-13 | 2012-02-08 | Ishida Co., Ltd. | X-ray inspection device |
US20110211667A1 (en) * | 2010-02-26 | 2011-09-01 | Abdelaziz Ikhlef | De-populated detector for computed tomography and method of making same |
US9482629B2 (en) | 2010-06-03 | 2016-11-01 | Hitachi, Ltd. | X-ray CT apparatus |
CN102917642A (en) * | 2010-06-03 | 2013-02-06 | 株式会社日立医疗器械 | X-ray ct device |
US20120166128A1 (en) * | 2010-12-28 | 2012-06-28 | Ikhlef Abdelaziz | Method and apparatus for detector calibration |
US9168008B2 (en) | 2011-11-03 | 2015-10-27 | General Electric Company | Coarse segmented detector architecture and method of making same |
US10168438B2 (en) | 2014-05-19 | 2019-01-01 | Diamond Light Source Limited | Analysis of signals from pixellated detectors of ionizing radiation |
GB2528168A (en) * | 2014-05-19 | 2016-01-13 | Diamond Light Source Ltd | Analysis of signals from pixellated detectors of ionizing radiation |
GB2528168B (en) * | 2014-05-19 | 2021-09-22 | Diamond Light Source Ltd | Analysis of signals from pixellated detectors of ionizing radiation |
US20160231441A1 (en) * | 2015-02-05 | 2016-08-11 | Kabushiki Kaisha Toshiba | Radiation detection device, radiation detection method, and computer program product |
US9864079B2 (en) * | 2015-02-05 | 2018-01-09 | Kabushiki Kaisha Toshiba | Radiation detection device, radiation detection method, and computer program product |
US10473796B2 (en) | 2018-02-07 | 2019-11-12 | FMI Medical Systems Co., Ltd. | Scintillating array with alignment features |
US10531848B2 (en) * | 2018-05-16 | 2020-01-14 | FMI Medical Systems Co., Ltd. | Packaging for CT detector |
Also Published As
Publication number | Publication date |
---|---|
US7825370B2 (en) | 2010-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7825370B2 (en) | Method and system for calibrating a computed tomography system | |
CA2786110C (en) | Method and apparatus for acquiring radiation data | |
JP5220617B2 (en) | Random reduction via TOFFOV | |
US7872221B2 (en) | Apparatus and methods for calibrating pixelated detectors | |
CN103998952B (en) | X-ray detector | |
US6624422B2 (en) | Method for dynamic stabilization of PET detector gains | |
US9029786B2 (en) | Nuclear medicine imaging apparatus, and nuclear medicine imaging method | |
EP2867701B1 (en) | Digital positron emission tomography (dpet) energy calibration method | |
US7632016B1 (en) | Digital detector calibration with known exposure | |
US8699659B2 (en) | Systems and methods for focal spot motion correction | |
EP1151322B1 (en) | X-ray imaging apparatus and method with scatter correction | |
JPS62161348A (en) | X-ray photographing method and apparatus | |
US5999588A (en) | Deadtime correction in a nuclear medicine imaging system | |
JP2011503535A (en) | Indirect radiation detector | |
US6100531A (en) | Dual-purpose radiation transmission source for nuclear medicine imaging system | |
US6259766B1 (en) | Computer tomography device | |
US6839401B2 (en) | X-ray computed tomography apparatus | |
US7056020B2 (en) | Alignment systems and methods for radiographic imaging systems | |
US7362844B2 (en) | Tomography appliance, and method for a tomography appliance | |
US5998793A (en) | Method and apparatus for correcting for random coincidences in a nuclear medicine imaging system | |
US6140649A (en) | Imaging attenuation correction employing simultaneous transmission/emission scanning | |
JPH10234724A (en) | X-ray computed tomograph | |
US6249003B1 (en) | Imaging attenuation correction method employing multiple energy scan masks and windows | |
US8415632B2 (en) | Methods and systems for calibrating a nuclear medicine imaging system | |
JP2013009874A (en) | X-ray ct apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: GENERAL ELECTRIC COMPANY, NEW YORK Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:IKHLEF, ABDELAZIZ;IMAI, YASUHIRO;SIGNING DATES FROM 20050613 TO 20050621;REEL/FRAME:016724/0766 Owner name: GENERAL ELECTRIC COMPANY, NEW YORK Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:IKHLEF, ABDELAZIZ;IMAI, YASUHIRO;REEL/FRAME:016724/0766;SIGNING DATES FROM 20050613 TO 20050621 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1552) Year of fee payment: 8 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 12TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1553); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 12 |