Radiation imaging processing device, method, and program

Abstract

To provide a radiation image processing device, a method and a program which accurately derive bone density.SOLUTION: A processor derives a bone part image of a subject on the basis of the first radiation image and the second radiation image acquired by imaging the subject including a bone part with the radiation having different energy distributions, derives a fat ratio distribution or muscle ratio distribution of the subject on the basis of the first radiation image and the second radiation image, and derives bone density in a bone part region of the subject on the basis of the bone part image, the fat ratio distribution or the muscle ratio distribution.SELECTED DRAWING: Figure 3

Description

Technical Field

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

Background Art

[0002] In bone diseases such as osteoporosis, the DXA method (Dual X-ray Absorptiometry) is known as one of the typical bone mineral quantification methods used for diagnosing bone density. The DXA method utilizes the fact that the radiation incident on and transmitted through the human body undergoes attenuation characterized by the attenuation coefficient μ (cm2 / g), density (g / cm3), and thickness t (cm) of the substances (such as bone) constituting the human body, and calculates the bone density from the pixel values of the radiation images obtained by imaging with radiation of two types of energies.

[0003] In addition, various methods for evaluating bone density using radiation images acquired by imaging a subject have been proposed. For example, in Patent Document 1, a bone part image in which a bone part of a subject is extracted is generated from a plurality of radiation images acquired by radiation having different energy distributions that have passed through the subject, the density of the region without bone in the radiation image is calculated as a correction value representing the fat content, and the pixel values of the bone part image are corrected by the correction value, thereby proposing a method for more accurately calculating the density of the bone part region and further the bone density.

Prior Art Documents

Patent Documents

[0004]

Patent Document 1

Summary of the Invention

Problems to be Solved by the Invention

[0005] Incidentally, the attenuation coefficient of radiation changes depending on the proportion of fat contained in the soft tissue of the subject, i.e., the fat percentage. Therefore, even with the same body thickness, the effect of radiation beam hardening changes depending on the fat percentage. Specifically, the higher the fat percentage, the lower the energy of the transmitted radiation, so the contrast of the bone region in the acquired radiation image increases. As a result, the higher the fat percentage, the higher the bone density value derived. Therefore, even if the bone density is the same, there is a problem that errors in bone density will occur depending on the amount of fat percentage.

[0006] This disclosure is made in view of the above circumstances and aims to enable the accurate derivation of bone density. [Means for solving the problem]

[0007] The radiation image processing device according to this disclosure comprises at least one processor, The processor derives a bone image of the subject based on a first and second radiation image obtained by photographing the subject, including the bone, with radiation that has a different energy distribution. Based on the first and second radiographic images, the fat percentage distribution or muscle percentage distribution of the subject is derived. Based on bone images and fat percentage distribution or muscle percentage distribution, bone density in the bone region of the subject is derived.

[0008] Furthermore, in the radiation image processing apparatus according to this disclosure, the processor derives the thickness distribution of the subject based on the first radiation image and the second radiation image. By referring to the relationship between fat percentage or muscle percentage and the body thickness distribution of the subject to a conversion coefficient for converting it to the standard soft tissue thickness distribution, a conversion coefficient corresponding to the fat percentage distribution or muscle percentage distribution is obtained. By transforming the thickness distribution of the subject using a conversion coefficient, the thickness distribution of the subject is converted to the thickness distribution of standard soft tissue. By referring to the relationship between body thickness distribution and correction coefficients for correcting pixel values ​​in bone images to bone density, a correction coefficient corresponding to the standard soft tissue thickness distribution is obtained. The bone density may be derived by correcting the bone image using a correction coefficient.

[0009] "Standard soft tissue" refers to soft tissue having a reduction coefficient and fat percentage or muscle percentage equivalent to that of standard human soft tissue. A material equivalent to that of standard human soft tissue can be used as the standard soft tissue. This material can be acrylic or urethane, for example. Furthermore, the correction coefficient for correcting the pixel values ​​of bone images to bone density is derived from a material equivalent to standard soft tissue and can be derived using radiographic images of phantoms of various thicknesses.

[0010] Furthermore, in the radiation image processing apparatus according to this disclosure, the processor removes scattered radiation components based on radiation scattered by the subject from the first radiation image and the second radiation image. The first and second radiographic images, from which scattered radiation components have been removed, may be used to derive the fat percentage distribution or muscle percentage distribution, and the bone density.

[0011] Furthermore, in the radiation image processing apparatus according to this disclosure, the processor may display bone density on a display.

[0012] The radiation image processing method disclosed herein derives a bone image of a subject based on a first radiation image and a second radiation image obtained by photographing the subject, including the bone, with radiation having different energy distributions. Based on the first and second radiographic images, the fat percentage distribution or muscle percentage distribution of the subject is derived. Based on bone images and fat percentage distribution or muscle percentage distribution, bone density in the bone region of the subject is derived.

