Radiation image processing device, method and program

Abstract

To provide a radiation image processing apparatus, a method and a program which can match contrasts of two radiation images having different imaging conditions and characteristics of devices that have acquired the images with each other.SOLUTION: A processor acquires two radiation images having a contrast based on a first characteristic about a first imaging device, derives the body thickness distribution of a subject on the basis of at least one of the two radiation images, removes a scattered ray component from the two radiation images on the basis of the first characteristic, derives a first bone part image and a first soft part image expressing a bone part and a soft tissue of the subject from the two radiation images, converts the first bone part image and the first soft part image into a second bone part image and a second soft part image having a contrast based on a second characteristic on the basis of the first characteristic, the second characteristic about a second imaging device and the body thickness distribution, and derives a processed radiation image having a contrast based on the second characteristic by adding the second bone part image and the second soft part image.SELECTED DRAWING: Figure 3

Description

[Technical Field]

[0001] The present disclosure relates to a radiation image processing device, method, and program. [Background technology]

[0002] Conventionally, when making a diagnosis using a radiological image, comparative interpretation is performed using a patient's past radiological images. For example, by displaying and comparatively interpreting the patient's radiological images acquired in the latest examination and in the past examinations, it is possible to confirm the degree of progression of a lesion or to detect an abnormality at an early stage.

[0003] Also proposed is a computer-aided image diagnosis system (CAD, hereinafter referred to as CAD) that automatically detects structures such as abnormal shadows in an image and highlights the detected structures.

[0004] Furthermore, when capturing a radiographic image of a subject, particularly if the subject is thick, there is a problem in that radiation is scattered within the subject, generating scattered rays, which reduces the contrast of the acquired radiographic image. For this reason, scattered ray removal processing is performed to remove scattered ray components contained in the radiographic image (see, for example, Patent Document 1). Specifically, the scattered ray components of the radiographic image are derived based on the radiation attenuation coefficient of the subject, and the derived scattered ray components are subtracted from the radiographic image, thereby performing the scattered ray removal processing. [Prior art documents] [Patent documents]

[0005] [Patent Document 1] Japanese Patent Application Laid-Open No. 2015-043959 Summary of the Invention [Problem to be solved by the invention]

[0006] In order to perform the above-mentioned comparative interpretation with high accuracy, it is important to reproduce the contrast of the two radiographic images to be compared. The contrast of a radiographic image varies depending on the characteristics of the imaging device, such as the energy of the radiation irradiated onto the subject, the top surface of the imaging table on which the subject is placed, and the anti-scatter grid for removing scattered radiation components contained in the radiation that has passed through the subject. Therefore, if the two radiographic images to be compared are acquired using the same imaging device under the same imaging conditions, the contrast of the two radiographic images will match, allowing for accurate comparative interpretation.

[0007] However, if two radiographic images for comparative interpretation are acquired using different imaging devices, the contrast between the two will differ. Furthermore, if the imaging conditions (tube voltage of the radiation source, imaging distance, tube current, etc.) are different, the contrast between the two radiographic images will also differ. Furthermore, the contrast of a radiographic image also changes depending on the scattered radiation contained in the radiographic image.

[0008] As described above, when the contrast of two radiographic images differs, accurate comparative interpretation is not possible. Experienced physicians can interpret radiographic images taking into account the difference in contrast of radiographic images due to differences in equipment, but this places a heavy burden on the physician. Furthermore, when detecting abnormal shadows from radiographic images using the above-mentioned CAD, the accuracy of detecting abnormal shadows may decrease if the contrast of the radiographic images input into the CAD differs.

[0009] The present disclosure has been made in consideration of the above circumstances, and is intended to be illustrative of the principles of the present invention, given that the imaging conditions and the characteristics of the device that acquired the image differ. The present invention aims to make it possible to match the contrast of two different radiographic images. [Means for solving the problem]

[0010] A radiological image processing method according to the present disclosure includes at least one processor, The processor acquires two radiographic images having contrast based on a first characteristic related to the first imaging device, the two radiographic images being acquired by imaging an object including a soft tissue and a bone portion with a first imaging device using radiation having different energy distributions; deriving a body thickness distribution of the subject based on at least one of the two radiographic images; removing, from the two radiation images, a scattered radiation component that is included in the radiation that has passed through the subject and that is scattered by the subject, based on the first characteristic; deriving a first bone image representing the bone tissue of the subject and a first soft tissue image representing the soft tissue of the subject by performing weighted subtraction on the two radiographic images from which the scattered radiation components have been removed; converting the first bone image and the first soft tissue image into a second bone image and a second soft tissue image having contrast based on the second characteristic, based on the first characteristic, a second characteristic related to a second imaging device different from the first imaging device, and a body thickness distribution; The second bone image and the second soft tissue image are added together to derive a processed radiographic image having contrast based on the second characteristic.

[0011] In the radiological image processing device according to the present disclosure, the processor derives scattered ray components according to the second characteristic based on the second characteristic and the body thickness distribution, Furthermore, the derived scattered radiation component may be used to derive a processed radiographic image.

[0012] Furthermore, in the radiation image processing device according to the present disclosure, the first characteristic includes the energy of radiation used in the first imaging device, a radiation attenuation coefficient corresponding to the body thickness distribution of an object interposed between the subject and a radiation detector that detects the radiation transmitted through the subject in the first imaging device, a ratio of scattered ray components contained in the radiation transmitted through the subject corresponding to the body thickness distribution, and a point spread function corresponding to the body thickness distribution; The second characteristic may include the energy of the radiation used in the second imaging device, a radiation attenuation coefficient depending on the body thickness distribution of an object interposed between the subject and the radiation detector that detects the radiation that has passed through the subject in the second imaging device, a ratio depending on the body thickness distribution of scattered radiation components contained in the radiation that has passed through the subject, and a point spread function depending on the body thickness distribution.

[0013] In the radiation image processing device according to the present disclosure, the processor may display the processed radiation image and the radiation image of the subject acquired by the second imaging device.

