Radiation image processing device, radiation image processing method, and radiation image processing program
An imaging condition acquisition unit acquires an imaging condition in a case in which a subject is imaged in a state in which an object is interposed between the subject and a radiation detector. A body thickness derivation unit derives a body thickness distribution of the subject based on the radiation image and the imaging condition. A characteristic acquisition unit acquires a radiation characteristic of the object in accordance with the body thickness distribution. A ray distribution unit derives a primary ray distribution and a scattered ray distribution of the radiation detected by the radiation detector by using the imaging condition, the body thickness distribution, and the radiation characteristic. A calculation unit updates the body thickness distribution such that an error between a sum of the primary ray distribution and the scattered ray distribution and a pixel value at each position of the radiation image is smaller than a predetermined threshold value, and repeats the derivation of the radiation characteristic based on the updated body thickness distribution and the derivation of the primary ray distribution and the scattered ray distribution.
1. A radiation image processing device comprising:
at least one processor,
wherein the processor is configured to
acquire an imaging condition in a case in which a radiation image of a subject is acquired by a radiation detector detecting radiation transmitted through the subject in a state in which an object is interposed between the subject and the radiation detector,
derive a body thickness distribution of the subject based on the radiation image and the imaging condition,
derive a radiation characteristic of the object based on the body thickness distribution,
derive a primary ray distribution and a scattered ray distribution of the radiation detected by the radiation detector by using the imaging condition, the body thickness distribution, and the radiation characteristic of the object, and
derive an error between a sum of the primary ray distribution and the scattered ray distribution and a pixel value at each position of the radiation image, update the body thickness distribution such that the error is smaller than a predetermined threshold value, and repeat derivation of the radiation characteristic based on the updated body thickness distribution and derivation of the primary ray distribution and the scattered ray distribution.
2. The radiation image processing device according to claim 1 ,
wherein the processor is configured to output a processed radiation image including, as a pixel value, the primary ray distribution derived based on the body thickness distribution of the subject in which the error is smaller than the threshold value.
3. The radiation image processing device according to claim 1 , further comprising:
a storage unit that stores the radiation characteristic measured in advance by using an imaging system that acquires the radiation image,
wherein the processor is configured to acquire the radiation characteristic stored in the storage unit.
4. The radiation image processing device according to claim 1 ,
wherein the object is at least one of an imaging table on which the subject is placed, a top plate, a scattered ray removal grid, or an air layer.
5. The radiation image processing device according to claim 1 ,
wherein the radiation characteristic is a primary ray transmittance and a scattered ray transmittance.
6. The radiation image processing device according to claim 1 ,
wherein the processor is configured to
acquire a processed radiation image for each of two radiation images based on radiation having different energy distributions and transmitted through the subject, and
derive a subtraction image obtained by extracting a specific structure of the subject by performing weighting subtraction between corresponding pixels of two processed radiation images.
7. The radiation image processing device according to claim 1 ,
wherein the processor is configured to
acquire a processed radiation image for each of two radiation images based on radiation having different energy distributions and transmitted through a subject including a bone part and a soft part,
derive a bone part image obtained by extracting the bone part of the subject by performing weighting subtraction between corresponding pixels of two processed radiation images, and
derive bone mineral information representing a bone mineral density in a bone region in the bone part image for each pixel of the bone region based on the imaging condition, the body thickness distribution, and a pixel value of the bone region.
8. The radiation image processing device according to claim 1 ,
wherein the processor is configured to
irradiate the subject with the radiation from a plurality of projection angles to acquire a processed radiation image for each of a plurality of radiation images corresponding to each of the plurality of projection angles, and
reconstruct a plurality of the processed radiation images to generate a tomographic image of the subject.
9. A radiation image processing method comprising:
acquiring an imaging condition in a case in which a radiation image of a subject is acquired by a radiation detector detecting radiation transmitted through the subject in a state in which an object is interposed between the subject and the radiation detector;
deriving a body thickness distribution of the subject based on the radiation image and the imaging condition;
deriving a radiation characteristic of the object based on the body thickness distribution;
deriving a primary ray distribution and a scattered ray distribution of the radiation detected by the radiation detector by using the imaging condition, the body thickness distribution, and the radiation characteristic of the object; and
deriving an error between a sum of the primary ray distribution and the scattered ray distribution and a pixel value at each position of the radiation image, updating the body thickness distribution such that the error is smaller than a predetermined threshold value, and repeating derivation of the radiation characteristic based on the updated body thickness distribution and derivation of the primary ray distribution and the scattered ray distribution.
10. A non-transitory computer-readable storage medium that stores a radiation image processing program causing a computer to execute a processing procedure comprising:
a procedure of acquiring an imaging condition in a case in which a radiation image of a subject is acquired by a radiation detector detecting radiation transmitted through the subject in a state in which an object is interposed between the subject and the radiation detector;
a procedure of deriving a body thickness distribution of the subject based on the radiation image and the imaging condition;
a procedure of deriving a radiation characteristic of the object based on the body thickness distribution;
a procedure of deriving a primary ray distribution and a scattered ray distribution of the radiation detected by the radiation detector by using the imaging condition, the body thickness distribution, and the radiation characteristic of the object; and
a procedure of deriving an error between a sum of the primary ray distribution and the scattered ray distribution and a pixel value at each position of the radiation image, updating the body thickness distribution such that the error is smaller than a predetermined threshold value, and repeating derivation of the radiation characteristic based on the updated body thickness distribution and derivation of the primary ray distribution and the scattered ray distribution.