[0013] The radiographic image processing program disclosed herein includes a procedure for deriving a bone image of a subject based on a first radiographic image and a second radiographic image obtained by photographing the subject, including the bone, with radiation having different energy distributions. Based on the first radiographic image and the second radiographic image, a procedure for deriving the fat percentage distribution or muscle percentage distribution of a subject, and Based on the bone image and the fat percentage distribution or muscle percentage distribution, causing a computer to execute a procedure for deriving the bone density in the bone region of the subject.

Effect of the Invention

[0014] According to the present disclosure, bone density can be accurately derived.

Brief Description of the Drawings

[0015] [Figure 1] Schematic block diagram showing the configuration of a radiographic imaging system to which a radiographic image processing apparatus according to an embodiment of the present disclosure is applied [Figure 2] Diagram showing the schematic configuration of a radiographic image processing apparatus according to an embodiment of the present disclosure [Figure 3] Diagram showing the functional configuration of a radiographic image processing apparatus according to an embodiment of the present disclosure [Figure 4] Diagram showing a soft tissue image [Figure 5] Diagram showing the relationship between the contrast of the bone part and the soft part with respect to the body thickness [Figure 6] Diagram showing a look-up table for obtaining a correction coefficient [Figure 7] Diagram showing the relationship between the fat percentage and the attenuation coefficient [Figure 8] Diagram showing the relationship between the body thickness and the attenuation coefficient [Figure 9] Diagram showing the relationship between the fat percentage and the body thickness conversion coefficient <000Embodiments of this disclosure will be described below with reference to the drawings. Figure 1 is a schematic block diagram showing the configuration of a radiation imaging system to which a radiation imaging processing device according to an embodiment of this disclosure is applied. As shown in Figure 1, the radiation imaging system according to this embodiment comprises a shooting device 1 and a radiation imaging processing device 10 according to this embodiment.

[0017] The imaging device 1 is an imaging device for performing energy subtraction using the so-called one-shot method, in which radiation such as X-rays emitted from the radiation source 3 and transmitted through the subject H is irradiated onto the first radiation detector 5 and the second radiation detector 6 with varying energies. During imaging, as shown in Figure 1, the first radiation detector 5, a radiation energy conversion filter 7 made of a copper plate or the like, and the second radiation detector 6 are arranged in 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 and the radiation energy conversion filter 7 are in close contact.

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

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

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

[0021] The storage 13 is implemented using an HDD (Hard Disk Drive), an SSD (Solid State Drive), and flash memory, etc. The storage 13, as a storage medium, stores the radiation image processing program 12 installed in the radiation image processing device 10. The CPU 11 reads the radiation image processing program 12 from the storage 13, expands it into memory 16, and executes the expanded radiation image processing program 12.

[0022] The radiation image processing program 12 is stored in a memory device of a server computer connected to the network, or in network storage, in a state that allows external access, and is downloaded and installed on the computers comprising the radiation image processing device 10 upon request. Alternatively, it is recorded on a recording medium such as a DVD (Digital Versatile Disc) or CD-ROM (Compact Disc Read Only Memory) and distributed, and then installed from that recording medium on the computers comprising the radiation image processing device 10.

[0023] Next, the functional configuration of the radiation image processing apparatus according to this embodiment will be described. Figure 3 is a diagram showing the functional configuration of the radiation image processing apparatus according to this embodiment. As shown in Figure 3, the radiation image processing apparatus 10 includes an image acquisition unit 21, a scattered radiation removal unit 22, a bone image output unit 23, a first output unit 24, a second output unit 25, and a display control unit 26. The CPU 11 then functions as the image acquisition unit 21, the scattered radiation removal unit 22, the bone image output unit 23, the first output unit 24, the second output unit 25, and the display control unit 26 by executing the radiation image processing program 12.

[0024] The image acquisition unit 21 causes the imaging device 1 to perform energy subtraction imaging of the subject H, thereby acquiring the first radiation image G1 and the second radiation image G2 of the subject H from the first and second radiation detectors 5 and 6. When acquiring the first radiation image G1 and the second radiation image G2, imaging conditions such as the imaging dose, beam 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 and 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 a scatter removal grid are set.

[0025] SOD and SID are used to calculate body thickness distribution, as described later. For SOD, it is preferable to acquire it using, for example, a TOF (Time Of Flight) camera. For SID, it is preferable to acquire it using, for example, a potentiometer, ultrasonic rangefinder, or laser rangefinder.

[0026] The shooting conditions can be set by the operator through input from the input device 15.

[0027] The scattered radiation removal unit 22 removes scattered radiation components from the first radiation image G1 and the second radiation image G2 acquired by the image acquisition unit 21. The removal of scattered radiation components will be described below. Any method can be used to remove scattered radiation components, such as the method described in Japanese Patent Application Publication No. 2015-043959. The scattered radiation removal process using the method described in Japanese Patent Application Publication No. 2015-043959 will be described below. In the following description, the first and second radiation images from which scattered radiation components have been removed will also be referred to as G1 and G2, respectively.

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