[0014] A radiological image processing method according to the present disclosure includes acquiring two radiological images having contrast based on a first characteristic of a first imaging device, the two radiological images being acquired by imaging a subject including soft tissue and bones using radiation with different energy distributions with a first imaging device; deriving a body thickness distribution of the subject based on at least one of the two radiographic images; removing, from the two radiation images, a scattered radiation component that is included in the radiation that has passed through the subject and that is scattered by the subject, based on the first characteristic; deriving a first bone image representing the bone tissue of the subject and a first soft tissue image representing the soft tissue of the subject by performing weighted subtraction on the two radiographic images from which the scattered radiation components have been removed; converting the first bone image and the first soft tissue image into a second bone image and a second soft tissue image having contrast based on the second characteristic, based on the first characteristic, a second characteristic related to a second imaging device different from the first imaging device, and a body thickness distribution; The second bone image and the second soft tissue image are added together to derive a processed radiographic image having contrast based on the second characteristic.

[0015] The radiation image processing method according to the present disclosure may be provided as a program for causing a computer to execute the method. [Effects of the Invention]

[0016] According to the present disclosure, it is possible to match the contrast of two radiographic images that are obtained under different imaging conditions and with different device characteristics. [Brief explanation of the drawings]

[0017] [Figure 1] FIG. 1 is a schematic block diagram illustrating a configuration of a radiographic image capturing system to which a radiographic image processing apparatus according to an embodiment of the present disclosure is applied. [Figure 2] FIG. 1 is a diagram showing a schematic configuration of a radiation image processing apparatus according to an embodiment of the present disclosure. [Figure 3] FIG. 1 is a diagram showing a functional configuration of a radiation image processing apparatus according to an embodiment of the present disclosure. [Figure 4] FIG. 1 is a diagram for explaining photographing of a reference object; [Figure 5] Diagram showing the spectrum of radiation [Figure 6] A diagram showing the radiation attenuation coefficients of soft tissue, bone tissue, and aluminum in the human body for radiation energy. [Figure 7] A diagram showing the relationship between the thickness of a reference object and the radiation attenuation coefficient [Figure 8] Diagram showing the scattered radiation model [Figure 9] FIG. 1 is a diagram schematically illustrating the processing performed by a conversion unit and a derivation unit. [Figure 10] Diagram showing the reading screen [Figure 11] A flowchart showing the processing performed in this embodiment DETAILED DESCRIPTION OF THE INVENTION

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

[0019] The imaging device 1A is an imaging device for performing energy subtraction using a so-called one-shot method, in which radiation such as X-rays emitted from a radiation source 2 and transmitted through a subject H lying supine on an imaging table 3 is irradiated at different energies onto a first radiation detector 5 and a second radiation detector 6. During imaging, as shown in FIG. 1 , an anti-scatter grid (hereinafter simply referred to as the grid) 4, a first radiation detector 5, a radiation energy converting filter 7 made of a copper plate or the like, and a second radiation detector 6 are arranged in this order from the side closest to the radiation source 2, and the radiation source 2 is driven. The first and second radiation detectors 5 and 6 are in close contact with the radiation energy converting filter 7. The grid 4, the first radiation detector 5, the radiation energy converting filter 7, and the second radiation detector 6 are removably attached to a bottom of a top plate 3A of the imaging table 3 by attachment parts 3B.

[0020] As a result, the first radiation detector 5 acquires a first radiographic image G1 of the subject H using low-energy radiation including so-called soft rays. The second radiation detector 6 acquires a second radiographic image G2 of the subject H using high-energy radiation from which the soft rays have been removed. The first and second radiographic images G1 and G2 are input to the radiographic image processing device 10.

[0021] 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 upon direct exposure to radiation, or so-called indirect type radiation detectors that convert radiation into visible light and then convert the visible light into electric charge signals. The radiation image signal readout method is preferably a TFT readout method in which the radiation image signal is read out by turning a TFT (thin film transistor) switch on and off, or an optical readout method in which the radiation image signal is read out by irradiating the detector with readout light, but is not limited to these, and other methods may also be used.

[0022] In the imaging device 1A, there are cases where imaging of the subject H is performed with only one radiation detector attached to the attachment portion 3B.

[0023] The grid 4 is made up of alternately arranged radiation-opaque materials such as lead and radiation-transmittable interspace materials such as aluminum or fiber at a fine grid density of, for example, 4.0 lines / mm. By using the grid 4, the scattered ray components of the radiation that has passed through the subject H can be removed, but cannot be completely removed. For this reason, the first and second radiographic images G1, G2 contain not only the primary ray components of the radiation that has passed through the subject H, but also the scattered ray components.

[0024] The primary ray component is a signal component of a pixel value represented by radiation that has passed through the subject H and reached the radiation detector without being scattered by the subject H. On the other hand, the scattered ray component is a signal component of a pixel value represented by radiation that has passed through the subject H and reached the radiation detector after being scattered by the subject H.

[0025] Since the imaging device 1B has the same configuration as the imaging device 1A, detailed description thereof will be omitted here. However, the imaging devices 1A and 1B differ in the material and thickness of the tabletop used and the characteristics of the grid used.

[0026] The image storage system 9 is a system that stores image data of radiation images acquired by the radiographic devices 1A and 1B. For example, the image storage system 9 stores multiple radiation images of the same patient that are captured on different dates and times. The image storage system 9 extracts an image from the stored radiation images in response to a request from the radiation image processing device 10 and transmits it to the device that issued the request. A specific example of the image storage system 9 is a PACS (Picture Archiving and Communication System).

[0027] In this embodiment, it is assumed that radiation images acquired in a past examination of the patient, who is the subject H, are stored in the image storage system 9. It is also assumed that the radiation images of the past examination (hereinafter referred to as past radiation images) were acquired by the imaging device 1B. In the radiation image processing device 10 according to this embodiment, the past radiation images of the patient are acquired from the image storage system 9, and a comparative interpretation is performed between the radiation image acquired in the latest examination and the past radiation images.

[0028] Next, the radiographic image processing apparatus according to this embodiment will be described. First, the hardware configuration of the radiographic image processing apparatus according to this embodiment will be described with reference to Fig. 2. As shown in Fig. 2, the radiographic image processing apparatus 10 is a computer such as a workstation, a server computer, or a personal computer, and includes a CPU (Central Processing Unit) 11 The radiation image processing apparatus 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 interface 16 (not shown) for connecting to a network. / F (Interface) 17. CPU 11, storage 13, display 14, input The output 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.

[0029] 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 radiographic image processing program 12 installed in the radiographic image processing device 10. The CPU 11 reads out the radiographic image processing program 12 from the storage 13, loads it into the memory 16, and executes the loaded radiographic image processing program 12.

