Estimation device, method and program

Abstract

To accurately predict a locomotorium disease in an estimation device, method, and program.SOLUTION: An estimation device includes at least one processor. The processor functions as a learned neural network for deriving an estimation result of at least one emphasis image that emphasizes a specific composition of a subject from a plain two-dimensional image acquired by executing plain radiography of the subject including a plurality of compositions. The learned neural network is learned by using, as teacher data, two radiation images acquired by imaging the subject by radioactive rays with different energy distribution, and an emphasis image for learning that emphasizes the specific composition of the subject derived from the two radiation images.SELECTED DRAWING: Figure 3

Description

[Technical field]

[0001] The present disclosure relates to an estimation device, a method, and a program. [Background technology]

[0002] Conventionally, energy subtraction processing has been known, which uses two radiation images obtained by irradiating a subject with two types of radiation with different energy distributions, taking advantage of the fact that the attenuation of transmitted radiation differs depending on the material that constitutes the subject. The energy subtraction processing is a method of obtaining an image in which a specific structure is emphasized by matching each pixel of the two radiation images obtained as described above, multiplying the pixels by an appropriate weighting coefficient, and then performing subtraction. In addition to bones and soft tissues, the composition of the human body, such as fat and muscle in the soft tissues, is also derived by energy subtraction processing (see Patent Document 1).

[0003] In addition, various methods have been proposed for deriving a radiological image different from the radiological image acquired by photographing a subject, using the radiological image acquired by photographing the subject. For example, Patent Document 2 proposes a method for deriving a bone image from a radiological image of a subject acquired by plain radiography, by using a trained model constructed by training a neural network using a radiological image of the subject acquired by plain radiography and a bone image of the same subject as training data.

[0004] Note that plain radiography is a method of radiography in which a subject is irradiated with radiation once to obtain a single two-dimensional image that is a transmission image of the subject. In the following description, a two-dimensional image obtained by plain radiography will be referred to as a plain two-dimensional image. [Prior art documents] [Patent documents]

[0005] [Patent Document 1] JP 2018-153605 A [Patent Document 2] U.S. Patent No. 7,545,965 Summary of the Invention [Problem to be solved by the invention]

[0006] However, it is desirable to estimate images that emphasize specific components of bones and the like with higher accuracy.

[0007] The present disclosure has been made in consideration of the above circumstances, and has an object to make it possible to estimate with high accuracy an image in which a specific composition is emphasized. [Means for solving the problem]

[0008] The present disclosure provides an estimation apparatus comprising at least one processor, The processor The present invention functions as a trained neural network that derives an estimation result of at least one enhanced image in which a specific composition of a subject is enhanced from a simple two-dimensional image obtained by simply photographing a subject including a plurality of compositions, The trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training enhanced image derived from the two radiation images, which enhances a specific composition of the subject.

[0009] In the estimation device according to the present disclosure, the learning weighted image may be derived by energy subtraction processing in which two radiation images are weighted and subtracted from each other.

[0010] In the estimation device according to the present disclosure, the emphasized image is at least one of a bone image in which a bone part of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized; The training images are Recognizing bones and soft tissues of the subject using at least one of the two radiation images; deriving attenuation coefficients for the bones and the soft tissues using the bone and soft tissue recognition results and the two radiological images; It may be derived by performing an energy subtraction process using the attenuation coefficient.

[0011] In the estimation device according to the present disclosure, the emphasized image is a bone image in which the bones of the subject are emphasized and a soft part image in which the soft parts of the subject are emphasized, The training images are deriving new weighting coefficients to be used in the weighted subtraction based on pixel values ​​of bone tissue included in the bone image and pixel values ​​of soft tissue included in the soft tissue image; A new bone image and a new soft tissue image are derived by performing weighted subtraction on the two radiographic images using new weighting coefficients; The weighting coefficients may be derived by repeatedly deriving a new weighting coefficient based on a new bone image, a new weighting coefficient based on a new soft tissue image, and a new bone image and a new soft tissue image based on the new weighting coefficients.

[0012] In the estimation device according to the present disclosure, the emphasized image is a bone image in which the bones of the subject are emphasized and a soft part image in which the soft parts of the subject are emphasized, The training images are While varying the attenuation coefficient for each different energy distribution for soft tissue, the thickness of the soft tissue, the attenuation coefficient for each different energy distribution for bone, and the thickness of the bone from their initial values, a difference between the value of (soft tissue attenuation coefficient × soft tissue thickness + bone attenuation coefficient × bone thickness) and each pixel value of the radiographic image is derived for each different energy distribution; Derive soft tissue attenuation coefficients and bone attenuation coefficients for different energy distributions that have a minimum difference or that have a difference less than a predetermined threshold; The weighting coefficient may be derived by performing an energy subtraction process using a weighting coefficient derived based on the attenuation coefficient of soft tissue and the attenuation coefficient of bone.

[0013] In the estimation device according to the present disclosure, the learning enhanced image is Deriving composition ratios of multiple components contained in the soft part of the subject; Derive the attenuation coefficient for each energy distribution that differs for each pixel of the two radiographic images according to the composition ratio; The weighting coefficient may be derived by performing an energy subtraction process using a weighting coefficient derived based on the derived attenuation coefficient.

[0014] In addition, in the estimation device according to the present disclosure, the composition ratio is deriving a body thickness of the subject as a first body thickness and a second body thickness for each pixel in each of the two radiographic images; It may be derived for each pixel of the radiographic image based on the first body thickness and the second body thickness.

[0015] In addition, in the estimation device according to the present disclosure, the composition ratio is Deriving a first body thickness and a second body thickness based on the attenuation coefficient for each of the different energy distributions for each of the plurality of compositions; The first body thickness and the second body thickness may be derived while changing the thickness of the composition and the attenuation coefficient for each composition, and may be derived based on the thickness of the composition at which the difference between the first body thickness and the second body thickness is equal to or less than a predetermined threshold value.

[0016] In the estimation device according to the present disclosure, the learning enhanced image is The energy subtraction process may be performed to remove from the two radiation images the scattered radiation components scattered by the subject out of the radiation irradiated to the subject, and the energy subtraction process may be performed on the two radiation images from which the scattered radiation components have been removed.

[0017] In the estimation device according to the present disclosure, the scattered radiation removal process includes acquiring radiation characteristics according to a body thickness distribution of an object interposed between the subject and a radiation detector that detects a radiation image, Deriving a primary ray distribution and a scattered ray distribution of radiation contained in each of the two radiographic images using the imaging conditions, the body thickness distribution, and the radiation characteristics of the object; deriving an error between the sum of the primary ray distribution and the scattered ray distribution for each of the two radiographic images and the pixel value at each position of the two radiographic images, updating the body thickness distribution so that the error is less than a predetermined threshold value, and repeatedly deriving radiation characteristics based on the updated body thickness distribution and deriving the primary ray distribution and the scattered ray distribution contained in each of the two radiographic images; This may be performed by subtracting the scattered radiation distribution when the error falls below a predetermined threshold value from each of the two radiation images.

[0018] In the estimation device according to the present disclosure, the scattered radiation removal process includes deriving a first primary ray distribution and a scattered ray distribution of radiation that has passed through the subject using two radiation images; deriving a second primary ray distribution and a scattered ray distribution of the radiation that has passed through the object, using the first primary ray distribution and the scattered ray distribution, and radiation characteristics of an object that is located between the subject and a radiation detector that detects a radiation image; This may be performed by deriving a radiation image after transmission through the subject and the object using the second primary ray distribution and the scattered ray distribution.

[0019] In the estimation device according to the present disclosure, the scattered radiation removal process derives a region detection image by detecting, in two radiation images, a subject region where radiation has passed through the subject and reached the radiation detection unit, and a direct radiation region where radiation has directly reached the radiation detection unit without passing through the subject; deriving a scattered radiation image relating to the scattered radiation component based on the region detection image and the scattered radiation spread information relating to the spread of the scattered radiation; This may be performed by subtracting the scattered radiation image from the two radiographic images.

[0020] In the estimation device according to the present disclosure, the learning enhanced image is deriving processing details of a first graininess reduction process for a first radiographic image having a higher S / N ratio of the two radiographic images; deriving processing details of a second graininess suppression process for a second radiographic image having a low S / N ratio based on processing details of the first graininess suppression process; performing graininess suppression processing on the first radiographic image based on the processing content of the first graininess suppression processing; performing graininess suppression processing on the second radiographic image based on the processing content of the second graininess suppression processing; It may be derived using two radiation images that have been subjected to graininess suppression processing.

[0021] In the estimation device according to the present disclosure, the processing content of the first graininess suppression processing is as follows: The image may be derived based on a physical quantity map of the subject that is derived based on at least one of the first radiographic image and the second radiographic image.

[0022] The estimation method according to the present disclosure is an estimation method for deriving an estimation result of at least one emphasized image in which a specific composition of a subject is emphasized from a plain radiographic image, using a trained neural network that derives an estimation result of at least one emphasized image in which a specific composition of the subject is emphasized from a plain radiographic image acquired by plain radiographing a subject including a plurality of compositions, the estimation method comprising: The trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training enhanced image derived from the two radiation images, which enhances a specific composition of the subject.

[0023] The estimation method according to the present disclosure may be provided as a program for causing a computer to execute the method. Effect of the Invention

[0024] According to the present disclosure, it is possible to estimate with high accuracy an image in which a specific composition is emphasized. [Brief description of the drawings]

[0025] [Figure 1] FIG. 1 is a schematic block diagram showing a configuration of a radiation image capturing system to which an estimation device according to a first embodiment of the present disclosure is applied; [Diagram 2] FIG. 1 is a diagram showing a schematic configuration of an estimation device according to a first embodiment. [Diagram 3] FIG. 1 is a diagram showing a functional configuration of an estimation device according to a first embodiment. [Figure 4] FIG. 1 is a diagram showing a schematic configuration of a neural network used in the present embodiment; [Diagram 5] Diagram showing training data [Figure 6] FIG. 1 is a diagram showing a schematic configuration of an information derivation device according to a first embodiment; [Figure 7] FIG. 1 is a diagram showing a functional configuration of an information derivation device according to a first embodiment; [Figure 8] Figure showing bone images [Figure 9] Soft tissue images [Figure 10] Diagram to explain neural network learning [Figure 11] Conceptual diagram of the processing performed by a trained neural network [Figure 12] A diagram showing the display screen of the estimation results [Figure 13] 1 is a flowchart of a learning process performed in the first embodiment. [Figure 14] 1 is a flowchart of an estimation process performed in the first embodiment. [Figure 15] FIG. 13 is a diagram showing a functional configuration of an information derivation device according to a second embodiment. [Figure 16] FIG. 13 is a diagram for explaining calculation of an index value representing attenuation [Figure 17] FIG. 13 is a diagram for explaining calculation of an index value representing attenuation [Figure 18] FIG. 13 is a diagram showing a functional configuration of an information derivation device according to a third embodiment. [Figure 19] FIG. 1 is a table showing the relationship between body thickness distribution and initial weighting coefficients. [Figure 20]FIG. 13 is a table defining the relationship between pixel values ​​of bones and thicknesses of bones. [Figure 21] FIG. 13 is a table showing the relationship between the thickness of the soft part and the thickness of the bone part and the weighting coefficient. [Figure 22] FIG. 13 is a diagram showing a functional configuration of an information derivation device according to a fourth embodiment. [Figure 23] FIG. 1 is a table showing the relationship between the initial value of the thickness of the soft part and the initial value of the attenuation coefficient of the soft part. [Figure 24] FIG. 1 shows a table defining the relationship between the thickness of soft tissue and bone and the attenuation coefficient. [Diagram 25] FIG. 13 is a diagram showing a functional configuration of an information derivation device according to a fifth embodiment. [Figure 26] FIG. 1 is a diagram for explaining the difference in body thickness derived from low-energy images and high-energy images. [Figure 27] FIG. 1 shows a table correlating the difference in body thickness derived from two radiographic images with the composition ratio of fat. [Figure 28] FIG. 13 is a diagram showing a functional configuration of a scattered radiation removal unit in an information derivation device according to a sixth embodiment. [Figure 29] FIG. 23 is a diagram for explaining photography in the sixth embodiment; [Diagram 30] FIG. 1 is a diagram for explaining measurement of scattered radiation transmittance according to the subject's body thickness. [Diagram 31] FIG. 1 is a diagram for explaining measurement of scattered radiation transmittance according to the subject's body thickness. [Diagram 32] A table showing the relationship between the subject's body thickness distribution and the scattered radiation transmittance of an object between the subject and the radiation detector. [Diagram 33] FIG. 1 is a diagram for explaining measurement of primary ray transmittance according to the subject's body thickness. [Diagram 34] FIG. 1 is a diagram for explaining measurement of primary ray transmittance according to the subject's body thickness. [Diagram 35] A table showing the relationship between the subject's body thickness distribution and the primary ray transmittance of an object between the subject and the radiation detector. [Diagram 36] A diagram showing the state where an air gap exists between the top plate and the grid. [Figure 37] FIG. 13 is a diagram showing a functional configuration of a scattered radiation removal unit in an information derivation device according to a seventh embodiment. [Figure 38] FIG. 13 is a diagram for explaining a function of a first lead-out portion in the seventh embodiment. [Figure 39] FIG. 1 is a diagram for explaining a method for estimating the body thickness of a subject. [Diagram 40] Diagram showing the point spread function [Diagram 41] Diagram showing the path of radiation [Diagram 42] FIG. 13 is a diagram showing a functional configuration of an information derivation device according to an eighth embodiment. [Diagram 43] FIG. 1 shows a bilateral filter for a first radiation image. [Diagram 44] FIG. 44 is a diagram showing a local region of a second radiographic image corresponding to the local region of the first radiographic image shown in FIG. 43; [Diagram 45] FIG. 13 is a diagram showing an example of a bilateral filter for a second radiation image; [Diagram 46] FIG. 13 is a diagram showing a functional configuration of an information derivation device according to a ninth embodiment. [Figure 47] FIG. 1 shows an example of a bilateral filter for a physical quantity map. [Figure 48] FIG. 23 is a diagram showing a functional configuration of a scattered radiation removal unit in an information derivation device according to a tenth embodiment. [Figure 49] FIG. 1 is a diagram for explaining detection of a subject region and a direct radiation region; [Figure 50] FIG. 13 is a diagram showing an area detection image of a specific line portion. [Figure 51] A diagram showing a scattered radiation image of a specific line portion. [Figure 52] A diagram showing the distribution of body thickness at specific line sections [Figure 53] FIG. 23 is a diagram showing a functional configuration of an information derivation device according to an eleventh embodiment. [Figure 54] FIG. 1 shows an example of the energy spectrum of radiation after passing through muscle tissue and radiation after passing through fat tissue. [Figure 55] Another example of training data DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0026] Hereinafter, an embodiment of the present disclosure will be described with reference to the drawings. Fig. 1 is a schematic block diagram showing a configuration of a radiographic image capturing system to which an estimation device according to a first embodiment of the present disclosure is applied. As shown in Fig. 1, the radiographic image capturing system according to the first embodiment includes an imaging device 1, an image storage system 9, an estimation device 10 according to the first embodiment, and an information derivation device 50. The imaging device 1, the estimation device 10, and the information derivation device 50 are connected to the image storage system 9 via a network (not shown).

[0027] The imaging device 1 is an imaging device capable of performing energy subtraction by a so-called one-shot method in which radiation such as X-rays emitted from a radiation source 3 and transmitted through a subject H is irradiated to a first radiation detector 5 and a second radiation detector 6 with different energies. During imaging, as shown in Fig. 1, the first radiation detector 5, a radiation energy converting filter 7 made of a copper plate or the like, and the second radiation detector 6 are arranged in this order from the side closest to the radiation source 3, and the radiation source 3 is driven. The first and second radiation detectors 5 and 6 are in close contact with the radiation energy converting filter 7.

[0028] As a result, the first radiation detector 5 acquires a first radiation image G1 of the subject H using low-energy radiation including so-called soft rays. Meanwhile, the second radiation detector 6 acquires a second radiation image G2 of the subject H using high-energy radiation from which soft rays have been removed. Therefore, the first radiation image G1 and the second radiation image G2 are acquired by capturing an image of the subject H using radiation with different energy distributions. The first and second radiation images G1, G2 are input to the estimation device 10. Both the first radiation image G1 and the second radiation image G1 are front images including the crotch area of ​​the subject H.

[0029] The first and second radiation detectors 5, 6 are capable of repeatedly recording and reading out radiation images, and may be so-called direct type radiation detectors that generate electric charges by directly receiving radiation, or so-called indirect type radiation detectors that convert radiation into visible light and then convert the visible light into an electric charge signal. As a method for reading out the radiation image signal, it is preferable to use a so-called TFT readout method in which the radiation image signal is read out by turning on and off a TFT (thin film transistor) switch, or a so-called optical readout method in which the radiation image signal is read out by irradiating the detector with readout light, but this is not limiting and other methods may also be used.

[0030] The imaging device 1 can also perform simple imaging of the subject H using only the first radiation detector 5 to obtain a simple radiation image G0, which is a simple two-dimensional image of the subject H. The imaging for obtaining the first and second radiation images G1, G2 is called energy subtraction imaging to distinguish it from simple imaging. In this embodiment, the first and second radiation images G1, G2 obtained by energy subtraction imaging are used as learning data, which will be described later. Furthermore, the simple radiation image G0 obtained by simple imaging is used to derive an estimation result of at least one emphasized image in which a specific composition of the subject H is emphasized, as will be described later.

[0031] The image storage system 9 is a system that stores image data of a radiation image acquired by the imaging device 1. The image storage system 9 extracts an image corresponding to a request from the estimation device 10 from the stored radiation images, and transmits the image to the device that has issued the request. A specific example of the image storage system 9 is a PACS (Picture Archiving and Communication System). In this embodiment, the image storage system 9 stores a large amount of teacher data for training a neural network, which will be described later.

[0032] Next, an estimation device according to a first embodiment will be described. First, a hardware configuration of the estimation device according to the first embodiment will be described with reference to FIG. 2. As shown in FIG. 2, the estimation device 10 is a computer such as a workstation, a server computer, or a personal computer, and includes a CPU (Central Processing Unit) 11, a non-volatile storage 13, and a memory 16 as a temporary storage area. The estimation device 10 also includes a display 14 such as a liquid crystal display, an input device 15 such as a keyboard and a mouse, and a network I / F (InterFace) 17 connected to a network (not shown). The CPU 11, the storage 13, the display 14, the input device 15, the memory 16, and the network I / F 17 are connected to a bus 18. The CPU 11 is an example of a processor in the present disclosure.

[0033] The storage 13 is realized by a hard disk drive (HDD), a solid state drive (SSD), a flash memory, etc. The storage 13 as a storage medium stores the estimation program 12A and the learning program 12B installed in the estimation device 10. The CPU 11 reads out the estimation program 12A and the learning program 12B from the storage 13, expands them in the memory 16, and executes the expanded estimation program 12A and the learning program 12B.

