ITERATIVE METHOD FOR MATERIAL DECOMPOSITION CALIBRATION IN A FULL SIZE PHOTON COUNTING COMPUTED TOMOGRAPHY SYSTEM

Abstract

A photon counting computed tomography (CT) method including, but not limited to, receiving a first forward model including a set of first parameters and a set of second parameters corresponding to a plurality of materials and path lengths by scanning a slab at plural tube voltages and plural tube currents of an X-ray tube; evaluating an image quality of a material decomposition image reconstructed by the set of first parameters and the set of second parameters; and updating at least one second parameters from the set of second parameters if the image quality of the material decomposition image does not satisfy a predetermined threshold, wherein the update of the at least one second parameter from the set of second parameters is achieved by updating the at least one second parameter from the set of second parameters to an energy dependent parameter from a constant value. Processing circuitry that is part of a photon counting computed tomography (CT) can implement the method.

Claims

1. A photon counting computed tomography (CT) apparatus, comprising: processing circuitry configured to: receive a first forward model including a set of first parameters and a set of second parameters corresponding to a plurality of materials and path lengths by scanning a slab at plural tube voltages and plural tube currents of an X-ray tube, evaluate an image quality of a material decomposition image reconstructed by the set of first parameters and the set of second parameters, and update at least one second parameters from the set of second parameters if the image quality of the material decomposition image does not satisfy a predetermined threshold, wherein the update of the at least one second parameter from the set of second parameters is achieved by updating the at least one second parameter from the set of second parameters to an energy dependent parameter from a constant value.

2. The photon counting computed tomography (CT) apparatus according to claim 1, wherein the at least one second parameter from the set of second parameters corresponds to a low tube current condition which is less influenced by pile-up than the at least one second parameter before being updated.

3. The photon counting computed tomography (CT) apparatus according to claim 2, wherein the processing circuitry is further configured to update remaining second parameters from the set of second parameters corresponding to higher tube current based on the at least one second parameter from the set of second parameters corresponding to low tube current.

4. The photon counting computed tomography (CT) apparatus according to claim 1, wherein the processing circuitry is further configured to update the forward model based on the updated at least one second parameter.

5. The photon counting computed tomography (CT) apparatus according to claim 4, wherein the processing circuitry is further configured to evaluate an image quality of the material decomposition image reconstructed based on the updated forward model.

6. The photon counting computed tomography (CT) apparatus according to claim 1, wherein the processing circuitry configured to evaluate the image quality of the material decomposition image reconstructed by the set of first parameters and the set of second parameters comprises processing circuitry configured to perform a pixel by pixel evaluation.

7. The photon counting computed tomography (CT) apparatus according to claim 1, wherein the set of second parameters are related to pileup correction.

8. The photon counting computed tomography (CT) apparatus according to claim 1, wherein the set of first parameters depend on a bin response function.

9. The photon counting computed tomography (CT) apparatus according to claim 1, wherein the processing circuitry configured to receive the first forward model comprises processing circuitry configured to scan a plurality of slabs at the plural tube voltages and the plural tube currents of an X-ray tube.

10. The photon counting computed tomography (CT) apparatus according to claim 1, wherein the energy dependent parameter updated from a constant value is updated by setting all count rate terms to zero.

11. A photon counting computed tomography (CT) method, comprising: receiving a first forward model including a set of first parameters and a set of second parameters corresponding to a plurality of materials and path lengths by scanning a slab at plural tube voltages and plural tube currents of an X-ray tube; evaluating an image quality of a material decomposition image reconstructed by the set of first parameters and the set of second parameters; and updating at least one second parameters from the set of second parameters if the image quality of the material decomposition image does not satisfy a predetermined threshold, wherein the update of the at least one second parameter from the set of second parameters is achieved by updating the at least one second parameter from the set of second parameters to an energy dependent parameter from a constant value.

12. The method according to claim 11, wherein the at least one second parameter from the set of second parameters corresponds to a low tube current condition which is less influenced by pile-up than the at least one second parameter before being updated.

13. The method according to claim 12, further comprising updating remaining second parameters from the set of second parameters corresponding to higher tube current based on the at least one second parameter from the set of second parameters corresponding to low tube current.

14. The method according to claim 11, further comprising updating the forward model based on the updated at least one second parameter.

15. The method according to claim 14, further comprising evaluating an image quality of the material decomposition image reconstructed based on the updated forward model.

16. The method according to claim 11, evaluating the image quality of the material decomposition image reconstructed by the set of first parameters and the set of second parameters comprises performing a pixel by pixel evaluation.

17. The method according to claim 11, wherein the set of second parameters are related to pileup correction.

18. The method according to claim 11, wherein the set of first parameters depend on a bin response function.

