The present disclosure relates to a radiographing system, a radiation display method, and a storage medium. In particular, to a technique for identifying constitutive substances of a subject.
There is known a radiographing system using a flat panel detector (FPD) formed from a semiconductor material. The radiographing system can shoot radiographic images, including still images and moving images. The FPD measures the total amount of electric charges generated by incident radiation.
There has been proposed a radiographing system that calculates the average value and variance value of pixel values in each predetermined area to estimate the radiation quantum number and the average value of energy (for example, refer to Japanese Patent Laid-Open No. 2009-285356). According to Japanese Patent Laid-Open No. 2009-285356, however, with, for example, an abrupt change in the density of an image, the estimated radiation quantum number and average value of energy can include major errors.
An aspect of the present disclosure provides a radiographing system that displays statistical images based on appropriate statistical processes, a radiation display method, and a storage medium. The radiographing system includes a radiation detection apparatus that detects radiation that passed through a subject and outputs image data, a spatial statistical processing unit that acquires an average pixel value and a variance value of a predetermined pixel based on the image data and performs a spatial statistical process, a time statistical processing unit that acquires an average pixel value and a variance value in a time-series manner based on a plurality of image data and performs a time statistical process, and a display unit that displays a spatial statistical image based on spatial statistical information that is a result of the spatial statistical process or a time statistical image based on time statistical information that is a result of the spatial statistical process.
Further features will become apparent from the following description of exemplary embodiments with reference to the attached drawings.
Exemplary embodiments will be described with reference to the attached drawings.
As illustrated in
The radiographing system includes an image storage unit 16 that stores radiation images based on the radiation detected by the radiation detection apparatus 12, a spatial statistical processing unit 18 that acquires an average pixel value and a variance value in each of a plurality of pixel areas and calculates radiation quantum number and average energy, and a time statistical processing unit 20 that acquires an average pixel value and a variance value in each pixel in a time-series manner and calculates radiation quantum number and average energy. That is, the radiographing system can generate a plurality of types of statistical images.
The radiographing system includes a display control unit 22 that controls display of a spatial statistical image based on spatial statistical information processed by the spatial statistical processing unit 18 and a time statistical image based on time statistical information processed by the time statistical processing unit 20. In addition, the radiographing system includes a display unit 24 that displays the spatial statistical image based on the spatial statistical information, the time statistical image based on the time statistical information, or a radiation image based on radiation.
The radiographing system includes an operation unit 30 for an operator to perform an operation and a control unit 26 that controls the radiation generation unit 10 and the radiation detection apparatus 12.
The control unit 26 is connected to the radiation generation unit 10 via a cable. The control unit 26 sets radiographing conditions (tube voltage, tube current, and irradiation time) in the radiation generation unit 10 to control the radiation generation unit 10. The radiation generation unit 10 acts as a radiation source that generates radiation. The radiation generation unit 10 is implemented by a radiation X-ray tube, for example, and emits radiation to the subject 14, e.g., a specific part of the subject.
The radiation generation unit 10 can emit radiation to a desired irradiation area. The radiation generation unit 10 is installed via a support unit (not illustrated) on a floor surface or a ceiling. The irradiation surface of the radiation generation unit 10 includes a diaphragm (not illustrated) that shields radiation. The operator can control the diaphragm shielding radiation to set the irradiation area of the radiation emitted from the radiation generation unit 10.
The radiographing system includes a radiation detection apparatus 12 that detects the radiation emitted from the radiation generation unit 10 and that penetrated the subject 14. The radiation detection apparatus 12 detects the radiation that passed through the subject 14 and outputs image data according to the radiation. The image data can also be called radiation image.
Specifically, the radiation detection apparatus 12 detects the radiation that passed through the subject 14 as electric charges equivalent to the penetrating radiation dose. For example, the radiation detection apparatus 12 is a direct conversion-type sensor, such as a-Se, that converts radiation directly into electric charges or an indirect-type sensor that uses a scintillator, such as CsI and a photoelectric conversion element such as a-Si. The radiation detection apparatus 12 also subjects the detected electric charges to A/D conversion to generate image data, and outputs the image data to the image storage unit 16.
The operation unit 30 operates a process in the radiographing system and transfers the contents of the operation to the control unit 26 and the display control unit 22. The display unit 24 is implemented by, for example, a liquid crystal display that displays various information for an operator, e.g., radiographing technician or doctor. The operation unit 30 is, for example, a mouse, operation buttons, or the like. The display unit 24 and the operation unit 30 can be implemented as a touch panel in which the display unit 24 and the operation unit 30 are integrated.
The control unit 26 is connected to the radiation detection apparatus 12 via a cable. The control unit 26 and the radiation detection apparatus 12 exchange power source, image data, control signals, and the like with each other via the cable. The radiation detection apparatus 12 detects the radiation that passed through the subject 14 and acquires a radiation image (image data) based on the subject 14. That is, the radiation generation unit 10 and the radiation detection apparatus 12 cooperate to perform radiographing.
The time statistical process by the time statistical processing unit 20 will now be described. The time statistical processing unit 20 acquires the average pixel value and variance value of each pixel in the radiation image in a time-series manner and estimates the radiation quantum number and average energy.
The radiation detection apparatus 12 includes a scintillator that converts radiation photons into visible light photons, a photoelectric conversion element that converts the visible light photons into electric charges, and an output circuit that converts a voltage with the converted electric charges into a digital signal and outputs the same. When the scintillator absorbs the radiation photons as an example of radiation quanta, the scintillator generates radiation photons. The number of the visible light photons generated at that time varies depending on the energy of the radiation photons absorbed by the scintillator. Specifically, the greater the energy of the radiation photons, the more the visible light photons are generated in the scintillator. The amount of electric charges depending on the number of the electric charges generated by the photoelectric conversion element is determined based on the number of the visible light photons absorbed by the photoelectric conversion element. The value of the digital signal finally output from the radiation detection apparatus is determined by subjecting the value of the analog voltage with the converted amount of electric charges to the analog/digital conversion. For example, the value of the digital signal output based on the radiation photons with specific energy is 30 LSB, and the value of the digital signal output based on the radiation photons with higher energy is 100 LSB. Therefore, each time the scintillator absorbs one radiation photon, the digital signal based on the amount of electric charges generated by the photoelectric conversion element can be acquired to identify the energy of the radiation photon from the value. In this case, LSB is the quantization unit of the analog/digital conversion, and 30 LSB means 30 quantization units.
The principles of estimating the average value of energy of radiation quanta will now be described.
The quantum noise results from variations in the quantum number of the radiation quanta, e.g., the number of radiation photons, per unit time. When being regarded as occurrence probability of discrete events per unit time, the variations in the quantum number of radiation quanta are pursuant to Poisson distribution as discrete probability distribution having a specific random variable for counting discrete events occurring at predetermined time intervals. In the Poisson distribution, when the random variable taking a value of a natural number satisfies a desired condition with a constant λ>0, the random variable is pursuant to the Poisson distribution of parameter λ.
For example, when the expected value of the quantum number of the radiation quanta applied to an arbitrary pixel per unit time is 10, the quantum number of the radiation quanta actually applied to an arbitrary pixel per unit time varies among values such as 12, 5, 13, 11, . . . . In such a case, as in the foregoing example, when the value of the digital signal output according to one radiation quantum with specific energy is 30 LSB, the actual pixel value varies among values such as 360 LSB, 150 LSB, 390 LSB, 330 LSB . . . . In the case of such an example, when the sample number as the number of the acquired pixel values is increased infinitely, the expected pixel value is 300 LSB, and the value of the variation (hereinafter, called variance) is 9000 LSB.
In addition, for example, it is assumed that the expected value of the quantum number of the radiation quanta applied to an arbitrary pixel per unit time is 3 and the value of the digital signal output according to one radiation quantum with specific energy is 100 LSB. In this case, when the sample number is increased infinitely, the expected pixel value is 300 LSB and the value of variation (hereinafter, called variance) is 30000 LSB.
That is, even with the same average of the pixel value, the image formed from the radiation quanta with larger energy has wider variance of pixel values than the image formed from the radiation quanta with smaller energy. With the use of this phenomenon, it is possible to estimate the energy of radiation quanta, such as radiation photons.
A method for estimating the energy of radiation quanta using equations will now be described. First, the radiation detection apparatus is exposed to radiation T (T is a natural number of greater than or equal to 2) times to acquire image data of T images from the radiation detection apparatus. When the pixel value of one pixel in the t-th (t is a natural number greater than or equal to 2 and less than or equal to T) image data is designated as I(t), the total number of the quantum numbers of the radiation quanta having reached and been absorbed by the pixel is designated as N, and the energy of the radiation quanta is designated as E, Equation (1) can be established as follows:
E×N=ΣI(t) (1)
When the arithmetic average of the quantum numbers of the radiation quanta having reached and been absorbed by the pixel of one pixel data calculated by Equation (1) is designated as nAve, nAve can be expressed by Equation (2) as follows:
n
Ave
=N/T=Σ/(t)/E/T (2)
When the sample variance of the quantum numbers of the radiation quanta having reached and been absorbed by the pixel of one pixel data calculated by Equation (1) is designated as nVar, nVar can be expressed by Equation (3) as follows:
n
Var
=Σ[{I(t)/E−nAve}2]/T (3)
In the Poisson distribution, the estimated value and the variance are equal to the parameter λ. In addition, as the number of samples is increased, the arithmetic average comes closer to the estimated value and the sample variance comes closer to the variance. Accordingly, when the number of samples is sufficiently large (preferably, infinity) and the arithmetic average nAve of the quantum numbers of the radiation quanta and the sample variance nVar of the quantum numbers of the radiation quanta are regarded as equal and approximated, Equation (4) can be derived based on the assumption that Equation (2) and Equation (3) are equal as follows:
E=τ{I(t)2}/I(t)t)}−τ/I(t)t)} (4)
In this way, the energy E of the radiation quanta having reached and been absorbed by an arbitrary pixel in the t-th image data can be estimated from the pixel value I(t) of the pixel.
When the arithmetic average of the pixel value I(t) is designated as IAve, the arithmetic average IAve of the pixel value I(t) can be expressed by Equation (5) using the arithmetic average nAve of the quantum numbers of the radiation quanta as follows:
I
Ave
=n
Ave
×E (5)
When the sample variance of the pixel value is designated as IVar, the sample variance IVar of the pixel value can be expressed from the sample variance nVar of the quantum numbers of the radiation quanta by Equation (6) as follows:
I
Var
=n
Var
×E
2 (6)
Therefore, the energy E of the radiation quanta having reached and been absorbed by the pixel can also be expressed by Equation (7) as follows:
E=I
Var
/I
Ave (7)
In actuality, the energy of the radiation quanta having reached and been absorbed by the pixel is not single. For example, when the radiation is generated by a general radiation generation apparatus under a tube voltage of 100 kV, radiation photons with various types of energy of 100 KeV or less can be generated. For such radiation, the average value of energy of the radiation quanta having reached and been absorbed by the pixel can be estimated by approximation on the assumption that Equation (4) can hold. In addition, the number of the radiation quanta can be estimated using Equation (1) from the average value of energy of the radiation quanta and the pixel value I(t) of an arbitrary pixel in the t-th image data. In this case, the estimation process by the time statistical processing unit 20 is called time statistical process.
The spatial statistical process performed by the spatial statistical processing unit 18 will now be described. The spatial statistical processing unit 18 estimates the radiation quantum number and the average energy using the average pixel value and variance value in each predetermined area of a radiation image. In this case, the predetermined area is set on the periphery of a pixel in which the average energy is to be estimated, and the radiation quantum number and the average energy are estimated on the assumption that the radiation energy is uniform.
For example, as illustrated in
The spatial statistical processing unit 18 can arbitrarily set the pixel range 44 in the predetermined area 42 to 5×5 pixels, 3×3 pixels, or the like for the spatial statistical process. The spatial statistical processing unit 18 can estimate the radiation quantum number and the average energy for one each pixel.
The display form on the display unit 24 will now be described with reference to
As illustrated in
In addition, as illustrated in
In this way, the operator can compare the radiation image based on the spatial statistical information or the time statistical image based on the time statistical information with the radiation image including morphologic information about the subject 14. Accordingly, the operator can add information about the constitutive substances of the subject 14 acquired from the spatial statistical image or the time statistical image to a hard-to-distinguish area of the radiation image acquired by radiographing.
The sequence of statistical process and display according to the first embodiment will be described with reference to
As illustrated in
Upon completion of the spatial statistical process performed by the spatial statistical processing unit 18, the spatial statistical processing unit 18 outputs the spatial statistical image based on the spatial statistical process. The display control unit 22 controls the display unit 24 to display the spatial statistical image output from the spatial statistical processing unit 18. The display unit 24 displays the spatial statistical image.
When the spatial statistical processing unit 18 repeatedly performs the spatial statistical process, the spatial statistical processing unit 18 outputs a spatial statistical image based on each spatial statistical process. The display control unit 22 controls the display unit 24 to display the latest spatial statistical image output from the spatial statistical processing unit 18. The display unit 24 displays the latest spatial statistical image output from the spatial statistical processing unit 18.
As illustrated in
At that time, as illustrated in
The spatial statistical process performed by the spatial statistical processing unit 18 and the time statistical process performed by the time statistical processing unit 20 are controlled based on the body motion of the subject 14. Information about the body motion of the subject 14 can be acquired by observing the subject 14 with an external camera or analyzing the image data. Specifically, based on subtraction images of a plurality of image data acquired in a time-series manner, the display control unit 22 determines that there is body motion of the subject 14 with a large difference, and determines that there is no body motion of the subject 14 with a small difference. Alternatively, the display control unit 22 can determine the presence or absence of body motion by creating an average image from the plurality of image data and calculating a motion vector between the average image and each one of the plurality of image data.
With the body motion of the subject 14, the spatial statistical processing unit 18 performs the spatial statistical process in many cases. At that time, the time statistical processing unit 20 does not perform the time statistical process. It is particularly effective in the case of observing lung cancer in breathing body motion. This is because the time statistical process is not stabilized with the body motion of the subject 14. With the body motion of the subject 14, the spatial statistical processing unit 18 outputs the spatial statistical image based on the spatial statistical process. The display control unit 22 controls the display unit 24 to display the spatial statistical image output from the spatial statistical processing unit 18. The display unit 24 displays the spatial statistical image.
Without the body motion of the subject 14, the time statistical processing unit 20 performs the time statistical process. It is particularly effective in the case where long-time breath-holding with lung cancer is possible. This is because the time statistical process is stabilized without the body motion of the subject 14. Without the body motion of the subject 14, the time statistical processing unit 20 outputs the time statistical image based on the time statistical process. The time statistical processing unit 20 can perform the time statistical process excluding the image data with the body motion.
The display control unit 22 controls the display unit 24 to display the time statistical image output from the time statistical processing unit 20. The display unit 24 displays the time statistical image.
According to the embodiment, the radiographing system includes the radiation detection apparatus that detects radiation that passed through the subject 14 and outputs image data, the spatial statistical processing unit that acquires the average pixel value and variance value in a predetermined area of the image data, calculates the radiation quantum number and average energy, and performs the spatial statistical process, the time statistical processing unit that acquires the average pixel value and variance value in a plurality of image data in a time-series manner, calculates the radiation quantum number and average energy and performs the time statistical process, and the display unit that displays the spatial statistical image based on the spatial statistical information or the time statistical image based on the time statistical information.
Accordingly, the radiographing system can display statistical images based on a plurality of kinds of statistical processes. That is, the radiographing system can display a statistical image based on an appropriate statistical process in any situation.
The sequence of statistical processes and display according to a second embodiment will be described with reference to
In the mode illustrated in
As illustrated in
At that time, as illustrated in
As described above, the radiographing system can display the statistical image based on the plurality of kinds of statistical processes. The statistical image based on the spatial statistical process and the time statistical process is higher in image resolution than the spatial statistical image based on the spatial statistical process. This enables the operator to identify the constitutive substances of the subject 14 in a quick and correct manner.
The sequence of statistical processes and display according to a third embodiment will be described with reference to
As illustrated in
Upon completion of the spatial statistical process by the spatial statistical processing unit 18, the spatial statistical processing unit 18 outputs the spatial statistical image based on the spatial statistical process. The display control unit 22 controls the display unit 24 to display the spatial statistical image output from the spatial statistical processing unit 18. The display unit 24 displays the spatial statistical image.
As illustrated in
At that time, as illustrated in
When the time statistical processing unit 20 repeatedly performs the time statistical process, the time statistical processing unit 20 outputs the time statistical image based on the time statistical process each time. Specifically, the time statistical processing unit 20 performs the time statistical process using the first output image data in a time range A, and outputs the time statistical image based on that time statistical process. Then, the time statistical processing unit 20 performs the time statistical process using the next output image data in a time range B, and outputs the time statistical image based on the time statistical process. The time range A and the time range B can overlap each other.
The display control unit 22 controls the display unit 24 to display the latest time statistical image output from the time statistical processing unit 20. The display unit 24 displays the latest time statistical image output from the time statistical processing unit 20.
Accordingly, the radiographing system can display the latest time statistical image. The operator can identify quickly the constitutive substances of the subject 14.
The sequence of statistical processes and display according to a fourth embodiment will be described with reference to
As illustrated in
Next, as illustrated in
Then, as illustrated in
Accordingly, the radiographing system can display the colored time statistical image and the black-and-white radiation image at the same time. The operator can observe the morphologic information and constitutive substances of the subject 14.
A fifth embodiment will be described with respect to
Although not illustrated, when parallel is selected as display mode and the spatial statistical image and the radiation image are selected as display targets, the display unit 24 displays the spatial statistical image and the radiation image in parallel as illustrated in
The operator presses an execute button 62 to enable the selection. The operator presses a cancel button to disable the setting. The number of image data necessary for the time statistical process can be set within a range of 10 to 200, for example. The repetition cycle of the time statistical process can be set within a range of 2 to 100, for example. In this case, the number of statistics is set to 100 and the repetition cycle is set to 5. That is, the display unit 24 displays the time statistical image having undergone the time statistical process using 100 image data. The display unit 24 updates and displays the time statistical image five times. Accordingly, the type of the image to be displayed can be set on the display unit 24, the number of image data to undergo the statistical process, and the repetition cycle of the statistical process.
When either the pixel range for the spatial statistical process or the number of statistics for the time statistical process is set, the display control unit can automatically set the other such that total predetermined pixels (for example, 100 pixels) can be achieved. Specifically, when the display control unit 22 sets the pixel range for the spatial statistical process to 10×10, the number of statistics for the time statistical process is set to 1. When the display control unit 22 sets the pixel range for the spatial statistical process to 5×5, the number of statistics for the time statistical process is set to 4. When the display control unit 22 sets the pixel range for the spatial statistical process to 2×2, the number of statistics for the time statistical process is set to 25.
Accordingly, the radiographing system can control the display mode of the image to be displayed on the display unit 24. The operator can observe the desired image.
Embodiments can also be realized by a computer of a system or apparatus that reads out and executes computer executable instructions recorded on a storage medium (e.g., non-transitory computer-readable storage medium) to perform the functions of one or more of the above-described embodiment(s) of the present invention, and by a method performed by the computer of the system or apparatus by, for example, reading out and executing the computer executable instructions from the storage medium to perform the functions of one or more of the above-described embodiment(s). The computer may comprise one or more of a central processing unit (CPU), micro processing unit (MPU), or other circuitry, and may include a network of separate computers or separate computer processors. The computer executable instructions may be provided to the computer, for example, from a network or the storage medium. The storage medium may include, for example, one or more of a hard disk, a random-access memory (RAM), a read only memory (ROM), a storage of distributed computing systems, an optical disk (such as a compact disc (CD), digital versatile disc (DVD), or Blu-ray Disc (BD)™), a flash memory device, a memory card, and the like.
While exemplary embodiments have been described, it is to be understood that the invention is not limited to the disclosed exemplary embodiments. The scope of the following claims is to be accorded the broadest interpretation so as to encompass all such modifications and equivalent structures and functions.
This application claims the benefit of Japanese Patent Application No. 2016-212127, filed Oct. 28, 2016, which is hereby incorporated by reference herein in its entirety.
Number | Date | Country | Kind |
---|---|---|---|
2016-212127 | Oct 2016 | JP | national |