[0030] The radiation image processing program 12 is stored in a state where it can be accessed from outside, in a storage device of a server computer connected to a network, or in a network storage, and is downloaded and installed in a computer constituting the radiation image processing apparatus 10 upon request. Alternatively, the program may be recorded on a recording medium such as a DVD (Digital Versatile Disc) or a CD-ROM (Compact Disc Read Only Memory) and distributed, and may be downloaded from the recording medium. It is installed in the computer that constitutes the radiation image processing apparatus 10.

[0031] Next, the functional configuration of the radiological image processing apparatus according to this embodiment will be described. Fig. 3 is a diagram showing the functional configuration of the radiological image processing apparatus according to this embodiment. As shown in Fig. 3, the radiological image processing apparatus 10 includes an image acquisition section 21, a scattered radiation removal section 22, a subtraction section 23, a conversion section 24, a derivation section 25, a display control section 26, and a characteristic derivation section 27. The CPU 11 executes the radiological image processing program 12 to function as the image acquisition section 21, the scattered radiation removal section 22, the subtraction section 23, the conversion section 24, the derivation section 25, the display control section 26, and the characteristic derivation section 27.

[0032] The image acquisition unit 21 causes the imaging device 1A to perform energy subtraction imaging of the subject H, thereby acquiring a first radiographic image G1 and a second radiographic image G2 of the subject H from the first and second radiation detectors 5, 6. When acquiring the first radiographic image G1 and the second radiographic image G2, imaging conditions are set, such as the imaging dose, radiation quality, tube voltage (kV), SID (Source Image Receptor Distance) which is the distance between the radiation source 2 and the surfaces of the first and second radiation detectors 5, 6, SOD (Source Object Distance) which is the distance between the radiation source 2 and the surface of the subject H, and the presence or absence of an anti-scatter grid. The imaging conditions may be set by an operator via the input device 15.

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

[0034] 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 imaging device 1A and the imaging conditions are transmitted to and stored in image storage system 9. The radiographic images acquired by imaging device 1B and the imaging conditions used when the radiographic images were acquired are also transmitted to and stored in image storage system 9.

[0035] The scattered radiation removal unit 22 removes scattered radiation components from each of the first radiographic image G1 and the second radiographic image G2 acquired by the image acquisition unit 21. The removal of scattered radiation components will be described below. As a method for removing scattered radiation components, any method can be used, for example, the method described in JP 2015-043959 A. The scattered radiation removal process using the method described in JP 2015-043959 A will be described below. In the following description, the first and second radiographic images from which scattered radiation components have been removed will also be referred to as G1 and G2, respectively.

[0036] First, the scattered radiation removal unit 22 acquires a virtual model of the subject H having an initial body thickness distribution Ts(x,y). The virtual model is data virtually representing the subject H, in which the body thickness according to the initial body thickness distribution Ts(x,y) is associated with the coordinate position of each pixel in the first radiographic image G1. The virtual model of the subject H having the initial body thickness distribution Ts(x,y) is assumed to be stored in advance in the storage 13, but may also be acquired from an external server on which the virtual model is saved.

[0037] Next, the scattered radiation removal unit 22 derives an estimated primary radiation image Ip(x,y) that is an estimate of a primary radiation image obtained by imaging the virtual model, and an estimated scattered radiation image Is(x,y) that is an estimate of a scattered radiation image obtained by imaging the virtual model, based on the virtual model, as shown in the following equations (1) and (2). Furthermore, the scattered radiation removal unit 22 derives an estimated image Im(x,y) that is an estimate of a first radiographic image G1 obtained by imaging the subject H, by combining the estimated primary radiation image Ip(x,y) and the estimated scattered radiation image Is(x,y), as shown in the following equation (3). Ip(x,y) = Io(x,y)×exp(-μ1Soft(T(x,y))×T(x,y)) (1) Is(x,y) = Io(x,y)×STPR1(T(x,y))*PSF1(T(x,y)) (2) Im(x,y) = Is(x,y)+Ip(x,y) (3)

[0038] where (x, y) are the coordinates of the pixel position in the first radiographic image G1, Io(x, y) is the pixel value of the first radiographic image G1 at the pixel position (x, y), Ip(x, y) is the primary ray component at the pixel position (x, y), and Is(x, y) is the scattered ray component at the pixel position (x, y). When deriving the first estimated image Im(x, y), the initial body thickness distribution Ts(x, y) is used as the body thickness distribution T(x, y) in equations (1) and (2).

[0039] In addition, μ1Soft(T(x,y)) in equation (1) is the human μSoft(T(x,y)) is the attenuation coefficient depending on the thickness distribution (x,y) of the soft tissue of the body. may be found in advance experimentally or by simulation and stored in the storage 13. STPR1(T(x,y)) in equation (2) is the scatter-to-primary ratio of the scattered radiation dose to the primary radiation dose contained in the radiation after passing through the subject H having a body thickness distribution T(x,y). STPR1(T(x,y)) may also be found in advance experimentally or by simulation and stored in the storage 13.

[0040] Furthermore, PSF1(T(x,y)) in equation (2) is a point spread function that represents the distribution of scattered rays spreading from one pixel according to the body thickness distribution T(x,y), and is defined according to the energy characteristics of the radiation. Furthermore, * is an operator indicating a convolution operation. PSF1 also varies depending on the distribution of the irradiation field in the imaging device 1, the distribution of the composition of the subject H, the exposure dose during imaging, the tube voltage, the imaging distance, and the characteristics of the radiation detectors 5 and 6. For this reason, PSF1 can be experimentally determined in advance for each energy characteristic of the radiation used by the imaging device 1A according to the irradiation field information, subject information, imaging conditions, etc., and stored in storage 13.

[0041] The attenuation coefficients μ1Soft, STPR1, and PSF1 are related to the first imaging device in this disclosure. This is an example of the first characteristic.

[0042] Next, the scattered radiation removal unit 22 corrects the initial body thickness distribution Ts(x, y) of the virtual model so that the difference between the estimated image Im and the first radiographic image G1 becomes smaller. The scattered radiation removal unit 22 continues to correct the initial body thickness distribution Ts(x, y) of the virtual model until the difference between the estimated image Im and the first radiographic image G1 satisfies a predetermined termination condition. By repeatedly deriving the thickness distribution T(x,y), scattered ray component Is(x,y), and primary ray component Ip(x,y), the body thickness distribution T(x,y), scattered ray component Is(x,y), and primary ray component Ip(x,y) are updated. When the termination condition is satisfied, the scattered ray removal unit 22 subtracts the scattered ray component Is(x,y) derived using equation (2) from the first radiographic image G1. This removes the scattered ray component contained in the first radiographic image G1. The body thickness distribution T(x,y) derived when the termination condition is satisfied is used in various calculations, which will be described later.