19. A non-transitory computer-readable storage medium storing computer-readable instructions that, when executed by a computer, cause the computer to perform a method for calibrating a photon counting computed tomography (CT)apparatus, the method comprising: receiving a first forward model including a set of first parameters and a set of second parameters corresponding to a plurality of materials and path lengths by scanning a slab at plural tube voltages and plural tube currents of an X-ray tube, evaluating an image quality of a material decomposition image reconstructed by the set of first parameters and the set of second parameters, and updating at least one second parameters from the set of second parameters if the image quality of the material decomposition image does not satisfy a predetermined threshold, wherein the update of the at least one second parameter from the set of second parameters is achieved by updating the at least one second parameter from the set of second parameters to an energy dependent parameter from a constant value.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0010] The application will be better understood in light of the description which is given in a non-limiting manner, accompanied by the attached drawings in which:

[0011] FIG. 1 shows an example of a PCD bin response function S.sub.b(E) for a photon counting detector. Each curve stands for an example function for each energy bin.

[0012] FIG. 2 shows an iterative material decomposition calibration and processing workflow.

[0013] FIG. 3 shows normalized linear attenuation coefficients for different materials.

[0014] FIG. 4 is a block diagram of an exemplary configuration of an X-ray CT apparatus imaging a person as a subject according to an exemplary embodiment described below

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0015] Reference throughout this specification to “one embodiment” or “an embodiment” means that a particular feature, structure, material, or characteristic described in connection with the embodiment is included in at least one embodiment of the application, but do not denote that they are present in every embodiment.

[0016] Thus, the appearances of the phrases “in one embodiment” or “in an embodiment” in various places throughout this specification are not necessarily referring to the same embodiment of the application. Furthermore, the particular features, structures, materials, or characteristics may be combined in any suitable manner in one or more embodiments.

[0017] In a transmission measurement using a photon counting energy-resolving detector (PCD), the forward model can be formulated as below:

N.sub.b(I.sub.1, . . . , M)=N.sub.0×∫dE w(E)S.sub.b(E)exp(−Σμ.sub.ml.sub.m),   (1)

where S.sub.b(E) denotes the bin response function defined as S.sub.b(E)=∫.sub.EbL.sup.EbHR(E, E′)dE′ where R(E,E′) is the detector response function, and E.sub.bL and E.sub.bH are the low and high energy threshold of each counting bin. FIG. 1 shows an example model of a typical S.sub.b(E) function for a PCD, where a long tail above the energy window is induced by charge sharing, k-escape and scattering effect. The low energy tail is mostly due to the finite energy resolution from the associated electronic noise. N.sub.0 is the total flux from an air scan, μ.sub.m and l.sub.m are the m.sup.t.sup.h basis material linear attenuation coefficient and path length. w(E) is the normalized incident X-ray spectrum. In practice, both w(E) and S.sub.b(E) are not exactly known, and they can be combined as one term, S.sub.wb(E)=w(E)S.sub.b(E), called thereafter the weighted bin response function. If S.sub.wb(E) can be calibrated through measurements, the decomposition problem at low flux condition can be well solved.

[0018] For a high flux scan condition (e.g., a few percent of pulse pileup), pulse pileup introduces additional spectral distortion in the measurement. One way to correct for the pileup effect is to introduce additional correction terms (e.g. see Dickmann above, who uses the measured count rate(s) as input). And this type of additional calibration is based on an accurate estimation of the flux independent weighted bin response S.sub.wb(E). This leads to the first two steps in the scheme using the following forward model:

N.sub.b(l.sub.1, . . . , M)=N.sub.0∫dE S.sub.wb(E)*P.sub.b(E, N.sub.b, N.sub.tot)exp(−Σμ.sub.ml.sub.m)

[0019] The concept of this two-step calibration method requires that the low flux response term S.sub.wb is flux independent. In a typical CT system, the minimum flux produced by the X-ray tube sometimes is still not low enough to meet the “pileup-free” requirement for estimation of S.sub.wb. Especailly, with the usage of bowtie filters, the center detector pixels will experience much higher flux exposure with air or thin path lengths compared to the peripheral pixels. With inaccurate S.sub.wb, it becomes more difficult to estimate the pileup term P.sub.b with high flux measurements, and introduces potential inconsistency of the calibration results across the detector channels, and at different flux conditions. Therefore, how to compensate the non-ideal measurements due to the system flux limit and improve the estimation of S.sub.wb, then improve the overall calibration quality and stability, is what we try to solve in this proposal.

[0020] In one non-limiting embodiment, an improved calibration builds upon the two-step calibration method for the PCD forward model for material decomposition. The method consists of two parts: 1) estimation of the flux independent weighted bin response function S.sub.wb(E) using the expectation maximization (EM) method, and 2) estimation of the pileup correction term P.sub.b(E, N.sub.b, N.sub.tot) which is a function of energy (E) and the measured bin counts (N.sub.b, N.sub.tot), where N.sub.b is the individual bin count and N.sub.tot is the total count of all the energy bins. The calibrated forward model can be expressed as:

N.sub.b(l.sub.1, . . . , M)=N.sub.0∫.sup.EmaxdES.sub.wb(E)*P.sub.b(E, N.sub.b, N.sub.tot)exp(−Σμ.sub.ml.sub.m)   (2)