[0029] Next, the scattered radiation removal unit 22 derives an estimated primary radiation image Ip(x,y) based on the virtual model, which is an estimated primary radiation image obtained by imaging the virtual model, and an estimated scattered radiation image Is(x,y) based on the scattered radiation image obtained by imaging the virtual model, as shown in equations (1) and (2) below. Furthermore, as shown in equation (3) below, the scattered radiation removal unit 22 derives an estimated image Im(x,y) by combining the estimated primary radiation image Ip(x,y) and the estimated scattered radiation image Is(x,y) as an estimated first radiation image G1 obtained by imaging the subject H. Ip(x,y) = Io(x,y)×exp(-μSoftT(x,y))×T(x,y)) (1) Is(x,y) = Io(x,y)×STPR(T(x,y))*PSF(T(x,y)) (2) Im(x,y) = Is(x,y)+Ip(x,y) (3)

[0030] Here, (x,y) is the coordinate of the pixel position of the first radiation image G1, Io(x,y) is the pixel value of the first radiation image G1 at pixel position (x,y), Ip(x,y) is the primary component at pixel position (x,y), and Is(x,y) is the scattered component at pixel position (x,y). Note that 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).

[0031] Furthermore, μSoft(T(x,y)) in equation (1) is an attenuation coefficient corresponding to the body thickness distribution (x,y) of the soft tissue of the human body at the pixel position (x,y). μSoft(T(x,y)) can be determined in advance experimentally or by simulation and stored in storage 13. In this embodiment, the attenuation coefficient μSoft(T(x,y)) is the attenuation coefficient of a standard soft tissue described later. Furthermore, STPR(T(x,y)) in equation (2) is the scatter-to-primary ratio of the scattered dose to the primary dose contained in the radiation after it has passed through an object H having a body thickness distribution T(x,y). STPR(T(x,y)) can also be determined in advance experimentally or by simulation and stored in storage 13.

[0032] Furthermore, PSF(T(x,y)) in equation (2) is a point spread function that represents the distribution of scattered rays spreading from a single pixel according to the body thickness distribution T(x,y), and is defined according to the energy characteristics of the radiation. Also, * is an operator that indicates a convolution operation. PSF also changes depending on the distribution of the irradiation field in the imaging device 1, the composition distribution 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, PSF should be experimentally determined in advance for each energy characteristic of the radiation used by the imaging device 1, according to the irradiation field information, subject information, and imaging conditions, and stored in storage 13.

[0033] Next, the scattered radiation removal unit 22 modifies the initial body thickness distribution Ts(x,y) of the virtual model so that the difference between the estimated image Im and the first radiation image G1 is small. The scattered radiation removal unit 22 updates the body thickness distribution T(x,y), scattered radiation component Is(x,y), and primary radiation component Ip(x,y) by repeatedly deriving the body thickness distribution T(x,y), scattered radiation component Is(x,y), and primary radiation component Ip(x,y) until the difference between the estimated image Im and the first radiation image G1 satisfies a predetermined termination condition. When the termination condition is met, the scattered radiation removal unit 22 subtracts the scattered radiation component Is(x,y) derived by equation (3) from the first radiation image G1. This removes the scattered radiation component contained in the first radiation image G1. The body thickness distribution T(x,y) derived when the termination condition is met is used when deriving bone density, which will be described later.

[0034] Meanwhile, the scattered radiation removal unit 22 performs the scattered radiation removal process on the second radiation image G2 in the same manner as on the first radiation image G1. In the following description, the reference numerals G1 and G2 will also be used for the first and second radiation images from which the scattered radiation component has been removed.

[0035] The bone image extraction unit 23 generates a bone image Gb representing the bone region of subject H from the first radiographic image G1 and the second radiographic image G2, from which the scattered radiation component has been removed by the scattered radiation removal unit 22. Figure 4 shows an example of a bone image Gb generated by the bone image extraction unit 23. Generally, the lumbar spine or femur is used to measure bone density. Therefore, the bone image Gb shown in Figure 4 is the bone image Gb generated from the first radiographic image G1 and the second radiographic image G2, which were obtained by taking images that include part of the femur and lumbar spine of subject H.

[0036] The bone image extraction unit 23 generates a bone image Gb in which only the bone portion of the subject H included in the first and second radiation images G1 and G2 is extracted by performing weighted subtraction between corresponding pixels on the first and second radiation images G1 and G2, as shown in equation (4) below. In equation (4) below, α1 is the weighting coefficient, and x and y are the coordinates of each pixel in the bone image Gb. Gb(x,y)=G1(x,y)-α1×G2(x,y) (4)

[0037] The first derivation unit 24 derives the fat percentage distribution of the subject H, as will be described later. Fat percentage refers to the proportion of fat in the soft tissues of the human body.

[0038] The second derivation unit 25 derives a bone density image representing the bone density in the bone region of the subject H based on the bone image Gb generated by the bone image derivation unit 23. The second derivation unit 25 derives the bone density image B by deriving the bone density B(x,y) for each pixel (x,y) of the bone image Gb. Note that the bone density B(x,y) may be derived for all bones included in the bone image Gb, or the bone density B(x,y) may be derived only for predetermined bones. As an example, in this embodiment, the predetermined bone included in the bone image Gb is the femur, and the second derivation unit 25 derives the bone density B(x,y) for each pixel of the bone region corresponding to the femur, and does not derive the bone density B(x,y) for bone regions corresponding to other bones included in the bone image Gb. Hereinafter, the bone region from which the second derivation unit 25 derives bone density B(x,y) refers to the bone region corresponding to the femur.