[0043] On the other hand, the scattered radiation removal unit 22 also performs the scattered radiation removal process on the second radiographic image G2 in the same manner as on the first radiographic image G1.

[0044] The derivation of the attenuation coefficients μ1Soft and STPR1 will be described below. The attenuation coefficients μ1Soft and STPR1 are derived by the characteristic derivation unit 27. Derivation of the attenuation coefficients μ1Soft and STPR1 In this case, the image acquisition unit 21 acquires a reference image K0 by having the imaging device 1A capture an image of a reference object simulating a human body. At this time, only one radiation detector may be used. Note that if the reference image K0 is stored in the image storage system 9, the image acquisition unit 21 acquires the reference image K0 from the image storage system 9. Also, in the following description, the reference number "1" is omitted for the sake of generality.

[0045] FIG. 4 is a diagram for explaining the imaging of a reference object. As shown in FIG. 4, the reference object 35 has portions with different thicknesses, such as 5 cm, 10 cm, and 20 cm, and is made of a material with a radiation transmittance similar to that of the soft tissues (fat and muscle) of the human body. Therefore, the reference object 35 simulates the radiation characteristics of the human body. Here, the soft tissue is made of a mixture of muscle and fat in a certain ratio. The mixture ratio of muscle and fat varies depending on gender, physique, etc., but can be defined by the average body fat percentage (25%). Therefore, a material such as acrylic, which corresponds to a composition where muscle and fat are mixed in a ratio of 0.75 and fat in a ratio of 0.25, is used as the reference object 35.

[0046] 4, when acquiring the reference image K0, the reference object 35 is placed on the top plate 3A of the imaging table 3, and the radiation source 2 is driven to irradiate the radiation detector (here, the first radiation detector 5) with radiation via the grid 4, thereby causing the image acquisition unit 21 to acquire the reference image K0. The pixel value of each pixel in the reference image K0 includes a primary ray component based on the radiation that has traveled straight through the reference object 35 and a scattered ray component based on the radiation scattered by the reference object 35.

[0047] Note that the reference object 35 is not limited to a single object having different thicknesses as shown in Fig. 4. Multiple reference objects each having different thicknesses may be used. In this case, the reference image K0 may be acquired by photographing the multiple reference objects at once, or the multiple reference objects may be photographed separately to acquire reference images corresponding to each of the reference objects.

[0048] When acquiring the reference image K0, imaging conditions such as the imaging dose, tube voltage, SID (Source Image Receptor Distance) which is the distance between the radiation source 2 and the surfaces of the first and second radiation detectors 5 and 6, and whether or not a grid 4 is used are also set.

[0049] The characteristic derivation unit 27 acquires the imaging conditions set during imaging. Furthermore, the characteristic derivation unit 27 acquires the radiation energy characteristics during imaging of the reference object 35 in order to derive the attenuation coefficients μSoft and STPR. The radiation energy characteristics may be acquired from the imaging device 1A, or may be stored in the image storage system 9 and acquired from the image storage system 9. Note that the nominal value of the imaging device 1A may be used for the energy characteristics, but since there are individual differences in characteristics between devices, it is preferable to measure them in advance using a semiconductor dosimeter.

[0050] Here, the energy characteristics are defined by one of the following: (i) the spectrum of the radiation emitted from the radiation source 2, (ii) the tube voltage [kV] and total filtration [mmAl equivalent], and (iii) the tube voltage [kV] and aluminum half-value layer [mmAl]. The radiation spectrum is a plot of the relationship between the radiation energy [keV] and the relative number of radiation photons. The tube voltage represents the maximum value of the emitted radiation energy distribution. The total filtration is the filtration of each component of the imaging device 1A, such as the radiation generator and collimator in the radiation source 2, converted into aluminum thickness. The greater the total filtration, the greater the effect of beam hardening in the imaging device 1, resulting in a greater concentration of high-energy components in the radiation wavelength distribution. The half-value layer is defined as the thickness of aluminum required to attenuate the dose by half for the emitted radiation energy distribution. The thicker the aluminum half-value layer, the greater the concentration of high-energy components in the radiation wavelength distribution.

[0051] Figure 5 shows the spectrum of radiation. In Figure 5, the spectrum corresponds to a tube voltage of 90 kV and a total filtration volume of 2.5 mmAl. The total filtration volume of 2.5 mmAl corresponds to a half-value layer of 2.96 mmAl.

[0052] The characteristic derivation unit 27 uses the energy characteristics of the radiation to derive the relationship between the thickness of the reference object 35 and the radiation attenuation coefficient of the reference object 35, reflecting the influence of beam hardening of objects present between the reference object 35 and the radiation detector 5.

[0053] The characteristic deriving unit 27 first derives the energy spectrum of the radiation from the acquired energy characteristics of the radiation using the well-known Tucker approximation formula, etc. If the acquired energy characteristics are the energy spectrum of the radiation, the acquired energy spectrum may be used as is.

[0054] Furthermore, the characteristic deriving unit 27 derives a radiation attenuation coefficient that depends on the thickness of the reference object 35 by simulating the radiation spectrum using the radiation attenuation characteristics of the soft tissue of the human body.

[0055] Here, when the energy spectrum of the radiation emitted from the radiation source 2 is Sin(E) and the thickness of the reference object 35 is t, the radiation dose Xbody(t) after passing through the reference object 35 can be calculated by the following formula (4) using the radiation attenuation coefficient μSoft(E) of the soft tissue of the human body. The radiation attenuation coefficients of the soft tissue, bone tissue, and aluminum of the human body with respect to the radiation energy are known, for example, as shown in FIG. 6. Aluminum is the interspace material of the grid 4. Here, FIG. 6 also shows the radiation attenuation coefficient of acrylic (Polymethyl methacrylate, PMMA), which is the material of the reference object 35. As shown in FIG. 6, The radiation attenuation coefficient of acrylic is almost identical to that of the soft tissues of the human body.

number