[0021] Plural different materials (e.g., 2-5 different materials such as polypropylene, water, aluminium, titanium/copper, and k-edge materials) can be used to calibrate the weighted bin response function S.sub.wb(E) at low flux. With more selective materials used in the calibration, the number of total path lengths is reduced to achieve equivalent or better results.

[0022] Step 1: With an appropriate tube spectrum modelling to capture the characteristic peaks in the incident spectrum, and a physical model to simulate the photon-counting detector spectral response, an initial guess of S.sub.wb(E) can be produced. By using the EM method (e.g., see Sidky et al.), S.sub.wb(E) can be reliably estimated for this very ill-conditioned problem based on a few transmission measurements.

[0023] Here, P.sub.b(E, N.sub.b, N.sub.tot) is assumed to be constant in Step 1. The calibrated forward model can be simplified to a system of linear equations

N.sub.b(l.sub.1, . . . , M)=N.sub.0∫.sup.EmaxdES.sub.wb(E)exp(−Σμ.sub.ml.sub.m)   (3)

[0024] Usually, the number of data measurements (M) is much smaller than the number of unknowns (E.sub.max). With the assumption of Poisson distribution of the data acquisition, an iterative EM algorithm can be derived to find the optimal estimation of the unknown energy bin response function S.sub.wb(E), as described below.

[0025] When estimating the bin response function using low flux data acquisition, the pileup effect correction P.sub.b is assumed to be a known term (e.g., constant). So, the model is simplified to

N.sub.b=N.sub.0∫dE S.sub.wb(E)[exp[−Σμ.sub.m(E)l.sub.m]]  (4)

[0026] Let A.sup.j(E)=exp[−Σμ.sub.m(E)l.sub.m.sup.j] represent the attenuated pathlength for j-th measurement. Thus, for each measurement j, we have

N.sub.b.sup.j=N.sub.0∫dE S.sub.wb(E)A.sup.j(E)=N.sub.0Σ.sub.ES.sub.wb(E)A.sup.j(E)   (5)

[0027] With M measurements, the data acquisition can be written in the matrix form below

[00001]${N_0\left(\begin{array}{ccc}{A}^1\left(1\right)& .Math.& A^1\left(E_{\max }\right)\\ .Math.& \ddots & .Math.\\ A^M\left(1\right)& .Math.& A^M\left(E_{\max }\right)\end{array}\right)}_{M\times E_{\max }}.Math.{\left(\begin{array}{c}\begin{array}{c}{S}_{\mathrm{wb}}\left(1\right)\\ .Math.\end{array}\\ S_{\mathrm{wb}}\left(E_{\max }\right)\end{array}\right)}_{E_{\max }\times 1}={\left(\begin{array}{c}\begin{array}{c}{N}_b^1\\ .Math.\end{array}\\ N_b^M\end{array}\right)}_{M\times 1}$

or A.Math.S.sub.wb=N.sub.b

[0028] By applying the EM iterative algorithm, the S.sub.wb can be estimated by

S.sub.wb.sup.(k+1)=S.sub.wb.sup.(k)⊙((A.sup.T.Math.(N.sub.b(A.Math.S.sub.wb.sup.(k))))(A.sup.T.Math.1))   (6)

where [0029] k: iteration number [0030] .Math.: matrix multiplication [0031] ⊙: element-wise multiplication [0032] : element-wise division [0033] 1: vector of ones with size of M×1
the updating formula for S.sub.wb(E) is given by

[00002]$\begin{array}{cc}{S}_{wb}^{\left(k+1\right)}\left(E\right)=S_{wb}^{\left(k\right)}\left(E\right)\frac{{\Sigma }_{j^{\prime }}{A}^{j^{\prime }}\left(E\right)\frac{N_b^{j^{\prime }}}{{\Sigma }_{E^{\prime }}{A}^{j^{\prime }}\left(E^{\prime }\right)s_{wb}^{\left(k\right)}\left(E^{\prime }\right)}}{{\Sigma }_{j^{\prime }}{A}^{j^{\prime }}\left(E\right)}& \left(7\right)\end{array}$

[0034] Step 2: Once S.sub.wb(E) is estimated from the calibration at each tube voltage (kVp) setting for each detector pixel, it is saved as a calibration table (e.g., in hardware or software) within the system. It will be used as an input to further estimate the pileup correction terms P.sub.b(E, N.sub.b, N.sub.tot) at higher flux scans. Both tables are then used for the material decomposition in object/patient scans to estimate the basis material path lengths.

[0035] The calibration tables are updated from time to time based on the system/detector performance variations. To further improve the calibration quality at high flux condition, an iterative calibration method is proposed based on the two-step calibration framework, by improving the estimation accuracy of the low flux response function S.sub.wb and achieve better overall estimation accuracy at higher flux as well.