[0034] The estimation program 12A and the learning program 12B are stored in a storage device of a server computer connected to a network or in a network storage in an externally accessible state, and are downloaded and installed in the computer constituting the estimation device 10 upon request. Alternatively, they are recorded on a recording medium such as a DVD (Digital Versatile Disc) or a CD-ROM (Compact Disc Read Only Memory) and distributed, and are installed in the computer constituting the estimation device 10 from the recording medium.

[0035] Next, a functional configuration of the estimation device according to the first embodiment will be described. FIG. 3 is a diagram showing the functional configuration of the estimation device according to the first embodiment. As shown in FIG. 3, the estimation device 10 includes an image acquisition unit 21, an information acquisition unit 22, an estimation unit 23, a learning unit 24, and a display control unit 25. The CPU 11 executes an estimation program 12A to function as the image acquisition unit 21, the information acquisition unit 22, the estimation unit 23, and the display control unit 25. The CPU 11 executes a learning program 12B to function as the learning unit 24.

[0036] The image acquisition unit 21 causes the imaging device 1 to perform energy subtraction imaging of the subject H, thereby acquiring a first radiation image G1 and a second radiation image G2, which are, for example, front images of the crotch area of ​​the subject H, from the first and second radiation detectors 5, 6. When acquiring the first radiation image G1 and the second radiation image G2, imaging conditions are set, such as the imaging dose, radiation quality, tube voltage, SID (Source Image Receptor Distance) which is the distance between the radiation source 3 and the surfaces of the first and second radiation detectors 5, 6, SOD (Source Object Distance) which is the distance between the radiation source 3 and the surface of the subject H, and the presence or absence of an anti-scatter grid.

[0037] The SOD and SID are used to calculate the body thickness distribution as described below. The SOD is preferably obtained by, for example, a TOF (Time Of Flight) camera. The SID is preferably obtained by, for example, a potentiometer, an ultrasonic range finder, a laser range finder, or the like.

[0038] The imaging conditions may be set by the operator through input device 15. The set imaging conditions are stored in storage 13. The first and second radiographic images G1, G2 acquired by energy subtraction imaging and the imaging conditions are also transmitted to and stored in image storage system 9. The first and second radiographic images G1, G2 are used to derive training data, which will be described later.

[0039] The image acquiring section 21 also acquires a simple radiographic image G0, which is a front image of the crotch area of ​​the subject H, by causing the imaging device 1 to perform simple imaging of the subject H using only the first radiation detector 5.

[0040] In this embodiment, the first and second radiographic images G1, G2 and the simple radiographic image G0 may be acquired by a program separate from the estimation program 12A and stored in the storage 13. In this case, the image acquisition unit 21 acquires the first and second radiographic images G1, G2 and the simple radiographic image G0 stored in the storage 13 by reading them out from the storage 13 for processing.

[0041] The information acquisition unit 22 acquires teacher data for training a neural network, which will be described later, from the image storage system 9 via the network I / F 17.

[0042] The estimation unit 23 derives estimation results of a bone image in which the bones contained in the subject H are emphasized and a soft tissue image in which the soft tissue is emphasized from the simple radiographic image G0. For this purpose, the estimation unit 23 derives estimation results of the bone image and the soft tissue image using a trained neural network 23A that outputs a bone image and a soft tissue image when the simple radiographic image G0 is input. Note that in this embodiment, the subject from which the estimation results of the bone image and the soft tissue image are derived is an image of the vicinity of the hip joint of the subject H, but is not limited thereto.

[0043] The learning unit 24 constructs a trained neural network 23A by machine learning the neural network using the teacher data. Examples of the neural network include a simple perceptron, a multilayer perceptron, a deep neural network, a convolutional neural network, a deep belief network, a recurrent neural network, and a probabilistic neural network. In this embodiment, a convolutional neural network is used as the neural network.

[0044] Fig. 4 is a diagram showing a neural network used in this embodiment. As shown in Fig. 4, the neural network 30 includes an input layer 31, an intermediate layer 32, and an output layer 33. The intermediate layer 32 includes, for example, a plurality of convolution layers 35, a plurality of pooling layers 36, and a fully connected layer 37. In the neural network 30, the fully connected layer 37 exists in front of the output layer 33. In the neural network 30, the convolution layers 35 and the pooling layers 36 are alternately arranged between the input layer 31 and the fully connected layer 37.

[0045] The configuration of the neural network 30 is not limited to the example shown in Fig. 4. For example, the neural network 30 may include one convolution layer 35 and one pooling layer 36 between the input layer 31 and the fully connected layer 37.

[0046] Fig. 5 is a diagram showing an example of training data used for training the neural network. As shown in Fig. 5, training data 40 is made up of training data 41 and correct answer data 42. In this embodiment, the data input to trained neural network 23A to obtain estimation results of bone images and soft tissue images is plain radiographic image G0, but training data 41 includes two radiographic images, a first radiographic image G1 and a second radiographic image G2, acquired by the above-mentioned energy subtraction imaging.

[0047] The correct answer data 42 is a bone image and a soft tissue image near the target bone (i.e., the femur) of the subject from which the learning data 41 was obtained. The bone image Gb and the soft tissue image Gs, which are the correct answer data 42, are derived from the first and second radiation images G1 and G2 by an information derivation device 50. The information derivation device 50 will be described below.

[0048] 6 is a schematic block diagram showing the configuration of the information derivation device according to the first embodiment. As shown in FIG. 6, the information derivation device 50 according to the first embodiment is a computer such as a workstation, a server computer, or a personal computer, and includes a CPU 51, a non-volatile storage 53, and a memory 56 as a temporary storage area. The information derivation device 50 also includes a display 54 such as a liquid crystal display, an input device 55 consisting of a pointing device such as a keyboard and a mouse, and a network I / F 57 connected to a network (not shown). The CPU 51, the storage 53, the display 54, the input device 55, the memory 56, and the network I / F 57 are connected to a bus 58.

[0049] The storage 53 is realized by an HDD, an SSD, a flash memory, etc., like the storage 53. The storage 53 as a storage medium stores an information derivation program 52. The CPU 51 reads out the information derivation program 52 from the storage 53, expands it in the memory 56, and executes the expanded information derivation program 52.

[0050] Next, a functional configuration of the information derivation device according to the first embodiment will be described. Fig. 7 is a diagram showing the functional configuration of the information derivation device according to the first embodiment. As shown in Fig. 7, the information derivation device 50 according to the first embodiment includes an image acquisition unit 61, a scattered radiation removal unit 62, and a subtraction unit 63. Then, the CPU 51 executes the information derivation program 52, so that the CPU 51 functions as the image acquisition unit 61, the scattered radiation removal unit 62, and the subtraction unit 63.

[0051] The image acquisition unit 61 acquires the first radiographic image G1 and the second radiographic image G2, which become the learning data 41 stored in the image storage system 9. Note that the image acquisition unit 61 may acquire the first radiographic image G1 and the second radiographic image G2 by causing the imaging device 1 to capture an image of the subject H, similar to the image acquisition unit 21 of the estimation device 10.

[0052] The image acquiring unit 61 also acquires the imaging conditions under which the first and second radiographic images were acquired and which are stored in the image storage system 9. The imaging conditions include the imaging dose, tube voltage, SID, SOD, and the presence or absence of an anti-scatter grid under which the first radiographic image G1 and the second radiographic image G2 were acquired.

[0053] Here, each of the first radiographic image G1 and the second radiographic image G2 contains scattered radiation components based on radiation scattered within the subject H in addition to primary radiation components of the radiation that has passed through the subject H. For this reason, the scattered radiation removal section 62 removes the scattered radiation components from the first radiographic image G1 and the second radiographic image G2. For example, the scattered radiation removal section 62 may remove the scattered radiation components from the first radiographic image G1 and the second radiographic image G2 by applying the method described in JP 2015-043959 A. When using the method described in JP 2015-043959 A, the body thickness distribution of the subject H and the scattered radiation components for removing the scattered radiation components are derived simultaneously.

[0054] The removal of the scattered radiation components from the first radiographic image G1 will be described below, but the removal of the scattered radiation components from the second radiographic image G2 can be performed in the same manner. First, the scattered radiation removal unit 62 acquires a virtual model of the subject H having an initial body thickness distribution T0(x,y). The virtual model is data virtually representing the subject H in which the body thickness according to the initial body thickness distribution T0(x,y) is associated with the coordinate position of each pixel of the first radiographic image G1. The virtual model of the subject H having the initial body thickness distribution T0(x,y) may be stored in advance in the storage 53 of the information derivation device 50. In addition, the scattered radiation removal unit 62 may calculate the body thickness distribution T(x,y) of the subject H based on the SID and SOD included in the imaging conditions. In this case, the initial body thickness distribution T0(x,y) can be obtained by subtracting the SOD from the SID.

[0055] Next, the scattered radiation removal unit 62 generates an estimated image based on the virtual model by combining an estimated primary ray image obtained by estimating a primary ray image obtained by photographing the virtual model and an estimated scattered ray image obtained by estimating a scattered ray image obtained by photographing the virtual model, as an estimated image based on the first radiation image G1 obtained by photographing the subject H.

[0056] Next, the scattered radiation removal unit 62 corrects the initial body thickness distribution T0(x,y) of the virtual model so that the difference between the estimated image and the first radiographic image G1 becomes smaller. The scattered radiation removal unit 62 repeatedly generates the estimated image and corrects the body thickness distribution until the difference between the estimated image and the first radiographic image G1 satisfies a predetermined end condition. The scattered radiation removal unit 62 derives the body thickness distribution when the end condition is satisfied as the body thickness distribution T(x,y) of the subject H. In addition, the scattered radiation removal unit 62 removes the scattered radiation components contained in the first radiographic image G1 by subtracting the scattered radiation components when the end condition is satisfied from the first radiographic image G1.

[0057] The subtraction unit 63 performs energy subtraction processing to derive a bone image Gb from which the bones of the subject H are extracted and a soft tissue image Gs from which the soft tissue is extracted from the first and second radiographic images G1, G2. The bone image Gb and the soft tissue image Gs derived by the subtraction unit 63 from the first and second radiographic images G1, G2 are an example of a learning-use emphasized image in the present disclosure. Note that the first and second radiographic images G1, G2 in the following processing are ones from which scattered radiation components have been removed. When deriving the bone image Gb, the subtraction unit 63 performs weighted subtraction between corresponding pixels on the first and second radiographic images G1, G2 as shown in the following formula (1), thereby generating a bone image Gb from which the bones of the subject H included in each of the radiographic images G1, G2 are extracted, as shown in FIG. 8. In formula (1), α is a weighting coefficient. The pixel value of each pixel in the bone region in the bone image Gb is the bone pixel value. Gb(x,y)=α·G2(x,y)-G1(x,y) (1)

[0058] On the other hand, when deriving a soft tissue image Gs, the subtraction unit 63 performs weighted subtraction between corresponding pixels of the first and second radiographic images G1, G2 as shown in the following equation (2) to generate a soft tissue image Gs from which the soft tissue of the subject H contained in each of the radiographic images G1, G2 has been extracted, as shown in Fig. 9. In equation (2), β is a weighting coefficient. Gs(x, y)=G1(x, y)-β×G2(x, y) (2)

[0059] The soft tissue image Gs represents a soft tissue region of the soft tissue of the subject H. In this embodiment, the "soft tissue" of the subject H refers to tissue other than bone tissue, and specifically includes muscle tissue, fat tissue, blood, and water.

[0060] The bone images Gb and soft tissue images Gs used as the correct answer data 42 are derived at the same time as the learning data 41 is acquired, and are transmitted to the image storage system 9. In the image storage system 9, the learning data 41 and the correct answer data 42 are associated with each other and stored as the teacher data 40. Note that, in order to improve the robustness of learning, teacher data 40 including images obtained by performing at least one of the following on the same image: enlargement / reduction, contrast change, movement, in-plane rotation, inversion, and noise addition, may be additionally created and stored as the learning data 41.

[0061] Returning to the estimation device 10, the learning unit 24 learns the neural network using a large number of teacher data 40. FIG. 10 is a diagram for explaining the learning of the neural network 30. When learning the neural network 30, the learning unit 24 inputs learning data 41, i.e., the first and second radiographic images G1 and G2, to the input layer 31 of the neural network 30. Then, the learning unit 24 causes the output layer 33 of the neural network 30 to output a bone image and a soft tissue image as output data 47. Then, the learning unit 24 derives the difference between the output data 47 and the correct answer data 42 as a loss L0. Note that the loss is derived between the bone image of the output data 47 and the bone image of the correct answer data 42, and between the soft tissue image of the output data and the soft tissue image of the correct answer data 42, but the loss is referred to as L0.

[0062] The learning unit 24 learns the neural network 30 based on the loss L0. Specifically, the learning unit 24 adjusts the kernel coefficients in the convolution layer 35, the weights of the connections between the layers, and the weights of the connections in the fully connected layer 37 (hereinafter referred to as parameters 48) so as to reduce the loss L0. For example, the backpropagation method can be used as a method for adjusting the parameters 48. The learning unit 24 repeats the adjustment of the parameters 48 until the loss L0 becomes equal to or less than a predetermined threshold value. As a result, when a simple radiographic image G0 is input, the parameters 48 are adjusted so as to output a bone image Gb and a soft tissue image Gs for the input simple radiographic image G0, and the trained neural network 23A is constructed. The constructed trained neural network 23A is stored in the storage 13.

[0063] Fig. 11 is a conceptual diagram of the process performed by the trained neural network 23A. As shown in Fig. 11, when a simple radiographic image G0 of a patient is input to the trained neural network 23A constructed as described above, the trained neural network 23A outputs a bone image Gb and a soft tissue image Gs for the input simple radiographic image G0.

[0064] The display control unit 25 displays the estimation results of the bone image Gb and the soft tissue image Gs estimated by the estimation unit 23 on the display 14. FIG. 12 is a diagram showing a display screen of the estimation results. As shown in FIG. 12, a display screen 70 has a first image display area 71 and a second image display area 72. A simple radiographic image G0 of the subject H is displayed in the first image display area 71. Furthermore, the bone image Gb and the soft tissue image Gs estimated by the estimation unit 23 are displayed in the second image display area 72.

[0065] Next, the process performed in the first embodiment will be described. FIG. 13 is a flowchart showing the learning process performed in the first embodiment. First, the information acquisition unit 22 acquires the teacher data 40 from the image storage system 9 (step ST1), and the learning unit 24 inputs the learning data 41 included in the teacher data 40 to the neural network 30 to output the bone image Gb and the soft tissue image Gs, and learns the neural network 30 using the loss L0 based on the difference from the correct answer data 42 (step ST2), and returns to step ST1. Then, the learning unit 24 repeats the processes of steps ST1 and ST2 until the loss L0 becomes a predetermined threshold value, and ends the learning process. Note that the learning unit 24 may end the learning process by repeating the learning a predetermined number of times. In this way, the learning unit 24 constructs the trained neural network 23A.

[0066] Next, the estimation process in the first embodiment will be described. FIG. 14 is a flowchart showing the estimation process in the first embodiment. It is assumed that the simple radiographic image G0 is acquired by imaging and stored in the storage 13. When an instruction to start the process is input from the input device 15, the image acquisition unit 21 acquires the simple radiographic image G0 from the storage 13 (step ST11). Next, the estimation unit 23 derives an estimation result of the bone image Gb and the soft tissue image Gs from the simple radiographic image G0 (step ST12). Then, the display control unit 25 displays the estimation result of the bone image Gb and the soft tissue image Gs derived by the estimation unit 23 on the display 14 together with the simple radiographic image G0 (step ST13), and the process ends.

[0067] In this embodiment, the trained neural network 23A constructed by learning using the first and second radiographic images G1 and G2 as teacher data is used to derive the estimation results of the bone image Gb and the soft tissue image Gs for the simple radiographic image G0. Here, in this embodiment, two radiographic images, the first and second radiographic images G1 and G2, are used for training the neural network. Therefore, compared with the case where one radiographic image and the bone image Gb and the soft tissue image Gs are used as teacher data, the trained neural network 23A can derive the estimation results of the bone image Gb and the soft tissue image Gs from the simple radiographic image G0 with higher accuracy. Therefore, according to this embodiment, the estimation results of the bone image Gb and the soft tissue image Gs can be derived with higher accuracy.

[0068] In the first embodiment, when deriving the bone image Gb and the soft tissue image Gs of the correct answer data 42, the bone and soft tissue of the subject H are recognized using at least one of the first and second radiographic images G1 and G2, and the bone and soft tissue attenuation coefficients for the bone and soft tissue are derived using the recognition results of the bone and soft tissue and the first and second radiographic images G1 and G2, and the bone image Gb and the soft tissue image Gs may be derived by performing energy subtraction processing using the derived attenuation coefficients. Hereinafter, this will be described as a second embodiment. Note that the energy subtraction processing in the second embodiment is described in, for example, International Publication No. 2020 / 175319.

[0069] Fig. 15 is a diagram showing the functional configuration of an information derivation device according to the second embodiment. In Fig. 15, the same components as those in Fig. 7 are given the same reference numerals, and detailed explanations are omitted. As shown in Fig. 15, the information derivation device 50A according to the second embodiment further includes a structure recognition unit 65 and a weighting coefficient derivation unit 66 in addition to the information derivation device 50 according to the first embodiment.

[0070] The structure recognition unit 65 recognizes structures included in the subject H using at least one of the first and second radiographic images G1, G2. In the second embodiment, the structures to be recognized are bones and soft tissues. In the second embodiment, the structure recognition unit 65 uses both the first and second radiographic images G1, G2 acquired by the image acquisition unit 61 for the recognition process, but may use either one of the first and second radiographic images G1, G2 for the recognition process. The structure recognition unit 65 may recognize structures using the first and second radiographic images G1, G2 after the scattered radiation removal process, or may recognize structures using the first and second radiographic images G1, G2 before the scattered radiation removal process.

[0071] The structure recognition unit 65 recognizes the position, size and / or shape, etc. of structures included in the subject H captured in the first and second radiographic images G1, G2. That is, the recognition process performed by the structure recognition unit 65 is a process for identifying the position, size and / or shape, etc. of structures having boundaries with other tissues, etc., in the subject H captured in the radiographic images. In the second embodiment, the bones and soft tissues of the subject H are recognized as structures.

[0072] The weighting coefficient derivation unit 66 derives index values ​​representing the attenuation of radiation for the structures recognized by the structure recognition unit 65, using the recognition result of the structure recognition unit 65 and the first and second radiographic images G1, G2, as weighting coefficients α, β to be used in subtraction processing. The attenuation coefficient is a so-called radiation source attenuation coefficient, and represents the degree (proportion) of attenuation of radiation due to absorption, scattering, etc. The attenuation coefficient varies depending on the specific composition (density, etc.) and thickness (mass) of the structure through which the radiation passes.