[0056] On the other hand, when imaging a reference object 35 in the imaging device 1A, the top plate 3A and grid 4 are present between the reference object 35 and the radiation detectors 5, 6. The material of the top plate 3A is acrylic, and the interspace material of the grid 4 is aluminum. If the radiation attenuation coefficient of acrylic is μPMMA(E), the thickness of the top plate 3A (i.e., the thickness of the acrylic) is tPMMA, the radiation attenuation characteristic of aluminum is μAl(E), and the thickness of the grid 4 (i.e., aluminum) is tAl, the X-ray dose Xout( t) is expressed by the following equation (5).

number

[0057] In addition, when the material of the tabletop 3A and the interspace material of the grid 4 are unknown, the X-ray dose Xout(t) after passing through the tabletop 3A and the grid 4 is calculated by the above formula (5). In this case, the energy characteristics (kV, TF0) of the radiation emitted from the radiation source 2 and the energy characteristics (kV, TF1) of the radiation after passing through the tabletop 3A and the grid 4 can be measured using a dosimeter, and the X-ray dose Xout(t) after passing through the tabletop 3A and the grid 4 can be calculated using the following equation (5-1) using the energy characteristics (kV, TF0) and the energy characteristics (kV, TF1). Energy characteristics represent the total filtration (mmAl equivalent) of radiation emitted by a certain tube voltage (kV).

number

[0058] The radiation attenuation coefficient of the reference object 35 in the imaging system including the tabletop 3A and the grid 4 is expressed as an exponential decay of the rate attenuation of radiation after passing through the reference object 35, as shown in the following equation (6), with the radiation dose when the reference object 35 is not present (i.e., when the thickness of the reference object 35 is 0) as the reference.

number

[0059] By solving equation (6) for the radiation attenuation coefficient μSoft(t) of soft tissue as shown in the following equation (7), the relationship between the thickness t of the reference object 35 and the radiation attenuation coefficient can be derived.

number

[0060] The reference object 35 has a plurality of different thicknesses. Therefore, the characteristic deriving unit 27 derives the radiation attenuation coefficient for each of the plurality of thicknesses of the reference object 35 using Equation (7). Then, for the radiation attenuation coefficients of thicknesses not present in the reference object 35, the characteristic deriving unit 27 derives the relationship between the thickness t of the reference object 35 and the radiation attenuation coefficient by performing an interpolation calculation using the radiation attenuation coefficients of thicknesses present in the reference object 35. FIG. 7 shows the relationship between the thickness t of the reference object 35 and the radiation attenuation coefficient when the tube voltage is 90 kV and the total filtration volume is 2.5 mmAl. The characteristic deriving unit 27 derives the relationship between the thickness of the reference object 35 and the radiation attenuation coefficient for each radiation energy characteristic and stores the derived relationship in the storage 13.

[0061] Furthermore, the characteristic derivation unit 27 derives the radiation attenuation coefficient corresponding to the thickness of the reference object 35 based on the derived relationship between the thickness of the reference object 35 and the radiation attenuation coefficient. Furthermore, based on the radiation attenuation coefficient corresponding to the thickness of the reference object 35, the characteristic derivation unit 27 derives the primary ray component included in the reference image K0.

[0062] Here, the pixel value of each pixel in the reference image K0 is I0o(x,y), and When the thickness of the corresponding reference object 35 is T0(x,y) and the radiation attenuation coefficient derived by the above formula (7) for the thickness T0(x,y) of each pixel of the reference image K0 is μSoft0(x,y), the characteristic deriving unit 27 derives the primary ray component I0p(x,y) contained in the pixel value of each pixel of the reference image K0 by the following formula (8). Note that the reference object 35 has a plurality of thicknesses in stages. Therefore, the characteristic deriving unit 27 calculates the primary ray component I0p(x, y) for each thickness of the reference object 35. Note that, for primary ray components corresponding to thicknesses not present in the reference object 35, the characteristic deriving unit 27 may derive the relationship between the thickness of the reference object 35 and the primary ray components by performing an interpolation calculation using primary ray components of thicknesses present in the reference object 35. I0p(x,y) = I0o(x,y)×exp(-μSoft0(x,y)×T0(x,y)) (8)

[0063] Furthermore, the characteristic deriving unit 27 derives the scattered ray component contained in the reference object 35 based on the difference between the pixel value of the reference image K0 and the primary ray component. That is, the characteristic deriving unit 27 derives the scattered ray component I0s(x, y) by the following equation (9). Note that the reference object 35 has a plurality of thicknesses. Since the reference object 35 has a stepped thickness, the scattered ray component I0s(x, y ) is derived. For scattered ray components corresponding to thicknesses not present in the reference object 35, the characteristic deriving unit 27 may derive the relationship between the thickness of the reference object 35 and the scattered ray components by performing an interpolation calculation using scattered ray components of thicknesses present in the reference object 35. I0s(x,y) = I0o(x,y)-I0p(x,y) (9)

[0064] The characteristic deriving unit 27 calculates the scattering coefficient for the primary ray component I0p(x, y) for each thickness of the reference object 35. The ratio of the random line component I0s(x,y) (i.e., I0s(x,y) / I0p(x,y)) is Since the thickness of the reference object 35 varies stepwise, the STPR for thicknesses not present in the reference object 35 can be derived by interpolation using the STPR for thicknesses present in the reference object 35.

[0065] FIG. 8 is a diagram showing the relationship between the thickness of the reference object and the STPR. FIG. 8 shows the relationship between the thickness of the reference object 35 and the STPR when the tube voltage is 90 kV and the total filtration volume is 2.5 mmAl. The characteristic derivation unit 27 stores the derived scattered radiation model in the storage 13. The reference object 35 simulates the radiation characteristics of the human body. Therefore, the relationship between the thickness of the reference object and the STPR shown in FIG. 8 represents the relationship between the thickness of the subject H and the STPR.

[0066] The relationship between the thickness of the reference object and the STPR may be derived for each energy characteristic of the radiation that can be emitted by the radiation source 2 of the imaging apparatus 1A and stored in the storage 13.

[0067] In this embodiment, in addition to the radiation attenuation coefficient of soft tissue, the radiation attenuation coefficient of bone tissue is also used. Therefore, the characteristic deriving unit 27 also derives the radiation attenuation coefficient of bone tissue. The derivation of the radiation attenuation coefficient μ1Bone of bone tissue will be explained below. The radiation attenuation coefficient μBone of the tissue is also an example of the first characteristic in the present disclosure. In the description, the reference number "1" has been omitted for the sake of generality.