[0036] The high-level workflow of the above process is demonstrated in FIG. 2. Steps 1) to 4) represent the calibration workflow, steps 5) to 8) demonstrate how the calibration tables are updated in an iterative process, and step 9) represents the output of the completed calibration tables used in the operational scans of patients/objects to produce spectral images.

[0037] First, Step 1), a series of low and high flux scans on rectangular slabs of different uniform known materials are collected in a photon counting CT system using stationary configuration at each tube kVp setting, which is the peak potential applied on the X-ray tube. Typical CT systems support a few kVp settings from 70 to 140 kVp which generate different energy spectrums from the X-ray tube for different scan protocol. For a CT scan, both mA and kVp need to be pre-selected before the tube is turned on. The photon counting detector generates b bins of counts Nb with pre-set energy thresholds for each of the slabs i. Then, Step 2), the low flux calibration data N.sub.b(1, . . . , i).sub.low is used to estimate a set of first parameters, the low flux response term S.sub.wb at each kVp setting for each detector pixel. Step 3), the low and high calibration data N.sub.b(1, . . . , i).sub.mA is used to estimate a set of second parameters, the pileup term P.sub.b at each mA, including the low mA scans. Step 4), with the estimation calibration tables of S.sub.wb, the set of first parameters, and P.sub.b, the set of second parameters, for each detector pixel, the quality of the calibration is checked on a quality phantom (e.g., a uniform water phantom, or phantom with multiple inserts with uniform known materials). Step 5), the image quality is assessed with predefined standards (e.g., HU at different ROIs, +/−5HU), and if it is passed, the current calibration tables are saved and then used for the following patient/object scans data processing. Otherwise, the procedure goes through the next pass of calibration analysis. Step 6), the second parameters in the set of second parameters is updated by setting the energy dependent parameter of one of the second parameters from the set of second parameters to a constant value. So the pileup reversing term P′.sub.b, the updated second parameter from the set of second parameters, in the low mA scans is calculated by setting the count rate related terms to 0 in P.sub.b(e.g., P.sub.b(E, N.sub.b, N.sub.tot), one will set N.sub.b and N.sub.tot to be ( ) with previously estimated parameters in P.sub.b at that low mA). Then, Step 7) the pseudo pileup-free low flux data N′.sub.b(1, . . . , i).sub.low is calculated based on the previous pass S.sub.wb and P′.sub.b. Step 8), N′.sub.b(1, . . , i).sub.low is then used to estimate the next pass of S.sub.wb. The process is repeated from the second step until the quality check passes for all the mA scans covered by the calibration. Step 9), these calibration tables are saved to be used in the operational scans of patients/objects to produce spectral images. Here, commonly examined standards are: image CT number accuracy, uniformity, spatial resolution, noise and artifacts. To check the quality of this calibration, these metrics should all be checked, especially the accuracy and artifacts like ring or bands in the image, which indicate the calibration is not good enough.

[0038] To choose the optimal materials and path lengths for this calibration, one can use the normalized linear attenuation coefficient vs. energy curves (FIG. 3) to choose the ones that are different from each other, e.g., polypropylene, water, aluminum, and titanium. Such a group of combinations allows calibrations which covers a large range of common materials present in human body.

[0039] In order to satisfy the low flux condition through the calibration measurement to minimize the pileup effect in the flow diagram, step 1, one can select to use nτ<x, where x ˜0.005-0.01 and n is the pixel count rate with the lowest tube flux setting, and τ is the effective dead time of the PCD Application Specific Integrated Circuit (ASIC). These ideal conditions are not always met with the x-ray tube used in the CT system. With bowtie filters, in the center detector region, the air and thin slab measurements can still have a significant portion (e.g. a few percent) of the pileup events that contaminate the estimation of low flux response function S.sub.wb, by forcing it to match the data with pileup effect. Hence, it increases the difficulty of the high flux pileup correction term estimation due to the inaccurate flux-independent term S.sub.wb. The iterative method modifies the two-step method by introducing a pileup-correction step utilizing the previous calibration results.

[0040] The calibration measurements should be taken with sufficient statistics to minimize the influence of the statistical fluctuation. One non-limiting example is to use >1000 times more statistics as the typical integration period for the calibration data sets to minimize the transferred statistical error in the calibration. Each energy bin b of the calibration measurements will be used to update the corresponding S.sub.wb(E) and P.sub.b(E, N.sub.b, N.sub.tot).

[0041] Since one can only do limited number of measurements with a few energy bins, the estimation is very ill-conditioned. In this case, a good initial guess is crucial for an accurate estimation as it provides additional constraints for the EM method. One of the design variations to accommodate non-ideal detectors is to allow a more flexible energy window for each bin in the initial guess of S.sub.b, especially with small variations in the actual energy threshold setting of the ASIC. By setting the low threshold x keV lower, and high threshold y keV higher, the initial S.sub.b becomes:

S.sub.b(E)=∫.sub.EbL−x.sup.EbH+yR(E, E′)dE′  (8)

where x, y can be chosen between 5 to 10 keV to allow certain variations in the ASIC performance, while still providing additional constraints for the EM problem.

