This disclosure relates to an X-ray CT apparatus and, in particular, relates to an improvement in an X-ray CT apparatus that has a function of aiding setting of many parameters to be set at a scanning planning stage prior to a scanning when performing an existence diagnosis of a lesioned portion, for example, a cancer in an organ of an object to be examined.
Since a contrast, amount of noises and a discrimination capacity of an image of a lesioned portion representing a target tissue with respect to a background tissue displayed on a display of an X-ray CT apparatus vary complexly in association with such as scanning conditions to be set (such as tube voltage, tube current, circumferential rotation speed, helical pitch and slice thickness), reconstruction conditions (such as reconstruction filters) and size of the target tissue, it was difficult to estimate a quality of images of the target tissue to be produced and to determine optimum parameters to be set at a scanning planning stage prior to a scanning.
JP-A-2004-97778 discloses an X-ray computer tomographic scanning apparatus comprising a plan assisting system for resolving such problem that includes a condition optimizing function as well as a scan planning screen building function having a graphical user interface (GUI). Wherein, the apparatus is constituted so as to assist a scanning planning in such a manner that when a dose represented by CTDI (CT Dose Index) is inputted on the scan planning screen having GUI, a DEI (Dose Efficiency Index) showing a target diameter having 50% detection rate (discrimination rate) on a reconstructed image and its phantom sample image corresponding thereto are displayed on the scan planning screen, and further, when a desired detection rate, target diameter or density difference between a target tissue defined by a target diameter and its background tissue are inputted, a target diameter, detection rate or DEI corresponding thereto is displayed on the scan planning screen.
However, when the energy (a tube voltage) for X-ray irradiation varies, the CT value of the target tissue for scanning and the CT value of the background tissue around the target tissue vary and the difference of the CT values, namely the contrast varies at the same time. Since the DEI shows a characteristic which varies in response to the contrast, when taking into account of such influences due to contrast, it is necessary to prepare a data base provided with many DEI characteristics for every contrasts. Since the data base provided with such DEI characteristics has to be prepared for every apparatus by making use of a phantom, many man hours are required therefor and at the same time since the DEI characteristics in the data bases are significantly affected by characters of individual evaluators on the target diameters, the level of user friendliness of the apparatus differs depending on experiences of users using the apparatus thereafter.
In an aspect of this disclosure, there is provided an X-ray CT apparatus having a versatile function in which only by setting limited numbers, of simple parameters to be set at a scanning planning stage prior to a scanning, the setting of other parameters to be set is aided.
More specifically, in an X-ray CT apparatus provided with an X-ray generating means, an X-ray detector disposed opposing to the X-ray generating means, a rotating means for rotating the X-ray generating means and the X-ray detector on a circular orbit of a same rotation center, a control means for controlling so as to irradiate X-rays from the X-ray generating means to an object laid along the rotation center and to detect by the X-ray detector X-ray dose penetrated through the object while rotating the X-ray generating means and the X-ray detector through the rotating means and a reconstruction computing means for performing reconstruction computation by making use of the data of the penetrated X-ray dose obtained under the control of the control means and acquiring a tomographic image,
the X-ray CT apparatus according to another aspect of this disclosure is further provided with an input means being inputted of information with respect to an target tissue of the object to be discriminated in the tomographic image and of an index with respect to correctness of the discrimination, and a scanning condition determining means for determining scanning conditions for discriminating the target tissue based on the index inputted by the input means. Further, a scanning condition determining method for an X-ray CT apparatus,
the scanning condition determining method for an X-ray CT apparatus according to another aspect of this disclosure is provided with,
(1) the step of inputting information with respect to an target tissue to be discriminated in a tomographic image to be acquired by the X-ray CT apparatus and an index with respect to correctness of the discrimination, and
(2) the step of determining scanning conditions for discriminating the target tissue with the index based on the information with respect to the target tissue, the index and apparatus characteristics of the X-ray CT apparatus.
According to another aspect of this disclosure, an X-ray CT apparatus having a versatile function is provided in which by setting limited numbers of simple parameters to be set at a scanning planning stage prior to a scanning, the setting of other parameters to be set is aided.
According to another aspect of this disclosure, by inputting simple input items through an input device at a scanning planning stage prior to a scanning, even a user with less experience of the X-ray CT apparatus can easily acquire an image having of a desired quality by the user with regard to a visibility (recognition capability).
a) and 3(b) are diagrams for explaining a principle applied in the present invention.
a), 4(b) and 4(c) are also diagrams for explaining a principle applied in the present invention.
Herein below, a preferred embodiment of an X-ray CT apparatus according to the present invention will be explained with reference to the accompanied drawings.
The scanner 10 is constituted by such as an X-ray generation device 11, a high voltage switching unit 18a, a high voltage generation device 18b, an X-ray control device 18c, the X-ray detector 12, a preamplifier 17, a collimator 13, a collimator control device 13a, a drive device 14 that rotatably drives the X-ray generation device 11 and the X-ray detector 12, a scanner control device 15 that controls the drive device 14 and the collimator control device 13a, a bed control device 21 that performs elevation of the bed 20 and movement control in a body axial direction, a bed movement measurement device 22 that measures the amount of movement of the bed 20 and a central control device 16 for controlling the scanner 10 and the operation unit 30.
The operation unit 30 is provided with an input and output device 31 and the computing device 32. The input and output device 31 is provided with the input device 311 constituted of a mouse and a keyboard, the display device 312 constituted of a monitor displaying such as reconstructed images and the memory device 313 constituted of a display memory, a main memory and a hard disk such as for storing display data such as the reconstructed images to be displayed on the display device 312 and for storing a reconstruction program. The computing device 32 is provided with an image processing means 321 that performs image processings such as a preprocessing, a filtering processing and a post processing based on projection data sent from the scanner 10 and a reconstruction computing device 322 for reconstructing the above projection data.
From the input device 311 in the operation unit 30 are inputted such as device conditions (such as a target angle, a target material and a target density of the X-ray tube, an aluminum equivalent of characteristic X-ray filtering of the X-ray tube, a scintillator material, a scintillator density, a scintillator thickness, a shape of bow tie filter, material of the bow tie filter and a shape of copper filter), information on inspection portion of scanning object (target tissue of inspection portion, size thereof and recognition rate thereof), scanning conditions (such as tube current, tube voltage, circumferential rotating speed and helical pitch), reconstruction conditions (image FOV, reconstruction filter, slice thickness of image and slice position to be reconstructed) and reference information (such as CT values or difference of the CT values of an target tissue and a background tissue of every inspection portions at a standard scanning condition or densities of an target tissue and a background tissue and X-ray attenuation characteristics of respective target tissue and background tissue, and based on the indications control signals necessary for the scanning are sent from the central control device 16 to the X-ray control device 18c, the bed moving device 21 and the scanner control device 15 and after receiving a scanning start signal a scanning operation begins. Further, such as the device conditions and the reference information are inputted in advance via the input device 311 and are stored in the memory device 313. When a scanning operation begins, a control signal is sent from the X-ray control device 18c to the high voltage generation device 18b, a high voltage is applied to the X-ray generation device 11 via the high voltage switching unit 18a and X-rays are irradiated from the X-ray generation device 11 to the object 40. At the same time, a control signal is sent from the scanner control device 15 to the drive device 14 and the X-ray generation device 11, the X-ray detector 12 and the preamplifier 17 are rotated around the object 40.
On the other hand, the bed 20 laying the object 40 is caused by the bed control device 21 to be standing still (at the time of circular scanning) or to move (at the time of helical scanning) in the circumferential rotation axis of such as the X-ray generation device 11. The X-rays irradiated are restricted by the collimator 13 of its irradiation region, are absorbed (attenuated, weakened), penetrate through the object 40 and are detected by the X-ray detector 12. The X-rays detected by the X-ray detector 12 are converted into current, are amplified by the preamplifier 17 and are inputted into the computing device 32 as a projection data signal. The projection data signal inputted to the computing device 32 is subjected to image reconstruction processing at the reconstruction computing device 322 in the computing device 32. The reconstructed image is stored in the memory device 313 within the input and output device 31 and is displayed on the display device 312 as a CT image.
When plotting distributions of pixel values (CT values) of an target tissue of an inspection portion, for example, a contrast study processed cancer in a liver and of a background tissue surrounding the same that are displayed on the display device in the X-ray CT apparatus, the plotting shows normal distributions as shown in
As seen from
D=f(SA,SB,C,σA,σB)
a),
a) is a discrimination capability characteristic with respect to contrast. Near a low contrast, the recognition rate is low due to the effect of noises, but when the contrast rises comparatively high, the recognition rate suddenly enhances, and when the contrast further rises, the recognition rate becomes constant (100%).
b) is a discrimination capability characteristic with respect to noises. When the noises are extremely low with respect to the contrast, the recognition rate is high even if the noises increase some, but when the noises increase to a certain level, the recognition rate suddenly lowers, and when the noises further increase, the recognition rate become constant (0%).
c) is a discrimination capability characteristic with respect to size. When the size is small, the recognition rate is low due to the effect of noises, but when the size increases to a certain degree, the recognition rate suddenly enhances, and when the size further increases, the recognition rate becomes constant (100%). When the contrast increases, the recognition rate is high even if the size is small, and when the noises increase, the recognition rate lowers when the size is small.
In the present invention, an X-ray CT apparatus having a function of aiding setting of parameters to be set by a user at a scanning planning stage prior to a scanning is constructed by making use of the relationship between the recognition rate D of the target tissue A, the sizes SA, SB of the target tissue A and the background tissue B, the contrast C between the target tissue A and the background tissue B and the standard deviation σA, σB corresponding to the amount of noises of the target tissue A and the background tissue B.
In the present embodiment, a sequence of scanning planning prior to a scanning will be explained herein below along with
In step S1, a user at first inputs a cancer exceeding 2 cm possibly formed in the liver as an target tissue A in the scanning target inspection portion from the input device 311, and further inputs a recognition rate, for example, 80% representing an index value of a desired image quality. When these are inputted, a corresponding background tissue B and the size of the background tissue B corresponding to the inputted size of the target tissue A are read out from the data base of the memory device 313.
In the data base of the memory device 313, sizes of background tissue corresponding to the sizes of target tissue for every inspection portions, X-ray attenuation characteristics and densities of both tissues and CT values of target tissue and background tissue for every inspection portions under the standard scanning condition are stored.
Further, in the data base of the memory device 313, are stored as shown in
In step S2, effective spectrum of irradiation X-rays and effective X-ray energy for every tube voltages irradiated from the X-ray generation device 11 to the target tissue A and the background tissue B set in step S1 and detected by the detector 12 are calculated.
Further, the calculation of the effective spectrum of irradiation X-rays for every tube voltages and respective X-ray energies can be performed in advance and the result thereof can be stored in advance in the data base of the memory device 313.
IE=Nρ/A∫ET0(l+T/m0c2)Q(dT/dl)−1exp(−μ(E)lt cot α)dT (1)
Ich∝(T0/TK,L)1.63 (2)
Wherein, IE is the energy intensity of bremsstahlung X-rays having photon energy E (=hμ), Ich is the energy intensity of the characteristic X-rays, N is Avogadro number, ρ, A are respectively density and atomic weight of the target, m0 is weight of an electron, c is velocity of light, T0 is an energy of an incident electron, Q is an energy intensity of X-rays emitted from one electron and is a value determined approximately by the ratio of photon energy E/electron energy T, dT/dl a theoretical formula with regard to stopping power proposed such as by Bethe, et al. lt is a penetration depth of electron and TK,L is an energy necessary for removing an electron from K, L electron orbits. Further, for the calculation of X-ray spectrum, already known other methods can be used.
As shown in
From the calculated effective spectra of the irradiated X-rays, ratios of every energies with respect to all of the spectra as contribution rates, and effective energies are calculated based on the calculated contribution rates.
Graphs 150 in
Graphs 155 and 156 in
Further, as an alternative, the apparatus can be constituted in such a manner that while causing to display the X-ray attenuation characteristic μA Of graph 155 and the respective effective energies 151, 152 and 153 and the X-ray attenuation characteristic μB of graph 156 and the respective effective energies 151, 152 and 153 on the display device 312, and when a user clicks a desired effective energy, a tube voltage corresponding thereto is selected. Thereby, after the user recognizes the difference of the X-ray attenuation coefficients, the user can select an effective energy under a desired condition such as putting priority on low exposure rather than recognition capability.
Step S4 is a step that is executed in a case where densities DA, DB of the target tissue A and the background tissue B in the inspection portion are not stored in the database within the memory 313, for example, as in the present embodiment where the target tissue is contrast imaged by making use of a contrast medium. When the densities DA, DB of the target tissue A and the background tissue B in the inspection portion are stored in the database within the memory 313, this step is omitted.
In step S4, densities DA, DB of the target tissue A of contrast study and the background tissue B corresponding thereto are calculated by making use of CT values PVA, PVB of the target tissue A of contrast scanning and the background tissue B corresponding thereto under a standard scanning condition (for example, tube voltage of 120 kv) which are stored in the database within the memory device 313, their attenuation coefficients μA, μB and an attenuation coefficient of water μW under the same condition and according to the following formulas 3-1 and 3-2;
PVA=(DAμA−μW)·1000/μW (3-1)
PVB=(DBμB−μW)·1000/μW (3-2)
In step S5, by making use of the densities DA, DB of the target tissue A of contrast scanning and the background tissue B calculated at step S4, the attenuation coefficients μA, μB corresponding to the tube voltage selected at step S3 and the attenuation coefficient of water μW under the same condition and according to the following formula a contrast C (difference of the CT values) between the target tissue A of contrast scanning and the background tissue B in the inspection object is calculated;
C=(DAμA−DBμB)·1000/μW (4)
In step S6, by making use of the size of the target tissue A set at step S1, its recognition rate, the size of the background tissue B read from the database in the memory device 313 based on the size of the target tissue A and the contrast C between the target tissue A and the background tissue B calculated at step S5 according to the tube voltage selected at step S3, and according to a formula relating to the recognition rate D of the target tissue A, the sizes SA, SB of the target tissue A and the background tissue B converted into pixel numbers thereof, the contrast C between the target tissue A and the background tissue B and amount of noises σA, σB of the target tissue A and the background tissue B or the approximation curves defining these relationships as exemplified in
In step S7, based on the amount of noises calculated at step S6, a product of tube current and time mAs defining X-ray irradiation dose irradiated from the X-ray generation device 11. The calculation sequence of mAs will be explained with reference of
A scanogram image 101 is acquired through a scanogram scanning 100. For every one line of 101a and 101b perpendicular to the body axis of the scanogram image 101, the maximum value b of projection value 102 and the integrated value RS of the projection value are determined, and an elliptical model approximating a human body 103 having major axis A and minor axis B (an elliptical body of a homogenous material having attenuation coefficient near a human body) is estimated. The attenuation coefficient μP of the elliptical model approximating human body is determined by taking into account of the effective spectra corresponding to the selected tube voltage. The irradiation dose mAs (tube current×circumferential rotating speed) is determined based on the lengths (A, B) of major axis and the minor axis of the elliptical model approximating human body 103, a helical pitch, a view weight for reconstruction, a reconstruction filter and the amount of noises calculated at step S7. More specifically, at first, a penetration distance lP (θ) of X-ray beams irradiated from circumferential rotating phase θ and passing through the center of the elliptical model approximating a human body 103 is calculated based on the following formula (5);
lP(θ)=A×B/sqrt(A2 cos2θ+B2 sin2θ) (5)
Subsequently, based on the calculated penetration distance lP (θ), a function (104) of modulating mAs (θ) in circumferential rotating direction is produced by modifying the equation shown in
σ2=σbase2×CST×CmAs×CL×CRW×CRF (6)
Wherein σbase2, CST, CmAs, CL, CRW and CRF are terms for taking into accounts of the influences when the slice thickness, irradiation dose (tube current×circumferential rotating speed), penetration distance view weight for reconstruction and reconstruction filter are respectively changed and are determined according to the following formulas (7);
CST=STbase/ST
CmAs=Σ{mAsbase/mAs(θ)}
CL=Σexp(μP(lP(θ)−lbase))
CRW=Σ{W(θ)2}/{ΣW(θ)}2
CRF=σRF2/σRFbase2 (7)
Wherein lbase, mAsbase and STbase are respectively a reference diameter of circular model approximating human body, a reference irradiation dose (tube current×circumferential rotating speed) and a reference slice thickness, W(θ) is a view weight for reconstruction, θ is a view phase (circumferential rotating phase), σbase2 is noises (variance value) acquired from an image obtained by reconstructing around 360 degrees projection data imaged under a condition of lbase, mAsbase and STbase without applying view weight, and σRF2, σRFbase2 are respectively noises of the reconstruction filter and noises of a reference reconstruction filter. Further, when no view weight for reconstruction is applied, CRW is 1, and when view weight for reconstruction is applied, the value increases. Although not specifically explained herein, the view weight for reconstruction relates to the helical pitch, and in case when the helical pitch is large, since the amount of data usable for the reconstruction decreases, CRW enlarges. Further, when defining a generally used reconstruction filter for stomach as a reference, CRW of the same is 1 and that of a filter of which high frequency components are intensified is larger than 1.
When scanning a plurality of portions once through a helical scanning and required to vary noises according to the portions, if a plurality of lines (for every portions) for defining noises are set and σbase2, lbase are varied at the lines set, the predetermined irradiation dose (tube current×circumferential rotating speed) can be achieved while maintaining the target recognition rate.
Further, with regard to the irradiation dose control based on the amount of noises, the present invention is not limited to the above method but can use any other methods such as conventional methods.
In step S8, from a plurality of standard projection data for every respective inspection portions of a human body imaged under a standard condition and stored in the database within the memory device 313, a standard data of a body type being close to the projection value of the scanogram of the object imaged previously is selected, the size of the selected projection data is corrected based on the integration value of the projection value of the scanogram, the pseudo projection data of the target tissue having the size inputted at step S1 is superposed on the size corrected standard projection data, the respective amounts of noises calculated previously at step S6 and converted into projection data are added on the standard projection data being superposed of the pseudo projection data and the standard projection data to which the noises converted into projection data are added and on which the pseudo projection data are superposed are displayed on the display device 312 as a presumed reconstruction image of the target tissue in recognition rate of 80% for serving the evaluation by an operator (herein, when converting the amount of noises calculated at step S6 into those in the projection data, it is sufficient if a relationship between projection data obtained by scanning such as water phantoms having various sizes in advance and values of image noises at those moments is prepared as a function and the amount of noises to be added to the projection data is determined based on the prepared function).
The points as has been explained hitherto will be explained specifically with reference to
The user evaluates the presumed reconstruction image 123 displayed on the screen 120, and when the user satisfies the image, the user pushes the execution button to move to the actual scanning operation according to the set condition. When the user dissatisfies the image, the user modifies such as the recognition rate, the size of the target tissue, the reconstruction condition and the contrast scanning condition that were inputted and set previously, and repeats steps S1 through S8.
As in the above manner, since the operator can recognize before the actual scanning the visibility of the target tissue from the presumed image presumed from the set scanning condition, erroneous settings of the scanning condition can be decreased, possible exposure of the object due to scanning once again is avoided and further diagnosis accuracy can be enhanced.
Further, different from a method that simply adds noises on a reconstructed image, in the present embodiment, since noises corresponding to the data value of the projection data are added, an image having quality further close to the actual one (with regard to such as noises and artifacts) can be obtained.
Further, herein, the data actually imaged of a human body are used as standard data, an virtual human body model produced from computer simulation and phantom data used for evaluating image quality and spatial resolution can be used.
As other usages of the present invention, the present invention can be used to determine an optimum scanner constitution by evaluating the effects of such as compensating filters and copper filters of different shapes and materials by means of the presumed image or the recognition rate. In this instance, the determination can be realized by preparing a plurality of compensating filters, copper filters and lead filters for the scanner, and exchanging the filter depending on the setting. Further, the present invention can be used to determine an optimum reconstruction condition by changing the shape of reconstruction filter depending on the set parameter and evaluating the effect of the reconstruction filter by means of the presumed image or the recognition rate. As will be understood from the above, the present invention can be used as a simulator for optimizing the scanning condition, the device condition, the contrast scanning condition and the reconstruction condition.
Further, the present invention is not limited to the above embodiments and can be carried out by modifying in a variety of manners within a range not exceeding the gist of the present invention. For example, in the above embodiment, as an index how accurately the target tissue is discriminated, the ratio of the portion of the normal distribution of the target tissue not overlapped with the normal distribution of the background tissue to the normal distribution of the target tissue is introduced which is defined as the discrimination capability, however, the present invention is not limited to such index as defined. For example, when data are obtained and stored in advance in what relationship between the normal distributions representing pixel distributions of the target tissue and the background tissue on the histograms as shown in
Number | Date | Country | Kind |
---|---|---|---|
2005-378521 | Dec 2005 | JP | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/JP2006/325777 | 12/25/2006 | WO | 00 | 5/27/2008 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2007/074772 | 7/5/2007 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
6778628 | Yamazaki et al. | Aug 2004 | B2 |
Number | Date | Country |
---|---|---|
2004-97778 | Apr 2004 | JP |
Number | Date | Country | |
---|---|---|---|
20090245459 A1 | Oct 2009 | US |