[0073] The weighting coefficient derivation unit 66 calculates an index value representing attenuation using a ratio or difference between pixel values ​​of corresponding pixels in the first and second radiographic images G1 and G2. FIG. 16 and FIG. 17 are diagrams for explaining the calculation of index values ​​representing attenuation. In FIG. 16 and FIG. 17, the calculation of index values ​​representing attenuation for three structures is explained. As shown in FIG. 16, it is assumed that there are three types of structures with compositions of "Ca", "Cb" and "Cc" in the first radiographic image G1, and the pixel values ​​of these in the first radiographic image G1 are all "V1". Meanwhile, the pixel values ​​of corresponding pixels in the second radiographic image G2 are "Va", "Vb" and "Vc", respectively. The degree of reduction in this pixel value corresponds to the degree of attenuation of radiation by each structure (each composition). 17, the weighting coefficient derivation unit 66 can calculate an index value μa representing the attenuation of a structure of composition "Ca", an index value μb representing the attenuation of a structure of composition "Cb", and an index value μc representing the attenuation of a structure of composition "Cc" by using a ratio or difference between a pixel value of a first radiographic image G1 and a corresponding pixel value of a second radiographic image G2. In the second embodiment, the weighting coefficient derivation unit 66 calculates an index value representing the attenuation of a bone portion and an index value representing the attenuation of a soft portion by defining the compositions as a bone portion and a soft portion.

[0074] Note that the index value μ representing attenuation can be calculated by knowing the ratio or difference between the pixel value in the first radiographic image G1 and the pixel value of the corresponding pixel in the second radiographic image G2. For this reason, although the pixel values ​​of the compositions "Ca", "Cb", and "Cc" are set to a common value "V1" in the first radiographic image G1 for the sake of explanation, the pixel values ​​of the compositions "Ca", "Cb", and "Cc" do not need to be common in the first radiographic image G1.

[0075] In this embodiment, the subtraction unit 63 derives the bone image Gb and the soft tissue image Gs. For this reason, the weighting coefficient derivation unit 66 derives the ratio Gs1(x,y) / Gs2(x,y) of the pixel value Gs1(x,y) of the first radiographic image G1, which is a low-energy image, to the pixel value Gs2(x,y) of the second radiographic image G2, which is a high-energy image, for the soft tissue regions in the first and second radiographic images G1 and G2 as an index value representing attenuation, i.e., the weighting coefficient α in the above formula (1). The ratio Gs1(x,y) / Gs2(x,y) represents the ratio μls / μhs of the attenuation coefficient μls for low-energy radiation to the attenuation coefficient μhs for high-energy radiation for soft tissue.

[0076] Furthermore, the weighting coefficient derivation unit 66 derives, for the bone regions in the first and second radiographic images G1, G2, a ratio Gb1(x,y) / Gb2(x,y) between the pixel value Gb1(x,y) of the first radiographic image G1, which is a low-energy image, and the pixel value Gb2(x,y) of the second radiographic image G2, which is a high-energy image, as an index value representing the attenuation coefficient, i.e., the weighting coefficient β in the above formula (2). Note that the ratio Gb1(x,y) / Gb2(x,y) represents the ratio μlb / μhb of the attenuation coefficient μlb for low-energy radiation to the attenuation coefficient μhb for high-energy radiation for the bone regions.

[0077] In the second embodiment, the subtraction unit 63 derives a bone image Gb and a soft tissue image Gs by performing energy subtraction processing according to the above formulas (1) and (2) using the weighting coefficients α and β derived by the weighting coefficient derivation unit 66. The weighting coefficients α and β in the formulas (1) and (2) are derived from the attenuation coefficients of the bones and the attenuation coefficients of the soft tissue derived by the weighting coefficient derivation unit 66.

[0078] Next, a third embodiment of the present disclosure will be described. FIG. 18 is a diagram showing a functional configuration of an information derivation device according to the third embodiment. In FIG. 18, the same components as those in FIG. 7 are given the same reference numbers, and detailed description will be omitted. In the third embodiment, when deriving the bone image Gb and the soft tissue image Gs that are the correct answer data 42, a new weighting factor to be used in weighted subtraction is derived based on the pixel values ​​of the bone included in the bone image Gb and the pixel values ​​of the soft tissue included in the soft tissue image Gs, and a new bone image and a new soft tissue image are derived by performing weighted subtraction using the new weighting factor on the first and second radiographic images G1 and G2, and a further new weighting factor is derived based on the new bone image and the new soft tissue image, and the derivation of a further new bone image and a further new soft tissue image based on the further new weighting factor is repeated, thereby deriving the bone image Gb and the soft tissue image Gs.

[0079] For this purpose, as shown in FIG. 18, the information derivation device 50B according to the third embodiment further includes an initial weighting coefficient setting unit 67 and a weighting coefficient derivation unit 68 in addition to the components of the information derivation device 50 according to the first embodiment.

[0080] The initial weighting coefficient setting unit 67 sets initial weighting coefficients, which are initial values ​​of weighting coefficients when the subtraction unit 63 performs subtraction processing, based on the body thickness distribution when the scattered radiation removal unit 62 satisfies the termination condition. Here, in the third embodiment, the subtraction unit 63 derives a bone image Gb and a soft tissue image Gs according to the above formulas (1) and (2) using the initial weighting coefficients set by the initial weighting coefficient setting unit 67 and the weighting coefficients α and β derived by the weighting coefficient derivation unit 68.

[0081] The initial weighting coefficient setting unit 67 sets initial weighting coefficients α0 and β0, which are initial values ​​of the weighting coefficients, based on the body thickness distribution when the scattered radiation removal unit 62 satisfies the termination condition. In the third embodiment, as shown in Fig. 19, a table LUT2 that specifies the relationship between the body thickness distribution and the initial weighting coefficients α0 and β0 is stored in the storage 53. The initial weighting coefficient setting unit 67 sets the initial weighting coefficients α0 and β0 based on the body thickness distribution with reference to the table LUT2.

[0082] Here, since the degree of the beam hardening described above depends on the thickness ts of the soft tissue and the thickness tb of the bone within the subject H, the attenuation coefficient μs of the soft tissue and the attenuation coefficient μb of the bone can be defined as μs(ts,tb) and μb(ts,tb) as functions of ts and tb.

[0083] In the energy subtraction process, since there are images having two different energy distributions, the attenuation coefficient of the soft tissue of the low-energy image (the first radiographic image G1 in this embodiment) can be expressed as μls(ts,tb), and the attenuation coefficient of the bone can be expressed as μlb(ts,tb). Also, the attenuation coefficient of the soft tissue of the high-energy image (the second radiographic image G2 in this embodiment) can be expressed as μhs(ts,tb), and the attenuation coefficient of the bone can be expressed as μhb(ts,tb).

[0084] In order to derive the bone image Gb, it is necessary to eliminate the contrast of the soft tissue contained in the radiographic image. Therefore, the weighting factor α can be calculated using the ratio of the attenuation coefficients of the soft tissue, α=μls(ts,tb) / μhs(ts,tb). In addition, in order to derive the soft tissue image Gs, it is necessary to eliminate the contrast of the bone contained in the radiographic image. Therefore, the weighting factor β can be calculated using the ratio of the attenuation coefficients of the bone, β=μlb(ts,tb) / μhb(ts,tb). In other words, the weighting factors α and β can be expressed as functions of the thickness ts of the soft tissue and the thickness tb of the bone.

[0085] In the third embodiment, the subtraction unit 63 first performs a subtraction process for performing weighted subtraction between corresponding pixels of the first and second radiographic images G1, G2 using the initial weighting coefficients α0, β0 set by the initial weighting coefficient setting unit 67. Thereafter, as will be described later, the subtraction process is performed using the weighting coefficients αnew, βnew derived by the weighting coefficient derivation unit 68.

[0086] The weighting factor derivation unit 68 derives new weighting factors αnew and βnew based on the pixel value Gb(x,y) of the bone contained in the bone image Gb. Here, the pixel value Gb(x,y) of the bone corresponds to the thickness of the bone of the subject H. For this reason, in the third embodiment, a reference object simulating bones having various thicknesses is photographed in advance to obtain a radiographic image of the reference object as a reference radiographic image. Then, using the relationship between the pixel value of the region of the reference object in the reference radiographic image and the thickness of the reference object, a table that defines the relationship between the pixel value and the thickness of the bone is derived in advance and stored in the storage 53. FIG. 20 is a diagram showing a table that defines the relationship between the pixel value and the thickness of the bone. The table LUT3 shown in FIG. 20 indicates that the lower the pixel value Gb(x,y) of the bone (i.e., the higher the brightness), the thicker the bone.

[0087] The weighting coefficient derivation unit 68 derives the bone thickness tb at each pixel of the bone image Gb from each pixel value Gb(x,y) of the bone image Gb with reference to table LUT3. Note that since regions in the bone image Gb where no bone exists consist only of soft tissue, the bone thickness tb is 0. On the other hand, for pixels in the bone image Gb where the bone thickness tb is not 0, the weighting coefficient derivation unit 68 derives the soft tissue thickness ts by subtracting the bone thickness tb from the body thickness distribution when the scattered radiation removal unit 62 satisfied the termination condition.

[0088] As described above, the weighting coefficients α and β can be expressed as functions of the thickness ts of the soft part and the thickness tb of the bone part. In the third embodiment, a table that specifies the relationship between the thickness ts of the soft part and the thickness tb of the bone part and the weighting coefficients α and β is stored in the storage 53. FIG. 21 is a diagram showing a table that specifies the relationship between the thickness ts of the soft part and the thickness tb of the bone part and the weighting coefficients α and β. As shown in FIG. 21, the table LUT4 three-dimensionally expresses the relationship between the thickness ts of the soft part and the thickness tb of the bone part and the weighting coefficient α (or β). Here, in the table LUT4, the weighting coefficient α (or β) becomes smaller as the thickness ts of the soft part and the thickness tb of the bone part become larger.

[0089] In the third embodiment, a plurality of tables LUT4 are prepared according to the energy distribution of the radiation used during imaging and stored in the storage of the information derivation device 50B. The weighting coefficient derivation unit 68 acquires information on the energy distribution of the radiation used during imaging based on the imaging conditions, reads out the table LUT4 corresponding to the acquired energy distribution information from the storage, and uses it to derive the weighting coefficients. Then, the weighting coefficient derivation unit 68 derives new weighting coefficients αnew and βnew based on the derived bone thickness tb and soft tissue thickness ts by referring to the table LUT4.

[0090] The subtraction unit 63 uses the new weighting coefficients αnew, βnew derived by the weighting coefficient derivation unit 68 to derive a new bone image Gbnew and a new soft tissue image Gsnew according to the above equations (1) and (2).

[0091] Note that the new bone image Gbnew and the new soft tissue image Gsnew may be used as the final bone image Gb and soft tissue image Gs, but in the third embodiment, the derivation of the weighting coefficients α and β and the subtraction process are repeated.

[0092] That is, the weighting coefficient derivation unit 68 derives a new bone thickness tbnew based on the pixel values ​​of the bone in the new bone image Gbnew with reference to table LUT3. Then, the weighting coefficient derivation unit 68 derives a difference Δtb between the new bone thickness tbnew and the bone thickness tb obtained in the previous process, and determines whether the difference Δtb is less than a predetermined threshold value. If the difference Δtb is equal to or greater than the threshold value, the weighting coefficient derivation unit 68 derives a new soft tissue thickness tsnew from the new bone thickness tbnew, and further derives new weighting coefficients αnew and βnew with reference to table LUT4 based on the new bone thickness tbnew and the new soft tissue thickness tsnew.

[0093] The subtraction unit 63 further performs subtraction processing using the new weighting factors αnew and βnew, and derives a further new bone image Gbnew and a further new soft tissue image Gsnew.

[0094] Then, the weighting coefficient derivation unit 68 derives a new bone thickness tbnew based on the new bone image Gbnew, and derives the difference Δtb between the new bone thickness tbnew and the bone thickness tb obtained in the previous process.

[0095] In the third embodiment, the subtraction unit 63 and the weighting factor derivation unit 68 repeat the subtraction process and the derivation of the weighting factors αnew and βnew until the difference Δtb derived by the weighting factor derivation unit 68 becomes less than a predetermined threshold value.

[0096] In the third embodiment, the bone image Gb and the soft tissue image Gs when the difference Δtb is less than the threshold value are used as the correct answer data.

[0097] Next, a fourth embodiment of the present disclosure will be described. FIG. 22 is a diagram showing a functional configuration of an information derivation device according to the fourth embodiment. In FIG. 22, the same reference numbers are given to the same configurations as those in FIG. 7, and detailed description will be omitted. In the fourth embodiment of the present disclosure, when deriving the bone image Gb and the soft tissue image Gs that are the correct answer data 42, the attenuation coefficient for each different energy distribution of the soft tissue, the thickness of the soft tissue, the attenuation coefficient for each different energy distribution of the bone, and the thickness of the bone are changed from the initial value, and the difference between the value of the attenuation coefficient of the soft tissue × the thickness of the soft tissue + the attenuation coefficient of the bone × the thickness of the bone and each pixel value of the radiographic image is derived for each different energy distribution, and the attenuation coefficient of the soft tissue and the attenuation coefficient of the bone for each different energy distribution that minimizes the difference or is less than a predetermined threshold value is derived, and the attenuation coefficient of the soft tissue and the attenuation coefficient of the bone for each different energy distribution are derived, and the weighting coefficient derived based on the attenuation coefficient of the soft tissue and the attenuation coefficient of the bone is used to perform energy subtraction processing, thereby deriving the bone image Gb and the soft tissue image Gs.

[0098] For this reason, as shown in FIG. 22, an information derivation device 50C according to the fourth embodiment further includes an initial value derivation unit 81, an attenuation coefficient derivation unit 82, and a weighting coefficient derivation unit 83 in addition to the components of the information derivation device 50 according to the first embodiment.

[0099] The initial value derivation unit 81 derives initial values ​​of attenuation coefficients for different energy distributions of soft tissue, soft tissue thickness, attenuation coefficients for different energy distributions of bone, and bone thickness for deriving weighting coefficients when performing energy subtraction processing. Specifically, it derives initial values ​​μls0, μhs0, ts0, μlb0, μhb0, tb0 of the soft-tissue attenuation coefficient μls for low-energy radiation, the soft-tissue attenuation coefficient μhs for high-energy radiation, the soft-tissue thickness ts, the bone attenuation coefficient μlb for low-energy radiation, the bone attenuation coefficient μhb for high-energy radiation, and the bone thickness tb.

[0100] In the fourth embodiment, the subtraction unit 63 performs a subtraction process using weighting coefficients derived by the weighting coefficient derivation unit 83, as described below, to perform weighted subtraction between corresponding pixels of the first and second radiographic images G1, G2, as shown in the above equations (1) and (2), thereby deriving a bone image Gb from which bones in the subject H are extracted and a soft tissue image Gs from which soft tissue is extracted.

[0101] Here, as described above, the attenuation coefficient μs of the soft tissue and the attenuation coefficient μb of the bone can be defined as μs(ts,tb) and μb(ts,tb) as functions of ts and tb. Therefore, in the fourth embodiment, as in the third embodiment, the attenuation coefficient of the soft tissue of the low-energy image (the first radiographic image G1 in this embodiment) can be expressed as μls(ts,tb) and the attenuation coefficient of the bone can be expressed as μlb(ts,tb). In addition, the attenuation coefficient of the soft tissue of the high-energy image (the second radiographic image G2 in this embodiment) can be expressed as μhs(ts,tb) and the attenuation coefficient of the bone can be expressed as μhb(ts,tb).

[0102] In order to derive the bone image Gb, it is necessary to eliminate the contrast of the soft tissue contained in the radiographic image. Therefore, the weighting factor α can be calculated by using the ratio of the attenuation coefficients of the soft tissue, α=μls(ts,tb) / μhs(ts,tb). In addition, in order to derive the soft tissue image Gs, it is necessary to eliminate the contrast of the bone contained in the radiographic image. Therefore, the weighting factor β can be calculated by using the ratio of the attenuation coefficients of the bone, β=μlb(ts,tb) / μhb(ts,tb). In the following description, the attenuation coefficients μls(ts,tb), μhs(ts,tb), μlb(ts,tb), and μhb(ts,tb) will be simply expressed as attenuation coefficients μls, μhs, μlb, and μhb, with (ts,tb) omitted.

[0103] The initial value derivation unit 81 uses the body thickness distribution when the scattered radiation removal unit 62 satisfies the end condition as the initial value ts0 of the soft tissue thickness ts. In the fourth embodiment, the body thickness distribution used when performing the scattered radiation removal process is the body thickness distribution assuming that the subject H is composed of only soft tissue. Therefore, the initial value tb0 of the bone thickness tb is 0. In addition, as the initial values ​​μls0, μhs0, μlb0, μhb0 of the attenuation coefficients, values ​​corresponding to the initial values ​​ts0, tb0 of the soft tissue and bone thicknesses ts, tb are derived. In the fourth embodiment, since the initial value tb0 of the bone thickness tb is 0, the attenuation coefficients μlb0, μhb0 of the bone are 0. As the soft tissue attenuation coefficients μls0, μhs0, values ​​corresponding to the initial value ts0 of the soft tissue thickness ts are derived. For this reason, in the fourth embodiment, a table is stored in the storage 53 that defines the relationship between the initial value ts0 of the soft-part thickness ts and the initial values ​​μls0 and μhs0 of the attenuation coefficients of the soft parts.

[0104] 23 is a diagram showing a table that specifies the relationship between the initial value ts0 of the soft-part thickness ts and the initial values ​​μls0 and μhs0 of the soft-part attenuation coefficients. The initial value derivation unit 81 derives the initial values ​​μls0 and μhs0 of the soft-part attenuation coefficients corresponding to the initial value ts0 of the soft-part thickness ts, with reference to the table LUT5 stored in the storage 53.

[0105] The attenuation coefficient deriving unit 82 derives the attenuation coefficients μls, μhs of the soft tissue and the attenuation coefficients μlb, μhb of the bone tissue for each different energy distribution. Here, for the energy subtraction process, a low energy image and a high energy image are obtained by capturing an image of the subject H with radiation having different energy distributions. In this embodiment, the first radiographic image G1 is a low energy image, and the second radiographic image G2 is a high energy image. The pixel values ​​G1(x,y) of each pixel of the first radiographic image G1, which is a low energy image, and the pixel values ​​G2(x,y) of each pixel of the second radiographic image G2, which is a high energy image, are expressed by the following formulas (3) and (4) using the thickness ts(x,y) of the soft tissue, the thickness tb(x,y) of the bone tissue, and the attenuation coefficients μls(x,y), μhs(x,y), μlb(x,y), and μhb(x,y) at the corresponding pixel positions. In addition, (x, y) is omitted in equations (3) and (4).

[0106] G1 = μls × ts + μlb × tb (3) G2 = μhs × ts + μhb × tb (4)

[0107] In order to derive the weighting coefficients α and β for performing the energy subtraction process, it is necessary to derive the attenuation coefficients μls(x,y), μhs(x,y), μlb(x,y), and μhb(x,y). As described above, the attenuation coefficients μls(x,y), μhs(x,y), μlb(x,y), and μhb(x,y) are expressed as functions of the thickness of the soft tissue ts and the thickness of the bone tb. Therefore, in order to derive the attenuation coefficients μls(x,y), μhs(x,y), μlb(x,y), and μhb(x,y), it is necessary to derive the thickness of the soft tissue ts and the thickness of the bone tb. Solving equations (3) and (4) for ts and tb gives the following equations (5) and (6).

[0108] ts={μhb×G1-μlb×G2} / {μls×μhb-μlb×μhs}(5) tb={μls×G2-μhs×G1} / {μls×μhb-μlb×μhs}(6)