[0042] To capture the spectrum variation across the fan beam after bowtie filter and detector response variation across different detector pixels, this calibration process is done pixel by pixel with each bowtie/filter configuration.

[0043] The design described in the present application employs more than two materials in the calibration, which provides more sensitivity to constraint the weighted bin response function estimation problem of the photon counting detectors.

[0044] In addition, the method utilizes the parameterization for the pileup correction terms P.sub.b which is a function of E, N.sub.b and N.sub.tot. The total count term N.sub.tot is introduced for a better approximation of the true pileup phenomena, and can significantly improve the model capability at higher flux condition with fewer parameters.

[0045] Furthermore, various calibration path length ranges are used at different fan angles to improve the calibration accuracy and efficiency. The slab scans used for the forward model calibration can be selected based on the imaging task to generate the best image quality.

[0046] An X-ray CT will be described as an example of a medical image diagnostic modality herein. In that context, an information processing method of information acquired by imaging performed by the X-ray CT will be described.

[0047] The X-ray CT is implemented, for example, in an X-ray CT apparatus 10 illustrated in FIG. 4. FIG. 4 is a block diagram illustrating an example of a configuration of the X-ray CT apparatus 10 according to a first embodiment. For example, the X-ray CT apparatus 10 has a gantry 110, a bed 130, and a console 140.

[0048] In FIG. 4, it is assumed that the longitudinal direction of a rotating shaft of a rotating frame 113 or a tabletop 133 of the bed 130 in a non-tilted state is a Z axis direction. Furthermore, it is assumed that an axial direction orthogonal to the Z axis direction and horizontal to a floor surface is an X axis direction. Furthermore, it is assumed that an axial direction orthogonal to the Z axis direction and perpendicular to the floor surface is a Y axis direction. Note that FIG. 4 illustrates the gantry 110 drawn from a plurality of directions for convenience of description and the X-ray CT apparatus 10 has one gantry 110.

[0049] The gantry 110 includes an X-ray tube 111, an X-ray detector 112, the rotating frame 113, an X-ray high voltage device 114, a control device 115, a wedge 116, a collimator 117, and a data acquisition system (DAS) 118.

[0050] The X-ray tube 111 is a vacuum tube having a cathode (filament) that generates thermoelectrons and an anode (target) that generates X-rays in response to a collision of thermoelectrons. The X-ray tube 111 emits the thermoelectrons toward the anode from the cathode by the application of a high voltage from the X-ray high voltage device 114, thereby generating the X-rays to be emitted to a subject P.

[0051] The X-ray detector 112 detects the X-rays emitted from the X-ray tube 111 and passed through the subject P, and outputs a signal corresponding to the dose of the detected X-rays to the DAS 118. The X-ray detector 112, for example, includes a plurality of detection element arrays in which a plurality of detection elements are arranged in a channel direction (channel direction) along one arc centered on a focal point of the X-ray tube 111. The X-ray detector 112, for example, has a structure in which the detection element arrays with the detection elements arranged in the channel direction are arranged in a row direction (slice direction and row direction).

[0052] For example, the X-ray detector 112 is an indirect conversion type detector having a grid, a scintillator array, and a photosensor array. The scintillator array has a plurality of scintillators. Each of the scintillators has a scintillator crystal that outputs light with a photon quantity corresponding to an incident X-ray dose. The grid has an X-ray shielding plate that is disposed on the surface of the scintillator array on an X-ray incident side and absorbs scatted X-rays. The grid may also be referred to as a collimator (a one-dimensional collimator or a two-dimensional collimator). The photosensor array has a function of converting light into an electrical signal corresponding to the amount of light from the scintillator, and has, for example, photosensors such as photodiodes. Note that the X-ray detector 112 may be a direct conversion type detector having a semiconductor element that converts the incident X-rays into electrical signals.

[0053] The rotating frame 113 is an annular frame that supports the X-ray tube 111 and the X-ray detector 112 so as to face each other and rotates the X-ray tube 111 and the X-ray detector 112 by the control device 115. For example, the rotating frame 113 is a casting made of aluminum. Note that the rotating frame 113 can further support the X-ray high voltage device 114, the wedge 116, the collimator 117, the DAS 118 and the like, in addition to the X-ray tube 111 and the X-ray detector 112. Moreover, the rotating frame 113 can further support various configurations not illustrated in FIG. 4. Hereinafter, in the gantry 110, the rotating frame 113 and a part, which rotationally moves with the rotating frame 113, are also referred to as a rotating part.

[0054] The X-ray high voltage device 114 has electric circuitry such as a transformer and a rectifier, and has a high voltage generation device that generates a high voltage to be applied to the X-ray tube 111 and an X-ray control device that controls an output voltage corresponding to the X-rays generated by the X-ray tube 111. The high voltage generation device may be a transformer type device or an inverter type device. Note that the X-ray high voltage device 114 may be provided on the rotating frame 113, or may also be provided on a fixed frame (not illustrated).

