JP6897801B2 - Ｘ線イメージング装置 - Google Patents

Ｘ線イメージング装置 Download PDF

Info

Publication number
JP6897801B2
JP6897801B2 JP2019565730A JP2019565730A JP6897801B2 JP 6897801 B2 JP6897801 B2 JP 6897801B2 JP 2019565730 A JP2019565730 A JP 2019565730A JP 2019565730 A JP2019565730 A JP 2019565730A JP 6897801 B2 JP6897801 B2 JP 6897801B2
Authority
JP
Japan
Prior art keywords
lattice
fiber bundle
ray
image
grid
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
JP2019565730A
Other languages
English (en)
Other versions
JPWO2019142451A1 (ja
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 JPWO2019142451A1 publication Critical patent/JPWO2019142451A1/ja
Application granted granted Critical
Publication of JP6897801B2 publication Critical patent/JP6897801B2/ja
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 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/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 or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements 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 or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4291—Arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating
    • 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/46—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 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/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 or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4035—Arrangements for generating radiation specially adapted for radiation diagnosis the source being combined with a filter or grating

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Molecular Biology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Optics & Photonics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Human Computer Interaction (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Toxicology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Description

本発明は、Ｘ線イメージング装置に関し、特に、被写体の暗視野像を生成する画像処理部を備えるＸ線イメージング装置に関する。
従来、被写体の暗視野像を生成する画像処理部を備えるＸ線イメージング装置が知られている。このようなＸ線イメージング装置は、たとえば、M Bech, et al., “Quantitative x-ray dark-field computed tomography” 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のＸ線イメージング装置は、位相格子と、検出器と、画像処理装置とを備えている。Ｘ線源からのＸ線は、被写体により散乱され、位相格子を通過して検出器に照射される。検出器により検出されるＸ線の干渉強度に基づいて、画像処理装置により暗視野像が生成される。また、被写体が３６０度回転された場合における各々の角度において生成された暗視野像を用いて画像の再構成が行われる。
上記M Bech, et al., “Quantitative x-ray dark-field computed tomography” 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
しかしながら、上記M Bech, et al., “Quantitative x-ray dark-field computed tomography” PHYSICS IN MEDICINE AND BIOLOGY 55(2010) P.5529-P.5539に記載されているような従来のＸ線イメージング装置では、Ｘ線の散乱角が大きい場合などにおいては、Ｘ線の干渉強度が透過長に対して指数関数的に変化しないため、ＦＢＰによる再構成を適切に行うことができない。すなわち、ＦＢＰにより再構成された暗視野像（３次元データ）においてアーチファクト（ノイズ）が発生するという不都合がある。この場合、再構成された被写体の像が不鮮明になるという問題点がある。
この発明は、上記のような課題を解決するためになされたものであり、この発明の１つの目的は、再構成された被写体の像が不鮮明になるのを抑制することが可能なＸ線イメージング装置を提供することである。
上記目的を達成するために、この発明の一の局面におけるＸ線イメージング装置は、Ｘ線源と、Ｘ線源から照射されたＸ線を検出する検出器と、Ｘ線源と検出器との間に設けられ、第１格子および第２格子を含む複数の格子と、被写体と、Ｘ線源と検出器と複数の格子とによって構成される撮像系とを相対的に回転させる回転機構と、回転機構により回転された場合における複数の回転角度の各々において、検出器により検出されたＸ線の強度分布に基づいて、Ｘ線の散乱に起因する暗視野像を生成する画像処理部とを備え、画像処理部は、暗視野像における複数の画素のうち、Ｘ線の散乱強度に基づく値が所定の閾値よりも大きい画素のＸ線の散乱強度に基づく値を、所定の設定値に低減させる散乱補正を実施するように構成されている。なお、「暗視野像」とは、物体の小角散乱に基づくＶｉｓｉｂｉｌｉｔｙの変化によって得られる、Ｖｉｓｉｂｉｌｉｔｙ像のことである。また、暗視野像は、小角散乱像とも呼ばれる。「Ｖｉｓｉｂｉｌｉｔｙ」とは、鮮明度のことである。
この発明の一の局面におけるＸ線イメージング装置では、上記のように、Ｘ線の散乱強度に基づく値が所定の閾値よりも大きい画素のＸ線の散乱強度に基づく値を、所定の設定値に低減させるように構成されている。これにより、散乱補正前において、Ｘ線の散乱強度に基づく値が、回転角度に対応して変化する（Ｘ線の）透過長に対して指数関数的に変化しない場合でも、散乱補正における所定の設定値を適切化してＸ線の散乱強度に基づく値を適切な値まで低減させることによって、Ｘ線の散乱強度に基づく値の回転角度（透過長）に対する変化を指数関数的な変化に近づけることができる。その結果、Ｘ線の散乱角が大きいことに起因するアーチファクト（ノイズ）の発生を抑制することができるとともに、再構成された被写体の像が不鮮明になるのを抑制することができる。
上記一の局面におけるＸ線イメージング装置において、好ましくは、画像処理部は、散乱補正において、Ｘ線の散乱強度に基づく値が所定の閾値よりも大きい画素のＸ線の散乱強度に基づく値を、所定の閾値に設定するように構成されている。このように構成すれば、所定の設定値と所定の閾値とを別個の値にする場合に比べて、ユーザによる所定の設定値の決定にかかる作業負担の増加を抑制することができる。
上記一の局面におけるＸ線イメージング装置において、好ましくは、被写体は、繊維束を含み、第１格子は、第１格子構成部分を有し、第２格子は、第１格子構成部分が延びる方向に沿って延びる第２格子構成部分を有し、画像処理部は、回転機構の回転軸が延びる方向と、第１格子構成部分および第２格子構成部分の各々が延びる方向とが揃っている場合に、繊維束を含む被写体の暗視野像に対して散乱補正を行うように構成されている。ここで、一般的に、暗視野像において、第１格子構成部分および第２格子構成部分の各々が延びる方向に沿って延びる繊維束が顕著に現われる。しかし、回転機構の回転軸が延びる方向と、第１格子構成部分および第２格子構成部分の各々が延びる方向とが揃っている場合では、第１格子構成部分および第２格子構成部分の各々が延びる方向に沿って延びる繊維束に加えて、第１格子構成部分および第２格子構成部分の各々が延びる方向と直交する方向に延びる繊維束が現われる。これは、回転機構により被写体が回転されることによって、回転角度によっては、第１格子構成部分および第２格子構成部分の各々が延びる方向と直交する方向に延びる繊維束に対するＸ線の透過長が長くなる場合（角度）があるからである。このため、第１格子構成部分および第２格子構成部分の各々が延びる方向と直交する方向に延びる繊維束によるＸ線の散乱角が大きくなる。したがって、回転機構の回転軸が延びる方向と、第１格子構成部分および第２格子構成部分の各々が延びる方向とが揃っている場合に、繊維束を含む被写体の暗視野像に対して散乱補正を行うことによって、第１格子構成部分および第２格子構成部分の各々が延びる方向と直交する方向に延びる繊維束によるＸ線の散乱に起因するアーチファクト（ノイズ）の発生をより効果的に抑制することができる。その結果、再構成された被写体の像が不鮮明になるのをより効果的に抑制することができる。
この場合、好ましくは、繊維束は、所定の方向に延びる第１繊維束と、所定の方向とは異なる方向に延びる第２繊維束とを含み、画像処理部は、散乱補正によって、第１繊維束に起因するＸ線の散乱強度に基づく値と、第２繊維束に起因するＸ線の散乱強度に基づく値との合計値を、所定の設定値に低減させるように構成されている。このように構成すれば、第１繊維束に起因するＸ線の散乱強度に基づく値、および、第２繊維束に起因するＸ線の散乱強度に基づく値をそれぞれ別個に散乱補正により低減させる場合と異なり、１回の散乱補正によってＸ線の散乱強度に基づく値の制御をすることができる。その結果、散乱補正に要する時間が長くなるのを抑制することができる。
上記第１繊維束に起因するＸ線の散乱強度に基づく値と第２繊維束に起因するＸ線の散乱強度に基づく値との合計値を所定の設定値に低減させるＸ線イメージング装置において、好ましくは、第１繊維束および第２繊維束のうちの一方は、第１格子構成部分および第２格子構成部分の各々が延びる方向に沿って延びるように設けられ、第１繊維束および第２繊維束のうちの他方は、第１格子および第２格子を正面から見て第１格子構成部分および第２格子構成部分の各々が延びる方向と直交する方向に沿って延びるように設けられている。このように構成すれば、散乱補正によって、第１格子構成部分および第２格子構成部分の各々が延びる方向と直交する方向に延びる第１繊維束および第２繊維束のうちの一方によるＸ線の散乱に起因するアーチファクト（ノイズ）の発生を抑制することができる。その結果、第１格子構成部分および第２格子構成部分の各々が延びる方向に沿って延びる第１繊維束および第２繊維束のうちの他方の像が不鮮明になるのを抑制することができる。
上記回転軸が延びる方向と第１格子構成部分および第２格子構成部分の各々が延びる方向とが揃っている場合に散乱補正を行うＸ線イメージング装置において、好ましくは、画像処理部は、回転機構の回転軸が延びる方向と、第１格子構成部分および第２格子構成部分の各々が延びる方向とが揃っている場合に、複数の回転角度の各々における暗視野像に対して散乱補正を行うとともに、散乱補正が行われた複数の回転角度の各々における暗視野像を再構成することにより第１立体暗視野データを生成するように構成されている。このように構成すれば、複数の回転角度の各々において、暗視野シグナル（Ｘ線の散乱強度に基づく値）を指数関数的減衰の変化に近づけることができるので、第１立体暗視野データが不鮮明になるのを抑制することができる。
この場合、好ましくは、画像処理部は、散乱補正の実施の有無の変更を受け付け可能に構成されており、第１格子および第２格子を正面から見て回転軸が延びる方向と第１格子構成部分および第２格子構成部分の各々が延びる方向とが略直交している場合に、複数の回転角度の各々における暗視野像に対して散乱補正を行わずに再構成することにより第２立体暗視野データを生成するとともに、第１立体暗視野データと第２立体暗視野データとを組み合わせて、被写体の立体像を取得するように構成されている。ここで、一般的に、暗視野像において、第１格子および第２格子を正面から見て回転軸が延びる方向と第１格子構成部分および第２格子構成部分の各々が延びる方向とが略直交している場合、第１格子構成部分および第２格子構成部分の各々が延びる方向に沿って延びる繊維束が暗視野像において顕著に現われる。これは、回転機構により被写体が回転される場合に、第１格子構成部分および第２格子構成部分の各々が延びる方向と直交する方向（回転軸に沿った方向）に延びる繊維束に対して、どの回転角度でも繊維束に対する暗視野感度があるためなのと、Ｘ線の透過長が回転角度によって大きく変化しないからである。これにより、第２立体暗視野データにおいて顕著に表示される繊維束と、第１立体暗視野データにおいて顕著に表示される繊維束とは互いに異なる（方向に延びる）繊維束となる。したがって、第１立体暗視野データと第２立体暗視野データとを組み合わせることによって、ユーザは、第１立体暗視野データに顕著に表示される繊維束と、第２立体暗視野データにおいて顕著に現われる繊維束との位置関係を立体像において視認することができる。
上記一の局面におけるＸ線イメージング装置において、好ましくは、画像処理部は、所定の設定値の変更を受け付け可能に構成されている。このように構成すれば、所定の設定値を被写体の種類に応じて変更することができるので、より適切に散乱補正を実施することができる。
本発明によれば、上記のように、再構成された被写体の像が不鮮明になるのを抑制することができる。
一実施形態によるＸ線イメージング装置の構成を示した模式図である。 図２（ａ）は、回転軸と格子方向とが直交している場合のＸ線イメージング装置の構成を示した構成図である。図２（ｂ）は、図２（ａ）の場合の立体暗視野データである。 図３（ａ）は、回転軸と格子方向とが揃っている場合のＸ線イメージング装置の構成を示した構成図である。図３（ｂ）は、図３（ａ）の場合（散乱補正なし）の立体暗視野データ（比較例）である。 一実施形態によるＸ線イメージング装置における、Ｘ線の光軸に対する被写体の配置角度とＸ線の散乱度合を説明するための模式図（ａ）〜（ｄ）である。 一実施形態によるＸ線イメージング装置によって被写体を回転させながら撮像している際の暗視野像の模式図（ａ）〜（ｄ）である。 一実施形態によるＸ線イメージング装置によって生成された暗視野像のサイノグラムである。（図６（ａ）は、図２（ａ）の場合の図である。図６（ｂ）は、図３（ａ）の場合（散乱補正なし）の図である。） 一実施形態によるＸ線イメージング装置による第１繊維束に起因する暗視野信号と第２繊維束に起因する暗視野信号との関係を説明するための図である。（図７（ａ）は、図２（ａ）の場合の図である。図７（ｂ）は、図３（ａ）の場合（散乱補正なし）の図である。） 図３（ａ）のＸ線イメージング装置によって生成された暗視野像に対して散乱補正を行うとともに再構成して生成された立体暗視野データである。 一実施形態によるＸ線イメージング装置による立体暗視野データの断面像である。（図９（ａ）は、図３（ｂ）の断面像（比較例）である。図９（ｂ）は、図８の断面像である。） 一実施形態によるＸ線イメージング装置による（図９とは異なる被写体を用いて生成された）立体暗視野データの断面像である。（図１０（ａ）は、図３（ｂ）の断面像（比較例）である。図１０（ｂ）は、図８の断面像である。） 一実施形態によるＸ線イメージング装置による立体暗視野データの（図１０とは異なる方向に沿った方向の）断面像である。（図１１（ａ）は、図３（ｂ）の断面像（比較例）である。図１１（ｂ）は、図８の断面像である。） 一実施形態によるＸ線イメージング装置による立体暗視野データを示した図である。（図１２（ａ）は、図３（ｂ）に示した立体暗視野データ（比較例）である。図１２（ｂ）は、図８に示した立体暗視野データである。） 一実施形態によるＸ線イメージング装置による図８に示した立体暗視野データと図２（ｂ）に示した立体暗視野データとを組み合わせて取得された立体像を示した図である。 一実施形態によるＸ線イメージング装置による撮像方法を説明するためのフローチャートである。
以下、本発明を具体化した実施形態を図面に基づいて説明する。
［本実施形態］
図１〜図１４を参照して、本実施形態によるＸ線イメージング装置１００の構成について説明する。
（Ｘ線イメージング装置の構成）
図１に示すように、Ｘ線イメージング装置１００は、被写体Ｔを通過したＸ線の拡散（散乱）を利用して、被写体Ｔの内部を画像化する装置である。また、Ｘ線イメージング装置１００は、タルボ（Ｔａｌｂｏｔ）効果を利用して、被写体Ｔの内部を画像化する装置である。Ｘ線イメージング装置１００は、たとえば、非破壊検査用途では、物体の内部の画像化に用いることが可能である。
被写体Ｔは、内部に繊維束１０（図２参照）を含む。被写体Ｔは、たとえば、繊維束１０として炭素繊維が用いられ、母材として樹脂１１（図２参照）が用いられる炭素繊維強化プラスチック（ＣＦＲＰ）である。なお、繊維束とは、繊維が多数集まって束状になったものである。本実施形態では、繊維束１０は、多数の繊維によって板状に形成されたものである。
図１は、Ｘ線イメージング装置１００をＸ方向から見た図である。図１に示すように、Ｘ線イメージング装置１００は、Ｘ線源１と、第１格子２と、第２格子３と、第３格子４と、検出器５と、画像処理部６と、制御部７と、回転機構８と、格子移動機構９とを備えている。なお、本明細書において、Ｘ線源１から第１格子２に向かう方向をＺ２方向、その逆向きの方向をＺ１方向とする。また、Ｚ方向と直交する面内の左右方向をＸ方向とし、紙面の奥に向かう方向をＸ２方向、紙面の手前側に向かう方向をＸ１方向とする。また、Ｚ方向と直交する面内の上下方向をＹ方向とし、上方向をＹ１方向、下方向をＹ２方向とする。
Ｘ線源１は、高電圧が印加されることにより、Ｘ線を発生させるとともに、発生されたＸ線をＺ２方向に向けて照射するように構成されている。
第１格子２は、Ｙ方向に所定の周期（ピッチ）ｄ１で配列される複数のスリット２ａおよび、Ｘ線位相変化部２ｂを有している。各スリット２ａおよびＸ線位相変化部２ｂはそれぞれ、直線状に延びるように形成されている。また、各スリット２ａおよびＸ線位相変化部２ｂはそれぞれ、平行に延びるように形成されている。第１格子２は、いわゆる位相格子である。なお、Ｘ線位相変化部２ｂは、請求の範囲の「第１格子構成部分」の一例である。
第１格子２は、Ｘ線源１と、第２格子３との間に配置されており、Ｘ線源１からＸ線が照射される。第１格子２は、タルボ効果により、第１格子２の自己像（図示せず）を形成するために設けられている。なお、可干渉性を有するＸ線が、スリットが形成された格子を通過すると、格子から所定の距離（タルボ距離）離れた位置に、格子の像（自己像）が形成される。これをタルボ効果という。
第２格子３は、Ｙ方向に所定の周期（ピッチ）ｄ２で配列される複数のＸ線透過部３ａおよびＸ線吸収部３ｂを有する。Ｘ線吸収部３ｂは、Ｘ線位相変化部２ｂが延びる方向に沿って延びている。各Ｘ線透過部３ａおよびＸ線吸収部３ｂはそれぞれ、直線状に延びるように形成されている。また、各Ｘ線透過部３ａおよびＸ線吸収部３ｂはそれぞれ、平行に延びるように形成されている。第２格子３は、いわゆる、吸収格子である。第１格子２、第２格子３はそれぞれ異なる役割を持つ格子であるが、スリット２ａおよびＸ線透過部３ａはそれぞれＸ線を透過させる。また、Ｘ線吸収部３ｂはＸ線を遮蔽する役割を担っており、Ｘ線位相変化部２ｂはスリット２ａとの屈折率の違いによってＸ線の位相を変化させる。なお、Ｘ線吸収部３ｂは、請求の範囲の「第２格子構成部分」の一例である。
第２格子３は、第１格子２と検出器５との間に配置されており、第１格子２を通過したＸ線が照射される。また、第２格子３は、第１格子２からタルボ距離離れた位置に配置される。第２格子３は、第１格子２の自己像と干渉して、検出器５の検出表面上にモアレ縞（図示せず）を形成する。
第３格子４は、Ｙ方向に所定の周期（ピッチ）ｄ３で配列される複数のＸ線透過部４ａおよびＸ線吸収部４ｂを有する。各Ｘ線透過部４ａおよびＸ線吸収部４ｂはそれぞれ、直線状に延びるように形成されている。また、各Ｘ線透過部４ａおよびＸ線吸収部４ｂはそれぞれ、平行に延びるように形成されている。第３格子４は、いわゆる、マルチスリットである。
第３格子４は、Ｘ線源１と第１格子２との間に配置されている。第３格子４は、各Ｘ線透過部４ａを通過したＸ線を線光源とすることにより、Ｘ線源１からのＸ線を多点光源化するように構成されている。３枚の格子（第１格子２、第２格子３、および、第３格子４）のピッチと格子間の距離とが一定の条件を満たすことにより、Ｘ線源１から照射されるＸ線の可干渉性を高めることが可能である。これにより、Ｘ線源１の管球の焦点サイズが大きくても干渉強度を保持できる。
検出器５は、Ｘ線を検出するとともに、検出されたＸ線を電気信号に変換し、変換された電気信号を画像信号として読み取るように構成されている。検出器５は、たとえば、ＦＰＤ（Ｆｌａｔ Ｐａｎｅｌ Ｄｅｔｅｃｔｏｒ）である。検出器５は、複数の変換素子（図示せず）と複数の変換素子上に配置された画素電極（図示せず）とにより構成されている。複数の変換素子および画素電極は、所定の周期（画素ピッチ）で、Ｘ方向およびＹ方向にアレイ状に配列されている。また、検出器５は、取得した画像信号を、画像処理部６に出力するように構成されている。
画像処理部６は、検出器５から出力された画像信号に基づいて、吸収像（図示せず）を生成するように構成されている。また、画像処理部６は、検出器５により検出されたＸ線の強度分布に基づいて暗視野像（図５参照）を生成するように構成されている。ここで、吸収像とは、被写体ＴによるＸ線の吸収の差によって生じるコントラストを画像化したものである。また、暗視野像とは、被写体Ｔの内部にある微細構造によるＸ線の屈折によって生じるコントラストを画像化したものである。また、画像処理部６は、回転機構８を回転させながら（複数の回転角度の各々において）撮像された複数の吸収像および複数の暗視野像をそれぞれ再構成することにより、３次元吸収像（図示せず）および３次元暗視野像（図２（ｂ）等参照）を生成する。また、画像処理部６は、たとえば、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）や画像処理用に構成されたＦＰＧＡ（Ｆｉｅｌｄ−Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）などのプロセッサを含む。
制御部７は、回転機構８を介して、被写体Ｔと、Ｘ線源１、検出器５、および、複数の格子（第１格子２、第２格子３、および、第３格子４）によって構成される撮像系２００を相対的に回転させるように構成されている。また、制御部７は、格子移動機構９を介して、第１格子２を格子面内において格子方向と直交する方向にステップ移動させるように構成されている。Ｘ線イメージング装置１００では、第１格子２を一定周期間隔に走査することにより得られた複数のモアレ縞（画像）から再構成画像を取得する手法（縞走査法）が用いられている。また、制御部７は、たとえば、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）などのプロセッサを含む。なお、格子方向は、請求の範囲の「第１格子構成部分および第２格子構成部分の各々が延びる方向」の一例である。
回転機構８は、制御部７からの信号に基づいて、被写体Ｔと撮像系２００とを相対的に回転させるように構成されている。具体的には、回転機構８は、被写体Ｔを軸線ＡＲ周りに回転させることにより、撮像系２００に対して被写体Ｔを相対的に回転させるように構成されている。図１では、軸線ＡＲが延びる方向（図１ではＹ方向）と、複数の格子の格子方向（図１ではＸ方向）とが直交している状態を図示しているが、複数の格子の格子方向は所定の方向（たとえばＹ方向）に変化させることが可能である。なお、格子方向とは、格子の格子パターンが延びる方向である。また、格子パターンとは、各格子のスリット２ａ、Ｘ線位相変化部２ｂ、Ｘ線透過部３ａ、および、Ｘ線吸収部３ｂなどのことである。また、回転機構８は、たとえば、モータなどによって駆動される回転ステージを含む。なお、軸線ＡＲは、請求の範囲の「回転軸」の一例である。
格子移動機構９は、制御部７からの信号に基づいて、第１格子２を格子面内（ＸＹ面内）において格子方向と直交する方向（図１ではＹ方向）にステップ移動させるように構成されている。具体的には、格子移動機構９は、第１格子２の周期ｄ１をｎ分割し、ｄ１／ｎずつ第１格子２をステップ移動させる。格子移動機構９は、少なくとも第１格子２の１周期ｄ１分、第１格子２をステップ移動させるように構成されている。なお、ｎは正の整数であり、たとえば、９などである。また、格子移動機構９は、たとえば、ステッピングモータやピエゾアクチュエータなどを含む。
図２（ａ）に示すように、繊維束１０は、被写体Ｔ内において、複数の異なる方向に延びる繊維束１０を含む。図２（ａ）に示す例では、被写体Ｔは、第１格子２および第２格子３を正面から見て（Ｚ方向から見て）格子方向と直交する方向（Ｙ方向）に延びる繊維束１０ａと、格子方向（Ｘ方向）に延びる繊維束１０ｂとが編み込まれた構造となっている。なお、図２（ａ）に示す例では、繊維束１０ａおよび繊維束１０ｂを区別して表示しているが、実際にはシート状に織り込まれた（または、編み込まれた）同種の炭素繊維である。また、Ｘ方向およびＹ方向は、それぞれ、請求の範囲の「所定の方向とは異なる方向」および「所定の方向」の一例である。また、繊維束１０ａおよび繊維束１０ｂは、それぞれ、請求の範囲の「第１繊維束」および「第２繊維束」の一例である。
また、被写体Ｔの樹脂１１には、クラック１２が形成されている場合がある。図２に示す例では、クラック１２は、Ｙ方向に延びている。なお、クラック１２は、樹脂１１に形成される空気の層などである。
図２に示す例では、軸線ＡＲが延びる方向（図２（ａ）ではＹ方向）と格子方向（図２（ａ）ではＸ方向）とが直交している。この場合、繊維束１０のうち格子方向に延びる繊維束（図２（ａ）では繊維束１０ｂ）の像が暗視野像において現われる。この場合に生成された複数の回転角度における暗視野像（図５の暗視野像２０ａ〜２０ｄ参照）を再構成することにより、繊維束１０ｂの立体像が鮮明に表示される３Ｄデータ２０（図２（ｂ）参照）が生成される。なお、３Ｄデータ２０は、請求の範囲の「第２立体暗視野像」の一例である。
図３に示す比較例では、軸線ＡＲが延びる方向（図３（ａ）ではＹ方向）と格子方向（図３（ａ）ではＹ方向）とが揃っている。この場合、繊維束１０のうち格子方向に延びている繊維束（図３（ａ）では繊維束１０ａ）およびクラック１２の像に加えて、格子方向と直交する方向（Ｘ方向）に延びている繊維束（図３（ａ）では繊維束１０ｂ）が暗視野像において現われる。この場合に生成された複数の回転角度における暗視野像（図５の暗視野像２１ｂ〜２１ｄ参照）を再構成することにより、繊維束１０ａおよびクラック１２の立体像に加えて、繊維束１０ｂの立体像が現われる３Ｄデータ２１（図３（ｂ）参照）が生成される。
ここで、図４を参照して、軸線ＡＲが延びる方向と格子方向とが揃っている場合（図３（ａ）参照）において、暗視野像（図５の暗視野像２１ｂ〜２１ｄ参照）に、格子方向と直交する方向（Ｘ方向）に延びている繊維束（図３（ａ）では繊維束１０ｂ）が現われる理由について説明する。
図４（ａ）は、ＸＺ平面において、被写体Ｔの繊維束１０ｂがＸ線の光軸と直交するように配置されている場合の例を示す模式図である。図４（ｂ）は、図４（ａ）における４００−４００線に沿う断面図である。図４（ｃ）は、ＸＺ平面において、被写体Ｔの繊維束１０ｂがＸ線の光軸に対して斜めになるように配置されている場合の例を示す模式図である。図４（ｄ）は、図４（ｃ）における５００−５００線に沿う断面図である。
ここで、Ｘ線の光軸と直交する方向に延びる繊維束１０にＸ線が入射した場合、繊維束１０によってＸ線が散乱（拡散）する。具体的には、被写体Ｔ内の繊維と樹脂１１（図２参照）との界面において、繊維と樹脂１１との屈折率の差によってＸ線が屈折する。被写体Ｔ内には、繊維束１０が積層されているため、また、繊維束１０は多数の繊維から構成されているため、多数の繊維を通過することによって多重の屈折が起こり、Ｘ線が散乱（拡散）する。
図４（ａ）に示す例では、繊維束１０ｂがＸ線の光軸と直交する配置となっている。すなわち、繊維束１０ｂがＸ方向に延びるように配置されている。したがって、図４（ｂ）に示すように、繊維の断面は略円形となっている。図４（ｃ）に示す例では、繊維束１０ｂがＸ線の光軸に対して斜めに配置されている。したがって、図４（ｄ）に示すように、繊維の断面は楕円形である。繊維の断面が楕円形となる場合、繊維の断面が略円形となる場合と比較して、Ｘ線が繊維と樹脂１１との界面を通過する距離が長くなる。また、繊維と樹脂１１との界面はざらざらした形状（微小な凹凸）となっており、Ｘ線が繊維と樹脂１１との界面を通過する距離が長くなる分だけ、Ｘ線の散乱（拡散）が増大する。
図５は、本実施形態によるＸ線イメージング装置１００において被写体Ｔを回転させながら撮像した際のある角度（第１角度θ１、第２角度θ２、第３角度θ３および第４角度θ４）における、軸線ＡＲが延びる方向と格子方向とが直交している場合（図２参照）の暗視野像２０ａ〜２０ｄと、軸線ＡＲが延びる方向と格子方向とが揃っている場合（図３参照）の暗視野像２１ａ〜２１ｄを示す模式図（ａ）〜（ｄ）である。
軸線ＡＲが延びる方向と格子方向とが直交している場合では、第１角度θ１の暗視野像（暗視野像２０ａ）（図５（ａ）参照）において、繊維束１０ｂが鮮明に現われている。また、角度が比較的大きい暗視野像（暗視野像２０ｂ〜２０ｄ）（図５（ｂ）〜（ｄ）参照）において、角度が大きくなるにつれて、繊維束１０ｂによるＸ線の散乱（拡散）（図４参照）が大きくなる。その結果、繊維束１０ｂの暗視野信号（図７（ａ）の実線参照）が大きくなる。なお、この場合、角度が大きくなっても繊維束１０ａによるＸ線の散乱（拡散）は大きく変化しないので、繊維束１０ａ（およびクラック１２）の暗視野信号（図７（ａ）の破線参照）は大きくならない。ここで、暗視野信号とは、暗視野像における所定の画素の画素値をＶとした場合にｆ（１／Ｖ）で求められる値である。なお、上記関数ｆ（ｘ）は、引数ｘの値が大きいほど値が大きくなる関数である。したがって、Ｘ線の散乱強度が大きく、暗視野像の画素値（明るさ）が小さい場合（すなわちＶが大きい場合）は、暗視野信号は大きくなる。なお、暗視野信号は、請求の範囲の「Ｘ線の散乱強度に基づく値」の一例である。
一方、軸線ＡＲが延びる方向と格子方向とが揃っている場合では、角度が比較的小さい第１角度θ１の暗視野像（暗視野像２１ａ）においては、繊維束１０ｂの暗視野信号（図７（ｂ）の実線参照）が比較的小さく、繊維束１０ａ（およびクラック１２）が鮮明に現われている。しかし、角度が大きくなるにつれて、繊維束１０ｂによるＸ線の散乱（拡散）が大きくなるため、繊維束１０ｂの暗視野信号（図７（ｂ）の実線参照）が繊維束１０ａの暗視野信号（図７（ｂ）の破線参照）に対して大きくなる。このため、暗視野像２１ｂ〜２１ｄにおいて、繊維束１０ｂの像が顕著に現われるため、繊維束１０ａ（およびクラック１２）の詳細な構造を把握することは難しくなる。
図６は、本実施形態によるＸ線イメージング装置１００において被写体Ｔを３６０度回転させながら撮像した際のサイノグラムである。図６のサイノグラムは、複数の回転角度の各々における検出器５の出力データを示したものである。軸線ＡＲが延びる方向と格子方向とが直交している場合のサイノグラム（図６（ａ）参照）、および、軸線ＡＲが延びる方向と格子方向とが揃っている場合のサイノグラム（図６（ｂ）参照）の各々に示されているように、所定の回転角度において、繊維束１０ｂによるＸ線の散乱（拡散）が大きくなっている（図６（ａ）および図６（ｂ）の斜線部分参照）ことが分かる。すなわち、図６（ａ）および図６（ｂ）の斜線部分は、暗視野信号が比較的大きい回転角度（図７（ａ）および（ｂ）の実線参照）に対応している。
ここで、本実施形態では、画像処理部６（図１参照）は、暗視野像における複数の画素のうち、暗視野信号が閾値Ｖ１（図７（ｂ）参照）よりも大きい画素の暗視野信号を、設定値Ｖ２に低減させる散乱補正を実施するように構成されている。なお、本実施形態では、設定値Ｖ２は、閾値Ｖ１と等しい。なお、閾値Ｖ１および設定値Ｖ２は、それぞれ、請求の範囲の「所定の閾値」および「所定の設定値」の一例である。
具体的には、画像処理部６は、暗視野像に含まれる複数の画素の各々に対して散乱補正を実施する。詳細には、画像処理部６は、暗視野像に含まれる複数の画素のうち、暗視野信号が閾値Ｖ１よりも大きい画素の暗視野信号を散乱補正により設定値Ｖ２（閾値Ｖ１）に低減させる。また、画像処理部６は、暗視野像に含まれる複数の画素のうち、暗視野信号が閾値Ｖ１以下の画素に対しては散乱補正を実施せず、暗視野信号を変化させない。
また、本実施形態では、画像処理部６は、軸線ＡＲが延びる方向と格子方向とが揃っている場合（図３（ａ）参照）に、繊維束１０を含む被写体Ｔの暗視野像に対して散乱補正を行うように構成されている。具体的には、画像処理部６は、軸線ＡＲが延びる方向と格子方向とが揃っている場合に、複数の回転角度の各々における暗視野像（図５の暗視野像２１ａ〜２１ｄ参照）に対して散乱補正を行う。
そして、図８に示すように、画像処理部６は、散乱補正が行われた複数の回転角度の各々における暗視野像を再構成することにより３Ｄデータ２２を生成する。なお、３Ｄデータ２２は、請求の範囲の「第１立体暗視野像」の一例である。
また、画像処理部６は、散乱補正の実施の有無の変更を受け付け可能に構成されている。本実施形態では、画像処理部６は、軸線ＡＲが延びる方向と格子方向とが直交している場合（図２（ａ）参照）には、散乱補正が行われないように構成されている。
また、本実施形態では、画像処理部６は、散乱補正によって、繊維束１０ａに起因する暗視野信号（図７（ｂ）の実線参照）と、繊維束１０ｂに起因する暗視野信号（図７（ｂ）の破線参照）との合計値を、設定値Ｖ２（閾値Ｖ１）に低減させるように構成されている。なお、設定値Ｖ２（閾値Ｖ１）は、ユーザにより予め設定されている値である。
また、本実施形態では、画像処理部６は、設定値Ｖ２の変更を受け付け可能に構成されている。具体的には、図示しない操作部においてユーザが設定値Ｖ２の変更を指示した場合、ユーザの指示に基づいて制御部７（図１参照）により画像処理部６の設定値Ｖ２が変更される。
図９（ａ）は、軸線ＡＲが延びる方向と格子方向とが揃っている場合（図３（ａ）参照）の、散乱補正が実施されていない場合の３Ｄデータ２１（図３（ｂ）参照）の所定の方向に沿った断面像（比較例）である。また、図９（ｂ）は、軸線ＡＲが延びる方向と格子方向とが揃っている場合（図３（ａ）参照）の、散乱補正が実施されている場合の３Ｄデータ２２（図８参照）の所定の方向に沿った断面像である。図９（ａ）に示す比較例では、断面像内において、繊維束１０ｂの像（斜線部分）が顕著に現われており、繊維束１０ａ（およびクラック１２）の像は現われていない。一方、図９（ｂ）に示す例では、図９（ａ）において現われていた繊維束１０ｂの像が現われず、繊維束１０ａの像（斜線部分）が顕著に現われている。なお、図９（ｂ）に示す例では、クラック１２の像は現われていない。
また、図１０に示す例では、図９とは異なる種類の被写体Ｔの画像を示す。図１０（ａ）は、軸線ＡＲが延びる方向と格子方向とが揃っている場合（図３（ａ）参照）の、散乱補正が実施されていない場合の３Ｄデータ２１（図３（ｂ）および図１２（ａ）参照）の所定の方向に沿った断面像（比較例）である。また、図１０（ｂ）は、軸線ＡＲが延びる方向と格子方向とが揃っている場合（図３（ａ）参照）の、散乱補正が実施されている場合の３Ｄデータ２２（図８および図１２（ｂ）参照）の所定の方向に沿った断面像である。図１０（ａ）に示す比較例では、断面像内において、クラック１２（筋状の部分）の像は現われてはいるが、図９（ａ）と同様に繊維束１０ｂの像（斜線部分）が顕著に現われており、クラック１２の像が視認しにくくなっている。一方、図１０（ｂ）に示す例では、図１０（ａ）において現われていた繊維束１０ｂの像が現われず、クラック１２の像（斜線部分）が顕著に現われている。なお、図１０（ｂ）に示す例では、繊維束１０ａの像は顕著に現われていない。
また、図１１では、図１０とは異なる方向に沿った断面像を示す。図１１（ａ）に示す比較例では、断面像内において、繊維束１０ｂによるＸ線の散乱に起因するアーチファクト（ノイズ）（被写体Ｔから外側に出ている放射状の線）が顕著に現われている。また、図１１（ｂ）では、図１１（ａ）の比較例に比べて、繊維束１０ｂによるＸ線の散乱に起因するアーチファクト（ノイズ）の発生が抑制されている。
また、図１２では、図１０および図１１の被写体Ｔの３Ｄデータを示す。図１２（ａ）に示すように、散乱補正が実施されていない場合の３Ｄデータ２１には、繊維束１０ｂの立体像が顕著に現われている。また、図１２（ｂ）に示すように、散乱補正が実施されている場合の３Ｄデータ２２には、クラック１２（および繊維束１０ａ）の立体像が顕著に現われている。
ここで、本実施形態では、図１３に示すように、画像処理部６は、第１格子２および第２格子３を正面から見て（Ｚ方向から見て）軸線ＡＲが延びる方向と格子方向とが略直交している場合（図２（ａ）参照）に、複数の回転角度の各々における暗視野像（図５の暗視野像２０ａ〜２０ｄ参照）に対して散乱補正を行わずに再構成することにより３Ｄデータ２０（図２（ｂ）参照）を生成するように構成されている。そして、画像処理部６は、３Ｄデータ２２と３Ｄデータ２０とを組み合わせて、被写体Ｔの立体像２３を取得するように構成されている。すなわち、画像処理部６は、３Ｄデータ２２と３Ｄデータ２０とを合成して、１つの立体像２３を取得（生成）するように構成されている。
（被写体の撮像）
次に、図１４を参照して、本実施形態によるＸ線イメージング装置１００による被写体Ｔを撮像する処理の流れについて説明する。
ステップＳ１において、格子方向と軸線ＡＲが延びる方向とを揃える。この場合、繊維束１０ａは、格子方向に沿って延びるように設けられ、繊維束１０ｂは、第１格子２および第２格子３を正面から見て格子方向と直交する方向に沿って延びるように設けられている。
次に、ステップＳ２において、制御部７により、回転機構８を介して、被写体Ｔが所定角度（たとえば１度）回転される。
次に、ステップＳ３において、制御部７により、格子移動機構９を介して第１格子２がステップ移動されながら、被写体Ｔが撮像される。次に、ステップＳ４において、画像処理部６により、被写体Ｔの暗視野像が生成される。
次に、ステップＳ５において、画像処理部６により、ステップＳ４において生成された暗視野像に対して散乱補正が実施される。
次に、ステップＳ６において、制御部７により、回転機構８が被写体Ｔを３６０度回転させたかどうかが判定される。被写体Ｔが３６０度回転していない場合、ステップＳ２に戻る。被写体Ｔが３６０度回転している場合、ステップＳ７へ進む。
ステップＳ７において、画像処理部６により、各回転角度において撮像された暗視野像が再構成され、被写体Ｔの３次元データ（３Ｄデータ２２（図８参照））が生成される。なお、この場合、画像処理部６により、生成された３Ｄデータ２２に対して少なくとも平滑化処理を含む補正処理が行われてもよい。平滑化処理は、たとえば、ガウシアンフィルタなどによるフィルター処理による平滑化を含む。
次に、ステップＳ８において、複数の格子が回動（回転）されることにより、格子方向と軸線ＡＲが延びる方向とが直交した状態にする。この場合、繊維束１０ｂは、格子方向に沿って延びるように設けられ、繊維束１０ａは、第１格子２および第２格子３を正面から見て格子方向と直交する方向に沿って延びるように設けられている。
次に、ステップＳ９において、制御部７により、回転機構８を介して、被写体Ｔが所定角度（たとえば１度）回転される。
次に、ステップＳ１０において、制御部７により、格子移動機構９を介して第１格子２がステップ移動されながら、被写体Ｔが撮像される。次に、ステップＳ１１において、画像処理部６により、被写体Ｔの暗視野像が生成される。
次に、ステップＳ１２において、制御部７により、回転機構８が被写体Ｔを３６０度回転させたかどうかが判定される。被写体Ｔが３６０度回転していない場合、ステップＳ９に戻る。被写体Ｔが３６０度回転している場合、ステップＳ１３へ進む。
ステップＳ１３において、画像処理部６により、各回転角度において撮像された暗視野像（図５の暗視野像２０ａ〜２０ｄ参照）が再構成され、被写体Ｔの３次元データ（３Ｄデータ２０（図２（ｂ）参照））が生成される。なお、この場合、画像処理部６により、生成された３Ｄデータ２０に対して平滑化処理を含む補正処理が行われてもよい。
ステップＳ１４において、画像処理部６により、ステップＳ７において生成された３Ｄデータ２２と、ステップ１３において生成された３Ｄデータ２０とが組み合わされて（合成されて）、被写体Ｔの立体像２３が取得される。
ステップＳ１５において、画像処理部６により、ステップＳ１４において取得された立体像２３から所定の特徴量が抽出される。この場合に抽出される特徴量としては、繊維束１０の編み込みの高さ、隣接する繊維束１０の間の隙間の大きさ、繊維束１０の幅、繊維束１０の長さ、繊維束１０の曲率、繊維束１０の肉厚、クラック１２の長さ、クラック１２の深さ、クラック１２の数量、および、クラック１２の曲率などである。また、繊維束１０の編み込みの規則性などを特徴量として抽出してもよい。
（本実施形態の効果）
本実施形態では、以下のような効果を得ることができる。
本実施形態では、上記のように、画像処理部６が、暗視野像における複数の画素のうち、暗視野信号が閾値Ｖ１よりも大きい画素の暗視野信号を、設定値Ｖ２（閾値Ｖ１）に低減させる散乱補正を実施するように、Ｘ線イメージング装置１００を構成する。これにより、散乱補正前において、暗視野信号が、回転角度に対応して変化する（Ｘ線の）透過長に対して指数関数的に変化しない場合でも、散乱補正における設定値Ｖ２を適切化して暗視野信号を適切な値まで低減させることによって、暗視野信号の回転角度（透過長）に対する変化を指数関数的な変化に近づけることができる。その結果、Ｘ線の散乱角が大きいことに起因するアーチファクト（ノイズ）の発生を抑制することができるとともに、再構成された被写体Ｔの像が不鮮明になるのを抑制することができる。
また、本実施形態では、上記のように、画像処理部６が、散乱補正において、暗視野信号が閾値Ｖ１よりも大きい画素の暗視野信号を、閾値Ｖ１に設定するように、Ｘ線イメージング装置１００を構成する。これにより、設定値Ｖ２と閾値Ｖ１とを別個の値にする場合に比べて、ユーザによる設定値Ｖ２の決定にかかる作業負担の増加を抑制することができる。
また、本実施形態では、上記のように、画像処理部６が、回転機構８の軸線ＡＲが延びる方向と、格子方向とが揃っている場合に、繊維束１０を含む被写体Ｔの暗視野像に対して散乱補正を行うように、Ｘ線イメージング装置１００を構成する。ここで、一般的に、暗視野像において、格子方向に沿って延びる繊維束１０が顕著に現われる。しかし、回転機構８の軸線ＡＲが延びる方向と、格子方向とが揃っている場合では、格子方向に沿って延びる繊維束１０に加えて、格子方向と直交する方向に延びる繊維束１０が現われる。これは、回転機構８により被写体Ｔが回転されることによって、回転角度によっては、格子方向と直交する方向に延びる繊維束１０に対するＸ線の透過長が長くなる場合（角度）があるからである。このため、格子方向と直交する方向に延びる繊維束１０によるＸ線の散乱角が大きくなる。したがって、回転機構８の軸線ＡＲが延びる方向と、格子方向とが揃っている場合に、繊維束１０を含む被写体Ｔの暗視野像に対して散乱補正を行うことによって、格子方向と直交する方向に延びる繊維束１０によるＸ線の散乱に起因するアーチファクト（ノイズ）の発生をより効果的に抑制することができる。その結果、再構成された被写体Ｔの像が不鮮明になるのをより効果的に抑制することができる。
また、本実施形態では、上記のように、画像処理部６が、散乱補正によって、繊維束１０ａに起因する暗視野信号と、繊維束１０ｂに起因する暗視野信号との合計値を、設定値Ｖ２（閾値Ｖ１）に低減させるように、Ｘ線イメージング装置１００を構成する。これにより、繊維束１０ａに起因する暗視野信号、および、繊維束１０ｂに起因する暗視野信号をそれぞれ別個に散乱補正により低減させる場合と異なり、１回の散乱補正によって暗視野信号の制御をすることができる。その結果、散乱補正に要する時間が長くなるのを抑制することができる。
また、本実施形態では、上記のように、繊維束１０ａおよび繊維束１０ｂのうちの一方が、格子方向に沿って延びるように設けられ、繊維束１０ａおよび繊維束１０ｂのうちの他方が、第１格子２および第２格子３を正面から見て格子方向と直交する方向に沿って延びるように設けられるように、Ｘ線イメージング装置１００を構成する。これにより、散乱補正によって、格子方向と直交する方向に延びる繊維束１０ａおよび繊維束１０ｂのうちの一方によるＸ線の散乱に起因するアーチファクト（ノイズ）の発生を抑制することができる。その結果、格子方向に沿って延びる繊維束１０ａおよび繊維束１０ｂのうちの他方の像が不鮮明になるのを抑制することができる。
また、本実施形態では、上記のように、画像処理部６が、回転機構８の軸線ＡＲが延びる方向と、格子方向とが揃っている場合に、複数の回転角度の各々における暗視野像に対して散乱補正を行うとともに、散乱補正が行われた複数の回転角度の各々における暗視野像を再構成することにより３Ｄデータ２２を生成するように、Ｘ線イメージング装置１００を構成する。これにより、複数の回転角度の各々において、暗視野信号を指数関数的減衰の変化に近づけることができるので、３Ｄデータ２２が不鮮明になるのを抑制することができる。
また、本実施形態では、上記のように、第１格子２および第２格子３を正面から見て軸線ＡＲが延びる方向と格子方向とが略直交している場合に、複数の回転角度の各々における暗視野像に対して散乱補正を行わずに再構成することにより３Ｄデータ２０を生成するとともに、３Ｄデータ２２と３Ｄデータ２０とを組み合わせて、被写体Ｔの立体像２３を取得するように、Ｘ線イメージング装置１００を構成する。ここで、一般的に、暗視野像において、第１格子２および第２格子３を正面から見て軸線ＡＲが延びる方向と格子方向とが略直交している場合、格子方向に沿って延びる繊維束１０が暗視野像において顕著に現われる。これは、回転機構８により被写体Ｔが回転される場合に、格子方向と直交する方向（軸線ＡＲに沿った方向）に延びる繊維束１０に対して、どの回転角度でも繊維束１０に対する暗視野感度があるためなのと、Ｘ線の透過長が回転角度によって大きく変化しないからである。これにより、３Ｄデータ２０において顕著に表示される繊維束１０と、３Ｄデータ２２において顕著に表示される繊維束１０とは互いに異なる（方向に延びる）繊維束１０となる。したがって、３Ｄデータ２０と３Ｄデータ２２とを組み合わせることによって、ユーザは、３Ｄデータ２０に顕著に表示される繊維束１０と、３Ｄデータ２２において顕著に表示される繊維束１０との位置関係を立体像２３において視認することができる。
また、本実施形態では、上記のように、画像処理部６が、設定値Ｖ２の変更を受け付け可能に構成されるように、Ｘ線イメージング装置１００を構成する。これにより、設定値Ｖ２を被写体Ｔの種類に応じて変更することができるので、より適切に散乱補正を実施することができる。
（変形例）
なお、今回開示された実施形態は、すべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した実施形態の説明ではなく、請求の範囲によって示され、さらに請求の範囲と均等の意味および範囲内でのすべての変更（変形例）が含まれる。
たとえば、上記実施形態では、暗視野像の暗視野信号の大きさに基づいて、散乱補正を実施する例を示したが、本発明はこれに限られない。たとえば、暗視野像の画素値に基づいて、散乱補正を行ってもよい。ここで、Ｘ線の散乱が大きい場合は暗視野像の画素値（明るさ）は小さくなる。したがって、画素値に基づいて散乱補正を行う場合、画像処理部６は、画素値が所定の閾値よりも小さい場合に散乱補正を行うように構成される。なお、この場合、画素値は、請求の範囲の「Ｘ線の散乱強度に基づく値」である。
また、上記実施形態では、所定の設定値（設定値Ｖ２）と所定の閾値（閾値Ｖ１）とが等しい例を示したが、本発明はこれに限られない。たとえば、所定の設定値（設定値Ｖ２）と所定の閾値（閾値Ｖ１）とが互いに異なっていてもよい。この場合、画像処理部６は、所定の設定値（設定値Ｖ２）として最適な値を、生成された暗視野像に基づいて推定するように構成されていてもよい。
また、上記実施形態では、回転機構８の回転軸（軸線ＡＲ）が延びる方向と、第１格子構成部分（Ｘ線位相変化部２ｂ）および第２格子構成部分（Ｘ線吸収部３ｂ）の各々が延びる方向（格子方向）とが揃っている場合に散乱補正を行う例を示したが、本発明はこれに限られない。たとえば、回転機構８の回転軸（軸線ＡＲ）が延びる方向と、第１格子構成部分（Ｘ線位相変化部２ｂ）および第２格子構成部分（Ｘ線吸収部３ｂ）の各々が延びる方向（格子方向）とが、所定の角度（たとえば１０度）ずれている場合に散乱補正を行ってもよい。
また、上記実施形態では、第１繊維束（繊維束１０ａ）と第２繊維束（繊維束１０ｂ）とが略直交している例を示したが、本発明はこれに限られない。たとえば、第１繊維束（繊維束１０ａ）と第２繊維束（繊維束１０ｂ）とが直交せずに交差していてもよい。
また、上記実施形態では、第２立体暗視野データ（３Ｄデータ２０）の生成において散乱補正を実施しない例を示したが、本発明はこれに限られない。たとえば、第２立体暗視野データ（３Ｄデータ２０）の生成において散乱補正を実施してもよい。
また、上記実施形態では、第１立体暗視野データ（３Ｄデータ２２）を生成した後に、第２立体暗視野データ（３Ｄデータ２０）を生成する例を示したが、本発明はこれに限られない。たとえば、第２立体暗視野データ（３Ｄデータ２０）を生成した後に、第１立体暗視野データ（３Ｄデータ２２）を生成してもよい。
また、上記実施形態では、画像処理部６が、第１立体暗視野データ（３Ｄデータ２２）と第２立体暗視野データ（３Ｄデータ２０）とを組み合わせる（合成する）例を示したが、本発明はこれに限られない。たとえば、第１立体暗視野データ（３Ｄデータ２２）および第２立体暗視野データ（３Ｄデータ２０）に加えて、回転軸（軸線ＡＲ）が延びる方向と第１格子構成部分（Ｘ線位相変化部２ｂ）および第２格子構成部分（Ｘ線吸収部３ｂ）の各々が延びる方向（格子方向）とが直交せずに交差している場合の立体暗視野データを組み合わせてもよい。
また、上記実施形態では、第３格子４が設けられている例を示したが、本発明はこれに限られない。第３格子４が設けられていなくてもよい。
また、上記実施形態では、縞走査法によって暗視野像を生成する例を示したが、本発明はこれに限られない。たとえば、第１格子２、第２格子３、または、第３格子４のうちのいずれかを、Ｘ線の光軸方向に直交する平面上において回転させる手法（いわゆる、モアレ１枚撮り手法）によって、暗視野像を生成してもよい。
また、上記実施形態では、第１格子２が位相格子である例を示したが、本発明はこれに限られない。たとえば、第１格子２は吸収格子であってもよい。
また、上記実施形態では、回転機構８は、被写体Ｔを回転させることにより、撮像系に対して被写体Ｔを相対的に回転させる例を示したが、本発明はこれに限られない。たとえば、撮像系を回転させることにより、撮像系に対して被写体Ｔを相対的に回転させるように構成されていてもよい。
また、上記実施形態では、複数の格子を回動（回転）させることにより、複数の格子の格子方向と被写体Ｔの向きとを相対的に変更する例を示したが、本発明はこれに限られない。たとえば、被写体Ｔを回動（回転）させることにより、複数の格子の格子方向と被写体Ｔの向きとを相対的に変更するように構成されていてもよい。
また、上記実施形態では、第１格子２を格子面内においてステップ移動させる例を示したが、本発明はこれに限られない。複数の格子のうち、いずれの格子をステップ移動させてもよい。
また、上記実施形態では、被写体Ｔとして、炭素繊維強化プラスチック（ＣＦＲＰ）を撮像する例を示したが、本発明はこれに限られない。たとえば、ガラス繊維強化プラスチック（ＧＦＲＰ）などを被写体Ｔとして用いてもよい。撮像する被写体内に繊維束が含まれていれば、どのような被写体であってもよい。
１ Ｘ線源
２ 第１格子
２ｂ Ｘ線位相変化部（第１格子構成部分）
３ 第２格子
３ｂ Ｘ線吸収部（第２格子構成部分）
５ 検出器
６ 画像処理部
８ 回転機構
１０ 繊維束
１０ａ 繊維束（第１繊維束）
１０ｂ 繊維束（第２繊維束）
２０ ３Ｄデータ（第２立体暗視野データ）
２０ａ、２０ｂ、２０ｃ、２０ｄ 暗視野像（回転軸が延びる方向と第１格子構成部分および第２格子構成部分の延びる方向とが略直交している場合の暗視野像）
２１ａ、２１ｂ、２１ｃ、２１ｄ 暗視野像（回転軸が延びる方向と第１格子構成部分および第２格子構成部分の延びる方向とが揃っている場合の暗視野像）
２２ ３Ｄデータ（第１立体暗視野データ）
２３ 立体像
１００ Ｘ線イメージング装置
２００ 撮像系
ＡＲ 軸線（回転軸）
Ｔ 被写体
Ｖ１ 閾値（所定の閾値）
Ｖ２ 設定値（所定の設定値）

Claims (8)

  1. Ｘ線源と、
    前記Ｘ線源から照射されたＸ線を検出する検出器と、
    前記Ｘ線源と前記検出器との間に設けられ、第１格子および第２格子を含む複数の格子と、
    被写体と、前記Ｘ線源と前記検出器と前記複数の格子とによって構成される撮像系とを相対的に回転させる回転機構と、
    前記回転機構により回転された場合における複数の回転角度の各々において、前記検出器により検出されたＸ線の強度分布に基づいて、Ｘ線の散乱に起因する暗視野像を生成する画像処理部とを備え、
    前記画像処理部は、前記暗視野像における複数の画素のうち、Ｘ線の散乱強度に基づく値が所定の閾値よりも大きい前記画素の前記Ｘ線の散乱強度に基づく値を、所定の設定値に低減させる散乱補正を実施するように構成されている、Ｘ線イメージング装置。
  2. 前記画像処理部は、前記散乱補正において、前記Ｘ線の散乱強度に基づく値が前記所定の閾値よりも大きい前記画素の前記Ｘ線の散乱強度に基づく値を、前記所定の閾値に設定するように構成されている、請求項１に記載のＸ線イメージング装置。
  3. 前記被写体は、繊維束を含み、
    前記第１格子は、第１格子構成部分を有し、
    前記第２格子は、前記第１格子構成部分が延びる方向に沿って延びる第２格子構成部分を有し、
    前記画像処理部は、前記回転機構の回転軸が延びる方向と、前記第１格子構成部分および前記第２格子構成部分の各々が延びる方向とが揃っている場合に、前記繊維束を含む前記被写体の前記暗視野像に対して前記散乱補正を行うように構成されている、請求項１に記載のＸ線イメージング装置。
  4. 前記繊維束は、所定の方向に延びる第１繊維束と、前記所定の方向とは異なる方向に延びる第２繊維束とを含み、
    前記画像処理部は、前記散乱補正によって、前記第１繊維束に起因する前記Ｘ線の散乱強度に基づく値と、前記第２繊維束に起因する前記Ｘ線の散乱強度に基づく値との合計値を、前記所定の設定値に低減させるように構成されている、請求項３に記載のＸ線イメージング装置。
  5. 前記第１繊維束および前記第２繊維束のうちの一方は、前記第１格子構成部分および前記第２格子構成部分の各々が延びる方向に沿って延びるように設けられ、前記第１繊維束および前記第２繊維束のうちの他方は、前記第１格子および前記第２格子を正面から見て前記第１格子構成部分および前記第２格子構成部分の各々が延びる方向と直交する方向に沿って延びるように設けられている、請求項４に記載のＸ線イメージング装置。
  6. 前記画像処理部は、前記回転機構の回転軸が延びる方向と、前記第１格子構成部分および前記第２格子構成部分の各々が延びる方向とが揃っている場合に、前記複数の回転角度の各々における前記暗視野像に対して前記散乱補正を行うとともに、前記散乱補正が行われた前記複数の回転角度の各々における前記暗視野像を再構成することにより第１立体暗視野データを生成するように構成されている、請求項３に記載のＸ線イメージング装置。
  7. 前記画像処理部は、前記散乱補正の実施の有無の変更を受け付け可能に構成されており、前記第１格子および前記第２格子を正面から見て前記回転軸が延びる方向と前記第１格子構成部分および前記第２格子構成部分の各々が延びる方向とが略直交している場合に、前記複数の回転角度の各々における前記暗視野像に対して前記散乱補正を行わずに再構成することにより第２立体暗視野データを生成するとともに、前記第１立体暗視野データと前記第２立体暗視野データとを組み合わせて、前記被写体の立体像を取得するように構成されている、請求項６に記載のＸ線イメージング装置。
  8. 前記画像処理部は、前記所定の設定値の変更を受け付け可能に構成されている、請求項１に記載のＸ線イメージング装置。
JP2019565730A 2018-01-19 2018-11-06 Ｘ線イメージング装置 Active JP6897801B2 (ja)

Applications Claiming Priority (3)

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

Publications (2)

Publication Number Publication Date
JPWO2019142451A1 JPWO2019142451A1 (ja) 2020-11-26
JP6897801B2 true JP6897801B2 (ja) 2021-07-07

Family

ID=67301445

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019565730A Active JP6897801B2 (ja) 2018-01-19 2018-11-06 Ｘ線イメージング装置

Country Status (4)

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

Families Citing this family (3)

* 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 株式会社リガク 構造情報取得方法及び構造情報取得装置
WO2024101253A1 (ja) * 2022-11-10 2024-05-16 コニカミノルタ株式会社 画像処理方法、画像処理装置、画像処理システム及びプログラム

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101510298B (zh) * 2009-03-17 2010-12-29 西北工业大学 一种ct伪影的综合校正方法
CN101943668B (zh) * 2009-07-07 2013-03-27 清华大学 X射线暗场成像系统和方法
CN102221565B (zh) * 2010-04-19 2013-06-12 清华大学 X射线源光栅步进成像系统与成像方法
JP6237001B2 (ja) * 2013-08-29 2017-11-29 コニカミノルタ株式会社 医用画像処理装置及び位相画像生成方法
JP2015118074A (ja) * 2013-12-20 2015-06-25 キヤノン株式会社 Ｘ線断層撮影装置およびｘ線断層撮影方法
JP6187298B2 (ja) * 2014-02-14 2017-08-30 コニカミノルタ株式会社 Ｘ線撮影システム及び画像処理方法
JP6626517B2 (ja) * 2015-06-26 2019-12-25 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 暗視野および位相コントラストｃｔのための堅牢な再構成
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

Also Published As

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

Similar Documents

Publication Publication Date Title
JP6897801B2 (ja) Ｘ線イメージング装置
JP6741080B2 (ja) Ｘ線位相イメージング装置
JP6943090B2 (ja) Ｘ線イメージング装置
EP2469472A1 (en) Computed tomography (ct) image reconstruction method and apparatus
JP2012200567A (ja) 放射線撮影システム及び放射線撮影方法
JP6673188B2 (ja) Ｘ線位相撮影装置
JP7006784B2 (ja) Ｘ線イメージング装置
WO2020095482A1 (ja) Ｘ線位相撮像システム
WO2019220689A1 (ja) Ｘ線イメージング装置
JP7180566B2 (ja) Ｘ線イメージング装置およびｘ線イメージング方法
JP6844461B2 (ja) Ｘ線位相イメージング装置および情報取得手法
JP7163969B2 (ja) Ｘ線位相差撮影システム
JP6798408B2 (ja) Ｘ線位相イメージング装置
JP6992902B2 (ja) Ｘ線位相イメージング装置
JP2011227028A (ja) 断層像再構成方法およびｘ線ｃｔ装置
JP2021143958A (ja) Ｘ線位相撮像システム、ｘ線画像処理装置、および、ｘ線画像処理方法
JP2007151849A (ja) Ｘ線ｃｔ撮影方法およびｘ線ｃｔ装置
KR101685005B1 (ko) 저선량 엑스선 콘빔 ct 영상 장치 및 이를 이용한 영상 생성 방법
JP2019072367A (ja) Ｘ線位相差撮影装置および位相コントラスト画像補正方法
JP2021085829A (ja) Ｘ線位相イメージング装置
JP2012220422A (ja) 断層像再構成方法およびｘ線ｃｔ装置
JP2005331527A (ja) コンピュータ断層撮影方法および装置
JP2012249787A (ja) Ｘ線光学素子及び該ｘ線光学素子を有するｘ線撮像装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200512

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20210511

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210524

R151 Written notification of patent or utility model registration

Ref document number: 6897801

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151