[0109] Here, the attenuation coefficients μls(x,y), μhs(x,y), μlb(x,y), and μhb(x,y) on the right-hand sides of equations (5) and (6) are expressed as functions of the soft tissue thickness ts and the bone thickness tb, so equations (5) and (6) cannot be solved algebraically.

[0110] For this reason, in the fourth embodiment, the error functions EL and EH shown in the following formulas (7) and (8) are set. The error functions EL and EH correspond to the difference between the value of (soft tissue attenuation coefficient x soft tissue thickness + bone attenuation coefficient x bone thickness) and each pixel value of the radiographic image for each different energy distribution. In order to simultaneously minimize the error functions EL and EH, the error function E0 shown in formula (9) is set in the fourth embodiment. Then, while changing the soft tissue thickness ts, the bone thickness tb, and the attenuation coefficients μls, μhs, μlb, and μhb from their initial values, a combination of the soft tissue thickness ts and the bone thickness tb that minimizes the error function E0 or makes the error function E0 less than a predetermined threshold value is derived. In this case, it is preferable to derive the soft tissue thickness ts and the bone thickness tb using an optimization algorithm such as the steepest descent method and the conjugate gradient method. The initial values ​​of the soft part thickness ts, the bone part thickness tb and the attenuation coefficients μls, μhs, μlb, and μhb used in this case are ts0, tb0, μls0, μhs0, μlb0, and μhb0 derived by the initial value derivation unit 81.

[0111] EL = G1 - {μls × ts + μlb × tb} (7) EH = G2 - {μhs × ts + μhb × tb} (8) E0=EL 2 +EH 2 (9)

[0112] The attenuation coefficients used in the process of deriving the soft tissue thickness ts and the bone thickness tb are derived by referring to a table that defines the relationship between the soft tissue thickness ts and the bone thickness tb and the attenuation coefficients. Such a table is stored in the storage 53 of the information derivation device.

[0113] Fig. 24 is a diagram showing a table that specifies the relationship between the soft tissue thickness ts and the bone thickness tb and the attenuation coefficient. As shown in Fig. 24, the table LUT6 three-dimensionally represents the relationship between the soft tissue thickness ts and the bone thickness tb and the attenuation coefficient μ. Note that, although Fig. 24 shows only one LUT6, a table is prepared for each of the attenuation coefficients μls, μhs, μlb, and μhb and is stored in the storage. Here, in the table LUT6, the greater the soft tissue thickness ts and the bone thickness tb, the smaller the attenuation coefficient μ.

[0114] The attenuation coefficient derivation unit 82 derives the soft tissue thickness ts and the bone thickness tb when the error function E0 is minimized or becomes less than a predetermined threshold value, and then derives the attenuation coefficients μls, μhs, μlb, and μhb by referring to table LUT6.

[0115] The weighting coefficient derivation unit 83 derives weighting coefficients α and β used when the subtraction unit 63 performs subtraction processing. That is, the weighting coefficient derivation unit 83 derives weighting coefficients α and β by performing the calculations of α=μls / μhs and β=μlb / μhb using the attenuation coefficients μls, μhs, μlb, and μhb derived by the attenuation coefficient derivation unit 82.

[0116] In the fourth embodiment, the subtraction unit 63 uses the weighting coefficients α and β derived by the weighting coefficient derivation unit 83 to derive the bone image Gb and the soft tissue image Gs according to the above equations (1) and (2).

[0117] Next, a fifth embodiment of the present disclosure will be described. FIG. 25 is a diagram showing a functional configuration of an information derivation device according to the fifth embodiment. In FIG. 25, the same components as those in FIG. 7 are given the same reference numbers, and detailed description will be omitted. In the fifth embodiment of the present disclosure, when deriving the bone image Gb and the soft tissue image Gs that are the correct answer data 42, the composition ratios of a plurality of compositions contained in the soft tissue of the subject H are derived, and the attenuation coefficients for different energy distributions in the soft tissue are derived for each pixel of the first and second radiographic images G1 and G2 according to the composition ratios, and subtraction processing is performed using the derived soft tissue attenuation coefficients and a weighting coefficient derived based on a predetermined bone attenuation coefficient, thereby deriving the bone image Gb and the soft tissue image Gs.

[0118] For this purpose, as shown in FIG. 25, an information derivation device 50D according to the fifth embodiment further includes a composition ratio derivation unit 84 and an attenuation coefficient setting unit 85 in addition to the components of the information derivation device 50 according to the first embodiment.

[0119] The composition ratio derivation unit 84 acquires the composition ratio of the subject H. In the fifth embodiment, the composition ratio derivation unit 84 acquires the composition ratio by deriving the composition ratio of the subject H based on the first and second radiographic images G1, G2. Note that in the fifth embodiment, the composition ratio of fat is derived as the composition ratio. Therefore, although the subject H includes bones, for the sake of explanation, the first and second radiographic images G1, G2 will be described as including only soft tissues and not bones.

[0120] In order to derive the composition ratio, the composition ratio deriving unit 84 first derives the body thickness of the subject H as the first body thickness and the second body thickness for each pixel of the first and second radiographic images G1 and G2, respectively. Specifically, the composition ratio deriving unit 84 assumes that the luminance distribution of the first radiographic image G1 matches the body thickness distribution of the subject H, and converts the pixel values ​​of the first radiographic image G1 into thicknesses using the attenuation coefficient of radiation in the muscles of the subject H, thereby deriving the first body thickness t1 of the subject H. In addition, the composition ratio deriving unit 84 assumes that the luminance distribution of the second radiographic image G2 matches the body thickness distribution of the subject H, and converts the pixel values ​​of the second radiographic image G2 into thicknesses using the attenuation coefficient of the muscles of the subject H, thereby deriving the second body thickness t2 of the subject H.

[0121] Here, since the degree of the beam hardening described above depends on the fat thickness tf and muscle thickness tm in the subject H, the fat attenuation coefficient μf and the muscle attenuation coefficient μm can be defined as μf(tf,tm) and μm(tf,tm) as nonlinear functions of the fat thickness tf and muscle thickness tm.

[0122] In the fifth embodiment, the attenuation coefficient of fat for the first radiation image G1, which is a low-energy image, can be expressed as μlf(tf,tm), and the attenuation coefficient of muscle can be expressed as μlm(tf,tm). Also, the attenuation coefficient of fat for the second radiation image G2, which is a high-energy image, can be expressed as μhf(tf,tm), and the attenuation coefficient of muscle can be expressed as μhm(tf,tm).

[0123] Furthermore, the pixel value G1(x,y) of each pixel in the soft tissue region of the first radiographic image G1, which is a low-energy image, and the pixel value G2(x,y) of each pixel in the soft tissue region of the second radiographic image G2, which is a high-energy image, are expressed by the following formulas (10) and (11), respectively, using the fat thickness tf(x,y), muscle thickness tm(x,y), and attenuation coefficients μlf(x,y), μhf(x,y), μlm(x,y), and μhm(x,y) at the corresponding pixel positions. Note that (x,y) is omitted in formulas (10) and (11).

[0124] G1 = μlf × tf + μlm × tm (10) G2 = μhf × tf + μhm × tm (11)

[0125] As described above, in the fifth embodiment, when deriving the first body thickness t1 and the second body thickness t2, the pixel values ​​of the first radiographic image G1 and the second radiographic image G2 are converted into thicknesses using the attenuation coefficient of the muscle in the subject H. Therefore, in the fifth embodiment, the composition ratio deriving unit 84 derives the first body thickness t1 and the second body thickness t2 using the following formulas (12) and (13). Note that (x, y) are omitted in formulas (12) and (13).

[0126] t1=G1 / μlm (12) t2=G2 / μhm (13)

[0127] If the subject H only contains muscle at the pixel positions where the first and second body thicknesses t1 and t2 are derived, the first body thickness t1 and the second body thickness t2 are the same. However, in the actual subject H, both muscle and fat are contained at the same pixel positions in the first and second radiation images G1 and G2. Therefore, the first and second body thicknesses t1 and t2 derived by the formulas (12) and (13) do not match the actual body thickness of the subject H. In addition, the first body thickness t1 derived from the first radiation image G1, which is a low-energy image, and the second body thickness t2 derived from the second radiation image G2, which is a high-energy image, the first body thickness t1 is greater than the second body thickness t2.

[0128] For example, as shown in Fig. 26, assume that the actual body thickness is 100 mm, and the thicknesses of fat and muscle are 30 mm and 70 mm, respectively. In this case, the first body thickness t1 derived from the first radiographic image G1 acquired using low-energy radiation is, for example, 80 mm, and the second body thickness t2 derived from the second radiographic image G2 acquired using high-energy radiation is, for example, 70 mm. The difference between the first body thickness t1 and the second body thickness t2 increases as the composition ratio of fat increases.

[0129] Here, the difference between the first body thickness t1 and the second body thickness t2 varies depending on the composition ratio of fat and muscle in the subject H. For this reason, in the fifth embodiment, subject models with variously changed composition ratios of fat are imaged using radiation with different energy distributions, and the body thickness is derived from each of the two radiographic images thus obtained. A table correlating the difference in body thickness derived from the two radiographic images with the composition ratio of fat is created in advance and stored in the storage 53.

[0130] Fig. 27 is a diagram showing a table associating the difference in body thickness derived from two radiographic images with the fat composition percentage. As shown in Fig. 27, the horizontal axis of table LUT7 represents the difference in body thickness derived from each of the two radiographic images, and the vertical axis represents the fat composition percentage. As shown in Fig. 27, the greater the difference in body thickness derived from each of the two radiographic images, the greater the fat composition percentage. Note that a table associating the difference in body thickness derived from each of the two radiographic images with the fat composition percentage is prepared for each energy distribution of radiation used during imaging, and is stored in storage 53.

[0131] The composition ratio derivation unit 84 derives the difference between the derived first body thickness t1 and the derived second body thickness t2, and derives the fat composition ratio R(x,y) by referring to the LUT7 stored in the storage of the information derivation device 50D. Note that the muscle composition ratio can be derived by subtracting the derived fat composition ratio R(x,y) from 100%.

[0132] The attenuation coefficient setting unit 85 sets the attenuation coefficient of radiation used when acquiring the first and second radiographic images G1, G2 for each pixel of the first and second radiographic images G1, G2 according to the composition ratio R(x, y) of fat. Specifically, the attenuation coefficient of the soft tissue of the subject H is set. Here, in the fifth embodiment, the first radiographic image G1 corresponds to a low-energy image, and the second radiographic image G2 corresponds to a high-energy image. For this reason, the attenuation coefficient setting unit 85 derives the soft tissue attenuation coefficient μls(x, y) for the low-energy image and the soft tissue attenuation coefficient μhs(x, y) for the high-energy image by the following formulas (14) and (15). Note that (x, y) is omitted in formulas (14) and (15). Furthermore, μlm is the attenuation coefficient of muscle in low-energy images, μlf is the attenuation coefficient of fat in low-energy images, μhm is the attenuation coefficient of muscle in high-energy images, and μhf is the attenuation coefficient of fat in high-energy images.

[0133] μls = (1-R) ​​× μlm + R × μlf (14) μhs = (1-R) ​​× μhm + R × μhf (15)

[0134] In the fifth embodiment, the subtraction unit 63 performs subtraction processing using the attenuation coefficients μls and μhs set by the attenuation coefficient setting unit 85. At this time, the subtraction unit 63 derives the weighting coefficients α and β in the above formulas (1) and (2). The weighting coefficient α used in the fifth embodiment is derived by α=μls / μhs using the attenuation coefficient set by the attenuation coefficient setting unit 85. Note that in the fifth embodiment, the attenuation coefficients μlb and μhb of the bones are prepared in advance and stored in the storage 53. Then, the subtraction unit 63 derives the bone image Gb and the soft tissue image Gs according to the above formulas (1) and (2) using the derived weighting coefficients α and β.

[0135] Next, an information derivation device according to a sixth embodiment of the present disclosure will be described. The configuration of the information derivation device according to the sixth embodiment is the same as that of the information derivation device 50 according to the first embodiment, and only the processing performed by the scattered radiation removal unit 62 is different, so detailed description will be omitted here. FIG. 28 is a diagram showing the functional configuration of the scattered radiation removal unit in the information derivation device according to the sixth embodiment. As shown in FIG. 28, the scattered radiation removal unit 62A in the information derivation device according to the sixth embodiment includes an imaging condition acquisition unit 91, a body thickness derivation unit 92, a characteristic acquisition unit 93, a line distribution derivation unit 94, and a calculation unit 95.

[0136] In the sixth embodiment, an imaging device 1A shown in Fig. 29 is used to image a subject H lying on his back on an imaging table 4, and an anti-scatter grid 8 is further used to obtain first and second radiation images G1, G2, and the first and second radiation images G1, G2 are used as teacher data 40. In Fig. 29, the imaging table 4 has a tabletop 4A, and below the tabletop 4A, an anti-scatter grid (hereinafter simply referred to as grid) 8, a first radiation detector 5, a filter 7, and a second radiation detector 6 are arranged in this order from the radiation source 3 side, and these are attached below the tabletop 4A by attachment parts 4B.

[0137] When generating the teacher data 40, the imaging condition acquisition unit 91 acquires from the image storage system 9 imaging conditions used in the energy subtraction imaging of the subject H to acquire the first and second radiation images G1, G2 as the learning data 41. In the sixth embodiment, the imaging conditions also include the radiation quality. The radiation quality is defined using one or more of the tube voltage [kV] of the radiation generator in the radiation source 3, the total filtration amount [mmAl equivalent], and the half value layer [mmAl]. The tube voltage means the maximum value of the generated radiation energy distribution. The total filtration amount is the filtration amount of each component of the imaging device 1, such as the radiation generator and collimator in the radiation source 3, converted into the thickness of aluminum. The larger the total filtration amount, the greater the influence of beam hardening in the imaging device 1, and the more high-energy components there are in the wavelength distribution of the radiation. The half value layer is defined by the thickness of aluminum required to attenuate the dose by half for the generated radiation energy distribution. The thicker the aluminum of the half value layer, the more high-energy components there are in the wavelength distribution of the radiation.

[0138] The body thickness derivation unit 92 derives the body thickness distribution of the subject H based on at least one of the first and second radiographic images G1, G2 and the imaging conditions. Hereinafter, the body thickness distribution derived by the body thickness derivation unit 92 is referred to as an initial body thickness distribution T0. Hereinafter, the derivation of the initial body thickness distribution T0 will be described. Note that, although the following description will be given of the derivation of the body thickness and the removal of scattered rays using the first radiographic image G1, the derivation of the body thickness and the removal of scattered rays can also be performed in a similar manner using the second radiographic image G2.

[0139] First, when the radiation source 3 is driven to irradiate the radiation detector 5 with radiation in the absence of the subject H, the reaching dose I0(x,y) of the radiation emitted from the radiation source 3 to the radiation detector 5 is expressed by the following formula (16). In formula (16), mAs included in the imaging conditions is the tube current time product, and kV is the tube voltage. When the half-value layer is also taken into consideration, the reaching dose I0(x,y) is expressed by the following formula (16-1). Here, F is a nonlinear function representing the amount of radiation that reaches the radiation detector 5 when a reference dose (e.g., 1 mAs) is irradiated to the radiation detector 5 in the absence of the subject H at a reference SID (e.g., 100 cm). F varies depending on the tube voltage or the tube voltage and the half-value layer. In addition, the reaching dose I0 is derived for each pixel of the radiation image acquired by the radiation detector 5, so (x,y) represents the pixel position of each pixel. In the following explanation, in order to include both cases where the half-value layer is taken into account and where it is not, the formulas will be expressed by including mmAl in parentheses, as shown in the following formula (16-2).

[0140] I0(x,y)=mAs×F(kV) / SID 2 (16) I0(x,y)=mAs×F(kV,mmAl) / SID 2 (16-1) I0(x,y)=mAs×F(kV(,mmAl)) / SID 2 (16-2)

[0141] In addition, if the initial body thickness distribution is T0, the attenuation coefficient of the subject H when the initial body thickness distribution is T0 is μ(T0), and the Scatter-to-Primary Ratio, which is the ratio of the scattered radiation amount and the primary radiation amount contained in the radiation after passing through the subject H having the initial body thickness distribution T0 when the spread of the scattered radiation is not considered, is STPR(T0), the dose I1 after passing through the subject H is expressed by the following formula (17). In addition, in formula (17), the initial body thickness distribution T0, the arriving dose I0, and the dose I1 are derived for each pixel of the simple radiographic image G0, but (x, y) are omitted. In addition, STPR is a nonlinear function that depends not only on the body thickness but also on the tube voltage [kV] and the half-value layer [mmAl], but in formula (17), the notations of kV and mmAl are omitted. I1=I0×exp{-μ(T0)×T0}×{1+STPR(T0)} (17)

[0142] In formula (17), the dose I1 is the pixel value of each pixel of the plain radiographic image G0, and the arrival dose I0 is derived by the above formulas (16) and (16-1). On the other hand, since F and STPR are nonlinear functions, formula (17) cannot be algebraically solved for T0. For this reason, the body thickness derivation unit 92 defines an error function E1 shown in the following formula (18) or formula (18-1). Then, T0 that minimizes the error function E1 or is less than a predetermined threshold value is derived as the initial body thickness distribution. At this time, the body thickness derivation unit 92 derives the initial body thickness distribution T0 using an optimization algorithm such as the steepest descent method and the conjugate gradient method.

[0143] E1=[I1-I0×exp{-μ(T0)×T0}×{1+STPR(T0)}] 2 (18) E1=|I1-I0×exp{-μ(T0)×T0}×{1+STPR(T0)}| (18-1)

[0144] The characteristic acquisition unit 93 acquires radiation characteristics of an object interposed between the subject H and the first and second radiation detectors 5, 6 during imaging. Here, when the radiation after passing through the subject H passes through an object interposed between the subject H and the radiation detector 5, the transmittance of the radiation changes depending on the radiation quality of the radiation after passing through the subject H. Furthermore, the primary rays and scattered rays contained in the radiation after passing through the subject H have different transmittances due to differences in the traveling direction and radiation quality of the radiation. For this reason, in the sixth embodiment, the primary ray transmittance and scattered ray transmittance of the object are used as the radiation characteristics of the object.

[0145] As described above, when the radiation after passing through the subject H passes through an object between the subject H and the radiation detector 5, the transmittance of the radiation changes depending on the radiation quality after passing through the subject H. In addition, the radiation quality after passing through the subject H depends on the body thickness distribution T of the subject H. For this reason, the primary ray transmittance and scattered ray transmittance can be expressed as functions of the body thickness distribution T of the subject H, as Tp(T) and Ts(T), respectively.

[0146] Furthermore, the radiation quality after passing through the subject H also depends on the radiation quality of the radiation source 3 included in the imaging conditions. The radiation quality depends on the tube voltage and the half-value layer. For this reason, the primary ray transmittance and the scattered ray transmittance are strictly expressed as Tp(kV(,mmAl),T) and Ts(kV(,mmAl),T), respectively. In the following explanation, the primary ray transmittance and the scattered ray transmittance may be simply expressed as Tp and Ts.