[0055] The control device 115 has processing circuitry having a central processing unit (CPU) and the like, and a driving mechanism such as a motor and an actuator. The control device 115 receives input signals from an input interface 143 and controls the operations of the gantry 110 and the bed 130. For example, the control device 115 controls the rotation of the rotating frame 113, the tilt of the gantry 110, the operation of the bed 130, and the like. As an example, as control for tilting the gantry 110, the control device 115 rotates the rotating frame 113 around an axis parallel to the X axis direction based on information on an input inclination angle (tilt angle). Note that the control device 115 may be provided in the gantry 110 or may also be provided in the console 140.

[0056] The wedge 116 is an X-ray filter for adjusting the dose of the X-rays emitted from the X-ray tube 111. Specifically, the wedge 116 is an X-ray filter that attenuates the X-rays emitted from the X-ray tube 111 such that the X-rays emitted from the X-ray tube 111 to the subject P have a predetermined distribution. For example, the wedge 116 is a wedge filter or a bow-tie filter and is manufactured by processing aluminum and the like to have a predetermined target angle and a predetermined thickness.

[0057] The collimator 117 is a lead plate and the like for narrowing down the emission range of the X-rays having transmitted through the wedge 116 and forms a slit by a combination of a plurality of lead plates and the like. Note that the collimator 117 may also be referred to as an X-ray diaphragm. Furthermore, although FIG. 4 illustrates a case where the wedge 116 is disposed between the X-ray tube 111 and the collimator 117, the collimator 117 may be disposed between the X-ray tube 111 and the wedge 116. In such a case, the wedge 116 attenuates the X-rays, which are emitted from the X-ray tube 111 and whose emission range is limited by the collimator 117, by allowing the X-rays to pass therethrough.

[0058] The DAS 118 acquires X-ray signals detected by each detector element included in the X-ray detector 112. For example, the DAS 118 has an amplifier that performs an amplification process on electrical signals output from each detector element and an A/D converter that converts the electrical signals to digital signals, and generates detection data. The DAS 118 is implemented by, for example, a processor.

[0059] The data generated by the DAS 118 is transmitted from a transmitter having a light emitting diode (LED) provided on the rotating frame 113 to a receiver having a photodiode provided on a non-rotating part (for example, a fixed frame and the like and not illustrated in FIG. 4) of the gantry 110 by optical communication, and is transmitted to the console 140. The non-rotating part is, for example, a fixed frame and the like that rotatably supports the rotating frame 113. Note that the data transmission method from the rotating frame 113 to the non-rotating part of the gantry 110 is not limited to the optical communication, and may adopt any non-contact type data transmission method or a contact type data transmission method.

[0060] The bed 130 is a device that places and moves the subject P to be scanned and includes a pedestal 131, a couch driving device 132, the tabletop 133, and a support frame 134. The pedestal 131 is a casing that supports the support frame 134 so as to be movable in a vertical direction. The couch driving device 132 is a driving mechanism that moves the tabletop 133, on which the subject P is placed, in a long axis direction of the tabletop 133 and includes a motor, an actuator and the like. The tabletop 133 provided on the upper surface of the support frame 134 is a plate on which the subject P is placed. Note that the couch driving device 132 may also move the support frame 134 in the long axis direction of the tabletop 133 in addition to the tabletop 133.

[0061] The console 140 has a memory 141, a display 142, the input interface 143, and processing circuitry 144. Although the console 140 is described as a separate body from the gantry 110, the gantry 110 may include the console 140 or a part of each component of the console 140.

[0062] The memory 141 is implemented by, for example, a semiconductor memory element such as a random access memory (RAM) and a flash memory, a hard disk, an optical disk, and the like. For example, the memory 141 stores a computer program for circuitry included in the X-ray CT apparatus 10 to perform its functions. Furthermore, the memory 141 stores various information obtained by imaging the subject P. Note that the memory 141 may be implemented by a server group (cloud) connected to the X-ray CT apparatus 10 via a network.

[0063] The display 142 displays various information. For example, the display 142 displays an image based on the imaging described herein. Furthermore, for example, the display 142 displays a graphical user interface (GUI) for receiving various instructions, settings, and the like from a user via the input interface 143. For example, the display 142 is a liquid crystal display or a cathode ray tube (CRT) display. The display 142 may be a desktop type display, or may be composed of a tablet terminal and the like capable of wirelessly communicating with the body of the X-ray CT apparatus 10.

[0064] Although the X-ray CT apparatus 10 is described as including the display 142 in FIG. 4, the X-ray CT apparatus 10 may include a projector instead of or in addition to the display 142. Under the control of the processing circuitry 144, the projector can perform projection onto a screen, a wall, a floor, the body surface of the subject P, and the like. As an example, the projector can also perform projection onto any plane, object, space, and the like by projection mapping.