[0068] Here, assuming that soft tissue and bony tissue overlap on the radiation transmission path, and the thickness of the soft tissue is tSoft, the radiation attenuation coefficient can be derived as a function dependent on the thickness of the soft tissue. If the energy spectrum of the radiation emitted from the radiation source 2 is Sin(E) and the thickness of the soft tissue of the subject H is tSoft, the radiation dose Xout1(tSoft) after passing through the subject H in the absence of bony tissue can be calculated for each thickness t of the subject H using the radiation attenuation coefficient μSoft(E) of the soft tissue of the human body according to the following equation (10). Note that, like the above equation (5), equation (10) takes into account the radiation attenuation coefficient of objects (i.e., the tabletop 3A and grid 4) present between the subject H and the radiation detectors 5 and 6.

number

[0069] The radiation dose Xout2(t) in the presence of bone tissue is further derived by the following equation (11) using the radiation attenuation coefficient μBone(E) of the bone tissue.

number

[0070] The radiation attenuation coefficient of bone tissue is expressed as an exponential decay of the radiation dose due to bone tissue, using the radiation dose when there is no bone tissue as a reference, as shown in the following equation (12).

number

[0071] By solving equation (12) for μBone(t) as shown in the following equation (13), the relationship between the thickness t of the subject H and the radiation attenuation coefficient of bone tissue can be derived, where tBone is the thickness of bone tissue.

number

[0072] As described above, the PSF also varies depending on the distribution of the irradiation field in the imaging device 1A, the distribution of the composition of the subject H, the irradiation dose during imaging, the tube voltage, the imaging distance, and the characteristics of the radiation detectors 5 and 6. For this reason, the characteristic derivation unit 27 may experimentally determine in advance a PSF for each energy characteristic of the radiation used by the imaging device 1A in accordance with the irradiation field information, the subject information, the imaging conditions, and the like, and store the PSF in the storage 13.

[0073] In this embodiment, the characteristic deriving unit 27 derives the radiation attenuation coefficient μ2Soft of soft tissue, the radiation attenuation coefficient μ2Bone of bone tissue, STPR2, and PSF2 of the imaging device 1B in the same manner as described above. The radiation attenuation coefficients μ2Bone, STPR2 and PSF2 are stored in Storage 13. For the imaging device 1B, the radiation attenuation coefficient μ2Soft of soft tissue and the radiation attenuation coefficient μ2Soft of bone tissue are The linear attenuation coefficients μBone, STPR2, and PSF2 are related to the second imaging device in this disclosure. This is an example of the second property.

[0074] The subtraction unit 23 performs energy subtraction processing to derive a bone image Gb from which bones of the subject H are extracted and a soft tissue image Gs from which soft tissue is extracted from the first and second radiological images G1, G2 that have been subjected to scattered radiation removal processing. The bone image Gb and soft tissue image Gs are examples of the first bone image and first soft tissue image according to the present disclosure. Note that the first and second radiological images G1, G2 in the subsequent processing are processed radiological images from which scattered radiation components have been removed.

[0075] When deriving the bone image Gb, the subtraction unit 23 performs weighted subtraction between corresponding pixels of the first and second radiographic images G1 and G2 as shown in the following equation (14), thereby extracting the bones of the subject H contained in each of the radiographic images G1 and G2. In equation (14), α1 is a weighting coefficient, which is set to a value that enables extraction of the bones of subject H contained in each radiographic image G1, G2 using equation (14), based on the radiation attenuation coefficients of bone tissue and soft tissue. Gb(x, y)=G1(x, y)-α1×G2(x, y) (14)

[0076] On the other hand, when deriving a soft tissue image Gs, the subtraction unit 23 performs weighted subtraction between corresponding pixels of the first and second radiographic images G1, G2 as shown in the following equation (15) to generate a soft tissue image Gs from which the soft tissue of the subject H contained in each radiographic image G1, G2 has been extracted. In equation (15), α2 is a weighting coefficient, and is set to a value that enables extraction of the soft tissue of the subject H contained in each radiographic image G1, G2 using equation (15), based on the radiation attenuation coefficients of bony tissue and soft tissue. Gs(x, y)=G1(x, y)-α2×G2(x, y) (15)

[0077] Next, a description will be given of the processing performed by the conversion unit 24 and the derivation unit 25. Fig. 9 is a diagram schematically showing the processing performed by the conversion unit 24 and the derivation unit 25. First, a description will be given of the processing performed by the conversion unit 24.

[0078] The bone image Gb and soft tissue image Gs derived by the subtraction unit 23 are derived from the first and second radiographic images G1 and G2 acquired by the imaging device 1A, and therefore have contrast based on the first characteristic (i.e., radiation energy, radiation attenuation coefficient, STPR1, and PSF1) of the imaging device 1A. The conversion unit 24 converts the bone image Gb and soft tissue image Gs derived by the subtraction unit 23 so that they have contrast based on the second characteristic of the imaging device 1B.

[0079] The converter 24 converts the contrast of the bone image Gb using the following equation (16) to derive a converted bone image Gbt. The converter 24 converts the contrast of the soft tissue image Gs using the following equation (17) to derive a converted soft tissue image Gst. Note that β1 in equation (16) is derived from β1 = μ2Soft(T(x,y)) / μ1Soft(T(x,y)), and β2 in equation (17) is derived from β2 = μ2Bone(T(x,y)) / μ1Bone(T(x,y)). The body thickness distribution T(x,y) derived by the scattered radiation removal unit 22 is used. The converted bone image Gbt and the converted soft tissue image Gst are examples of a second bone image and a second soft tissue image according to the present disclosure. Gbt(x,y)=β1×Gb(x,y) (16) Gst(x,y)=β2×Gs(x,y) (17)

[0080] The derivation unit 25 derives a composite radiographic image Gc by adding corresponding pixels of the converted bone image Gbt and the converted soft tissue image Gst derived by the conversion unit 24. The composite radiographic image Gc is an example of a processed radiographic image in the present disclosure.

[0081] Here, the bone image Gb and the soft tissue image Gs are derived from the first and second radiographic images G1 and G2 from which scattered radiation components have been removed, and therefore the converted bone image Gbt, the converted soft tissue image Gst, and the composite radiographic image Gc do not contain scattered radiation components. Therefore, although the composite radiographic image Gc may be used as is for comparative interpretation with the first and second radiographic images G1 and G2 or the bone image Gb and the soft tissue image Gs, in this embodiment, scattered radiation components according to the second characteristic are added to the composite radiographic image Gc.