[0147] Here, the primary ray transmittance Tp and scattered ray transmittance Ts of an object interposed between the subject H and the radiation detector 5 depend on the body thickness distribution T of the subject H as described above. For this reason, in the sixth embodiment, a phantom having various thicknesses simulating the body thickness distribution T of the subject is used to measure the primary ray transmittance Tp and scattered ray transmittance Ts of the object according to the body thickness distribution T of the subject H, and a table that specifies the relationship between the body thickness distribution T of the subject H and the primary ray transmittance Tp and scattered ray transmittance Ts of the object is generated based on the measurement result, and stored in the storage of the information derivation device according to the sixth embodiment. Hereinafter, the measurement of the primary ray transmittance Tp and scattered ray transmittance Ts of the object according to the body thickness distribution T of the subject H will be described.

[0148] First, the calculation of the scattered radiation transmittance Ts will be described. Fig. 30 and Fig. 31 are diagrams for explaining the measurement of the scattered radiation transmittance Ts according to the body thickness of the subject H. First, as shown in Fig. 30, a phantom 101 simulating a human body is placed on the surface of the radiation detector 5, and a lead plate 102 is further placed on the phantom 101. Here, the phantom 101 has various thicknesses such as 5 cm, 10 cm, and 20 cm, and is made of a material such as acrylic having a radiation transmittance similar to that of water. In this state, the radiation source 3 is driven to irradiate the radiation detector 5 with radiation, and the characteristic acquisition unit 93 acquires a radiation image K0 for measurement. The signal value of the radiation image K0 is large in the area of ​​the radiation detector 5 directly irradiated with radiation, and the signal value decreases in the order of the area of ​​the phantom 101 and the area of ​​the lead plate 102.

[0149] Since the lead plate 102 does not transmit radiation, the signal value in the area of ​​the lead plate 102 in the radiation image K0 should be 0. However, radiation scattered by the phantom 101 reaches the area of ​​the radiation detector 5 corresponding to the lead plate 102. For this reason, the area of ​​the lead plate 102 in the radiation image K0 has a signal value S0 corresponding to the scattered radiation component by the phantom 101.

[0150] Next, as shown in FIG. 31, a phantom 101 is placed on the tabletop 4A of the imaging device 1, and a lead plate 102 is further placed on the phantom 101. Then, in the same manner as in imaging the subject H, a radiation detector 5 and a grid 8 are disposed below the tabletop 4A, and the radiation source 3 is driven to irradiate the radiation detector 5 with radiation, whereby the characteristic acquisition unit 93 acquires a radiation image K1 for measurement. As with the radiation image K0, the signal value of the radiation image K1 is large in the area of ​​the radiation detector 5 directly irradiated with radiation, and the signal value decreases in the order of the area of ​​the phantom 101 and the area of ​​the lead plate 102. Here, as shown in FIG. 31, when imaging is performed with the tabletop 4A and the grid 8 interposed between the phantom 101 and the radiation detector 5, not only the radiation scattered by the phantom 101 but also the radiation scattered by the tabletop 4A and the grid 8 reaches the area of ​​the radiation detector 5 corresponding to the lead plate 102. Therefore, the region of the lead plate 102 in the radiation image K1 has a signal value S1 corresponding to the scattered radiation components caused by the phantom 101, the table top 4A, and the grid 8.

[0151] Note that signal value S1 includes scattered ray components due to the tabletop 4A and grid 8, and is therefore greater than signal value S0 shown in Fig. 30. Therefore, when imaging a phantom 101 having a thickness t, the scattered ray transmittance Ts of the objects interposed between the subject H and the radiation detector 5, i.e., the tabletop 4A and grid 8, can be calculated by S1 / S0.

[0152] In the sixth embodiment, the characteristic acquisition unit 93 uses at least two types of phantoms with different thicknesses to calculate the scattered ray transmittance Ts corresponding to each thickness as shown in Fig. 30 and Fig. 31. Furthermore, the characteristic acquisition unit 93 derives the scattered ray transmittance Ts of a thickness not included in the phantom 101 by interpolating the scattered ray transmittances Ts for a plurality of measured thicknesses. In this way, the characteristic acquisition unit 93 generates a table LUT8 that indicates the relationship between the body thickness of the subject H and the scattered ray transmittance Ts of an object interposed between the subject H and the radiation detector 5 as shown in Fig. 32 by interpolating the scattered ray transmittance for a thickness between each thickness.

[0153] Next, calculation of the primary ray transmittance will be described. Fig. 33 and Fig. 34 are diagrams for explaining measurement of the primary ray transmittance Tp according to the body thickness of the subject H. First, as shown in Fig. 33, a phantom 101 simulating a human body is placed on the surface of the radiation detector 5. Here, the phantom 101 used is the same as that used when the scattered ray transmittance Ts was derived. In this state, the radiation source 3 is driven to irradiate the radiation detector 5 with radiation, and the characteristic acquisition unit 93 acquires a radiation image K2 for measurement. A signal value S2 of an area corresponding to the phantom 101 in the radiation image K2 includes both the primary ray component and the scattered ray component of the radiation that has passed through the phantom 101. Here, the scattered ray component of the radiation that has passed through the phantom 101 is the signal value S0 in the radiation image K0 obtained by the method shown in Fig. 30 above. Therefore, the primary ray component of the radiation that has passed through the phantom 101 is derived by S2-S0.

[0154] Next, as shown in Fig. 34, a phantom 101 is placed on the tabletop 4A of the imaging device 1, and in a similar manner to the case of imaging a subject H, a radiation detector 5 and a grid 8 are arranged below the tabletop 4A. In this state, the radiation source 3 is driven to irradiate the radiation detector 5 with radiation, and the characteristic acquisition unit 93 acquires a radiation image K3 for measurement. A signal value S3 of a region corresponding to the phantom 101 in the radiation image K3 includes both a primary ray component and a scattered ray component of the radiation that has passed through the phantom 101, the tabletop 4A, and the grid 8. Here, the scattered ray component of the radiation that has passed through the phantom 101, the tabletop 4A, and the grid 8 is the signal value S1 in the radiation image K1 obtained by the method shown in Fig. 31 above. For this reason, the primary ray component of the radiation that has passed through the phantom 101, the tabletop 4A, and the grid 8 is derived by S3-S1.

[0155] Therefore, when the phantom 101 is imaged, the primary ray transmittance Tp of the tabletop 4A and the grid 8 interposed between the subject H and the radiation detector 5 can be calculated by (S3-S1) / (S2-S0). In the sixth embodiment, the characteristic acquisition unit 93 uses at least two types of phantoms with different thicknesses to calculate the primary ray transmittance Tp corresponding to each thickness as shown in Figs. 33 and 34. Furthermore, the characteristic acquisition unit 93 derives the primary ray transmittance Tp of a thickness not included in the phantom 101 by interpolating the primary ray transmittances Tp for the measured thicknesses. As a result, the characteristic acquisition unit 93 generates a table LUT9 that indicates the relationship between the body thickness of the subject H and the primary ray transmittance Tp of an object interposed between the subject H and the radiation detector 5 as shown in Fig. 35.

[0156] The tables LUT8 and LUT9 generated as described above are stored in the storage 53 of the information derivation device according to the sixth embodiment. Note that the tables are generated according to various imaging conditions (i.e., radiation quality, dose, and radiation source distance) and further according to the type of grid 8 to be used, and are stored in the storage 13.

[0157] The characteristic acquisition unit 93 refers to tables LUT8 and LUT9 stored in the storage 53 in accordance with the imaging conditions acquired by the imaging condition acquisition unit 91, and acquires the primary ray transmittance Tp(T0) and scattered ray transmittance Ts(T0) corresponding to the initial body thickness distribution T0 for an object interposed between the subject H and the radiation detector 5. Note that since the primary ray transmittance Tp and the scattered ray transmittance Ts also depend on the radiation quality, the primary ray transmittance Tp and the scattered ray transmittance Ts can be expressed as Tp(kV(,mmAl),T0) and Ts(kV(,mmAl),T0), respectively.

[0158] The radiation distribution deriving unit 94 derives the primary radiation distribution and the scattered radiation distribution of the radiation detected by the radiation detector 5 using the imaging conditions, the body thickness distribution, and the radiation characteristics of the object interposed between the subject H and the radiation detector 5. Here, the primary radiation distribution Ip0 and the scattered radiation distribution Is0 of the radiation after passing through the subject H are expressed by the following formulas (19) and (20) when the body thickness distribution is T. The PSF in formula (20) is a point spread function that represents the distribution of scattered radiation spreading from one pixel, and is defined according to the radiation quality and the body thickness. Also, * indicates a convolution operation. The primary radiation distribution Ip0 and the scattered radiation distribution Is0 are derived for each pixel of the simple radiation image G0, but (x, y) are omitted in formulas (19) and (20). In addition, in the sixth embodiment, as described below, the body thickness distribution, the primary ray distribution Ip0, and the scattered ray distribution Is0 are repeatedly derived, but when deriving the primary ray distribution Ip0 and the scattered ray distribution Is0 for the first time, the initial body thickness distribution T0 is used as the body thickness distribution T.

[0159] Ip0=I0×exp{-μ(T)×T} (19) Is0=Ip0×STPR(kV(,mmAl),T)*PSF(kV(,mmAl),T) (20)

[0160] Furthermore, the ray distribution derivation unit 94 derives the primary ray distribution Ip1 and scattered ray distribution Is1 that reach the radiation detector 5 according to the following equations (21) and (22) using the primary ray transmittance Tp and scattered ray transmittance Ts of an object interposed between the subject H and the radiation detector 5. Furthermore, the sum Iw1 of the primary ray distribution Ip1 and scattered ray distribution Is1 is derived according to the following equation (23). In equations (21) and (22), the initial body thickness distribution T0 is used as the body thickness distribution T when deriving the primary ray distribution Ip1 and scattered ray distribution Is1 the first time.

[0161] Ip1 = Ip0 × Tp (kV (, mmAl), T) (21) Is1 = Is0 × Ts (kV(,mmAl),T) (22) Iw1 = Ip1 + Is1 (23)

[0162] The calculation unit 95 derives an error E2 between the sum Iw1 of the primary ray distribution Ip1 and the scattered ray distribution Is1 and the dose at each pixel position of the first radiographic image G1, i.e., the pixel value I1. The error E2 is derived using the following formula (24) or formula (24-1). In formula (24) and formula (24-1), N represents the number of pixels in the first radiographic image G1, and Σ represents the sum for all pixels in the first radiographic image G1. Note that formula (24-1) calculates I1 / Iw1 within log, and therefore the error E2 can be derived independently of the dose irradiated to the subject H, i.e., the reaching dose I0.

[0163] E2=(1 / N)×Σ{I1-Iw1} 2 (twenty four) E2=(1 / N)×Σ|log{I1 / Iw1}| (24-1)

[0164] Then, the calculation unit 95 updates the body thickness distribution T so as to minimize the error E2 or make the error E2 less than a predetermined threshold value. Then, the calculation unit 95 repeats the acquisition of the primary ray transmittance Tp and the scattered ray transmittance Ts and the derivation of the primary ray distribution Ip1 and the scattered ray distribution Is1 based on the updated body thickness distribution. Here, the calculation performed by the calculation unit 95 is referred to as a repeated calculation. Also, in the sixth embodiment, the calculation unit 95 repeats the calculation so that the error E2 becomes less than a predetermined threshold value. Then, the calculation unit 95 outputs a processed first radiographic image Gc1 having a pixel value that is the primary ray distribution Ipc derived based on the body thickness distribution Tc of the subject H in which the error E2 becomes less than a predetermined threshold value.

[0165] The repeated acquisition of the primary ray transmittance Tp and the scattered ray transmittance Ts, and the repeated derivation of the primary ray distribution Ip1 and the scattered ray distribution Is1 are performed by the characteristic acquisition unit 93 and the line distribution derivation unit 94, respectively.

[0166] Meanwhile, a primary ray distribution Ipc is derived for the second radiographic image G2 in the same manner as for the first radiographic image G1. The primary ray distribution for the first radiographic image G1 is designated Ipc-1, and the primary ray distribution for the second radiographic image G2 is designated Ipc-2. The calculation unit 95 then outputs a processed second radiographic image Gc2 having the primary ray distribution Ipc-2 as pixel values.

[0167] In the sixth embodiment, a subtraction unit 63 derives a bone image Gb and a soft tissue image Gs using the processed first and second radiation images Gc1, Gc2.

[0168] The trained neural network 23A constructed by learning using the bone image Gb and soft tissue image Gs derived in the sixth embodiment as correct answer data will output a bone image Gb and a soft tissue image Gs that take into account objects interposed between the subject H and the radiation detector when a simple radiation image G0 obtained by simply photographing the subject H using the imaging device 1A shown in FIG. 29 is input.

[0169] In the sixth embodiment, the tabletop 4A and the grid 8 of the imaging table 4 are used as objects between the subject H and the radiation detector 5. However, as shown in FIG. 36, there are cases where an air layer 103 is present between the tabletop 4A and the grid 8. In such a case, it is preferable that the line distribution deriving unit 94 derives the primary ray distribution Ip1 and the scattered ray distribution Is1 by including the air layer 103 in the objects between the subject H and the radiation detector 5. In this case, the primary ray distribution Ip1 and the scattered ray distribution Is1 may be derived by convoluting the point spread function PSFair(kV(,mmAl),tair) according to the thickness tair of the air layer 103 with the above equations (21) and (22), as shown in the following equations (21-1) and (22-1). The thickness tair of the air layer 103 is the distance between the lower surface of the tabletop 4A and the surface of the grid 8 on the subject H side.

[0170] Ip1=Ip0×Tp(kV(,mmAl),T)*PSF air (kV(,mmAl),t air ) (21-1) Is1 = Is0 × Ts (kV (, mmAl), T) * PSF air (kV(,mmAl),t air ) (22-1)

[0171] Next, an information derivation device according to a seventh embodiment of the present disclosure will be described. The configuration of the information derivation device according to the seventh embodiment is the same as that of the information derivation device according to the first embodiment, and only the processing performed by the scattered radiation removal unit 62 is different, so detailed description will be omitted here. FIG. 37 is a diagram showing the functional configuration of the scattered radiation removal unit in the information derivation device according to the seventh embodiment. As shown in FIG. 37, the scattered radiation removal unit 62B in the information derivation device according to the seventh embodiment includes a first derivation unit 97, a second derivation unit 98, and an image generation unit 99.

[0172] In the information derivation device according to the seventh embodiment, the first derivation unit 97 derives a first primary ray distribution and a scattered ray distribution of the radiation transmitted through the subject H using the first and second radiation images G1 and G2, the second derivation unit 98 derives a second primary ray distribution and a scattered ray distribution of the radiation transmitted through the object using the first primary ray distribution and the scattered ray distribution, and the radiation characteristics of the object interposed between the subject H and the radiation detector that detects the radiation image, and the image generation unit 99 performs a scattered ray removal process by deriving a radiation image after transmission through the subject and the object using the second primary ray distribution and the scattered ray distribution. The scattered ray removal process performed by the information derivation device according to the seventh embodiment is described in International Publication No. 2020 / 241664. Hereinafter, the scattered ray removal process from the first radiation image G1 will be described, but the scattered ray removal process can also be performed on the second radiation image G2 in the same manner.

[0173] The first derivation unit 97 uses the first radiographic image G1 to estimate the components (primary ray components and scattered ray components) of the radiation that has passed through the subject H. In the derivation process performed by the first derivation unit 97 (hereinafter referred to as the first derivation process), the "radiation that has passed through the subject H" refers to the radiation that has passed through the subject H and has not yet passed through objects such as the tabletop 4A and the grid 8.

[0174] Moreover, the component of radiation that has passed through the subject H specifically refers to the component of radiation that has passed through the subject H and / or the component of radiation that has been scattered by the subject H. In other words, the component of radiation that has passed through the subject H is the primary ray component after passing through the subject H. The component of radiation that has been scattered by the subject H is the scattered ray component after passing through the subject H. Regarding radiation that is incident on the subject H toward an arbitrary position X, if the subject H is considered to be an operator g1 that generates primary rays and an operator h1 that generates scattered ray components, then as shown in FIG. 38, the primary ray component after passing through the subject H is g1(X), and the scattered ray component after passing through the subject H is h1(X).

[0175] The first derivation unit 97 uses the first radiation image G1 to estimate the primary ray component g1(X) after passing through the subject H, the scattered ray component h1(X) after passing through the subject H, or both. In the seventh embodiment, the first derivation unit 97 uses the first radiation image G1 to estimate each of the primary ray component g1(X) after passing through the subject H and the scattered ray component h1(X) after passing through the subject H.

[0176] When estimating the primary ray component g1(X) after passing through the subject H from the first radiographic image G1, the first derivation unit 97 estimates the scattered ray component h1(X) after passing through the subject H by subtracting the estimated primary ray component g1(X) from the first radiographic image G1. When estimating the scattered ray component h1(X) after passing through the subject H from the first radiographic image G1, the first derivation unit 97 estimates the primary ray component g1(X) after passing through the subject H by subtracting the estimated scattered ray component h1(X) from the first radiographic image G1.

[0177] The first derivation process performed by the first derivation unit 97 can be performed, for example, by estimating the body thickness of the subject H using the first radiographic image G1, and then estimating the components of the radiation that have passed through the subject H using the estimated body thickness of the subject H. In this case, the first derivation unit 97 estimates a primary ray component g1(X) of radiation that has passed through the subject H and a scattered ray component h1(X) of radiation that has been scattered by the subject H for each pixel of the first radiographic image G1 (or for each predetermined section consisting of a plurality of pixels) based on the estimated body thickness of the subject H.

[0178] For example, as shown in FIG. 39, the pixel value V2 when the subject H is present ("subject present") is smaller than the pixel value V1 when the subject H is not present ("no subject"). This is due to absorption by the subject H, etc. Therefore, the difference Δ (=V1-V2) between them is related to the body thickness of the subject H. On the other hand, the pixel value V1 when the subject H is not present can be known from the pixel value of the region (direct region) where the radiation reaches the radiation detector 5 without passing through the subject H, or from a previous experiment (imaging without the subject H being present), etc. Therefore, the first derivation unit 97 can estimate the body thickness of the subject H from the pixel value V2 of the first radiographic image G1 captured with the subject H present.

[0179] Furthermore, both the primary ray component g1(X) and the scattered ray component h1(X) after passing through the subject H are related to the body thickness of the subject H. For example, the thicker the body of the subject H, the smaller the primary ray component g1(X) due to absorption, etc. of the subject H, and the larger the scattered ray component h1(X) relative to the incident radiation. Such properties of the subject H, that is, the amount of radiation that has a specific energy transmitted through and scattered by the subject H, can be determined in advance by experiments, etc., prior to radiography.

[0180] For this reason, the first derivation unit 97 holds, for example, characteristics related to the amount of transmission and the amount of scattering (hereinafter referred to as subject scattering characteristics) in the form of a function, a table, or the like for each subject H or each imaging region of the subject H. Then, by using the energy of the radiation used in imaging and the estimated actual body thickness of the subject H to obtain the amount of transmission and the amount of scattering of the radiation, the primary ray component g1(X) and the scattered ray component h1(X) after transmission through the subject H are estimated.