[0039] Specifically, the second derivation unit 25 derives the bone density B(x,y) corresponding to each pixel by converting the pixel values ​​Gb(x,y) of the bone region in the bone image Gb to the pixel values ​​of the bone image acquired under standard imaging conditions. More specifically, the second derivation unit 25 derives the bone density B(x,y) for each pixel by correcting the pixel values ​​Gb(x,y) of the bone image Gb using a correction coefficient obtained from a lookup table described later.

[0040] Here, the higher the tube voltage at radiation source 3 and the higher the energy of the radiation emitted from radiation source 3, the smaller the contrast (i.e., the difference in pixel values) between soft tissue and bone in the radiation image. Also, during the process of radiation passing through subject H, the low-energy component of the radiation is absorbed by subject H, resulting in beam hardening, which increases the energy of the radiation. The increase in radiation energy due to beam hardening increases with the thickness of subject H.

[0041] Figure 5 shows the relationship between the contrast between bone and soft tissue and the thickness of subject H. Figure 5 shows the relationship between the contrast between bone and soft tissue and the thickness of subject H at three tube voltages: 80kV, 90kV, and 100kV. As shown in Figure 5, the higher the tube voltage, the lower the contrast. Furthermore, beyond a certain thickness of subject H, the greater the thickness, the lower the contrast. Note that the larger the pixel value Gb(x,y) in the bone region of the bone image Gb, the greater the contrast between bone and soft tissue. Therefore, the relationship shown in Figure 5 shifts towards higher contrast as the pixel value Gb(x,y) in the bone region of the bone image Gb increases.

[0042] In this embodiment, a lookup table (not shown) for obtaining correction coefficients to compensate for the difference in contrast in the bone image Gb depending on the tube voltage during acquisition, and for the decrease in contrast due to the effect of beam hardening, is stored in storage 13. The correction coefficient is a coefficient for correcting each pixel value Gb(x,y) of the bone image Gb.

[0043] Figure 6 shows an example of a lookup table stored in storage 13. In Figure 6, a lookup table LUT1 is shown as an example, with the reference imaging condition set to a tube voltage of 90kV. As shown in Figure 6, in lookup table LUT1, a larger correction coefficient is set as the tube voltage increases and the thickness of the subject increases. In the example shown in Figure 6, since the reference imaging condition is a tube voltage of 90kV, the correction coefficient is 1 when the tube voltage is 90kV and the body thickness is 0. Note that although lookup table LUT1 is shown in two dimensions in Figure 6, the correction coefficient differs depending on the pixel value of the bone region. Therefore, lookup table LUT1 is actually a three-dimensional table with an axis representing the pixel value of the bone region.

[0044] Here, the correction coefficient is derived by taking radiographic images using a phantom that simulates the human body, containing materials corresponding to the soft tissues and bone tissues of the human body. For example, acrylic or urethane can be used as the material corresponding to the soft tissues. For example, hydroxyapatite can be used as the material corresponding to the bone tissues of the human body. The radiation attenuation coefficient and fat content of the material corresponding to the soft tissues in the phantom are predetermined values ​​depending on the material. As mentioned above, the greater the thickness of the subject's body, the smaller the contrast between bone and soft tissue. As a result, even if the bone density is the same, the bone density derived from the radiographic image will be smaller as the body thickness increases, given the same body thickness.

[0045] Therefore, by imaging phantoms of various thicknesses, pixel values ​​for the bone region corresponding to the body thickness are derived from the acquired radiographic images, and a correction coefficient corresponding to the body thickness is derived so that the derived pixel values ​​correspond to the same bone density. Furthermore, the correction coefficient is derived according to various tube voltages. This allows the lookup table LUT1 shown in Figure 6 to be derived.

[0046] The second derivation unit 25 extracts a pixel-by-pixel correction coefficient C0(x,y) from the lookup table LUT1, which corresponds to the imaging conditions including the body thickness distribution T(x,y) of the subject H and the tube voltage setting value stored in the storage 13. Then, as shown in equation (5) below, the second derivation unit 25 derives the bone density B(x,y) for each pixel of the bone region in the bone region image Gb by multiplying it by the correction coefficient C0(x,y). This derives a bone density image B in which the bone density B(x,y) is the pixel value. The bone density B(x,y) thus derived represents the pixel value of the bone region in the bone region included in the radiographic image obtained by imaging the subject H with a tube voltage of 90kV, which is the standard imaging condition, and from which the effects of beam hardening have been removed. In this embodiment, the unit of bone density is g / cm 2 Let it be so. B(x,y) = C0(x,y) × Gb(x,y) (5)

[0047] Incidentally, the radiation attenuation coefficient changes depending on the fat content of the soft tissue. Therefore, even with the same body thickness, the effect of radiation beam hardening changes depending on the fat content. Specifically, the higher the fat content, the larger the proportion of relatively low-energy radiation in the transmitted radiation, resulting in a higher contrast in the bone region of the acquired radiation image. As a result, the higher the fat content, the larger the pixel value Gb(x,y) of the bone image Gb. Furthermore, the lookup table LUT1 shown in Figure 6 is derived using a phantom. The fat content of the material corresponding to the soft tissue in the phantom is a constant value. Therefore, when the pixel value Gb(x,y) of the bone image Gb is corrected using the lookup table shown in Figure 6, a bone density different from the actual bone density is derived depending on the magnitude of the fat content.

[0048] Therefore, in this embodiment, it is conceivable to derive the fat percentage for each pixel from the first and second radiographic images G1 and G2, pre-derive a correction coefficient that reduces the pixel value Gb(x,y) as the fat percentage increases, and then use the correction coefficient to correct the pixel value Gb(x,y) of the bone image Gb.

[0049] However, if the lookup table LUT1 shown in Figure 6 is further modified to include fat percentages, the computational load for obtaining correction coefficients increases. Therefore, in this embodiment, changes in fat percentage are converted into changes in the thickness of a standard soft tissue (hereinafter referred to as standard soft tissue), and bone density B(x,y) is derived from the pixel values ​​Gb(x,y) of the bone image Gb using only the thickness, without using fat percentage. In this embodiment, the soft tissue used in the phantom used to derive the correction coefficients in the lookup table LUT1 is used as the standard soft tissue. Therefore, the standard soft tissue has predetermined attenuation coefficients and fat percentages according to the material corresponding to the soft tissue of the phantom. The conversion of changes in fat percentage to changes in the thickness of the standard soft tissue will be explained below.

[0050] The fat percentage is derived by the first derivation unit 24. First, the first derivation unit 24 generates a soft tissue image Gs from the first and second radiation images G1 and G2, respectively, by extracting the soft tissue of the subject H contained in each radiation image G1 and G2, using the following equation (6). In equation (6), α2 is the weighting coefficient. Gs(x, y)=G1(x, y)-α2×G2(x, y) (6)

[0051] Furthermore, the first derivation unit 24 derives the fat percentage rf(x,y) at each pixel position (x,y) in the soft tissue image Gs using the following equation (7). The fat percentage rf(x,y) at each pixel position (x,y) in the soft tissue image Gs is an example of the fat percentage distribution according to this disclosure. In equation (7), αm is a weighting coefficient corresponding to the attenuation coefficient of muscle tissue, and αf is a weighting coefficient corresponding to the attenuation coefficient of adipose tissue. Also, Δ(x,y) represents the density difference distribution. The density difference distribution is the distribution of density changes on the image as seen from the density obtained when radiation reaches the first radiation detector 5 and the second radiation detector 6 without passing through the subject H. The distribution of density changes on the image is calculated by subtracting the density of each pixel in the subject H region from the density in the pass-through region obtained when radiation in the soft tissue image Gs directly irradiates the first radiation detector 5 and the second radiation detector 6. rf(x,y)={αm-Δ(x,y) / T(x,y)} / (αm-αf) (7)

[0052] The second derivation unit 25 determines the body thickness at which the transmitted radiation dose is the same as the attenuation coefficient for standard soft tissue, and converts the change in fat percentage into a change in body thickness. Here, since the attenuation coefficient for the soft tissue of subject H depends on the fat percentage rf and the body thickness T, the attenuation coefficient for soft tissue is expressed as μ(rf,T). When the body thickness T is constant, the relationship between fat percentage and the attenuation coefficient is shown in Figure 7. Also, when the fat percentage is constant, the relationship between body thickness and the attenuation coefficient is shown in Figure 8. If the incident X-ray dose is I0 and the X-ray dose transmitted through subject H is Ih, then Ih is expressed by the following equation (8) using the attenuation coefficient μ(rf,T). Note that in equation (8) and equations (9) and (10) described later, values ​​are derived for each pixel, but (x,y) is omitted. Ih = I0·exp(-μ(rf,T)·T) (8)

[0053] On the other hand, in the radiographic image obtained by photographing the phantom, let rf1(x,y) be the fat content of each pixel in the standard soft tissue (i.e., the region of material corresponding to soft tissue), let T1(x,y) be the thickness of the standard soft tissue, and let μ(rf1,T1) be the attenuation coefficient of the standard soft tissue. In this case, if the incident X-ray dose is I0 and the X-ray dose transmitted through the region of material corresponding to soft tissue in the phantom is I1, then I1 is expressed by the following equation (9). I1=I0·exp(-μ(rf1,T1)·T1) (9)

[0054] Since the fat content of the standard soft tissue is constant, the attenuation coefficient μ(rf1,T1) depends only on the body thickness T1. Therefore, the attenuation coefficient of the standard soft tissue can be expressed as μ(T1).

[0055] When the transmitted dose to subject H is equal to that to subject with standard soft tissue, Ih = I1. Therefore, from equations (8) and (9), we derive T1 such that μ(rf,T)·T = μ(T1)·T1. Here, if we set ΔI = I1 - Ih, ΔI=μ(T1)·T1-μ(rf,T)·T (10) Therefore, in equation (10), when T1 changes, μ(T1) also changes, so (ΔI) 2 We derive the value of T1 that minimizes (ΔI). Specifically, we gradually increase or decrease T1 from an initial value (e.g., T) to find the value of T1 that minimizes (ΔI). 2 By changing T1 so that (ΔI) becomes smaller, 2 We derive the value of T1 that minimizes [the expression].

[0056] Then, by determining the relationship between the body fat percentage and the derived ratio of T1 to the body thickness T of subject H (T1 / T), it is possible to convert the body thickness of subject H with a certain body fat percentage to the standard soft tissue thickness T1. If the derived ratio of T1 to the body thickness T of subject H (T1 / T) is taken as the body thickness conversion coefficient K(x,y), the relationship between the body fat percentage and the body thickness conversion coefficient is shown in Figure 9.

[0057] In this embodiment, the relationship between body fat percentage and body thickness conversion coefficient shown in Figure 9 is stored in storage 13 as a lookup table LUT2. The first derivation unit 24 derives the soft tissue image Gs of the subject H using equation (6) and derives the body fat percentage rf(x,y) for each pixel of the soft tissue image Gs using equation (7). The second derivation unit 25 then refers to LUT2 to derive a body thickness conversion coefficient K(x,y) corresponding to the body fat percentage rf(x,y) and converts the body thickness T(x,y) of the subject H to the standard soft tissue body thickness T1(x,y). The second derivation unit 25 then refers to the lookup table LUT1 shown in Figure 6 to derive a correction coefficient C0(x,y) corresponding to the derived body thickness T1(x,y) and derives the bone density B(x,y) using equation (5).

[0058] In the scattered radiation removal section 22, the body thickness distribution T(x,y) and scattered radiation components are derived using the above equations (1) to (3). The attenuation coefficient μSoft(T(x,y)) used in this process is the attenuation coefficient for a standard soft section.

[0059] The display control unit 26 displays the bone density estimated by the second derivation unit 25 on the display 14. Figure 10 shows the bone density display screen. As shown in Figure 10, the display screen 30 has an image display area 31. The image display area 31 displays a bone density image Gd representing the bone density distribution in the radiographic image of the subject H. In the bone density image Gd, patterns are applied to the bone regions according to the magnitude of the derived bone density. In Figure 10, for the sake of simplicity, patterns representing bone density are applied only to the femur. Below the image display area 31, a reference 32 is displayed indicating the magnitude of the bone density for the applied patterns. The operator can easily recognize the patient's bone density by reading the bone density image Gd while referring to the reference 32. Alternatively, different colors may be applied to the bone density image Gd according to the bone density instead of patterns.

[0060] Next, the processing performed in this embodiment will be described. Figure 11 is a flowchart showing the processing performed in this embodiment. First, the image acquisition unit 21 causes the imaging device 1 to perform energy subtraction imaging of the subject H, thereby acquiring first and second radiation images G1 and G2 (radiation image acquisition; step ST1). Next, the scattered radiation removal unit 22 derives the body thickness distribution T(x,y) of the subject H (step ST2) and removes the scattered radiation component from the first and second radiation images G1 and G2 (step ST3). Steps ST2 and ST3 are performed in parallel. Furthermore, the bone image extraction unit 23 derives a bone image Gb from which the bone portion of the subject H has been extracted from the first and second radiation images G1 and G2 from which the scattered radiation component has been removed (step ST4).

[0061] Next, the first derivation unit 24 derives a soft tissue image Gs from the first and second radiation images G1 and G2, from which scattered radiation components have been removed, using the above formula (6) to extract the soft tissue of subject H, and further derives the fat percentage rf(x,y) of subject H from the soft tissue image Gs (step ST5). Then, the second derivation unit 25 derives a body thickness conversion coefficient K(x,y) corresponding to the fat percentage rf(x,y) by referring to the lookup table LUT2 shown in Figure 9 (step ST6). Furthermore, the second derivation unit 25 converts the body thickness of subject H to the body thickness of standard soft tissue using the body thickness conversion coefficient K(x,y) (body thickness conversion; step ST7). Next, the second derivation unit 25 refers to the lookup table LUT1 shown in Figure 6 to obtain a correction coefficient C0(x,y) corresponding to the body thickness of the converted subject H (step ST8), derives bone density B(x,y) using the obtained correction coefficient C0(x,y) (step ST9), and terminates the process.

[0062] Thus, in this embodiment, the bone density in the bone region of the subject H is derived based on the bone image and fat percentage distribution of the subject H. Therefore, since the derived bone density takes fat percentage into account, the bone density can be derived with high accuracy according to this embodiment.

[0063] In particular, in this embodiment, since the body thickness distribution T(x,y) of the subject H is converted to the standard soft tissue thickness distribution T1(x,y) according to the fat percentage distribution of the subject H, it is no longer necessary to prepare a lookup table LUT1 for determining the correction coefficient C0(x,y) for converting the bone image Gb pixel value Gb(x,y) to bone density, according to the fat percentage. Therefore, the computational load when deriving bone density can be reduced.

[0064] Furthermore, by using the first and second radiographic images G1 and G2, from which scattered radiation components have been removed, the fat percentage distribution and bone density can be derived with greater accuracy.

[0065] Furthermore, as shown in Figure 1, bone density can be derived using a radiographic image obtained by an imaging method (hereinafter referred to as simple radiography) that irradiates the subject H with radiation and obtains a two-dimensional image which is a transmission image of the subject H. Therefore, bone density can be derived more easily compared to when bone density is derived using a device that uses the DXA method or a CT (Computed Tomography) device. For this reason, the radiographic image processing device according to this embodiment can be applied to continuous use such as health checkups or monitoring the progress of treatment.

[0066] In the above embodiment, bone density of the femur near the hip joint is derived and displayed, but the bones covered are not limited to the femur. The technology of this disclosure can also be applied to estimate information related to bone density of any bone, such as the vertebrae, femur and tibia near the knee joint, lumbar vertebrae, calcaneus, and metacarpals.

[0067] Furthermore, although the above embodiment displays a bone density image on the display screen 30, it is not limited to this. The second derivation unit 25 may derive a representative value of the bone density of the target area (for example, around the joint of the femur) and display the derived representative value of the bone density. As the representative value, the mean, median, maximum, and minimum values ​​of bone density can be used.

[0068] Furthermore, in the above embodiment, bone density is derived based on bone image Gb and fat percentage distribution, but bone density may also be derived based on bone image Gb and muscle percentage distribution. Here, the soft tissues of the human body include muscle tissue, adipose tissue, internal organs, blood, and water. Since muscle tissue, internal organs, blood, and water have similar radiation absorption characteristics, non-adipose tissue, including muscle tissue, internal organs, blood, and water, can be treated as muscle tissue. For this reason, the muscle percentage rm(x,y) can be derived from the fat percentage rf(x,y) derived by equation (7) using the following equation (11). The muscle percentage rm(x,y) is an example of the muscle percentage distribution of this disclosure. Note that in the thigh of the human body, the main components of soft tissue are muscle tissue and adipose tissue, so the muscle percentage may be derived by approximating non-adipose tissue as muscle tissue. rm(x,y)=1-rf(x,y) (11)

[0069] Furthermore, when the body thickness T is constant, the relationship between muscle percentage and the attenuation coefficient is monotonically increasing, unlike the relationship between fat percentage and the attenuation coefficient shown in Figure 7, which is monotonically decreasing. Therefore, the relationship between muscle percentage and the body thickness conversion coefficient is also monotonically increasing, unlike the relationship between fat percentage and the body thickness conversion coefficient shown in Figure 9, which is monotonically decreasing. Accordingly, the relationship between muscle percentage and the body thickness conversion coefficient can be derived and stored in storage 13, and by referring to this relationship, a body thickness conversion coefficient corresponding to the muscle percentage rm(x,y) can be derived, and the body thickness T(x,y) of the subject H can be converted to the standard soft tissue body thickness T1(x,y) using the derived body thickness conversion coefficient. Thus, similar to the above embodiment using fat percentage, a correction coefficient C0(x,y) corresponding to the derived body thickness T1(x,y) can be derived, and bone density B(x,y) can be derived using the above equation (5).

[0070] Furthermore, in the above embodiment, when performing energy subtraction processing to derive bone density, the first and second radiographic images G1 and G2 are acquired by the one-shot method, but the method is not limited to this. The first and second radiographic images G1 and G2 may also be acquired by the 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 subject H included in the first radiographic image G1 and the second radiographic image G2 may shift due to the movement of subject H. For this reason, it is preferable to perform alignment of the subject in the first radiographic image G1 and the second radiographic image G2 before performing the processing of this embodiment.

[0071] Furthermore, in the above embodiment, bone density is derived using radiation images acquired in a system that photographs subject H using first and second radiation detectors 5 and 6. However, the technology of this disclosure can also be applied when acquiring first and second radiation images G1 and G2 using a accumulative phosphor sheet instead of radiation detectors. In this case, two accumulative phosphor sheets are stacked and radiation transmitted through subject H is irradiated onto them to accumulate and record radiation image information of subject H on each accumulative phosphor sheet. The first and second radiation images G1 and G2 can then be acquired by photoelectrically reading the radiation image information from each accumulative phosphor sheet. Note that even when acquiring first and second radiation images G1 and G2 using accumulative phosphor sheets, the two-shot method may be used.

[0072] Furthermore, the radiation in the above embodiment is not particularly limited, and in addition to X-rays, alpha rays or gamma rays can be used.

[0073] Furthermore, in the above embodiment, the hardware structure of the Processing Unit that performs various processes such as the image acquisition unit 21, the scattered radiation removal unit 22, the bone image output unit 23, the first output unit 24, the second output unit 25, and the display control unit 26 of the radiation image processing device 10 can be any of the following types of processors. As mentioned above, these types of 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 manufacturing, such as an FPGA (Field Programmable Gate Array), and a dedicated electrical circuit, which is a processor with a circuit configuration specifically designed to perform a particular process, such as an ASIC (Application Specific Integrated Circuit).