[0065] The input interface 143 receives various input operations from a user, converts the received input operations into electrical signals, and outputs the electrical signals to the processing circuitry 144. For example, the input interface 143 is implemented by a mouse, a keyboard, a trackball, a switch, a button, a joystick, a touch pad for performing an input operation by touching an operation surface, a touch screen in which a display screen and a touch pad are integrated, non-contact input circuitry using an optical sensor, voice input circuitry, and the like. Note that the input interface 143 may be composed of a tablet terminal and the like capable of wirelessly communicating with the body of the X-ray CT apparatus 10. Furthermore, the input interface 143 may be circuitry that receives an input operation from a user by motion capture. As an example, the input interface 143 can receive a user's body movement, line of sight, and the like as an input operation by processing a signal acquired via a tracker or an image collected for a user. Furthermore, the input interface 143 is not limited to one including physical operation parts such as a mouse and a keyboard. For example, an example of the input interface 143 includes electric signal processing circuitry which receives an electric signal corresponding to an input operation from an external input device separately provided from the X-ray CT apparatus 10 and outputs the electric signal to the processing circuitry 144.

[0066] The processing circuitry 144 controls the overall operation of the X-ray CT apparatus 10 by performing a control function 144a, an imaging function 144b, an acquisition function 144c, a model generation function 144d, a noise reduction processing function 144e, and an output function 144f.

[0067] For example, the processing circuitry 144 reads a computer program corresponding to the control function 144a from the memory 141 and executes the read computer program, thereby controlling various functions, such as the imaging function 144b, the acquisition function 144c, the model generation function 144d, the noise reduction processing function 144e, and the output function 144f, based on various input operations received from a user via the input interface 143.

[0068] Furthermore, for example, the processing circuitry 144 reads a computer program corresponding to the imaging function 144b from the memory 141 and executes the read computer program, thereby imaging the subject P. For example, the imaging function 144b controls the X-ray high voltage device 114 to supply the X-ray tube 111 with a high voltage. With this, the X-ray tube 111 generates X-rays to be emitted to the subject P. Furthermore, the imaging function 144b controls the couch driving device 132 to move the subject P into an imaging port of the gantry 110. Furthermore, the imaging function 144b adjusts the position of the wedge 116 and the opening degree and position of the collimator 117, thereby controlling the distribution of the X-rays emitted to the subject P. Furthermore, the imaging function 144b controls the control device 115 to rotate the rotating part. Furthermore, while the imaging is performed by the imaging function 144b, the DAS 118 acquires X-ray signals from the respective detection elements in the X-ray detector 112 and generates detection data.

[0069] Furthermore, the imaging function 144b performs pre-processing on the detection data output from the DAS 118. For example, the imaging function 144b performs pre-processing, such as logarithmic transformation processing, offset correction processing, inter-channel sensitivity correction processing, and beam hardening correction, on the detection data output from the DAS 118. Note that the data subjected to the pre-processing is also described as raw data. Furthermore, the detection data before the pre-processing and the raw data subjected to the pre-processing are also collectively described as projection data.

[0070] Furthermore, for example, the processing circuitry 144 reads a computer program corresponding to the acquisition function 144c from the memory 141 and executes the read computer program to acquire image data. Furthermore, for example, the processing circuitry 144 reads a computer program corresponding to the model generation function 144d from the memory 141 and executes the read computer program, thereby creating and/or using the energy dependent forward model. Furthermore, for example, the processing circuitry 144 reads a computer program corresponding to the noise reduction processing function 144e from the memory 141 and executes the read computer program, thereby reducing noise in input subject data to obtain denoised data. Furthermore, for example, the processing circuitry 144 reads a computer program corresponding to the output function 144f from the memory 141 and executes the read computer program, thereby outputting an image based on the denoised data.

[0071] In the X-ray CT apparatus 10 illustrated in FIG. 4, the respective processing functions are stored in the memory 141 in the form of the computer programs executable by a computer. The processing circuitry 144 is a processor that performs a function corresponding to each computer program by reading and executing the computer program from the memory 141. In other words, the processing circuitry 144 having read the computer program has a function corresponding to the read computer program.

[0072] Note that, in FIG. 4, it has been described that the control function 144a, the imaging function 144b, the acquisition function 144c, the model generation function 144d, the noise reduction processing function 144e, and the output function 144f are implemented by the single processing circuitry 144, but the processing circuitry 144 may be configured by combining a plurality of independent processors, and each processor may be configured to perform each function by executing each computer program. Furthermore, each processing function of the processing circuitry 144 may be performed by being appropriately distributed or integrated into a single circuit or a plurality of processing circuits.

[0073] Furthermore, the processing circuitry 144 may also perform the functions by using a processor of an external device connected via the network. For example, the processing circuitry 144 reads and executes the computer program corresponding to each function from the memory 141 and uses, as computation resources, a server group (cloud) connected to the X-ray CT apparatus 10 via the network, thereby performing each function illustrated in FIG. 4.

