Multi-energy CT imaging system and imaging method

1. A multi-energy CT imaging system comprising:

a stage for carrying an object to be inspected;

a voltage-regulatable X-ray generation device for emitting X-rays for radiating said object to be inspected;

a detector for receiving X-rays emitted from said X-ray generation device and penetrating said object to be inspected, and for outputting detection data;

a rack having said X-ray generation device and said detector mounted thereon; and

a data processing and control device for controlling said stage, said X-ray generation device, said detector and said rack, processing detected data from said detector, and during one rotation of scan of said X-ray generation device, evenly dividing a circular track of said X-ray generation device into angle intervals of a predetermined number according to a predetermined number of energies and setting a different high voltage of said X-ray generation device for each angle interval, and when said rack rotates from the current angle interval into a next angle interval, controlling said X-ray generation device according to an instruction from said data processing and control device to switch said X-ray generation device to a voltage set in the next angle interval.

2. A method for imaging with a multi-energy CT imaging system, comprising:

(a) providing a multi-energy CT imaging system of claim 1 ;

(b) evenly dividing the circular track of the X-ray generation device of said multi-energy CT imaging system into angle intervals of a predetermined number according to the predetermined number of energy and setting a different high voltage of said X-ray generation device for each angle interval;

(c) rotating said X-ray generation device along the circular track and working in each angle interval at a high voltage set for said angle interval to emit X-rays to irradiate the object to be inspected;

(d) acquiring projection data of different energies using the detector of said multi-energy CT imaging system;

(e) obtaining line integral data of the projection data of said different energies using the data processing and control device of said multi-energy CT imaging system;

(f) splicing said line integral data of different energies one by one to obtain a complete line integral data covering 360°;

(g) reconstructing according to said complete line integral data to obtain a hybrid reconstruction result;

(h) using said hybrid reconstruction result as a priori image to obtain attenuation coefficient reconstruction results in different energies;

(i) obtaining decomposition coefficients from said attenuation coefficient reconstruction results and obtaining an atomic number distribution image and electron density distribution image using said decomposition coefficient.

3. The method of claim 2 , wherein the circular track of step (b) is divided into N angle intervals and the line integral data in different energies obtained in step (e) are made to be P 1 , P 2 , . . . P n , and in step (g), reconstruction is performed using a common circular track filtering back projection method so as to make the obtained hybrid reconstruction result to be f 0 , N and n being integers greater than or equal to 2;

in step (h), an attenuation coefficient reconstruction result fi that satisfies the condition of ∥H i f t −p i ∥ 2 <ε is obtained by means of arg min G(f i , f 0 ), wherein, G (f i f 0 ) is a function describing the difference between the priori image and the reconstruction result, H i , is a projection operator in the i th energy, εis a parameter related to measurement of data noise variance, i=1, 2, . . . , n.

4. The method of claim 3 , characterized in that in step (h), a decomposition coefficient is obtained using a post-processing base function decomposition method, then said decomposition coefficient is used to obtain the atomic number distribution image and electron density distribution image.

5. The method of claim 3 , characterized in that in step (i), a decomposition coefficient is obtained using a post-processing base material decomposition method, then said decomposition coefficient is used to obtain the atomic number distribution image and electron density distribution image.

6. The method of claim 3 , characterized in that in place of said step (i), there is such a step (j) as performing orthographic projection to the attenuation coefficient reconstruction result f i , obtained in step (h) to obtain the completed projection data in which the scanning plane in different angles have been missing, then obtaining coefficients of different base functions by means of pre-processing, thus obtaining the atomic number distribution image and electron density distribution image.