[0181] The derivation result output by the first derivation unit 97 (hereinafter referred to as the first derivation result) is the primary ray component g1(X) at the position P1 after passing through the subject H, the scattered ray component h1(X) at the position P1 after passing through the subject H, or the intensity distribution f1(X) of the radiation at the position P1 after passing through the subject H. The intensity distribution f1(X) of the radiation at the position P1 is, for example, the sum or weighted sum of the primary ray component g1(X) and the scattered ray component h1(X). In the seventh embodiment, the first derivation unit 97 outputs the intensity distribution f1(X) of the radiation at the position P1 after passing through the subject H as the first derivation result, for example, in the form of an image or in the form of a collection of data from which an image can be constructed. Note that the first derivation unit 97 can also output one of the primary ray component g1(X) or the scattered ray component h1(X) after passing through the subject H as the derivation result.

[0182] The second derivation unit 98 estimates the components of the radiation that have passed through an object, using the derivation results of the first derivation unit 97 and the scattering characteristics of objects such as the tabletop 4A and the grid 8, through which the radiation further passes after passing through the subject H. In the derivation process of the second derivation unit 98 (hereinafter referred to as the second derivation process), "passing through an object" means that the subject H has passed through a certain position and then passed through the object. Therefore, it also includes passing directly through an object without passing through the subject H, depending on the specific shape of the subject H, etc.

[0183] Specifically, the second derivation unit 98 estimates a component of radiation that has passed through the subject H and an object, or a component of radiation that has been scattered by at least one of the subject H and the object. The component of radiation that has passed through the subject H and an object is a primary ray component after passing through the object. The component of radiation that has been scattered by at least one of the subject H and the object is a scattered ray component after passing through the object.

[0184] The scattering characteristic of an object determines the distribution of the amount of radiation transmitted through the object and / or the amount of radiation scattered by the object. In the seventh embodiment, the scattering characteristic f2(X) includes a first characteristic g2(X) that determines the distribution of the amount of radiation transmitted through the object and a second characteristic h2(X) that determines the distribution of the amount of radiation scattered by the object. Specifically, the scattering characteristic f2(X) is the sum or weighted sum of the first characteristic g2(X) and the second characteristic h2(X), e.g., f2(X)=g2(X)+h2(X).

[0185] The first characteristic g2(X) is a function or a table that determines the amount of radiation penetrating the object toward an arbitrary position X and directly without passing through the subject H. The second characteristic h2(X) is a function or a table that determines the amount of radiation penetrating the object toward an arbitrary position X and directly without passing through the subject H. For example, when the object is only the top plate 4A of the imaging table 4, the first characteristic g2(X) determines the distribution of the amount of radiation penetrating the top plate 4A, and the second characteristic h2(X) determines the distribution of the amount of scattered radiation of the top plate 4A. The specific configuration of the object (whether the imaging table 4 is in use or not, etc.) is known before radiography. For this reason, the first characteristic g2(X) and the second characteristic h2(X) can be obtained in advance by experiments, etc., for each specific configuration of the object or each combination of objects. Furthermore, if an object is considered to generate a primary ray component and a scattered ray component from incident radiation, the first characteristic g2(X) is an operator that generates a primary ray component in response to the incident radiation, and the second characteristic h2(X) is an operator that generates a scattered ray component in response to the incident radiation.

[0186] In the seventh embodiment, the second derivation unit 98 holds in advance the first characteristic g2(X) and the second characteristic h2(X), for example, for each specific configuration of the object. As a result, the second derivation unit 98 holds in advance the scattering characteristic f2(X) of the object. However, the second derivation unit 98 can acquire the first characteristic g2(X), the second characteristic h2(X) and / or the scattering characteristic f2(X) as necessary.

[0187] As shown in Fig. 40, the intensity distribution of radiation (X) incident on an object toward an arbitrary point X0 after it has passed through the object can be approximated by a PSF (point spread function) 120. The PSF 120 is, for example, a Gaussian function. Of the radiation (X) incident on an object toward an arbitrary point X0, the component that reaches the arbitrary point X0 and its vicinity is a primary ray component distribution 121, and the portion of the PSF 120 excluding the primary ray component distribution 121 is a scattered ray component distribution 122.

[0188] Furthermore, since the energy of the radiation used for imaging, and the material (density, etc.) and thickness (mass) of the object such as the tabletop 4A are known, the specific shape of the PSF 120, such as the peak height and half-width, are determined in advance. Therefore, for example, the second characteristic h2(X) can be obtained in advance by performing a deconvolution operation of the above-mentioned scattered ray component distribution 122 on a radiation image obtained by imaging without placing the subject H. Also, the first characteristic g2(X) can be obtained in advance by subtracting the second characteristic h2(X) from the radiation image obtained by imaging without placing the subject H, or by performing a deconvolution operation of the primary ray component distribution 121.

[0189] The second derivation unit 98 estimates the components of the radiation transmitted through the object by applying the scattering characteristics of the object to the first derivation result, which is the derivation result of the first derivation unit 97. Specifically, the object is irradiated with radiation having a distribution represented by the first derivation result. For this reason, the second derivation unit 98 estimates the components of the radiation transmitted through the object by taking the first derivation result (f1(X)) as the argument of the scattering characteristics f2(X) of the object. That is, the second derivation unit 98 estimates the components of the radiation transmitted through the object by calculation based on the following formula (25). In the seventh embodiment, the first derivation result is f1(X)=g1(X)+h1(X), so the above formula (25) can be expressed as the following formula (26) and expanded as the following formula (27).

[0190] f2(f1(X))=g2(f1(X))+h2(f1(X)) (25) f2(f1(X))=g2(g1(X)+h1(X))+h2(g1(X)+h1(X)) (26) f2(f1(X))=g2g1(X)+g2h1(X)+h2g1(X)+h2h1(X) (27)

[0191] As shown in FIG. 41, "g2g1(X)" in the first term on the right side of the above formula (27) represents ray Ra1, which is used for imaging and passes through the subject H and the object (tabletop 4A is shown in FIG. 41) to reach pixel P(X0) at any point X0. "g2h1(X)" in the second term on the right side of the formula (27) represents ray Ra2, which is used for imaging and is scattered by scatterer D1 included in the subject H, passes through the object, and reaches pixel P(X0) at any point X0. "h2g1(X)" in the third term on the right side of the formula (27) represents ray Ra3, which is used for imaging and passes through subject H, passes through scatterer D3 included in the object, and reaches pixel P(X0) at any point X0. Additionally, “h2h1(X)” in the fourth term on the right-hand side of equation (27) represents the ray Ra4 used for imaging, which is scattered by the scatterer D2 contained in the subject H, and then further scattered by the scatterer D4 contained in the object, to reach the pixel P(X0) at an arbitrary point X0.

[0192] From the above, the second derivation unit 98 obtains the first term "g2g1(X)" and / or the sum of the second to fourth terms "g2h1(X)+h2g1(X)+h2h1(X)" of equation (27). This is because the first term "g2g1(X)" of equation (27) represents the distribution of the primary ray components after passing through an object, and the sum of the second to fourth terms "g2h1(X)+h2g1(X)+h2h1(X)" represents the distribution of the scattered ray components after passing through an object. In the seventh embodiment, the second derivation unit 98 obtains the distribution g2g1(X) of the primary ray components after passing through an object, and outputs this as a derivation result.

[0193] The image generation unit 99 uses the derivation result of the second derivation unit 98 to generate a processed radiographic image that forms an image of the subject H using the radiation that has passed through the subject H and the object. When the second derivation unit 98 estimates a primary ray component of radiation that has passed through the subject H and the object, the image generation unit 99 generates a processed radiographic image by imaging the second derivation result that is the derivation result of the second derivation unit 98. When the second derivation unit 98 estimates a scattered ray component of radiation that has been scattered by the subject H or the object, the image generation unit 99 generates a processed radiographic image by subtracting the second derivation result that is the derivation result of the second derivation unit 98 from the first radiographic image G1.

[0194] In the seventh embodiment, the second derivation unit 98 outputs the distribution of the primary ray components after passing through an object, and the image generation unit 99 generates a processed radiographic image by imaging this. Therefore, the distribution g2g1(X) of the primary ray components output by the second derivation unit 98 is a processed radiographic image from which scattered ray components have essentially been removed. Note that the image generation unit 99 can perform various types of image processing (for example, contrast adjustment processing or structure emphasis processing) on ​​the generated processed radiographic image as necessary.

[0195] The scattered radiation removal processing according to the sixth and seventh embodiments may be performed in the first to fifth embodiments. In this case, the information derivation devices 50A to 50D have the scattered radiation removal unit 62A or the scattered radiation removal unit 62B instead of the scattered radiation removal unit 62.

[0196] Next, an eighth embodiment of the present disclosure will be described. FIG. 42 is a diagram showing a functional configuration of an information derivation device according to the eighth embodiment. In FIG. 42, the same components as those in FIG. 7 are given the same reference numerals, and detailed description will be omitted. In the seventh embodiment of the present disclosure, graininess suppression processing is performed on the first and second radiographic images G1, G2, and energy subtraction processing is performed using the first and second radiographic images G1, G2 on which graininess suppression processing has been performed. For this reason, as shown in FIG. 42, an information derivation device 50E according to the eighth embodiment further includes a processing content derivation unit 86 and a graininess suppression processing unit 87 in addition to the information derivation device 50 according to the first embodiment.

[0197] The processing content derivation section 86 derives the processing content of the first graininess suppression processing for the first radiographic image G1, and derives the processing content of the second graininess suppression processing for the second radiographic image G2 based on the processing content of the first graininess suppression processing.

[0198] In this embodiment, the subject H is imaged by the one-shot method to obtain the first and second radiographic images G1, G2. In the case of the one-shot method, the radiation transmitted through the subject H is irradiated onto two radiation detectors 5, 6 that are stacked with a radiation energy conversion filter 7 interposed therebetween. Therefore, the second radiation detector 6 located away from the radiation source 3 is irradiated with a smaller dose of radiation than the first radiation detector 5 located closer to the radiation source 3. As a result, the second radiographic image G2 has more quantum noise of radiation and a lower S / N ratio than the first radiographic image G1. For this reason, it is necessary to perform graininess suppression processing on the second radiographic image G2 in particular to suppress graininess caused by quantum noise.

[0199] The processing content derivation unit 86 derives the processing content of the first graininess suppression processing for the first radiographic image G1, which has a higher S / N ratio, out of the first radiographic image G1 and the second radiographic image G1. The processing content derivation unit 86 then derives the processing content of the second graininess suppression processing for the second radiographic image G2 based on the processing content of the first graininess suppression processing. The derivation of the processing content will be described below.

[0200] The graininess suppression process may be a filtering process using a smoothing filter such as a Gaussian filter having a predetermined size such as 3×3 or 5×5 centered on the pixel of interest. However, when a Gaussian filter is used, there is a possibility that the edges of structures included in the first and second radiographic images G1 and G2 may become blurred. For this reason, in the eighth embodiment, the graininess suppression process is performed using an edge-preserving smoothing filter that suppresses graininess while preventing edge blurring. As the edge-preserving smoothing filter, a bilateral filter is used that uses a normal distribution weighting for pixels adjacent to the pixel of interest such that the weighting decreases the farther away from the pixel of interest and the weighting decreases the greater the difference in pixel value between the pixel of interest and the pixel of interest.

[0201] Fig. 43 is a diagram showing an example of a bilateral filter for a first radiographic image G1. Fig. 43 shows two local regions A1 of 5 x 5 pixels near an edge included in the first radiographic image G1. The two local regions A1 shown in Fig. 43 are the same, but the positions of the target pixels are different. In the local region A1 on the left side of Fig. 43, a low density pixel adjacent to the edge boundary is the target pixel P11. As described above, the bilateral filter weights pixels adjacent to the target pixel according to a normal distribution such that the weighting decreases as the pixel moves away from the target pixel and the weighting decreases as the difference in pixel value between the pixel and the target pixel increases.

[0202] For this reason, the processing content derivation unit 86 determines the filter size of the bilateral filter based on the difference between the pixel value of the pixel of interest and the pixel values ​​of the pixels surrounding the pixel of interest. For example, the smaller the difference between the pixel value of the pixel of interest and the pixel values ​​of the pixels surrounding the pixel of interest, the larger the filter size is made. In addition, the processing content derivation unit 86 determines the weight of the bilateral filter based on the difference between the pixel value of the pixel of interest and the pixel values ​​of the pixels surrounding the pixel of interest. For example, the smaller the difference between the pixel value of the pixel of interest and the pixel values ​​of the pixels surrounding the pixel of interest, the larger the weight for the pixels close to the pixel of interest is made than the weight for the pixels farther from the pixel of interest.

[0203] As a result, for the pixel of interest P11 in the local area A1 on the left side of Figure 43, the processing content derivation unit 86 derives a 3x3 bilateral filter F11 having weights as shown on the left side of Figure 43 as the processing content of the first graininess reduction processing.

[0204] In the local area A1 on the right side of Fig. 43, a high-density pixel adjacent to the edge boundary is the pixel of interest P12. For this reason, the processing content derivation unit 86 derives a 3 × 3 bilateral filter F12 having weights as shown on the right side of Fig. 43 as the processing content of the first graininess reduction processing for the pixel of interest P12 in the local area A1 on the right side of Fig. 43.

[0205] FIG. 44 is a diagram showing a local area A2 of a second radiographic image corresponding to the local area A1 of the first radiographic image shown in FIG. 43. FIG. 44 shows a local area A2 of 5×5 pixels near an edge included in the second radiographic image G2. The local area A2 is an area at the same position as the local area A1 in the first radiographic image G1 shown in FIG. 43. The two local areas A2 shown in FIG. 44 are the same, but the positions of the target pixels are different. In the local area A2 on the left side of FIG. 44, the pixel corresponding to the target pixel P11 shown in the local area A1 on the left side of FIG. 43 is the target pixel P21. In the local area A2 on the right side of FIG. 44, the pixel corresponding to the target pixel P12 shown in the local area A1 on the right side of FIG. 43 is the target pixel P22.

[0206] Here, the second radiation detector 6 from which the second radiation image G2 is acquired is irradiated with a lower dose of radiation than the first radiation detector 5 from which the first radiation image G1 is acquired. For this reason, the second radiation image G2 has more quantum noise of radiation and poor graininess compared to the first radiation image G1, and edges do not appear clearly. In addition, due to the influence of quantum noise, low density pixels are included in high density areas near edge boundaries, and high density pixels are included in low density areas. For this reason, it is not possible to appropriately determine a bilateral filter that suppresses graininess while preserving edges from the second radiation image G2, as in the case of the first radiation image G1.

[0207] In this case, it is conceivable to use a smoothing filter such as a Gaussian filter, but such a smoothing filter cannot suppress noise while preserving edges at the same time, and as a result, the edges are buried in noise, making it impossible to restore the edges of the structures included in the second radiographic image G2.

[0208] For this reason, in the eighth embodiment, the processing content derivation unit 86 derives the processing content of the second graininess suppression processing for the second radiographic image based on the processing content of the first graininess suppression processing for the first radiographic image G1. That is, the processing content derivation unit 86 derives the processing content of the second graininess suppression processing so that the processing content of the second graininess suppression processing performed for each pixel of the second radiographic image G2 is the same as the processing content of the first graininess suppression processing performed for pixels of the first radiographic image G1 corresponding to each pixel of the second radiographic image G2. Specifically, the processing content derivation unit 86 derives a bilateral filter having the same size and the same weight as the bilateral filter determined for each pixel of the first radiographic image G1 as the processing content of the second graininess suppression processing for the second radiographic image G2.

[0209] Fig. 45 is a diagram showing an example of a bilateral filter for the second radiographic image G2. Fig. 45 shows a local area A2 of 5 × 5 pixels near an edge included in the second radiographic image G2, similar to Fig. 43. As shown in Fig. 45, for a pixel of interest P21 in the local area A2 of the second radiographic image G2, the processing content derivation unit 86 derives, as processing content of the second graininess suppression processing, a bilateral filter F21 having the same size and the same weight as the bilateral filter F11 derived for the pixel of interest P11 in the local area A1 of the first radiographic image G1.

[0210] Furthermore, for the pixel of interest P22 in the local region A2 of the second radiographic image G2, the processing content derivation unit 86 derives, as the processing content of the second graininess suppression processing, a bilateral filter F22 having the same size and weight as the bilateral filter F12 derived for the pixel of interest P12 in the local region A1 of the first radiographic image G1.

[0211] In addition, when deriving the processing content of the second graininess reduction processing, it is necessary to associate pixel positions between the first radiographic image G1 and the second radiographic image G2. For this reason, it is preferable to align the positions of the first radiographic image G1 and the second radiographic image G2.

[0212] The graininess suppression processing unit 87 performs graininess suppression processing on the first radiographic image G1 and the second radiographic image G2. That is, the graininess suppression processing is performed on the first radiographic image G1 and the second radiographic image G2 based on the processing content derived by the processing content derivation unit 86. Specifically, the first radiographic image G1 is subjected to filtering processing using a bilateral filter derived for the first radiographic image G1. Furthermore, the second radiographic image G2 is subjected to filtering processing using a bilateral filter derived based on the first radiographic image G1.

[0213] In the eighth embodiment, the subtraction unit 63 derives a bone image Gb and a soft tissue image Gs by using the first and second radiographic images G1 and G2 that have been subjected to the graininess suppression process, according to the above formulas (1) and (2).

[0214] Next, a ninth embodiment of the present disclosure will be described. FIG. 46 is a diagram showing a functional configuration of an information derivation device according to the ninth embodiment. In FIG. 46, the same components as those in FIG. 42 are given the same reference numerals, and detailed description will be omitted. The information derivation device 50F according to the ninth embodiment of the present disclosure is different from the information derivation device 50E according to the eighth embodiment in that the information derivation device 50F further includes a map derivation unit 88 that derives a physical quantity map of the subject H based on at least one of the first radiographic image G1 and the second radiographic image G2, and the processing content derivation unit 86 derives the processing content of the second graininess suppression processing for the second radiographic image G2 based on the physical quantity map.

[0215] The map derivation unit 88 derives a physical quantity map for the subject H. Examples of physical quantities include the body thickness and bone density of the subject H. The body thickness can be the body thickness distribution when the scattered radiation removal unit 62 satisfies the termination condition. The bone density can be the bone density derived from the bone image Gb by a predetermined method.

[0216] In the first and second radiographic images G1 and G2 to be subjected to the graininess suppression process, the contrast of the structures contained in the images varies depending on the imaging conditions. Therefore, when performing edge-preserving smoothing process using a bilateral filter, it is necessary to control the strength of the edges to be preserved depending on the imaging conditions. On the other hand, in the body thickness map representing the body thickness of the subject H, the contrast of the structures contained in the map is represented by thickness (mm) that does not depend on the imaging conditions.

[0217] Therefore, in the ninth embodiment, the processing content derivation unit 86 derives the processing content of the first graininess suppression processing for the first radiographic image G1 based on the physical quantity map. FIG. 47 is a diagram showing an example of a bilateral filter for the physical quantity map. Note that FIG. 47 shows two local regions A3 of 5×5 pixels near the edge included in the physical quantity map. The two local regions A3 shown in FIG. 47 are the same, but the positions of the target pixels are different. In the local region A3 on the left side of FIG. 47, a high density pixel on the edge is the target pixel P31. Therefore, for the target pixel P31 of the local region A3 on the left side of FIG. 47, a 3×3 bilateral filter F31 having a weight as shown on the left side of FIG. 47 is derived as the processing content of the first graininess suppression processing.

