CN111566471B - X射线成像装置 - Google Patents

X射线成像装置 Download PDF

Info

Publication number
CN111566471B
CN111566471B CN201880082441.7A CN201880082441A CN111566471B CN 111566471 B CN111566471 B CN 111566471B CN 201880082441 A CN201880082441 A CN 201880082441A CN 111566471 B CN111566471 B CN 111566471B
Authority
CN
China
Prior art keywords
grating
dark field
fiber bundle
image
ray
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201880082441.7A
Other languages
English (en)
Other versions
CN111566471A (zh
Inventor
佐野哲
木村健士
白井太郎
土岐贵弘
堀场日明
森本直树
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shimadzu Corp
Original Assignee
Shimadzu Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shimadzu Corp filed Critical Shimadzu Corp
Publication of CN111566471A publication Critical patent/CN111566471A/zh
Application granted granted Critical
Publication of CN111566471B publication Critical patent/CN111566471B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. 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 for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • A61B6/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/40—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4035—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis the source being combined with a filter or grating
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/46—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A61B6/466—Displaying means of special interest adapted to display 3D data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. 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 for radiation diagnosis, e.g. 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/041—Phase-contrast imaging, e.g. using grating interferometers
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/42—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4291—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating

Abstract

该X射线成像装置(100)具备：旋转机构(8)，其使由X射线源(1)、检测器(5)、第一光栅(2)以及第二光栅(3)构成的摄像系统(200)相对地旋转；以及图像处理部(6)，其在多个旋转角度的各个旋转角度下基于X射线的强度分布来生成暗场像。而且，图像处理部(6)构成为实施使暗场像的多个像素中的、暗场信号比阈值(V1)大的像素的暗场信号降低到设定值(V2)的散射校正。

Description

X射线成像装置
技术领域
本发明涉及一种X射线成像装置，特别是涉及一种具备生成被摄体的暗场像的图像处理部的X射线成像装置。
背景技术
以往，已知一种具备生成被摄体的暗场像的图像处理部的X射线成像装置。这种X射线成像装置例如在M Bech,et al.,“Quantitative x-ray dark-field computedtomography”PHYSICS IN MEDICINE AND BIOLOGY 55(2010)P.5529-P.5539中被公开。
上述M Bech,et al.,“Quantitative x-ray dark-field computed tomography”PHYSICS IN MEDICINE AND BIOLOGY 55(2010)P.5529-P.5539的X射线成像装置具备相位光栅、检测器以及图像处理装置。来自X射线源的X射线由于被摄体而发生散射，通过相位光栅后被照射到检测器。由图像处理装置基于由检测器检测到的X射线的干涉强度来生成暗场像。另外，使用在被摄体旋转了360度的情况下的各个角度下生成的暗场像来进行图像的重构。
在如上述M Bech,et al.,“Quantitative x-ray dark-field computedtomography”PHYSICS IN MEDICINE AND BIOLOGY 55(2010)P.5529-P.5539所记载那样的以往的X射线成像装置中，在X射线的散射十分微小的情况下，能够近似为X射线的干涉强度相对于散射体的透射长度进行指数函数式的衰减。在该情况下，能够通过普通的FBP(Filtered Back Projection：滤波背投影)来适当地进行重构。
现有技术文献
专利文献
非专利文献1：M Bech,et al.,“Quantitative x-ray dark-field computedtomography”PHYSICS IN MEDICINE AND BIOLOGY 55(2010)P.5529-P.5539
发明内容
发明要解决的问题
然而，在如上述M Bech,et al.,“Quantitative x-ray dark-field computedtomography”PHYSICS IN MEDICINE AND BIOLOGY 55(2010)P.5529-P.5539所记载那样的以往的X射线成像装置中，在X射线的散射角大等情况下，由于X射线的干涉强度并非相对于透射长度进行指数函数式的变化，因此无法适当地进行基于FBP的重构。即，存在以下问题：在通过FBP重构出的暗场像(三维数据)中产生伪影(噪声)。在该情况下，存在重构出的被摄体的像变得不清晰这样的问题。
本发明是为了解决如上所述的问题而完成的，本发明的一个目的在于提供一种能够抑制重构出的被摄体的像变得不清晰的X射线成像装置。
用于解决问题的方案
为了实现上述目的，本发明的一个方面的X射线成像装置具备：X射线源；检测器，其检测从X射线源照射的X射线；多个光栅，所述多个光栅设置在X射线源与检测器之间，包括第一光栅和第二光栅；旋转机构，其使被摄体与摄像系统相对地旋转，摄像系统由X射线源、检测器以及多个光栅构成；以及图像处理部，其基于在通过旋转机构进行了旋转的情况下的多个旋转角度的各个旋转角度下由检测器检测出的X射线的强度分布，来生成由X射线的散射引起的暗场像，其中，图像处理部构成为实施使以下像素的基于X射线的散射强度的值降低到规定的设定值的散射校正：暗场像的多个像素中的、基于X射线的散射强度的值比规定的阈值大的像素。此外，“暗场像”是指根据基于物体的小角散射的Visibility的变化而获得的Visibility像。另外，暗场像也被称为小角散射像。“Visibility”是指清晰度。
在本发明的一个方面的X射线成像装置中，如上所述那样构成为使基于X射线的散射强度的值比规定的阈值大的像素的基于X射线的散射强度的值降低到规定的设定值。由此，在散射校正前，即使在基于X射线的散射强度的值并非相对于与旋转角度对应地变化的(X射线的)透射长度进行指数函数式的变化的情况下，也能够通过优化散射校正中的规定的设定值并使基于X射线的散射强度的值降低到适当的值来使基于X射线的散射强度的值相对于旋转角度(透射长度)的变化接近于指数函数式的变化。其结果，能够抑制由X射线的散射角大引起的伪影(噪声)的产生，并且能够抑制重构出的被摄体的像变得不清晰。
在上述一个方面的X射线成像装置中，优选的是，图像处理部构成为：在散射校正中，将基于X射线的散射强度的值比规定的阈值大的像素的基于X射线的散射强度的值设定为规定的阈值。如果像这样构成，则与将规定的设定值和规定的阈值设为不同的值的情况相比，能够抑制用户决定规定的设定值所耗费的作业负担的增加。
在上述一个方面的X射线成像装置中，优选的是，被摄体包括纤维束，第一光栅具有第一光栅构成部分，第二光栅具有沿着第一光栅构成部分所延伸的方向延伸的第二光栅构成部分，图像处理部构成为：在旋转机构的旋转轴延伸的方向与第一光栅构成部分及第二光栅构成部分各自延伸的方向一致的情况下，对包括纤维束的被摄体的暗场像进行散射校正。在此，一般情况下，在暗场像中显著地显现沿着第一光栅构成部分和第二光栅构成部分各自延伸的方向延伸的纤维束。但是，在旋转机构的旋转轴延伸的方向与第一光栅构成部分及第二光栅构成部分各自延伸的方向一致的情况下，除了显现沿着第一光栅构成部分及第二光栅构成部分各自延伸的方向延伸的纤维束以外，还显现沿着与第一光栅构成部分及第二光栅构成部分各自延伸的方向正交的方向延伸的纤维束。这是因为，由于通过旋转机构使被摄体旋转，因此根据旋转角度不同，存在X射线对于沿着与第一光栅构成部分及第二光栅构成部分各自延伸的方向正交的方向延伸的纤维束的透射长度变长的情况(角度)。因此，由沿着与第一光栅构成部分及第二光栅构成部分各自延伸的方向正交的方向延伸的纤维束引起的X射线的散射角变大。因而，在旋转机构的旋转轴延伸的方向与第一光栅构成部分及第二光栅构成部分各自延伸的方向一致的情况下，通过对包括纤维束的被摄体的暗场像进行散射校正，能够更有效地抑制由沿着与第一光栅构成部分及第二光栅构成部分各自延伸的方向正交的方向延伸的纤维束导致的X射线的散射所引起的伪影(噪声)的产生。其结果，能够更有效地抑制重构出的被摄体的像变得不清晰。
在该情况下，优选的是，纤维束包括沿规定方向延伸的第一纤维束以及沿着与规定方向不同的方向延伸的第二纤维束，图像处理部构成为：通过散射校正，使基于由第一纤维束引起的X射线的散射强度的值和基于由第二纤维束引起的X射线的散射强度的值的合计值降低到规定的设定值。如果像这样构成，则与分别独立地通过散射校正使基于由第一纤维束引起的X射线的散射强度的值和基于由第二纤维束引起的X射线的散射强度的值降低的情况不同，能够通过一次的散射校正来控制基于X射线的散射强度的值。其结果，能够抑制散射校正所需的时间变长。
在上述使基于由第一纤维束引起的X射线的散射强度的值和基于由第二纤维束引起的X射线的散射强度的值的合计值降低到规定的设定值的X射线成像装置中，优选的是，第一纤维束和第二纤维束中的一方被设置为沿着第一光栅构成部分和第二光栅构成部分各自延伸的方向延伸，第一纤维束和第二纤维束中的另一方被设置为沿着从正面观察第一光栅和第二光栅时与第一光栅构成部分及第二光栅构成部分各自延伸的方向正交的方向延伸。如果像这样构成，则通过散射校正，能够抑制由第一纤维束和第二纤维束中的沿着与第一光栅构成部分及第二光栅构成部分各自延伸的方向正交的方向延伸的一方导致的X射线的散射所引起的伪影(噪声)的产生。其结果，能够抑制第一纤维束和第二纤维束中的沿着第一光栅构成部分和第二光栅构成部分各自延伸的方向延伸的另一方的像变得不清晰。
在上述在旋转轴延伸的方向与第一光栅构成部分及第二光栅构成部分各自延伸的方向一致的情况下进行散射校正的X射线成像装置中，优选的是，图像处理部构成为：在旋转机构的旋转轴延伸的方向与第一光栅构成部分及第二光栅构成部分各自延伸的方向一致的情况下，对多个旋转角度的各个旋转角度下的暗场像进行散射校正，并且对多个旋转角度的各个旋转角度下的被进行了散射校正的暗场像进行重构，由此生成第一立体暗场数据。如果像这样构成，则能够在多个旋转角度的各个旋转角度下使暗场信号(基于X射线的散射强度的值)的变化接近指数函数式的衰减的变化，因此能够抑制第一立体暗场数据变得不清晰。
在该情况下，优选的是，图像处理部构成为能够受理是否实施散射校正的变更，并且图像处理部构成为：在从正面观察第一光栅和第二光栅时旋转轴延伸的方向与第一光栅构成部分及第二光栅构成部分各自延伸的方向大致正交的情况下，对多个旋转角度的各个旋转角度下的暗场像不进行散射校正而进行重构，由此生成第二立体暗场数据，并且将第一立体暗场数据和第二立体暗场数据进行组合来获取被摄体的立体像。在此，一般情况下，在暗场像中，在从正面观察第一光栅和第二光栅时旋转轴延伸的方向与第一光栅构成部分及第二光栅构成部分各自延伸的方向大致正交的情况下，在暗场像中显著地显现沿着第一光栅构成部分和第二光栅构成部分各自延伸的方向延伸的纤维束。这是由于，在通过旋转机构使被摄体旋转的情况下，对于沿着与第一光栅构成部分及第二光栅构成部分各自延伸的方向正交的方向(沿着旋转轴的方向)延伸的纤维束来说，无论在哪个旋转角度，都具有针对纤维束的暗场灵敏度，并且X射线的透射长度不会根据旋转角度而大幅地变化。由此，在第二立体暗场数据中显著地显示的纤维束和在第一立体暗场数据中显著地显示的纤维束成为互不相同(沿互不相同的方向延伸)的纤维束。因而，通过将第一立体暗场数据和第二立体暗场数据进行组合，用户能够在立体像中视觉识别在第一立体暗场数据中显著地显示的纤维束与在第二立体暗场数据中显著地显示的纤维束的位置关系。
在上述一个方面的X射线成像装置中，优选的是，图像处理部构成为能够受理规定的设定值的变更。如果像这样构成，则能够根据被摄体的种类来变更规定的设定值，因此能够更适当地实施散射校正。
发明的效果
根据本发明，如上所述，能够抑制重构出的被摄体的像变得不清晰。
附图说明
图1是示出一个实施方式的X射线成像装置的结构的示意图。
图2的(a)是示出旋转轴与光栅方向正交的情况下的X射线成像装置的结构的结构图。图2的(b)是图2的(a)的情况下的立体暗场数据。
图3的(a)是示出旋转轴与光栅方向一致的情况下的X射线成像装置的结构的结构图。图3的(b)是图3的(a)的情况下(无散射校正)的立体暗场数据(比较例)。
图4是用于说明一个实施方式的X射线成像装置的、被摄体相对于X射线的光轴的配置角度以及X射线的散射程度的示意图(a)～(d)。
图5是一边通过一个实施方式的X射线成像装置使被摄体旋转一边进行摄像时的暗场像的示意图(a)～(d)。
图6是由一个实施方式的X射线成像装置生成的暗场像的正弦图。(图6的(a)是图2的(a)的情况下的图。图6的(b)是图3的(a)的情况下(无散射校正)的图。)
图7是用于说明一个实施方式的X射线成像装置的由第一纤维束引起的暗场信号与由第二纤维束引起的暗场信号的关系的图。(图7的(a)是图2的(a)的情况下的图。图7的(b)是图3的(a)的情况下(无散射校正)的图。)
图8是对由图3的(a)的X射线成像装置生成的暗场像进行散射校正并进行重构而生成的立体暗场数据。
图9是由一个实施方式的X射线成像装置得到的立体暗场数据的截面像。(图9的(a)是图3的(b)的截面像(比较例)。图9的(b)是图8的截面像。)
图10是由一个实施方式的X射线成像装置得到的(使用与图9不同的被摄体生成的)立体暗场数据的截面像。(图10的(a)是图3的(b)的截面像(比较例)。图10的(b)是图8的截面像。)
图11是由一个实施方式的X射线成像装置得到的立体暗场数据的(沿着与图10不同的方向的方向的)截面像。(图11的(a)是图3的(b)的截面像(比较例)。图11的(b)是图8的截面像。)
图12是示出由一个实施方式的X射线成像装置得到的立体暗场数据的图。(图12的(a)是图3的(b)示出的立体暗场数据(比较例)。图12的(b)是图8示出的立体暗场数据。)
图13是示出将由一个实施方式的X射线成像装置得到的在图8中示出的立体暗场数据与在图2的(b)中示出的立体暗场数据进行组合而获取到的立体像的图。
图14是用于说明一个实施方式的X射线成像装置的摄像方法的流程图。
具体实施方式
以下，基于附图来说明将本发明具体化的实施方式。
[本实施方式]
参照图1～图14对本实施方式的X射线成像装置100的结构进行说明。
(X射线成像装置的结构)
如图1所示，X射线成像装置100是利用通过了被摄体T的X射线的扩散(散射)来将被摄体T的内部图像化的装置。另外，X射线成像装置100是利用塔尔博特(Talbot)效应将被摄体T的内部图像化的装置。X射线成像装置100例如在无损检查用途中能够用于物体的内部的图像化。
被摄体T在内部包括纤维束10(参照图2)。被摄体T例如是将碳纤维用作纤维束10、将树脂11(参照图2)用作母材的碳纤维强化塑料(CFRP)。此外，所谓纤维束是将多个纤维集合成束状而得到的。在本实施方式中，纤维束10是利用多个纤维形成为板状而得到的。
图1是从X方向观察X射线成像装置100的图。如图1所示，X射线成像装置100具备X射线源1、第一光栅2、第二光栅3、第三光栅4、检测器5、图像处理部6、控制部7、旋转机构8以及光栅移动机构9。此外，在本说明书中，将从X射线源1朝向第一光栅2的方向设为Z2方向，将其逆向的方向设为Z1方向。另外，将与Z方向正交的面内的左右方向设为X方向，将朝向纸面的内部的方向设为X2方向，将朝向纸面的近前侧的方向设为X1方向。另外，将与Z方向正交的面内的上下方向设为Y方向，将向上方向设为Y1方向，将向下方向设为Y2方向。
X射线源1构成为通过被施加高电压来产生X射线，并且朝向Z2方向照射所产生的X射线。
第一光栅2具有沿Y方向以规定的周期(间距)d1排列的多个狭缝2a和X射线相位变化部2b。各狭缝2a和X射线相位变化部2b分别形成为直线状地延伸。另外，各狭缝2a和X射线相位变化部2b分别形成为平行地延伸。第一光栅2是所谓的相位光栅。此外，X射线相位变化部2b是本发明的“第一光栅构成部分”的一例。
第一光栅2被配置在X射线源1与第二光栅3之间，被X射线源1照射X射线。第一光栅2是为了通过塔尔博特效应形成第一光栅2的自身像(未图示)而设置的。此外，当具有相干性的X射线通过形成有狭缝的光栅时，在与光栅相距规定的距离(塔尔博特距离)的位置处形成光栅的像(自身像)。将该情况称为塔尔博特效应。
第二光栅3具有沿Y方向以规定的周期(间距)d2排列的多个X射线透过部3a和X射线吸收部3b。X射线吸收部3b沿着X射线相位变化部2b所延伸的方向延伸。各X射线透过部3a和X射线吸收部3b分别形成为直线状地延伸。另外，各X射线透过部3a和X射线吸收部3b分别形成为平行地延伸。第二光栅3是所谓的吸收光栅。第一光栅2、第二光栅3分别是具有不同的作用的光栅，但狭缝2a和X射线透过部3a分别使X射线透过。另外，X射线吸收部3b发挥屏蔽X射线的作用，X射线相位变化部2b由于折射率与狭缝2a的折射率不同而使X射线的相位变化。此外，X射线吸收部3b是本发明的“第二光栅构成部分”的一例。
第二光栅3被配置在第一光栅2与检测器5之间，被照射通过了第一光栅2的X射线。另外，第二光栅3被配置在与第一光栅2相距塔尔博特距离的位置。第二光栅3与第一光栅2的自身像发生干涉，从而在检测器5的检测表面上形成莫尔条纹(未图示)。
第三光栅4具有沿Y方向以规定的周期(间距)d3排列的多个X射线透过部4a和X射线吸收部4b。各X射线透过部4a和X射线吸收部4b分别形成为直线状地延伸。另外，各X射线透过部4a和X射线吸收部4b分别形成为平行地延伸。第三光栅4是所谓的多狭缝。
第三光栅4被配置在X射线源1与第一光栅2之间。第三光栅4构成为：将通过了各X射线透过部4a的X射线作为线光源，由此使来自X射线源1的X射线多点光源化。通过使三个光栅(第一光栅2、第二光栅3以及第三光栅4)的间距以及光栅间的距离满足固定的条件，能够提高从X射线源1照射的X射线的相干性。由此，即使X射线源1的管球的焦点尺寸大，也能够保持干涉强度。
检测器5构成为检测X射线，并且将检测到的X射线转换为电信号，读取转换得到的电信号来作为图像信号。检测器5例如是FPD(Flat Panel Detector：平板检测器)。检测器5由多个转换元件(未图示)和配置在多个转换元件上的像素电极(未图示)构成。多个转换元件和像素电极以规定的周期(像素间距)沿X方向和Y方向阵列状地排列。另外，检测器5构成为将获取到的图像信号输出到图像处理部6。
图像处理部6构成为基于从检测器5输出的图像信号来生成吸收像(未图示)。另外，图像处理部6构成为基于由检测器5检测出的X射线的强度分布来生成暗场像(参照图5)。在此，吸收像是将由于由被摄体T引起的X射线的吸收差异而产生的对比度进行图像化而得到的。另外，暗场像是将由于由处于被摄体T的内部的细微构造引起的X射线的折射而产生的对比度进行图像化而得到的。另外，图像处理部6对一边使旋转机构8旋转一边(在多个旋转角度的各个旋转角度下)进行拍摄所得到的多个吸收像和多个暗场像分别进行重构，由此生成三维吸收像(未图示)和三维暗场像(参照图2的(b)等)。另外，图像处理部6例如包括GPU(Graphics Processing Unit：图形处理单元)、构成为用于图像处理的FPGA(Field-Programmable Gate Array：现场可编程门阵列)等处理器。
控制部7构成为借助旋转机构8使被摄体T与由X射线源1、检测器5及多个光栅(第一光栅2、第二光栅3及第三光栅4)构成的摄像系统200相对地旋转。另外，控制部7构成为借助光栅移动机构9使第一光栅2在光栅面内沿着与光栅方向正交的方向进行步进移动。在X射线成像装置100中，使用了根据通过以固定周期间隔扫描第一光栅2所得到的多个莫尔条纹(图像)来获取重构图像的方法(条纹扫描法)。另外，控制部7例如包括CPU(CentralProcessing Unit：中央处理单元)等处理器。此外，光栅方向是本发明的“第一光栅构成部分和第二光栅构成部分各自延伸的方向”的一例。
旋转机构8构成为基于来自控制部7的信号使被摄体T与摄像系统200相对地旋转。具体地说，旋转机构8构成为通过使被摄体T绕着轴线AR进行旋转来使被摄体T相对于摄像系统200相对地旋转。在图1中，图示了轴线AR延伸的方向(在图1中为Y方向)与多个光栅的光栅方向(在图1中为X方向)正交的状态，但多个光栅的光栅方向能够变化为规定方向(例如Y方向)。此外，光栅方向是指光栅的光栅图案延伸的方向。另外，光栅图案是指各光栅的狭缝2a、X射线相位变化部2b、X射线透过部3a以及X射线吸收部3b等。另外，旋转机构8例如包括由电动机等驱动的旋转台。此外，轴线AR是本发明的“旋转轴”的一例。
光栅移动机构9构成为基于来自控制部7的信号使第一光栅2在光栅面内(XY面内)沿着与光栅方向正交的方向(在图1中为Y方向)进行步进移动。具体地说，光栅移动机构9将第一光栅2的周期d1分割为n份，使第一光栅2每次步进移动d1/n。光栅移动机构9构成为使第一光栅2至少步进移动与第一光栅2的一个周期d1相当的量。此外，n是正的整数，例如是9等。另外，光栅移动机构9例如包括步进电动机、压电致动器等。
如图2的(a)所示，纤维束10在被摄体T内包括沿不同的多个方向延伸的纤维束10。在图2的(a)所示的例子中，被摄体T为编织了纤维束10a和纤维束10b而得到的构造，其中，该纤维束10a沿着从正面观察(从Z方向观察)第一光栅2和第二光栅3时与光栅方向正交的方向(Y方向)延伸，该纤维束10b沿着光栅方向(X方向)延伸。此外，在图2的(a)所示的例子中，区分地显示了纤维束10a和纤维束10b，但实际上是织成(或者，编织为)片状的同种的碳纤维。另外，X方向和Y方向分别是本发明的“与规定方向不同的方向”和“规定方向”的一例。另外，纤维束10a和纤维束10b分别是本发明的“第一纤维束”和“第二纤维束”的一例。
另外，有时在被摄体T的树脂11中形成有裂纹12。在图2所示的例子中，裂纹12沿Y方向延伸。此外，裂纹12是在树脂11中形成的空气层等。
在图2所示的例子中，轴线AR延伸的方向(在图2的(a)中为Y方向)与光栅方向(在图2的(a)中为X方向)正交。在该情况下，在暗场像中显现纤维束10中的沿光栅方向延伸的纤维束(图2的(a)中为纤维束10b)的像。通过对在该情况下生成的多个旋转角度下的暗场像(参照图5的暗场像20a～20d)进行重构，来生成清晰地显示纤维束10b的立体像的3D数据20(参照图2的(b))。此外，3D数据20是本发明的“第二立体暗场像”的一例。
在图3所示的比较例中，轴线AR延伸的方向(在图3的(a)中为Y方向)与光栅方向(在图3的(a)中为Y方向)一致。在该情况下，在暗场像中除了显现纤维束10中的沿光栅方向延伸的纤维束(在图3的(a)中为纤维束10a)和裂纹12的像以外，还显现沿着与光栅方向正交的方向(X方向)延伸的纤维束(在图3的(a)中为纤维束10b)。通过对在该情况下生成的多个旋转角度下的暗场像(参照图5的暗场像21b～21d)进行重构，来生成3D数据21(参照图3的(b))，该3D数据21除了显现纤维束10a和裂纹12的立体像以外，还显现纤维束10b的立体像。
在此，参照图4，对在轴线AR延伸的方向与光栅方向一致的情况下(参照图3的(a))在暗场像(参照图5的暗场像21b～21d)中显现沿着与光栅方向正交的方向(X方向)延伸的纤维束(在图3的(a)中为纤维束10b)的理由进行说明。
图4的(a)是表示在XZ平面内配置为被摄体T的纤维束10b与X射线的光轴正交的情况下的例子的示意图。图4的(b)是沿图4的(a)中的400-400线的截面图。图4的(c)是表示在XZ平面内配置为被摄体T的纤维束10b相对于X射线的光轴倾斜的情况下的例子的示意图。图4的(d)是沿图4的(c)中的500-500线的截面图。
在此，在向沿着与X射线的光轴正交的方向延伸的纤维束10入射了X射线的情况下，X射线因纤维束10而散射(扩散)。具体地说，在被摄体T内的纤维与树脂11(参照图2)的界面，由于纤维与树脂11的折射率的差导致X射线发生折射。在被摄体T内，层叠有纤维束10，另外纤维束10由多个纤维构成，因此由于通过多个纤维而发生多重的折射，X射线发生散射(扩散)。
在图4的(a)所示的例子中为纤维束10b与X射线的光轴正交的配置。即，纤维束10b被配置为沿X方向延伸。因而，如图4的(b)所示，纤维的截面大致为圆形。在图4的(c)所示的例子中，纤维束10b相对于X射线的光轴倾斜地配置。因而，如图4的(d)所示，纤维的截面为椭圆形。在纤维的截面为椭圆形的情况下，与纤维的截面大致为圆形的情况相比，X射线通过纤维与树脂11的界面的距离变长。另外，纤维与树脂11的界面为粗糙形状(微小的凹凸)，X射线的散射(扩散)增大与X射线通过纤维与树脂11的界面的距离变长的部分相当的量。
图5是示出在本实施方式的X射线成像装置100中一边使被摄体T旋转一边进行摄像时的某角度(第一角度θ1、第二角度θ2、第三角度θ3以及第四角度θ4)下的、轴线AR延伸的方向与光栅方向正交的情况(参照图2)下的暗场像20a～20d以及轴线AR延伸的方向与光栅方向一致的情况(参照图3)下的暗场像21a～21d的示意图(a)～(d)。
在轴线AR延伸的方向与光栅方向正交的情况下，在第一角度θ1的暗场像(暗场像20a)(参照图5的(a))中清晰地显现出纤维束10b。另外，在角度比较大的暗场像(暗场像20b～20d)(参照图5的(b)～(d))中，随着角度变大，由纤维束10b引起的X射线的散射(扩散)(参照图4)变大。其结果，纤维束10b的暗场信号(参照图7的(a)的实线)变大。此外，在该情况下，即使角度变大，由纤维束10a引起的X射线的散射(扩散)也不会大幅地变化，因此纤维束10a(以及裂纹12)的暗场信号(参照图7的(a)的虚线)不会变大。在此，暗场信号是在将暗场像中的规定像素的像素值设为V的情况下通过f(1/V)求出的值。此外，上述函数f(x)是自变量x的值越大则值越大的函数。因而，在X射线的散射强度大而暗场像的像素值(明亮度)小的情况下(即，在V大的情况下)，暗场信号变大。此外，暗场信号是本发明的“基于X射线的散射强度的值”的一例。
另一方面，在轴线AR延伸的方向与光栅方向一致的情况下，在角度较小的第一角度θ1的暗场像(暗场像21a)中，纤维束10b的暗场信号(参照图7的(b)的实线)比较小，清晰地显现出纤维束10a(以及裂纹12)。但是，随着角度变大，由纤维束10b引起的X射线的散射(扩散)变大，因此纤维束10b的暗场信号(参照图7的(b)的实线)相对于纤维束10a的暗场信号(参照图7的(b)的虚线)变大。因此，在暗场像21b～21d中显著地显现纤维束10b的像，因此难以掌握纤维束10a(以及裂纹12)的详细的构造。
图6是在本实施方式的X射线成像装置100中一边使被摄体T旋转360度一边进行摄像时的正弦图。图6的正弦图示出了多个旋转角度的各个旋转角度下的检测器5的输出数据。如轴线AR延伸的方向与光栅方向正交的情况下的正弦图(参照图6的(a))以及轴线AR延伸的方向与光栅方向一致的情况下的正弦图(参照图6的(b))分别示出的那样，可知在规定的旋转角度下，由纤维束10b引起的X射线的散射(扩散)变大(参照图6的(a)和图6的(b)的斜线部分)。即，图6的(a)和图6的(b)的斜线部分对应于暗场信号比较大的旋转角度(参照图7的(a)和图7的(b)的实线)。
在此，在本实施方式中，图像处理部6(参照图1)构成为实施使暗场像的多个像素中的、暗场信号比阈值V1(参照图7的(b))大的像素的暗场信号降低到设定值V2的散射校正。此外，在本实施方式中，设定值V2与阈值V1相等。此外，阈值V1和设定值V2分别是本发明的“规定的阈值”和“规定的设定值”的一例。
具体地说，图像处理部6对暗场像所包含的多个像素的每个像素实施散射校正。详细地说，图像处理部6通过散射校正使暗场像所包含的多个像素中的、暗场信号比阈值V1大的像素的暗场信号降低到设定值V2(阈值V1)。另外，图像处理部6不对暗场像所包含的多个像素中的、暗场信号为阈值V1以下的像素实施散射校正，不使暗场信号变化。
另外，在本实施方式中，图像处理部6构成为：在轴线AR延伸的方向与光栅方向一致的情况下(参照图3的(a))，对包括纤维束10的被摄体T的暗场像进行散射校正。具体地说，在轴线AR延伸的方向与光栅方向一致的情况下，图像处理部6对多个旋转角度中的各个旋转角度下的暗场像(参照图5的暗场像21a～21d)进行散射校正。
然后，如图8所示，图像处理部6通过对多个旋转角度的各个旋转角度下的被进行了散射校正的暗场像进行重构，来生成3D数据22。此外，3D数据22是本发明的“第一立体暗场像”的一例。
另外，图像处理部6构成为能够受理是否实施散射校正的变更。在本实施方式中，图像处理部6构成为：在轴线AR延伸的方向与光栅方向正交的情况下(参照图2的(a))不进行散射校正。
另外，在本实施方式中，图像处理部6构成为：通过散射校正，使由纤维束10a引起的暗场信号(参照图7的(b)的实线)和由纤维束10b引起的暗场信号(参照图7的(b)的虚线)的合计值降低到设定值V2(阈值V1)。此外，设定值V2(阈值V1)是由用户预先设定的值。
另外，在本实施方式中，图像处理部6构成为能够受理设定值V2的变更。具体地说，在用户在未图示的操作部中指示了设定值V2的变更的情况下，由控制部7(参照图1)基于用户的指示来变更图像处理部6的设定值V2。
图9的(a)是轴线AR延伸的方向与光栅方向一致的情况(参照图3的(a))下的、未实施散射校正时的3D数据21(参照图3的(b))的沿规定方向的截面像(比较例)。另外，图9的(b)是轴线AR延伸的方向与光栅方向一致的情况(参照图3的(a))下的、实施了散射校正时的3D数据22(参照图8)的沿规定方向的截面像。在图9的(a)所示的比较例中，在截面像内显著地显现出纤维束10b的像(斜线部分)，纤维束10a(以及裂纹12)的像没有显现。另一方面，在图9的(b)所示的例子中，在图9的(a)中显现出的纤维束10b的像没有显现，而显著地显现出纤维束10a的像(斜线部分)。此外，在图9的(b)所示的例子中，没有显现裂纹12的像。
另外，在图10所示的例子中示出与图9不同种类的被摄体T的图像。图10的(a)是轴线AR延伸的方向与光栅方向一致的情况(参照图3的(a))下的、未实施散射校正时的3D数据21(参照图3的(b)和图12的(a))的沿规定方向的截面像(比较例)。另外，图10的(b)是轴线AR延伸的方向与光栅方向一致的情况(参照图3的(a))下的、实施了散射校正时的3D数据22(参照图8和图12的(b))的沿规定方向的截面像。在图10的(a)所示的比较例中，在截面像内显现出裂纹12(条纹状的部分)的像，但与图9的(a)同样地，显著地显现出纤维束10b的像(斜线部分)，难以视觉识别裂纹12的像。另一方面，在图10的(b)所示的例子中，在图10的(a)中显现出的纤维束10b的像未显现，而显著地显现出裂纹12的像(斜线部分)。此外，在图10的(b)所示的例子中，没有显著地显现纤维束10a的像。
另外，在图11中示出沿着与图10不同的方向的截面像。在图11的(a)所示的比较例中，在截面像内显著地显现出由纤维束10b导致的X射线的散射所引起的伪影(噪声)(从被摄体T出到外侧的放射状的线)。另外，在图11的(b)中，与图11的(a)的比较例相比，抑制了由纤维束10b导致的X射线的散射所引起的伪影(噪声)的产生。
另外，在图12中示出图10和图11的被摄体T的3D数据。如图12的(a)所示，在未实施散射校正的情况下的3D数据21中显著地显现出纤维束10b的立体像。另外，如图12的(b)所示，在实施了散射校正的情况下的3D数据22中显著地显现出裂纹12(以及纤维束10a)的立体像。
在此，在本实施方式中，如图13所示，图像处理部6构成为：在从正面观察(从Z方向观察)第一光栅2和第二光栅3时、轴线AR延伸的方向与光栅方向大致正交的情况下(参照图2的(a))，对多个旋转角度的各个旋转角度下的暗场像(参照图5的暗场像20a～20d)不进行散射校正而进行重构，由此生成3D数据20(参照图2的(b))。而且，图像处理部6构成为将3D数据22和3D数据20进行组合来获取被摄体T的立体像23。即，图像处理部6构成为将3D数据22和3D数据20进行合成来获取(生成)一个立体像23。
(被摄体的摄像)
接着，参照图14对本实施方式的X射线成像装置100拍摄被摄体T的处理的流程进行说明。
在步骤S1中，使光栅方向与轴线AR延伸的方向一致。在该情况下，纤维束10a被设置为沿光栅方向延伸，纤维束10b被设置为沿着从正面观察第一光栅2和第二光栅3时与光栅方向正交的方向延伸。
接着，在步骤S2中，控制部7借助旋转机构8使被摄体T旋转规定角度(例如1度)。
接着，在步骤S3中，控制部7一边借助光栅移动机构9使第一光栅2进行步进移动一边拍摄被摄体T。接着，在步骤S4中，图像处理部6生成被摄体T的暗场像。
接着，在步骤S5中，由图像处理部6对在步骤S4中生成的暗场像实施散射校正。
接着，在步骤S6中，控制部7判定旋转机构8是否使被摄体T旋转了360度。在被摄体T没有旋转360度的情况下，返回到步骤S2。在被摄体T旋转了360度的情况下，进入步骤S7。
在步骤S7中，图像处理部6对在各旋转角度下拍摄到的暗场像进行重构，来生成被摄体T的三维数据(3D数据22(参照图8))。此外，在该情况下，图像处理部6也可以对所生成的3D数据22进行至少包括平滑化处理的校正处理。平滑化处理例如包括基于由高斯滤波器等进行的滤波处理的平滑化。
接着，在步骤S8中，通过使多个光栅进行转动(旋转)，来成为光栅方向与轴线AR延伸的方向正交的状态。在该情况下，纤维束10b被设置为沿着光栅方向延伸，纤维束10a被设置为沿着从正面观察第一光栅2和第二光栅3时与光栅方向正交的方向延伸。
接着，在步骤S9中，控制部7借助旋转机构8使被摄体T旋转规定角度(例如1度)。
接着，在步骤S10中，控制部7一边借助光栅移动机构9使第一光栅2进行步进移动一边拍摄被摄体T。接着，在步骤S11中，图像处理部6生成被摄体T的暗场像。
接着，在步骤S12中，控制部7判定旋转机构8是否使被摄体T旋转了360度。在被摄体T没有旋转360度的情况下，返回到步骤S9。在被摄体T旋转了360度的情况下，进入步骤S13。
在步骤S13中，由图像处理部6对在各旋转角度下拍摄到的暗场像(参照图5的暗场像20a～20d)进行重构，来生成被摄体T的三维数据(3D数据20(参照图2的(b)))。此外，在该情况下，图像处理部6也可以对所生成的3D数据20进行包括平滑化处理的校正处理。
在步骤S14中，图像处理部6将在步骤S7中生成的3D数据22和在步骤13中生成的3D数据20进行组合(合成)，来获取被摄体T的立体像23。
在步骤S15中，图像处理部6从在步骤S14中获取到的立体像23提取规定的特征量。作为在该情况下提取出的特征量，是纤维束10的编织高度、相邻的纤维束10之间的间隙的大小、纤维束10的宽度、纤维束10的长度、纤维束10的曲率、纤维束10的壁厚、裂纹12的长度、裂纹12的深度、裂纹12的数量以及裂纹12的曲率等。另外，也可以将纤维束10的编织的规则性等作为特征量来提取。
(本实施方式的效果)
在本实施方式中，能够获得如下的效果。
在本实施方式中，如上所述，X射线成像装置100构成为：图像处理部6实施使暗场像的多个像素中的、暗场信号比阈值V1大的像素的暗场信号降低到设定值V2(阈值V1)的散射校正。由此，在散射校正前，即使在暗场信号并非相对于与旋转角度对应地变化的(X射线的)透射长度进行指数函数式的变化的情况下，也能够通过优化散射校正中的设定值V2并使暗场信号降低到适当的值来使暗场信号相对于旋转角度(透射长度)的变化接近于指数函数式的变化。其结果，能够抑制由X射线的散射角大引起的伪影(噪声)的产生，并且能够抑制重构出的被摄体T的像变得不清晰。
另外，在本实施方式中，如上所述，X射线成像装置100构成为：图像处理部6在散射校正中将暗场信号比阈值V1大的像素的暗场信号设定为阈值V1。由此，与将设定值V2和阈值V1设为不同的值的情况相比，能够抑制由用户决定设定值V2所耗费的作业负担的增加。
另外，在本实施方式中，如上所述，X射线成像装置100构成为：在旋转机构8的轴线AR延伸的方向与光栅方向一致的情况下，图像处理部6对包括纤维束10的被摄体T的暗场像进行散射校正。在此，一般情况下，在暗场像中显著地显现沿光栅方向延伸的纤维束10。但是，在旋转机构8的轴线AR延伸的方向与光栅方向一致的情况下，除了显现沿光栅方向延伸的纤维束10以外，还显现沿着与光栅方向正交的方向延伸的纤维束10。这是因为，由于通过旋转机构8使被摄体T旋转，因此根据旋转角度不同，存在X射线对于沿着与光栅方向正交的方向延伸的纤维束10的透射长度变长的情况(角度)。因此，由沿着与光栅方向正交的方向延伸的纤维束10引起的X射线的散射角变大。因而，在旋转机构8的轴线AR延伸的方向与光栅方向一致的情况下，通过对包括纤维束10的被摄体T的暗场像进行散射校正，能够更有效地抑制由沿着与光栅方向正交的方向延伸的纤维束10导致的X射线的散射所引起的伪影(噪声)的产生。其结果，能够更有效地抑制重构出的被摄体T的像变得不清晰。
另外，在本实施方式中，如上所述，X射线成像装置100构成为：图像处理部6通过散射校正使由纤维束10a引起的暗场信号和由纤维束10b引起的暗场信号的合计值降低到设定值V2(阈值V1)。由此，与分别独立地通过散射校正使由纤维束10a引起的暗场信号和由纤维束10b引起的暗场信号降低的情况不同，能够通过一次的散射校正来控制暗场信号。其结果，能够抑制散射校正所需的时间变长。
另外，在本实施方式中，如上所述，X射线成像装置100构成为：纤维束10a和纤维束10b中的一方被设置为沿光栅方向延伸，纤维束10a和纤维束10b中的另一方被设置为沿着从正面观察第一光栅2和第二光栅3时与光栅方向正交的方向延伸。由此，通过散射校正，能够抑制由纤维束10a和纤维束10b中的沿着与光栅方向正交的方向延伸的一方导致的X射线的散射所引起的伪影(噪声)的产生。其结果，能够抑制纤维束10a和纤维束10b中的沿光栅方向延伸的另一方的像变得不清晰。
另外，在本实施方式中，如上所述，X射线成像装置100构成为：在旋转机构8的轴线AR延伸的方向与光栅方向一致的情况下，图像处理部6对多个旋转角度的各个旋转角度下的暗场像进行散射校正，并且对多个旋转角度的各个旋转角度下的被进行了散射校正的暗场像进行重构，由此生成3D数据22。由此，能够在多个旋转角度的各个旋转角度下使暗场信号的变化接近指数函数式的衰减的变化，因此能够抑制3D数据22变得不清晰。
另外，在本实施方式中，如上所述，X射线成像装置100构成为：在从正面观察第一光栅2和第二光栅3时轴线AR延伸的方向与光栅方向大致正交的情况下，对多个旋转角度的各个旋转角度下的暗场像不进行散射校正而进行重构，由此生成3D数据20，并且将3D数据22和3D数据20进行组合来获取被摄体T的立体像23。在此，一般情况下，在暗场像中，在从正面观察第一光栅2和第二光栅3时轴线AR延伸的方向与光栅方向大致正交的情况下，在暗场像中显著地显现沿光栅方向延伸的纤维束10。这是由于，在通过旋转机构8使被摄体T旋转的情况下，对于沿着与光栅方向正交的方向(沿着轴线AR的方向)延伸的纤维束10来说，无论在哪个旋转角度下都具有针对纤维束10的暗场灵敏度，并且X射线的透射长度不会根据旋转角度而大幅地变化。由此，在3D数据20中显著地显示的纤维束10与在3D数据22中显著地显示的纤维束10成为互不相同(沿互不相同的方向延伸)的纤维束10。因而，通过将3D数据20和3D数据22进行组合，用户能够在立体像23中视觉识别在3D数据20中显著地显示的纤维束10与在3D数据22中显著地显示的纤维束10的位置关系。
另外，在本实施方式中，如上所述，X射线成像装置100构成为：图像处理部6能够受理设定值V2的变更。由此，能够根据被摄体T的种类来变更设定值V2，因此能够更适当地实施散射校正。
(变形例)
此外，应该认为此次公开的实施方式在所有方面均为例示性的而非限制性的。本发明的范围由权利要求书示出，并非由上述实施方式的说明示出，还包含与权利要求书同等的含义和范围内的所有变更(变形例)。
例如，在上述实施方式中示出了基于暗场像的暗场信号的大小来实施散射校正的例子，但本发明并不限于此。例如，也可以基于暗场像的像素值进行散射校正。在此，在X射线的散射大的情况下，暗场像的像素值(明亮度)变小。因而，在基于像素值进行散射校正的情况下，图像处理部6构成为在像素值小于规定的阈值的情况下进行散射校正。此外，在该情况下，像素值是本发明的“基于X射线的散射强度的值”。
另外，在上述实施方式中示出了规定的设定值(设定值V2)与规定的阈值(阈值V1)相等的例子，但本发明并不限于此。例如，规定的设定值(设定值V2)与规定的阈值(阈值V1)也可以互不相同。在该情况下，图像处理部6也可以构成为基于所生成的暗场像来估计最适合作为规定的设定值(设定值V2)的值。
另外，在上述实施方式中示出了在旋转机构8的旋转轴(轴线AR)延伸的方向与第一光栅构成部分(X射线相位变化部2b)及第二光栅构成部分(X射线吸收部3b)各自延伸的方向(光栅方向)一致的情况下进行散射校正的例子，但本发明并不限于此。例如，也可以在旋转机构8的旋转轴(轴线AR)延伸的方向与第一光栅构成部分(X射线相位变化部2b)及第二光栅构成部分(X射线吸收部3b)各自延伸的方向(光栅方向)错开规定的角度(例如10度)的情况下进行散射校正。
另外，在上述实施方式中示出了第一纤维束(纤维束10a)与第二纤维束(纤维束10b)大致正交的例子，但本发明并不限于此。例如，第一纤维束(纤维束10a)与第二纤维束(纤维束10b)也可以不正交而交叉。
另外，在上述实施方式中示出了在第二立体暗场数据(3D数据20)的生成中不实施散射校正的例子，但本发明并不限于此。例如，也可以在第二立体暗场数据(3D数据20)的生成中实施散射校正。
另外，在上述实施方式中示出了在生成第一立体暗场数据(3D数据22)之后生成第二立体暗场数据(3D数据20)的例子，但本发明并不限于此。例如，也可以在生成了第二立体暗场数据(3D数据20)之后生成第一立体暗场数据(3D数据22)。
另外，在上述实施方式中示出了图像处理部6将第一立体暗场数据(3D数据22)和第二立体暗场数据(3D数据20)进行组合(合成)的例子，但本发明并不限于此。例如，也可以除了第一立体暗场数据(3D数据22)和第二立体暗场数据(3D数据20)以外，还组合旋转轴(轴线AR)延伸的方向与第一光栅构成部分(X射线相位变化部2b)及第二光栅构成部分(X射线吸收部3b)各自延伸的方向(光栅方向)不正交而交叉的情况下的立体暗场数据。
另外，在上述实施方式中示出了设置有第三光栅4的例子，但本发明并不限于此。也可以不设置第三光栅4。
另外，在上述实施方式中示出了通过条纹扫描法生成暗场像的例子，但本发明并不限于此。例如，也可以通过使第一光栅2、第二光栅3以及第三光栅4中的任一个光栅在与X射线的光轴方向正交的平面上旋转的方法(所谓的莫尔单张拍摄方法)来生成暗场像。
另外，在上述实施方式中示出了第一光栅2为相位光栅的例子，但本发明并不限于此。例如，第一光栅2也可以是吸收光栅。
另外，在上述实施方式中示出了以下例子：旋转机构8通过使被摄体T旋转来使被摄体T相对于摄像系统相对地旋转，但本发明并不限于此。例如，也可以构成为通过使摄像系统旋转来使被摄体T相对于摄像系统相对地旋转。
另外，在上述实施方式中示出了以下例子：通过使多个光栅转动(旋转)来使多个光栅的光栅方向与被摄体T的朝向相对地变更，但本发明并不限于此。例如，也可以构成为通过使被摄体T转动(旋转)来使多个光栅的光栅方向与被摄体T的朝向相对地变更。
另外，在上述实施方式中示出了使第一光栅2在光栅面内进行步进移动的例子，但本发明并不限于此。也可以使多个光栅中的任一个光栅进行步进移动。
另外，在上述实施方式中示出了将碳纤维强化塑料(CFRP)作为被摄体T进行拍摄的例子，但本发明并不限于此。例如，也可以将玻璃纤维强化塑料(GFRP)等用作被摄体T。如果在要拍摄的被摄体内包括纤维束，则可以是任何被摄体。
附图标记说明
1：X射线源；2：第一光栅；2b：X射线相位变化部(第一光栅构成部分)；3：第二光栅；3b：X射线吸收部(第二光栅构成部分)；5：检测器；6：图像处理部；8：旋转机构；10：纤维束；10a：纤维束(第一纤维束)；10b：纤维束(第二纤维束)；20：3D数据(第二立体暗场数据)；20a、20b、20c、20d：暗场像(旋转轴延伸的方向与第一光栅构成部分及第二光栅构成部分延伸的方向大致正交的情况下的暗场像)；21a、21b、21c、21d：暗场像(旋转轴延伸的方向与第一光栅构成部分及第二光栅构成部分延伸的方向一致的情况下的暗场像)；22：3D数据(第一立体暗场数据)；23：立体像；100：X射线成像装置；200：摄像系统；AR：轴线(旋转轴)；T：被摄体；V1：阈值(规定的阈值)；V2：设定值(规定的设定值)。

Claims (8)

1.一种X射线成像装置，具备：
X射线源；
检测器，其检测从所述X射线源照射的X射线；
多个光栅，所述多个光栅设置在所述X射线源与所述检测器之间，包括第一光栅和第二光栅；
旋转机构，其使被摄体与摄像系统相对地旋转，所述摄像系统由所述X射线源、所述检测器以及所述多个光栅构成；以及
图像处理部，其基于在通过所述旋转机构进行了旋转的情况下的多个旋转角度的各个旋转角度下由所述检测器检测出的X射线的强度分布，来生成由X射线的散射引起的暗场像，
其中，所述图像处理部构成为实施使以下像素的基于X射线的散射强度的值降低到规定的设定值的散射校正：所述暗场像的多个像素中的、所述基于X射线的散射强度的值比规定的阈值大的所述像素。
2.根据权利要求1所述的X射线成像装置，其特征在于，
所述图像处理部构成为：在所述散射校正中，将所述基于X射线的散射强度的值比所述规定的阈值大的所述像素的所述基于X射线的散射强度的值设定为所述规定的阈值。
3.根据权利要求1所述的X射线成像装置，其特征在于，
所述被摄体包括纤维束，
所述第一光栅具有第一光栅构成部分，
所述第二光栅具有沿着所述第一光栅构成部分所延伸的方向延伸的第二光栅构成部分，
所述图像处理部构成为：在所述旋转机构的旋转轴延伸的方向与所述第一光栅构成部分及所述第二光栅构成部分各自延伸的方向一致的情况下，对包括所述纤维束的所述被摄体的所述暗场像进行所述散射校正。
4.根据权利要求3所述的X射线成像装置，其特征在于，
所述纤维束包括沿规定方向延伸的第一纤维束以及沿着与所述规定方向不同的方向延伸的第二纤维束，
所述图像处理部构成为：通过所述散射校正，使基于由所述第一纤维束引起的所述X射线的散射强度的值和基于由所述第二纤维束引起的所述X射线的散射强度的值的合计值降低到所述规定的设定值。
5.根据权利要求4所述的X射线成像装置，其特征在于，
所述第一纤维束和所述第二纤维束中的一方被设置为沿着所述第一光栅构成部分和所述第二光栅构成部分各自延伸的方向延伸，所述第一纤维束和所述第二纤维束中的另一方被设置为沿着从正面观察所述第一光栅和所述第二光栅时与所述第一光栅构成部分及所述第二光栅构成部分各自延伸的方向正交的方向延伸。
6.根据权利要求3所述的X射线成像装置，其特征在于，
所述图像处理部构成为：在所述旋转机构的旋转轴延伸的方向与所述第一光栅构成部分及所述第二光栅构成部分各自延伸的方向一致的情况下，对所述多个旋转角度的各个旋转角度下的所述暗场像进行所述散射校正，并且对所述多个旋转角度的各个旋转角度下的被进行了所述散射校正的所述暗场像进行重构，由此生成第一立体暗场数据。
7.根据权利要求6所述的X射线成像装置，其特征在于，
所述图像处理部构成为能够受理是否实施所述散射校正的变更，并且所述图像处理部构成为：在从正面观察所述第一光栅和所述第二光栅时、所述旋转轴延伸的方向与所述第一光栅构成部分及所述第二光栅构成部分各自延伸的方向大致正交的情况下，对所述多个旋转角度的各个旋转角度下的所述暗场像不进行所述散射校正而进行重构，由此生成第二立体暗场数据，并且将所述第一立体暗场数据和所述第二立体暗场数据进行组合来获取所述被摄体的立体像。
8.根据权利要求1所述的X射线成像装置，其特征在于，
所述图像处理部构成为能够受理所述规定的设定值的变更。
CN201880082441.7A 2018-01-19 2018-11-06 X射线成像装置 Active CN111566471B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2018-007624 2018-01-19
JP2018007624 2018-01-19
PCT/JP2018/041152 WO2019142451A1 (ja) 2018-01-19 2018-11-06 Ｘ線イメージング装置

Publications (2)

Publication Number Publication Date
CN111566471A CN111566471A (zh) 2020-08-21
CN111566471B true CN111566471B (zh) 2023-06-16

Family

ID=67301445

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201880082441.7A Active CN111566471B (zh) 2018-01-19 2018-11-06 X射线成像装置

Country Status (4)

Country Link
US (1) US11154270B2 (zh)
JP (1) JP6897801B2 (zh)
CN (1) CN111566471B (zh)
WO (1) WO2019142451A1 (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7136033B2 (ja) * 2019-07-26 2022-09-13 株式会社島津製作所 Ｘ線位相差撮影装置
JP7454541B2 (ja) 2021-11-24 2024-03-22 株式会社リガク 構造情報取得方法及び構造情報取得装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101510298A (zh) * 2009-03-17 2009-08-19 西北工业大学 一种ct伪影的综合校正方法
CN101943668A (zh) * 2009-07-07 2011-01-12 清华大学 X射线暗场成像系统和方法
JP2015043931A (ja) * 2013-08-29 2015-03-12 コニカミノルタ株式会社 医用画像処理装置及び位相画像生成方法
JP2015118074A (ja) * 2013-12-20 2015-06-25 キヤノン株式会社 Ｘ線断層撮影装置およびｘ線断層撮影方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102221565B (zh) * 2010-04-19 2013-06-12 清华大学 X射线源光栅步进成像系统与成像方法
JP6187298B2 (ja) * 2014-02-14 2017-08-30 コニカミノルタ株式会社 Ｘ線撮影システム及び画像処理方法
EP3314576B1 (en) * 2015-06-26 2019-11-27 Koninklijke Philips N.V. Robust reconstruction for dark-field and phase contrast ct
JP6602630B2 (ja) * 2015-10-05 2019-11-06 株式会社日立ハイテクサイエンス Ｘ線検査装置及びｘ線検査方法
JP6613988B2 (ja) * 2016-03-30 2019-12-04 コニカミノルタ株式会社 放射線撮影システム
DE102016223797A1 (de) * 2016-11-30 2018-05-30 Technische Universität München Röntgen-CT-Verfahren, Probenhalter und Röntgen-CT-Vorrichtung

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101510298A (zh) * 2009-03-17 2009-08-19 西北工业大学 一种ct伪影的综合校正方法
CN101943668A (zh) * 2009-07-07 2011-01-12 清华大学 X射线暗场成像系统和方法
JP2015043931A (ja) * 2013-08-29 2015-03-12 コニカミノルタ株式会社 医用画像処理装置及び位相画像生成方法
JP2015118074A (ja) * 2013-12-20 2015-06-25 キヤノン株式会社 Ｘ線断層撮影装置およびｘ線断層撮影方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Design of Acquisition Schemes and Setup Geometry for Anisotropic X-ray Dark-Field Tomography (AXDT);Sharma Y.等;《Scientific Reports》;20170609;第1-10页 *
Laminate fibre structure characterisation of carbon fibre-reinforced polymers by X-ray scatter dark field imaging with a grating interferometer;Revol Vincent等;《NDT&E International》;20130508;第64-71页 *

Also Published As

Publication number Publication date
WO2019142451A1 (ja) 2019-07-25
JPWO2019142451A1 (ja) 2020-11-26
JP6897801B2 (ja) 2021-07-07
US20210059630A1 (en) 2021-03-04
CN111566471A (zh) 2020-08-21
US11154270B2 (en) 2021-10-26

Similar Documents

Publication Publication Date Title
CN110049725B (zh) X射线相位成像装置
CN109425624B (zh) X射线成像装置
JP2019144177A (ja) Ｘ線撮影システム
DE102010016997B4 (de) Inspektionssystem und Verfahren mit Mehrfachbildphasenverschiebungsanalyse
JP2012200567A (ja) 放射線撮影システム及び放射線撮影方法
CN111566471B (zh) X射线成像装置
JP7006784B2 (ja) Ｘ線イメージング装置
JP6673188B2 (ja) Ｘ線位相撮影装置
WO2020066135A1 (ja) Ｘ線位相イメージング装置
JP6881682B2 (ja) Ｘ線イメージング装置
CN112304986A (zh) X射线成像装置和x射线成像方法
CN109283200B (zh) X射线相位成像装置和含有纤维的材料的缺陷检测方法
CN108720855B (zh) Ｘ射线相位成像装置
WO2020246220A1 (ja) 放射線撮影システム及び拡大吸収コントラスト画像生成方法
WO2020090168A1 (ja) Ｘ線位相差撮影システム
WO2020039660A1 (ja) Ｘ線位相イメージング装置
JP2021143958A (ja) Ｘ線位相撮像システム、ｘ線画像処理装置、および、ｘ線画像処理方法
CN109646027A (zh) 乳腺成像系统及其光路装置
WO2013051647A1 (ja) 放射線撮影装置及び画像処理方法
JP2021085829A (ja) Ｘ線位相イメージング装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant