Source: https://patents.google.com/patent/JP6451400B2/en
Timestamp: 2020-02-23 00:05:00
Document Index: 599572856

Matched Legal Cases: ['art 17', 'art 18', 'art 53', 'art 51', 'art 13', 'art 17', 'art 181', 'art 182', 'art 183', 'art 184']

JP6451400B2 - Image processing system and image processing apparatus - Google Patents
JP6451400B2
JP6451400B2 JP2015036006A JP2015036006A JP6451400B2 JP 6451400 B2 JP6451400 B2 JP 6451400B2 JP 2015036006 A JP2015036006 A JP 2015036006A JP 2015036006 A JP2015036006 A JP 2015036006A JP 6451400 B2 JP6451400 B2 JP 6451400B2
absorption image
JP2015036006A
JP2016154766A (en
淳子 吉田
2015-02-26 Application filed by コニカミノルタ株式会社 filed Critical コニカミノルタ株式会社
2015-02-26 Priority to JP2015036006A priority Critical patent/JP6451400B2/en
2016-09-01 Publication of JP2016154766A publication Critical patent/JP2016154766A/en
2019-01-16 Publication of JP6451400B2 publication Critical patent/JP6451400B2/en
238000003384 imaging method Methods 0 claims description 94
238000003702 image correction Methods 0 claims description 50
230000001678 irradiating Effects 0 claims description 32
238000002083 X-ray spectrum Methods 0 description 13
The present invention relates to an image processing system and an image processing apparatus.
Conventionally, X-ray imaging apparatuses using the Talbot effect, such as a Talbot interferometer and a Talbot-Lau interferometer, are known. The Talbot effect is a phenomenon in which, when coherent light is transmitted through a first grating provided with slits at a certain period, the grating image is formed at a certain period in the light traveling direction. This lattice image is called a self-image, and the Talbot interferometer and the Talbot-Lau interferometer place a second grating at a position connecting the self images, and measure interference fringes generated by slightly shifting the second grating. If an object is placed in front of the second grating, moire is disturbed. Therefore, if X-ray imaging is performed with a Talbot interferometer or a Talbot-Lau interferometer, a subject is placed in front of the first grating and coherent X-rays are taken. Can be obtained, and a reconstructed image (absorption image, differential phase image, small angle scattered image) of the subject can be obtained by calculating the obtained moire fringe image.
By the way, when a moire fringe image is taken with an X-ray imaging apparatus using a Talbot interferometer or a Talbot-low interferometer as described above, and it is simply reconstructed, image unevenness due to the grating used is reflected. .
Therefore, in order to reduce this image unevenness, a subject reconstructed image generated based on the subject moire fringe image obtained by photographing the subject is photographed without placing the subject (hereinafter referred to as background photographing (BG photographing). )) To subtract or divide by the background reconstructed image (BG reconstructed image) generated based on the background moire fringe image (hereinafter referred to as BG moire fringe image) obtained by Is done).
However, since the X-ray spectrum changes when passing through the subject, even in the corrected absorption image / small angle scattered image subjected to the BG process, image unevenness due to the first grating or the second grating remains.
Therefore, for example, in Patent Document 1, during BG imaging, imaging is performed with a homogeneous member that causes an X-ray spectrum change equivalent to that when X-rays pass through the subject being interposed, and the obtained BG moire is obtained. It is described that by performing BG processing using a BG reconstructed image created based on a fringe image, image unevenness remaining in the subject reconstructed image due to the change in the X-ray spectrum described above is described.
Japanese Patent Laid-Open No. 2014-135989
However, the method of Patent Document 1 needs to perform imaging in a state in which a member that causes an X-ray spectrum change equivalent to that of a subject is interposed during BG imaging, or needs to be previously captured. Compared to this, it is troublesome and complicated.
An object of the present invention is to make it possible to easily reduce image unevenness without having to intervene a member that causes an X-ray spectrum change equivalent to that of a subject during BG imaging.
A plurality of slits arranged in a direction orthogonal to the X-ray irradiation direction, and a plurality of gratings arranged side by side in the X-ray irradiation direction;
A conversion element that accumulates electric charges and generates an electric signal according to the X-rays irradiated by the X-ray source and transmitted through the plurality of gratings is arranged in a two-dimensional manner, and the electric signal generated by the conversion element is imaged. An X-ray detector that reads as a signal and obtains a moire fringe image;
An X-ray imaging apparatus using a Talbot interferometer or a Talbot-Lau interferometer,
A subject absorption image is generated based on a moire fringe image with a subject obtained by placing a subject at a subject placement position provided on the X-ray irradiation path and irradiating the X-ray from the X-ray source. Subject absorption image generation means;
A no-subject absorption image based on a no-subject moire fringe image obtained by irradiating an X-ray from the X-ray source without placing a subject at a subject placement position provided on the X-ray irradiation path. A subject-free absorption image generation means for generating;
Image unevenness correction means for correcting image unevenness of the subject absorption image using the non-subject absorption image;
An image processing system comprising :
For each pixel of the non-subject absorption image, a high-frequency component is extracted by subtracting 1 from the value obtained by dividing the pixel value by the low-frequency component value obtained by blurring the non-subject absorption image. And a non- subject absorption image correction means for correcting only the ratio of the high-frequency component in the non-subject absorption image by multiplying the high-frequency component by a correction coefficient for matching the high-frequency component with the high-frequency component of the subject absorption image. ,
The subject-free absorbed image correction means sets a correction coefficient corresponding to a value of each pixel of an image obtained by correcting the subject-absorbed image using the subject-free absorbed image before correction to a corresponding pixel of the subject-free absorbed image. And make the correction,
The image unevenness correction unit corrects image unevenness of the subject absorption image using the non-subject absorption image corrected by the non-subject absorption image correction unit.
Invention according to claim 2,
For each pixel of the non-subject absorption image, a high-frequency component is extracted by subtracting 1 from the value obtained by dividing the pixel value by the low-frequency component value obtained by blurring the non-subject absorption image. And a non-subject absorption image correction means for correcting only the ratio of the high-frequency component in the non-subject absorption image by multiplying the high-frequency component by a correction coefficient for matching the high-frequency component with the high-frequency component of the subject absorption image. ,
The non-subject absorbed image correction means performs the correction by setting different correction coefficients for the region corresponding to the bone region and the region corresponding to the soft region in the subject absorption image,
The invention according to claim 3,
The non-subject absorption image correction means performs the correction by setting different correction coefficients for the region corresponding to the metal region and the region corresponding to the resin region in the subject absorption image,
The invention according to claim 4,
The subject-free absorbed image correction means performs the correction by setting a correction coefficient corresponding to the value of the high-frequency component in the pixel for each pixel in the subject-free absorbed image,
The invention according to claim 5 is the invention according to any one of claims 1 to 4 ,
The image unevenness correcting means corrects the image unevenness of the subject absorption image by dividing the subject absorption image by the corrected subject-free absorption image.
The invention according to claim 6,
An image processing apparatus that performs image processing on an image acquired in an X-ray imaging apparatus using a Talbot interferometer or a Talbot-Lau interferometer,
In the X-ray imaging apparatus, a subject absorption image is generated based on a moire fringe image with a subject acquired by placing a subject at a subject placement position provided on an X-ray irradiation path and irradiating the subject with X-rays. Subject absorption image generation means;
In the X-ray imaging apparatus, a subject-free absorption image is obtained based on a moiré fringe image without a subject acquired by irradiating an X-ray without placing a subject at a subject placement position provided on an X-ray irradiation path. A subject-free absorption image generation means for generating;
For each pixel of the non-subject absorption image, a high-frequency component is extracted by subtracting 1 from the value obtained by dividing the pixel value by the low-frequency component value obtained by blurring the non-subject absorption image. and a subject without absorption image correction means to correct only of the high frequency component in the object without absorption image correction factor by multiplying the high frequency component for adjusting the high-frequency component to the high frequency component of the object absorption image ,
Image unevenness correction means for correcting image unevenness of the subject absorption image using the no-subject absorption image corrected by the no-subject absorption image correction means;
The subject-free absorbed image correction means sets a correction coefficient corresponding to a value of each pixel of an image obtained by correcting the subject-absorbed image using the subject-free absorbed image before correction to a corresponding pixel of the subject-free absorbed image. Then, the correction is performed.
Invention according to claim 7,
For each pixel of the non-subject absorption image, a high-frequency component is extracted by subtracting 1 from the value obtained by dividing the pixel value by the low-frequency component value obtained by blurring the non-subject absorption image. A non-subject absorption image correction unit that corrects only the ratio of the high-frequency component in the non-subject absorption image by multiplying the high-frequency component by a correction coefficient for matching the high-frequency component to the high-frequency component of the subject absorption image;
The non-subject absorbed image correction means performs the correction by setting different correction coefficients for the region corresponding to the bone region and the soft region in the subject absorption image .
The subject-free absorbed image correction means performs the correction by setting different correction coefficients for the region corresponding to the metal region and the region corresponding to the resin region in the subject absorbed image .
The subject-free absorbed image correction means performs the correction by setting a correction coefficient corresponding to the value of the high-frequency component in each pixel in the subject-free absorbed image .
The invention according to claim 10 is the invention according to any one of claims 6 to 9, wherein
According to the present invention, it is possible to easily reduce image unevenness without the need for interposing a member that causes an X-ray spectrum change equivalent to that of a subject during BG imaging.
1 is a diagram illustrating an overall configuration of a medical image system according to an embodiment. It is a top view of a multi slit. It is a block diagram which shows the functional structure of the main-body part of FIG. It is a block diagram which shows the functional structure of the controller of FIG. It is a figure explaining the principle of a Talbot interferometer. It is a flowchart which shows the imaging | photography control processing performed by the control part of FIG. It is a figure which shows the moire fringe image obtained by imaging | photography of 5 steps. It is a subject absorption image obtained by photographing a phantom as a subject using a Talbot-Lau interferometer. FIG. 9 is a BG absorption image obtained by shooting without placing a phantom under the same shooting conditions as in FIG. 8. FIG. (A) is an image obtained by taking the natural logarithm of the absorption image I AB obtained by dividing the subject absorption image of FIG. 8 by the BG absorption image of FIG. 9, and (b) is a differential absorption image of the image of (a). It is. It is a flowchart of the absorption image generation process performed by the control part of FIG. (A) is a figure which shows the correction coefficient for streak unevenness correction, (b) is a figure which shows the correction coefficient for granular unevenness correction. (A) is the image which took the natural logarithm of the corrected absorption image produced | generated using the correction coefficient for soft tissues to the whole image, (b) is a differential absorption image of the image of (a). (A) is an image obtained by taking a natural logarithm of a corrected absorption image generated using a correction coefficient for a bone portion in a bone portion region and a soft tissue region in a soft tissue region, and (b) is (a) It is a differential absorption image of the image. It is an image obtained by taking the natural logarithm of a corrected absorption image generated using a correction coefficient for correcting streak unevenness over the entire image. This is an image obtained by taking the natural logarithm of a corrected absorption image generated using a correction coefficient for correcting streak unevenness for streaky unevenness and for a granular unevenness.
FIG. 1 shows a medical image system according to an embodiment of the present invention. The medical image system includes an X-ray imaging apparatus 1 and a controller 5. The X-ray imaging apparatus 1 performs X-ray imaging using a Talbot-Lau interferometer, and the controller 5 generates a reconstructed image of the subject using a plurality of moire fringe images obtained by the X-ray imaging.
As shown in FIG. 1, the X-ray imaging apparatus 1 includes an X-ray source 11, a multi-slit 12, a subject table 13, a first grating 14, a second grating 15, an X-ray detector 16, a holding part 17, and a body part 18. Etc.
The X-ray imaging apparatus 1 is a vertical type, and an X-ray source 11, a multi slit 12, a subject table 13, a first grating 14, a second grating 15, and an X-ray detector 16 are arranged in this order in the z direction, which is the gravitational direction. Placed in. The distance between the focal point of the X-ray source 11 and the multi-slit 12 is d 1 (mm), the distance between the focal point of the X-ray source 11 and the X-ray detector 16 is d 2 (mm), and the distance between the multi-slit 12 and the first grating 14. The distance is represented by d3 (mm), and the distance between the first grating 14 and the second grating 15 is represented by d4 (mm). The position of the subject table 13 may be provided between the first grid 14 and the second grid 15.
The distance d1 is preferably 5 to 500 (mm), more preferably 5 to 300 (mm).
The distance d2 is generally at least 3000 (mm) or less because the height of the photographing room is generally about 3 (m) or less. Especially, the distance d2 is preferably 400 to 3000 (mm), more preferably 500 to 2000 (mm).
The distance (d1 + d3) between the focal point of the X-ray source 11 and the first grating 14 is preferably 300 to 3000 (mm), and more preferably 400 to 1800 (mm).
The distance (d1 + d3 + d4) between the focal point of the X-ray source 11 and the second grating 15 is preferably 400 to 3000 (mm), and more preferably 500 to 2000 (mm).
Each distance may be set by calculating an optimum distance at which the lattice image (self-image) by the first lattice 14 overlaps the second lattice 15 from the wavelength of the X-rays emitted from the X-ray source 11.
The X-ray source 11, the multi slit 12, the subject table 13, the first grating 14, the second grating 15, and the X-ray detector 16 are integrally held by the same holding unit 17 and the positional relationship in the z direction is fixed. ing. The holding portion 17 is formed in an arm shape, and is attached to the main body portion 18 so as to be movable in the z direction by a driving portion 18 a provided in the main body portion 18.
The X-ray source 11 is held via a buffer member 17a. Any material may be used for the buffer member 17a as long as it can absorb shocks and vibrations, and examples thereof include an elastomer. Since the X-ray source 11 generates heat upon irradiation with X-rays, it is preferable that the buffer member 17a on the X-ray source 11 side is additionally a heat insulating material.
The X-ray source 11 includes an X-ray tube, generates X-rays by the X-ray tube, and irradiates the X-rays in the gravity direction (z direction). As the X-ray tube, for example, a Coolidge X-ray tube or a rotary anode X-ray tube widely used in the medical field can be used. As the anode, tungsten or molybdenum can be used.
The focal diameter of the X-ray is preferably 0.03 to 3 (mm), more preferably 0.1 to 1 (mm).
In the X-ray irradiation direction of the X-ray source 11, an irradiation field stop (not shown) for narrowing the X-ray irradiation range is provided.
The multi slit 12 (third grating) is a diffraction grating, and a plurality of slits are provided at predetermined intervals in the x direction as shown in FIG. The multi-slit 12 is formed on a substrate having a low X-ray absorption rate such as silicon or glass by using a material having a high X-ray shielding power such as tungsten, lead, or gold, that is, a high X-ray absorption rate. For example, the resist layer is masked in a slit shape by photolithography, and UV is irradiated to transfer the slit pattern to the resist layer. A slit structure having the same shape as the pattern is obtained by exposure, and a metal is embedded between the slit structures by electroforming to form a multi-slit 12.
The slit period of the multi slit 12 is 1 to 60 (μm). As shown in FIG. 2, the slit period is defined as a period between adjacent slits. The width (length in the x direction) of the slit is 1 to 60 (%) of the slit period, and more preferably 10 to 40 (%). The height (length in the z direction) of the slit is 1 to 500 (μm), preferably 1 to 150 (μm).
When the slit period of the multi slit 12 is w 0 (μm) and the slit period of the first grating 14 is w 1 (μm), the slit period w 0 can be obtained by the following equation.
w 0 = w 1 · (d3 + d4) / d4
By determining the period w 0 so as to satisfy the equation, the self-images formed by the X-rays that have passed through the slits of the multi-slit 12 and the first grating 14 overlap each other on the second grating 15. Can be in a suitable state.
As shown in FIG. 1, a drive unit 12 a that moves the multi slit 12 in the x direction orthogonal to the z direction is provided adjacent to the multi slit 12. As the drive unit 12a, for example, a relatively large reduction ratio type drive mechanism such as a worm reducer can be used alone or in combination.
The subject table 13 is a table for placing a subject provided at a subject placement position on the X-ray irradiation path from the X-ray source 11.
The first grating 14 is a diffraction grating in which a plurality of slits having a predetermined period are provided in the x direction as in the multi-slit 12 (see FIG. 2). The first lattice 14 can be formed by photolithography using UV as in the case of the multi-slit 12, or a silicon substrate is deeply digged with a fine fine line by a so-called ICP method to form a lattice structure only with silicon. It is good as well. The slit period of the first grating 14 is 1 to 20 (μm). The width of the slit is 20 to 70 (%) of the slit period, and preferably 35 to 60 (%). The height of the slit is 1 to 100 (μm).
When a phase type is used as the first grating 14, the height of the slit is π / 8 to 15 × π / phase difference between the two materials forming the slit period, that is, the material of the X-ray transmitting part and the X-ray shielding part. The height is 8. The height is preferably π / 2 or π. When an absorption type is used as the first grating 14, the height of the slit is set to a height at which X-rays are sufficiently absorbed by the X-ray shielding part.
When the first grating 14 is a phase type having a phase difference of π / 2 due to the materials of the X-ray transmitting part and the X-ray shielding part, the distance d4 between the first grating 14 and the second grating 15 is approximately the following condition: It is necessary to satisfy.
d4 = (m + 1/2) · w 1 2 / λ
Note that m is an integer, and λ is the wavelength of X-rays.
The second grating 15 is a diffraction grating in which a plurality of slits having a predetermined period are provided in the x direction as in the multi-slit 12 (see FIG. 2). The second grating 15 can also be formed by photolithography. The slit period of the second grating 15 is 1 to 20 (μm). The width of the slit is 30 to 70 (%) of the slit cycle, and preferably 35 to 60 (%). The height of the slit is 1 to 100 (μm).
In the present embodiment, each of the first grating 14 and the second grating 15 has a lattice plane perpendicular to the z direction (parallel in the xy plane), and the slit direction of the first grating 14 and the second grating 15. The slit direction is arranged with a predetermined angle (slightly) in the xy plane, but both may be arranged in parallel.
The multi-slit 12, the first grating 14, and the second grating 15 can be configured as follows, for example.
Focal diameter of X-ray source 11: 300 (μm), tube voltage: 40 (kVp), additional filter: aluminum 1.6 (mm)
Distance d1 from the focal point of the X-ray source 11 to the multi slit 12: 240 (mm)
Distance d3 from the multi slit 12 to the first grating 14: 1110 (mm)
Distance d3 + d4: 1370 (mm) from the multi slit 12 to the second grating 15
Multi slit 12 size: 10 (mm square), slit period: 22.8 (μm)
Size of the first grating 14: 50 (mm square), slit period: 4.3 (μm)
Size of the second grating 15: 50 (mm square), slit period: 5.3 (μm)
The X-ray detector 16 has two-dimensionally arranged conversion elements that generate electric signals in accordance with the irradiated X-rays, and reads the electric signals generated by the conversion elements as image signals.
The pixel size of the X-ray detector 16 is 10 to 300 (μm), more preferably 50 to 200 (μm).
It is preferable that the position of the X-ray detector 16 is fixed to the holding unit 17 so as to contact the second grating 15. This is because the moire fringe image obtained by the X-ray detector 16 becomes blurred as the distance between the second grating 15 and the X-ray detector 16 increases.
As the X-ray detector 16, an FPD (Flat Panel Detector) can be used. F
There are two types of PDs: an indirect conversion type that converts detected X-rays into electrical signals via a photoelectric conversion element, and a direct conversion type that converts detected X-rays directly into electrical signals. Also good.
In the indirect conversion type, photoelectric conversion elements are two-dimensionally arranged together with TFTs (thin film transistors) under a scintillator plate such as CsI or Gd 2 O 2 S to constitute each pixel. When the X-rays incident on the X-ray detector 16 are absorbed by the scintillator plate, the scintillator plate emits light. Charges are accumulated in each photoelectric conversion element by the emitted light, and the accumulated charges are read as an image signal.
In the direct conversion type, an amorphous selenium film having a film pressure of 100 to 1000 (μm) is formed on glass by thermal vapor deposition of amorphous selenium, and the amorphous selenium film and the electrode are arranged on a two-dimensionally arranged TFT array. Vapor deposited. When the amorphous selenium film absorbs X-rays, a voltage is released in the material in the form of electron-hole pairs, and a voltage signal between the electrodes is read by the TFT.
Note that imaging means such as a CCD (Charge Coupled Device), an X-ray camera or the like is used as the X-ray detector 16.
It may be used as
As shown in FIG. 3, the main body unit 18 includes a control unit 181, an operation unit 182, a display unit 183, a communication unit 184, a storage unit 185, and the like.
The control unit 181 includes a CPU (Central Processing Unit), a RAM (Random Access Memory), and the like, and executes various processes in cooperation with a program stored in the storage unit 185. The control unit 181 is connected to each unit such as the X-ray source 11, the drive unit 12 a, the drive unit 18 a, and the X-ray detector 16. The timing of X-ray irradiation from the X-ray source 11, the X-ray irradiation conditions, the reading timing of the image signal by the X-ray detector 16, the movement of the multi-slit 12 and the like are controlled in accordance with the imaging condition setting information.
The operation unit 182 includes an exposure switch and the like, generates an operation signal corresponding to these operations, and outputs the operation signal to the control unit 181.
The display unit 183 displays the operation screen, the operation status of the X-ray imaging apparatus 1 and the like on the display according to the display control of the control unit 181.
The communication unit 184 includes a communication interface and communicates with the controller 5 on the network. For example, the communication unit 184 transmits the moire fringe image read by the X-ray detector 16 and stored in the storage unit 185 to the controller 5.
The storage unit 185 stores a program executed by the control unit 181 and data necessary for executing the program. The storage unit 185 stores the moire fringe image obtained by the X-ray detector 16.
The controller 5 controls the imaging operation of the X-ray imaging apparatus 1 according to the operation by the operator. The controller 5 performs image processing on a series of moire fringe images obtained by the X-ray imaging apparatus 1 as an image processing apparatus. For example, a reconstructed image (absorption image, small angle scattered image, differential phase image) of a subject is generated using a series of moire fringe images obtained by the X-ray imaging apparatus 1.
As illustrated in FIG. 4, the controller 5 includes a control unit 51, an operation unit 52, a display unit 53, a communication unit 54, and a storage unit 55.
The control unit 51 includes a CPU (Central Processing Unit) and a RAM (Random Access Memory).
) And the like, and various processes including an absorption image generation process to be described later are executed in cooperation with a program stored in the storage unit 55. The control unit 51 functions as a subject absorption image generation unit, a non-subject absorption image generation unit, a non-subject absorption image correction unit, and an image unevenness correction unit.
The operation unit 52 includes a keyboard having cursor keys, numeric input keys, various function keys, and the like, and a pointing device such as a mouse, and includes a key pressing signal pressed by the keyboard and an operation signal by the mouse. Is output to the control unit 51 as an input signal. It is good also as a structure provided with the touchscreen comprised integrally with the display of the display part 53, and producing | generating the operation signal according to these operation to the control part 51. FIG.
The display unit 53 includes, for example, a monitor such as a CRT (Cathode Ray Tube) or an LCD (Liquid Crystal Display), and displays an operation screen, a generated reconstructed image, and the like according to display control of the control unit 51. indicate.
The communication unit 54 includes a communication interface and communicates with the X-ray imaging apparatus 1 and the X-ray detector 16 on the network by wire or wirelessly. For example, the communication unit 54 transmits imaging conditions and control signals to the X-ray imaging apparatus 1 and receives a moire fringe image from the X-ray imaging apparatus 1 or the X-ray detector 16.
The storage unit 55 stores a program executed by the control unit 51 and data necessary for executing the program. For example, the storage unit 55 stores imaging order information that is imaging information reserved by a RIS (Radiology Information System), a HIS (Hospital Information System), or the like (not shown). The imaging order information includes patient information such as a patient ID and a patient name, imaging part (subject part) information, and the like.
In addition, the storage unit 55 stores an imaging condition table in which a subject part is associated with an imaging condition suitable for photographing the subject part.
Further, the storage unit 55 stores the moire fringe image acquired by the X-ray imaging apparatus 1 based on the imaging order information, the reconstructed image generated based on the moire fringe image, and the like in association with the imaging order information.
The storage unit 55 stores in advance gain correction data corresponding to the X-ray detector 16, a defective pixel map, and the like. The defective pixel map is position information (coordinates) of defective pixels (including those without pixels) of the X-ray detector 16.
<Operation of medical imaging system>
Here, an X-ray imaging method using the Talbot-Lau interferometer of the X-ray imaging apparatus 1 will be described.
As shown in FIG. 5, when X-rays emitted from the X-ray source 11 pass through the first grating 14, the transmitted X-rays form an image at a constant interval in the z direction. This image is called a self-image, and the phenomenon in which a self-image is formed is called the Talbot effect. The second grating 15 is arranged substantially parallel to the self-image at a position connecting the self-images, and a moire fringe image (indicated by Mo in FIG. 5) is obtained by X-rays transmitted through the second grating 15. That is, the first grating 14 forms a periodic pattern, and the second grating 15 converts the periodic pattern into moire fringes. If a subject (indicated by H in FIG. 5) exists between the X-ray source 11 and the first grating 14, the phase of the X-ray is shifted depending on the subject. Therefore, as shown in FIG. Disturbed by borders. The disturbance of the moire fringes can be detected by processing the moire fringe image, and the subject image can be imaged. This is the principle of the Talbot interferometer.
In the X-ray imaging apparatus 1, a multi-slit 12 is disposed near the X-ray source 11 between the X-ray source 11 and the first grating 14, and X-ray imaging using a Talbot-Lau interferometer is performed. The Talbot interferometer is based on the premise that the X-ray source 11 is an ideal point source. However, since a focal point having a large focal diameter is used for actual imaging, it is as if a plurality of point sources are connected by the multi slit 12. Multiple light sources are used as if they were irradiated with X-rays. This is an X-ray imaging method using a Talbot-Lau interferometer, and a Talbot effect similar to that of a Talbot interferometer can be obtained even when the focal diameter is somewhat large.
In the present embodiment, photographing is performed by a fringe scanning method. In the fringe scanning, in general, any one of the gratings (the multi slit 12, the first grating 14, and the second grating 15) (in this embodiment, the multi slit 12) is formed in the slit period direction ( The image is moved M times (M is a positive integer, M> 2) while moving relatively in the x direction), and M moire necessary to generate one reconstructed image is obtained. It means obtaining a fringe image. Specifically, assuming that the slit period of the grating to be moved is d (μm), the imaging is repeated by moving the grating in the slit period direction by d / M (μm) to obtain M moiré fringe images. .
In the conventional Talbot-Lau interferometer, the multi-slit 12 is used for the purpose of increasing the number of light sources and increasing the irradiation dose as described above, and in order to obtain a plurality of moire fringe images by the fringe scanning method, The second grating 15 was moved relative to each other. However, in the present embodiment, the first grating 14 or the second grating 15 is not moved relatively, but the positions of the first grating 14 and the second grating 15 are fixed and the first grating 14 and the second grating 15 are fixed. On the other hand, by moving the multi slit 12, a plurality of moire fringe images having a constant periodic interval are obtained.
FIG. 6 is a flowchart showing an imaging control process executed by the control unit 181 of the X-ray imaging apparatus 1. The flow of the imaging control process will be described with reference to FIG.
First, when the exposure switch of the operation unit 182 is operated by the operator (step S1; YES), the control unit 181 controls the X-ray source 11, the X-ray detector 16, and the drive unit 12a to perform a series of steps. And a series of moire fringe images having different moire fringe phases are acquired (step S2).
In a series of imaging, first, X-ray irradiation by the X-ray source 11 is started with the multi-slit 12 stopped. In the X-ray detector 16, after resetting to remove unnecessary charges remaining in the previous imaging, charges are accumulated at the timing of X-ray irradiation, and the charges accumulated at the timing of stopping X-ray irradiation are converted into image signals. Is read as This is one step of shooting. The movement of the multi-slit 12 is started at the timing when the photographing for one step is finished, stopped when the predetermined amount is moved, and the photographing of the next step is performed. In this manner, the movement and stop of the multi-slit 12 are repeated for a predetermined number of steps, and when the multi-slit 12 is stopped, X-ray irradiation and image signal reading are performed. When shooting with the multi-slit 12 moving by one slit period is completed, a series of shooting for acquiring a plurality of moire fringe images necessary to generate one reconstructed image is completed.
The number M of steps in a series of photographing is 2 to 20, more preferably 3 to 10. From the viewpoint of obtaining a reconstructed image with high visibility in a short time, 5 steps are preferable (Reference 1: K. Hibino, BFOreb and DIFarrant, Phase shifting for nonsinusoidal wave forms with phase-shift errors, J. Opt. Soc. Am. A, Vol. 12, 761-768 (1995), Reference 2: A. Momose, W. Yashiro, Y. Takeda, Y. Suzuki and T. Hattori, Phase Tomography by X-ray Talbot Interferometetry for biological imaging, Jpn. J. Appl. Phys., Vol. 45, 5254-5262 (2006)). Here, description will be made on the assumption that shooting is performed in five steps.
For example, assume that the slit period of the multi-slit 12 is 22.8 (μm), and five-step imaging is performed in 10 seconds. Shooting is performed every time the multi slit 12 moves and stops 4.56 (μm) corresponding to 1/5 of the slit period. In terms of shooting time, shooting is performed 2, 4, 6, 8, and 10 seconds after the exposure switch is turned on. When the multi-slit 12 can be moved at a constant feed amount with ideal feed accuracy, five moire fringe images corresponding to one slit period of the multi-slit 12 can be obtained by photographing in five steps as shown in FIG. It is done.
When the photographing of each series of steps is completed, the control unit 181 causes the controller 5 to transmit the moire fringe image of each step by the communication unit 184 (step S3). The communication unit 184 may send the image to the controller 5 one by one every time the shooting of each step is completed. After the shooting of each step is completed and all the moire fringe images are obtained, the images are summarized. May be transmitted.
In the present embodiment, shooting with the subject placed on the subject table 13 (subject shooting) and BG shooting without placing the subject on the subject table 13 are performed, and the subject moire fringe image (moire fringe image with subject) and BG are taken. A moire fringe image (a moire fringe image without a subject) is acquired.
In the controller 5, when the communication unit 54 receives a series of subject moire fringe images and BG moire fringe images from the main body 18, the control unit 51 converts the received subject moire fringe images and BG moire fringe images into a series. Based on this, reconstructed images such as an absorption image, a differential phase image, and a small angle scattered image are generated.
Of the reconstructed image, the absorption image I AB (x, y) is generally following [Expression 1] are generated by (reference 3: A.Momose, W.Yashiro, H.Kuwabara and K. Kawabata, Grating-Based X-ray Phase Imaging Using Multiline X-ray Source, Jpn. J. Appl. Phys., Vol. 48, 076512 (2009)).
I k SAMPLE (x, y) and I k BG (x, y) represent pixel values of the subject moire fringe image and BG moire fringe image of the k-th shooting, respectively. (x, y) represents the two-dimensional coordinates of the pixel in each image. I AB_SAMPLE (x, y) and I AB_BG (x, y) represent the pixel values of the subject absorption image and the BG absorption image (absence without subject absorption image), respectively.
FIG. 8 shows an example of a subject absorption image I AB_SAMPLE (x, y) obtained by photographing a phantom as a subject using a Talbot-Lau interferometer with a tungsten tube and 40 kVp (additional filter, aluminum 2.3 mm). . FIG. 9 shows an example of a BG absorption image I AB_BG (x, y) obtained by photographing without placing a phantom under the same photographing conditions as in FIG. Note that the arrows and broken lines shown in FIG. 9 and subsequent figures indicate image unevenness. In each image of FIG. 8 and FIG. 9, fine streaks are included in the vertical direction on the paper surface. This is mainly due to image unevenness caused by unevenness in the height and period of the lattice generated when the second lattice 15 is manufactured. Conceivable.
Image unevenness due to the grid should appear in the same manner in the subject absorption image I AB_SAMPLE (x, y) and the BG absorption image I AB_BG (x, y), and the subject absorption image I AB_SAMPLE (x, y) and the BG absorption image By taking the ratio of I AB — BG (x, y), it should be possible to obtain an absorption image I AB (x, y) representing the attenuation of X-rays by the subject without the influence of unevenness due to the lattice. However, in practice, the absorption image I AB (x, y) generated by taking the ratio of the subject absorption image I AB_SAMPLE (x, y) and the BG absorption image I AB_BG (x, y) It is known that unevenness remains.
10A shows an absorption image I AB obtained by dividing the subject absorption image I AB_SAMPLE (x, y) of FIG. 8 by the BG absorption image I AB_BG (x, y) of FIG. 9 according to [Equation 1]. This is an image obtained by taking the natural logarithm of (x, y). FIG. 10B is an image for visualizing the image unevenness remaining in the image of FIG. 10A in an easy-to-understand manner. The x direction of each pixel in the image of FIG. It is a differential absorption image calculated by subtracting pixel values of both adjacent pixels.
In the images shown in FIGS. 10A and 10B, vertical streaky irregularities (indicated by arrows in FIGS. 8 and 10A) and granular irregularities (see FIG. 8 and FIG. 9). 8, the region surrounded by the broken line in FIG. 10A remains. These are image unevenness commonly seen in the subject absorption image I AB_SAMPLE (x, y) and the BG absorption image I AB_BG (x, y), but the ratio of the unevenness to the pixel value is the subject absorption image I AB_SAMPLE (x, y). ) And the BG absorption image I AB_BG (x, y), it is considered that image unevenness also remains in the absorption image I AB (x, y) after BG processing. The ratio of unevenness to the pixel value differs between the subject absorption image I AB_SAMPLE (x, y) and the BG absorption image I AB_BG (x, y) because the X-rays on the low energy side are transmitted when the X-rays pass through the subject. This is because the X-ray spectrum distribution changes due to absorption, and the higher energy side increases. Generally, unevenness is less likely to occur when the energy of X-rays is higher. In other words, the subject absorption image I AB_SAMPLE (x, y) has a smaller ratio of unevenness to the pixel value than the BG absorption image I AB_BG (x, y).
Therefore, the inventor of the present application extracts a high frequency component corresponding to image unevenness from the BG absorption image I AB — BG (x, y), and adds a correction coefficient a (x, y) (hereinafter referred to as coefficient a (x, y, As shown in FIGS. 10 (a) and 10 (b), correction is performed by multiplying a ≦ 1) and attenuating the high-frequency component ratio in accordance with the subject absorption image I AB_SAMPLE (x, y). The present inventors have found that it is possible to reduce (remove) image unevenness.
In Patent Document 1, image unevenness is reduced (attenuated) on the premise that imaging is performed in a state where a member that causes an X-ray spectrum change equivalent to that of a subject is interposed during BG imaging. In this case, since the amount of X-rays reaching the X-ray detector is reduced by interposing the member, the pixel value of the entire BG absorption image is reduced, and the BG absorption image I AB_BG ( There is also the disadvantage that noise increases compared to x, y). In the present invention, since the image unevenness is removed by correcting only the high-frequency component of the BG absorption image I AB — BG (x, y) obtained by the conventional BG imaging, the image unevenness can be removed without degrading the image quality.
In FIGS. 10A and 10B, an image obtained by taking the natural logarithm of the absorption image I AB (x, y) is shown as an example. However, in the X-ray image, the transmittance of X-rays by the subject is shown. An image to be represented (here, corresponding to an absorption image) is generally displayed according to a value proportional to the thickness of the subject by taking a natural logarithm. In the following description, absorption image I AB (x, y) corrected absorption image I AB2 with reduced image unevenness remaining in (x, y) is described process for obtaining a corrected absorption image I AB2 (x, If the image unevenness of y) is reduced, even the image having the natural logarithm reduces the image unevenness. Therefore, as an image showing the effect of the processing of the present embodiment, the natural logarithm of the absorption image which is a general display method is used. An example of an image taken is shown (see FIGS. 13A, 14A, 15, and 16).
Hereinafter, the absorption image generation process including the correction process of the ratio of the high frequency component of the BG absorption image I AB_BG (x, y) will be described. FIG. 11 shows a flowchart of the absorption image generation process executed by the controller 51 of the controller 5. The absorption image generation process is executed in cooperation with the control unit 51 and a program stored in the storage unit 55.
First, the control unit 51 generates a subject absorption image I AB_SAMPLE (x, y) based on a series of subject moiré fringe images received by the communication unit 54, and generates a series of BG moiré fringe images received by the communication unit 54. Based on this, a BG absorption image I AB_BG (x, y) is generated (step S11).
In step S11, the subject absorption image I AB_SAMPLE (x, y) is generated using [Formula 2], and the BG absorption image I AB_BG (x, y) is generated using [Formula 3]. The moire fringe image is preferably subjected to an offset correction process, a gain correction process, a defective pixel correction process, and the like in advance.
Next, the control unit 51 extracts a high frequency component of the BG absorption image I AB_BG (x, y) (step S12).
BG absorption image I AB_BG (x, y) and the high frequency component of, represents the ratio of the mean value of the peripheral pixels of the fine irregularity structure of the BG absorption image I AB_BG (x, y). In step S2, first, for example, the BG absorption image I AB_BG (x, y) is averaged with surrounding pixels, or a filtering process for giving a blur such as a Gaussian filter is performed to obtain the BG absorption image I AB_BG (x, y). ) Low frequency component L AB_BG (x, y) is obtained. Next, the high frequency component H AB_BG (x, y) is extracted based on the BG absorption image I AB_BG (x, y) and the low frequency component L AB_BG (x, y). Here, by dividing the BG absorption image I AB — BG (x, y) by the low frequency component L AB — BG (x, y), the pixel value is 1 + high frequency component (high frequency component represented with the pixel value 1 as a reference). Is obtained. That is, the high frequency component H AB_BG (x, y) of the BG absorption image I AB_BG (x, y) can be extracted by the following [Equation 4].
High frequency component H AB — BG (x, y) = (BG absorption image I AB — BG (x, y) / low frequency component L AB — BG (x, y)) − 1 −1 [Expression 4]
Next, the control unit 51 corrects the ratio of the high frequency component H AB_BG (x, y) of the BG absorption image I AB_BG (x, y) (step S13).
To obtain a corrected BG absorption image I AB_BG2 (x, y) in which the ratio of the high frequency component of the BG absorption image I AB_BG (x, y) is matched with the high frequency component of the subject absorption image I AB_SAMPLE (x, y). For the high frequency component H AB_BG (x, y) on the left side of [Equation 4], a coefficient a (x, y) for adjusting the ratio of this high frequency component to the high frequency component of the subject absorption image I AB_SAMPLE (x, y) is A BG absorption image may be obtained from the multiplied expression. That is, the corrected BG absorption image I AB — BG2 (x, y) in which only the ratio of the high frequency component is corrected can be obtained by [Equation 5] below.
Corrected BG absorption image I AB_BG2 (x, y) = low frequency component L AB_BG (x, y) × (1 + a (x, y) × high frequency component H AB_BG (x, y)) [Equation 5]
Next, the control unit 51 performs BG processing using the corrected BG absorption image I AB_BG2 (x, y) (step S14). That is, by dividing the subject absorption image I AB_SAMPLE (x, y) to a corrected BG absorption image I AB_BG2 (x, y), to correct the image unevenness of the object absorption image I AB_SAMPLE (x, y). As a result, a corrected absorption image I AB2 (x, y) with reduced image unevenness is generated.
Here, the coefficient a (x, y) used in step S13 is preferably set to a value corresponding to the pixel value of the absorption image I AB (x, y) after BG processing. If the subject material is the same, the thicker the subject, the more X-ray absorption by the subject increases, so the change in the X-ray spectrum increases. As described above, the ratio of high-frequency components caused by image unevenness on the subject absorption image I AB_SAMPLE (x, y) varies depending on the X-ray spectrum. Therefore, a value corresponding to the pixel value of the absorption image I AB (x, y) after BG processing indicating the value corresponding to the thickness of the subject when the subject is the same material is set as the coefficient a (x, y). It is preferable. Since influence on the pixel values of the image unevenness in the absorption image I AB after BG treatment (x, y) is small, the subject absorption image I AB_SAMPLE (x, y) the uncorrected BG absorption image I AB_BG ( There is no problem even if a (x, y) is determined based on the pixel value of the absorption image I AB (x, y) after the conventional BG processing divided by x, y).
The coefficient a (x, y) is preferably set for each subject material that changes in the same X-ray spectrum. If the subject material is different, the X-ray spectrum after passing through the subject will be different even if the pixel value of the absorption image I AB (x, y) is the same, resulting in image unevenness on the subject absorption image I AB_SAMPLE (x, y). This is because the ratio of the resulting high frequency component is different.
In addition, as described above, the image unevenness due to the lattice of the subject absorption image I AB_SAMPLE (x, y) and the BG absorption image I AB_BG (x, y) includes streak (vertical streak) unevenness and granular unevenness. Although existing, the same coefficient a (x, y) cannot be corrected with high accuracy. For this reason, the coefficient a (x, y) is preferably set separately for streak unevenness correction and granular unevenness correction.
Therefore, in the present embodiment, as shown in FIGS. 12A and 12B, the coefficient a is used for correcting streak unevenness (see FIG. 12A) and granular unevenness correction (see FIG. 12B). ) For each subject material (here, bone part, soft part (soft tissue)), the pixel value of the absorption image I AB (x, y) after BG processing (or a value obtained by taking the natural logarithm thereof) ) Is converted into a function, or a table is created and stored in the storage unit 55 in advance. In step S13 of FIG. 11, the control unit 51 calculates the pixel value (or its natural logarithm) of the absorption image I AB (x, y) after the BG processing based on the function or table stored in the storage unit 55. The coefficient a (x, y) corresponding to the (taken value) is obtained and set. The coefficient a in 12 (a) and 12 (b) is obtained experimentally and empirically. Specifically, first, for a material such as acrylic or aluminum, the relationship between the coefficient a at which the image unevenness of the absorption image I AB (x, y) is most attenuated and the absorption image pixel value I AB (x, y) Ask for. On the basis of these data, a coefficient a is created so that the unevenness remaining most in the absorption image I AB (x, y) obtained by photographing the phantom of the bone part or the soft part or the actual photographing target is eliminated. For example, as for the coefficients shown in FIGS. 12A and 12B, the coefficient a obtained with acrylic is applied as the coefficient a for the soft part, and the coefficient a obtained with acrylic as the average value of the coefficient a obtained with acrylic / aluminum. Yes. Note that such a method for obtaining the coefficient a is not limited to the bone part and the soft part, but is particularly useful in the human body.
FIG. 13A shows a corrected absorption image I AB2 (x, y) based on the subject absorption image I AB_SAMPLE (x, y) shown in FIG. 8 and the BG absorption image I AB_BG (x, y) shown in FIG. Is a natural logarithm of the corrected absorption image I AB2 (x, y) generated using the soft tissue coefficient a (x, y) for the entire image. FIG. 13B is a differential absorption image calculated by subtracting pixel values of pixels adjacent to each other in the x direction (left and right direction in the drawing) of each pixel of the image in FIG.
FIG. 14A shows a corrected absorption image I AB2 (x, y) based on the subject absorption image I AB_SAMPLE (x, y) shown in FIG. 8 and the BG absorption image I AB_BG (x, y) shown in FIG. , The corrected absorption image I AB_BG (x, y) generated using the coefficient a (x, y) for the bone part B1 and the soft tissue B2 for the soft part area B2 is generated. It is an image that takes the natural logarithm of. FIG. 14B is a differential absorption image calculated by subtracting pixel values of pixels adjacent to each other in the x direction (left and right direction in the drawing) of each pixel of the image in FIG.
13 (a) and 14 (a) use the coefficient a for correcting streak unevenness shown in FIG. 12 (a) for streaky unevenness, and FIG. 12 (b) for granular unevenness. It is generated using the coefficient a for correcting graininess irregularity shown in FIG.
As shown in FIGS. 13A and 13B, when the entire image is corrected with the coefficient a (x, y) for soft tissue, FIGS. 10A and 10B generated by the conventional method are shown in FIGS. It can be seen that although the stripe-like unevenness and the granular unevenness are both attenuated, the stripe-like unevenness remains in the bone (for example, see the arrow in FIG. 13B). On the other hand, when the bone region B1 and the soft region B2 are corrected using the bone and soft tissue coefficients a (x, y), respectively, in FIGS. 14 (a) and 14 (b). As shown, it can be seen that the unevenness of the longitudinal muscles can be further reduced for the bone.
Regarding the discrimination between the bone region B1 and the soft region B2, a threshold value is provided for the pixel value of the small-angle scattered image Isc (x, y) generated from the same moire fringe image, and the bone portion is equal to or greater than the threshold value. It can be identified as a soft part. The small-angle scattered image I sc (x, y) can be obtained by the following [Equation 6].
I sc_SAMPLE (x, y) and I sc_BG (x, y) represent pixel values of the subject small-angle scattered image and the BG small-angle scattered image, respectively.
Further, FIG. 15, generates an object absorption image I AB_SAMPLE shown in FIG. 8 (x, y), BG absorption image I AB_BG shown in FIG. 9 (x, y) corrected absorption image I AB2 from (x, y) In this case, it is an image obtained by taking the natural logarithm of the corrected absorption image I AB2 (x, y) generated using the coefficient a (x, y) for correcting streak unevenness over the entire image. FIG. 16 shows a case where the corrected absorption image I AB2 (x, y) is generated from the subject absorption image I AB_SAMPLE (x, y) shown in FIG. 8 and the BG absorption image I AB_BG (x, y) shown in FIG. This is an image obtained by taking the natural logarithm of the corrected absorption image generated by using the coefficient a (x, y) for correcting the streak unevenness for the streaky unevenness and for the granular unevenness for the granular unevenness correction. .
15 and 16 are generated using the coefficient a for the bone part for the bone part and the coefficient a for the soft tissue for the soft part.
In the image obtained by taking the natural logarithm of the BG-processed absorption image I AB (x, y) shown in FIG. 10 (a), the granular structure seen in the area surrounded by the broken line is drawn white compared to the surrounding area. However, in FIG. 15, it is drawn in black. This is because the coefficient a (x, y) appropriate for the stripe-like unevenness is smaller than the coefficient a (x, y) appropriate for the granular unevenness (see FIGS. 12A and 12B). ), Indicating that the correction has been applied excessively. Therefore, by setting different coefficients a (x, y) for the stripe-like unevenness and the granular unevenness, the granular image unevenness can be appropriately corrected as shown in FIG.
The discrimination between the stripe-like unevenness and the granular unevenness can be performed, for example, by comparing the value of the high-frequency component H AB_BG (x, y) calculated by [Equation 4] with a predetermined threshold value. As can be seen from the stripe-like unevenness and the density of the granular unevenness depicted in the BG absorption image I AB — BG (x, y) in FIG. 9, the change in pixel value is larger in the granular unevenness than in the surrounding area. That is, the value of the high frequency component H AB_BG (x, y) is small (the absolute value is large, but the value is small with respect to the pixel values of the surrounding pixels). Therefore, the value of the high frequency component H AB_BG (x, y) calculated by [Equation 4] is compared with a predetermined threshold value. Can be determined.
As described above, the control unit 51 of the controller 5 generates the subject absorption image I AB_SAMPLE (x, y) based on the subject moire fringe image acquired by the X-ray imaging apparatus 1 and based on the BG moire fringe image. Te BG absorption image I AB_BG (x, y) generates, BG absorption image I AB_BG (x, y) to correct the proportion of the high-frequency components in the resulting corrected BG absorption image I AB_BG2 the (x, y) The corrected absorption image I AB2 (x, y) is generated by correcting the subject absorption image. The correction of the ratio of the high frequency component in the BG absorption image I AB_BG (x, y) is performed, for example, by extracting the high frequency component H AB_BG (x, y) for each pixel of the BG absorption image I AB_BG (x, y). This is performed by multiplying the high frequency component by a coefficient a (x, y) for adjusting the component H AB_BG (x, y) to the subject absorption image I AB_SAMPLE (x, y).
Accordingly, the BG processing is performed after the ratio of the image unevenness with respect to the pixel value of the BG absorption image I AB_BG (x, y) is matched with the ratio of the image unevenness of the subject absorption image I AB_SAMPLE (x, y). It is possible to reduce image unevenness with high accuracy from AB_SAMPLE (x, y). Further, unlike the prior art, it is not necessary to interpose a member that causes an X-ray spectrum change equivalent to that of the subject at the time of BG imaging, so that it is possible to easily reduce image unevenness without trouble.
The coefficient a (x, y) is preferably set based on the pixel value of the absorption image I AB (x, y) after the conventional BG processing is generated. Thereby, it is possible to correct the BG absorption image I AB_BG (x, y) with high accuracy according to the thickness of the subject.
The coefficient a (x, y) is preferably set differently for the region corresponding to the bone region and the region corresponding to the soft region in the subject absorption image I AB_SAMPLE (x, y). Thereby, it becomes possible to perform correction suitable for each of the bone region and the soft region.
The coefficient a (x, y) is preferably set for each pixel in the BG absorption image I AB_BG (x, y) according to the value of the high-frequency component at that pixel. This makes it possible to accurately correct both streaky unevenness and granular unevenness.
In addition, the description in this embodiment mentioned above is a suitable example which concerns on this invention, and is not limited to this.
For example, in the above embodiment, the X-ray imaging apparatus using the Talbot-Lau interferometer that moves the multi slit 12 with respect to the first grating 14 and the second grating 15 during imaging has been described as an example. May be applied to an X-ray imaging apparatus using a Talbot-Lau interferometer of a type that moves either the multi-slit 12 or the first grating 14 or the second grating 15 or two of them. The present invention may also be applied to an X-ray imaging apparatus using a Talbot interferometer that moves either the first grating 14 or the second grating 15 relative to another grating.
In the above embodiment, the medical image system in which the multi-slit 12, the first grating 14, and the second grating 15 are one-dimensional gratings has been described as an example. However, the present invention uses a two-dimensional grating to form a two-dimensional image. It is also possible to apply to a medical image system that performs fringe scanning.
In the above embodiment, the case where three types of reconstructed images are generated has been described. However, the present invention can be applied to any medical image system that generates at least an absorption image.
In addition, the detailed configuration and detailed operation of each device constituting the medical image system can be changed as appropriate without departing from the spirit of the invention.
DESCRIPTION OF SYMBOLS 1 X-ray imaging apparatus 11 X-ray source 12 Multi-slit 12a Drive part 13 Subject stand 14 1st grating | lattice 15 2nd grating | lattice 16 X-ray detector 17 Holding | maintenance part 17a Buffer member 18 Main body part 181 Control part 182 Operation part 183 Display part 184 Communication unit 185 Storage unit 18a Drive unit 5 Controller 51 Control unit 52 Operation unit 53 Display unit 54 Communication unit 55 Storage unit
The image non-uniformity correction means, an image processing system for correcting the image unevenness of the object absorption image using the corrected subject without absorption image by the subject without absorption image correction means.
The image processing according to any one of claims 1 to 4 , wherein the image unevenness correction unit corrects image unevenness of the subject absorption image by dividing the subject absorption image by the corrected subject-free absorption image. System .
For each pixel of the non-subject absorption image, a high-frequency component is extracted by subtracting 1 from the value obtained by dividing the pixel value by the low-frequency component value obtained by blurring the non-subject absorption image. A non- subject absorption image correction unit that corrects only the ratio of the high-frequency component in the non-subject absorption image by multiplying the high-frequency component by a correction coefficient for matching the high-frequency component to the high-frequency component of the subject absorption image;
The subject-free absorbed image correction means sets a correction coefficient corresponding to a value of each pixel of an image obtained by correcting the subject-absorbed image using the subject-free absorbed image before correction to a corresponding pixel of the subject-free absorbed image. An image processing apparatus that performs the correction .
The subject-free absorbed image correction means is an image processing apparatus that performs the correction by setting different correction coefficients for a region corresponding to a bone region and a region corresponding to a soft region in the subject absorbed image .
The subject-absorbed absorption image correction means performs the correction by setting different correction coefficients for a region corresponding to a metal region and a region corresponding to a resin region in the subject absorption image .
The non-subject absorbed image correction means is an image processing apparatus that performs the correction for each pixel in the non-subject absorbed image by setting a correction coefficient according to the value of a high frequency component in the pixel .
The image processing according to any one of claims 6 to 9 , wherein the image unevenness correction unit corrects image unevenness of the subject absorption image by dividing the subject absorption image by the corrected subject-free absorption image. Equipment .
JP2015036006A 2015-02-26 2015-02-26 Image processing system and image processing apparatus Active JP6451400B2 (en)
JP2015036006A JP6451400B2 (en) 2015-02-26 2015-02-26 Image processing system and image processing apparatus
US15/008,284 US10076298B2 (en) 2015-02-26 2016-01-27 Image processing system and image processing device
JP2016154766A JP2016154766A (en) 2016-09-01
JP6451400B2 true JP6451400B2 (en) 2019-01-16
ID=56799032
JP2015036006A Active JP6451400B2 (en) 2015-02-26 2015-02-26 Image processing system and image processing apparatus
US (1) US10076298B2 (en)
JP (1) JP6451400B2 (en)
JP4377571B2 (en) * 2001-09-11 2009-12-02 株式会社東芝 X-ray flat panel detector, X-ray image diagnostic apparatus, and X-ray image correction method
AT524056T (en) * 2007-11-15 2011-09-15 Suisse Electronique Microtech Interferometer device and method
JP2010249533A (en) * 2009-04-10 2010-11-04 Canon Inc Source grating for talbot-lau-type interferometer
JP5302838B2 (en) * 2009-09-30 2013-10-02 富士フイルム株式会社 Image processing apparatus, image processing method, image processing program, and X-ray image photographing apparatus
RU2545319C2 (en) * 2009-12-10 2015-03-27 Конинклейке Филипс Электроникс Н.В. Phase-contrast image formation
US8989474B2 (en) * 2010-03-18 2015-03-24 Konica Minolta Medical & Graphic, Inc. X-ray image capturing system
JP5708652B2 (en) * 2010-08-31 2015-04-30 コニカミノルタ株式会社 X-ray imaging system
EP2633814A4 (en) * 2010-10-27 2014-03-19 Fujifilm Corp Radiography system and radiograph generation method
JP5475925B2 (en) * 2011-04-20 2014-04-16 富士フイルム株式会社 Radiation imaging apparatus and image processing method
US8767915B2 (en) * 2011-07-29 2014-07-01 The Johns Hopkins University Differential phase contrast X-ray imaging system and components
JP6076981B2 (en) * 2011-08-19 2017-02-08 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． Frequency dependent composite of X-ray images of different modalities
JP5857800B2 (en) * 2012-03-01 2016-02-10 コニカミノルタ株式会社 Joint imaging apparatus and imaging object fixing unit
EP2867655B1 (en) * 2012-06-27 2019-05-08 Koninklijke Philips N.V. Dark-field imaging
US8989347B2 (en) * 2012-12-19 2015-03-24 General Electric Company Image reconstruction method for differential phase contrast X-ray imaging
JP2014135989A (en) 2013-01-16 2014-07-28 Konica Minolta Inc Medical image system
WO2014167901A1 (en) * 2013-04-08 2014-10-16 コニカミノルタ株式会社 Diagnostic medical image system and method for introducing talbot imaging device system into diagnostic medical image system used for general imaging
JP6291812B2 (en) * 2013-11-29 2018-03-14 コニカミノルタ株式会社 Medical imaging system
JP6187298B2 (en) * 2014-02-14 2017-08-30 コニカミノルタ株式会社 X-ray imaging system and image processing method
US9330456B2 (en) * 2014-04-29 2016-05-03 General Electric Company Systems and methods for regularized Fourier analysis in x-ray phase contrast imaging
2015-02-26 JP JP2015036006A patent/JP6451400B2/en active Active
2016-01-27 US US15/008,284 patent/US10076298B2/en active Active
US10076298B2 (en) 2018-09-18
JP2016154766A (en) 2016-09-01
US20160253785A1 (en) 2016-09-01
Marshall 2006 A comparison between objective and subjective image quality measurements for a full field digital mammography system
CN101453954B (en) 2011-08-17 Radiation image pick-up device
US8824629B2 (en) 2014-09-02 Radiation imaging system and image processing method
JP4118535B2 (en) 2008-07-16 X-ray inspection equipment
Ref document number: 6451400