[0218] In the local area A3 on the right side of Fig. 47, a low density pixel adjacent to the edge boundary is the pixel of interest P32. Therefore, for the pixel of interest P32 in the local area A3 on the right side of Fig. 47, a 3 × 3 bilateral filter F32 having weights as shown on the right side of Fig. 47 is derived as the processing content of the first graininess reduction processing.

[0219] In the ninth embodiment, the processing content derivation unit 86 also derives a bilateral filter identical to the bilateral filter determined for each pixel of the first radiographic image G1 as the processing content of the second graininess suppression processing for the second radiographic image G2. That is, in the second radiographic image G2, in a local region corresponding to the local region A3 of the first radiographic image G1, a bilateral filter having the same size and the same weight as the bilateral filter F31 is derived as the processing content of the second graininess suppression processing for a pixel corresponding to the pixel of interest P31 of the local region A3. Also, in the second radiographic image G2, a bilateral filter having the same size and the same weight as the bilateral filter F32 is derived as the processing content of the second graininess suppression processing for a pixel corresponding to the pixel of interest P32 of the local region A3 of the first radiographic image G1.

[0220] In the ninth embodiment, the graininess suppression processing unit 87 performs graininess suppression processing on the first radiographic image G1 and the second radiographic image G2 based on the processing content derived by the processing content derivation unit 86. The subtraction unit 63 derives a bone image Gb and a soft tissue image Gs according to the above formulas (1) and (2) using the first and second radiographic images G1 and G2 that have been subjected to the graininess suppression processing.

[0221] The graininess suppression processing performed in the eighth and ninth embodiments may be performed in the first to seventh embodiments.

[0222] Next, an information derivation device according to a tenth embodiment of the present disclosure will be described. The configuration of the information derivation device according to the tenth embodiment is the same as that of the information derivation device according to the first embodiment, and only the process performed by the scattered radiation removal unit 62 is different, so detailed description will be omitted here. In the tenth embodiment, the first and second radiographic images G1 and G2 acquired by photographing the subject H using the imaging device 1A shown in FIG. 29 as in the sixth embodiment are used as the teacher data 40.

[0223] Fig. 48 is a diagram showing the functional configuration of a scattered radiation removal unit in the information derivation device according to the tenth embodiment. As shown in Fig. 48, the scattered radiation removal unit 62C in the information derivation device according to the tenth embodiment includes an area detection unit 150, a scattered radiation image derivation unit 151, a pixel value calculation unit 152, a pixel value replacement unit 153, a boundary position adjustment unit 154, and a calculation unit 155.

[0224] The area detection section 150 obtains area detection images by detecting, in the first and second radiation images G1, G2, subject areas where radiation has passed through the subject H to reach the first and second radiation detectors 5, 6, and direct radiation areas where radiation has not passed through the subject H but has passed only through the tabletop 4A to directly reach the first and second radiation detectors 5, 6.

[0225] Specifically, as shown in Fig. 49, the region detection section 150 detects a region in the first radiographic image G1 where the pixel value is less than the region threshold as the subject region 160, and detects a region in the first radiographic image G1 where the pixel value is equal to or greater than the region threshold as the direct radiation region 161. For example, when the region detection image is a binarized image in which the subject region 160 is set to "0" and the direct radiation region 161 is set to "1", a distribution as shown in Fig. 50 is obtained on a specific line 162 of the first radiographic image G1. Note that the region detection section 150 may use a trained neural network that has been trained to detect the subject region and the direct radiation region, instead of using threshold processing for region detection.

[0226] The region threshold value is determined for each imaging condition, and it is preferable to use a region threshold value corresponding to the imaging condition when imaging the subject H. The region detection process may be performed on both the first radiographic image G1 and the second radiographic image G2, or on either one of them.

[0227] The scattered radiation image acquiring section 151 acquires a scattered radiation image relating to the scattered radiation components based on the region detection image and the scattered radiation spread information relating to the spread of the scattered radiation. The calculation section 155 acquires a radiographic image from which the scattered radiation components have been removed by subtracting the scattered radiation image from the first and second radiographic images G1, G2. Note that the scattered radiation spread information used by the scattered radiation image acquiring section 151 may be, for example, that shown in FIG. 40.

[0228] The scattered radiation image derivation section 151 then obtains a scattered radiation image containing scattered radiation components by convolving the region detection image with the PSF, ie, "region detection image*PSF (* is a convolution operator)." The scattered radiation image derivation section 151 derives scattered radiation images for each of the first radiation image G1 and the second radiation image G2.

[0229] In the scattered radiation image, the pixel values ​​of the scattered radiation components in the direct radiation region may be directly used as the pixel values ​​of the direct radiation region 161 of the first and second radiation images G1, G2, but the pixel values ​​of the direct radiation region 161 are often saturated (exceeding the pixel values ​​that the imaging sensors of the first and second radiation detectors 5, 6 can receive).For this reason, it is preferable to theoretically calculate, as the pixel values ​​of the direct radiation region of the scattered radiation image, non-saturated scattered radiation pixel values ​​that are not saturated in pixel value and that take into account the influence of scattered rays, and to replace the theoretically calculated non-saturated scattered radiation pixel values ​​with the pixel values ​​of the direct radiation region of the scattered radiation image.

[0230] The pixel value calculation unit 152 calculates a non-saturated scattered radiation pixel value corresponding to the imaging dose by referring to a relationship for non-saturated scattered radiation pixel values ​​indicating the relationship between the radiation dose and non-saturated scattered radiation pixel values, which are non-saturated pixel values ​​obtained when there is no pixel value saturation and the effect of scattered radiation is taken into consideration.Then, the pixel value replacement unit 153 replaces the pixel values ​​of the direct radiation region with the non-saturated scattered radiation pixel values.Specifically, as shown in FIG. 51, in the scattered radiation image derived by the scattered radiation image derivation unit 151, the pixel values ​​of the direct radiation region 161 are replaced with the non-saturated scattered radiation pixel values ​​calculated by the pixel value calculation unit 152.

[0231] The non-saturated scattered radiation pixel value relations are determined in advance in a table (not shown) for the direct radiation region for each imaging condition and stored in the storage 53. The pixel value calculation unit 152 refers to the table for the direct radiation region and uses any one of the non-saturated scattered radiation pixel value relations corresponding to the imaging conditions when imaging the subject H, or uses a combination of multiple non-saturated scattered radiation pixel value relations that satisfy the imaging conditions when imaging the subject H. In addition, it is preferable to calculate the non-saturated scattered radiation pixel value by, for example, "non-saturated scattered radiation pixel value = non-saturated pixel value obtained when there is no saturation of pixel value × scattered radiation content rate of the direct radiation region".

[0232] As described above, the region detection section 150 distinguishes between the subject region 160 and the direct radiation region 161 using the region threshold value, but the direct radiation region 161 detected using the region threshold value may also include a region that has passed through a soft tissue region close to the skin of the subject H. If scattered radiation is assumed to be emitted from such a soft tissue region, subtracting the scattered radiation image from the first and second radiation images G1, G2 will result in excessive removal of scattered radiation components. For this reason, the boundary position between the subject region 160 and the direct radiation region 161 can be adjusted by a specific width, and a scattered radiation image is obtained based on the region detection image after the boundary position has been adjusted and the scattered radiation spread information.

[0233] The boundary position adjustment unit 154 adjusts the boundary position between the subject region 160 and the direct radiation region 161 by a specific width. The adjustment of the boundary position is preferably performed when the pixel value of the direct radiation region exceeds a pixel value threshold value, or when the radiation dose exceeds a dose threshold value. In addition, the specific width is preferably determined based on at least one of the part of the subject H, the imaging method, and the imaging conditions of the subject H. In addition, the specific width is preferably determined based on pixel values ​​in the radiation image near the boundary position.

[0234] Specifically, when the specific width is determined based on the shape of the subject H among the parts of the subject H, that is, the body thickness distribution of the subject H, the specific width is a portion Px of the direct radiation region 161 detected by the region threshold value where the body thickness exceeds 0, as shown in the body thickness distribution of the specific line 162 portion in Fig. 52. By adjusting the boundary position BP by the specific width corresponding to this portion Px, the subject region 160 expands by the specific width, while the direct radiation region 161 narrows by the specific width. In addition, when the specific width is determined based on pixel values ​​in the vicinity of the boundary position, it is preferable to determine the specific width based on the inter-pixel distance between the pixel at the boundary position and a neighboring pixel that is within a specific range from the boundary position and has a pixel value in a specific range.

[0235] In the tenth embodiment, a subtraction unit 63 derives a bone image Gb and a soft tissue image Gs using the processed first and second radiation images G1, G2.

[0236] The scattered radiation removal process according to the tenth embodiment may be performed in the first to fifth, eighth and ninth embodiments.

[0237] In addition, in each of the above-mentioned embodiments, the estimation results of the bone image Gb and the soft tissue image are derived from the simple radiographic image G0, but the present invention is not limited to this. For example, the estimation results of the muscle image and the fat image may be derived. This will be described below as the eleventh embodiment.

[0238] FIG. 53 is a diagram showing a functional configuration of an information derivation device according to an eleventh embodiment. In FIG. 53, the same configuration as in FIG. 7 is given the same reference numerals, and detailed description is omitted. In the eleventh embodiment of the present disclosure, instead of deriving a bone image Gb and a soft tissue image Gs as the correct answer data 42, a muscle image Gm and a fat image Gf are derived. For this reason, as shown in FIG. 53, an information derivation device 50G according to the eleventh embodiment further includes a muscle image derivation unit 110 and a fat image derivation unit 111 in addition to the information derivation device 50 according to the first embodiment.

[0239] The muscle image derivation unit 110 derives the muscle image Gm using the soft tissue image Gs. For this purpose, the muscle image derivation unit 110 derives muscle mass based on the pixel value for each pixel in the soft tissue region in the soft tissue image Gs. The soft tissue includes muscle tissue, fat tissue, blood, and water. In the muscle image derivation unit 110 of the eleventh embodiment, tissue other than fat tissue in the soft tissue is considered as muscle tissue. That is, in the muscle image derivation unit 110 of the eleventh embodiment, non-fatty tissue including blood and water is treated as muscle tissue.

[0240] The muscle image derivation unit 110 separates muscle and fat from the soft tissue image Gs by utilizing the difference in energy characteristics between muscle tissue and fat tissue. As shown in FIG. 54, the dose of radiation after passing through the subject H, which is a human body, is lower than that of radiation before entering the subject H. In addition, since muscle tissue and fat tissue have different absorbed energy and different attenuation coefficients, the energy spectrum of the radiation after passing through the subject H is different from that of the radiation after passing through the muscle tissue. As shown in FIG. 54, the energy spectrum of the radiation that passes through the subject H and is irradiated to each of the first radiation detector 5 and the second radiation detector 6 depends on the body composition of the subject H, specifically, the ratio of muscle tissue to fat tissue. Since fat tissue is more easily permeable to radiation than muscle tissue, the dose of radiation after passing through the human body is lower when the ratio of muscle tissue is higher than that of fat tissue.

[0241] For this reason, the muscle image derivation unit 110 separates muscle and fat from the soft tissue image Gs by utilizing the difference in the energy characteristics between the muscle tissue and the fat tissue described above. That is, the muscle image derivation unit 110 generates a muscle image and a fat image from the soft tissue image Gs.

[0242] Note that the specific method by which the muscle image derivation unit 110 separates muscle and fat from the soft tissue image Gs is not limited. As an example, the muscle image derivation unit 110 of the eleventh embodiment derives a muscle image from the soft tissue image Gs by the following formula (28) and formula (29). Specifically, first, the muscle image derivation unit 110 derives the muscle ratio rm(x, y) at each pixel position (x, y) in the soft tissue image Gs by formula (28). Note that μm in formula (28) is a weighting coefficient according to the attenuation coefficient of muscle tissue, and μf is a weighting coefficient according to the attenuation coefficient of fat tissue. Also, Δ(x, y) represents the concentration difference distribution. The concentration difference distribution is the distribution on the image of the concentration change seen from the concentration obtained by the radiation reaching the first radiation detector 5 and the second radiation detector 6 without penetrating the subject H. The distribution of density changes in the image is calculated by subtracting the density of each pixel in the area of ​​the subject H from the density in the direct-passage area obtained by irradiating the first radiation detector 5 and the second radiation detector 6 with radiation in the soft tissue image Gs directly. rm(x,y)={μf-Δ(x,y) / T(x,y)} / (μf-μm) (28)

[0243] Furthermore, the muscle image derivation unit 110 generates a muscle image Gm from the soft part image Gs by the following formula (29): In addition, in formula (29), x and y are the coordinates of each pixel of the muscle image Gm. Gm(x,y)=rm(x,y)×Gs(x,y) (29)

[0244] The fat image derivation unit 111 generates a fat image Gf from the soft tissue image Gs using the muscle ratio rm according to the following formula (30). Gf(x,y)=(1-rm(x,y))×Gs(x,y) (30)

[0245] In the eleventh embodiment, the muscle image Gm and fat image Gf derived by the information derivation device 50G are used as correct answer data of the teacher data. Fig. 55 is a diagram showing the teacher data derived in the eleventh embodiment. As shown in Fig. 55, the teacher data 40A is composed of learning data 41 including the first and second radiation images G1 and G2, and correct answer data 42A including the muscle image Gm and fat image Gf.

[0246] By training a neural network using the training data 40A shown in FIG. 55, it is possible to construct a trained neural network 23A that outputs a muscle image Gm and a fat image Gf as estimation results when a plain radiation image G0 is input.

[0247] In the first to tenth embodiments, the estimation results of the bone image Gb and the soft tissue image Gs are derived, but the present invention is not limited to this. The estimation result of either the bone image Gb or the soft tissue image Gs may be derived. In this case, the trained neural network 23A may be constructed by performing training using teacher data in which the correct answer data is either the bone image Gb or the soft tissue image Gs.

[0248] In the eleventh embodiment, the estimation results of the muscle image Gm and the fat image Gf are derived, but the present invention is not limited to this. The estimation results of either the muscle image Gm or the fat image Gf may be derived. In this case, the trained neural network 23A may be constructed by performing training using teacher data in which the correct answer data is either the muscle image Gm or the fat image Gf.

[0249] In each of the above embodiments, the trained neural network 23A is constructed by training the neural network in the estimation device 10, but the present invention is not limited to this. The trained neural network 23A constructed in a device other than the estimation device 10 may be used in the estimation unit 23 of the estimation device 10 in this embodiment.

[0250] In addition, in each of the above-described embodiments, the first and second radiographic images G1, G2 are acquired by a one-shot method when performing the energy subtraction process, but the present invention is not limited to this. The first and second radiographic images G1, G2 may be acquired by a so-called two-shot method in which imaging is performed twice using only one radiation detector. In the case of the two-shot method, the position of the subject H contained in the first radiographic image G1 and the second radiographic image G2 may be shifted due to the body movement of the subject H. For this reason, it is preferable to perform the process of this embodiment after aligning the position of the subject in the first radiographic image G1 and the second radiographic image G2.

[0251] The alignment process may be performed using, for example, the method described in JP 2011-255060 A. The method described in JP 2011-255060 A generates a plurality of first band images and a plurality of second band images representing structures having different frequency bands for each of the first and second radiographic images G1 and G2, obtains the amount of positional deviation between corresponding positions in the first band image and the second band image of the corresponding frequency band, and aligns the first radiographic image G1 and the second radiographic image G2 based on the amount of positional deviation.

[0252] In the above-described embodiments, the bone image Gb and soft tissue image Gs as the correct answer data of the teacher data, and the bone image Gb and soft tissue image Gs from the simple radiographic image G0 are estimated using radiographic images acquired in a system that captures the image of the subject H using the first and second radiation detectors 5 and 6. However, the technology of the present disclosure can be applied to a case where the first and second radiographic images G1 and G2 are acquired using stimulable phosphor sheets instead of the radiation detectors. In this case, two stimulable phosphor sheets are stacked and irradiated with radiation that has passed through the subject H, and the radiographic image information of the subject H is accumulatively recorded on each stimulable phosphor sheet, and the radiographic image information is photoelectrically read from each stimulable phosphor sheet to acquire the first and second radiographic images G1 and G2. The two-shot method may also be used when the first and second radiographic images G1 and G2 are acquired using stimulable phosphor sheets.

[0253] Moreover, the radiation in the above embodiment is not particularly limited, and in addition to X-rays, α-rays, γ-rays, etc. can be used.

[0254] In the above embodiment, for example, the image acquisition unit 21, information acquisition unit 22, estimation unit 23, learning unit 24, and display control unit 25 of the estimation device 10, and the image acquisition unit 61, scattered radiation removal unit 62, and subtraction unit 63 of the information derivation device 50, etc. may have a hardware structure of a processing unit that performs various processes, such as the image acquisition unit 21, information acquisition unit 22, estimation unit 23, learning unit 24, and display control unit 25 of the estimation device 10, and the image acquisition unit 61, scattered radiation removal unit 62, and subtraction unit 63 of the information derivation device 50, etc. As described above, the above various processors include a CPU, which is a general-purpose processor that executes software (programs) and functions as various processing units, as well as a programmable logic device (PLD), which is a processor whose circuit configuration can be changed after manufacture, such as an FPGA (Field Programmable Gate Array), a dedicated electric circuit, etc., which is a processor having a circuit configuration designed specifically to perform specific processes, such as an ASIC (Application Specific Integrated Circuit), etc.

[0255] A single processing unit may be configured with one of these various processors, or may be configured with a combination of two or more processors of the same or different types (e.g., a combination of multiple FPGAs or a combination of a CPU and an FPGA). Also, multiple processing units may be configured with a single processor.

[0256] As an example of configuring multiple processing units with one processor, first, there is a form in which one processor is configured with a combination of one or more CPUs and software, as typified by computers such as client and server, and this processor functions as multiple processing units. Second, there is a form in which a processor is used that realizes the functions of the entire system including multiple processing units with one IC (Integrated Circuit) chip, as typified by System On Chip (SoC), etc. In this way, the various processing units are configured using one or more of the above various processors as a hardware structure.

[0257] Furthermore, more specifically, the hardware structure of these various processors can be an electric circuit that combines circuit elements such as semiconductor elements. [Explanation of symbols]

