WO2012128335A1 - 医用画像表示システム - Google Patents
医用画像表示システム Download PDFInfo
- Publication number
- WO2012128335A1 WO2012128335A1 PCT/JP2012/057404 JP2012057404W WO2012128335A1 WO 2012128335 A1 WO2012128335 A1 WO 2012128335A1 JP 2012057404 W JP2012057404 W JP 2012057404W WO 2012128335 A1 WO2012128335 A1 WO 2012128335A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- image
- grating
- slit
- ray
- subject
- Prior art date
Links
- 238000003384 imaging method Methods 0.000 claims abstract description 158
- 238000010521 absorption reaction Methods 0.000 claims abstract description 122
- 238000012545 processing Methods 0.000 claims description 47
- 238000006243 chemical reaction Methods 0.000 claims description 13
- LFEUVBZXUFMACD-UHFFFAOYSA-H lead(2+);trioxido(oxo)-$l^{5}-arsane Chemical compound [Pb+2].[Pb+2].[Pb+2].[O-][As]([O-])([O-])=O.[O-][As]([O-])([O-])=O LFEUVBZXUFMACD-UHFFFAOYSA-H 0.000 abstract description 21
- 238000013399 early diagnosis Methods 0.000 abstract description 7
- 238000000034 method Methods 0.000 description 143
- 230000002159 abnormal effect Effects 0.000 description 94
- 238000012937 correction Methods 0.000 description 45
- 238000003745 diagnosis Methods 0.000 description 43
- 238000001514 detection method Methods 0.000 description 41
- 230000008569 process Effects 0.000 description 40
- 238000003860 storage Methods 0.000 description 29
- 238000004891 communication Methods 0.000 description 26
- 210000000481 breast Anatomy 0.000 description 16
- 208000004434 Calcinosis Diseases 0.000 description 14
- 210000001519 tissue Anatomy 0.000 description 14
- 230000008859 change Effects 0.000 description 11
- 206010028980 Neoplasm Diseases 0.000 description 10
- 206010006187 Breast cancer Diseases 0.000 description 8
- 208000026310 Breast neoplasm Diseases 0.000 description 8
- 230000000694 effects Effects 0.000 description 8
- 230000002308 calcification Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 239000000463 material Substances 0.000 description 7
- 230000002950 deficient Effects 0.000 description 6
- 230000003902 lesion Effects 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 210000000845 cartilage Anatomy 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 5
- BUGBHKTXTAQXES-UHFFFAOYSA-N Selenium Chemical compound [Se] BUGBHKTXTAQXES-UHFFFAOYSA-N 0.000 description 4
- 230000005856 abnormality Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 239000010408 film Substances 0.000 description 4
- 230000003211 malignant effect Effects 0.000 description 4
- 238000002360 preparation method Methods 0.000 description 4
- 229910052711 selenium Inorganic materials 0.000 description 4
- 239000011669 selenium Substances 0.000 description 4
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 3
- 201000011510 cancer Diseases 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 3
- 238000005520 cutting process Methods 0.000 description 3
- 210000003414 extremity Anatomy 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 238000000206 photolithography Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 239000010703 silicon Substances 0.000 description 3
- 229910052710 silicon Inorganic materials 0.000 description 3
- 210000002435 tendon Anatomy 0.000 description 3
- 206010043248 Tendon rupture Diseases 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000012984 biological imaging Methods 0.000 description 2
- 238000001574 biopsy Methods 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 210000005075 mammary gland Anatomy 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000003825 pressing Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 125000006850 spacer group Chemical group 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 210000000779 thoracic wall Anatomy 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 2
- 229910052721 tungsten Inorganic materials 0.000 description 2
- 239000010937 tungsten Substances 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 241000167854 Bourreria succulenta Species 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 125000002066 L-histidyl group Chemical group [H]N1C([H])=NC(C([H])([H])[C@](C(=O)[*])([H])N([H])[H])=C1[H] 0.000 description 1
- 206010065433 Ligament rupture Diseases 0.000 description 1
- ZOKXTWBITQBERF-UHFFFAOYSA-N Molybdenum Chemical compound [Mo] ZOKXTWBITQBERF-UHFFFAOYSA-N 0.000 description 1
- 206010047571 Visual impairment Diseases 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 210000001188 articular cartilage Anatomy 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 235000019693 cherries Nutrition 0.000 description 1
- 238000004195 computer-aided diagnosis Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 229920001971 elastomer Polymers 0.000 description 1
- 239000000806 elastomer Substances 0.000 description 1
- 238000005323 electroforming Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 210000004247 hand Anatomy 0.000 description 1
- 239000011810 insulating material Substances 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 239000011133 lead Substances 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000009607 mammography Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 229910052750 molybdenum Inorganic materials 0.000 description 1
- 239000011733 molybdenum Substances 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 238000007740 vapor deposition Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
- 210000000707 wrist Anatomy 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/06—Diaphragms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/483—Diagnostic techniques involving scattered radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/484—Diagnostic techniques involving phase contrast X-ray imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/502—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of breast, i.e. mammography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5247—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from an ionising-radiation diagnostic technique and a non-ionising radiation diagnostic technique, e.g. X-ray and ultrasound
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21K—TECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
- G21K2207/00—Particular details of imaging devices or methods using ionizing electromagnetic radiation such as X-rays or gamma rays
- G21K2207/005—Methods and devices obtaining contrast from non-absorbing interaction of the radiation with matter, e.g. phase contrast
Definitions
- the present invention relates to a medical image display system.
- phase contrast imaging is performed in which X-ray images with high visibility are obtained by edge enhancement using X-ray refraction during magnified imaging (see, for example, Patent Documents 1 and 2).
- the absorption contrast method is effective for photographing a subject with large X-ray absorption such as bone.
- the phase contrast method has a small X-ray absorption difference and can image breast tissue, articular cartilage, and soft tissue around the joint, which are difficult to appear as an image by the absorption contrast method. Application to diagnosis is expected.
- 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 arranges the second lattice at a position connecting the self-images, and measures interference fringes generated by slightly shifting the second lattice. If an object is placed in front of the second grating, the moire is disturbed.
- the object is placed in front of the first grating and irradiated with coherent X-rays. It is possible to obtain a reconstructed image of the subject by calculating a moire image.
- a Talbot-Lau interferometer has been proposed in which a multi-slit is installed between the X-ray source and the first grating to increase the X-ray irradiation dose (see, for example, Patent Document 6).
- a conventional Talbot-Lau interferometer takes a plurality of moire images at a constant interval while moving the first grating or the second grating (relatively moving both gratings). Is provided for the increase of
- the applicant of the present application is a system that enables scanning with high mechanical accuracy by moving the multi slit with respect to the first grating and the second grating in the Talbot-Lau interferometer and obtains a high-definition image.
- the applicant of the present application has applied for a system that can obtain a high-definition image in a Talbot-Lau interferometer (see Patent Document 8).
- a reconstructed image is created by a fringe scanning method from a plurality of moire images with a constant periodic interval obtained by a Talbot interferometer and a Talbot-low interferometer.
- a method of creating a reconstructed image from a single moire image using a Fourier transform method is also known (see, for example, Non-Patent Document 1).
- the reconstructed image obtained by the Fourier transform method has a lower spatial resolution than the fringe scanning method, but does not require a plurality of moire images unlike the fringe scanning method. Therefore, it is possible to shorten the shooting time and to suppress the influence of the body movement of the subject during a plurality of shootings. Further, since the scanning mechanism operation at the time of photographing is eliminated, there is no false image due to an error of the grating or multi-slit feeding mechanism.
- CAD Computer-Aided
- the CAD detection result is used as diagnosis support information
- the CAD detection result is used as the second opinion for the doctor's interpretation result based on the X-ray image (absorption image).
- the abnormal shadow candidate detection algorithm by increasing the threshold setting used to determine whether or not it is an abnormal shadow candidate, at least obviously positive candidates are detected, or by lowering the threshold setting
- the interpretation time varies greatly depending on the detection error, etc., the fatigue level of the doctor is different, and the diagnostic accuracy is not stable.
- An object of the present invention is to effectively use a reconstructed image created from a moire image generated by a Talbot interferometer or a Talbot-Lau interferometer to realize early diagnosis and improve diagnostic accuracy. .
- a medical image display system includes: An X-ray source that emits X-rays; A first grating and a second grating configured by arranging a plurality of slits in a direction orthogonal to the X-ray irradiation axis direction; Subject table, An X-ray detector that two-dimensionally arranges a conversion element that generates an electrical signal according to the irradiated X-ray, and reads the electrical signal generated by the conversion element as an image signal; A fringe scanning type imaging apparatus or a Fourier transform type imaging apparatus, An image processing unit that generates at least two of an X-ray absorption image, a differential phase image, and a small-angle scattered image of the subject based on the image signal of the subject imaged by any of the imaging devices; A display unit for displaying an image generated by the image processing unit; A control unit that controls display on the display unit of the image generated by the image processing unit; Is provided.
- control unit sequentially switches and displays at least two images generated by the image processing unit on the display unit every predetermined time.
- the fringe scanning imaging apparatus may be a Talbot-Lau interferometer having a multi-slit disposed in the vicinity of the X-ray source and moving the multi-slit relative to the first grating and the second grating. preferable.
- the present invention it is possible to realize early diagnosis and improve diagnostic accuracy by effectively using a reconstructed image created from a moire image generated by a Talbot interferometer or a Talbot-low interferometer. It becomes.
- FIG. 1 is a diagram illustrating a medical image display system (including a side view of an X-ray imaging apparatus) according to the present embodiment. It is a top view of a multi slit. It is the top view and side view of the state which hold
- FIG. 7A It is the top view which expanded and showed the holding
- FIG. 7A It is a figure which shows the state which hold
- FIG. 7A It is sectional drawing which shows the rotation tray which can rotate 1st grating
- FIG. 1 shows a medical image display system according to the first embodiment.
- the medical image display system includes an X-ray imaging apparatus 1 and a controller 5.
- the X-ray imaging apparatus 1 is an apparatus having a first imaging mode that functions as a fringe scanning imaging apparatus and a second imaging mode that functions as a Fourier transform imaging apparatus.
- the fringe scanning type imaging apparatus captures images in a plurality of steps with a Talbot-Lau interferometer for a reconstructed image by the fringe scanning method, and generates a plurality of moire images.
- the Fourier transform type imaging device captures images in one or two directions for a reconstructed image by the Fourier transform method, and generates one or two moire images.
- the configuration of the X-ray imaging apparatus 1 will be described taking an example of an apparatus that performs imaging using a finger as a subject, but is not limited thereto.
- 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)
- the distance between the multi-slit 12 and the first grating 14 is represented by d3 (mm)
- the distance between the first grating 14 and the second grating 15 is represented by d4 (mm).
- the distance d1 is preferably 5 to 500 (mm), more preferably 5 to 300 (mm).
- the distance d2 is preferably at least 3000 (mm) or less since the height of the radiology room is generally about 3 (m) or less.
- the distance d2 is preferably 400 to 2500 (mm), and 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 2800 (mm), 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
- 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 a C-shaped arm shape, and is attached to the main body portion 18 so as to be movable (up and down) 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 from the X-ray tube, and irradiates the X-rays in the z direction (gravity direction).
- 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.
- 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).
- the multi slit 12 is a diffraction grating, and a plurality of slits are arranged at predetermined intervals as shown in FIG. 2A.
- the plurality of slits are arranged in a direction (indicated by white arrows in FIG. 2A) orthogonal to the X-ray irradiation axis direction (z direction in FIG. 1).
- 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.
- 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. 2A, the slit period is defined as a period between adjacent slits.
- the width of the slit (the length of each slit in the slit arrangement direction) is 1 to 60 (%) of the slit period, and more preferably 10 to 40 (%).
- the height of the slit (the height in the z direction) is 1 to 500 ( ⁇ m), preferably 1 to 150 ( ⁇ m).
- the multi slit 12 is held by a holder 12b having a rack 12a as shown in FIG. 2B.
- the rack 12 a is provided in the slit arrangement direction of the multi slit 12.
- the rack 12a is engaged with a pinion 122c of the driving unit 122 described later, and moves the multi slit 12 held by the holder 12b in the slit arrangement direction according to the rotation (phase angle) of the pinion 122c.
- the X-ray imaging apparatus 1 is provided with a multi-slit rotating unit 121 and a driving unit 122.
- the multi-slit rotator 121 is a mechanism for rotating the multi-slit 12 around the X-ray irradiation axis according to the rotation of the first grating 14 and the second grating 15 around the X-ray irradiation axis.
- the drive unit 122 is a mechanism for moving the multi-slit 12 in the slit arrangement direction for photographing a plurality of moire images.
- FIG. 3 shows a plan view and a cross-sectional view taken along line AA ′ of the multi-slit rotating unit 121 and the driving unit 122.
- the multi-slit rotating unit 121 includes a motor unit 121a, a gear unit 121b, a gear unit 121c, a support unit 121d, and the like.
- the motor part 121a, the gear part 121b, and the gear part 121c are held by the holding part 17 via the support part 121d.
- the motor unit 121a is a pulse motor that can be switched to microstep driving, is driven in accordance with control from the control unit 181 (see FIG. 8), and drives the gear unit 121c to the X-ray irradiation axis (see FIG. 8) via the gear unit 121b. 3) (represented by a dashed line R in FIG. 3).
- the gear part 121c has an opening part 121e for mounting the multi slit 12 held by the holder 12b.
- the multi slit 12 mounted on the opening 121e can be rotated around the X-ray irradiation axis, and the slit arrangement direction of the multi slit 12 can be varied.
- the gear portion 121c does not have to be on the entire circumference, and the range indicated by the two-dot chain line in FIG. ).
- the opening 121e has a shape and size that allows the multi slit 12 held by the holder 12b to be fitted from above.
- the size w4 of the opening 121e in the slit arrangement direction is slightly larger than the size W2 of the holder 12b in the slit arrangement direction, and the multi-slit 12 can be slid in the slit arrangement direction.
- the size w3 in the direction orthogonal to the slit arrangement direction in the opening 121e is a dimension that allows precise fitting with the size W1 in the direction orthogonal to the slit arrangement direction in the holder 12b, and the holder 12b is attached to the opening 121e. Then, the rack 12a provided in the holder 12b is disposed outside the opening 121e so as to be engageable with a pinion 122c described later.
- the drive unit 122 includes a precision reduction gear that moves the multi slit 12 in the slit arrangement direction in units of 0.1 ⁇ m to several tens of ⁇ m according to the multi slit period.
- the drive unit 122 includes a motor unit 122a, a gear unit 122b, a pinion 122c, and the like, and is fixed to the gear unit 121c of the multi-slit rotating unit 121 by an L-shaped sheet metal (not shown).
- the multi slit 12 and the drive part 122 are rotated integrally.
- the motor unit 122a is driven in accordance with control from the control unit 181 and rotates the pinion 122c via the gear unit 122b.
- the pinion 122c engages with the rack 12a of the multi slit 12 and rotates to move the multi slit 12 in the slit arrangement direction.
- the subject table 13 is a table for placing a finger as a subject.
- the subject table 13 is preferably provided at a height at which the patient's elbow can be placed. In this way, by being configured to be placed up to the patient's elbow, the patient can have a comfortable posture, and the movement of the imaging part of the fingertip can be reduced during imaging for a relatively long time.
- the subject table 13 is provided with a subject holder 130 for fixing the subject.
- the subject holder 130 is detachable depending on the subject.
- the subject holder 130 is a plate-like member having an elliptical shape 131 such as a mouse that can be easily grasped by the palm.
- the palm shape is a gentle convex curved surface, and the subject is less likely to get tired by grasping the elliptical shape 131 with the palm.
- the downward movement of the subject can be suppressed.
- the X-ray complex refractive index of the subject holder 130 is not uniform for the X-ray dose reaching the X-ray detector 16.
- An inter-finger spacer 133 is preferably provided on the subject holder 130 in order to further stabilize the subject posture.
- the subject holder 130 is created according to the shape of the palm for each patient, and the subject holder 130 for the patient is magnetized on the subject table 13 at the time of photographing. It is preferable to attach by etc.
- the subject holder 130 Since the load from the arm to the wrist is supported by the subject table 13, the subject holder 130 only needs to be able to withstand the weight of the fingertip and the force pressed by the patient from above, and can be made of plastic that can be mass-produced at low cost. It is.
- the first grating 14 is a diffraction grating provided with a plurality of slits arranged in a direction orthogonal to the z direction, which is the X-ray irradiation axis direction, like the multi-slit 12.
- 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).
- the slit height (z-direction height) is a phase difference of ⁇ between the two materials forming the slit period, that is, the materials of the X-ray transmitting portion and the X-ray shielding portion. / 8 to 15 ⁇ ⁇ / 8.
- the height is preferably ⁇ / 4 to 3 ⁇ ⁇ / 4.
- the height of the slit is set to a height at which X-rays are sufficiently absorbed by the X-ray shielding part.
- the above-mentioned condition is that the first grating 14 is a ⁇ / 2 type grating, that is, the case where the phase difference due to the materials of the X-ray shielding part and the X-ray transmitting part of the first grating is in the vicinity of ⁇ / 2.
- the first lattice 14 may be a ⁇ -type, and a condition corresponding to the type of the lattice to be used may be calculated.
- the second grating 15 is a diffraction grating provided with a plurality of slits arranged in a direction orthogonal to the z direction, which is the X-ray irradiation axis direction, like the multi-slit 12.
- 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 period, and preferably 35 to 60 (%).
- the height of the slit is 1 to 100 ( ⁇ m).
- each of the first grating 14 and the second grating 15 has a grating plane perpendicular to the z direction (parallel in the xy plane), and the slit arrangement direction of the first grating 14 and the second grating 15.
- the slit arrangement direction is inclined at a predetermined angle in the xy plane, but both may be arranged in parallel.
- lattice 15 are disk shape.
- 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 tube 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)
- 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)
- FIG. 5 shows a plan view of the lattice rotation unit 210.
- FIG. 6 shows a plan view and a DD ′ sectional view of the grating rotating unit 210 in a state where the first grating 14 and the second grating 15 are mounted.
- the lattice rotation unit 210 includes a handle 211, a relative angle adjustment unit 213, a stopper 214, and the like on a rotation tray 212.
- the rotating tray 212 has an opening 212 a for holding the first lattice 14 and the second lattice 15.
- the first lattice 14 includes a circular lattice portion 140 in which a plurality of slits are arranged, and a first holder portion 141 and a second holder for attaching the lattice portion 140 to the opening 212a. And a holder part 142 (see FIG. 6).
- the first holder portion 141 is a member attached to the outer periphery of the lattice portion 140 and having the same radius (peripheral radius) as the opening portion 212a.
- the first holder portion 141 is fitted to the opening portion 212a when the first lattice 14 is attached.
- the second holder part 142 is a member that is attached to the outer side of the first holder part 141 and has a slightly larger radius (periphery radius) than the opening 212a.
- the second holder portion 142 has a part of the outer periphery that is gear processed.
- a protrusion 142 a is provided at a predetermined position on the outer periphery of the second holder portion 142.
- the second grid 15 includes a circular grid section 150 in which a plurality of slits are arranged, and a holder section 151 for attaching the grid section 150 to the opening 212a.
- the holder portion 151 is a disk-shaped member having a radius substantially the same as the radius of the opening 212a.
- the lattice part 150 is held on the upper surface of the central part of the holder part 151 (see FIG. 6).
- the second grid 15 is fitted into the bottom surface of the opening 212a.
- the first grid 14 is fitted into the opening 212 a from above the second grid 15. Accordingly, the first grating 14 and the second grating 15 are held on the rotating tray 212 in the state shown in FIG.
- the relative angle in the slit direction of the first grating 14 and the second grating 15 held in the opening 212a is adjusted by the relative angle adjusting unit 213 according to the imaging mode.
- the X-ray imaging apparatus 1 performs imaging in one or two directions for a first imaging mode in which imaging is performed in a plurality of steps for a reconstructed image by a fringe scanning method, and for a reconstructed image by a Fourier transform method. 2 shooting modes.
- the relative angle between the slit direction of the first grating 14 and the slit direction of the second grating 15 required in photographing for the fringe scanning method depends on the period of the second grating, the image size, and the number of fringes.
- the relative angle adjustment unit 213 can automatically adjust the relative angle between the first grating 14 and the second grating 15 according to the imaging mode set by the operation unit 182. It has become.
- the relative angle adjustment unit 213 includes a motor unit 213a, a first gear 213b, a second gear 213c, and a lever 213d.
- the motor unit 213a engages with the second gear 213c, and rotates the second gear 213c according to control from the control unit 181.
- the center of the second gear 213c is connected to the center of the first gear 213b via the lever 213d, and the circumference thereof is engaged with the first gear 213b.
- the first gear 213b rotates around the second gear 213c with the center of the second gear 213c as the rotation axis, and the first gear 14
- the first grating 14 can be rotated around the X-ray irradiation axis without engaging the gear portion of the two holder part 142 and rotating the second grating 15.
- the first grating 14 and the second grating 15 are formed when the protrusion 142a of the second holder 142 abuts against a stopper (convex protrusion) 214 provided on the rotating tray 212 at the time of factory shipment.
- the position of the stopper 214 and the relative angle between the first grating 14 and the second grating 15 are set so that the relative angle in the slit direction becomes an optimum relative angle in the first imaging mode (imaging mode for the fringe scanning method). It is adjusted in advance and attached to the opening 212a.
- the relative angle adjustment unit is controlled by the control unit 181 so that the relative angle between the first grating 14 and the second grating is optimal for the second shooting mode.
- the motor unit 213a employing the pulse motor 213 is driven (energization control). Accordingly, the first gear 213b rotates through the second gear 213c and engages with the gear portion of the second holder part 142, and the relative angle in the slit direction between the first grating 14 and the second grating 15 is the second angle.
- the second holder 142 is rotated so as to be optimal for the shooting mode.
- the energization state of the pulse motor After that, by changing the energization state of the pulse motor, the energization state (less than 50% of the rated current at the time of driving, etc.) to the extent that the motor self-holding force (excitation force) that overcomes the spring force described later is exhibited.
- the second holder part 142 can be maintained in this phase. Since the rotation angle at this time is as small as about 1 degree, first, the second holder part 142 is rotated counterclockwise by the pulse motor of the motor part 213a until the projection part 142a reaches the reference position 215.
- the second holder part 142 is rotated by microstep driving by switching the rotation direction of the second holder part 142 clockwise.
- the second holder part 142 is biased by a spring (not shown).
- the projecting part 142a is caused by the biasing force of the spring. Returns to the position of the stopper 214. That is, the first grating 14 and the second grating 15 return to the optimum relative angle for the first imaging mode.
- the relative angle between the first grating 14 and the second grating 15 is adjusted to an angle corresponding to the photographing mode.
- the grid rotating unit 210 can also integrally rotate the first grid 14 and the second grid 15 whose relative angles are adjusted around the X-ray irradiation axis (indicated by a dotted line R in FIG. 6) with respect to the subject. it can.
- the structure extending linearly in parallel with the slit direction of the first grating 14 and the second grating 15 can be clearly photographed. There is a characteristic that it cannot.
- the grating rotating unit 210 rotates the first grating 14 and the second grating 15 integrally around the X-ray irradiation axis while maintaining the relative angle by the following mechanism, so that the object is focused on the arrangement direction of the structure to be noted.
- the angle in the slit direction of the first grating 14 and the second grating 15 can be adjusted.
- the rotary tray 212 is provided with the handle 211 as described above.
- the handle 211 is a protrusion for an operator such as a radiographer to manually rotate the rotating tray 212 around the X-ray irradiation axis (indicated by a dotted line R in FIG. 6).
- the rotating tray 212 has recesses 212b to 212e for fixing the rotation angle of the rotating tray 212.
- the recesses 212b to 212e are positions at a predetermined rotation angle (here, the position where the recess 212b faces the ball of the tray fixing member 171b is a 0 ° position) (here, the position where the recess 212b faces the ball of the tray fixing member 171b).
- Angle detection sensors SE1 to SE4 are provided in the recesses 212b to 212e, respectively, and detect that they are engaged with the tray fixing member 171b and output the detection signal to the control unit 181.
- the rotation tray 212 since the rotation tray 212 is manually rotated, it is not necessary to provide an electric cord or the like for integrally rotating the first grid 14 and the second grid 15 in a range touched by the patient, thereby ensuring safety. be able to.
- the position (angle) of the first grid 14 and the second grid 15 when the rotating tray 212 is set to 0 ° is set as the home position.
- the position (angle) at which the slit direction of the first grating 14 and the slit direction of the multi-slit 12 are parallel when the first grating 14 and the second grating 15 are at the home position is defined as the home position of the multi-slit 12.
- FIG. 7A is an enlarged plan view showing the holding portion 171 of the lattice rotating unit 210 in the holding unit 17, and FIG. 7B is a cross-sectional view taken along the line EE ′ in FIG. 7A.
- FIG. 7C is a diagram illustrating a state in which the lattice rotation unit 210 is held by the holding unit 17.
- the holding portion 171 has a size that fits precisely with the rotating tray 212 of the lattice rotating unit 210, and has an opening 171a that rotatably holds the rotating tray 212, and the rotating tray 212.
- a tray fixing member 171b for fixing the rotation angle.
- the space between the bottom of the opening 171a and the mounting portion of the X-ray detector 16 is preferably hollow or made of aluminum or carbon having a high X-ray transmittance so as not to prevent the transmission of X-rays.
- the tray fixing member 171b When the tray fixing member 171b is positioned so that any of the recesses 212b to 212e faces the tray fixing member 171b, the tray engaging member 171b guides the ball in the direction of the arrow in FIGS. 7A and 7B.
- a slide guide (a guide of a pressing spring) (not shown) is used.
- the slide guide of the tray fixing member 171b causes the ball to engage with the facing recess and
- the angle detection sensor (any of SE1 to SE4) provided detects the engagement of the ball and outputs a detection signal to the control unit 181.
- the control unit 181 can detect the rotation angle of the rotating tray 212, that is, the rotation angles of the first grating 14 and the second grating 15.
- a mounting portion 212f of the X-ray detector 16 is provided below the opening 212a of the rotating tray 212, and the first grating 14, the second grating 15, and the X-ray detector 16 are integrated. It may be possible to rotate as follows. In this way, since the vertical and horizontal sharpness anisotropy of the X-ray detector 16 is not affected, the vertical and horizontal sharpnesses of the reconstructed image are rotated by the first and second gratings 14 and 15. It can be generally constant regardless of the angle.
- 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).
- 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 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.
- an FPD Felat Panel Detector
- the FPD includes an indirect conversion type in which X-rays are converted into electric signals by a photoelectric conversion element via a scintillator, and a direct conversion type in which X-rays are directly converted into electric signals. Any of these may be used.
- photoelectric conversion elements are two-dimensionally arranged with TFTs (thin film transistors) under a scintillator plate such as CsI or Gd 2 O 2 to constitute each pixel.
- TFTs thin film transistors
- 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.
- an amorphous selenium film having a film thickness of 100 to 1000 ( ⁇ m) is formed on glass by thermal vapor deposition of amorphous selenium, and the amorphous selenium film and electrodes are arranged on a two-dimensionally arranged TFT array. Vapor deposited.
- 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.
- imaging means such as a CCD (Charge Coupled Device) or an X-ray camera may be used as the X-ray detector 16.
- the FPD is reset to remove unnecessary charges remaining after the previous photographing (reading). Thereafter, charges are accumulated at the timing when the X-ray irradiation starts, and the charges accumulated at the timing when the X-ray irradiation ends are read as an image signal. Note that dark reading for offset correction is performed immediately after resetting or after reading an image signal.
- the main body 18 includes a control unit 181, an operation unit 182, a display unit 183, a communication unit 184, and a storage unit 185.
- the control unit 181 includes a CPU (Central Processing Unit), a RAM (Random Access Memory), and the like, and controls each unit of the X-ray imaging apparatus 1 in cooperation with a program stored in the storage unit 185 and various types. Execute the process. For example, the control unit 181 executes various processes including an imaging control process described later.
- the operation unit 182 includes a touch panel configured integrally with the display of the display unit 183 in addition to a key group used for input operations such as an exposure switch and an imaging condition, and generates an operation signal corresponding to these operations to generate a control unit. It outputs to 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 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 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 functions as an image processing unit that creates a diagnostic reconstructed image using the moire image obtained by the X-ray imaging apparatus 1.
- 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), a RAM (Random Access Memory), and the like, and creates and displays a reconstructed image by a fringe scanning method to be described later in cooperation with a program stored in the storage unit 55.
- Various processes including a process and a reconstructed image creation / display process by a Fourier transform method are executed.
- the control unit 51 performs the first imaging mode or the second in the X-ray imaging apparatus 1.
- At least two of the X-ray absorption image, the differential phase image, and the small angle scattering image are created based on the moire image obtained by the imaging in the imaging mode, and the display unit 53 controls the created image.
- the control unit 51 functions as an image processing unit and a control 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
- the display unit 53 includes, for example, a monitor such as a CRT (Cathode Ray Tube) or an LCD (Liquid Crystal Display), and the operation screen and the operation status of the X-ray imaging apparatus 1 according to the display control of the control unit 51.
- a monitor such as a CRT (Cathode Ray Tube) or an LCD (Liquid Crystal Display)
- the created reconstructed image is displayed.
- 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 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.
- the storage unit 55 stores imaging order information indicating an order reserved by a RIS, HIS or the like or a reservation device (not shown).
- the imaging order information is information such as a patient name, an imaging region, and an imaging mode.
- the storage unit 55 also displays setting information set by the operation unit 52, for example, the type of image displayed in step S22 in FIG. 14 or the switching timing of image display in step S22 in FIG. Is stored in association with.
- the storage unit 55 stores the moire image obtained by the X-ray detector 16 and the diagnostic reconstructed image created based on the moire image in association with the imaging order information.
- the storage unit 55 displays a reference image (details will be described later) indicating a typical case of a lesion as a lesion name and an image type (a fringe scanning method or a Fourier transform method, an absorption image or a differential phase image or a small angle scattered image), and the like. Store in association with each other. Further, 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 controller 5 when a list display of imaging order information is instructed by operating the operation unit 52, the imaging unit information is read from the storage unit 55 by the control unit 51 and displayed on the display unit 53.
- radiographing order information is designated by the operation unit 52
- setting information of radiographing conditions (including radiographing mode) according to the designated radiographing order information, an instruction for warming up the X-ray source 11, and the like are transmitted by the communication unit 54. It is transmitted to the line imaging apparatus 1.
- an imaging mode is set in the X-ray imaging apparatus 1. That is, the controller 5 functions as a setting unit that sets the shooting mode.
- the control unit 51 starts from the sleep state for preventing the internal battery consumption to the imaging ready state.
- the communication unit 184 receives imaging condition setting information from the controller 5, preparation for X-ray imaging is executed.
- An X-ray imaging method (imaging method in the first imaging mode) using the Talbot-Lau interferometer of the X-ray imaging apparatus 1 will be described.
- FIG. 10 when the X-rays irradiated 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 in parallel at a position connecting the self-images, and the second grating 15 is slightly tilted from a position parallel to the grating direction of the first grating 14, and therefore the second grating 15.
- a moire image M is obtained by the X-rays transmitted through.
- the phase of the X-ray is shifted by the subject H, so that the interference fringes on the moire image M are bordered on the edge of the subject H as shown in FIG. 10. Disturbed.
- the disturbance of the interference fringes can be detected by processing the moire image M, and the subject image can be imaged. This is the principle of the Talbot interferometer and Talbot low interferometer.
- 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.
- 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 images, the first grating 14 or the second grating 15 is used. It was moved relative. 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 images having a constant cycle interval are obtained. Note that when a moire image is obtained in the second shooting mode, the multi-slit 12 is not moved, and shooting is performed once or twice while rotating the subject and the slit direction by 90 degrees.
- FIG. 11 is a flowchart showing an imaging control process executed by the control unit 181 of the X-ray imaging apparatus 1.
- the imaging control process is executed by the cooperation of programs stored in the control unit 181 and the storage unit 185.
- Step S1 based on the setting information received from the controller 5, it is determined which of the first imaging mode (for fringe scanning method) or the second imaging mode (for Fourier transform method) is set.
- step S2 If it is determined that the first shooting mode is set (step S1; first shooting mode), the first shooting mode process is executed (step S2).
- step S1; second shooting mode the second shooting mode process is executed (step S3).
- FIG. 12 is a flowchart showing a first imaging mode process executed by the control unit 181 of the X-ray imaging apparatus 1 in step S2 of FIG.
- the first shooting mode process is executed in cooperation with a program stored in the control unit 181 and the storage unit 185.
- the X-ray imaging method using the Talbot-Lau interferometer is used for X-ray imaging in the first imaging mode, and the fringe scanning method is used for reconstruction of the subject image.
- the drive unit 122 is driven and stopped by the control of the control unit 181, so that the multi-slit 12 is moved by a plurality of steps at regular intervals, and imaging is performed for each step. An image is obtained.
- the number of steps 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)).
- the X-ray source 11 is switched to the warm-up state by the control unit 181 (step S101).
- the relative angle adjusting unit 213 of the grating rotating unit 210 is controlled so that the relative angle between the first grating 14 and the second grating 15 is optimal for the first imaging mode (the protrusion 142a contacts the stopper 214).
- the first grid 14 is rotated (so that it is in position). Thereby, the relative angle of the 1st grating
- the first grid 14 and the second grid 15 are integrally rotated in accordance with the operation of the operator, and the slit directions of the first grid 14 and the second grid 15 with respect to the subject are set (step S103). That is, an operator such as a photographic engineer rotates the handle 211 of the grid rotation unit 210, and the first grid 14 and the second grid 15 according to the arrangement direction of the structure to be noted of the subject placed on the subject table 13. Set the slit direction.
- a detection signal is output from any of the angle detection sensors SE1 to SE4 to the control unit 181 and the control is performed.
- the unit 181 the rotation angle from the home position of the rotary tray 212 (that is, the first grid 14 and the second grid 15) of the grid rotating unit 210 corresponding to the set slit direction is acquired.
- the motor unit 121a of the multi-slit rotating unit 121 is controlled by a pulse according to the rotation angle of the first grating 14 and the second grating 15, and the multi-slit according to the rotation angle of the first grating 14 and the second grating 15. 12 is rotated (step S104).
- the pulse motor of the motor unit 121a is controlled, and the rotation angle of the multi-slit 12 from the home position is rapidly close to the rotation angle of the rotation tray 212 (for example, about 29 ° when the rotation tray 212 is set to 30 °). It is rotated.
- step S105 the motor unit 121a is switched to microstep precision control, and shooting is performed at a plurality of rotation angles while rotating the multi-slit 12 little by little, and a plurality of moire images for adjustment are generated (step S105).
- the rotation tray 212 is set to 30 °
- the multi-slit 12 is set to three rotation angles of 29.5 °, 30 °, and 30.5 °, and low-dose X-rays are emitted from the X-ray source 11. Irradiated and photographed. Thereby, three moire images for adjustment are acquired.
- shooting is performed without placing the subject on the subject table 13.
- the plurality of adjustment moire images that have been photographed are displayed side by side on the display unit 183 in association with the rotation angle of the multi slit 12 (step S106).
- step S103 the rotation tray 212 is rotated.
- the first grating 14 and the second grating 15 are rotated while maintaining the relative angle.
- interference fringes that is, moire
- the relative angle between the multi slit 12 and the first and second gratings 14 and 15, that is, the rotary tray 212 on which these are placed In general, the smaller the relative angle between the multi-slit 12 and the first grating 14, the more moire images with the clearer fringes are obtained.
- the multi slit 12 is disposed in the vicinity of the X-ray source 11 which is a heat generating portion, it is easily affected by heat. Therefore, in consideration of deformation of the multi-slit 12, etc., not only the multi-slit 12 is rotated by the same angle as the rotary tray 212, but also the motor unit 121a is micro-step driven to perform fine adjustment in steps S105 to S108. It is valid.
- the operator observes the moire image displayed on the display unit 183 in step S106, and selects the rotation angle at which the interference fringes are clear as the rotation angle used for photographing.
- the sharpness of the interference fringes is observed by the operator's visual observation, the sharpness indicating the degree of the sharpness of the interference fringes is set to the maximum value in the sine curve (see FIG. 17) described later, and the minimum value. In the case of MIN, it can be expressed by the following formula.
- a rotation angle that automatically reaches a maximum value may be set by a program instead of an operator.
- step S107 When the rotation angle of the multi slit 12 is input by the operation unit 182 (step S107; YES), the motor unit 121a is re-driven so that the rotation angle from the home position of the multi slit 12 becomes the input rotation angle. Then, the position of the multi slit 12 is finely adjusted (step S108).
- step S109 After the rotation angle of the multi-slit 12 is adjusted, when the subject is placed on the subject table 13 and the exposure switch is turned on by the operator (step S109; YES), the multi-slit 12 is arranged in the slit arrangement direction by the drive unit 122. And a plurality of steps of photographing are executed, and a plurality of moire images with a subject are generated (step S110).
- X-ray irradiation by the X-ray source 11 is started with the multi-slit 12 stopped.
- the X-ray detector 16 accumulates charges in accordance with the timing of X-ray irradiation, and reads the accumulated charges as image signals in accordance with the timing of X-ray irradiation stop.
- the drive unit 122 is activated by the control of the control unit 181 and the movement of the multi-slit 12 is started.
- the driving unit 122 is stopped to stop the movement of the multi-slit, and the next step photographing is performed.
- 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.
- the read image signal is output to the main body 18 as a moire image.
- the slit period of the multi-slit 12 is 22.8 ( ⁇ m), and five-step shooting 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.
- the slit period of the second grating 15 is relatively small and the movement amount of each step is small, but the slit period of the multi slit 12 is the first. It is relatively larger than the two grids 15, and the movement amount of each step is also large.
- the amount of movement of the second grating 15 with a slit period of 5.3 ( ⁇ m) per step is 1.06 ( ⁇ m)
- the amount of movement of the multi-slit 12 with a slit period of 22.8 ( ⁇ m) is It is 4.56 ( ⁇ m), about four times as large.
- a moving pulse motor (drive source)
- the method of moving is smaller. This indicates that it is easy to obtain a moire image along a sine curve, which will be described later, and that a high-definition reconstructed image can be obtained even when the activation and the stop are repeated.
- the image based on the conventional method is sufficiently suitable for diagnosis, the accuracy of the entire drive transmission system including the motor (drive source) (particularly the start characteristics and stop characteristics) is relaxed, and the components of the drive transmission system are reduced. This shows that the cost can be reduced.
- the moire image of each step is transmitted from the communication unit 184 of the main body unit 18 to the controller 5 (step S111).
- a moire image with a subject is transmitted from the main body 18 to the controller 5 one by one every time photographing of each step is completed.
- dark reading is performed in the X-ray detector 16 to obtain a dark image (offset correction data) for correcting image data with a subject (step S112).
- the dark reading is performed at least once.
- the average value may be acquired as a dark image by performing multiple dark readings.
- the dark image is transmitted from the communication unit 184 to the controller 5 (step S113).
- the offset correction data based on the dark reading is commonly used for correcting each moire image signal.
- the acquisition of the dark image may be performed by performing dark reading of the corresponding step after generating the moire image of each step and generating offset correction data dedicated to each step.
- step S114 the operator enters an ON switch waiting state for the exposure switch.
- the operator removes the subject from the subject table 13 and retracts the patient so that a moire image without the subject can be created.
- the exposure switch is pressed.
- step S114 When the exposure switch is pressed (step S114; YES), the multi-slit 12 is moved in the slit arrangement direction by the driving unit 122, and shooting in a plurality of steps is performed without a subject, and a plurality of moire images without a subject are obtained. It is generated (step S115).
- the moire image of each step is transmitted from the communication unit 184 of the main body unit 18 to the controller 5 (step S116).
- a moire image without a subject is transmitted from the main body 18 to the controller 5 one by one by the communication unit 184 every time photographing of each step is completed.
- dark reading is performed in the X-ray detector 16, and a dark image without a subject is acquired (step S117).
- the dark reading is performed at least once.
- the average value may be acquired as a dark image by performing multiple dark readings.
- the dark image is transmitted from the communication unit 184 to the controller 5 (step S118), and a series of imaging for one imaging order is completed.
- the acquisition of the dark image may be performed by performing dark reading of the corresponding step after generating the moire image of each step and generating offset correction data dedicated to each step.
- the controller 5 when the moire image is received by the communication unit 54, the received moire image is stored in the storage unit 55 in association with the shooting order information specified at the start of shooting.
- FIG. 13 is a flowchart showing a second imaging mode process executed by the control unit 181 of the X-ray imaging apparatus 1 in step S3 of FIG.
- the second shooting mode process is executed by the cooperation of the program stored in the control unit 181 and the storage unit 185.
- the X-ray source 11 is first switched to the warm-up state by the control unit 181 (step S201).
- the relative angle adjustment unit 213 of the grid rotation unit 210 is controlled so that the relative angle between the first grid 14 and the second grid 15 is optimal for the second imaging mode (the projection 142a is at a predetermined angle from the home position). Adjustment is performed so that the rotation position is reached (step S202).
- step S203 to step S208 is performed.
- the processing in steps S203 to S208 is the same as that described in steps S103 to 108 in FIG.
- step S210 When the subject is placed on the subject table 13 and the exposure switch is turned on by the operator (step S209; YES), shooting is performed and a moire image with the subject is generated (step S210). That is, radiation is emitted from the X-ray source 11 and reading is performed by the X-ray detector 16. In the second shooting mode, only one image is shot without moving the multi slit 12 while the drive unit 122 is stopped. When shooting is completed, a moire image obtained by shooting is transmitted from the communication unit 184 of the main body 18 to the controller 5 (step S211).
- dark reading is performed in the X-ray detector 16, and a dark image (offset correction data) for correcting image data with a subject is acquired (step S212).
- the dark reading is performed at least once. Alternatively, the average value may be acquired as a dark image by performing multiple dark readings.
- the dark image is transmitted from the communication unit 184 to the controller 5 (step S213).
- the offset correction data based on the dark reading is commonly used for correcting the moire image signal.
- step S214 the operator enters an ON switch waiting state for the exposure switch.
- the operator removes the subject from the subject table 13 and retracts the patient so that a moire image without the subject can be created.
- the exposure switch is pressed.
- step S214 When the exposure switch is pressed (step S214; YES), shooting is performed without a subject, and a moire image without a subject is generated (step S215).
- step S215 as in step S210, only one image is taken without moving the multi-slit 12 while the driving unit 122 is stopped.
- a moire image is transmitted from the communication unit 184 of the main body unit 18 to the controller 5 (step S216).
- dark reading is performed in the X-ray detector 16, and a dark image without a subject is acquired (step S217).
- the dark reading is performed at least once.
- the average value may be acquired as a dark image by performing multiple dark readings.
- the dark image is transmitted from the communication unit 184 to the controller 5 (step S218), and a series of shooting for one shooting order is completed.
- control unit 51 of the controller 5 when a moire image is received by the communication unit 54, when the shooting mode set in the shooting order information currently being processed is the first shooting mode, fringe scanning is performed. Reconstructed image creation / display processing by the method is executed, and in the second imaging mode, reconstructed image creation / display processing by the Fourier transform method is executed.
- FIG. 14 is a flowchart showing the reconstructed image creation / display process by the fringe scanning method executed by the control unit 51.
- the reconstructed image creation / display process by the fringe scanning method is executed in cooperation with the control unit 51 and a program stored in the storage unit 55.
- steps S11 to S13 correction processing for correcting variation of each pixel of the X-ray detector 16 is executed for each of a plurality of moire images with a subject and a plurality of moire images without a subject. Specifically, an offset correction process (step S11), a gain correction process (step S12), and a defective pixel correction process (step S13) are executed.
- step S11 offset correction is performed on each moire image with a subject based on the dark image for correcting the image data with a subject. Based on the dark image for correcting the image data without a subject, an offset process is performed on each moire image without the subject.
- step S12 gain correction data corresponding to the X-ray detector 16 used for imaging is read from the storage unit 55, and gain correction is performed on each moire image based on the read gain correction data.
- the defective pixel map data indicating the defective pixel position
- the position indicated by the defective pixel position map in each moire image is read. Pixel values (signal values) are interpolated and calculated by surrounding pixels.
- step S14 X-ray intensity fluctuation correction (trend correction) is performed between the plurality of moire images (step S14).
- X-ray intensity fluctuation correction trend correction
- one reconstructed image is created based on a plurality of moire images. For this reason, if there is fluctuation (variation) in the intensity of X-rays irradiated in capturing each moiré image, an elaborate reconstructed image cannot be obtained, and fine signal changes may be overlooked. Therefore, in step S14, processing for correcting a signal value difference due to X-ray intensity fluctuations at the time of imaging in a plurality of moire images is performed.
- a correction method using a signal value of a predetermined pixel of each moire image, a signal value difference in a predetermined direction of the X-ray detector 16 between each moire image is corrected ( Any one of a method for correcting one-dimensionally and a method for correcting a signal value difference in a two-dimensional direction between each moire image (two-dimensional correction) may be used.
- a direct X outside the moire image area (subject placement area) 161 of the X-ray detector 16 is used for each of a plurality of moire images.
- a signal value of a pixel at a predetermined position P corresponding to the line area is acquired.
- the first moire image (for example, the first moire image taken with a subject) is normalized by the average signal value of the pixels at the acquired position P for the second and subsequent images, and the normalized position
- correction coefficients for the second and subsequent moire images are calculated.
- the X-ray intensity fluctuation is corrected by multiplying the second and subsequent moire images by a correction coefficient.
- detection means such as a sensor for detecting the X-ray irradiation amount is provided on the back side of the X-ray detector 16, and based on the X-ray irradiation amount at the time of capturing each moire image output from the detection means, It is also possible to correct the signal value difference caused by the X-ray intensity fluctuation at the time of imaging.
- an average signal value of pixels in a predetermined row L1 (a row indicates a reading line direction in the X-ray detector 16) is calculated for each of a plurality of moire images.
- the first moire image is normalized by the average signal value of the second and subsequent pixels, and two images are obtained based on the signal value of each pixel in the normalized row L1 and the second and subsequent rows L1.
- a correction coefficient in the row direction of each moire image after the eye is calculated.
- the X-ray intensity fluctuation in the row direction is corrected by multiplying the second and subsequent moire images by a correction coefficient corresponding to the position in the row direction.
- the fluctuation of the X-ray intensity in the one-dimensional direction between each imaging can be easily corrected.
- the X-ray intensity fluctuation in the reading line direction of the X-ray detector 16 and the like caused by this difference It can be corrected.
- the two-dimensional correction In the two-dimensional correction, first, for each of a plurality of moire images, an average of pixels in each of a predetermined row L1 and column L2 (the column indicates a direction orthogonal to the reading line direction in the X-ray detector 16). A signal value is calculated. Next, the first moiré image is normalized by the average signal value of the pixels in the second and subsequent rows L1, and based on the signal values of the pixels in the normalized row L1 and the second and subsequent rows L1. The correction coefficient in the row direction of each of the second and subsequent moire images is calculated.
- the first moire image is normalized by the average signal value of the pixels in the second and subsequent columns L2, and is based on the signal values of the respective pixels in the normalized column L2 and the second and subsequent columns L2.
- the correction coefficients in the column direction of the second and subsequent moire images are calculated.
- the correction coefficients for the pixels in the second and subsequent moire images are calculated by multiplying the correction coefficients in the row direction and the column direction.
- the correction coefficients for the pixels in the second and subsequent moire images are calculated by multiplying the correction coefficients in the row direction and the column direction.
- by multiplying each pixel by a correction coefficient in the row direction and the column direction fluctuations in X-ray intensity in the two-dimensional direction are corrected.
- the fluctuation of the X-ray intensity in the two-dimensional direction between each imaging can be easily corrected.
- step S15 the analysis of the moire image is performed (step S15), and it is determined whether or not it can be used to create a reconstructed image (step S16).
- the multi-slit 12 can be moved at a constant feed amount with ideal feed accuracy, five moire images corresponding to one slit period of the multi-slit 12 can be obtained in five steps as shown in FIG. Since the moire image of each step is a result of stripe scanning at a constant cycle interval of 0.2 cycles, when attention is paid to any one pixel of each moire image, the X-ray relative intensity obtained by normalizing the signal value Draws a sine curve as shown in FIG. Therefore, the controller 5 obtains the X-ray relative intensity by paying attention to the pixel having the moire image obtained in each step.
- the X-ray relative intensity obtained from each moiré image forms a sine curve as shown in FIG. 17, it is determined that a moiré image having a constant periodic interval is obtained and can be used to create a reconstructed image. be able to.
- the sine curve shape depends on the opening width of the multi-slit 12, the period of the first grating 14 and the second grating 15, and the distance between the gratings of the first and second gratings. In the case of coherent light, it has a triangular wave shape, but the X-rays act as quasi-coherent light due to the multi-slit effect, thereby drawing a sine curve.
- the analysis in step S15 is performed for each of the moire image with the subject and the moire image without the subject.
- step S16 If there is a moire image in which a sine curve cannot be formed in the moire image at each step, it is determined that it cannot be used to create a reconstructed image (step S16; NO), and an instruction is given to change the shooting timing and reshoot.
- Control information is transmitted from the controller 5 to the X-ray imaging apparatus 1 (step S17). For example, as shown in FIG. 17, if the third step is originally 0.4 cycle, and the cycle is shifted and a moiré image of 0.35 cycle is obtained, the feeding accuracy of the drive unit 122 is reduced. This is considered to be caused (for example, noise superimposition on the drive pulse of the pulse motor). Therefore, it suffices to instruct to re-shoot only the third step by advancing the shooting timing by 0.05 cycles.
- the imaging timing is adjusted according to the control information, and re-imaging is executed.
- a reconstructed image with a subject and no subject are used by using a plurality of moire images with and without a subject, respectively.
- step S21 using the reconstructed image without the subject, correction processing for removing the phase of the interference fringes and removing image unevenness (artifact) is performed from the reconstructed image with the subject (step S21).
- step S21 the X-ray dose distribution unevenness due to the slit direction change of the multi-slit 12 and the first and second gratings 14 and 15 at the time of imaging, the dose distribution unevenness due to the manufacturing variation of the slit.
- processing for removing image unevenness (artifact) including mainly unevenness due to the image of the subject holder 130 in the image is included.
- the signal value of the corresponding differential phase image without the subject is calculated from the signal value of each pixel of the differential phase image with the subject.
- Subtraction processing is performed (public literature (A); Timm Weitkamp, Ana Diazand, Christian David, franz Pfeiffer and Marco Stampanoni, Peter Cloetens and Eric Ziegler, X-ray Phase Imaging with a grating interferometer, OPTICS EXPRESS, Vol. 13, No.
- FIGS. 18A to 18C show an example of a reconstructed image created by the fringe scanning method based on a moire image obtained by shooting a cherry as a subject.
- 18A is an absorption image
- FIG. 18B is a differential phase image
- FIG. 18C is a small angle scattering image.
- the absorption image has a characteristic of representing a large structural change of the subject.
- the differential phase image has a feature that it represents a phase change of the tissue edge of the subject.
- the small-angle scattered image has a characteristic of representing scattering in the tissue of the subject.
- step S22 When the processing in step 21 is completed, the created reconstructed image is displayed on the display unit 53 (step S22).
- the display mode of the reconstructed image in step S22 will be described later.
- FIG. 19 is a flowchart showing a reconstructed image creation / display process by the Fourier transform method executed by the control unit 51.
- the reconstructed image creation / display process by the Fourier transform method is executed in cooperation with the control unit 51 and the program stored in the storage unit 55.
- steps S31 to S33 a correction process for correcting variations in each pixel of the X-ray detector 16 is executed for each of a plurality of moire images with a subject and a plurality of moire images without a subject. Specifically, offset correction processing (step S31), gain correction processing (step S32), and defective pixel correction processing (step S33) are executed. The contents of each process are the same as those described in steps S11 to S13 in FIG.
- X-ray intensity fluctuation correction (trend correction) between the moire image with the subject and the moire image without the subject is performed (step S34).
- the specific processing content of the X-ray intensity fluctuation correction is the same as that described in step S14 in FIG.
- a reconstructed image of the subject is created by the Fourier transform method.
- the creation of the reconstructed image by the Fourier transform method can be performed by a known method (see Non-Patent Document 1).
- FIG. 20A shows an example of a moiré image with a subject photographed in the second photographing mode.
- H1 is a magic pen and H2 is a USB memory.
- FIG. 20B shows the result of two-dimensional Fourier transform of the moire image of FIG. 20A.
- FIG. 21A shows an example of a moiré image without a subject photographed in the second photographing mode.
- FIG. 21B shows the result of two-dimensional Fourier transform of the moire image of FIG. 21A. Since the calculation result after the Fourier transform is a complex number, in FIG. 20B and FIG.
- the norm (amplitude) of the real part and the imaginary part is displayed.
- a low frequency component referred to as a 0th order component
- a component near the interference fringe frequency referred to as a primary component
- a high frequency component is obtained side by side.
- the direction in which the zeroth-order component and the first-order component are arranged is related to the direction of stripes in the moire image, and is substantially perpendicular to the direction of stripes in the moire image.
- the fringes of the moiré image for Fourier transform are fine horizontal stripes as compared with A2 in FIG.
- An image obtained by performing a Fourier transform on a moiré image for Fourier transform is an image in which a zeroth-order component and a first-order component are arranged vertically as indicated by A4 in FIG.
- the moire image for the fringe scanning method (second to the first grating 14).
- the stripes of the image obtained by slightly tilting the grid are vertical as shown by B2 in FIG.
- the stripes of the Fourier transform moire image are finer vertical stripes than B2 of FIG.
- An image obtained by performing a Fourier transform on a moiré image for Fourier transform is an image in which a zeroth-order component and a first-order component are arranged side by side as indicated by B4 in FIG.
- the moire image for the fringe scanning method when the orientation of the multi slit 12, the first grating 14, and the second grating 15 is 45 ° obliquely, the moire image for the fringe scanning method (slightly relative to the first grating 14).
- the stripes of the image obtained by slightly tilting the second grating are inclined 45 ° (oblique in the direction opposite to the slit direction) as indicated by C2 in FIG.
- the stripes of the moiré image for Fourier transform (the image obtained by further tilting the second grating) are finer diagonal stripes in the same direction as C2 in FIG. 24, as indicated by C3 in FIG.
- the image obtained by performing Fourier transform on the moiré image for Fourier transform is an image in which the zeroth-order component and the first-order component are arranged at an angle of 45 ° opposite to the fringe direction, as indicated by C4 in FIG.
- the zero-order component is cut out by the Hanning window W shown in FIG. 25 (step S36).
- the peripheral part of the Hanning window W is dropped to 0, and the central part of the Hanning window W is passed as it is.
- the cutting window function is not limited to the Hanning window, and a Hamming window, a Gaussian window, or the like may be used depending on the application.
- each of the extracted 0th order component and 1st order component is subjected to inverse Fourier transform (step S38).
- step S42 using the reconstructed image without the subject, the phase of the interference fringes is removed from the reconstructed image with the subject, and correction processing for removing image unevenness (artifact) is performed (step S42).
- the processing in step S42 is the same as that described in step S21 in FIG.
- the inventors of the present application have the information on the differential phase image and the small-angle scattered image of the Talbot interferometer and the Talbot-Lau interferometer as the grating (multi-slit 12, first grating).
- the window w used in steps S36 and S37 is not the conventional square but the slit direction of the grating as shown in FIG. If the rectangle extends in the orthogonal direction, the high-frequency component of the signal in the direction orthogonal to the slit direction of the grating including image information can be extracted without dropping, and blur in the direction orthogonal to the slit direction of the grating is reduced. It was found that this can be done (referred to as an improved Fourier transform method).
- the present invention is characterized in that the direction parallel to the slit direction of the lattice, in which a reduction in spatial resolution is unavoidable in principle, is not included in the image originally.
- a great advantage can be obtained for the imaging of the Fourier transform method using a dimensional lattice.
- FIG. 27 shows an example in which a rectangular window W is set in A4 of FIG. 22.
- FIG. 28A shows an example of a reconstructed image of a subject obtained by the fringe scanning method.
- FIG. 28B shows an example of a reconstructed image obtained by the improved Fourier transform method.
- FIG. 28C shows an example of a reconstructed image obtained by a conventional Fourier transform method.
- the reconstructed images in FIGS. 28A to 28C are differential phase images obtained by photographing with the slit direction of the grating being vertical.
- an image obtained by the fringe scanning method has little blur in both the vertical and horizontal directions.
- the image obtained by the improved Fourier transform method is blurred only in the vertical direction and is not blurred in the horizontal direction.
- FIG. 28A shows an example of a reconstructed image of a subject obtained by the fringe scanning method.
- FIG. 28B shows an example of a reconstructed image obtained by the improved Fourier transform method.
- FIG. 28C shows an example of a reconstructed image obtained by a conventional Fourier transform method.
- the image obtained by the conventional Fourier transform method is blurred both in the vertical direction and in the horizontal direction.
- differential phase images are shown, but the blur directions by the respective methods are the same in the absorption image and the small-angle scattering image.
- the improved Fourier transform method only the signal component in the direction parallel to the slit direction of the grating is blurred, so the first photographing is performed by arranging the subject so that the longitudinal direction of the subject is perpendicular to the slit direction of the grating.
- the relative angle between the subject and the grid is rotated by 90 °, the second image is taken, and reconstructed images are generated from the respective moire images obtained by the first and second times, By synthesizing the two reconstructed images generated, a two-dimensional image with little blur in both the vertical and horizontal directions of the subject can be acquired (in the case of a differential phase image and a small angle scattered image).
- the window w to be cut out is limited to a narrow range in both the vertical and horizontal directions. Therefore, it is inevitable that the spatial resolution is greatly reduced.
- the present method it is possible to generate a two-dimensional image using a one-dimensional lattice without greatly reducing the resolution.
- subject information at the four corners of the composite image is lost.
- the region of interest of the subject is generally located at the center of the imaging region.
- the above-mentioned subject information deficiency is less likely to be a problem. Further, since the photographing itself only needs to be performed twice, the influence of the subject's body movement can be suppressed.
- changing the shooting direction changing the slit direction relative to the subject
- the control unit 51 of the controller 5 determines the first shot image and the second shot image. After the reconstructed image creation / display processing by the Fourier transform method shown in FIG. 19 is executed for each, the two images are synthesized. When the same part of the subject is not drawn on the same pixel in the first image and the second image (when the subject is deformed or moved), either one of the images is moved in parallel or rotated. Then, after aligning at the position where the error between the two images becomes the smallest, the composition is performed. Various methods can be used as the synthesis method.
- the pixel of the first captured image is f1 (x, y)
- the pixel of the second captured image is f2 (x, y)
- the pixel of the composite image is g (x, y).
- Calculate and take the average power take the square root of the sum of squares).
- g (x, y) ⁇ (f1 (x, y) ⁇ 2 + f2 (x, y) ⁇ 2)
- the first captured image may be displayed in red
- the second captured image may be displayed in color.
- the reconstructed image obtained by the fringe scanning method is clearer and less blurred than the reconstructed image obtained by the Fourier transform method.
- the scanning for the fringe scanning method multiple images are taken continuously, so the imaging time becomes longer depending on the detector capture time, processing time before and after X-ray exposure, mechanism operation time, etc. 1 minute), body movement is likely to occur.
- the Fourier transform method one image is obtained by one image, so the image taking time depends only on the X-ray exposure time and can be suppressed to about 5 seconds. I can expect.
- the improved Fourier transform method can suppress the degradation of spatial resolution.
- the object that can fix the subject is the fringe scanning method, and the Fourier transform method is used to suppress body movement, and (2) The Fourier transform method is used for the simple inspection, which is a more precise inspection.
- a fringe scanning method it is possible to acquire an image according to the purpose and to perform imaging with less burden on the patient and less imaging.
- the relative angle between the first grating 14 and the second grating 15 can be easily adjusted to switch between the imaging for the fringe scanning method and the imaging for the Fourier transform method. Optimal shooting can be performed.
- step S43 the created reconstructed image is displayed on the display unit 53 (step S43).
- the display mode of the reconstructed image in step S22 of FIG. 14 and step S43 of FIG. 19 will be described.
- FIG. 29 shows an example of a display mode when the reconstructed image is displayed on the display unit 53 in step S22 of FIG. 14 and step S43 of FIG.
- steps S22 and S43 the three reconstructed images of the absorption image, the differential phase image, and the small angle scattered image are sequentially switched to the same position (region R0) of the display unit 53 at predetermined time intervals.
- Cycle display In addition, a stop button or a pause button may be provided on a screen outside the image area of the display unit 53 so that any image can be continuously displayed in a stationary state according to the operation of the operation unit 52. preferable.
- FIG. 18A to FIG. 18C the positions of the subjects in the three images are the same, and each represents information on different characteristics of the subject. For example, information on a large structural change appears in the absorption image.
- the differential phase image information on the phase change of the tissue edge appears.
- Information on scattering in the tissue appears in the small-angle scattering image. Therefore, as shown in FIG.
- the doctor who performs the interpretation does not need to move the line of sight and is tired. Can be interpreted while maintaining a high degree of concentration.
- the afterimage effect sin-called subliminal effect
- the afterimage effect when images are switched every predetermined time makes it possible to reconstruct a plurality of information (features) on the subject in his / her head, enabling high-accuracy diagnosis. Become.
- the three reconstructed images of the absorption image ⁇ the differential phase image ⁇ the small angle scattered image are sequentially switched in this order to be circulated, but the number of images to be displayed may be two or more.
- the order of switching display and the type of image are not particularly limited, and can be set for each part and for each user.
- the order of switching display is absorption image ⁇ small angle scattered image ⁇ differential phase image because a doctor can efficiently diagnose without taking extra time.
- a large structure change information appears in the absorption image. Therefore, by displaying the absorption image first, the doctor can grasp the distribution of fat and mammary glands and large lesions in the entire breast.
- a small-angle scattered image it is possible to detect calcification and an accumulation portion of cancer tissue without calcification with high sensitivity, and to determine the presence or absence of breast cancer.
- the differential phase image it is possible to detect the mass, the margin between the breast cancer tissue and the normal tissue, and identify the spread range of the breast cancer.
- a differential phase image is displayed to detect the mass and the marginal part of breast cancer tissue and normal tissue, to identify the extent of breast cancer spread, and then to display a small angle scattered image. Therefore, it is good to detect the accumulation part of the cancer tissue without calcification and calcification, and which display order should be decided for each doctor so that the abnormal shadow detection ability of the doctor finally becomes high. It may be determined.
- Absorption images have been conventionally used as diagnostic breast images, and are images that doctors are most familiar with (and are familiar with) diagnosis. Therefore, based on the diagnostic resolution that doctors have cultivated over many years of diagnosis 1 It is preferable to perform the next diagnosis.
- a small-angle scattered image and a differential phase image are displayed, and re-interpretation based on these images enables self-assessment of the results of primary diagnosis such as the presence or absence of abnormal shadows and benign / malignant abnormal shadows. Can be modified.
- abnormal shadows that were initially invisible and benign / malignant differences will become visible.
- the doctor will eventually establish a new step-up diagnostic resolution, so that even a diagnosis based only on the absorption image can make a diagnosis with higher accuracy than before. It is preferable.
- the subject part when the subject part is a limb, it is preferable to switch and display in the order of absorption image ⁇ differential phase image.
- the doctor can determine where the cartilage or tendon is located, and then displaying the differential phase image to determine the presence or absence of cartilage wear or tendon rupture. Is possible.
- an abnormal shadow candidate is detected from a medical image by an abnormal shadow candidate detection device (CAD: Computer-Aided Diagnosis), and the detection result is displayed together with the medical image for diagnosis. It is done.
- an absorption image is used as a medical image for diagnosis as described above, and detection of abnormal shadow candidates by CAD is conventionally performed only on the absorption image. Therefore, in the second embodiment, as in the conventional absorption image-based diagnosis system, among the reconstructed images, first, abnormal shadow candidates are detected with respect to the absorption image by interpretation and CAD, and the abnormal shadow candidates are detected.
- a small angle scattered image or a differential phase image is used in the secondary diagnosis for discriminating true positive / false positive of a detected region or the like will be described.
- the configuration of the X-ray imaging apparatus 1 and the controller 5 in the medical image display system according to the second embodiment and the operations from imaging to creation of a reconstructed image are the same as those described in the first embodiment, and therefore will be described. Incorporate.
- an absorption image of the reconstructed image is displayed on the display unit 53.
- the doctor interprets the image by observing the displayed absorption image.
- drawing of an image is omitted, but it is assumed that a breast image is drawn.
- an abnormal shadow candidate detection process is executed on the absorption image, and an abnormal shadow candidate is detected from the absorption image.
- the abnormal shadow candidate detection program is stored in the storage unit 55 of the controller 5, and the abnormal shadow candidate detection program stored in the control unit 51 and the storage unit 55 of the controller 5 As a result, the abnormal shadow candidate detection process is executed on the absorption image.
- a known algorithm can be applied as an abnormal shadow candidate detection algorithm.
- a method using an iris filter or a method using a Laplacian filter disclosed in Japanese Patent Application Laid-Open No. 10-91758 (Journal of the Institute of Electronics, Information and Communication Engineers (D-II) , Vol. J76-D-II, no. 2, pp. 241-249, 1993), etc. are applicable.
- micro-calcification cluster shadow candidates for example, Morphology filter (The Institute of Electronics, Information and Communication Engineers (D-II), Vol.J71-D-II, no.7, pp.1170-1176, 1992), Laplacian filter (The Institute of Electronics, Information and Communication Engineers Journal (D-II), Vol. J71-D-II, no. 10, pp. 1994-2001, 1998) Applicable.
- Morphology filter The Institute of Electronics, Information and Communication Engineers (D-II), Vol.J71-D-II, no.7, pp.1170-1176, 1992
- Laplacian filter The Institute of Electronics, Information and Communication Engineers Journal (D-II), Vol. J71-D-II, no. 10, pp. 1994-2001, 1998) Applicable.
- an annotation indicating the position where the abnormal shadow candidate is detected is displayed on the absorption image displayed on the display unit 53 as shown in FIG.
- the elliptical annotation indicates the position of the microcalcification cluster shadow candidate
- the rectangular annotation indicates the position of the tumor shadow candidate.
- the dotted annotation indicates a candidate suspected of being false positive.
- the doctor refers to the medical image and the annotation displayed on the display unit 53 and designates a candidate area as a target of the secondary diagnosis (for example, double-click) by the operation unit 52.
- a candidate region to be subjected to the secondary diagnosis is designated by the operation unit 52, display according to the designated type of abnormal shadow candidate (for example, a mass, a microcalcification cluster,...) Is performed.
- a tumor may appear in any of an absorption image, a small angle scattering image, and a differential phase image. Therefore, when a rectangular region, that is, a tumor candidate region is designated by the operation unit 52, as shown in FIGS. 30C to 30E, the small-angle scattered image ⁇ the differential phase image ⁇ the absorption image is switched every predetermined time. Are cycled.
- a small-angle scattered image is displayed on the display unit 53.
- areas other than the designated tumor candidate area are blackened (in a low luminance state converted to black). )Is displayed.
- the differential phase image is switched and displayed.
- FIG. 30D when the differential phase image is displayed, the area other than the designated area is also displayed after being blackened.
- E of FIG. 30 an absorption image in which the area other than the designated area is blackened is displayed. This is repeated.
- microcalcification clusters tend to appear in the differential phase image. Therefore, when an elliptical area, that is, a microcalcification cluster candidate area is designated by the operation unit 52, the differential phase image ⁇ the small angle scattered image ⁇ the absorption image is switched and displayed cyclically every predetermined time. At this time, in order to prevent dazzling due to unnecessary portions other than the designated area, areas other than the designated microcalcification cluster candidate areas are blackened and displayed.
- the absorption image-based abnormal shadow candidate detection result is true positive.
- False positive can be determined based on a small-angle scattered image and / or a differential phase image that reproduces features that do not appear in the absorption image of the subject.
- an abnormal shadow candidate is detected for an absorption image obtained by X-ray mammography first, and if there is an abnormality, for example, an ultrasonic diagnostic apparatus is used. Additional imaging was performed with other modalities to improve diagnostic accuracy, and finally connected to biopsy (biopsy test). As a result, the patient needed two visits and was burdensome. Ultrasound diagnosis may be performed in parallel in advance, but in this case, it was useless for a patient for whom an abnormal shadow candidate was not detected. Further, since the association between the abnormal shadow candidate area detected by CAD and the area to be imaged by the ultrasonic diagnostic apparatus or the like depends on the operation of the photographer himself, there is a risk of misdiagnosis due to an association error.
- a differential representing a feature different from an absorption image obtained by the same one imaging set as the absorption image used for CAD detection for an area where an abnormal shadow candidate is detected in CAD Since the phase image and the small-angle scattered image are provided, it is not necessary to perform another imaging such as ultrasonic diagnosis, and the burden on the patient can be reduced. Moreover, early diagnosis by the same doctor is possible. Further, since the subject and the arrangement thereof are the same as those of the absorption image, it is not necessary to align the images for the second diagnosis and is accurate, so that the diagnostic accuracy can be improved. Furthermore, since only the blackening process is performed without performing the alignment, the processing time can be shortened. By displaying the area other than the target area of the secondary diagnosis after performing the blackening process, it is possible to improve the diagnostic performance without requiring an interpretation time.
- a small-angle scattering image and a differential phase image are displayed, and re-reading based on these images is performed, so that the result of the primary diagnosis of the presence / absence of abnormal shadows and the benign / malignant abnormal shadows, etc. Can be self-correcting.
- the absorption image is read again after reading the small-angle scattered image and differential phase image, abnormal shadows that were not initially visible, benign / malignant differences, etc. will eventually become visible, and this diagnosis will be made.
- the doctor will eventually establish a new step-up diagnostic resolution so that even a diagnosis based only on the absorption image can make a diagnosis with higher accuracy than before. It is preferable.
- the operation unit 52 it is possible to set in advance by the operation unit 52 for each user and for each part, which timing of which kind of image to be switched and displayed is a small angle scattered image, a differential phase image, or an absorption image. Also, whether to perform blackening processing or display the entire image area may be set in advance.
- abnormal shadow candidates may be detected in the absorption image, and a small angle scattered image and a differential phase image may be created only for the detected region. Thereby, the processing time can be further shortened.
- the processing and display of a modified example described below may be performed.
- the abnormal shadow candidate detection process is executed on the absorption image in cooperation with the control unit 51 and the program stored in the storage unit 55, and the absorption image Detection of abnormal shadow candidates from.
- a reduced image 531a in which an annotation is added to the position of the detected abnormal shadow candidate on the absorption image is created.
- a main image 531b is created by combining the left breast image and the right breast image of the small angle scattered image or differential phase image with the chest wall. Then, as shown in FIG. 31, a diagnostic screen 531 in which the main image 531b is arranged at the center and the reduced image 531a is arranged outside the subject area of the main image 531b is displayed on the display unit 53.
- an abnormal shadow candidate is detected for an absorption image as in the prior art, and the result is obtained as a main life size consisting of a differential phase image or a small angle scattered image.
- the reduced image 531a is displayed in a place that does not interfere with the observation of the image 531b. Therefore, with regard to abnormal shadow candidates detected on the basis of the same absorption image as in the past, while confirming the position, a small-angle scattered image or differential phase image that reproduces features different from the absorption image is observed, and the abnormal shadow candidate is true. It becomes possible to make a secondary diagnosis of positive or false positive.
- an image representing characteristics different from the absorption image which is created by processing an image obtained by the same single imaging as the absorption image used for CAD detection. Therefore, the patient does not need to go to the hospital twice and the burden on the patient can be reduced. Early diagnosis by the same doctor is also possible. Further, since the position of the abnormal shadow candidate can be confirmed with the reduced image having the same arrangement of the main image and the subject, the position of the abnormal shadow candidate can be searched from the main image with high accuracy for diagnosis.
- both the small angle scattered image and the differential phase image may be switched and displayed as the main image 531b.
- the small-angle scattered image and the differential phase image having the same object arrangement with respect to the detector, which are created from the image obtained by one imaging set, are used, the first image of the left breast image and the right breast image are used. If the chest wall is aligned, the main image 531b in which the left and right breasts are aligned can be provided by simply switching the image data without performing alignment processing on the other side.
- the controller 51 and the abnormal shadow candidate detection program stored in the storage unit 55 cooperate to detect abnormalities from the absorption image, the small-angle scattered image, and the differential phase image.
- a shadow candidate is detected.
- the algorithm of the abnormal shadow candidate detection program applied to each image is common. After detection, as an image to be displayed first as an image from which an abnormal shadow candidate is detected, an absorption image, a small-angle scattered image, and a differential phase image are cyclically displayed on the display unit 53 at predetermined intervals as shown in FIG.
- the following modes are conceivable as the detection results of the abnormal shadow candidates of the absorption image, the small angle scattered image, and the differential phase image and the display order of each image.
- an abnormality detection candidate is detected from only one image
- images are displayed in the order of an image where an abnormal shadow candidate is detected ⁇ an absorption image ⁇ the remaining image.
- an abnormal shadow candidate is detected only from the absorption image
- the absorption image familiar to the doctor is displayed.
- the priorities of the small angle scattered image and the differential phase image may be set in advance from the operation unit 52.
- the abnormal shadow candidate detection program primarily detects a candidate area that may be an abnormal shadow candidate using a predetermined detection algorithm, and is calculated from the feature amount (or feature amount) of each detected primary candidate. Based on whether or not the index value) exceeds a predetermined threshold value, it is finally determined whether or not it is an abnormal shadow candidate (see, for example, Japanese Patent Application Laid-Open No. 2007-151465).
- TP apparent true positive zone
- FP apparent false positive zone
- the abnormal shadow candidates between the threshold values are gray zones in which the judgment of whether they are true positive or false positive differs depending on the sensitivity of the abnormal shadow candidate detection program. For a candidate in this gray zone, it is necessary to separately determine whether or not it is abnormal by a method other than CAD. Conventionally, the gray zone is determined by another modality. In the above (2a), when no abnormal shadow candidate is detected in all of the absorption image, the small-angle scattered image, and the differential phase image, the threshold value is lowered to approach the second threshold value, and the abnormal shadow in the gray zone closer to FP Allow candidates to be detected. Then, by first displaying an image from which an abnormal shadow candidate is detected, a doctor can concentrately interpret abnormal shadow candidates in a gray zone close to TP.
- the threshold value is increased to approach the first threshold value, and an abnormality in a gray zone that is closer to TP Enable to detect shadow candidates. Then, even if the threshold value is raised, an image in which abnormal shadow candidates are detected is displayed first, so that the doctor can concentrate on reading the abnormal shadow candidates in the gray zone close to the FP.
- two or more kinds of switching display of the absorption image, the small-angle scattering image, and the differential phase image may be performed independently as shown in FIG. 29, or as shown in FIGS. 33A, 33B, and 34.
- the small angle scattered image and the differential phase image may be displayed side by side on the screen of the display unit 53.
- FIG. 33A shows an example of a diagnostic screen 532 in which an absorption image 532a, a small-angle scattered image 532b, a differential phase image 532c, and a switching display image 532d are arranged.
- the switching display image 532d is an image obtained by switching and displaying two or more of the absorption image 532a, the small-angle scattered image 532b, and the differential phase image 532c every predetermined time.
- the diagnostic screen 532 shown in FIG. 33A when an ROI is set on any of the absorption image 532a, the small-angle scattered image 532b, and the differential phase image 532c by the operation unit 52, the ROI on the image is surrounded by a rectangle. The ROI is displayed in an identifiable manner on other reconstructed images as shown in FIG. 33A. Then, the area other than the ROI of the switching display image 532d is blackened. As shown in FIG. 33A, the diagnosis screen 532 is provided with a speed adjustment lever 532e. By operating the speed adjustment lever by the operation unit 52, the image displayed in the switching display image 532d is switched. The interval can be adjusted.
- the doctor examines the familiar absorption image 532 a and selects an area that seems to be suspicious by the operation unit 52 through an inquiry. Then, the corresponding region of the small-angle scattered image and the differential phase image is displayed in an identifiable manner as the ROI together with the selected region of the absorption image. Therefore, the wear and loss of the cartilage in that portion are diagnosed by the small-angle scattered image and the differential phase image. Or can be diagnosed for tendon or ligament tears. Then, it is possible to perform a diagnosis in which information of each image is integrated by a circular display focused on a place where the doctor wants to check. Further, according to the diagnostic screen 532, in addition to the switching display, three images representing different characteristics of the subject are displayed at the same time, so that the final confirmation can be performed with the most characteristic image without any operation. .
- an abnormal shadow candidate may be detected from the absorption image by CAD, and an annotation indicating the position of the abnormal shadow candidate may be displayed on the absorption image 532a as shown in FIG. 33B.
- an area detected as an abnormal shadow candidate by CAD may be automatically set as an ROI.
- the absorption image 533a that the doctor is used to making a diagnosis and the switching display image 533b may be displayed side by side on one screen.
- the switching display image 533b may display an absorption image at the initial display stage, and perform switching display for the region set as the ROI by the operation unit 52 in the absorption image 533a.
- the switching display images 532d and 533b may be switched among three images of absorption image ⁇ small angle scattered image (differential phase image) ⁇ differential phase image (small angle scattered image). It is good also as switching the image according to an abnormal shadow candidate detection result.
- the control unit 51 of the controller 5 in the X-ray imaging apparatus 1 uses the first imaging mode by the fringe scanning imaging apparatus or the second imaging system by the Fourier transform imaging apparatus. At least two of an X-ray absorption image, a differential phase image, and a small angle scattering image are created based on a moire image obtained by imaging in the imaging mode. Then, the display unit 53 of the created image is controlled.
- the control unit 51 sequentially switches and displays at least two of the generated X-ray absorption image, differential phase image, and small angle scattered image at the same position on the display unit 53.
- control is performed so as to sequentially switch and display an absorption image ⁇ a small-angle scattered image ⁇ a differential phase image.
- the absorption image information on a large structural change appears. Therefore, by displaying the absorption image first, the doctor can grasp the distribution of fat and mammary glands and large lesions in the entire breast.
- the X-ray absorption image, the differential phase image, and the small angle scattered image are sequentially switched and displayed at the same position on the display unit 53 in the order corresponding to the subject region, and thus created from the moire image.
- the reconstructed image can be used effectively, early diagnosis can be realized, and the diagnostic accuracy can be improved.
- abnormal shadow candidates are detected by CAD from the X-ray absorption image, differential phase image, and small angle scattered image, and switching is performed based on the detection result from each image.
- the type of image to be displayed and the display order it is possible to perform display according to the detection result of the abnormal shadow candidate. Specifically, by displaying in order from the image in which the abnormal shadow candidate is detected, it is possible for the doctor to focus on and interpret the portion detected as the abnormal shadow candidate, and to improve the diagnostic accuracy. Become.
- each image of the X-ray absorption image, the differential phase image, and the small angle scattered image on the same screen as the switching display, three images representing different characteristics of the subject are displayed. Since the images are displayed at the same time as the switching display, it is possible to perform the final confirmation of the diagnosis with the image in which the feature of the abnormal shadow candidate appears most after performing the diagnosis in which the features of each image are integrated by the switching display.
- the control unit 51 first displays an absorption image on the display unit 53, and then, from the absorption image as in the conventional case, CAD.
- the display unit 53 is controlled to perform switching display in the order of small-angle scattered image ⁇ differential phase image ⁇ absorption image.
- the display unit 53 is controlled to perform switching display in the order of differential phase image ⁇ small angle scattered image ⁇ absorption image.
- a differential phase image or a small angle scattered image is arranged as the main image in the center of the screen, and a reduced image in which the detection position of the abnormal shadow candidate is displayed on the absorption image is displayed outside the subject area of the main image.
- Control as follows. In this way, for a region where an abnormal shadow candidate is detected in CAD, a differential phase image representing a feature different from the absorption image obtained by the same imaging set as the absorption image used for CAD detection, By displaying the small-angle scattered image on the display unit 53, it is not necessary to perform another imaging such as ultrasonic diagnosis, and the burden on the patient can be reduced. Moreover, early diagnosis by the same doctor is possible.
- the subject and the arrangement thereof are the same as those of the absorption image, it is not necessary to align the images for the second diagnosis and is accurate, so that the diagnostic accuracy can be improved. Furthermore, since it is not necessary to align each image, the processing time can be shortened.
- the said embodiment is a suitable example of this invention, and is not limited to this.
- the X-ray imaging apparatus 1 has a multi-slit, and a plurality of moire images for the fringe scanning method are obtained by moving the multi-slit relative to the first grating and the second grating.
- the Talbot-Lau interferometer has been described as an example of generating a Talbot-Lau interferometer. However, the first grating and the second grating are moved relative to each other at regular intervals, and the X-ray source is irradiated every movement at regular intervals.
- a Talbot interferometer that generates a plurality of moire images for the fringe scanning method by repeating a process in which a radiation detector reads an image signal in accordance with X-rays may be used.
- the absorption image, the differential phase image, and the small angle scattered image obtained by reconstructing a plurality of moire images generated by the Talbot interferometer may be sequentially switched and displayed at the same position on the display unit 53 as described above.
- a reconstructed image based on one-dimensional image data photographed by an apparatus capable of both the fringe scanning method and the Fourier transform method (including the improved type) is used.
- the present invention is not limited to this. Instead, an apparatus dedicated to the Fourier transform method (including the improved type) may be used.
- the image was taken with a dedicated Fourier transform imaging device that uses the first and second gratings as a two-dimensional grating, and a dedicated Fourier transform imaging device that uses a multi-grating (two-dimensional grating) near the focal position. You may apply to the reconstructed image based on two-dimensional image data.
- the display method of the present invention has been described by taking the case where the absorption image, the differential phase image, and the small angle scattered image are sequentially switched and displayed at the same position on the display unit 53 as an example.
- the present invention is not limited to this, and can be applied to the case where a plurality of images created by performing different image processing on the same captured image are displayed.
- the saturation (color) of each display screen may be switched to, for example, black monotone, red monotone, or blue monotone so that it can be easily recognized that the image has been switched. .
- 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 arranged in this order (hereinafter referred to as the first arrangement).
- the arrangement of the X-ray source 11, the multi-slit 12, the first grating 14, the subject table 13, the second grating 15, and the X-ray detector 16 also includes the first grating 14 and A reconstructed image can be obtained by moving the multi slit 12 while the second grating 15 is fixed.
- the subject center and the first grid 14 are separated from each other by the thickness of the subject, which is slightly inferior in sensitivity compared to the above-described embodiment.
- the arrangement effectively uses X-rays by the amount of X-ray absorption in the first grating 14.
- the effective spatial resolution at the subject position depends on the X-ray focal spot diameter, the spatial resolution of the detector, the magnification of the subject, the thickness of the subject, and the like. When the resolution is 120 ⁇ m (Gauss half width) or less, the effective spatial resolution is smaller in the second arrangement than in the first arrangement. It is preferable to determine the order of arrangement of the first grating 14 and the object table 13 in consideration of sensitivity, spatial resolution, the amount of X-ray absorption in the first grating 14, and the like.
- the order of shooting with a subject and shooting without a subject is not limited to the above embodiment, and any order may be used. The same applies to the order of creating a reconstructed image with a subject and creating a reconstructed image without a subject.
- the medical field it may be used as a medical image display system that displays X-ray images.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Toxicology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
X線を照射するX線源と、
前記X線の照射軸方向と直交する方向に複数のスリットが配列されて構成された第1格子及び第2格子と、
被写体台と、
照射されたX線に応じて電気信号を生成する変換素子が2次元状に配置され、当該変換素子により生成された電気信号を画像信号として読み取るX線検出器と、
を備える縞走査型撮影装置又はフーリエ変換型撮影装置と、
前記何れかの撮影装置で撮影された被写体の画像信号に基づき、当該被写体のX線吸収画像、微分位相画像、及び小角散乱画像のうち少なくとも2つを生成する画像処理部と、
前記画像処理部で生成された画像を表示する表示部と、
前記画像処理部で生成された画像の前記表示部への表示を制御する制御部と、
を備える。
以下、図面を参照して本発明の第1の実施形態について説明する。
図1に、第1の実施形態に係る医用画像表示システムを示す。医用画像表示システムは、X線撮影装置1とコントローラ5を備える。X線撮影装置1は、縞走査型撮影装置として機能する第1の撮影モードと、フーリエ変換型撮影装置として機能する第2の撮影モードとを有する装置である。縞走査型撮影装置は、縞走査法による再構成画像用にタルボ・ロー干渉計により複数ステップで撮影を行い、複数のモアレ画像を生成するものである。フーリエ変換型撮影装置は、フーリエ変換法による再構成画像用に1又は2方向で撮影を行い、1又は2のモアレ画像を生成するものである。
本実施形態においては、X線撮影装置1の構成は、手指を被写体として撮影する装置を例にとり説明するが、これに限定されるものではない。
距離d2は、一般的に放射線科の撮影室の高さは3(m)程度又はそれ以下であることから、少なくとも3000(mm)以下であることが好ましい。なかでも、距離d2は400~2500(mm)が好ましく、さらに好ましくは500~2000(mm)である。
X線源11の焦点と第1格子14間の距離(d1+d3)は、好ましくは300~2800(mm)であり、さらに好ましくは400~1800(mm)である。
X線源11の焦点と第2格子15間の距離(d1+d3+d4)は、好ましくは400~3000(mm)であり、さらに好ましくは500~2000(mm)である。
それぞれの距離は、X線源11から照射されるX線の波長から、第2格子15上に第1格子14による格子像(自己像)が重なる最適な距離を算出し、設定すればよい。
X線源11は、緩衝部材17aを介して保持されている。緩衝部材17aは、衝撃や振動を吸収できる材料であれば何れの材料を用いてもよいが、例えばエラストマー等が挙げられる。X線源11はX線の照射によって発熱するため、X線源11側の緩衝部材17aは加えて断熱素材であることが好ましい。
X線の焦点径は、0.03~3(mm)が好ましく、さらに好ましくは0.1~1(mm)である。
マルチスリット12のスリット周期をw0(μm)、第1格子14のスリット周期をw1(μm)とすると、スリット周期w0は下記式により求めることができる。
w0=w1・(d3+d4)/d4
当該式を満たすように周期w0を決定することにより、マルチスリット12及び第1格子14の各スリットを通過したX線により形成される自己像が、それぞれ第2格子15上で重なり合い、いわばピントが合った状態とすることができる。
モータ部121aは、マイクロステップ駆動に切り替え可能なパルスモータであり、制御部181(図8参照)からの制御に応じて駆動され、ギア部121bを介してギア部121cをX線照射軸(図3に一点鎖線Rで示す)を中心として回転させる。ギア部121cは、ホルダー12bに保持されたマルチスリット12を装着するための開口部121eを有している。ギア部121cを回転させることにより、開口部121eに装着されたマルチスリット12をX線照射軸周りに回転させ、マルチスリット12のスリット配列方向を可変することができる。なお、撮影において、マルチスリット12は0°~90°程度回転できればよいので、ギア部121cは全周にある必要はなく、図3に2点鎖線で示す範囲(正逆回転方向にそれぞれ90°)で回転できればよい。
開口部121eは、ホルダー12bに保持されたマルチスリット12を上部から嵌め込むことが可能な形状及びサイズとなっている。ここでは、開口部121eにおけるスリット配列方向のサイズw4はホルダー12bにおけるスリット配列方向のサイズW2より若干大きくなっており、マルチスリット12をスリット配列方向にスライドさせることが可能となっている。なお、開口部121eにおけるスリット配列方向に直交する方向のサイズw3は、ホルダー12bにおけるスリット配列方向に直交する方向のサイズW1との精密嵌合可能な寸法としており、ホルダー12bを開口部121eに装着すると、ホルダー12bに設けられたラック12aは開口部121eの外に、後述するピニオン122cと係合可能に配置される。
モータ部122aは、例えば、制御部181からの制御に応じて駆動され、ギア部122bを介してピニオン122cを回転させる。ピニオン122cは、マルチスリット12のラック12aと係合して回転することで、マルチスリット12をスリット配列方向に移動させる。
被写体ホルダー130上には、更に被写体姿勢を安定させるため、指間スペーサ133を備えることが好ましい。また、患者毎に手や指間の大きさは異なるので、患者毎の手のひらの形状に合わせて被写体ホルダー130を作成しておき、撮影時には、その患者用の被写体ホルダー130を被写体台13にマグネット等で取り付けることが好ましい。腕から手首までの荷重は被写体台13が支えるので、被写体ホルダー130は指先部分の加重と患者が上方から押さえる力に耐えるものであればよく、安価で量産が可能なプラスチック成形とすることが可能である。
d4=(m+1/2)・w1 2/λ
なお、mは整数であり、λはX線の波長である。
X線源11のX線管の焦点径;300(μm)、管電圧:40(kVp)、付加フィルタ:アルミ1.6(mm)
X線源11の焦点からマルチスリット12までの距離d1 : 240(mm)
マルチスリット12から第1格子14までの距離d3 :1110(mm)
マルチスリット12から第2格子15までの距離d3+d4:1370(mm)
マルチスリット12のサイズ:10(mm四方)、スリット周期:22.8(μm)
第1格子14のサイズ:50(mm四方)、スリット周期:4.3(μm)
第2格子15のサイズ:50(mm四方)、スリット周期:5.3(μm)
ここで、本実施形態において、第1格子14は、複数のスリットが配列されてなる円形状の格子部140と、この格子部140を開口部212aに取り付けるための第1ホルダー部141及び第2ホルダー部142と、により構成されている(図6参照)。第1ホルダー部141は、格子部140の外周に取り付けられた、開口部212aと略同じ半径(外周の半径)の部材であり、第1格子14の装着時に開口部212aと嵌合する。第2ホルダー部142は、第1ホルダー部141より更に外側に取り付けられた、開口部212aよりやや半径(外周の半径)の大きい部材である。第2ホルダー部142は、その外周の一部がギア加工されている。また、第2ホルダー部142の外周の所定位置には、突起部142aが設けられている。
第2格子15は、複数のスリットが配列されてなる円形状の格子部150と、この格子部150を開口部212aに取り付けるためのホルダー部151と、により構成されている。ホルダー部151は、開口部212aの半径と略同様の半径を有する円盤状の部材である。格子部150はホルダー部151の中央部上面に保持されている(図6参照)。
ここで、X線撮影装置1は、縞走査法による再構成画像用に複数ステップで撮影を行う第1の撮影モードと、フーリエ変換法による再構成画像用に1又は2方向で撮影を行う第2の撮影モードとを有している。縞走査法用の撮影において要求される第1格子14のスリット方向と第2格子15のスリット方向との相対角度は第2格子の周期と画像サイズと縞の本数に依存する。縞走査法においては、モアレ画像における干渉縞の本数が少ないほど、また、干渉縞が鮮明なほど、このモアレ画像に基づいて作成される再構成画像が鮮明となることが知られている(非特許文献2参照)。そこで、第2格子の周期を5.3μm、60mm角の画像の中に干渉縞が0~3本程度にすると仮定すると、相対角度は0度~±0.015度とする必要がある。一方、フーリエ変換法用の撮影において要求される第1格子14のスリット方向と第2格子15のスリット方向との相対角度はX線検出器16の画素ピッチと空間分解能に依存する。一般的に使用される検出器(空間分解能30μm~200μm)とすると、その相対角度は0.4度~3度とする必要がある。よって、第1の撮影モードと第2の撮影モードを切り替えて撮影を行うためには、撮影モードに応じて第1格子と第2格子の相対角度を調整する必要がある。しかし、例えば、上記構成において、縞走査法の角度の0.005度のずれは、縞1周期分に相当する。縞走査法で常に縞が広がった状態を維持するためには、ミリ度の精度での調整が要求されるため、手動で第1格子14及び第2格子15のスリット方向の相対角度を調整するのは困難である。
そこで、X線撮影装置1においては、相対角調整部213により、操作部182で設定された撮影モードに応じて第1格子14と第2格子15の相対角度を自動的に調整することが可能となっている。
なお、このときの回転角度は1度程度と小さいため、まず、モータ部213aのパルスモータにより、突起部142aが基準位置215にくるまで反時計周りに第2ホルダー部142を回転させ、突起部142aが基準位置215に来たことが図示しないセンサにより検知されたら、第2ホルダー部142の回転方向を時計回りに切り替えてマイクロステップ駆動より第2ホルダー部142を回転させることが好ましい。
第2ホルダー部142は図示しないバネに付勢されており、モータ部213aの駆動により第1ギア213bと第2ホルダー部142との係合が解除されると、バネの付勢力により突起部142aがストッパー214の位置に戻る。即ち、第1格子14と第2格子15が第1の撮影モードに最適な相対角度に戻る。
格子回転部210は、また、相対角度が調整された第1格子14及び第2格子15を被写体に対してX線照射軸周り(図6に点線Rで示す)に一体的に回転させることができる。
ここで、1次元格子(スリット)を用いるタルボ干渉計及びタルボ・ロー干渉計では、第1格子14及び第2格子15のスリット方向と平行に線状に延びる構造物は鮮明に撮影することができないという特性がある。よって、被写体の注目すべき構造物の配置方向に応じて、第1格子14及び第2格子15のスリット方向の角度を調整する必要がある。格子回転部210は、以下の機構により第1格子14及び第2格子15をその相対角度を維持したまま一体的にX線照射軸周りに回転させ、被写体の注目すべき構造物の配置方向に対する第1格子14及び第2格子15のスリット方向の角度を調整することができる。
このように、回転トレイ212を手動で回転させるので、患者が触れる範囲に第1格子14及び第2格子15を一体的に回転させるための電気コード等を設ける必要がなく、安全性を確保することができる。
図7A、図7Bに示すように、保持部分171には、格子回転部210の回転トレイ212と精密嵌合するサイズであり、回転トレイ212を回転可能に保持する開口部171aと、回転トレイ212の回転角度を固定するためのトレイ固定部材171bと、が設けられている。開口部171aの底部とX線検出器16の載置部の間は、X線の透過を妨げないように、中空とするか又はX線透過率の高いアルミやカーボン等とすることが好ましい。トレイ固定部材171bは、凹部212b~212eの何れかがトレイ固定部材171bと対向するように位置したときにその対向する凹部に係合するボールと、ボールを図7A、図7Bの矢印方向に誘導するための図示しないスライドガイド(押圧バネのガイド)により構成されている。凹部212b~212eの何れかがトレイ固定部材171bと対向する位置で回転トレイ212の回転が停止すると、トレイ固定部材171bのスライドガイドにより、対向している凹部にボールが係合するとともに、凹部に設けられた角度検知センサ(SE1~SE4の何れか)によりボールの係合が検知されて制御部181に検知信号が出力される。これにより、制御部181は、回転トレイ212の回転角度、即ち、第1格子14及び第2格子15の回転角度を検知できるようになっている。
X線検出器16としては、FPD(Flat Panel Detector)を用いることができる。FPDにはX線をシンチレータを介して光電変換素子により電気信号に変換する間接変換型、X線を直接的に電気信号に変換する直接変換型があるが、何れを用いてもよい。
なお、CCD(Charge Coupled Device)、X線カメラ等の撮影手段をX線検出器16として用いてもよい。
まずFPDはリセットを行い、前回の撮影(読取)以降に残存する不要な電荷を取り除く。その後、X線の照射が開始するタイミングで電荷の蓄積が行われ、X線の照射が終了するタイミングで蓄積された電荷が画像信号として読み取られる。なお、リセットの直後や画像信号の読み取り後等に、オフセット補正用のダーク読み取りを行う。
制御部181は、CPU(Central Processing Unit)やRAM(Random Access Memory)等から構成され、記憶部185に記憶されているプログラムとの協働により、X線撮影装置1の各部を制御するとともに各種処理を実行する。例えば、制御部181は、後述する撮影制御処理をはじめとする各種処理を実行する。
表示部183は制御部181の表示制御に従って、ディスプレイに操作画面やX線撮影装置1の動作状況等を表示する。
記憶部185は、制御部181により実行されるプログラム、プログラムの実行に必要なデータを記憶している。また、記憶部185はX線検出器16によって得られたモアレ画像を記憶する。
制御部51は、CPU(Central Processing Unit)やRAM(Random Access Memory)等から構成され、記憶部55に記憶されているプログラムとの協働により、後述する縞走査法による再構成画像作成・表示処理、フーリエ変換法による再構成画像作成・表示処理をはじめとする各種処理を実行する。この縞走査法による再構成画像作成・表示処理、フーリエ変換法による再構成画像作成・表示処理を実行することにより、制御部51は、X線撮影装置1において第1の撮影モード、又は第2の撮影モードによる撮影によって得られたモアレ画像に基づいて、X線吸収画像、微分位相画像、及び小角散乱画像のうち少なくとも2つを作成し、作成された画像の表示部53への制御を行う。即ち、制御部51は、画像処理部、制御部として機能する。
また、記憶部55は、操作部52により設定された設定情報、例えば、図14のステップS22や図19のステップS43において表示する画像の種類や画像の表示の切替タイミング等を部位情報やユーザIDに対応付けて記憶する。
また、記憶部55は、X線検出器16によって得られたモアレ画像、モアレ画像に基づき作成された診断用の再構成画像を撮影オーダ情報に対応付けて記憶する。
また、記憶部55は、病変の典型的な症例を示す参照画像(詳細後述)を病変名及び画像の種類(縞走査法orフーリエ変換法、吸収画像or微分位相画像or小角散乱画像)等と対応付けて記憶する。
更に、記憶部55は、X線検出器16に対応するゲイン補正データ、欠陥画素マップ等を予め記憶する。
X線撮影装置1においては、通信部184によりコントローラ5から撮影条件の設定情報等が受信されると、X線撮影準備が実行される。
図10に示すように、X線源11から照射されたX線が第1格子14を透過すると、透過したX線がz方向に一定の間隔で像を結ぶ。この像を自己像といい、自己像が形成される現象をタルボ効果という。自己像を結ぶ位置に第2格子15が平行に配置され、当該第2格子15はその格子方向が第1格子14の格子方向と平行な位置からわずかに傾けられているので、第2格子15を透過したX線によりモアレ画像Mが得られる。X線源11と第1格子14間に被写体Hが存在すると、被写体HによってX線の位相がずれるため、図10に示すようにモアレ画像M上の干渉縞は被写体Hの辺縁を境界に乱れる。この干渉縞の乱れを、モアレ画像Mを処理することによって検出し、被写体像を画像化することができる。これがタルボ干渉計及びタルボ・ロー干渉計の原理である。
なお、第2の撮影モードによりモアレ画像を得る場合には、マルチスリット12の移動は行わず、1回撮影又は被写体とスリット方向を90度回転させて2回撮影を行う。
次いで、格子回転部210の相対角調整部213が制御され、第1格子14と第2格子15の相対角度が第1の撮影モードに最適となるように(突起部142aがストッパー214に接触する位置にくるように)第1格子14が回転される。これにより、第1格子14と第2格子15の相対角度が調整される(ステップS102)。
一般的には、マルチスリット12と第1格子14との相対角度が少ないほど、干渉縞の鮮明性の高いモアレ画像が得られる。しかし、マルチスリット12は発熱部であるX線源11近傍に配置されるので熱影響を受けやすい。そのため、マルチスリット12の変形等を考慮して、マルチスリット12を回転トレイ212と同じ角度だけ回転させるだけでなく、モータ部121aをマイクロステップ駆動させてステップS105~S108における微調整を行うことが有効である。
干渉縞の鮮明度=(MAX-MIN)/(MAX+MIN)=振幅/平均値
まず、マルチスリット12が停止した状態でX線源11によるX線の照射が開始される。X線検出器16ではリセット後、X線照射のタイミングに合わせて電荷が蓄積され、X線の照射停止のタイミングに合わせて蓄積された電荷が画像信号として読み取られる。これが1ステップ分の撮影である。1ステップ分の撮影が終了するタイミングで、制御部181の制御により駆動部122が起動され、マルチスリット12の移動が開始される。所定量移動すると駆動部122が停止されることによりマルチスリットの移動が停止され、次のステップの撮影が行われる。このようにして、マルチスリット12の移動と停止が所定のステップ数分だけ繰り返され、マルチスリット12が停止したときにX線の照射と画像信号の読み取りが行われる。読み取られた画像信号はモアレ画像として本体部18に出力される。
尚、ダーク画像の取得は、各ステップのモアレ画像取得後に、当該ステップのダーク読取を行って、各ステップ専用のオフセット補正データを生成することとしても良い。
尚、ダーク画像の取得は、各ステップのモアレ画像取得後に、当該ステップのダーク読取を行って、各ステップ専用のオフセット補正データを生成することとしても良い。
なお、被写体無しの複数のモアレ画像及びダーク読取は被写体有りの撮影直後に実施されることが精度的に最も好ましいが、被写体画像の再構成までの時間短縮をはかるために、始業開始等に事前に取得しておいたデータを用いることとしても良い。
コントローラ5においては、通信部54によりモアレ画像が受信されると、受信されたモアレ画像が撮影開始時に指定された撮影オーダ情報と対応付けて記憶部55に記憶される。
次いで、格子回転部210の相対角調整部213が制御され、第1格子14と第2格子15の相対角度が第2の撮影モードに最適となるように(突起部142aがホームポジションから所定角度回転した位置にくるように)調整される(ステップS202)。
撮影が終了すると、本体部18の通信部184からコントローラ5に撮影により得られたモアレ画像が送信される(ステップS211)。
撮影が終了すると、本体部18の通信部184からコントローラ5に、モアレ画像が送信される(ステップS216)。
なお、被写体無しの複数のモアレ画像及びダーク読取は被写体有りの撮影直後に実施されることが精度的に最も好ましいが、被写体画像の再構成までの時間短縮をはかるために、始業開始等に事前に取得しておいたデータを用いることとしても良い。
ステップS13においては、撮影に用いられたX線検出器16に対応する欠陥画素マップ(欠陥画素位置を示すデータ)が記憶部55から読み出され、各モアレ画像における欠陥画素位置マップで示す位置の画素値(信号値)が周辺画素により補間算出される。
具体的な処理としては、各モアレ画像の予め定められた1点の画素の信号値を用いて補正する方法、各モアレ画像間におけるX線検出器16の所定方向の信号値差を補正する(一次元補正する)方法、各モアレ画像間における2次元方向の信号値差を補正する(二次元補正する)方法、の何れであってもよい。
なお、上記サインカーブ形状は、マルチスリット12の開口幅、第1格子14及び第2格子15の周期、及び第1格子及び第2格子の格子間距離に依存し、また、放射光のようなコヒーレント光の場合には三角波形状となるが、マルチスリット効果によりX線が準コヒーレント光として作用する為、サインカーブを描くものとなる。ステップS15の解析は、被写体有りのモアレ画像と被写体無しのモアレ画像についてそれぞれ行われる。
X線撮影装置1では、当該制御情報に従って撮影のタイミングが調整され、再撮影が実行される。
具体的には、複数のモアレ画像の干渉縞を加算することにより吸収画像(X線吸収画像)が作成される(ステップS18)。また、縞走査法の原理を用いて干渉縞の位相が計算され、微分位相画像が作成される(ステップS19)。また、縞走査法の原理を用いて干渉縞のVisibilityが計算され(Visibility=2×振幅÷平均値)、小角散乱画像が作成される(ステップS20)。
例えば、被写体有りの再構成画像が微分位相画像である場合には、被写体有りの微分位相画像の各画素の信号値から被写体無しの微分位相画像の対応する(同じ位置の画素)の信号値を減算する処理が行われる(公知文献(A);Timm Weitkamp,Ana Diazand,Christian David, franz Pfeiffer and Marco Stampanoni, Peter Cloetens and Eric Ziegler, X-ray Phase Imaging with a grating interferometer,OPTICSEXPRESS,Vol.13, No.16,6296-6004(2005)、公知文献(B);Atsushi Momose, Wataru Yashiro, Yoshihiro Takeda, Yoshio Suzuki and Tadashi Hattori, Phase Tomography by X-ray Talbot Interferometry for Biological Imaging, Japanese Journal of Applied Physics, Vol.45, No.6A, 2006, pp.5254-5262(2006)参照)。
被写体有りの再構成画像が吸収画像、小角散乱画像である場合には、公知文献(C)に記載されているように、被写体有りの再構成画像の各画素の信号値を被写体無しの再構成画像の対応する画素の信号値で除算する割り算処理が行われる(公知文献(C);F.Pfeiffer, M.Bech,O.Bunk, P.Kraft, E.F.Eikenberry, CH.Broennimann,C.Grunzweig, and C.David,Hard-X-ray dark-field imaging using a grating interferometer, nature materials Vol.7,134-137(2008))。
図18Aに示すように、吸収画像は、被写体の大きな構造の変化を表すという特徴がある。図18Bに示すように、微分位相画像は、被写体の組織辺縁の位相変化を表すという特徴がある。図18Cに示すように、小角散乱画像は、被写体の組織内の散乱を表すという特徴がある。
図19は、制御部51により実行されるフーリエ変換法による再構成画像作成・表示処理を示すフローチャートである。フーリエ変換法による再構成画像作成・表示処理は、制御部51と記憶部55に記憶されているプログラムとの協働により実行される。
図20B、図21Bに示すように、1枚のモアレ画像をフーリエ変換すると、低周波成分(0次成分と呼ぶ)と干渉縞周波数付近の成分(1次成分と呼ぶ)、又は、0次成分と1次成分に加えさらに高周波成分(X線撮影装置1の干渉性に依存)が並んで得られる。0次成分と1次成分の並ぶ方向はモアレ画像の縞の方向に関係しており、モアレ画像の縞の方向に対して略直角になる。
例えば、図22のA1に示すように、マルチスリット12、第1格子14、第2格子15の格子の向きが縦であった場合、縞走査法用のモアレ画像(第1格子14に対して第2格子をわずかに傾けて得られた画像)の縞は、図22のA2に示すように横となる。フーリエ変換用のモアレ画像(第2格子を更に傾けて得られた画像)の縞は、図22のA3に示すように、図22のA2に比べて細かな横縞となる。フーリエ変換用のモアレ画像をフーリエ変換した画像は、図22のA4に示すように、0次成分、1次成分が縦に並んだ画像となる。
図23のB1に示すように、マルチスリット12、第1格子14、第2格子15の格子の向きが横であった場合、縞走査法用のモアレ画像(第1格子14に対して第2格子をわずかに傾けて得られた画像)の縞は、図23のB2に示すように縦となる。フーリエ変換用のモアレ画像(第2格子を更に傾けて得られた画像)の縞は、図23のB3に示すように、図23のB2に比べて細かな縦縞となる。フーリエ変換用のモアレ画像をフーリエ変換した画像は、図23のB4に示すように、0次成分、1次成分が横に並んだ画像となる。
図24のC1に示すように、マルチスリット12、第1格子14、第2格子15の向きが斜め45°であった場合、縞走査法用のモアレ画像(第1格子14に対してわずかに第2格子を少し傾けて得られた画像)の縞は、図24のC2に示すように斜め45°(スリット方向とは逆方向の斜め)となる。フーリエ変換用のモアレ画像(第2格子を更に傾けて得られた画像)の縞は、図24のC3に示すように、図24のC2と同方向の、より細かな斜め縞となる。フーリエ変換用のモアレ画像をフーリエ変換した画像は、図24のC4に示すように、0次成分、1次成分が縞方向とは逆の斜め45°に並んだ画像となる。
次いで、フーリエ変換により得られた画像において、1次成分が図26に示すようにキャリア周波数(=モアレ周波数)分シフトされ、Hanning窓Wで切り出される(ステップS37)。切り出しの窓関数はHanning窓に限定されず、用途に応じてHamming窓、ガウス窓等を使用しても良い。
次いで、切り出された0次成分、1次成分のそれぞれが逆フーリエ変換される(ステップS38)。
なお、図28A~図28Cにおいては、微分位相画像を示しているが、吸収画像、小角散乱画像においても各手法によるぼけの方向は同様である。
2次元格子を用いてフーリエ変換方式で撮影することも可能だが、縦横ともに1次成分が存在するため、切り出す窓wは縦横ともに狭い範囲に限定される。そのため空間分解能が大きく低下することを避けられない。一方、本方式によれば、1次元格子を使用して、分解能を大きく低下させずに2次元画像を生成することが可能となる。なお、2方向の撮影に基づく合成画像においては、合成画像の四隅部の被写体情報が欠損するが、医用現場の放射線撮影に於いては、一般的に被写体の関心領域を撮影領域の中央部に位置させることが多いので、上記の被写体情報欠損も問題なることが少ない。また、撮影自体も2回で済むため、被写体の体動影響も抑制可能である。
なお、撮影方向の変更(被写体に対するスリット方向の変更)する際には、縞走査法方式と同様に、マルチスリット12、第1格子14及び第2格子15を同時に90°回転させる必要がある。
g(x、y)=√(f1(x、y)^2+f2(x、y)^2)
また、例えば、1回目の撮影画像を赤、2回目の撮影画像を青にする等、カラー表示することとしてもよい。
以下、図14のステップS22及び図19のステップS43における再構成画像の表示態様について説明する。
よって、図29に示すように、3つの再構成画像を表示部53の同一位置(領域R0)に所定時間毎に順次切り替えて表示すれば、読影を行う医師は、視線の移動が不要で疲れを誘発しないので高い集中度合いを維持して読影を行うことができる。また所定時間毎に画像を切り替えた際の残像効果(所謂サブリミナル効果)により、被写体に関する複数の情報(特徴)を自身の頭の中で再構築可能となり、高い精度の診断を行うことが可能となる。
上述のように、吸収画像には、大きな構造の変化の情報が表れる。よって、吸収画像を最初に表示することで、医師は、***全体における脂肪や乳腺の分布や大きな病変を把握することができる。次いで、小角散乱画像を表示することで、石灰化と、石灰化を伴わない癌組織の集積部を感度高く検出し、乳癌の有無を判別することができる。次いで、微分位相画像を表示することで、腫瘤や乳癌組織と正常組織の辺縁部を検出し、乳癌の拡がり範囲を識別することができる。
なお、吸収画像に次いで、まず微分位相画像を表示することで、腫瘤や乳癌組織と正常組織の辺縁部を検出し、乳癌の拡がり範囲を識別し、然る後、小角散乱画像を表示することで、石灰化と、石灰化を伴わない癌組織の集積部を検出することとしても良く、どちらの表示順にするかは、最終的に医師の異常陰影検出能が高くなるように医師毎に定めることとしても良い。
また、吸収画像は、従来から診断用の***画像として用いられており、医師が最も診断に慣れた(見慣れた)画像であるため、各医師が長年の診断で培ってきた診断分解能に基づく1次診断を行うのに好ましい。吸収画像の表示後、小角散乱画像、微分位相画像を表示し、これらの画像に基づく再読影を行うことで、異常陰影の有無自体や異常陰影の良性/悪性等の1次診断の結果を自己修正可能となる。また、小角散乱画像、微分位相画像、の読影後に、再度、吸収画像の読影を行うと、最初は見えなかった異常陰影や、良性/悪性の差異等が、やがて、視認できるようになり、この診断方法を繰り返すことで、医師は、やがて、新たなステップアップした診断分解能を確立することになり、吸収画像のみに基づく診断であっても、従前よりも高精度の診断を行うことができるようになり、好ましい。
以下、本発明の第2の実施形態について説明する。
画像診断を支援するために、従来より、医用画像から異常陰影候補検出装置(CAD:Computer-Aided Diagnosis)による異常陰影候補の検出を行い、診断用の医用画像とともにその検出結果を表示することが行なわれている。従来、診断用の医用画像としては、上述のように吸収画像が用いられており、CADによる異常陰影候補の検出についても、従来、吸収画像に対してのみ行われている。
そこで、第2の実施形態においては、従来の吸収画像ベースの診断体系と同様に、再構成画像のうち、まず吸収画像に対して読影及びCADで異常陰影候補の検出を行い、異常陰影候補として検出された領域等の真陽性/偽陽性を判別する2次診断において小角散乱画像や微分位相画像を用いる例について説明する。
まず、図30のAに示すように、表示部53に、再構成画像のうち吸収画像が表示される。医師は、この表示された吸収画像を観察して読影を行う。なお、図30のAにおいては、画像の描画を省略しているが、***画像が描画されているものとする。
本実施形態の方式では、CADにおいて異常陰影候補の検出があった領域に対し、CADの検出に用いた吸収画像と同じ一回の撮影セットによって得られた、吸収画像とは異なる特徴を表す微分位相画像や小角散乱画像を提供するので、超音波診断等の再度の撮影が不要となり、患者の負担を低減することができる。また、同一医師による早期診断が可能である。また、被写体及びその配置が吸収画像と同一であるので、2回の診断の画像同士の位置合わせが不要で、かつ正確であり、診断精度の向上を図ることができる。更に、位置合わせを行わず、黒化処理を行うだけであるので、処理時間の短縮を図ることができる。2次診断の対象となる領域以外は黒化処理をして表示することで、いたずらに読影時間を必要とすることなく、診断性能を高めることができる。
異常陰影候補が検出された場合、吸収画像上の検出された異常陰影候補の位置にアノテーションを付した縮小画像531aが作成される。また、小角散乱画像又は微分位相画像の左***画像と右***画像を胸壁で合わせた主画像531bが作成される。そして、図31に示すように、主画像531bが中央に配置され、主画像531bの被写体領域外に縮小画像531aが配置された診断用画面531が表示部53に表示される。
この変形例においても、第2の実施形態と同様に、CADの検出に用いた吸収画像と同じ1回の撮影により得られた画像を処理して作成した、吸収画像とは異なる特徴を表す画像を提供するので、患者は2回の通院は不要で患者負担を低減することができる。同一医師による早期診断も可能となる。また、主画像と被写体の配置が同一の縮小画像で異常陰影候補の位置を確認できるので、異常陰影候補の位置を主画像から精度良く捜して診断を行うことが可能となる。
以下、本発明の第3の実施形態について説明する。
第3の実施形態においては、CADの検出結果を用いて第1の実施形態における切替表示の表示順を決定する例について説明する。
(1)1つの画像からのみ異常検出候補が検出された場合
この場合、異常陰影候補が検出された画像→吸収画像→残りの画像の順に表示する。吸収画像のみから異常陰影候補が検出された場合、まず、医師が見慣れた吸収画像から表示する。小角散乱画像と微分位相画像の優先順位は、操作部52から予め設定しておくこととしてもよい。
この場合、まず、医師が見慣れた吸収画像を表示し、次いで、操作部52からの切替指示に応じて、他の画像を表示する。切替指示に応じて表示する画像の種類は、操作部52から予め設定しておくこととする。或いは、医師の診断し易い表示順を予め設定しておくこととしてもよい。
(2a)3つの画像の全てから異常陰影候補が検出されなかった場合
異常陰影候補検出プログラムにおいて異常陰影候補の検出に用いられる閾値を下げて再度3つの画像のそれぞれに対して異常陰影候補検出処理を行うことにより、偽陽性に近いほうの候補を検出し、(1)と同様に、異常陰影候補が検出された画像から順に表示する。
(2b)3つの画像の全てから異常陰影候補が検出された場合
異常陰影候補検出プログラムにおいて異常陰影候補の検出に用いられる閾値を上げて再度3つの画像のそれぞれに対して異常陰影候補検出処理を行い、より真陽性に近いほうの候補まで検出し、(1)と同様に、異常陰影候補が検出された画像から順に表示する。
一般的に、異常陰影候補検出プログラムは、異常陰影候補の可能性のある候補領域を所定の検出アルゴリズムで一次検出し、検出された各一次候補についての特徴量(又は、特徴量から算出された指標値)が予め定められた閾値を超えたか否かに基づいて、最終的に異常陰影候補であるか否かの判定を行っている(例えば、特開2007-151465号公報参照)。図32に示す第1の閾値を超えている場合は明白な真陽性ゾーン(TP)、第2の閾値を下回る場合は明白な偽陽性ゾーン(FP)とすると、第1の閾値と第2の閾値の間にある異常陰影候補は、真陽性であるか偽陽性であるかの判断が異常陰影候補検出プログラムの感度によって異なるグレーゾーンとなる。このグレーゾーンにある候補については、CAD以外の方法で異常であるか否かを別途判定する必要がある。従来は、グレーゾーンについては他のモダリティで判別を行っていた。上述の(2a)では、吸収画像、小角散乱画像、微分位相画像の全てにおいて異常陰影候補が検出されなかった場合、閾値を下げて第2の閾値に近づけ、よりFPに近いグレーゾーンの異常陰影候補を検出できるようにする。そして、異常陰影候補が検出された画像を最初に表示することで、医師がTPに近いグレーゾーンの異常陰影候補を集中して読影できるようにする。また、上述の(2b)では、吸収画像、小角散乱画像、微分位相画像の全てにおいて異常陰影候補が検出された場合、閾値を上げて第1の閾値に近づけ、よりTPに近いグレーゾーンの異常陰影候補を検出できるようにする。そして、閾値を上げても異常陰影候補が検出された画像を最初に表示することで、医師がFPに近いグレーゾーンの異常陰影候補を集中して読影できるようにする。
検出された数が多いほうの画像と吸収画像を交互に切替表示、同時2画面表示、又は、重ねて1画面表示する。
(3a)吸収画像から異常陰影候補が検出されず、微分位相画像と小角散乱画像から異常陰影候補が検出された場合
微分位相画像用、小角散乱画像用の異常陰影候補検出プログラムの閾値を上げて再度異常陰影候補検出処理を行い、より真陽性に近い候補のみを検出し、(1)と同様に、異常陰影候補が検出された画像から順に表示する。
具体的には、被写体部位が***画像である場合には、吸収画像→小角散乱画像→微分位相画像の順次切替表示するよう制御する。吸収画像には、大きな構造の変化の情報が表れる。よって、吸収画像を最初に表示することで、医師は、***全体における脂肪や乳腺の分布や大きな病変を把握することができる。次いで、小角散乱画像を表示することで、石灰化と、石灰化を伴わない癌組織の集積部を感度高く検出し、乳癌の有無を判別することができる。次いで、微分位相画像を表示することで、腫瘤や乳癌組織と正常組織の辺縁部を検出し、乳癌の拡がり範囲を識別することができる。
また、被写体部位が手足等である場合には、吸収画像→微分位相画像の順に切替表示するよう制御する。まず吸収画像を表示することで、医師は、軟骨や腱のある場所の見当をつけることができ、次いで、微分位相画像を表示することで、軟骨の摩耗や腱の破裂などの有無を判別することができる。
このように、CADにおいて異常陰影候補の検出があった領域に対し、CADの検出に用いた吸収画像と同じ一回の撮影セットによって得られた、吸収画像とは異なる特徴を表す微分位相画像や小角散乱画像を表示部53に表示することで、超音波診断等の再度の撮影が不要となり、患者の負担を低減することができる。また、同一医師による早期診断が可能である。また、被写体及びその配置が吸収画像と同一であるので、2回の診断の画像同士の位置合わせが不要で、かつ正確であり、診断精度の向上を図ることができる。更に、各画像の位置合わせを行う必要がないので、処理時間の短縮を図ることができる。
例えば、上記実施形態においては、X線撮影装置1を、マルチスリットを有し、マルチスリットを第1格子及び第2格子に対して相対的に移動させることで縞走査法用の複数のモアレ画像を生成するタルボ・ロー干渉計の構成とした例を説明したが、第1格子と第2格子とを一定周期間隔で相対移動させ、一定周期間隔での移動毎にX線源により照射されたX線に応じて放射線検出器が画像信号を読み取る処理を繰り返すことで縞走査法用の複数のモアレ画像を生成するタルボ干渉計としてもよい。そして、タルボ干渉計で生成された複数のモアレ画像を再構成することにより得られる吸収画像、微分位相画像、小角散乱画像を上述のように表示部53の同一位置に順次切り替え表示することとしてもよい。
また、上記実施例では、縞走査方式とフーリエ変換方式(改良型含む)の両方式が可能な装置により撮影された1次元画像データに基づく再構成画像を使用したが、これに限定されるものでは無く、フーリエ変換方式(改良型含む)専用の装置を用いても良い。
また、第1格子及び第2格子を2次元格子としたフーリエ変換方式専用の撮影装置や、更に焦点位置近傍にマルチ格子(2次元格子)を併用したフーリエ変換方式専用の撮影装置で撮影された2次元画像データに基づく再構成画像に適用しても良い。
また、切替表示に際しては、画像が切替られたことを容易に視認可能なように、各表示画面の彩度(色)を、例えば、黒モノトーン、赤モノトーン、青モノトーンのように切替えても良い。
第2の配置においては、被写体の厚み分だけ、被写体中心と第1格子14は離れることになり、上記の実施形態に比べ感度の点でやや劣ることになるが、一方で、被写体への被曝線量低減を考慮すると、当該配置の方が第1格子14でのX線吸収分だけX線を有効に活用していることになる。
また、被写***置での実効的な空間分解能は、X線の焦点径、検出器の空間分解能、被写体の拡大率、被写体の厚さ等に依存するが、上記実施例に於ける検出器の空間分解能が120μm(ガウスの半値幅)以下の場合には、第1の配置よりも第2の配置の方が実効的な空間分解能は小さくなる。
感度、空間分解能、及び、第1格子14でのX線吸収量等を考慮して、第1格子14、被写体台13の配置順をきめることが好ましい。
11 X線源
12 マルチスリット
12a ラック
12b ホルダー
121 マルチスリット回転部
121a モータ部
121b ギア部
121c ギア部
121d 支持部
121e 開口部
122 駆動部
122a モータ部
122b ギア部
122c ピニオン
13 被写体台
130 被写体ホルダー
131 楕円形状
133 指間スペーサ
14 第1格子
140 格子部
141 第1ホルダー部
142 第2ホルダー部
142a 突起部
15 第2格子
150 格子部
151 ホルダー部
16 X線検出器
17 保持部
17a 緩衝部材
171a 開口部
171b トレイ固定部材
18 本体部
181 制御部
182 操作部
183 表示部
184 通信部
185 記憶部
18a 駆動部
210 格子回転部
211 ハンドル
212 回転トレイ
212a 開口部
212b~212e 凹部
213 相対角調整部
213a モータ部
213b 第1ギア
213c 第2ギア
213d レバー
214 ストッパー
5 コントローラ
51 制御部
52 操作部
53 表示部
54 通信部
55 記憶部
Claims (3)
- X線を照射するX線源と、
前記X線の照射軸方向と直交する方向に複数のスリットが配列されて構成された第1格子及び第2格子と、
被写体台と、
照射されたX線に応じて電気信号を生成する変換素子が2次元状に配置され、当該変換素子により生成された電気信号を画像信号として読み取るX線検出器と、
を備える縞走査型撮影装置又はフーリエ変換型撮影装置と、
前記何れかの撮影装置で撮影された被写体の画像信号に基づき、当該被写体のX線吸収画像、微分位相画像、及び小角散乱画像のうち少なくとも2つを生成する画像処理部と、
前記画像処理部で生成された画像を表示する表示部と、
前記画像処理部で生成された画像の前記表示部への表示を制御する制御部と、
を備える医用画像表示システム。 - 前記制御部は、前記画像処理部で生成された少なくとも2つの画像を前記表示部に所定時間毎に順次切替表示させる請求項1に記載の医用画像表示システム。
- 前記縞走査型撮影装置は、前記X線源近傍に配置されたマルチスリットを有し、当該マルチスリットを前記第1格子及び第2格子に対して相対移動させるタルボ・ロー干渉計である請求項1又は2に記載の医用画像表示システム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013506010A JP5915645B2 (ja) | 2011-03-23 | 2012-03-22 | 医用画像表示システム |
US14/006,838 US20140010344A1 (en) | 2011-03-23 | 2012-03-22 | Medical image display system |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011063470 | 2011-03-23 | ||
JP2011-063470 | 2011-03-23 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2012128335A1 true WO2012128335A1 (ja) | 2012-09-27 |
Family
ID=46879475
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2012/057404 WO2012128335A1 (ja) | 2011-03-23 | 2012-03-22 | 医用画像表示システム |
Country Status (3)
Country | Link |
---|---|
US (1) | US20140010344A1 (ja) |
JP (1) | JP5915645B2 (ja) |
WO (1) | WO2012128335A1 (ja) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015015851A1 (ja) * | 2013-07-30 | 2015-02-05 | コニカミノルタ株式会社 | 医用画像システム及び関節軟骨状態のスコア判定方法 |
JP2015024068A (ja) * | 2013-07-29 | 2015-02-05 | コニカミノルタ株式会社 | 医用画像処理装置 |
JP2015507984A (ja) * | 2012-02-24 | 2015-03-16 | パウル・シェラー・インスティトゥート | 人体組織の微小石灰化の複数の異なるタイプの非侵襲的分類システム |
JP2017524139A (ja) * | 2014-08-13 | 2017-08-24 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | トモグラフィにおける暗視野イメージング |
CN108236473A (zh) * | 2016-12-26 | 2018-07-03 | 株式会社岛津制作所 | X射线相位摄影装置 |
JP2019005047A (ja) * | 2017-06-22 | 2019-01-17 | 株式会社島津製作所 | X線イメージング装置およびx線イメージング画像の合成方法 |
JP2019045471A (ja) * | 2017-09-01 | 2019-03-22 | 株式会社島津製作所 | X線イメージング装置 |
WO2020054158A1 (ja) * | 2018-09-11 | 2020-03-19 | 株式会社島津製作所 | X線位相撮像システム |
WO2020066135A1 (ja) * | 2018-09-25 | 2020-04-02 | 株式会社島津製作所 | X線位相イメージング装置 |
WO2020188970A1 (ja) * | 2019-03-15 | 2020-09-24 | コニカミノルタ株式会社 | 表示制御装置、放射線撮影システム及びプログラム |
US10837922B2 (en) | 2017-09-01 | 2020-11-17 | Shimadzu Corporation | X-ray imaging apparatus |
WO2023171726A1 (ja) * | 2022-03-10 | 2023-09-14 | コニカミノルタ株式会社 | 評価方法、評価装置及びプログラム |
WO2024070776A1 (ja) * | 2022-09-30 | 2024-04-04 | コニカミノルタ株式会社 | 画像表示装置、画像表示システム、画像表示方法、画像処理装置及びプログラム |
WO2024070777A1 (ja) * | 2022-09-30 | 2024-04-04 | コニカミノルタ株式会社 | 画像表示装置、画像表示システム、画像表示方法、画像処理装置及びプログラム |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008102685A1 (ja) * | 2007-02-21 | 2008-08-28 | Konica Minolta Medical & Graphic, Inc. | 放射線画像撮影装置及び放射線画像撮影システム |
AU2012268876A1 (en) * | 2012-12-24 | 2014-07-10 | Canon Kabushiki Kaisha | Non-linear solution for 2D phase shifting |
JP6369206B2 (ja) * | 2014-08-06 | 2018-08-08 | コニカミノルタ株式会社 | X線撮影システム及び画像処理装置 |
JP2016050891A (ja) * | 2014-09-01 | 2016-04-11 | キヤノン株式会社 | X線撮像装置 |
JP2016106721A (ja) * | 2014-12-03 | 2016-06-20 | キヤノン株式会社 | 画像処理装置および画像処理方法 |
JP6613988B2 (ja) * | 2016-03-30 | 2019-12-04 | コニカミノルタ株式会社 | 放射線撮影システム |
JP2017205138A (ja) * | 2016-05-16 | 2017-11-24 | コニカミノルタ株式会社 | 放射線撮影装置 |
JP6750310B2 (ja) * | 2016-05-30 | 2020-09-02 | コニカミノルタ株式会社 | タルボ撮影装置 |
EP3391821B1 (en) * | 2017-04-20 | 2024-05-08 | Shimadzu Corporation | X-ray phase contrast imaging system |
KR102081869B1 (ko) * | 2018-01-08 | 2020-02-26 | 연세대학교 원주산학협력단 | 엑스선 영상 생성 장치 및 방법 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010050483A1 (ja) * | 2008-10-29 | 2010-05-06 | キヤノン株式会社 | X線撮像装置およびx線撮像方法 |
JP2010253157A (ja) * | 2009-04-28 | 2010-11-11 | Konica Minolta Medical & Graphic Inc | X線干渉計撮影装置及びx線干渉計撮影方法 |
JP2011045655A (ja) * | 2009-08-28 | 2011-03-10 | Konica Minolta Medical & Graphic Inc | X線撮影装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4218347B2 (ja) * | 2003-01-17 | 2009-02-04 | コニカミノルタホールディングス株式会社 | 画像診断支援装置 |
JP4820680B2 (ja) * | 2006-04-12 | 2011-11-24 | 株式会社東芝 | 医用画像表示装置 |
US20120014578A1 (en) * | 2010-07-19 | 2012-01-19 | Qview Medical, Inc. | Computer Aided Detection Of Abnormalities In Volumetric Breast Ultrasound Scans And User Interface |
JP2011224329A (ja) * | 2010-03-30 | 2011-11-10 | Fujifilm Corp | 放射線撮影システム及び方法 |
-
2012
- 2012-03-22 US US14/006,838 patent/US20140010344A1/en not_active Abandoned
- 2012-03-22 JP JP2013506010A patent/JP5915645B2/ja active Active
- 2012-03-22 WO PCT/JP2012/057404 patent/WO2012128335A1/ja active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010050483A1 (ja) * | 2008-10-29 | 2010-05-06 | キヤノン株式会社 | X線撮像装置およびx線撮像方法 |
JP2010253157A (ja) * | 2009-04-28 | 2010-11-11 | Konica Minolta Medical & Graphic Inc | X線干渉計撮影装置及びx線干渉計撮影方法 |
JP2011045655A (ja) * | 2009-08-28 | 2011-03-10 | Konica Minolta Medical & Graphic Inc | X線撮影装置 |
Non-Patent Citations (1)
Title |
---|
W. YASHIRO: "On the origin of visibility contrast in x-ray Talbot interferometry", OPTICS EXPRESS, vol. 18, no. 16, 2 August 2010 (2010-08-02), pages 16890 - 16901 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015507984A (ja) * | 2012-02-24 | 2015-03-16 | パウル・シェラー・インスティトゥート | 人体組織の微小石灰化の複数の異なるタイプの非侵襲的分類システム |
JP2015024068A (ja) * | 2013-07-29 | 2015-02-05 | コニカミノルタ株式会社 | 医用画像処理装置 |
US10085701B2 (en) | 2013-07-30 | 2018-10-02 | Konica Minolta, Inc. | Medical image system and joint cartilage state score determination method |
WO2015015851A1 (ja) * | 2013-07-30 | 2015-02-05 | コニカミノルタ株式会社 | 医用画像システム及び関節軟骨状態のスコア判定方法 |
JP2017524139A (ja) * | 2014-08-13 | 2017-08-24 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | トモグラフィにおける暗視野イメージング |
CN108236473A (zh) * | 2016-12-26 | 2018-07-03 | 株式会社岛津制作所 | X射线相位摄影装置 |
JP2019005047A (ja) * | 2017-06-22 | 2019-01-17 | 株式会社島津製作所 | X線イメージング装置およびx線イメージング画像の合成方法 |
US10837922B2 (en) | 2017-09-01 | 2020-11-17 | Shimadzu Corporation | X-ray imaging apparatus |
JP2019045471A (ja) * | 2017-09-01 | 2019-03-22 | 株式会社島津製作所 | X線イメージング装置 |
JP7110697B2 (ja) | 2017-09-01 | 2022-08-02 | 株式会社島津製作所 | X線イメージング装置 |
WO2020054158A1 (ja) * | 2018-09-11 | 2020-03-19 | 株式会社島津製作所 | X線位相撮像システム |
JPWO2020054158A1 (ja) * | 2018-09-11 | 2021-09-24 | 株式会社島津製作所 | X線位相撮像システム |
JP7040625B2 (ja) | 2018-09-11 | 2022-03-23 | 株式会社島津製作所 | X線位相撮像システム |
WO2020066135A1 (ja) * | 2018-09-25 | 2020-04-02 | 株式会社島津製作所 | X線位相イメージング装置 |
WO2020188970A1 (ja) * | 2019-03-15 | 2020-09-24 | コニカミノルタ株式会社 | 表示制御装置、放射線撮影システム及びプログラム |
WO2023171726A1 (ja) * | 2022-03-10 | 2023-09-14 | コニカミノルタ株式会社 | 評価方法、評価装置及びプログラム |
WO2024070776A1 (ja) * | 2022-09-30 | 2024-04-04 | コニカミノルタ株式会社 | 画像表示装置、画像表示システム、画像表示方法、画像処理装置及びプログラム |
WO2024070777A1 (ja) * | 2022-09-30 | 2024-04-04 | コニカミノルタ株式会社 | 画像表示装置、画像表示システム、画像表示方法、画像処理装置及びプログラム |
Also Published As
Publication number | Publication date |
---|---|
US20140010344A1 (en) | 2014-01-09 |
JPWO2012128335A1 (ja) | 2014-07-24 |
JP5915645B2 (ja) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5915645B2 (ja) | 医用画像表示システム | |
JP5935693B2 (ja) | 医用画像表示方法 | |
JP5652245B2 (ja) | X線撮影システム | |
JP5900324B2 (ja) | X線撮影システム | |
JP5331940B2 (ja) | 放射線撮影システム及び放射線画像生成方法 | |
JP5857800B2 (ja) | 関節撮影装置及び撮影対象固定ユニット | |
JP5708652B2 (ja) | X線撮影システム | |
JP5238786B2 (ja) | 放射線撮影装置及び放射線撮影システム | |
JP2012090944A (ja) | 放射線撮影システム及び放射線撮影方法 | |
JP2012115576A (ja) | 放射線画像検出装置、放射線撮影装置、放射線撮影システム | |
JP2012148068A (ja) | 放射線画像取得方法および放射線画像撮影装置 | |
JP2012200567A (ja) | 放射線撮影システム及び放射線撮影方法 | |
JP2012095865A (ja) | 放射線撮影装置、放射線撮影システム | |
JP2012143550A (ja) | 放射線画像撮影装置および放射線画像取得方法 | |
JP5831614B2 (ja) | X線撮影システム | |
JP2013085631A (ja) | 関節撮影装置 | |
WO2012169426A1 (ja) | 放射線撮影システム | |
JP6237001B2 (ja) | 医用画像処理装置及び位相画像生成方法 | |
JP6413685B2 (ja) | X線撮影システム | |
JP6451400B2 (ja) | 画像処理システム及び画像処理装置 | |
WO2012056992A1 (ja) | 放射線画像検出装置、放射線撮影装置、放射線撮影システム | |
WO2012147749A1 (ja) | 放射線撮影システム及び放射線撮影方法 | |
WO2012057046A1 (ja) | 放射線撮影装置及び放射線撮影システム | |
JP2012135561A (ja) | X線撮影装置 | |
JP2011206280A (ja) | 放射線撮影方法及びシステム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 12759925 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2013506010 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14006838 Country of ref document: US |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 12759925 Country of ref document: EP Kind code of ref document: A1 |