JP2012170645A - X-ray imaging apparatus, and x-ray imaging method - Google Patents
X-ray imaging apparatus, and x-ray imaging method Download PDFInfo
- Publication number
- JP2012170645A JP2012170645A JP2011036114A JP2011036114A JP2012170645A JP 2012170645 A JP2012170645 A JP 2012170645A JP 2011036114 A JP2011036114 A JP 2011036114A JP 2011036114 A JP2011036114 A JP 2011036114A JP 2012170645 A JP2012170645 A JP 2012170645A
- Authority
- JP
- Japan
- Prior art keywords
- ray
- detector
- subject
- distance
- dividing element
- 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.)
- Pending
Links
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Measurement Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
【課題】位相イメージングにより被検体の位相像を取得するに当たり、被検体における吸収情報に対する位相情報の割合が最大限取得可能となるX線撮像装置を提供する。
【解決手段】X線撮像装置であって、X線源と、
前記X線源から照射された発散X線を複数のX線に分割するための複数の開口を具備する、前記X線源と被検体との間に配置された分割素子と、
前記被検体を透過し屈折した前記発散X線によるX線ビームを検出するため、該X線ビームの照射方向の被検体よりも下流に配置された検出器と、
前記検出器により検出されたX線ビームの強度変化の情報を演算処理する演算処理手段と、を備え、
前記X線源と前記分割素子の距離をL1、前記分割素子と前記検出器の距離をL2とするとき、L2/L1の値が1.0より大きく2.4未満の範囲となるように構成されている。
【選択図】 図1An object of the present invention is to provide an X-ray imaging apparatus capable of acquiring the maximum ratio of phase information to absorption information in a subject when acquiring a phase image of the subject by phase imaging.
An X-ray imaging apparatus comprising: an X-ray source;
A splitting element disposed between the X-ray source and the subject, comprising a plurality of openings for splitting the divergent X-rays irradiated from the X-ray source into a plurality of X-rays;
A detector disposed downstream of the subject in the irradiation direction of the X-ray beam in order to detect an X-ray beam by the divergent X-ray transmitted through the subject and refracted;
An arithmetic processing means for arithmetically processing information on an intensity change of the X-ray beam detected by the detector,
When the distance between the X-ray source and the dividing element is L1, and the distance between the dividing element and the detector is L2, the value of L2 / L1 is in the range of more than 1.0 and less than 2.4. Has been.
[Selection] Figure 1
Description
本発明はX線を用いたX線撮像装置、およびX線撮像方法に関する。 The present invention relates to an X-ray imaging apparatus and an X-ray imaging method using X-rays.
放射線を用いた非破壊検査法は工業利用から医療利用まで幅広い分野で用いられている。
例えば、X線は波長が約1pm以上10nm以下の電磁波であり、各物質の質量減弱係数差を用いた吸収コントラスト法はX線の透過能の高さを利用し、金属系材料などの内部亀裂検査、食品に混入された異物や手荷物検査などのセキュリティ分野の用途として実用化されている。
しかしながら、プラスチックやアクリル樹脂などを代表とする高分子材料系は、質量減弱係数差が小さいため吸収コントラスト法では検出・検知することが困難である。
Nondestructive inspection methods using radiation are used in a wide range of fields from industrial use to medical use.
For example, X-rays are electromagnetic waves having a wavelength of about 1 pm to 10 nm, and the absorption contrast method using the difference in mass attenuation coefficient of each substance makes use of the high X-ray transmissivity, so that internal cracks in metallic materials, etc. It has been put to practical use in security fields such as inspection, inspection of foreign matters mixed in food and baggage inspection.
However, polymer material systems such as plastic and acrylic resin are difficult to detect and detect by the absorption contrast method because the difference in mass attenuation coefficient is small.
放射線の吸収によるコントラストが形成され難い軽元素で構成されている被検体や低密度差の境界面を有する被検体に対しては、被検体によるX線の位相変化を検出する位相イメージングが有効である。
このような位相イメージングを用いた方法は、高分子材料系のイメージングや軟組織などの医療系への応用が有力である。
各種X線位相イメージングにおいて、特許文献1や特許文献2ではX線の被検知物による位相変化による屈折効果を利用したX線撮像方法が提案されている。
これらの方法では、X線光学手段によって空間的に分割したX線を被検体に照射し、被検体で屈折したX線を二次元検出器に入射するように設計されている。
このような手法によれば、屈折したX線の屈折量を二次元検出器の強度分布から得ることができる。
Phase imaging that detects X-ray phase changes by the subject is effective for subjects that are composed of light elements that do not easily form contrast due to radiation absorption or that have a low density difference interface. is there.
Such a method using phase imaging is likely to be applied to imaging of polymer materials and medical systems such as soft tissues.
In various X-ray phase imaging,
These methods are designed to irradiate a subject with X-rays spatially divided by X-ray optical means and to irradiate the two-dimensional detector with X-rays refracted by the subject.
According to such a method, the amount of refraction of the refracted X-ray can be obtained from the intensity distribution of the two-dimensional detector.
ところで、特許文献1や特許文献2に記載された方法に限らず、被検体の屈折率を可視化する際、屈折情報から吸収情報を完全に分離して屈折情報だけから被検体の位相像を得ることは困難である。
特許文献1や特許文献2では、被検体の特定空間の平均屈折率を検出することができるが、特定空間の吸収は一様であることが前提条件となっている。
しかしながら、特定空間内の屈折率が一様でない場合は、必ずしも吸収率が一様ではない。
そのため、検出器上で形成するX線ビームの位置ずれから被検体の屈折量を求める際には、被検体で屈折したX線ビームの微小な変化量と被検体で減衰したX線の強度むらによる変化量を切り分けることができない。
By the way, not only the methods described in
In
However, when the refractive index in a specific space is not uniform, the absorptance is not necessarily uniform.
Therefore, when obtaining the amount of refraction of the subject from the positional deviation of the X-ray beam formed on the detector, the minute change amount of the X-ray beam refracted by the subject and the unevenness in the intensity of the X-ray attenuated by the subject. The amount of change due to cannot be separated.
本発明は、上記課題に鑑み、位相イメージングにより被検体の位相像を取得するに当たり、被検体における吸収情報に対する位相情報の割合が最大限取得可能となるX線撮像装置およびX線撮像方法の提供を目的とする。 In view of the above problems, the present invention provides an X-ray imaging apparatus and an X-ray imaging method capable of acquiring the maximum ratio of phase information to absorption information in a subject when acquiring a phase image of the subject by phase imaging. With the goal.
本発明のX線撮像装置は、
X線源と、
前記X線源から照射された発散X線を複数のX線に分割するための複数の開口を具備する、前記X線源と被検体との間に配置された分割素子と、
前記被検体を透過し屈折した前記発散X線によるX線ビームを検出するため、該X線ビームの照射方向の被検体よりも下流に配置された検出器と、
前記検出器により検出されたX線ビームの強度変化の情報を演算処理する演算処理手段と、を備え、
前記X線源と前記分割素子の距離をL1、前記分割素子と前記検出器の距離をL2とするとき、L2/L1の値が1.0より大きく2.4未満の範囲となるように構成されていることを特徴とする。
また、本発明のX線撮像方法は、X線源から照射された発散X線を分割素子により複数のX線に分割し、該分割された前記発散X線によるX線ビームを被検体に透過させた際に生じるX線ビームの強度変化を、検出器により検出するX線撮像方法であって、
前記分割素子と前記被検体と前記検出器とを、この順に前記X線源による前記発散X線の照射方向の光路上に配置し、
前記X線源と前記分割素子の距離と、前記分割素子と前記検出器の距離の関係を、前記X線ビームが検出器上で重ならない距離の関係とすることによって、
前記被検体における吸収情報に対する位相情報の割合を最大限取得可能とすることを特徴とする。
The X-ray imaging apparatus of the present invention
An X-ray source;
A splitting element disposed between the X-ray source and the subject, comprising a plurality of openings for splitting the divergent X-rays irradiated from the X-ray source into a plurality of X-rays;
A detector disposed downstream of the subject in the irradiation direction of the X-ray beam in order to detect an X-ray beam by the divergent X-ray transmitted through the subject and refracted;
An arithmetic processing means for arithmetically processing information on an intensity change of the X-ray beam detected by the detector,
When the distance between the X-ray source and the dividing element is L1, and the distance between the dividing element and the detector is L2, the value of L2 / L1 is in the range of more than 1.0 and less than 2.4. It is characterized by being.
In the X-ray imaging method of the present invention, divergent X-rays irradiated from an X-ray source are divided into a plurality of X-rays by a dividing element, and an X-ray beam by the divided divergent X-rays is transmitted to a subject. An X-ray imaging method for detecting a change in the intensity of an X-ray beam generated by a detector using a detector,
The splitting element, the subject, and the detector are arranged in this order on the optical path in the irradiation direction of the divergent X-ray by the X-ray source,
By making the relationship between the distance between the X-ray source and the splitting element and the distance between the splitting element and the detector into the relationship of the distance where the X-ray beam does not overlap on the detector,
It is possible to obtain the maximum ratio of phase information to absorption information in the subject.
本発明によれば、位相イメージングにより被検体の位相像を取得するに当たり、被検体における吸収情報に対する位相情報の割合が最大限取得可能となるX線撮像装置およびX線撮像方法を実現することができる。 According to the present invention, it is possible to realize an X-ray imaging apparatus and an X-ray imaging method capable of acquiring the maximum ratio of phase information to absorption information in a subject when acquiring a phase image of the subject by phase imaging. it can.
本発明は、X線源と分割素子と被検体と検出器とが、この順にX線源による発散X線の放射方向の光路上に配置し、X線源と分割素子の距離をL1とし、分割素子と検出器の距離をL2とするとき、
これらのL1とL2とを、X線ビームが検出器上で重ならない距離の関係とすることによって、被検体における吸収情報に対する位相情報の割合が最大限取得可能となることを見出したことによるものである。
このようなL2/L1の値は、X線源の焦点サイズf、分割素子の開口幅Ga、検出器の画素サイズPsから求めることができるが、これらの詳細については以下の本発明の実施形態の説明により明らかにする。
図1に、本実施形態における球面波を発生するX線源と発生したX線を空間的に分割する分割素子と二次元X線検出器を用いた、X線の位相変化から像を得るX線撮像装置の構成例を説明する図を示す。
図1に示すように、本実施形態のX線撮像装置では、発散X線102を発生させて被検体に対して照射するX線源101を備える。
そして、分割素子103、被検体104、検出器106を、この順にX線源101による発散X線102の照射方向の光路上に配置する。
発散X線102を空間的に分割する分割素子103は、シート状、ドット状の形状のX線ビームを形成する。これは、医療用レントゲン装置で、被検体を撮像中に発生する散乱X線を除去するために使用する集束型グリッドと同様に、X線を透過し易い軽元素で作られた箔112とX線を遮光する重元素で作られた箔113を交互に重ねて形成される。
X線を透過するための箔112はアルミニウム、紙、合成樹脂などが用いられる。
一方、X線を遮光するための箔113は白金、金、鉛、タンタル、タングステン等が用いられる。
また、X線を遮光するための箔113を一定の間隔を維持して両端を固定することで、X線を透過する部分を空気としても良く、本件では空気を含めて箔と呼ぶ。
In the present invention, the X-ray source, the dividing element, the subject, and the detector are arranged in this order on the optical path in the radiation direction of the divergent X-ray from the X-ray source, and the distance between the X-ray source and the dividing element is L1. When the distance between the dividing element and the detector is L2,
By finding that the ratio of the phase information to the absorption information in the subject can be obtained to the maximum by making these L1 and L2 a distance relationship in which the X-ray beams do not overlap on the detector. It is.
Such a value of L2 / L1 can be obtained from the focal point size f of the X-ray source, the aperture width Ga of the dividing element, and the pixel size Ps of the detector. Details of these will be described in the following embodiments of the present invention. It will be clarified by explanation of
FIG. 1 shows an X image obtained from an X-ray phase change using an X-ray source for generating a spherical wave, a dividing element for spatially dividing the generated X-ray, and a two-dimensional X-ray detector in this embodiment. The figure explaining the structural example of a linear imaging device is shown.
As shown in FIG. 1, the X-ray imaging apparatus of this embodiment includes an
Then, the dividing
The dividing
The
On the other hand, platinum, gold, lead, tantalum, tungsten, or the like is used for the
Further, by fixing the both ends of the
検出器106は、被検体を透過し屈折した発散X線によるX線ビームを検出するため、該X線ビームの照射方向において被検体よりも下流に配置されている。そして、検出器106は、複数の画素111が二次元的に配置されて構成されている。
図2は検出器106の一部分を図1におけるY方向から見た図である。
X線ビーム105の強度は、画素111エリア内に入射したX線のフォトン数を積算することで得ているため、検出器106の空間分解能は画素111の空間サイズで決まっている。
図2に示すように、被検体104を透過したX線ビーム105を複数の画素111に跨って照射することで、被検体が無い時の基準X線ビーム114に対する位置変化量を得ることが可能である。
一般的には、検出器106で形成されるX線ビーム105の幅が細くなるほど、各画素111で検出するX線強度に対するX線ビーム105の位置ずれで生じる強度差が相対的に向上する。
よって、X線ビーム105の位置ずれで生じる強度差がノイズに埋もれる可能性が低減し、被検体104の屈折量を明確に得ることが可能になるため位相検出感度が向上する。
本実施形態では、被検体を透過中に屈折したX線ビーム105と、被検体が無い時の基準X線ビーム114の検出器106上での位置ずれ量dxpと、X線ビーム105が検出器106上で形成する幅Gdを用いて位相検出感度をdxp/Gdで表す。
図2はX線ビーム105が隣接する2つの画素111に跨っているが、跨っている画素111の数は2以上跨っていれば本質的にいくつでも良い。
The
FIG. 2 is a view of a part of the
Since the intensity of the
As shown in FIG. 2, by irradiating the
In general, as the width of the
Therefore, the possibility that the intensity difference caused by the positional deviation of the
In the present embodiment, the
In FIG. 2, the
図2に示すように、X線ビーム105が検出器106上で形成する幅Gdは、X線源から発生する発散X線102の輝度は一定、分割素子103の開口幅をGaとし被検体104がない場合は、次の式1で表すことができる。
Gd=Ga×(L1+L2)/L1+f×L2/L1…式1
また、被検体を透過中に屈折したX線ビーム105と、被検体が無い時の基準X線ビーム114の検出器106上での位置ずれ量dxpは次の式2で表すことができる。
dxp=L2×tan(α)…式2
但し、αはX線ビーム105が被検体104で屈折する角度である。
被検体104を透過中に減衰したX線ビーム105は、被検体104の形状に応じて強度が変化するため、X線ビーム105の強度の重心位置と被検体が無い時の基準X線ビーム114の強度の重心位置にずれが生じる。
被検体104の吸収差で生じるX線ビーム105の位置ずれdxaは、被検体104の形状や厚みのみに依存する。
被検体104の屈折で生じるX線ビーム105の位置ずれdxpと吸収成分で生じるX線ビームの位置ずれdxaを完全に分離することは困難である。
しかし、dxaはL2に依存しないため、式2より、L2を大きくすることで吸収情報に対する位相情報の割合を高くすることが可能である。
また、一般的には、検出器106の画素111のサイズよりもdxpの方が小さいため、dxpを検出するためにもL2を大きく取ることが重要である。
As shown in FIG. 2, the width Gd formed on the
Gd = Ga × (L1 + L2) / L1 + f × L2 /
Further, the positional deviation amount dxp on the
dxp = L2 × tan (α) (2)
Here, α is an angle at which the
Since the intensity of the
The positional deviation dxa of the
It is difficult to completely separate the positional deviation dxp of the
However, since dxa does not depend on L2, it is possible to increase the ratio of the phase information to the absorption information by increasing L2 from
In general, since dxp is smaller than the size of the
しかしながら、有限の焦点サイズから発生する発散X線102を用いる場合は、L2を大きく取ると隣り合うX線ビーム105が検出器106上で重なるため、個々のX線ビーム105の正確な位置情報を得ることが難しくなる。
隣り合うX線ビーム105が検出器106上で重なり合わない条件は、一つのX線ビーム105の一次元方向の位置情報を検出するために必要な画素111の数nとしたとき次の式3で表すことができる。
Gd+dxp+dxa<n×Ps…式3
式3に式1を代入すると次の式4になる。
Ga×(L1+L2)/L1+f×L2/L1+dxp+dxa<n×Ps…式4
また、式4をL2/L1で括ると次の式5で表すことができる。
L2/L1+(dxp+dxa)/(Ga+f)<(n×Ps−Ga)/(Ga+f)…式5
Ga+fは100μm前後なのに対してdxp+dxaは最大で数ミクロン程度なためGa+f≫dxp+dxaとなり、式5は次の式6に変換することができる。
L2/L1<(n×Ps−Ga)/(Ga+f)…式6
図3に式1、式2、式6から、横軸に装置のジオメトリL2/L1、縦軸に位相検出感度dxp/Gdをプロットした図を示す。
図3の細線は、f=100μm、Ps=100μm、Ga=33μm、n=2の計算結果である。
式6よりL2/L1=1.25のとき位相検出感度dxp/Gdは最大値を取り、L2/L1>1.25では隣り合うX線ビーム105が検出器106上で重なるため、個々のX線ビーム105の正確な位置情報を得ることが難しくなる。
同様に、図3の太線はf=50μm、Ps=100μm、Ga=33μm、n=2の計算結果である。
式6よりL2/L1=2.0のとき位相検出感度dxp/Gdは最大値を取り、L2/L1>2.0では隣り合うX線ビーム105が検出器106上で重なり合う。
図3には示さないが、X線源101の焦点サイズを50〜150μm、画素111のサイズを50〜100μm、分割素子103の開口幅を20〜75μmとする。これにより、吸収情報に対する位相情報の割合を最大限得られる適切なL2/L1の範囲は1.0より大きく2.4未満(1.0<L2/L1<2.4)となる。
また、X線の輝度と位相検出感度を考慮した場合は、X線源の焦点サイズを80以上100μm以下とし、画素のサイズを50以上150μm以下とし、分割素子の開口の幅を20以上50μm以下とするとき、
前記X線源と前記分割素子の距離であるL1に対して、前記分割素子と前記検出器の距離であるL2が1.1〜1.5倍長とする(1.1<L2/L1<1.5とする)ことにより、吸収情報に対する位相情報の割合を最大限取得することが可能となる。
However, in the case of using
The condition that adjacent X-ray beams 105 do not overlap on the
Gd + dxp + dxa <n ×
Substituting
Ga × (L1 + L2) / L1 + f × L2 / L1 + dxp + dxa <n ×
Further, when
L2 / L1 + (dxp + dxa) / (Ga + f) <(n × Ps−Ga) / (Ga + f)
Ga + f is around 100 μm, whereas dxp + dxa is about several microns at the maximum, so Ga + f >> dxp + dxa, and
L2 / L1 <(n × Ps−Ga) / (Ga + f) (Formula 6)
FIG. 3 is a diagram in which the geometry L2 / L1 of the apparatus is plotted on the horizontal axis and the phase detection sensitivity dxp / Gd is plotted on the vertical axis from
The thin line in FIG. 3 is the calculation result when f = 100 μm, Ps = 100 μm, Ga = 33 μm, and n = 2.
From Equation 6, when L2 / L1 = 1.25, the phase detection sensitivity dxp / Gd takes the maximum value, and when L2 / L1> 1.25, the adjacent X-ray beams 105 overlap on the
Similarly, the thick lines in FIG. 3 are the calculation results when f = 50 μm, Ps = 100 μm, Ga = 33 μm, and n = 2.
From Equation 6, when L2 / L1 = 2.0, the phase detection sensitivity dxp / Gd takes the maximum value, and when L2 / L1> 2.0, adjacent X-ray beams 105 overlap on the
Although not shown in FIG. 3, the focus size of the
In consideration of X-ray luminance and phase detection sensitivity, the focal point size of the X-ray source is set to 80 to 100 μm, the pixel size is set to 50 to 150 μm, and the opening width of the dividing element is set to 20 to 50 μm. And when
L2 which is the distance between the dividing element and the detector is 1.1 to 1.5 times longer than L1 which is the distance between the X-ray source and the dividing element (1.1 <L2 / L1 < 1.5), it is possible to obtain the maximum ratio of the phase information to the absorption information.
本実施形態の分割素子103は一次元を元に明記しているが、二次元でも良い。
分割素子103が二次元の場合、X線ビーム105が被検体104で屈折する角度はx方向、y方向で独立している。
検出器106の画素111が二次元に配置されている場合は、y方向の屈折量に関しても図2で示したx方向の屈折量の検出方法と同様の方法で求めることが可能である。
The dividing
When the dividing
When the
以下に、本発明の実施例について説明する。
[実施例1]
実施例1においては、図4を用いて、本発明を適用したX線撮像装置の構成例について説明する。
図4において、101は発散X線102を発生するX線源、103は分割素子、104は被検体、106は二次元検出器である。
また、108、109、110はそれぞれX線源101、分割素子103、検出器106の移動・回転手段、107は検出器により検出されたX線ビームの強度変化の情報を演算処理する演算装置(演算処理手段)である。
本実施例において、X線源101としては、モリブデン、銀やタングステンターゲットの回転対陰極型のX線発生装置を用いる。
X線源101から発生した発散X線102の発生源における輝度の分布は正規分布を仮定し、焦点サイズは正規分布関数の半値全幅とする。X線源101から発生した発散X線102はX線源101から100cm離れた位置に配置した分割素子103により空間的に分割される。
Examples of the present invention will be described below.
[Example 1]
In the first embodiment, a configuration example of an X-ray imaging apparatus to which the present invention is applied will be described with reference to FIG.
In FIG. 4, 101 is an X-ray source for generating
In this embodiment, as the
The distribution of luminance at the source of
分割素子103は、X線を透過するための箔112とX線を遮光するための箔113を交互に重ね合わせ、集束位置は分割素子103からX線源方向100cmの位置にある。
この分割素子103にはD25、D33、D50、D75の4種類があり、それぞれX線を透過するための箔112の膜厚(Ga)とX線を遮光するための箔113の膜厚(Gb)が、次の表1に示すようにそれぞれ異なっている。
[表1]
The dividing
There are four types of
[Table 1]
被検体104は分割素子103の直後に設置し、分割素子103と検出器106の距離を10〜300cmとする。
図5は検出器106の一部分を図4におけるY方向から見た図である。
図5に示すように、被検体104を透過したX線ビーム105を複数の画素111に跨って照射することでX線ビーム105の検出器106上での位置及び屈折量を得ることが可能である。
本実施例の画素サイズは100μmである。
被検体を透過したX線ビーム105が照射されている場合の2つの画素のX線検出強度をI01、I02、被検体が無い時の基準X線ビーム114が照射されている2つの画素のX線検出強度をI11、I12とすると、
検出器106で検出可能なX線ビーム105の基準X線ビーム114に対する強度の変化量dIは、次の式7によって表される。
dI=(I01−I02)/(I01+I02)−(I11−I12)/(I11+I12)…式7
図6にD25の分割素子を用いてL1=1m,L2=1m,f=100μmの条件下で円柱状のアクリル樹脂を透過した個々のX線ビーム105の式7のdIをマッピングしたシミュレーション結果を示す。
本シミュレーションのX線源101は、17.5keV(0.071nm)の単色光で焦点サイズは100μmの正規分布関数を使用している。
また、検出器106のMTFの広がり関数には10lp/mmで0.1となる正規分布関数を使用している。
X線ビーム105が検出器106及び画素111にランダムに到達することによる空間的揺らぎで生じるポアソンノイズに代表されるノイズ系は含んでいない。 図6(a)の計算には、アクリル樹脂のパラメータに複素屈折率1−δ―iβを使用した。
また、図6(b)の計算には、アクリル樹脂のパラメータに複素屈折率1−δ―iβの実部を使用した。
また、図6(c)の計算には、アクリル樹脂のパラメータに複素屈折率1−δ―iβの虚部を使用した。
このことから、図6(a)の微分像にはアクリル樹脂の吸収成分と位相成分、図6(b)の微分像にはアクリル樹脂の位相成分、図6(c)の微分像にはアクリル樹脂の吸収成分をそれぞれ計算している。
つまり、図6(b)と図6(c)を比較することで、図6(a)に含まれる吸収情報に対する位相情報の割合が分かる。
The subject 104 is installed immediately after the dividing
FIG. 5 is a view of a part of the
As shown in FIG. 5, the position and the amount of refraction of the
The pixel size in this embodiment is 100 μm.
The X-ray detection intensities of the two pixels when the
An intensity change dI of the
dI = (I01−I02) / (I01 + I02) − (I11−I12) / (I11 + I12)...
FIG. 6 shows a simulation result obtained by mapping dI of Equation 7 of each
The
The MTF spread function of the
It does not include a noise system typified by Poisson noise caused by spatial fluctuation caused by the
In the calculation of FIG. 6B, the real part of the complex refractive index 1-δ-iβ was used as the parameter of the acrylic resin.
In the calculation of FIG. 6C, the imaginary part of the complex refractive index 1-δ-iβ was used as the parameter of the acrylic resin.
From this, the differential image of FIG. 6A is the absorption component and phase component of acrylic resin, the differential image of FIG. 6B is the phase component of acrylic resin, and the differential image of FIG. Each of the resin absorption components is calculated.
That is, by comparing FIG. 6B and FIG. 6C, the ratio of the phase information to the absorption information included in FIG.
図6(b)と図6(c)の破線部のラインプロファイルを図7の破線でそれぞれ示す。
図7(a)と図7(b)の破線のラインプロファイルから微分位相像の信号強度は、微分吸収像の信号強度よりも4倍程度大きいため、図6(a)に含まれる位相成分の割合は約80%であることが分かる。
同じシミュレーションをD25、D33、D50、D75の各分割素子103に対して、X線源101と分割素子103の距離を100cm、分割素子103と検出器106の距離を10〜300cmまで変化させた。そのときの微分像に含まれる位相成分の割合をプロットした図を図8に示す。
図8からD25の分割素子103ではL2=130cm、D33の分割素子103ではL2=120cm、D50の分割素子103ではL2=120cm、D75の分割素子103ではL2=100cmのとき、それぞれの位相成分の割合の最大値が得られる。
The line profiles of the broken line parts in FIG. 6B and FIG. 6C are indicated by broken lines in FIG.
Since the signal intensity of the differential phase image is about four times larger than the signal intensity of the differential absorption image based on the broken line profiles in FIGS. 7A and 7B, the phase component included in FIG. It can be seen that the percentage is about 80%.
In the same simulation, the distance between the
8 to D25, L2 = 130 cm,
吸収成分は被検体104の形状や厚みに依存するのに対して、位相成分は式2より分割素子103と検出器106の距離に依存する。
このことから、同じ分割素子103で同じ被検体104を撮像する場合、得られる微分像に含まれる位相成分の割合が最大になるL2/L1の値はL1及びL2の値によらず一定である。
図9に、D25の分割素子を用いて、X線源101と分割素子103までの距離を25cm、50cm、100cm、分割素子103と検出器106の距離を10〜300cmまで変化させた。そのときのアクリル樹脂の微分像に含まれる位相成分の割合をプロットした図を示す。
X線源101と分割素子103までの距離を増加すると微分像に含まれる位相成分の割合が増加するが、位相成分が最大になるL2/L1の値は常に1.3である。
同様の計算をD25の分割素子に対してf=90μmの条件下でシミュレーションを行うと位相成分が最大になるL2/L1の値は常に1.4、f=80μmの条件下で位相成分が最大になるL2/L1の値は常に1.5である。
The absorption component depends on the shape and thickness of the subject 104, whereas the phase component depends on the distance between the splitting
Therefore, when the
In FIG. 9, the distance between the
Increasing the distance between the
When the same calculation is performed for the D25 splitting element under the condition of f = 90 μm, the phase component becomes maximum L2 / L1 is always 1.4, and the phase component is maximum under the condition of f = 80 μm. The value of L2 / L1 is always 1.5.
本シミュレーションで得られる位相成分が最大になるL2/L1の値は、式6で得られる位相成分が最大になるL2/L1の値よりも1割程度小さい。
これは、式6を導出する過程でGa+f≫dxp+dxaの近似を用いていることと、検出器106のMTFによるX線ビーム105のボケを考慮していなことに起因している。
しかしながら、図9より、本シミュレーションで得られる位相成分が最大になるL2/L1と、式6で得られる位相成分が最大になるL2/L1のそれぞれの吸収成分に対する位相成分の割合は、
1〜2%程度の差しかないため、式6で得られるL2/L1を用いて被検体104を撮像しても良い。
The value of L2 / L1 that maximizes the phase component obtained in this simulation is about 10% smaller than the value of L2 / L1 that maximizes the phase component obtained by Equation 6.
This is because the approximation of Ga + f >> dxp + dxa is used in the process of deriving Equation 6, and the blur of the
However, from FIG. 9, the ratio of the phase component to the respective absorption components of L2 / L1 where the phase component obtained in this simulation is maximized and L2 / L1 where the phase component obtained by Equation 6 is maximized is
Since the difference is only about 1 to 2%, the subject 104 may be imaged using L2 / L1 obtained by Expression 6.
図9から100μmの焦点サイズのX線源101、D=25の分割素子103を用いて微分像に含まれる位相成分の割合が最大にするX線撮像装置の実際のジオメトリは、
L1=0.3mのときL2=0.4m、
L1=1mのときL2=1.3m、
L1=1.55mのときL2=2m、
L1=2.0mのときL2=2.6m、
などが選択肢として考えられる。
同様に、図8からD=50の分割素子103では、
L1=0.3mのときL2=0.36m、
L1=1.0mのときL2=1.2m、
L1=1.67mのときL2=2m、
L1=2mのときL2=2.4m、
などが選択肢として考えられる。
図9より、L2/L1が一定下では、L1及びL2が長くなるほど微分像に含まれる位相成分の割合が増加するが、L1+L2は装置の大きさに直接影響を及ぼしたり、X線ビーム105一本当たりのX線強度が減少したりする課題が残っている。
The actual geometry of the X-ray imaging apparatus that maximizes the ratio of the phase component included in the differential image using the
L2 = 0.4 m when L1 = 0.3 m,
L2 = 1.3 m when L1 = 1 m,
L2 = 2m when L1 = 1.55m,
L2 = 2.6 m when L1 = 2.0 m,
Etc. are considered as options.
Similarly, from FIG. 8, in the dividing
L2 = 0.36m when L1 = 0.3m,
When L1 = 1.0m, L2 = 1.2m,
L2 = 2m when L1 = 1.67m,
L2 = 2.4m when L1 = 2m,
Etc. are considered as options.
From FIG. 9, when L2 / L1 is constant, the ratio of the phase component included in the differential image increases as L1 and L2 become longer. However, L1 + L2 directly affects the size of the apparatus, and the
[実施例2]
実施例2においては、実施例1を説明した一次元方向のX線位置変化検出に一次元加えて、二次元方向のX線位置変化検出を同時に測定する方法について説明する。
図10を用いて、X線撮像装置の構成例について説明する。図10(a)において、101は発散X線を発生するX線源、103と201は一次元分割素子、104は被検体、106は二次元検出器である。
108、109、117、202はそれぞれX線源101、一枚目の分割素子103、二枚目の分割素子201、検出器106の移動・回転手段である。
[Example 2]
In the second embodiment, a method of simultaneously measuring two-dimensional X-ray position change detection in addition to the one-dimensional X-ray position change detection described in the first embodiment will be described.
A configuration example of the X-ray imaging apparatus will be described with reference to FIG. In FIG. 10A, 101 is an X-ray source that generates divergent X-rays, 103 and 201 are one-dimensional division elements, 104 is a subject, and 106 is a two-dimensional detector.
X線源101から発生したX線は、一枚目の分割素子103のX線が透過する箔112と二枚目の分割素子201のX線が透過する箔112が空間的に重なっている部分のみ透過することができる。このため、分割素子103、201で分割したX線ビーム105は二次元ドットアレイ状を形成している。
一般的な検出器106における二次元情報は、x方向、y方向で90度ずれているため、一枚目の分割素子103と二枚目の分割素子201の箔が積層されている方向が、互いに90度ずれているとき最も都合がよい。
検出器106における二次元情報のx方向、y方向の角度が90度以外の場合は、その角度に合わせて一枚目の分割素子103と二枚目の分割素子201の角度を調整することが望ましい。
また、検出器106のx方向、y方向の角度と、一枚目の分割素子103と二枚目の分割素子201の箔が積層されている方向の角度が異なっていても、式7を用いてdIのマッピングを作成する際に補正を加えることで屈折情報を明確にすることができる。
The X-ray generated from the
Since the two-dimensional information in the
When the angle in the x direction and y direction of the two-dimensional information in the
Even if the angles of the
X線ビーム105が被検体104で屈折する角度は、x方向とy方向で独立しているため、被検体104の位相成分の割合を最大にするL2/L1の比は、二次元の場合でも一次元と同じである。
このことから、図8及び図9よりD25の一次元分割素子103とD25の一次元分割素子201の2枚用いて二次元状の分割素子とした場合では、
L2/L1=1.3のときx方向y方向共に位相成分の割合が最大値を得ることができる。
また、図8より1枚目の分割素子103と2枚目の分割素子201の距離は10cm程度離れていても、x方向とy方向の微分像に含まれる位相成分の割合の誤差は数%程度ですむため画質に大きな問題を生じにくい。
図9に示すように、1枚目の分割素子103と2枚目の分割素子201の位置関係で生じる微分像に含まれる位相成分の割合の誤差は、X線源101と検出器106の距離が長くなるほど少なくなる。
Since the angle at which the
From this, in the case where a two-dimensional dividing element is used by using two pieces of the one-
When L2 / L1 = 1.3, the maximum ratio of the phase components can be obtained in both the x and y directions.
Further, from FIG. 8, even if the distance between the
As shown in FIG. 9, the error of the ratio of the phase component included in the differential image generated by the positional relationship between the
二次元状の分割素子は、図10(b)に示すように二枚の一次元分割素子103を互いに面内で90度維持した状態で貼り合わせて形成されたものでもよい。また、X線を遮光するタングステンなどの板に放電加工で穴をアレイ状に空けたものや、X線の透過するシリコン基板上にX線を遮光する物質を微細加工にてアレイ状に形成されたものでもよい。
As shown in FIG. 10B, the two-dimensional dividing element may be formed by bonding two one-
101:X線源
102:発散X線
103:分割素子
104:被検体
105:X線ビーム
106:検出器
107:演算装置
108:X線源を移動・回転する手段
109:分割素子を移動・回転する手段
110:二次元検出器を移動・回転する手段
111:画素
112:X線を透過するための箔
113:X線を遮光するための箔
114:被検体が無い時の基準X線ビーム
101: X-ray source 102: Divergent X-ray 103: Splitting element 104: Subject 105: X-ray beam 106: Detector 107: Computing device 108: Means for moving / rotating the X-ray source 109: Moving / rotating the splitting element Means 110: means for moving / rotating the two-dimensional detector 111: pixel 112: foil for transmitting X-ray 113: foil for shielding X-ray 114: reference X-ray beam when there is no subject
Claims (6)
X線源と、
前記X線源から照射された発散X線を複数のX線に分割するための複数の開口を具備する、前記X線源と被検体との間に配置された分割素子と、
前記被検体を透過し屈折した前記発散X線によるX線ビームを検出するため、該X線ビームの照射方向の被検体よりも下流に配置された検出器と、
前記検出器により検出されたX線ビームの強度変化の情報を演算処理する演算処理手段と、を備え、
前記X線源と前記分割素子の距離をL1、前記分割素子と前記検出器の距離をL2とするとき、L2/L1の値が1.0より大きく2.4未満の範囲となるように構成されていることを特徴とするX線撮像装置。 An X-ray imaging apparatus,
An X-ray source;
A splitting element disposed between the X-ray source and the subject, comprising a plurality of openings for splitting the divergent X-rays irradiated from the X-ray source into a plurality of X-rays;
A detector disposed downstream of the subject in the irradiation direction of the X-ray beam in order to detect an X-ray beam by the divergent X-ray transmitted through the subject and refracted;
An arithmetic processing means for arithmetically processing information on an intensity change of the X-ray beam detected by the detector,
When the distance between the X-ray source and the dividing element is L1, and the distance between the dividing element and the detector is L2, the value of L2 / L1 is in the range of more than 1.0 and less than 2.4. An X-ray imaging apparatus characterized by that.
前記画素のサイズを50以上150μm以下とし、
前記分割素子の前記開口の幅を20以上50μm以下とするとき、
前記L2/L1の値が1.1より大きく1.5未満の範囲となるように構成されていることを特徴とする請求項1に記載のX線撮像装置。 The focal spot size of the X-ray source is 80 to 100 μm,
The size of the pixel is 50 to 150 μm,
When the width of the opening of the dividing element is 20 to 50 μm,
The X-ray imaging apparatus according to claim 1, wherein the value of L2 / L1 is configured to be in a range greater than 1.1 and less than 1.5.
0.3m≦L1≦2m、0.3m≦L2≦4.8m
The distance L1 which is the distance between the X-ray source and the dividing element and the distance L2 which is the distance between the dividing element and the detector satisfy the following ranges, respectively. The X-ray imaging apparatus according to Item.
0.3m ≦ L1 ≦ 2m, 0.3m ≦ L2 ≦ 4.8m
前記分割素子と前記被検体と前記検出器とを、この順に前記X線源による前記発散X線の照射方向の光路上に配置し、
前記X線源と前記分割素子の距離と、前記分割素子と前記検出器の距離の関係を、前記X線ビームが検出器上で重ならない距離の関係とすることによって、
前記被検体における吸収情報に対する位相情報の割合を最大限取得可能とすることを特徴とするX線撮像方法。 Diversity X-rays irradiated from an X-ray source are divided into a plurality of X-rays by a dividing element, and the intensity change of the X-ray beam generated when the X-ray beam generated by the divided X-rays is transmitted through the subject. Is an X-ray imaging method for detecting
The splitting element, the subject, and the detector are arranged in this order on the optical path in the irradiation direction of the divergent X-ray by the X-ray source,
By making the relationship between the distance between the X-ray source and the splitting element and the distance between the splitting element and the detector into the relationship of the distance where the X-ray beam does not overlap on the detector,
An X-ray imaging method characterized in that a maximum ratio of phase information to absorption information in the subject can be acquired.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011036114A JP2012170645A (en) | 2011-02-22 | 2011-02-22 | X-ray imaging apparatus, and x-ray imaging method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011036114A JP2012170645A (en) | 2011-02-22 | 2011-02-22 | X-ray imaging apparatus, and x-ray imaging method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2012170645A true JP2012170645A (en) | 2012-09-10 |
| JP2012170645A5 JP2012170645A5 (en) | 2014-04-03 |
Family
ID=46974114
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2011036114A Pending JP2012170645A (en) | 2011-02-22 | 2011-02-22 | X-ray imaging apparatus, and x-ray imaging method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2012170645A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2016077904A (en) * | 2014-10-17 | 2016-05-16 | キヤノン株式会社 | Imaging method, image processing apparatus, computer-readable medium, method, apparatus and system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH04241841A (en) * | 1991-01-10 | 1992-08-28 | Fuji Photo Film Co Ltd | Radiograph photographing device |
| US5812629A (en) * | 1997-04-30 | 1998-09-22 | Clauser; John F. | Ultrahigh resolution interferometric x-ray imaging |
| JP2000512764A (en) * | 1997-01-24 | 2000-09-26 | クウォンタ・ビジョン・インコーポレイテッド | Inspection apparatus using small-angle topographic method for determining internal structure and composition of object |
| WO2007125833A1 (en) * | 2006-04-24 | 2007-11-08 | The University Of Tokyo | X-ray image picking-up device and x-ray image picking-up method |
-
2011
- 2011-02-22 JP JP2011036114A patent/JP2012170645A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH04241841A (en) * | 1991-01-10 | 1992-08-28 | Fuji Photo Film Co Ltd | Radiograph photographing device |
| JP2000512764A (en) * | 1997-01-24 | 2000-09-26 | クウォンタ・ビジョン・インコーポレイテッド | Inspection apparatus using small-angle topographic method for determining internal structure and composition of object |
| US5812629A (en) * | 1997-04-30 | 1998-09-22 | Clauser; John F. | Ultrahigh resolution interferometric x-ray imaging |
| WO2007125833A1 (en) * | 2006-04-24 | 2007-11-08 | The University Of Tokyo | X-ray image picking-up device and x-ray image picking-up method |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2016077904A (en) * | 2014-10-17 | 2016-05-16 | キヤノン株式会社 | Imaging method, image processing apparatus, computer-readable medium, method, apparatus and system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5269041B2 (en) | X-ray imaging apparatus and X-ray imaging method | |
| US9234856B2 (en) | X-ray apparatus and X-ray measuring method | |
| US20190250113A1 (en) | X-ray interrogation system using multiple x-ray beams | |
| US9031189B2 (en) | X-ray imaging apparatus and X-ray imaging method | |
| US20150055745A1 (en) | Phase Contrast Imaging Using Patterned Illumination/Detector and Phase Mask | |
| KR101331542B1 (en) | X-ray imaging apparatus and x-ray imaging method | |
| EP2557437A2 (en) | Beam Sensing | |
| JP2012090945A (en) | Radiation detection device, radiographic apparatus, and radiographic system | |
| JP2011224329A (en) | Radiation imaging system and method | |
| CN106706157A (en) | Quasi-concentric visual axis-based ICF (inertial confinement fusion) hot spot electronic temperature detection device | |
| US9239304B2 (en) | X-ray imaging apparatus | |
| WO2017045657A1 (en) | A method of detection of defects in materials with internal directional structure and a device for performance of the method | |
| JP2012170645A (en) | X-ray imaging apparatus, and x-ray imaging method | |
| US20150153290A1 (en) | X-ray apparatus and method of measuring x-rays | |
| US11000249B2 (en) | X-ray detector for grating-based phase-contrast imaging | |
| JP2008256587A (en) | X-ray inspection apparatus and X-ray inspection method | |
| JP2014012132A (en) | X-ray apparatus and method for measuring x-ray | |
| JP5733908B2 (en) | X-ray imaging device | |
| JP5610885B2 (en) | X-ray imaging apparatus and imaging method | |
| JP2013150780A (en) | Radiographic apparatus and radiographic method | |
| JP2014014670A (en) | X-ray apparatus and x-ray measurement method | |
| JP2017044588A (en) | X-ray detector and X-ray imaging system | |
| JP2023166197A (en) | X-ray detection device and X-ray detection method | |
| JP5930614B2 (en) | X-ray imaging device | |
| Tkaczyk et al. | Response of a CsI/amorphous-Si flat panel detector as function of incident x-ray angle |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20131212 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140214 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20140214 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20141017 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20141202 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150130 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150609 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20151020 |