[0258] 1. 1A Imaging Device 3 Radiation source 4. Shooting stand 4A Tabletop 4B Mounting part 5, 6 Radiation detectors 7 Radiation Energy Conversion Filter 8 Anti-scatter grid 9. Image Storage System 10 Estimation device 11, 51 CPU 12 Estimation Processing Program 12B Study Program 13, 53 Storage 14,54 Display 15, 55 Input Devices 16,56 memory 17, 57 Network I / F 18, 58 Bus 21 Image acquisition unit 22 Information Acquisition Department 23 Estimation part 23A Trained Neural Networks 24 Learning Department 25 Display control section 30 Neural Networks 31 Input Layer 32 Middle Class 33 Output layer 35 Convolutional Layer 36 Pooling Layer 37 Fully connected layer 40, 40A Teacher Data 41 Training Data 42, 42A Correct data 47 Output Data 48 Parameters 50,50A, 50B, 50C, 50D, 50E, 50G information derivation device 52 Information derivation program 61 Image acquisition unit 62, 62A, 62B Scattered radiation removal section 63 Subtraction Section 65 Structure recognition part 66 Weighting coefficient derivation part 67 Initial weighting coefficient setting unit 68 Weighting coefficient derivation part 70 display screen 71 First image display area 72 Second image display area 81 Initial value derivation part 82 Attenuation coefficient derivation part 83 Weighting coefficient derivation part 84 Composition ratio derivation part 85 Attenuation coefficient setting section 86 Processing content derivation part 87 Grain suppression processing unit 88 Map Derivation Part 91 Shooting condition acquisition unit 92 Body thickness derivation part 93 Characteristics Acquisition Unit 94 Line distribution derivation part 95 Arithmetic section 97 First derivation part 98 Second derivation part 99 Image Generation Unit 101 Phantom 102 Lead plate 103 Air Layer 110 Muscle image extraction unit 111 Fat image extraction unit 120 PSF 121 Distribution of Primary Radiation Components 122 Distribution of scattered radiation components 150 Area detection unit 151 Scattered Radiation Image Extraction Unit 152 Pixel value calculation unit 153 Pixel value replacement unit 154 Boundary position adjustment section 155 Arithmetic section 160 Subject area 161 Direct radiation area 162 Specific Line A1~A3 Local area BP boundary position D1~D4 Scatterers F11, F12, F21, F22, F31, F32 Bilateral Filters G0 Plain radiographic image G1, G2 radiological images Gb Bone image Gf Fat Images Gm muscle images Gs Soft tissue image H Subject K1~K3 Radiation images LUT1~LUT9 tables P11, P12, P21, P22, P31, P32 Areas of interest Ra radiation

Claims

1. At least one processor; The processor, The present invention functions as a trained neural network that derives an estimation result of at least one enhanced image in which a specific composition of a subject is enhanced from a simple two-dimensional image acquired by simply photographing the subject including a plurality of compositions, the trained neural network is constructed by using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions and a training emphasized image derived from the two radiation images in which a specific composition of the subject is emphasized, inputting the two radiation images into the neural network to output an emphasized image from the neural network, and training the neural network so that a difference between the emphasized image output from the neural network and the training emphasized image becomes small.

2. 2. The estimation device according to claim 1, wherein the learning weighted image is derived by energy subtraction processing that performs weighted subtraction on the two radiation images.

3. A system comprising at least one processor, The processor, The present invention functions as a trained neural network that derives an estimation result of at least one enhanced image in which a specific composition of a subject is enhanced from a simple two-dimensional image acquired by simply photographing the subject including a plurality of compositions, the trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training emphasized image in which a specific composition of the subject is emphasized and which is derived by an energy subtraction process of performing weighted subtraction between the two radiation images; the emphasized image is at least one of a bone image in which a bone portion of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, The learning emphasis image is recognizing the bones and the soft tissue of the subject using at least one of the two radiation images; deriving attenuation coefficients for the bones and the soft tissue using the recognition results of the bones and the soft tissue and the two radiological images; The attenuation coefficient is derived by performing the energy subtraction process.

4. A method for detecting a signal from a computer comprising: The processor, The present invention functions as a trained neural network that derives an estimation result of at least one enhanced image in which a specific composition of a subject is enhanced from a simple two-dimensional image acquired by simply photographing the subject including a plurality of compositions, the trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training emphasized image in which a specific composition of the subject is emphasized and which is derived by an energy subtraction process of performing weighted subtraction between the two radiation images; the emphasized image is a bone image in which a bone portion of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, The learning emphasis image is deriving a new weighting coefficient to be used in the weighted subtraction based on pixel values ​​of the bone portion included in the bone image and pixel values ​​of the soft tissue included in the soft tissue image; performing the weighted subtraction on the two radiographic images using the new weighting coefficients to derive a new bone image and a new soft tissue image; An estimation device that derives a further new weighting coefficient based on the new bone image, a further new weighting coefficient based on the new soft tissue image, and a further new bone image and a further new soft tissue image based on the further new weighting coefficient.

5. A method for detecting a signal from a computer comprising: The processor, The present invention functions as a trained neural network that derives an estimation result of at least one enhanced image in which a specific composition of a subject is enhanced from a simple two-dimensional image acquired by simply photographing the subject including a plurality of compositions, the trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training emphasized image in which a specific composition of the subject is emphasized and which is derived by an energy subtraction process of performing weighted subtraction between the two radiation images; the emphasized image is a bone image in which a bone portion of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, The learning emphasis image is deriving a difference between the value of (attenuation coefficient of the soft part)×(thickness of the soft part)+(attenuation coefficient of the bone part)×(thickness of the bone part) for each different energy distribution and each pixel value of the radiographic image while changing the attenuation coefficient for each different energy distribution for the soft part, the thickness of the soft part, the attenuation coefficient for each different energy distribution for the bone part, and the thickness of the bone part from initial values; Deriving an attenuation coefficient of the soft tissue and an attenuation coefficient of the bone for each of the different energy distributions, in which the difference is minimized or the difference is less than a predetermined threshold value; The estimation device derives the attenuation coefficient of the soft tissue and the attenuation coefficient of the bone by performing the energy subtraction process using a weighting coefficient derived based on the attenuation coefficient of the soft tissue and the attenuation coefficient of the bone.

6. A method for detecting a signal from a computer comprising: The processor, The present invention functions as a trained neural network that derives an estimation result of at least one enhanced image in which a specific composition of a subject is enhanced from a simple two-dimensional image acquired by simply photographing the subject including a plurality of compositions, the trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training emphasized image in which a specific composition of the subject is emphasized and which is derived by an energy subtraction process of performing weighted subtraction between the two radiation images; The learning emphasis image is Deriving composition ratios of a plurality of compositions contained in the soft part of the subject; deriving an attenuation coefficient for each energy distribution that differs for each pixel of the two radiological images according to the composition ratio; The estimation device derives the energy subtraction process using weighting coefficients derived based on the derived attenuation coefficients.

7. The composition ratio is deriving a body thickness of the subject as a first body thickness and a second body thickness for each pixel of each of the two radiographic images; The estimation device according to claim 6 , wherein the thickness is derived for each pixel of the radiation image based on the first body thickness and the second body thickness.

8. The composition ratio is Deriving the first body thickness and the second body thickness based on the attenuation coefficient for each of the different energy distributions for each of the plurality of compositions; The estimation device according to claim 7, wherein the first body thickness and the second body thickness are derived while changing the thickness of the composition and the attenuation coefficient for each composition, and the first body thickness and the second body thickness are derived based on the thickness of the composition at which the difference between the first body thickness and the second body thickness is equal to or less than a predetermined threshold value.

9. The learning emphasis image is 9. The estimation device according to claim 2, wherein a scattered radiation removal process is performed to remove from the two radiographic images a scattered radiation component that is scattered by the subject out of the radiation irradiated to the subject, and the energy subtraction process is performed on the two radiographic images from which the scattered radiation component has been removed.

10. the scattered radiation removal process acquires radiation characteristics of an object interposed between the subject and a radiation detector that detects the radiation image, the radiation characteristics corresponding to a body thickness distribution of the subject; deriving a primary ray distribution and a scattered ray distribution of radiation contained in each of the two radiographic images using imaging conditions, the body thickness distribution, and radiation characteristics of the object; deriving an error between the sum of the primary ray distribution and the scattered ray distribution for each of the two radiological images and the pixel value at each position of the two radiological images, updating the body thickness distribution so that the error becomes less than a predetermined threshold value, and repeatedly deriving the radiation characteristics based on the updated body thickness distribution and deriving the primary ray distribution and the scattered ray distribution included in each of the two radiological images; 10. The estimation device according to claim 9, wherein the estimation is performed by subtracting the scattered radiation distribution when the error becomes less than a predetermined threshold value from each of the two radiation images.

11. the scattered radiation removal process deriving a first primary ray distribution and a scattered ray distribution of the radiation that has passed through the subject using the two radiation images; deriving a second primary ray distribution and a scattered ray distribution of the radiation that has passed through an object, using the first primary ray distribution and the scattered ray distribution, and radiation characteristics of an object that is located between the subject and a radiation detector that detects the radiation image; The estimation device according to claim 9 , wherein the estimation is performed by deriving a radiation image after transmission through the subject and the object using the second primary ray distribution and the scattered ray distribution.

12. the scattered radiation removal process derives a region detection image by detecting, in the two radiation images, a subject region where the radiation has passed through the subject and reached a radiation detection unit, and a direct radiation region where the radiation has directly reached the radiation detection unit without passing through the subject; deriving a scattered radiation image relating to the scattered radiation component based on the region detection image and scattered radiation spread information relating to the spread of the scattered radiation; The estimation device according to claim 9 , wherein the estimation is performed by subtracting the scattered radiation image from the two radiation images.

13. The learning emphasis image is deriving processing details of a first graininess suppression processing for a first radiographic image having a higher S / N ratio of the two radiographic images; deriving processing details of a second graininess suppression process for a second radiographic image having a low S / N ratio based on processing details of the first graininess suppression process; performing graininess suppression processing on the first radiographic image based on processing details of the first graininess suppression processing; performing graininess suppression processing on the second radiographic image based on processing details of the second graininess suppression processing; The estimation device according to claim 2 , wherein the estimation is derived using the two radiation images that have been subjected to the graininess suppression processing.

14. The processing content of the first graininess suppression processing is as follows: The estimation device according to claim 13 , wherein the estimation value is derived based on a physical quantity map of the subject that is derived based on at least one of the first radiographic image and the second radiographic image.

15. An estimation method for deriving an estimation result of at least one emphasized image in which a specific composition of a subject is emphasized from a plain radiographic image obtained by plain radiographing the subject including a plurality of compositions, using a trained neural network that derives an estimation result of at least one emphasized image in which a specific composition of the subject is emphasized from the plain radiographic image, the method comprising: the trained neural network is constructed by using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions and a training emphasized image derived from the two radiation images in which a specific composition of the subject is emphasized, inputting the two radiation images into the neural network to output an emphasized image from the neural network, and training the neural network so that a difference between the emphasized image output from the neural network and the training emphasized image becomes small.

16. An estimation method for deriving an estimation result of at least one emphasized image in which a specific composition of a subject is emphasized from a plain radiographic image obtained by plain radiographing the subject including a plurality of compositions, using a trained neural network that derives an estimation result of at least one emphasized image in which a specific composition of the subject is emphasized from the plain radiographic image, the method comprising: the trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training emphasized image in which a specific composition of the subject is emphasized and which is derived by an energy subtraction process of performing weighted subtraction between the two radiation images; the emphasized image is at least one of a bone image in which a bone portion of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, The learning emphasis image is recognizing the bones and the soft tissue of the subject using at least one of the two radiation images; deriving attenuation coefficients for the bones and the soft tissue using the recognition results of the bones and the soft tissue and the two radiological images; The attenuation coefficient is derived by performing the energy subtraction process using the attenuation coefficient.

17. An estimation method for deriving an estimation result of at least one emphasized image in which a specific composition of a subject is emphasized from a plain radiographic image obtained by plain radiographing the subject including a plurality of compositions, using a trained neural network that derives an estimation result of at least one emphasized image in which a specific composition of the subject is emphasized from the plain radiographic image, the method comprising: the trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training emphasized image in which a specific composition of the subject is emphasized and which is derived by an energy subtraction process of performing weighted subtraction between the two radiation images; the emphasized image is a bone image in which a bone portion of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, The learning emphasis image is deriving a new weighting coefficient to be used in the weighted subtraction based on pixel values ​​of the bone portion included in the bone image and pixel values ​​of the soft tissue included in the soft tissue image; performing the weighted subtraction on the two radiographic images using the new weighting coefficients to derive a new bone image and a new soft tissue image; An estimation method in which the weighting coefficients are derived by repeatedly deriving a new weighting coefficient based on the new bone image, a new weighting coefficient based on the new soft tissue image, and a new bone image and a new soft tissue image based on the new weighting coefficients.

18. An estimation method for deriving an estimation result of at least one emphasized image in which a specific composition of a subject is emphasized from a plain radiographic image obtained by plain radiographing the subject including a plurality of compositions, using a trained neural network that derives an estimation result of at least one emphasized image in which a specific composition of the subject is emphasized from the plain radiographic image, the method comprising: the trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training emphasized image in which a specific composition of the subject is emphasized and which is derived by an energy subtraction process of performing weighted subtraction between the two radiation images; the emphasized image is a bone image in which a bone portion of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, The learning emphasis image is deriving a difference between the value of (attenuation coefficient of the soft part)×(thickness of the soft part)+(attenuation coefficient of the bone part)×(thickness of the bone part) for each different energy distribution and each pixel value of the radiographic image while changing the attenuation coefficient for each different energy distribution for the soft part, the thickness of the soft part, the attenuation coefficient for each different energy distribution for the bone part, and the thickness of the bone part from initial values; Deriving an attenuation coefficient of the soft tissue and an attenuation coefficient of the bone for each of the different energy distributions, in which the difference is minimized or the difference is less than a predetermined threshold value; The attenuation coefficient of the soft tissue and the attenuation coefficient of the bone are derived by performing the energy subtraction process using a weighting coefficient derived based on the attenuation coefficient of the soft tissue and the attenuation coefficient of the bone.

19. An estimation method for deriving an estimation result of at least one emphasized image in which a specific composition of a subject is emphasized from a plain radiographic image obtained by plain radiographing the subject including a plurality of compositions, using a trained neural network that derives an estimation result of at least one emphasized image in which a specific composition of the subject is emphasized from the plain radiographic image, the method comprising: the trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training emphasized image in which a specific composition of the subject is emphasized and which is derived by an energy subtraction process of performing weighted subtraction between the two radiation images; The learning emphasis image is Deriving composition ratios of a plurality of compositions contained in the soft part of the subject; deriving an attenuation coefficient for each energy distribution that differs for each pixel of the two radiological images according to the composition ratio; The estimation method includes performing the energy subtraction process using weighting coefficients derived based on the derived attenuation coefficients.

20. An estimation program that causes a computer to execute a procedure for deriving an estimation result of at least one emphasized image in which a specific composition of a subject is emphasized from a plain radiographic image obtained by plain radiographing a subject including a plurality of compositions, using a trained neural network that derives an estimation result of at least one emphasized image in which a specific composition of the subject is emphasized from the plain radiographic image, the estimation program comprising: the trained neural network is constructed and trained by using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions and a training emphasized image derived from the two radiation images in which a specific composition of the subject is emphasized, inputting the two radiation images into the neural network to output an emphasized image from the neural network, and training the neural network so that a difference between the emphasized image output from the neural network and the training emphasized image becomes small.

21. An estimation program that causes a computer to execute a procedure for deriving an estimation result of at least one emphasized image in which a specific composition of a subject is emphasized from a plain radiographic image obtained by plain radiographing a subject including a plurality of compositions, using a trained neural network that derives an estimation result of at least one emphasized image in which a specific composition of the subject is emphasized from the plain radiographic image, the estimation program comprising: the trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training emphasized image in which a specific composition of the subject is emphasized and which is derived by an energy subtraction process of performing weighted subtraction between the two radiation images; the emphasized image is at least one of a bone image in which a bone portion of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, The learning emphasis image is recognizing the bones and the soft tissue of the subject using at least one of the two radiation images; deriving attenuation coefficients for the bones and the soft tissue using the recognition results of the bones and the soft tissue and the two radiological images; An estimation program that derives the attenuation coefficient by performing the energy subtraction process using the attenuation coefficient.

22. An estimation program that causes a computer to execute a procedure for deriving an estimation result of at least one emphasized image in which a specific composition of a subject is emphasized from a plain radiographic image obtained by plain radiographing a subject including a plurality of compositions, using a trained neural network that derives an estimation result of at least one emphasized image in which a specific composition of the subject is emphasized from the plain radiographic image, the estimation program comprising: the trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training emphasized image in which a specific composition of the subject is emphasized and which is derived by an energy subtraction process of performing weighted subtraction between the two radiation images; the emphasized image is a bone image in which a bone portion of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, The learning emphasis image is deriving a new weighting coefficient to be used in the weighted subtraction based on pixel values ​​of the bone portion included in the bone image and pixel values ​​of the soft tissue included in the soft tissue image; performing the weighted subtraction on the two radiographic images using the new weighting coefficients to derive a new bone image and a new soft tissue image; An estimation program in which the weighting coefficients are derived by repeatedly deriving a new weighting coefficient based on the new bone image, a new weighting coefficient based on the new soft tissue image, and a new bone image and a new soft tissue image based on the new weighting coefficients.

23. An estimation program that causes a computer to execute a procedure for deriving an estimation result of at least one emphasized image in which a specific composition of a subject is emphasized from a plain radiographic image obtained by plain radiographing a subject including a plurality of compositions, using a trained neural network that derives an estimation result of at least one emphasized image in which a specific composition of the subject is emphasized from the plain radiographic image, the estimation program comprising: the trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training emphasized image in which a specific composition of the subject is emphasized and which is derived by an energy subtraction process of performing weighted subtraction between the two radiation images; the emphasized image is a bone image in which a bone portion of the subject is emphasized and a soft part image in which a soft part of the subject is emphasized, The learning emphasis image is deriving a difference between the value of (attenuation coefficient of the soft part)×(thickness of the soft part)+(attenuation coefficient of the bone part)×(thickness of the bone part) for each different energy distribution and each pixel value of the radiographic image while changing the attenuation coefficient for each different energy distribution for the soft part, the thickness of the soft part, the attenuation coefficient for each different energy distribution for the bone part, and the thickness of the bone part from initial values; Deriving an attenuation coefficient of the soft tissue and an attenuation coefficient of the bone for each of the different energy distributions, in which the difference is minimized or the difference is less than a predetermined threshold value; an estimation program for deriving the attenuation coefficient of the soft tissue and the attenuation coefficient of the bone by performing the energy subtraction process using a weighting coefficient derived based on the attenuation coefficient of the soft tissue and the attenuation coefficient of the bone.

24. An estimation program that causes a computer to execute a procedure for deriving an estimation result of at least one emphasized image in which a specific composition of a subject is emphasized from a plain radiographic image obtained by plain radiographing a subject including a plurality of compositions, using a trained neural network that derives an estimation result of at least one emphasized image in which a specific composition of the subject is emphasized from the plain radiographic image, the estimation program comprising: the trained neural network is trained using, as training data, two radiation images obtained by photographing a subject with radiation having different energy distributions, and a training emphasized image in which a specific composition of the subject is emphasized and which is derived by an energy subtraction process of performing weighted subtraction between the two radiation images; The learning emphasis image is Deriving composition ratios of a plurality of compositions contained in the soft part of the subject; deriving an attenuation coefficient for each energy distribution that differs for each pixel of the two radiological images according to the composition ratio; The estimation program derives the attenuation coefficient by performing the energy subtraction process using a weighting coefficient derived based on the derived attenuation coefficient.

Images