[0074] A single processing unit may consist of one of these various processors, or it may consist of 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). Alternatively, multiple processing units may be composed of a single processor.

[0075] Examples of configuring multiple processing units with a single processor include, firstly, a configuration where one or more CPUs and software combine to form a single processor, as exemplified by client and server computers, and this processor functions as multiple processing units. Secondly, a configuration using a processor that realizes the functions of the entire system, including multiple processing units, on a single IC (Integrated Circuit) chip, as exemplified by System-on-a-Chip (SoC). Thus, various processing units are configured, in terms of hardware structure, using one or more of the above-mentioned processors.

[0076] Furthermore, the hardware structure of these various processors can more specifically utilize electrical circuits (Circuitry) that combine circuit elements such as semiconductor devices. [Explanation of Symbols]

[0077] 1. Imaging device 3 Radiation source 5, 6 Radiation detectors 7. Radiation energy conversion filter 10. Radiation image processing device 11 CPU 12. Radiation Image Processing Program 13 Storage 14 displays 15 Input Devices 16 memory 17 Network Interface 18 bus 21 Image acquisition unit 22 Bone Image Output Unit 23 Scattered radiation removal section 24 First derivation part 25 Second derivation part 26 Display Control Unit 30 display screen 31 Image display area 32 References Gb bone image Gd bone density image H Subject

Claims

1. Equipped with at least one processor, The aforementioned processor, Based on the first and second radiation images obtained by photographing a subject including bone with radiation having different energy distributions, an image of the bone portion of the subject is derived. Based on the first and second radiographic images, the fat percentage distribution or muscle percentage distribution of the subject is derived. Based on the first and second radiation images, the thickness distribution of the subject is derived. By referring to the relationship between the fat percentage or muscle percentage and a conversion coefficient for converting the body thickness distribution of the subject to a standard soft tissue thickness distribution, the conversion coefficient corresponding to the fat percentage distribution or muscle percentage distribution of the subject is obtained. By converting the thickness distribution of the subject using the conversion coefficient, the thickness distribution of the subject is converted to the thickness distribution of the standard soft tissue. By referring to the relationship between the body thickness distribution and the correction coefficient for correcting the pixel values ​​of the bone image to bone density, the correction coefficient corresponding to the body thickness distribution of the standard soft tissue is obtained. A radiographic image processing apparatus that derives bone density in the bone region of a subject by correcting the bone image using the correction coefficient.

2. The processor removes the scattered radiation component based on the radiation scattered by the subject from the first radiation image and the second radiation image. The radiation image processing apparatus according to claim 1, which uses the first radiation image and the second radiation image from which the scattered radiation component has been removed to derive the fat percentage distribution or the muscle percentage distribution and the bone density.

3. The radiographic image processing apparatus according to claim 1 or 2, wherein the processor displays the bone density on a display.

4. Based on the first and second radiation images obtained by photographing a subject including bone with radiation having different energy distributions, an image of the bone portion of the subject is derived. Based on the first and second radiographic images, the fat percentage distribution or muscle percentage distribution of the subject is derived. Based on the first and second radiation images, the thickness distribution of the subject is derived. By referring to the relationship between the fat percentage or muscle percentage and a conversion coefficient for converting the body thickness distribution of the subject to a standard soft tissue thickness distribution, the conversion coefficient corresponding to the fat percentage distribution or muscle percentage distribution of the subject is obtained. By converting the thickness distribution of the subject using the conversion coefficient, the thickness distribution of the subject is converted to the thickness distribution of the standard soft tissue. By referring to the relationship between the body thickness distribution and the correction coefficient for correcting the pixel values ​​of the bone image to bone density, the correction coefficient corresponding to the body thickness distribution of the standard soft tissue is obtained. A radiographic image processing method for deriving bone density in the bone region of a subject by correcting the bone image using the correction coefficient.

5. A procedure for deriving a bone image of a subject based on a first and second radiation image obtained by photographing the subject, including the bone, with radiation having different energy distributions, A procedure for deriving the fat percentage distribution or muscle percentage distribution of the subject based on the first and second radiographic images, A procedure for deriving the thickness distribution of the subject based on the first and second radiation images, A procedure for obtaining the conversion coefficient corresponding to the fat percentage distribution or muscle percentage distribution of the subject by referring to the relationship between the fat percentage or muscle percentage and the body thickness distribution of the subject to convert it to the body thickness distribution of standard soft tissue, A procedure for converting the thickness distribution of the subject to the thickness distribution of the standard soft tissue by converting the thickness distribution of the subject using the conversion coefficient, A procedure for obtaining the correction coefficient corresponding to the standard soft tissue thickness distribution by referring to the relationship between the body thickness distribution and the correction coefficient for correcting the pixel values ​​of the bone image to bone density, A radiographic image processing program that causes a computer to perform a procedure for deriving the bone density in the bone region of the subject by correcting the bone image using the correction coefficient.

Images