[0074] Furthermore, although FIG. 4 illustrates only the single memory 141, the X-ray CT apparatus 10 may include a plurality of physically separated memories. For example, the X-ray CT apparatus 10 may separately include, as the memory 141, a memory that stores a computer program required when circuitry included in the X-ray CT apparatus 10 performs its function, a memory that stores various information obtained by imaging the subject P, and a memory that stores the model described herein.

[0075] Additional embodiments are set forth in the listing of parentheticals listed below. [0076] (1) A photon counting computed tomography (PCCT) apparatus including, but not limited to, processing circuitry configured to: receive a first forward model including a set of first parameters and a set of second parameters corresponding to a plurality of materials and path lengths by scanning a slab at plural tube voltages and plural tube currents of an X-ray tube, evaluate an image quality of a material decomposition image reconstructed by the set of first parameters and the set of second parameters, and update at least one second parameters from the set of second parameters if the image quality of the material decomposition image does not satisfy a predetermined threshold, wherein the update of the at least one second parameter from the set of second parameters is achieved by updating the at least one second parameter from the set of second parameters to an energy dependent parameter from a constant value. [0077] (2) The PCCT apparatus according to (1), wherein the at least one second parameter from the set of second parameters corresponds to a low tube current condition which is less influenced by pile-up than the at least one second parameter before being updated. [0078] (3) The PCCT apparatus according to either of (1) and (2), wherein the processing circuitry is further configured to update remaining second parameters from the set of second parameters corresponding to higher tube current based on the at least one second parameter from the set of second parameters corresponding to low tube current. [0079] (4) The PCCT apparatus according to any one of (1)-(3), wherein the processing circuitry is further configured to update the forward model based on the updated at least one second parameter. [0080] (5) The PCCT apparatus according to (4), wherein the processing circuitry is further configured to evaluate an image quality of the material decomposition image reconstructed based on the updated forward model. [0081] (6) The PCCT apparatus according to any one of (1)-(5), wherein the processing circuitry configured to evaluate the image quality of the material decomposition image reconstructed by the set of first parameters and the set of second parameters comprises processing circuitry configured to perform a pixel by pixel evaluation. [0082] (7) The PCCT apparatus according to any one of (1)-(6), wherein the set of second parameters are related to pileup correction. [0083] (8) The PCCT apparatus according to any one of (1)-(7), wherein the set of first parameters depend on a bin response function. [0084] (9) The PCCT apparatus according to any one of (1)-(8), wherein the processing circuitry configured to receive the first forward model comprises processing circuitry configured to scan a plurality of slabs at the plural tube voltages and the plural tube currents of an X-ray tube. [0085] (10) The PCCT apparatus according to any one of (1)-(9), wherein the energy dependent parameter updated from a constant value is updated by setting all count rate terms to zero. [0086] (11) A photon counting computed tomography (CT) method, including, but not limited to: receiving a first forward model including a set of first parameters and a set of second parameters corresponding to a plurality of materials and path lengths by scanning a slab at plural tube voltages and plural tube currents of an X-ray tube; evaluating an image quality of a material decomposition image reconstructed by the set of first parameters and the set of second parameters; and updating at least one second parameters from the set of second parameters if the image quality of the material decomposition image does not satisfy a predetermined threshold, wherein the update of the at least one second parameter from the set of second parameters is achieved by updating the at least one second parameter from the set of second parameters to an energy dependent parameter from a constant value. [0087] (12) The method according to (11), wherein the at least one second parameter from the set of second parameters corresponds to a low tube current condition which is less influenced by pile-up than the at least one second parameter before being updated. [0088] (13) The method according to either of (11) and (12), further comprising updating remaining second parameters from the set of second parameters corresponding to higher tube current based on the at least one second parameter from the set of second parameters corresponding to low tube current. [0089] (14) The method according to any one of (11)-(13), further comprising updating the forward model based on the updated at least one second parameter. [0090] (15) The method according to (14), further comprising evaluating an image quality of the material decomposition image reconstructed based on the updated forward model. [0091] (16) The method according to any one of (11)-(15), evaluating the image quality of the material decomposition image reconstructed by the set of first parameters and the set of second parameters comprises performing a pixel by pixel evaluation. [0092] (17) The method according to any one of (11)-(16), wherein the set of second parameters are related to pileup correction. [0093] (18) The method according to any one of (11)-(17), wherein the set of first parameters depend on a bin response function. [0094] (19) The method according to any one of (11)-(18), wherein receiving the first forward model includes, but is not limited to, scanning a plurality of slabs at the plural tube voltages and the plural tube currents of an X-ray tube. [0095] (20) The method according to any one of (11)-(19), wherein the energy dependent parameter updated from a constant value is updated by setting all count rate terms to zero. [0096] (21) A non-transitory computer-readable storage medium storing computer-readable instructions that, when executed by a computer, cause the computer to perform a method according to any one of (11)-(20).

[0097] Numerous modifications and variations of the embodiments presented herein are possible in light of the above teachings. It is therefore to be understood that within the scope of the claims, the disclosure may be practiced otherwise than as specifically described herein.