[0082] For this purpose, the derivation unit 25 uses the second characteristics, namely STPR2 and PSF2, derived by the characteristic derivation unit 27 for the imaging device 1B to derive a scattered ray image Isc representing the scattered ray components according to the second characteristics, according to the following equation (18). are scattered ray components corresponding to the scattered ray components contained in the radiographic image acquired when the subject H is imaged by the radiography device 1B. In equation (18), Gc(x, y) is the pixel value of each pixel of the composite radiographic image Gc. The body thickness distribution T(x, y) is derived by the scattered ray removal unit 22. Isc(x,y) = Gc(x,y)×STPR2(T(x,y))*PSF2(T(x,y)) (18)

[0083] The derivation unit 25 then adds up corresponding pixels of the combined radiation image Gc and the scattered radiation image Isc to derive the processed radiation image Gp.

[0084] The display control unit 26 displays the interpretation screen on the display 14. FIG. 10 is a diagram showing the interpretation screen. As shown in FIG. 10, the interpretation screen 40 has a first display area 41 that displays the radiographic image of the latest examination, and a second display area 42 that displays the past radiographic image. The first display area displays the processed radiographic image Gp acquired by the imaging device 1A and derived by the radiographic image processing device 10 according to this embodiment. The second display area displays the past radiographic image Gm acquired from the image storage system 9.

[0085] Next, the processing performed in this embodiment will be described. FIG. 11 is a flowchart showing the processing performed in this embodiment. It is assumed that the first and second radiographic images G1, G2 are acquired by the radiography device 1A and stored in the storage 13. It is also assumed that the first characteristic and the second characteristic are acquired by the characteristic derivation unit 27 and stored in the storage 13. When an instruction to start processing is input from the input device 15, the image acquisition unit 21 acquires the first and second radiographic images G1, G2 from the storage 13 (radiographic image acquisition; step ST1). Next, the scattered radiation removal unit 22 derives the body thickness distribution of the subject H from the first and second radiographic images G1, G2 (step ST2) and removes scattered radiation components from each of the first and second radiographic images G1, G2 (step ST3).

[0086] Next, the subtraction unit 23 derives a bone image Gb in which the bones of the subject H are extracted and a soft tissue image Gs in which the soft tissues are extracted from the first and second radiographic images G1, G2 from which the scattered radiation components have been removed (subtraction; step ST4).

[0087] Next, the conversion unit 24 converts the contrast of the bone image Gb and the soft tissue image Gs to derive a converted bone image Gbt and a converted soft tissue image Gst (conversion; step ST5). Next, the derivation unit 25 derives a composite radiographic image Gc by adding together the converted bone image Gbt and the converted soft tissue image Gst derived by the conversion unit 24 (step ST6). The derivation unit 25 also derives a scattered radiation image Isc representing the scattered radiation component according to the second characteristic (step ST7). The derivation unit 25 then derives a processed radiographic image Gp by adding together the composite radiographic image Gc and the scattered radiation image Isc (step ST8). Furthermore, the display control unit 26 displays the processed radiographic image Gp and the previous radiographic image Gm on the display 14 (image display; step ST9), and the process ends.

[0088] Here, the radiographic image acquired by the radiography apparatus 1A contains scattered ray components according to the characteristics of the energy of the radiation used in the radiography apparatus 1A, the top plate of the radiography table on which the subject H is placed in the radiography apparatus 1A, and an anti-scatter grid for removing scattered ray components contained in the radiation that has passed through the subject H in the radiography apparatus 1A. Therefore, the radiographic image acquired by the radiography apparatus 1A has contrast according to the first characteristic of the radiography apparatus 1A.

[0089] On the other hand, the radiation image acquired by the imaging device 1B is determined based on the energy of the radiation used in the imaging device 1B, the top plate of the imaging table on which the subject H is placed in the imaging device 1B, and the scattered radiation removal device for removing scattered radiation components contained in the radiation that has passed through the subject H in the imaging device 1B. It contains scattered radiation components according to the characteristics of the grid, etc. Therefore, the radiation image acquired by the imaging device 1B has contrast according to the second characteristic of the imaging device 1B.

[0090] In this embodiment, scattered radiation components are removed from the first and second radiographic images G1 and G2 acquired by the radiography device 1A, and a bone image Gb and a soft tissue image Gs are derived from the first and second radiographic images G1 and G2 from which scattered radiation components have been removed. Furthermore, the bone image Gb and the soft tissue image Gs are converted based on the first characteristic, the second characteristic, and the body thickness distribution T(x,y) of the subject H so that they have contrast based on the second characteristic of the radiography device 1B, and the converted bone image Gbt and the converted soft tissue image Gst are combined to derive a combined radiographic image Gc. Therefore, in this embodiment, the radiographic image acquired by the radiography device 1A can be converted to have a contrast corresponding to that of the radiographic image acquired by the radiography device 1B. This allows the contrast of the radiographic images used for comparative interpretation to match, even if they are acquired by different radiography devices 1A and 1B. Therefore, this embodiment enables accurate comparative interpretation.

[0091] Furthermore, by deriving a scattered radiation image representing the scattered radiation component according to the second characteristic based on the second characteristic and the body thickness distribution, and deriving a processed radiation image Gp using the derived scattered radiation image, it is possible to match the contrast of the radiation images acquired by the different imaging devices 1A and 1B, including the scattered radiation component.

[0092] In 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 this is not limiting. The first and second radiographic images G1, G2 may also 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 included in the first radiographic image G1 and the second radiographic image G2 may be shifted due to bodily movement of the subject H. For this reason, it is preferable to align the position of the subject in the first radiographic image G1 and the second radiographic image G2 before performing the process of this embodiment.

[0093] Furthermore, in the above embodiment, the bone disease prediction process is performed using radiographic images acquired by a system that captures the first and second radiographic images G1, G2 of the subject H using the first and second radiation detectors 5, 6. However, the technology of the present disclosure can also be applied to cases where the first and second radiographic images G1, G2 are acquired using stimulable phosphor sheets instead of radiation detectors. In this case, two stimulable phosphor sheets are placed one on top of the other and irradiated with radiation that has passed through the subject H, and radiographic image information of the subject H is stored and recorded on each stimulable phosphor sheet. The radiographic image information is then photoelectrically read from each stimulable phosphor sheet to acquire the first and second radiographic images G1, G2. The two-shot method may also be used when acquiring the first and second radiographic images G1, G2 using stimulable phosphor sheets.

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

