US10213177B2 - X-ray CT apparatus, data processing device, and projection data generation method - Google Patents
X-ray CT apparatus, data processing device, and projection data generation method Download PDFInfo
- Publication number
- US10213177B2 US10213177B2 US15/321,903 US201515321903A US10213177B2 US 10213177 B2 US10213177 B2 US 10213177B2 US 201515321903 A US201515321903 A US 201515321903A US 10213177 B2 US10213177 B2 US 10213177B2
- Authority
- US
- United States
- Prior art keywords
- data
- sensitivity variation
- ray
- noise
- measurement data
- 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.)
- Active, expires
Links
- 238000000034 method Methods 0.000 title claims abstract description 172
- 238000012545 processing Methods 0.000 title claims description 70
- 230000035945 sensitivity Effects 0.000 claims abstract description 182
- 230000008569 process Effects 0.000 claims abstract description 133
- 238000005259 measurement Methods 0.000 claims abstract description 111
- 238000012937 correction Methods 0.000 claims abstract description 76
- 238000001514 detection method Methods 0.000 claims abstract description 44
- 239000000284 extract Substances 0.000 claims abstract description 21
- 238000000605 extraction Methods 0.000 claims abstract description 21
- 238000011946 reduction process Methods 0.000 claims description 42
- 238000009499 grossing Methods 0.000 claims description 13
- 230000009467 reduction Effects 0.000 claims description 9
- 238000011156 evaluation Methods 0.000 claims description 8
- 238000004846 x-ray emission Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000004904 shortening Methods 0.000 claims description 2
- 238000002591 computed tomography Methods 0.000 description 30
- 238000003860 storage Methods 0.000 description 21
- 230000009466 transformation Effects 0.000 description 13
- 230000006870 function Effects 0.000 description 12
- 238000013480 data collection Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 3
- NRNCYVBFPDDJNE-UHFFFAOYSA-N pemoline Chemical compound O1C(N)=NC(=O)C1C1=CC=CC=C1 NRNCYVBFPDDJNE-UHFFFAOYSA-N 0.000 description 3
- 239000004243 E-number Substances 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 101001109518 Homo sapiens N-acetylneuraminate lyase Proteins 0.000 description 1
- 102100022686 N-acetylneuraminate lyase Human genes 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000011895 specific detection Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
Definitions
- the present invention relates to an X-ray computed tomography (CT) apparatus, a data processing device, and a projection data generation method.
- CT computed tomography
- the present invention relates to a noise reduction process of measurement data obtained by the X-ray CT apparatus.
- the X-ray CT apparatus emits X-ray using multiple angles around an object, and causes an X-ray detector to acquire the X-ray attenuated by the object.
- An X-ray beam is shaped in such a way that the X-ray generated in an X-ray tube is transmitted through a bowtie filter, and thereafter, the X-ray beam is emitted to the object.
- the X-ray detector is configured so that multiple detector blocks having multiple X-ray detection elements are arrayed in parallel in a rotation direction (hereinafter, referred to as a channel direction) and a rotation axis direction (column direction) of the apparatus.
- the X-ray detector having the multiple detector blocks In a case where the X-ray detector having the multiple detector blocks is used in this way, a clearance is generated between the adjacent detector blocks. Due to this clearance, sensitivity variations occur between the detection element on an inner side of the detector blocks and the detection element on an, end side of the detector blocks. In another case, due to variations in a light receiving area which occur during a manufacturing process of the X-ray detector, the sensitivity variations occur between the respective detection elements. In particular, the sensitivity variations caused by the former factor always occur in a specific detection element, thereby causing a ring-like artifact to appear on an image. In order to restrain these sensitivity variations, the X-ray CT apparatus in the related art performs air correction on measurement data acquired by the X-ray detector.
- air data scanned in a state where the object is not present is measured in advance, and the air data is referred to as sensitivity of the detection element, thereby correcting a sensitivity variation component included in the measurement data scanned in a state where the object is present.
- a representative method thereof includes a noise reduction method disclosed in PTL 1.
- PTL 1 discloses a method of reducing noise included in measurement data by modelling statistical fluctuations of the measurement data which are caused by an X-ray detection process.
- the present invention is made in view of the above-described problems, and an object thereof is to provide an X-ray CT apparatus, a data processing device, and a projection data generation method, which can generate projection data for obtaining an image with satisfactory image quality by accurately removing the influence of sensitivity variations in an X-ray detector which are included in measurement data obtained by the X-ray CT apparatus scanning an object.
- an X-ray CT apparatus includes a scan gantry unit that emits X-ray to an object, and that detects the X-ray transmitted through the object, a bed on which the object is laid, and in which the object is loaded into or unloaded from an X-ray emission range of the scan gantry unit, and an operation desk including a data processing device that controls each portion of the scan gantry unit, that acquires measurement data measured by the scan gantry unit, and that generates an image including an object tomographic image from the measurement data.
- the data processing device includes an extraction unit that acquires air data measured using the X-ray CT apparatus and the measurement data obtained by scanning the object, that extracts sensitivity variation data which is a sensitivity variation component of a detection element from the air data, and that extracts blank data from which the sensitivity variation component is removed, and a projection data generation unit that removes the sensitivity variation component and noise which are included in the measurement data, based on the sensitivity variation data, and that uses the blank data so as to perform a correction process of the measurement data from which the sensitivity variation component and the noise are removed.
- a data processing device acquires air data measured using an X-ray CT apparatus and measurement data obtained by scanning an object.
- the data processing device includes an extraction unit that that extracts sensitivity variation data which is a sensitivity variation component of a detection element from the air data, and that extracts blank data from which the sensitivity variation component is removed, and a projection data generation unit that removes the sensitivity variation component and noise which are included in the measurement data, based on the sensitivity variation data, and that uses the blank data so as to perform a correction process of the measurement data from which the sensitivity variation component and the noise are removed.
- a projection data generation method includes a step of causing a data processing device to acquire air data measured using an X-ray CT apparatus and measurement data obtained by scanning an object, to extract sensitivity variation data which is a sensitivity variation component of a detection element from the air data, and to extract blank data from which the sensitivity variation component is removed, a step of removing the sensitivity variation component and noise which are included in the measurement data, based on the sensitivity variation data, and a step of performing a correction process using the blank data on the measurement data from which the sensitivity variation component and the noise are removed.
- an X-ray CT apparatus a data processing device, and a projection data generation method, which can generate projection data for obtaining an image with satisfactory image quality by accurately removing the influence of sensitivity variations in an X-ray detector which are included in measurement data obtained by the X-ray CT apparatus scanning an object.
- FIG. 1 is an overall configuration diagram of an X-ray CT apparatus 1 .
- FIG. 2 is a layout view when adjacent detector blocks B 1 and B 2 are viewed in a column direction.
- FIG. 3 is a view when the adjacent detector blocks B 1 and B 2 are viewed in an X-ray emission direction.
- FIG. 4 is a view illustrating a function of a sensitivity variation extraction unit 31 .
- FIG. 5 is a view illustrating a function of a projection data generation unit 32 .
- FIG. 6 is a flowchart illustrating a procedure in a sensitivity variation extraction process.
- FIG. 7( a ) illustrates air data 40
- FIG. 7( b ) illustrates sensitivity variation data 41
- FIG. 7( c ) illustrates blank data 42 obtained by performing an end element interpolation process and a smoothing process.
- FIG. 8 is a flowchart illustrating a procedure in a projection data generation process.
- FIG. 9( a ) illustrates measurement data 50
- FIG. 9( b ) illustrates sensitivity variation correction-processed measurement data 51 .
- FIG. 10 is a flowchart illustrating a procedure in a projection data generation process according to a second embodiment.
- FIG. 11 is a flowchart illustrating a procedure in a projection data generation process according to a third embodiment.
- FIG. 1 is an overall configuration diagram of the X-ray CT apparatus 1 .
- the X-ray CT apparatus 1 includes a scan gantry unit 100 , a bed 106 , and an operation desk 120 .
- the scan gantry unit 100 is a device that emits X-ray to an object, and that detects the X-ray transmitted through the object.
- the operation desk 120 is a device that controls each portion of the scan gantry unit 100 , that acquires transmitted X-ray data (measurement data) measured by the scan gantry unit 100 , and that generates an image such as an object tomographic image from the measurement data.
- the bed 106 is a device on which the object is laid, and in which the object is loaded into or unloaded from an X-ray emission range of the scan gantry unit 100 .
- the scan gantry unit 100 includes an X-ray source 101 , a rotary disc 102 , a bowtie filter 103 , a collimator 104 , an X-ray detector 107 , a data collection device 108 , a gantry control device 109 , a bed control device 110 , and an X-ray control device 111 .
- the operation desk 120 includes an input device 121 , a data processing device (image processing device) 122 , a storage device 123 , a system control device 124 , and a display device 125 .
- An opening portion 105 is disposed in the rotary disc 102 of the scan gantry unit 100 .
- the X-ray source 101 and the X-ray detector 107 are arranged to face each other via the opening portion 105 .
- the object laid on the bed 106 is inserted into the opening portion 105 .
- the rotary disc 102 is rotated around the object by a drive force transmitted through a drive transmission system from a rotary disc drive device.
- the rotary disc drive device is controlled by the gantry control device 109 .
- the X-ray source 101 is controlled by the X-ray control device 111 , and continuously or intermittently emits X-ray having predetermined intensity.
- the X-ray control device 111 controls an X-ray tube voltage and an X-ray tube current to be applied or supplied to the X-ray source 101 , in accordance with the X-ray tube voltage and the X-ray tube current which are determined by the system control device 124 of the operation desk 120 .
- the bowtie filter 103 and the collimator 104 are disposed in an X-ray emission port of the X-ray source 101 .
- the bowtie filter 103 is a device for adjusting the X-ray emitted from the X-ray source 101 so as to be incident on the X-ray detector 107 with uniform intensity.
- the collimator 104 is a device for limiting an emission range of the X-ray which is emitted from an X-ray tube 101 and which is adjusted by the bowtie filter 103 .
- the collimator 104 is shaped in a cone beam (conical or pyramidal beam). An opening width of the collimator 104 is controlled by the system control device 124 .
- the X-ray which is emitted from the X-ray source 101 , which passes through the bowtie filter 103 and the collimator 104 , and which is transmitted through the object is absorbed (attenuated) in each tissue inside the object, passes through the object, and is incident on the X-ray detector 107 .
- an X-ray detection element group configured to have a combination of a scintillator and a photodiode is two-dimensionally arrayed in a channel direction (circumferential direction) and a column direction (body axis direction).
- a structure of the X-ray detector 107 will be described later (refer to FIGS. 2 and 3 ).
- the X-ray detector 107 is disposed so as to face the X-ray source 101 via the object.
- the X-ray detector 107 detects the amount of the X-ray transmitted through the object, and outputs the amount of the X-ray to the data collection device 108 .
- the data collection device 108 collects the amount of the X-ray detected by each X-ray detection element of the X-ray detector 107 at a predetermined sampling interval, and converts the amount of the X-ray into a digital signal. As the measurement data, the data collection device 108 sequentially outputs the digital signal to the image processing device 122 of the operation desk 120 . In order to correct the sensitivity of the X-ray detector 107 and the like, the data collection device 108 uses air data. The air data is the measurement data scanned in a state where the object is not present.
- the image processing device 122 acquires the measurement data input from the data collection device 108 , and performs a predetermined preprocessing on the measurement data, thereby generating projection data required for reconfiguration.
- the image processing device 122 uses the generated projection data so as to reconfigure an object image such as a tomographic image.
- the system control device 124 causes the storage device 123 to store object image data reconfigured by the image processing device 122 , and causes the display device 125 to display the object image data.
- the image processing device 122 may cause the storage device 123 to store the measurement data input from the data collection device 108 .
- the image processing device 122 may perform a projection data generation process after reading the measurement data at any optional timing different from scanning timing.
- output offset correction is generally performed using a dark current of the X-ray detector 107 .
- a noise reduction process is performed using an inter-proximity element smoothing process or the like.
- An air correction process is performed using the air data on noise-reduced data.
- a reference correction process, a logarithm transformation process, and a phantom correction process for restraining a beam hardening effect, and the like are performed by using a value of a reference detector, thereby generating the projection data.
- the image processing device 122 extracts the sensitivity variations between the detection elements which are included in the air data, as sensitivity variation data.
- the image processing device 122 extracts and stores blank data obtained by removing a sensitivity variation component from the air data.
- the image processing device 122 performs the noise reduction process. Thereafter, the image processing device 122 performs a correction process using the above-described blank data on the measurement data from which the sensitivity variation component and the noise are removed. Thereafter, similarly to the projection data generation process in the related art, the reference correction process, the logarithm transformation process, and the phantom correction process and the like are performed so as to generate the projection data. A process of separating and extracting the sensitivity variation data and the blank data from the air data, and the projection data generation process will be described in detail later.
- the system control device 124 is a computer including a central processing unit (CPU), a read only memory (ROM), and random access memory (RAM) and the like.
- the storage device 123 is a data recording device such as a hard disk, and stores in advance a program or data or the like for realizing a function of the X-ray CT apparatus 1 .
- the display device 125 is configured to include a display device such as a liquid crystal panel and a CRT monitor, and a logic circuit for performing a display process in cooperation with the display device, and is connected to the system control device 124 .
- the display device 125 displays the object image output from the image processing device 122 , and various information items handled by the system control device 124 .
- the input device 121 is configured to include a keyboard, a pointing device such as a mouse, ten keys, and various switch buttons and the like, and outputs various instructions or information items input by an operator to the system control device 124 .
- the operator interactively operates the X-ray CT apparatus 1 by using the display device 125 and the input device 121 .
- the input device 121 may be a touch panel-type input device configured to be integral with a display screen of the display device 125 .
- the bed 106 includes a top plate on which the object is laid, a vertical moving device, and a top plate drive device.
- a height of the top plate is vertically raised and lowered under the control of the bed control device 110 .
- the top plate is longitudinally moved in a body axis direction, or laterally moved in a direction perpendicular to the body axis and in a direction parallel to a floor surface (in the rightward and leftward direction in the drawing).
- the bed control device 110 moves the top plate at a bed moving speed and in a bed moving direction which are determined by the system control device 124 .
- FIG. 2 is a view when an arrangement of two adjacent detector blocks B 1 and B 2 in the X-ray detector 107 are viewed in the column direction.
- FIG. 3 is a view when viewed in the X-ray emission direction.
- the respective detector blocks B 1 and B 2 are configured so that multiple detection elements d having a scintillator element 22 and a photoelectric conversion element 23 are arrayed in the channel direction and the column direction.
- scattered X-ray collimators 21 are disposed parallel to each other in the channel direction and the column direction so that the X-ray scattered in the object is less likely to be incident on the detection elements d.
- the detector blocks B 1 and B 2 are arranged on an arc along the rotary disc 102 . Therefore, as illustrated in FIGS. 2 and 3 , a clearance C is generated between the adjacent detector blocks B 1 and B 2 . Due to the clearance C, sensitivity variations occur in an end element and an internal element of the respective detector blocks B 1 and B 2 in the X-ray detector 107 . These cause the ring-like artifact on the image.
- the image processing device 122 of the X-ray CT apparatus 1 includes a sensitivity variation extraction unit 31 and a projection data generation unit 32 .
- the sensitivity variation extraction unit 31 and the projection data generation unit 32 in the image processing device 122 according to the present invention will be described.
- the sensitivity variation extraction unit 31 acquires air data 40 which is the measurement data obtained by scanning in a state where the object is not present, and extracts a sensitivity variation component between the detection elements of the X-ray detector 107 which is included in the air data 40 , as sensitivity variation data 41 .
- the sensitivity variation extraction unit 31 removes the sensitivity variation component from the air data 40 , and extracts blank data 42 .
- the sensitivity variation data 41 and the blank data 42 which are extracted are stored in the storage device 123 (refer to FIG. 7 ).
- the air data 40 is data prior to the logarithm transformation, in which an offset correction process and a reference correction process are applied to the measurement data obtained in a state where the object is not present. It is desirable to periodically scan the air data 40 .
- the air data 40 is scanned under each scanning condition such as a scanner's rotation speed and a collimator's opening width, and each air data 40 is stored in the storage device 123 .
- the projection data generation unit 32 acquires measurement data 50 obtained by scanning the object, and performs an offset correction process 32 a . Thereafter, referring to the above-described sensitivity variation data 41 , the projection data generation unit 32 performs a sensitivity variation correction process 32 b . Subsequently, the projection data generation unit 32 applies a noise reduction process 32 c to the measurement data 51 (refer to FIG. 9 ) from which the sensitivity variation component is removed by performing the sensitivity variation correction process 32 b . Furthermore, a blank correction process 32 d is performed using the above-described blank data 42 on the measurement data from which the sensitivity variation component and the noise are removed.
- the projection data generation unit 32 applies a reference correction process 32 e , a logarithm transformation process 32 f , and a phantom correction process 32 g to the object attenuation data, thereby preparing projection data 53 required for obtaining a tomographic image.
- the sensitivity variation correction process 32 b , the noise reduction process 32 c , and the blank correction process 32 d will be described later.
- the image processing device 122 of the X-ray CT apparatus performs a sensitivity variation extraction process in accordance with a procedure illustrated in a flowchart in FIG. 6 .
- the image processing device 122 uses the sensitivity variation data 41 and the blank data 42 which are extracted by the sensitivity variation extraction process so as to perform the projection data generation process illustrated in FIG. 8 .
- the air data 40 Prior to the sensitivity variation extraction process, the air data 40 is scanned under various scanning conditions, and is stored in the storage device 123 .
- an end element map is prepared in advance, and is stored in the storage device 123 .
- the end element map is prepared by mapping positions (numbers in the channel direction and the column direction) of the detection elements located in end portions of the respective detector blocks B 1 , B 2 , . . . of the X-ray detector 107 .
- numbers of i ⁇ 1, . . . , i, . . . , I ⁇ and j ⁇ 1, . . . , j, . . . , I ⁇ are respectively given in the channel direction and the column direction of the detection element.
- a value of the end element map of the detection element of i-th channel and j-th column is set to b i, j .
- the element located in the end portion is determined in advance in accordance with a clearance size and a channel width of the detector block.
- information of the end element map can be changed in accordance with these changes.
- the image processing device 122 first performs a data reading process (Step S 101 ).
- Step S 101 the image processing device 122 reads the air data 40 and the end element map corresponding to the scanning condition from the storage device 123 .
- a value of the air data 40 in the detection element of i-th channel and j-th column is expressed by p i, j .
- the image processing device 122 performs interpolation and extrapolation on the end element of the detector block (Step S 102 ).
- Step S 102 a method of independently performing the interpolation in the channel direction and the column direction will be described. However, it is also possible to simultaneously performing the interpolation on two dimensions in the channel direction and the column direction.
- the interpolation data p′ i, j _ ch in the channel direction obtained by Expression (1) and Expression (2) above and the original air data 40 are combined with each other, thereby calculating interpolation-purpose data p′ i, j _ sl in the column direction as expressed in Expression (3) below.
- an interpolation coefficient a′ k (i, j) is expressed by Expression (5) below.
- the interpolated data can be obtained in which the sensitivity variations are corrected between the detection elements of the inner portion and the end portion of the detector block.
- an example using the Lagrange interpolation has been described.
- a known interpolation method such as spline interpolation may be used.
- Step S 103 a Gaussian filter is applied in the channel direction and the column direction of the interpolated data (smoothing process; Step S 103 ) so as to obtain the blank data 42 .
- the sensitivity variations between the detection elements which occur due to a size of a light receiving area generated mainly in a detector manufacturing process are corrected.
- FIG. 7( a ) is a view of the air data 40 read in Step S 101 .
- FIG. 7( c ) is a view illustrating an example of the blank data 42 obtained by the end element interpolation process and the smoothing process in Step S 102 and Step S 103 .
- the blank data 42 is data from which the sensitivity variation component is removed and which reflects only the X-ray attenuation obtained using the bowtie filter 103 .
- Step S 103 without being limited to the Gaussian filter, a known smoothing filter such as a movement averaging filter may be used.
- the image processing device 122 obtains a difference between the air data 40 and the blank data 42 in each detection element, and extracts the sensitivity variation data 41 (Step S 104 ).
- FIG. 7( b ) is an example of the sensitivity variation data 41 obtained by the process in Step S 104 .
- the sensitivity variation data 41 is data obtained by extracting only the sensitivity variation component removed from the air data 40 by the process in Step S 102 and Step S 103 .
- the image processing device 122 causes the storage device 123 to store the sensitivity variation data 41 and the blank data 42 (Step S 105 ), and completes the sensitivity variation extraction process.
- the above-described sensitivity variation extraction process may be performed at any timing as long as the process is performed from when the object starts to be scanned until the sensitivity variation correction process (Step S 203 ) is performed in the projection data generation process to be described later (refer to FIG. 8 ).
- the sensitivity variation extraction process may be performed immediately after the air data is scanned, and the data may be stored in the storage device 123 .
- the image processing device 122 acquires the measurement data 50 obtained by scanning the object under a predetermined scanning condition, from the data collection device 108 or the storage device 123 (Step S 201 ). If the measurement data 50 is acquired, the image processing device 122 applies the offset correction process to the acquired measurement data 50 (Step S 202 ). Thereafter, the image processing device 122 performs the sensitivity variation correction process (Step S 203 ).
- the image processing device 122 acquires the sensitivity variation data 41 from the storage device 123 , and divides the measurement data 50 by the sensitivity variation data 41 . In this manner, the image processing device 122 calculates sensitivity variation correction-processed measurement data 51 . In this case, the same sensitivity variation data 41 in a viewing direction is applied to the measurement data 50 of the detection element in any optional channel and column.
- FIG. 9( a ) illustrates an example of the offset correction-processed measurement data 50 in Step S 202 .
- FIG. 9( b ) illustrates an example of the sensitivity variation correction-processed measurement data 51 in Step S 203 .
- the image processing device 122 applies the noise reduction process to the sensitivity variation correction-processed measurement data 51 as illustrated in FIG. 9( b ) (Step S 204 ).
- the noise reduction process an inter-proximity element smoothing process may be performed.
- the image processing device 122 applies the blank correction process to noise reduction-processed measurement data (Step S 205 ).
- the image processing device 122 acquires the blank data 42 from the storage device 123 , and divides the noise reduction-processed measurement data (processed data in Step S 204 ) by the blank data 42 . In this manner, it is possible to obtain X-ray attenuation data of only the object (hereinafter, referred to as object attenuation data) which has no attenuation using the bowtie filter 103 .
- the image processing device 122 applies the reference correction process (Step S 206 ), the logarithm transformation process (Step S 207 ), and the phantom correction process (Step S 208 ) in the related art to the above-described object attenuation data, and obtains the projection data 53 .
- the image processing device 122 causes the storage device 123 to store the generated projection data 53 (Step S 209 ), and completes a series of the projection data generation processes.
- the X-ray CT apparatus 1 includes the scan gantry unit 100 that emits the X-ray to the object, and that detects the X-ray transmitted through the object, the bed 106 on which the object is laid, and in which the object is loaded into or unloaded from the X-ray emission range of the scan gantry unit 100 , and the operation desk 120 including the image processing device 122 that controls each portion of the scan gantry unit 100 , that acquires the measurement data measured by the scan gantry unit 100 , and that generates the image including the object tomographic image from the measurement data.
- the image processing device 122 includes the extraction unit 31 that acquires the air data measured using the X-ray CT apparatus 1 and the measurement data obtained by scanning the object, that extracts the sensitivity variation data which is the sensitivity variation component of the detection element from the air data, and that extracts the blank data from which the sensitivity variation component is removed, and the projection data generation unit 32 that removes the sensitivity variation component and the noise which are included in the measurement data, based on the sensitivity variation data, and that uses the blank data so as to perform the correction process of the measurement data from which the sensitivity variation component and the noise are removed.
- the image processing device 122 includes the sensitivity variation extraction unit 31 that acquires the air data 40 measured using the X-ray CT apparatus 1 and the measurement data obtained by scanning the object, that extracts the sensitivity variation data 41 which is the sensitivity variation component of the detection element from the air data 40 , and that extracts the blank data 40 from which the sensitivity variation component is removed, and the projection data generation unit 32 that removes the sensitivity variation component and the noise which are included in the measurement data, based on the sensitivity variation data 41 , and that uses the blank data 40 so as to perform the correction process of the measurement data from which the sensitivity variation component and the noise are removed.
- the projection data generation method includes the step of causing the data processing device to acquire the air data measured using the X-ray CT apparatus and the measurement data obtained by scanning the object, to extract the sensitivity variation data which is the sensitivity variation component of the detection element from the air data, and to extract the blank data from which the sensitivity variation component is removed, the step of removing the sensitivity variation component and the noise which are included in the measurement data, based on the sensitivity variation data, and the step of performing the correction process using the blank data on the measurement data from which the sensitivity variation component and the noise are removed.
- the image processing device 122 extracts the sensitivity variation component between the detection elements which are included in the air data 40 obtained by scanning of the X-ray CT apparatus 1 , and extracts the blank data 42 from which the sensitivity variation component is removed, from the air data 40 . Then, after removing the sensitivity variation component included in the measurement data 50 obtained by scanning the object, based on the above-described sensitivity variation data 41 , the image processing device 122 performs the noise removal process (for example, the inter-proximity element smoothing process).
- the noise removal process for example, the inter-proximity element smoothing process
- the correction process using the blank data 42 is performed on the measurement data 50 from which the sensitivity variation component and the noise are removed.
- the object attenuation data from which the influence of the X-ray attenuation is removed by the bowtie filter 103 and which indicates the X-ray attenuation of only the object is extracted.
- the reference correction process, the logarithm transformation process, and the phantom correction process and the like are performed on the object attenuation data, thereby preparing the projection data 53 required for preparing the object tomographic image.
- an image is reconfigured using this projection data 53 , it is possible to prepare an image having no influence of the sensitivity variations in the detection element. Accordingly, it is possible to restrain a ring-like artifact from appearing on the image.
- the blank correction process is performed after the noise reduction process. Accordingly, it is possible to consider the X-ray attenuation in the bowtie filter 103 , and it is possible to perform a proper noise reduction process.
- the sensitivity variation correction or the blank correction is performed on the measurement data before the logarithm transformation.
- the sensitivity variation correction or the blank correction may be performed after the logarithm transformation. That is, the measurement data after the logarithm transformation is subtracted from the sensitivity variation data obtained from the air data after the logarithm transformation. In this manner, it is possible to obtain the measurement data after the sensitivity variation correction.
- the projection data generation unit may perform the inter-proximity element smoothing process for removing the noise, and may perform the correction process.
- the inter-proximity element smoothing process for removing the noise may be performed.
- a configuration of the X-ray CT apparatus 1 according to the second embodiment is the same as that according to the first embodiment.
- repeated description will be omitted, and description will be made by giving the same reference numerals to each unit which is the same as that according to the first embodiment.
- the noise reduction process is performed on the data 51 obtained by removing the sensitivity variation component from the offset correction-processed measurement data.
- a sensitivity variation component removal process can be included in the noise reduction process. That is, according to the second embodiment, a successive approximation noise reduction method included in the evaluation function established by using the sensitivity variation data 41 extracted from the air data 40 is employed. In this manner, it is possible to simultaneously perform the removal of the sensitivity variation component and the removal of the noise from the offset correction-processed measurement data 50 .
- FIG. 10 is a flowchart illustrating a procedure in a projection data generation process according to the second embodiment.
- the sensitivity variation data 41 and the blank data 42 which are used in the projection data generation process are separated and extracted from the air data 40 in the same procedure (refer to FIG. 6 ) as that according to the first embodiment, and are stored in the storage device 123 .
- the image processing device 122 acquires the measurement data 50 obtained by scanning the object under a predetermined scanning condition, from the data collection device 108 or the storage device 123 (Step S 301 ). If the measurement data 50 is acquired, the image processing device 122 applies the offset correction process (Step S 302 ). Thereafter, the image processing device 122 performs the noise reduction process in which the sensitivity variations are considered (Step S 303 ).
- Step S 303 The noise reduction process in which the sensitivity variations are considered in Step S 303 will be described.
- the image processing device 122 acquires the offset correction-processed measurement data and the sensitivity variation data 41 .
- the offset processed-measurement data of i-th channel and j-th column is expressed as y i, j .
- noise-reduced measurement data targeted for estimation is expressed as x i, j .
- a probability of detecting the y i, j number of X-ray photons can be modeled using Poisson distribution.
- the simultaneous probability of all detection elements of the X-ray detector 107 is expressed by Expression (7) below.
- Expression (7) is called a likelihood function, and y i, j which maximizes the likelihood function is obtained by using a successive approximation method. In this manner, statistical noise can be effectively reduced.
- NPL 1 discloses a method of maximizing a posterior probability obtained in such a way that the likelihood function is multiplied by a prior probability.
- the sensitivity variations in the detection element are not considered.
- the sensitivity variation data 41 of i-th channel and j-th column is set to z i,j
- the noise-reduced measurement data in a state having no sensitivity variation is set to x′ i, j .
- Expression (8) is for optimization similarly to NPL1, for example.
- the sensitivity variation and noise-reduced data x′ i, j is successively and approximately estimated.
- Step S 303 the blank correction process is performed using the blank data 42 on the sensitivity variation component and noise-reduced measurement data (Step S 304 ).
- the image processing device 122 divides the sensitivity variation component and noise reduction-processed measurement data by the blank data 42 , thereby obtaining the object attenuation data.
- the image processing device 122 applies the reference correction process (Step S 305 ), the logarithm transformation process (Step S 306 ), and the phantom correction process (Step S 307 ) in the related art, to the object attenuation data obtained in the process until Step S 304 , thereby obtaining the projection data 53 .
- the image processing device 122 causes the storage device 123 to store the generated projection data 53 (Step S 308 ), and completes a series of the projection data generation processes.
- the projection data generation unit 32 removes the sensitivity variation component and the noise from the measurement data by employing the successive approximation noise reduction method using the evaluation function including the sensitivity variation data, and perform the correction process.
- the sensitivity variation component and the noise are removed from the measurement data by employing the successive approximation noise reduction method using the evaluation function including the sensitivity variation data.
- the process of removing the sensitivity variation component can be included in the noise reduction process. If this process procedure is employed, even in a case where the sensitivity variation data cannot be accurately extracted from the air data, the noise reduction process can be more accurately performed, compared to the process procedure according to the first embodiment.
- operator's selection is received with regard to whether or not the noise reduction process considering the sensitivity variations is to be performed on the measurement data.
- FIG. 11 is a flowchart illustrating a process flow according to the third embodiment.
- the image processing device 122 acquires the measurement data 50 obtained by scanning the object under a predetermined scanning condition (Step S 401 ), and applies the offset correction process (Step S 402 ). Thereafter, the selection whether or not to perform the noise reduction process on the measurement data is received (Step S 403 ).
- the operator inputs an instruction indicating whether or not to perform the noise reduction process via the input device 121 .
- the selection process in Step S 403 may be performed prior to Step S 401 .
- Step S 404 If an instruction indicating that the noise reduction process is not performed (Step S 403 ; No), the process proceeds to Step S 404 .
- the processes in Step S 404 to Step S 408 are the same as those in the procedure of the projection data generation process in the related art. That is, the projection data is obtained by applying the air correction process using the air data 40 (Step S 404 ), the reference correction process (Step S 405 ), the logarithm transformation process (Step S 406 ), and the phantom correction process (Step S 407 ) to the measurement data 50 obtained after the offset correction process is performed in Step S 402 .
- the image processing device 122 causes the storage device 123 to store the prepared projection data (Step S 408 ), and completes the projection data generation process.
- Step S 403 in a case where an instruction indicating that the noise reduction process is performed is input in the selection process in Step S 403 (Step S 403 ; Yes), the process proceeds to Step S 409 .
- the processes in Step S 409 to Step S 415 are the same as those in the procedure of the projection data generation process (refer to FIG. 8 ) according to the first embodiment. That is, the sensitivity variation correction process is performed on the offset correction-processed measurement data (Step S 409 ). The noise reduction process is performed on the sensitivity variation corrected measurement data 51 (Step S 410 ). Thereafter, the blank correction process is applied (Step S 411 ).
- the image processing device 122 applies the reference correction process (Step S 412 ), the logarithm transformation process (Step S 413 ), and the phantom correction process (Step S 414 ) to the object attenuation data obtained by the blank correction process, thereby obtaining the projection data 53 .
- the image processing device 122 causes the storage device 123 to store the generated projection data 53 (Step S 415 ), and completes a series of the projection data generation processes.
- Step S 409 to Step S 415 may be substituted with the projection data generation process (Step S 303 to Step S 308 in FIG. 10 ) according to the second embodiment.
- the noise reduction process using the evaluation function in which the sensitivity variations are considered is performed.
- the selection unit (input device 121 ) is further provided which receives selection whether or not to perform the noise reduction process from the measurement data.
- the projection data generation unit performs the noise reduction process in a case where it is selected to perform the noise reduction process via the input device 121 .
- the step of selecting whether or not to perform the noise reduction process from the measurement data is further provided.
- the noise reduction process is performed in a case where it is selected to perform the noise reduction process.
- the third embodiment a user's operation can select whether to improve image quality by applying the noise reduction process which needs a calculation time or whether to give priority to shortening the calculation time. Therefore, the third embodiment is practically adopted.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Human Computer Interaction (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
Description
- PTL 1: U.S. Pat. No. 7,676,074
- NPL 1: “Penalized-Likelihood Sinogram Restoration for Computed Tomography”, P. J. La Riviere, J. Bian and P. A. Vargas, IEEE. Trans. Med. Imag., vol. 25, No. 8, 2006
-
- 1 X-RAY CT APPARATUS,
- 100 SCAN GANTRY UNIT,
- 101 X-RAY SOURCE,
- 102 ROTARY DISC,
- 103 BOWTIE FILTER,
- 104 COLLIMATOR,
- 107 X-RAY DETECTOR,
- 120 OPERATION DESK,
- 121 INPUT DEVICE,
- 122 IMAGE PROCESSING DEVICE (DATA PROCESSING DEVICE),
- 123 STORAGE DEVICE,
- 124 SYSTEM CONTROL DEVICE,
- 125 DISPLAY DEVICE,
- B1, B2 DETECTOR BLOCK,
- C CLEARANCE,
- 31 SENSITIVITY VARIATION EXTRACTION UNIT,
- 32 PROJECTION DATA GENERATION UNIT,
- 40 AIR DATA,
- 41 SENSITIVITY VARIATION DATA,
- 42 BLANK DATA,
- 50 MEASUREMENT DATA,
- 51 SENSITIVITY VARIATION CORRECTION-PROCESSED MEASUREMENT DATA,
- 53 PROJECTION DATA
Claims (12)
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014158472 | 2014-08-04 | ||
JP2014-158472 | 2014-08-04 | ||
PCT/JP2015/070954 WO2016021417A1 (en) | 2014-08-04 | 2015-07-23 | X-ray ct device, data processing device, and projection data generation method |
Publications (2)
Publication Number | Publication Date |
---|---|
US20170135664A1 US20170135664A1 (en) | 2017-05-18 |
US10213177B2 true US10213177B2 (en) | 2019-02-26 |
Family
ID=55263689
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/321,903 Active 2036-04-26 US10213177B2 (en) | 2014-08-04 | 2015-07-23 | X-ray CT apparatus, data processing device, and projection data generation method |
Country Status (4)
Country | Link |
---|---|
US (1) | US10213177B2 (en) |
JP (1) | JP6571659B2 (en) |
CN (1) | CN106572833B (en) |
WO (1) | WO2016021417A1 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018227523A1 (en) | 2017-06-16 | 2018-12-20 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for image data processing in computerized tomography |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4654796A (en) * | 1983-08-29 | 1987-03-31 | Hitachi Medical Corporation | System for computing correction factors for sensitivity differences in detector channels of an X-ray tomography apparatus |
JP2004065287A (en) | 2002-08-01 | 2004-03-04 | Hitachi Medical Corp | X-ray ct device |
JP2008113792A (en) | 2006-11-02 | 2008-05-22 | Hitachi Medical Corp | Radiation measuring apparatus |
US7676074B2 (en) * | 2005-09-12 | 2010-03-09 | Purdue Research Foundation | Systems and methods for filtering data in a medical imaging system |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6307912B1 (en) * | 1999-11-29 | 2001-10-23 | General Electric Company | Methods and apparatus for optimizing CT image quality with optimized data acquisition |
JP4175791B2 (en) * | 2001-08-20 | 2008-11-05 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Image generation method and X-ray CT apparatus |
CN102376097A (en) * | 2010-08-25 | 2012-03-14 | 东软飞利浦医疗设备***有限责任公司 | Computed tomography (CT) iteration reconstruction method for asymmetrical detector |
-
2015
- 2015-07-23 JP JP2016540150A patent/JP6571659B2/en active Active
- 2015-07-23 US US15/321,903 patent/US10213177B2/en active Active
- 2015-07-23 WO PCT/JP2015/070954 patent/WO2016021417A1/en active Application Filing
- 2015-07-23 CN CN201580039240.5A patent/CN106572833B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4654796A (en) * | 1983-08-29 | 1987-03-31 | Hitachi Medical Corporation | System for computing correction factors for sensitivity differences in detector channels of an X-ray tomography apparatus |
JP2004065287A (en) | 2002-08-01 | 2004-03-04 | Hitachi Medical Corp | X-ray ct device |
US7676074B2 (en) * | 2005-09-12 | 2010-03-09 | Purdue Research Foundation | Systems and methods for filtering data in a medical imaging system |
JP2008113792A (en) | 2006-11-02 | 2008-05-22 | Hitachi Medical Corp | Radiation measuring apparatus |
Non-Patent Citations (2)
Title |
---|
International Search Report in connection with PCT.JP2015/070954. |
LaRiviere et al. "Penalized-Likelihood Sinogram Restoration for Computed Tomography" IEEE Transactions on Medical Imaging, vol. 25, No. 8 (2006) (Year: 2006). * |
Also Published As
Publication number | Publication date |
---|---|
JP6571659B2 (en) | 2019-09-04 |
CN106572833A (en) | 2017-04-19 |
JPWO2016021417A1 (en) | 2017-05-25 |
WO2016021417A1 (en) | 2016-02-11 |
US20170135664A1 (en) | 2017-05-18 |
CN106572833B (en) | 2019-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11327031B2 (en) | Photon counting X-ray CT apparatus | |
US8494244B2 (en) | System and method for blood vessel stenosis visualization and quantification using spectral CT analysis | |
JP6313168B2 (en) | X-ray CT apparatus, image processing apparatus, and image processing program | |
US10878602B2 (en) | Image reconstruction device, X-ray CT device, and image reconstruction method | |
JP5918374B2 (en) | X-ray CT apparatus and X-ray CT image processing method | |
US7801264B2 (en) | Method for calibrating a dual -spectral computed tomography (CT) system | |
JP6325256B2 (en) | X-ray computed tomography apparatus and medical image processing apparatus | |
JP6603233B2 (en) | Data processing apparatus, X-ray CT apparatus, and reference correction method | |
JP6176828B2 (en) | Image reconstruction apparatus, image reconstruction method, and X-ray computed tomography apparatus | |
JP4474304B2 (en) | Ring artifact removal method and X-ray CT apparatus | |
JP2006026410A (en) | X-ray computed tomographic system and its data correction method | |
US9980694B2 (en) | X-ray CT apparatus and image calculating device for X-ray CT apparatus | |
KR20140142803A (en) | Radiography imaging apparatus and method for generating an radiographic image | |
US10213177B2 (en) | X-ray CT apparatus, data processing device, and projection data generation method | |
JP7106392B2 (en) | Sensitivity correction method and photon counting detector | |
JP5942216B2 (en) | X-ray CT apparatus and image processing apparatus | |
JP4703221B2 (en) | X-ray CT system | |
JP2015167860A (en) | X-ray computer tomography apparatus, positional deviation specification method and positional deviation specification program | |
JP2022013686A (en) | X-ray ct apparatus, data processing method, and program | |
JP6510203B2 (en) | Image processing apparatus, X-ray CT apparatus and program | |
US20230145523A1 (en) | Medical image processing apparatus, x-ray ct apparatus, medical image processing method and non-volatile storage medium storing program | |
US11324470B2 (en) | X-ray diagnostic apparatus, medical image processing apparatus, and medical image processing method | |
US11944484B2 (en) | Material decomposition calibration method and apparatus for a full size photon counting CT system | |
US20230351647A1 (en) | X-ray diagnostic apparatus, medical image processing apparatus, and storage medium | |
JP2016154839A (en) | X-ray ct apparatus and image processing device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: HITACHI, LTD., JAPAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:TAKAHASHI, HISASHI;GOTO, TAIGA;HIROKAWA, KOICHI;REEL/FRAME:040759/0452 Effective date: 20161109 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
AS | Assignment |
Owner name: FUJIFILM HEALTHCARE CORPORATION, JAPAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:HITACHI, LTD.;REEL/FRAME:058026/0559 Effective date: 20211013 |
|
AS | Assignment |
Owner name: FUJIFILM HEALTHCARE CORPORATION, JAPAN Free format text: CORRECTIVE ASSIGNMENT TO CORRECT THE THE PROPERTY AND APPLICATION NUMBERS PREVIOUSLY RECORDED AT REEL: 058026 FRAME: 0559. ASSIGNOR(S) HEREBY CONFIRMS THE ASSIGNMENT;ASSIGNOR:HITACHI LTD.;REEL/FRAME:058917/0853 Effective date: 20211013 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1551); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 4 |