[0095] In the above embodiment, the hardware structure of the processing unit that executes various processes, such as the image acquisition unit 21, the scattered radiation removal unit 22, the subtraction unit 23, the conversion unit 24, the derivation unit 25, the display control unit 26, and the characteristic derivation unit 27, is as follows: The following various processors can be used: As mentioned above, in addition to the CPU, which is a general-purpose processor that executes software (programs) and functions as various processing units, FPGAs (Field Programmable Gate Arrays) Programmable logic devices (PLDs), which are processors whose circuit configuration can be changed after manufacturing, such as programmable logic devices (PLCs), and application-specific integrated circuits (ASICs), which are processors with circuit configurations designed specifically to execute specific processes. It includes dedicated electrical circuits that are processors.

[0096] 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 (for example, 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.

[0097] Examples of configuring multiple processing units with a single processor include, first, a form in which one processor is configured with a combination of one or more CPUs and software, and this processor functions as multiple processing units, as typified by computers such as client and server. Second, a form in which a processor is used to realize the functions of an entire system including multiple processing units with a single IC (Integrated Circuit) chip, as typified by systems on chips (SoCs). In this way, various processing units are configured using one or more of the above-mentioned various processors as a hardware structure.

[0098] Furthermore, as the hardware structure of these various processors, more specifically, an electric circuit (Circuitry) combining circuit elements such as semiconductor elements can be used. . [Explanation of symbols]

[0099] 1A, 1B Imaging device 2 Radiation source 3. Photo stand 3A Top plate 3B Mounting part 4 Anti-scatter grid 5,6 Radiation detectors 7 Radiation Energy Conversion Filter 9. Image Storage System 10 Radiation image processing device 11 CPU 12 Radiation Image Processing Program 13. Storage 14 Display 15 Input Devices 16 memory 17 Network I / F 18 Bus 21 Image acquisition unit 22 Scattered radiation removal section 23 Subtraction section 24 Conversion unit 25 Derivation part 26 Display control unit 35 Reference object 40 Reading screen 41,42 Image display area Gb Bone image Gm Past radiological images Gp converted image Gs Soft tissue images

Claims

1. at least one processor; The processor: two radiological images having contrast based on a first characteristic of a first imaging device, the two radiological images being obtained by imaging a subject including soft tissue and bones using radiation with different energy distributions, the two radiological images being obtained by an imaging method in which the radiation with different energy distributions is irradiated twice onto the subject at different timings; deriving a body thickness distribution of the subject based on at least one of the two radiographic images; removing, from the two radiation images, a scattered radiation component that is included in the radiation that has passed through the subject and that is scattered by the subject, based on the first characteristic; deriving a first bone image representing the bone tissue of the subject and a first soft tissue image representing the soft tissue of the subject by performing weighted subtraction on the two radiographic images in which the subject has been aligned and the scattered radiation components have been removed; converting the first bone image and the first soft tissue image into a second bone image and a second soft tissue image having contrast based on the second characteristic, based on the first characteristic, a second characteristic related to a second imaging device different from the first imaging device, and the body thickness distribution; a radiation image processing device that derives a processed radiation image having contrast based on the second characteristic by adding the second bone image and the second soft tissue image;

2. the processor derives a scattered ray component according to the second characteristic based on the second characteristic and the body thickness distribution; The radiation image processing apparatus according to claim 1 , further comprising: a processing unit for deriving the processed radiation image using the derived scattered radiation component.

3. the first characteristics include energy of the radiation used in the first imaging device, a radiation attenuation coefficient corresponding to the body thickness distribution of an object interposed between the subject and a radiation detector that detects the radiation transmitted through the subject in the first imaging device, a ratio of the scattered ray components contained in the radiation transmitted through the subject that corresponds to the body thickness distribution, and a point spread function corresponding to the body thickness distribution; 3. The radiographic image processing device according to claim 1, wherein the second characteristics include: energy of the radiation used in the second imaging device; a radiation attenuation coefficient corresponding to the body thickness distribution of an object interposed between the subject and a radiation detector that detects the radiation that has passed through the subject in the second imaging device; a ratio of the scattered ray component contained in the radiation that has passed through the subject that corresponds to the body thickness distribution; and a point spread function corresponding to the body thickness distribution.

4. The radiation image processing apparatus according to claim 1 , wherein the processor displays the processed radiation image and a radiation image of the subject acquired by the second imaging device.

5. two radiological images having contrast based on a first characteristic of a first imaging device, the two radiological images being obtained by imaging a subject including soft tissue and bones using radiation with different energy distributions, the two radiological images being obtained by an imaging method in which the radiation with different energy distributions is irradiated twice onto the subject at different timings; deriving a body thickness distribution of the subject based on at least one of the two radiographic images; removing, from the two radiation images, a scattered radiation component that is included in the radiation that has passed through the subject and that is scattered by the subject, based on the first characteristic; deriving a first bone image representing the bone tissue of the subject and a first soft tissue image representing the soft tissue of the subject by performing weighted subtraction on the two radiographic images in which the subject has been aligned and the scattered radiation components have been removed; a second characteristic relating to a second image capture device different from the first characteristic and the first image capture device; and converting the first bone image and the first soft tissue image into a second bone image and a second soft tissue image having contrast based on the second characteristic based on the body thickness distribution; A radiological image processing method for deriving a processed radiological image having contrast based on the second characteristic by adding the second bone image and the second soft tissue image.

6. a step of acquiring two radiological images having contrast based on a first characteristic of a first imaging device, the two radiological images being obtained by imaging a subject including soft tissue and bones using radiation having different energy distributions with a first imaging device, the two radiological images being obtained by an imaging method in which the subject is irradiated twice with the radiation having different energy distributions at different times; deriving a body thickness distribution of the subject based on at least one of the two radiographic images; a step of removing, from the two radiation images, a scattered radiation component that is included in the radiation that has passed through the subject and that is scattered by the subject, based on the first characteristic; deriving a first bone image representing the bone tissue of the subject and a first soft tissue image representing the soft tissue of the subject by performing weighted subtraction on two radiographic images obtained by aligning the subject and removing the scattered radiation component; converting the first bone image and the first soft tissue image into a second bone image and a second soft tissue image having contrast based on the second characteristic, based on the first characteristic, a second characteristic related to a second imaging device different from the first imaging device, and the body thickness distribution; and a procedure of deriving a processed radiographic image having contrast based on the second characteristic by adding the second bone image and the second soft tissue image.

Images