JP2013123460A - 放射線撮影装置及び放射線撮影方法 - Google Patents
放射線撮影装置及び放射線撮影方法 Download PDFInfo
- Publication number
- JP2013123460A JP2013123460A JP2011272142A JP2011272142A JP2013123460A JP 2013123460 A JP2013123460 A JP 2013123460A JP 2011272142 A JP2011272142 A JP 2011272142A JP 2011272142 A JP2011272142 A JP 2011272142A JP 2013123460 A JP2013123460 A JP 2013123460A
- Authority
- JP
- Japan
- Prior art keywords
- image
- grating
- group
- pixel
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title description 27
- 238000003384 imaging method Methods 0.000 claims abstract description 52
- 230000005855 radiation Effects 0.000 claims abstract description 44
- 230000010363 phase shift Effects 0.000 claims abstract description 30
- 230000000694 effects Effects 0.000 claims abstract description 12
- 230000007547 defect Effects 0.000 claims description 52
- 230000002159 abnormal effect Effects 0.000 claims description 32
- 230000000737 periodic effect Effects 0.000 claims description 25
- LFEUVBZXUFMACD-UHFFFAOYSA-H lead(2+);trioxido(oxo)-$l^{5}-arsane Chemical compound [Pb+2].[Pb+2].[Pb+2].[O-][As]([O-])([O-])=O.[O-][As]([O-])([O-])=O LFEUVBZXUFMACD-UHFFFAOYSA-H 0.000 claims description 20
- 238000003860 storage Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 abstract description 5
- 230000008859 change Effects 0.000 description 18
- 239000011295 pitch Substances 0.000 description 16
- 238000010521 absorption reaction Methods 0.000 description 15
- 238000012937 correction Methods 0.000 description 10
- 230000005540 biological transmission Effects 0.000 description 8
- 238000009826 distribution Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 230000005856 abnormality Effects 0.000 description 4
- 238000007689 inspection Methods 0.000 description 4
- 239000004065 semiconductor Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 239000010408 film Substances 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 3
- MCVAAHQLXUXWLC-UHFFFAOYSA-N [O-2].[O-2].[S-2].[Gd+3].[Gd+3] Chemical compound [O-2].[O-2].[S-2].[Gd+3].[Gd+3] MCVAAHQLXUXWLC-UHFFFAOYSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 2
- XQPRBTXUXXVTKB-UHFFFAOYSA-M caesium iodide Chemical compound [I-].[Cs+] XQPRBTXUXXVTKB-UHFFFAOYSA-M 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000002601 radiography Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- BUGBHKTXTAQXES-UHFFFAOYSA-N Selenium Chemical compound [Se] BUGBHKTXTAQXES-UHFFFAOYSA-N 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 210000000845 cartilage Anatomy 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 229910052697 platinum Inorganic materials 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 229910052711 selenium Inorganic materials 0.000 description 1
- 239000011669 selenium Substances 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 239000007779 soft material Substances 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 1
- 229910052721 tungsten Inorganic materials 0.000 description 1
- 239000010937 tungsten Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/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/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
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21K—TECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
- G21K2207/00—Particular details of imaging devices or methods using ionizing electromagnetic radiation such as X-rays or gamma rays
- G21K2207/005—Methods and devices obtaining contrast from non-absorbing interaction of the radiation with matter, e.g. phase contrast
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Measurement Of Radiation (AREA)
Abstract
【課題】画像データに生じた線状欠陥の影響を低減し、精度よく位相微分画像を生成可能とする。
【解決手段】X線撮影装置は、X線源から放射された放射線を通過させて第1の周期パターン像を生成する第1の格子と、第1の周期パターン像を部分的に遮蔽してモアレ縞が生じた第2の周期パターン像を生成する第2の格子と、行方向とこれに直交する列方向に沿って2次元配列された複数の画素部により第2の周期パターン像を検出して画像データを生成するX線画像検出器と、画像データに画像処理を施す画像処理部とを備える。行方向及び列方向のうち、線状欠陥が少ない方向に並ぶ所定数の画素部を1グループGr(x,n)とし、このグループGr(x,n)を該方向に所定画素ずつ変更しながら順に設定する。各グループGr(x,n)内に含まれる画素値により構成された強度変調信号の位相ズレ量を算出することにより位相微分画像を生成する。
【選択図】図7
【解決手段】X線撮影装置は、X線源から放射された放射線を通過させて第1の周期パターン像を生成する第1の格子と、第1の周期パターン像を部分的に遮蔽してモアレ縞が生じた第2の周期パターン像を生成する第2の格子と、行方向とこれに直交する列方向に沿って2次元配列された複数の画素部により第2の周期パターン像を検出して画像データを生成するX線画像検出器と、画像データに画像処理を施す画像処理部とを備える。行方向及び列方向のうち、線状欠陥が少ない方向に並ぶ所定数の画素部を1グループGr(x,n)とし、このグループGr(x,n)を該方向に所定画素ずつ変更しながら順に設定する。各グループGr(x,n)内に含まれる画素値により構成された強度変調信号の位相ズレ量を算出することにより位相微分画像を生成する。
【選択図】図7
Description
本発明は、被検体による放射線の位相変化に基づく画像を検出する放射線撮影装置及び放射線撮影方法に関する。
放射線、例えば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線位相イメージングの一種として、2枚の格子とX線画像検出器とを用いてX線タルボ干渉計を構成することにより、X線の位相変化を検出するX線撮影装置が知られている(例えば、特許文献1、2参照)。
このX線撮影装置は、X線源から見て被検体の背後に第1の格子を配置し、第1の格子からタルボ距離だけ離れた位置に第2の格子を配置し、その背後にX線画像検出器を配置したものである。タルボ距離は、第1の格子を通過したX線が、タルボ効果によって第1の格子の自己像(縞画像)を形成する距離である。この自己像は、被検体でのX線の位相変化で屈折が生じることにより変調される。この変調量を検出することにより、X線の位相変化が画像化される。
上記変調量の検出方法として縞走査法が知られている。縞走査法とは、第1の格子に対して第2の格子を、所定のピッチずつ間欠的に移動させながら、その各停止中にX線源からX線を放射し、被検体、第1及び第2の格子を通過したX線を、X線画像検出器を構成する複数の画素部により検出する方法である。第2の格子の移動方向は、第1の格子の面に平行で、かつ第1の格子の格子方向に垂直な方向である。
縞走査法では、X線画像検出器の各画素部について、第2の格子の移動に対する画素値の強度変化を表す強度変調信号を生成し、この強度変調信号の位相ズレ量(被検体が存在しない場合の初期位置からの位相差)を算出することにより、上記変調量に関連する画像が得られる。この画像は、被検体の屈折率を反映した画像であり、X線の位相変化(位相シフト)の微分量に対応するため、位相微分画像と呼ばれる。この位相微分画像を積分処理することにより、位相シフト分布を表す位相コントラスト画像が生成される。
しかし、縞走査法では、1枚の位相コントラスト画像を取得するために複数回(例えば、5回)の撮影を行う必要があるため、装置の振動や被検体の体動により画質が低下するという問題がある。また、複数回の撮影により被検体の被曝量が大きいことも問題である。特許文献1には、第1及び第2の格子を移動させずに一度の撮影で得られる単一の画像データから位相微分画像を生成することが記載されているが、その具体的な方法については述べられていない。
そこで、本出願人は、特願2010−265241号において、第1及び第2の格子を移動させずに一度の撮影で得られる単一の画像データから位相微分画像を生成する具体的方法を提案している。この方法では、第1の格子に対して第2の格子22を格子面内方向に相対的に僅かに傾斜させて、モアレ縞を発生させた状態で、X線画像検出器により単一の画像データを取得する。そして、モアレ縞に交わる方向に並ぶ所定数の画素部を1グループとし、このグループを変更しながら、各グループに含まれる画素部の画素値により強度変調信号を構成し、この強度変調信号の位相ズレ量を算出することにより位相微分画像が生成される。
しかしながら、特許文献2に記載のように、X線画像検出器により得られる画像データには、異常画素値が線状に連続した線状欠陥が生じることがある。特願2010−265241号に記載の方法では、強度変調信号を構成するためのグループが一方向に延伸しているため、線状欠陥の方向とグループの延伸方向とが一致した場合には、グループを構成する画素値がすべて異常画素値となり、正常な強度変調信号が得られない。この結果、位相微分画像には、画像データの線状欠陥に対応した線状欠陥が生じてしまう。
このように位相微分画像に線状欠陥が生じた場合には、これに積分処理を施すことにより得られる位相コントラスト画像には、線状欠陥の一部を積分領域に含む画素すべてに欠陥の影響が生じるため、広範囲に画質低下が生じてしまう。
本発明は、画像データに生じた線状欠陥の影響を低減し、精度よく位相微分画像を生成可能とする放射線撮影装置及び放射線撮影方法を提供することを目的とする。
上記目的を達成するために、本発明の放射線撮影装置は、放射線源から放射された放射線を通過させて第1の周期パターン像を生成する第1の格子と、前記第1の周期パターン像を部分的に遮蔽してモアレ縞が生じた第2の周期パターン像を生成する第2の格子と、行方向とこれに直交する列方向に沿って2次元配列された複数の画素部により、前記第2の周期パターン像を検出して画像データを生成する放射線画像検出器と、前記行方向及び前記列方向のうち、線状欠陥が少ない第1の方向に並ぶ所定数の前記画素部を1グループとし、このグループを前記第1の方向に所定画素ずつ変更しながら順に設定するグループ設定部と、前記各グループ内に含まれる画素値により構成された強度変調信号の位相ズレ量を算出することにより位相微分画像を生成する位相微分画像生成部と、を備えるものである。
前記モアレ縞は、前記第2の格子を、第1の格子に対して格子面内方向に相対的に傾斜して配置することにより生成されたものであり、前記第1及び第2の格子の格子方向にほぼ直交していることが好ましい。
前記モアレ縞は、前記第1及び第2の格子の対向方向の位置関係、または、前記第1及び第2の格子の格子ピッチを調整することにより生成されたものであり、前記第1及び第2の格子の格子方向にほぼ平行であってもよい。
前記第1の方向は、前記モアレ縞にほぼ直交していることが好ましい。
前記モアレ縞は、前記第2の格子を、第1の格子に対して格子面内方向に相対的に傾斜して配置し、かつ、前記第1及び第2の格子の対向方向の位置関係、または、前記第1及び第2の格子の格子ピッチを調整することにより生成されたものであり、前記第1及び第2の格子の格子方向に直交せず、かつ平行でなくてもよい。
前記放射線画像検出器の画素欠陥のマップ情報を記憶する画素欠陥情報記憶部と、前記マップ情報に基づき、前記各グループ内の画素値から異常画素値を除去する異常画素値除去部と、を備え、前記位相微分画像生成部は、前記各グループ内に含まれる前記異常画素値以外の正常画素値に基づいて前記強度変調信号の位相ズレ量を算出し、前記位相微分画像を生成することが好ましい。
前記各グループ内の前記正常画素値の数をカウントする正常画素値計数部を備えることが好ましい。この場合、前記正常画素値の数が基準値より少ないグループが存在する場合にエラー報知を行うエラー報知部を備えることが好ましい。
前記グループを構成する前記画素部の数は、前記モアレ縞の周期の整数倍に含まれる前記画素部の数に等しいことが好ましい。
前記位相微分画像を積分処理して位相コントラスト画像を生成する位相コントラスト画像生成部を備えることが好ましい。
前記第1の格子は、吸収型格子であり、入射した放射線を幾何光学的に投影することにより前記第1の周期パターン像を生成することが好ましい。
前記第1の格子は、吸収型格子または位相型格子であり、入射した放射線にタルボ効果を生じさせて前記第1の周期パターン像を生成するものであってもよい。
前記放射線源から放射された放射線を部分的に遮蔽して焦点を分散化するマルチスリットを備えることが好ましい。
本発明の放射線撮影方法は、放射線源から放射された放射線を第1の格子を通過させて第1の周期パターン像を生成するステップと、前記第1の周期パターン像を第2の格子により部分的に遮蔽してモアレ縞が生じた第2の周期パターン像を生成するステップと、行方向とこれに直交する列方向に沿って2次元配列された複数の画素部を有する放射線画像検出器により前記第2の周期パターン像を検出して画像データを生成するステップと、前記行方向及び前記列方向のうち、線状欠陥が少ない第1の方向に並ぶ所定数の前記画素部を1グループとし、このグループを前記第1の方向に所定画素ずつ変更しながら順に設定するステップと、前記各グループ内に含まれる画素値により構成された強度変調信号の位相ズレ量を算出することにより位相微分画像を生成するステップと、を備える。
本発明によれば、放射線画像検出器の行方向及び前記列方向のうち、線状欠陥が少ない第1の方向に並ぶ所定数の前記画素部を1グループとし、このグループを前記第1の方向に所定画素ずつ変更しながら順に設定するので、画像データに生じた線状欠陥の影響を低減し、精度よく位相微分画像を生成することができる。
(第1実施形態)
図1において、X線撮影装置10は、X線源11、格子部12、X線画像検出器13、メモリ14、画像処理部15、画像記録部16、撮影制御部17、コンソール18、及びシステム制御部19を備える。X線源11は、周知のように、回転陽極型のX線管(図示せず)と、X線の照射野を制限するコリメータ(図示せず)とを有し、撮影制御部17の制御に基づき、被検体Hに向けてX線を放射する。
図1において、X線撮影装置10は、X線源11、格子部12、X線画像検出器13、メモリ14、画像処理部15、画像記録部16、撮影制御部17、コンソール18、及びシステム制御部19を備える。X線源11は、周知のように、回転陽極型のX線管(図示せず)と、X線の照射野を制限するコリメータ(図示せず)とを有し、撮影制御部17の制御に基づき、被検体Hに向けてX線を放射する。
格子部12は、第1の格子21及び第2の格子22を備える。第1及び第2の格子21,22は、X線照射方向であるZ方向に関してX線源11に対向配置されている。X線源11と第1の格子21との間には、被検体Hが配置可能な間隔が設けられている。X線画像検出器13は、周知のように、半導体回路を用いたフラットパネル検出器であり、第2の格子22の背後に近接して配置されている。X線画像検出器13の検出面13aは、Z方向に直交している。
第1の格子21は、Z方向に直交する格子面(XY面)におけるY方向に延伸された複数のX線吸収部21a及びX線透過部21bを備えた吸収型格子である。X線吸収部21a及びX線透過部21bは、X方向(Z方向及びY方向に直交する方向)に交互に配列されており、縞状のパターンを形成している。第2の格子22は、第1の格子21と同様にY方向に延伸され、かつX方向に交互に配列された複数のX線吸収部22a及びX線透過部22bを備えた吸収型格子である。X線吸収部21a,22aは、金(Au)、白金(Pt)等のX線吸収性を有する材料により形成されている。X線透過部21b,22bは、シリコン(Si)や樹脂等のX線透過性を有する材料や空隙により形成されている。
第1の格子21は、X線源11から射出されたX線を部分的に通過させて第1の周期パターン像(以下、G1像という)を生成する。第2の格子22は、第1の格子21により生成されたG1像を部分的に透過させて第2の周期パターン像(以下、G2像という)を生成する。被検体Hが配置されていない場合において、G1像は、第2の格子22の格子パターンとほぼ一致する。詳しくは後述するが、第2の格子22に対して第1の格子21がZ軸周り(格子面内方向)に僅かに傾斜しており、G2像には、その傾斜角に応じた周期を有するモアレ縞が生じている。
X線画像検出器13は、G2像を検出して画像データを生成する。メモリ14は、X線画像検出器13から読み出された画像データを一時的に記憶する。画像処理部15は、メモリ14に記憶された画像データに基づいて位相微分画像を生成し、この位相微分画像に基づいて位相コントラスト画像を生成する。画像記録部16は、位相微分画像と位相コントラスト画像とを記録する。
コンソール18は、操作部18a及びモニタ18bを備えている。操作部18aは、キーボードやマウス等により構成され、X線源11の管電圧、管電流、照射時間等の撮影条件の設定や、撮影実行指示等の操作入力を可能とする。
モニタ18bは、撮影条件等の撮影情報や、画像記録部16に記録された位相微分画像及び位相コントラスト画像の表示を行う。システム制御部19は、操作部18aから入力される信号に応じて各部を統括的に制御する。
図2において、X線画像検出器13は、周知のように、入射X線により半導体膜(図示せず)に生じた電荷を収集する画素電極31と、画素電極31によって収集された電荷を読み出すためのTFT(Thin Film Transistor)32とを備えた画素部30が2次元格子状に多数配列されたものである。半導体膜は、例えば、アモルファスセレンにより形成されている。
また、X線画像検出器13は、ゲート走査線33、走査回路34、信号線35、及び読み出し回路36を備える。ゲート走査線33は、画素部30の行ごとに設けられている。走査回路34は、TFT32をオン/オフするための走査信号を各ゲート走査線33に付与する。信号線35は、画素部30の列ごとに設けられている。読み出し回路36は、各信号線35を介して画素部30から電荷を読み出し、画像データに変換して出力する。各画素部30の詳細な層構成については、例えば、特開2002−26300号公報に記載された層構成と同様である。
読み出し回路36は、積分アンプ、A/D変換器、補正回路(いずれも図示せず)等を備える。積分アンプは、各画素部30から信号線35を介して出力された電荷を積分して画像信号を生成する。A/D変換器は、積分アンプにより生成された画像信号を、デジタル形式の画像データに変換する。補正回路は、画像データに対して、暗電流補正、ゲイン補正、リニアリティ補正等を行う。この補正後の画像データがメモリ14に記憶される。
X線画像検出器13は、入射X線を半導体膜で直接電荷に変換する直接変換型に限られず、ヨウ化セシウム(CsI)やガドリウムオキシサルファイド(GOS)等のシンチレータで入射X線を可視光に変換し、可視光をフォトダイオードで電荷に変換する間接変換型であってもよい。さらに、X線画像検出器13を、シンチレータとCMOSセンサ等の固体撮像素子とを組み合わせて構成してもよい。
X線画像検出器13により生成される画像データには、信号線35やゲート走査線33の異常により、信号線35に沿う方向(以下、列方向という)やゲート走査線33に沿う方向(以下、行方向という)に複数の異常画素値が並んだ線状欠陥が生じることがある。この線状欠陥を含む画素欠陥のマップ情報は、X線画像検出器13の製造最終段階の検査工程において取得される。
X線画像検出器13は、第1及び第2の格子21,22の格子方向(Y方向)に線状欠陥の方向が一致しないように配置されている。例えば、X線画像検出器13の信号線35に異常が生じ、列方向に沿って線状欠陥が生じた場合には、行方向がY方向に沿うようにX線画像検出器13を配置する。なお、列方向と行方向との双方に線状欠陥が存在する場合には、列方向と行方向とのうち線状欠陥が少ない方向がY方向に沿うようにX線画像検出器13を配置する。
図3において、X線源11から放射されるX線は、X線焦点11aを発光点としたコーンビームである。第1の格子21は、タルボ効果が生じず、X線透過部21bを通過したX線をほぼ幾何光学的に投影するように構成される。具体的には、X方向へのX線透過部21bの幅を、X線源11から照射されるX線の実効波長より十分大きな値とし、X線の大部分がX線透過部21bで回折しないようにすることで実現される。X線源11の回転陽極としてタングステンを用い、管電圧を50kVとした場合には、X線の実効波長は約0.4Åである。この場合には、X線透過部21bの幅を1〜10μm程度とすればよい。
これにより、G1像は、第1の格子21からZ方向下流への距離に依らず、常に第1の格子21の自己像となる。G1像は、X線焦点11aからZ方向下流への距離に比例して拡大される。
第2の格子22の格子ピッチp2は、前述のように、第2の格子22の格子パターンが第2の格子22の位置におけるG1像に一致するように設定されている。具体的には、第2の格子22の格子ピッチp2は、第1の格子21の格子ピッチp1、X線焦点11aと第1の格子21との間の距離L1、第1の格子21と第2の格子22との間の距離L2と、下式(1)をほぼ満たすように設定されている。
G1像は、被検体HでX線に位相変化が生じて屈折することにより変調される。この変調量には、被検体HでのX線の屈折角φ(x)が反映される。同図には、被検体HでのX線の位相変化を表す位相シフト分布Φ(x)に応じて屈折するX線の経路が例示されている。符号X1は、被検体Hが存在しない場合にX線が直進する経路を示し、符号X2は、被検体Hにより屈折したX線の経路を示している。
位相シフト分布Φ(x)は、X線の波長をλ、被検体Hの屈折率分布をn(x,z)として、下式(2)で表される。
屈折角φ(x)は、位相シフト分布Φ(x)と、下式(3)の関係にある。
第2の格子22の位置において、X線は、屈折角φ(x)に応じた量だけX方向に変位する。この変位量Δxは、X線の屈折角φ(x)が微小であることに基づいて、近似的に下式(4)で表される。
このように、変位量Δxは、位相シフト分布Φ(x)の微分値に比例する。そして、変位量Δxは、後述する強度変調信号の位相ズレ量ψ(x)(被検体Hがある場合とない場合とでの強度変調信号の位相ズレ量)と、下式(5)に示すように関連している。
上式(3)〜(5)から、強度変調信号の位相ズレ量ψ(x)を求めることにより、位相微分画像が得られることが分かる。なお、位相ズレ量ψ(x)に定数を乗じたり加算したりしたものを位相微分画像と定義してもよい。
図4において、G1像が第2の格子22に対してZ軸周りに僅かに傾斜するように、第1の格子21は、第2の格子22に対してZ軸周りに角度θだけ傾斜して配置される。この配置により、G2像には、Y方向に下式(6)で表される周期T(以下、モアレ周期Tという)を有するモアレ縞MSが生じる。
X線画像検出器13の画素部30のX方向の大きさをDx(以下、主画素サイズDxという)、Y方向の大きさをDy(以下、副画素サイズDyという)とする。第2の格子22の傾斜角θは、モアレ周期Tが副画素サイズDyのほぼ整数倍となるように設定されている。主画素サイズDxは、副画素サイズDyとほぼ同一とする。
図5において、Y方向に沿って並ぶM個の画素部30を1グループGr(x,n)とする。ここで、M,nは、正の整数である。nは、1グループGr(x,n)内の先頭の画素部30のy座標を表している。本実施形態では、1グループGr(x,n)内の画素数Mを、1モアレ周期Tに含まれる画素数ν(図4の例では、ν=5)と同一とする。
I(x,y)は、座標x,yの画素部30の画素値を示している。各画素値I(x,y)は、メモリ14に記憶された画像データから取得される。1グループGr(x,n)内の画素値I(x,n)〜I(x,n+M−1)は、図6に示すように、1周期分の強度変調信号を構成する。これは、1モアレ周期T内における画素部30のy座標に応じて、G1像と第2の格子22とのX方向への重なり度合いが異なり、強度変調量が変化することによる。したがって、1グループGr(x,n)内の画素値I(x,n)〜I(x,n+M−1)が構成する強度変調信号は、従来の縞走査法において、第1または第2の格子を、その格子方向にほぼ垂直な方向(X方向)に間欠移動させながら取得した1周期分の強度変調信号に相当する。
同図において、ψ(x,n)は、被検体Hが配置されていない場合の強度変調信号(破線)に対する被検体Hが配置された場合の強度変調信号(実線)の位相ズレ量を示している。
前述のように、X線画像検出器13は、列方向と行方向とのうち線状欠陥が少ない方向がY方向に沿うように配置されているため、線状欠陥は、図7に示すように、X方向に沿って生じやすい。グループGr(x,n)は、Y方向に延伸するように設定されるため、グループGr(x,n)内の画素値がすべて異常画素値となる確率は低い。
図8において、画像処理部15は、グループ設定部40、画素欠陥情報記憶部41、異常画素値除去部42、位相微分画像生成部43、位相コントラスト画像生成部44、及び正常画素値計数部45を有する。
グループ設定部40は、図9に示すように、X方向に並ぶ画素部30の各列について、グループGr(x,n)の先頭の座標をY方向に1画素ずつ変更しながら(nを1ずつ変更しながら)、設定可能な最終グループに達するまで、グループGr(x,n)を順に設定する。
画素欠陥情報記憶部41は、前述のように、X線画像検出器13の製造最終段階の検査工程で取得された線状欠陥を含む画素欠陥のマップ情報を記憶する。異常画素値除去部42は、画素欠陥情報記憶部41に記憶された画素欠陥のマップ情報に基づき、各グループGr(x,n)内の画素値から異常画素値を除去する。
位相微分画像生成部43は、異常画素値が除去された各グループGr(x,n)内の正常画素値のみを用いて強度変調信号を構成し、各強度変調信号の位相ズレ量ψ(x,n)を算出することにより位相微分画像を生成する。具体的には、位相微分画像生成部43は、縞走査法と同様な方法により位相微分画像を生成する。この方法は、例えば、「応用光学 光計測入門 谷田貝豊彦著 丸善株式会社 136〜138頁」により知られている。
位相微分画像生成部43は、下記の連立方程式(7)〜(9)を演算し、演算結果を次式(10)に適用することにより、位相ズレ量ψ(x,n)を算出する。
ここで、参照位相δkは、下式(11)で表される。
なお、位相微分画像生成部43は、連立方程式(7)〜(9)において、各左辺の和から異常画素値に対応する走査位置kの項を除外して演算を行う。ここで、異常画素値に対応する走査位置kの項が除外されることにより、参照位相δkは非等間隔となる。例えば、図10に示すように、画素値I(x,n+1)が異常画素値である場合には、k=1の項が連立方程式(7)〜(9)の各左辺の和から除外され、非等間隔の参照位相δ0,δ2〜δM−1に対応する正常画素値I(x,n),I(x,n+2)〜I(x,n+M−1)を用いて演算が行われる。これは、同図に示すように、正常画素値I(x,n),I(x,n+2)〜I(x,n+M−1)のみを正弦波でフィッティングし、このフィッティング波形の位相ズレ量ψ(x,n)を算出することに相当する。
位相コントラスト画像生成部44は、位相微分画像生成部43により生成された位相微分画像をX方向に沿って積分処理することにより、位相シフト分布を表す位相コントラスト画像を生成する。
正常画素値計数部45は、異常画素値除去部42により除去された異常画素値以外の正常画素値の数を、各グループGr(x,n)ごとにカウントする。システム制御部19は、正常画素値のカウント数Cが所定の基準値S(例えば3)より少ないグループGr(x,n)が存在するか否かを判定する。カウント数Cが基準値Sより少ないグループGr(x,n)が存在する場合には、そのグループGr(x,n)に対応する強度変調信号の位相ズレ量ψ(x,n)が精度よく算出されないため、システム制御部19は、モニタ18bを制御して、異常が生じている旨のエラーメッセージの表示(エラー報知)を行う。ここで、システム制御部19及びモニタ18bがエラー報知部を構成する。なお、このエラー報知を、音声やランプの点灯等により行ってもよい。
なお、上記の連立方程式(7)〜(9)でパラメータがa0,a1,a2の3個であることから分かるように、位相ズレ量ψ(x,y)の算出には、正常画素値の数が最低3個必要である。このため、基準値Sを“3”とすることが好ましい。また、いわゆるサンプリング定理を満たすためには、正常画素値の数は、最低4個必要であるため、より厳しく、基準値Sを“4”とすることも好ましい。
次に、図11に示すフローチャートに沿ってX線撮影装置10の作用を説明する。被検体Hが配置され、操作部18aにより撮影指示がなされると(ステップS10でYES)、X線源11によりX線が放射され、被検体H、第1及び第2の格子21,22を介して生成されたG2像がX線画像検出器13により検出され画像データが生成される(ステップS11)。
画像データがメモリ14に格納されると、グループ設定部40により、前述のグループGr(x,n)が設定され(ステップS13)、異常画素値除去部42により、画素欠陥情報記憶部41に記憶された画素欠陥のマップ情報に基づき、グループGr(x,n)内の画素値から異常画素値が除去される(ステップS14)。異常画素値が除去されると、正常画素値計数部45により、グループGr(x,n)内の正常画素値の数がカウントされる(ステップS15)。
システム制御部19により、グループGr(x,n)内の正常画素値のカウント数Cが所定の基準値Sより少ないか否かが判定され(ステップS16)、カウント数Cが基準値Sより少ない場合には(ステップS16でYES)、モニタ18bにエラーメッセージの表示が行われる(ステップS17)。この場合には、位相微分画像及び位相コントラスト画像の生成が行われずに動作が終了する。
一方、カウント数Cが基準値S以上の場合には(ステップS16でNO)、位相微分画像生成部43により、異常画素値が除去された各グループGr(x,n)内の正常画素値のみを用いて強度変調信号が構成され、その位相ズレ量ψ(x,n)が算出される(ステップS18)。
そして、ステップS13で設定したグループGr(x,n)が最終グループであるか否かが判定され(ステップS19)、最終グループでない場合には(ステップS19でNO)、グループ設定部40により、グループGr(x,n)の先頭の座標がY方向に1画素分変更され(ステップS20)、新たなグループGr(x,n)が設定される(ステップS13)。なお、ステップS20において、グループGr(x,n)をY方向に変更し、Y方向に関して設定可能な最終グループに達した場合には、X方向に隣接する画素列について同様にグループGr(x,n)が設定される。
このように、グループGr(x,n)が変更されながら、ステップS13〜S19が繰り返し実行される。そして、グループGr(x,n)が最終グループに達した場合には(ステップS19でYES)、1画面分の位相ズレ量ψ(x,n)の算出結果により位相微分画像が生成される(ステップS21)。
そして、位相コントラスト画像生成部44により、位相微分画像に対して積分処理が行われて、位相コントラスト画像が生成される(ステップS22)。モニタ18bには、位相微分画像及び位相コントラスト画像が画像表示される(ステップS23)。このように、画像データから異常画素値が除去されたうえで、位相ズレ量ψ(x,n)が算出されるため、精度よく位相微分画像及び位相コントラスト画像が生成される。
なお、上記実施形態では、図5に示すように、1グループGr(x,n)内の画素数Mを、1モアレ周期Tに含まれる画素数νと同一としているが、図12に示すように、1グループGr(x,n)内の画素数Mを、1モアレ周期Tに含まれる画素数νのN倍(ここで、Nは2以上の整数)と同一としてもよい。
また、図13に示すように、1グループGr(x,n)内の画素数Mは、1モアレ周期Tに含まれる画素数νまたはそのN倍と一致しなくてもよい。さらに、図14に示すように、1グループGr(x,n)内の画素数Mは、1モアレ周期Tに含まれる画素数νより少なくてもよい。
また、上記実施形態では、グループ設定部40は、図9に示すように、グループGr(x,n)をY方向に1画素ずつ変更しているが、グループGr(x,n)を、Y方向に2以上の画素数を単位として変更してもよい。
また、上記実施形態では、第2の格子22のX線吸収部22aの延伸方向をY方向とし、これに対して第1の格子21のX線吸収部21aの延伸方向を角度θだけ傾斜させているが、逆に、第1の格子21のX線吸収部21aの延伸方向をY方向とし、これに対して第2の格子22のX線吸収部22aの延伸方向を角度θだけ傾斜させてもよい。さらには、y方向に対して、第1の格子21のX線吸収部21aの延伸方向と、第2の格子22のX線吸収部22aの延伸方向とを逆方向に傾斜させ、両者が角度θをなすようにしてもよい。
また、上記実施形態では、カウント数Cが基準値Sより少ない場合に、位相微分画像及び位相コントラスト画像の生成を行わずに動作を終了しているが、カウント数Cが基準値Sより少ないグループが所定数に達した場合に、位相微分画像及び位相コントラスト画像の生成を行わずに動作を終了してもよい。また、これらの終了動作を行わずに、カウント数Cが基準値Sより少ないグループの検出のみを行い、位相微分画像及び位相コントラスト画像の生成を行ってもよい。この場合、位相微分画像を生成した後、カウント数Cが基準値Sより少ないグループに対応する画素を周辺画素の画素値を用いて補間すればよい。
また、上記実施形態では、X線画像検出器13は、第2の格子22の背後に近接して配置され、第2の格子22により生成されるG2像をほぼ等倍率で検出しているが、X線画像検出器13と第2の格子22との間に間隔を設けてもよい。X線画像検出器13と第2の格子22とのZ方向への間隔をL3とすると、下式(12)の倍率Rで拡大されたG2像がX線画像検出器13により検出される。
この場合には、X線画像検出器13により検出されるモアレ縞の周期T’は、上式(6)で表されるモアレ周期TのR倍(すなわちT’=RT)となる。このため、モアレ周期T’に基づいて、同様にグループGr(x,n)の設定を行えばよい。
(第2実施形態)
上記第1実施形態では、第1及び第2の格子21,22の格子面内方向への相対的な傾斜によりG2像にモアレ縞を生じさせているが、本実施形態では、第1及び第2の格子21,22を傾斜させずに、上式(1)の関係を僅かに崩すように第1及び第2の格子21,22の位置関係(距離L1,L2)、若しくは第1及び第2の格子21,22の格子ピッチp1,p2を調整することで、図15に示すように、X方向に周期を有するモアレ縞MSをG2像に生じさせる。
上記第1実施形態では、第1及び第2の格子21,22の格子面内方向への相対的な傾斜によりG2像にモアレ縞を生じさせているが、本実施形態では、第1及び第2の格子21,22を傾斜させずに、上式(1)の関係を僅かに崩すように第1及び第2の格子21,22の位置関係(距離L1,L2)、若しくは第1及び第2の格子21,22の格子ピッチp1,p2を調整することで、図15に示すように、X方向に周期を有するモアレ縞MSをG2像に生じさせる。
第2の格子22の位置でのG1像のX方向へのパターン周期p3は、第2の格子22の格子ピッチp2とは僅かにずれている。モアレ縞MSは、X方向に、下式(13)で表される周期Tを有する。
本実施形態では、X線画像検出器13は、第1及び第2の格子21,22の格子方向に直交する方向(X方向)に線状欠陥の方向が一致しないように配置されている。例えば、X線画像検出器13の信号線35に異常が生じ、列方向に沿って線状欠陥が生じた場合には、行方向がX方向に沿うようにX線画像検出器13を配置する。なお、列方向と行方向との双方に線状欠陥が存在する場合には、列方向と行方向とのうち線状欠陥が少ない方向がX方向に沿うようにX線画像検出器13を配置する。
本実施形態では、グループ設定部40は、図16に示すように、Y方向に並ぶ画素部30の各行について、グループGr(n,y)の先頭の座標をX方向に1画素ずつ変更しながら(nを1ずつ変更しながら)、設定可能な最終グループに達するまで、グループGr(n,y)を順に設定する。異常画素値除去部42は、画素欠陥情報記憶部41に記憶された画素欠陥のマップ情報に基づき、各グループGr(n,y)内の画素値から異常画素値を除去する。
前述のように、X線画像検出器13は、列方向と行方向とのうち線状欠陥が少ない方向がX方向に沿うように配置されているため、線状欠陥は、Y方向に沿って生じやすい。グループGr(n,y)は、X方向に延伸するように設定されるため、グループGr(n,y)内の画素値がすべて異常画素値となる確率は低い。
位相微分画像生成部43は、異常画素値が除去された各グループGr(n,y)内の正常画素値のみを用いて強度変調信号を構成し、各強度変調信号の位相ズレ量ψ(n,y)を算出することにより位相微分画像を生成する。位相ズレ量ψ(n,y)の算出方法は、第1実施形態と同様である。具体的には、下記の連立方程式(14)〜(16)を演算し、演算結果を次式(17)に適用すればよい。
位相微分画像生成部43は、連立方程式(14)〜(16)において、各左辺の和から異常画素値に対応する走査位置kの項を除外して演算を行う。ここで、異常画素値に対応する走査位置kの項が除外されることにより、参照位相δkは非等間隔となる。
なお、本実施形態においてもX線画像検出器13と第2の格子22との間に間隔L3を設けてもよい。この場合には、上式(13)で表されるモアレ周期Tに、上式(12)で表される倍率Rを乗じたモアレ周期T’に基づいてグループGr(n,y)を設定すればよい。
本実施形態においても第1実施形態と同様に、1グループGr(n,y)内の画素数Mは、1モアレ周期Tに含まれる画素数νまたはそのN倍と一致しなくてもよい。さらに、グループ再設定部61は、グループGr(n,y)を、X方向に2以上の画素数を単位として変更してもよい。本実施形態のX線撮影装置のその他の構成及び作用は、第1実施形態と同一であるため説明は省略する。
なお、第1実施形態で示した第1及び第2の格子21,22の格子面内方向への相対的な傾斜と、第2実施形態で示した第1及び第2の格子21,22の位置関係や格子ピッチのずれが同時に生じることにより、X方向とY方向とのいずれにも平行でない方向に周期を有するモアレ縞がG2像に生じることがある。この場合には、モアレ縞はX方向及びY方向に成分を有するため、第1または第2実施形態のいずれかの方法を用いることにより位相微分画像を生成することが可能である。また、この場合には、X方向とY方向とのいずれにも平行でない斜め方向に並ぶ複数の画素部30によりグループを構成し、同様に位相微分画像を生成することも可能である。
(第3実施形態)
上記第1及び第2実施形態では、第1及び第2の格子21,22を設定した後、モアレ縞の周期方向(グループの延伸方向)に一致する線状欠陥が少なくなるように、X線画像検出器13の方向を設定している。本実施形態では、逆に、X線画像検出器13の方向を設定した後、第1及び第2の格子21,22の格子面内方向への相対的な傾斜角度θや、及び第1及び第2の格子21,22のZ方向への位置関係を調整することにより、モアレ縞の周期方向(グループの延伸方向)に一致する線状欠陥が少なくなるように、モアレ縞の周期方向を設定する。本実施形態のX線撮影装置のその他の構成及び作用は、第1または第2実施形態と同一であるため説明は省略する。
上記第1及び第2実施形態では、第1及び第2の格子21,22を設定した後、モアレ縞の周期方向(グループの延伸方向)に一致する線状欠陥が少なくなるように、X線画像検出器13の方向を設定している。本実施形態では、逆に、X線画像検出器13の方向を設定した後、第1及び第2の格子21,22の格子面内方向への相対的な傾斜角度θや、及び第1及び第2の格子21,22のZ方向への位置関係を調整することにより、モアレ縞の周期方向(グループの延伸方向)に一致する線状欠陥が少なくなるように、モアレ縞の周期方向を設定する。本実施形態のX線撮影装置のその他の構成及び作用は、第1または第2実施形態と同一であるため説明は省略する。
なお、モアレ縞の周期方向グループの延伸方向は必ずしも一致する必要はなく、グループの延伸方向と線状欠陥が可能な限り一致しないように調整すればよい。こうすることで、グループ内の全ての画素値が異常画素値となることがなくなり、位相微分画像を精度よく生成することができる。
(その他の実施形態)
上記各実施形態では、グループ設定部40により設定されるグループの大きさは一定であるが、グループに含まれる正常画素値の数が所定数より少ない場合に、グループに含まれる正常画素値の数が所定数を超えるように、グループの大きさを拡大してもよい。
上記各実施形態では、グループ設定部40により設定されるグループの大きさは一定であるが、グループに含まれる正常画素値の数が所定数より少ない場合に、グループに含まれる正常画素値の数が所定数を超えるように、グループの大きさを拡大してもよい。
上記各実施形態では、画素欠陥情報記憶部41に記憶させる画素欠陥のマップ情報を、X線画像検出器13の製造最終段階の検査工程で取得しているが、X線撮影装置10にマップ情報の取得機能を設け、X線画像検出器13の製造後にマップ情報を取得可能としてもよい。また、X線撮影装置10の使用時にマップ情報を適宜更新可能としてもよい。
上記各実施形態では、被検体Hを配置した状態でのみ撮影を行なっているが、被検体Hを配置せずに同様に撮影を行うことにより、位相微分画像を生成して、これをオフセット画像として記憶しておき、被検体Hを配置して撮影が行われた場合に生成された位相微分画像からオフセット画像を減算するオフセット補正を行ってもよい。このように被検体Hを配置せずに行う撮影(プレ撮影)についても本発明を適用可能である。
マップ情報の更新及びオフセット補正を実施する場合には、例えば、図17のフローチャートに示す手順に沿って各処理を実行する。まず、画素欠陥情報記憶部41のマップ情報の更新を行い(ステップS30)、グループの延伸方向と線状欠陥が可能な限り一致しないようにX線画像検出器13や第1及び第2の格子21,22の位置調整を行う(ステップS31)。この後、被検体Hを配置せずに行うプレ撮影(ステップS32)と、被検体Hを配置して行う本撮影(ステップS33)とを順に実行する。
本撮影が完了すると、プレ撮影で得られたオフセット画像を、本撮影で得られた位相微分画像から減算するオフセット補正が行われる(ステップS34)。そして、オフセット補正後の位相微分画像が位相コントラスト画像生成部44により積分処理され、位相コントラスト画像が生成される(ステップS35)。このように、X線画像検出器13や第1及び第2の格子21,22の位置調整を行った場合には、その後にプレ撮影を行い、オフセット画像を取得することが好ましい。なお、プレ撮影と本撮影の順序は逆であってもよい。
また、上記各実施形態では、被検体HをX線源11と第1の格子21との間に配置しているが、被検体Hを第1の格子21と第2の格子22との間に配置してもよい。
また、上記各第実施形態では、X線源11から射出されるコーンビーム状のX線を射出するX線源11を用いているが、平行ビーム状のX線を射出するX線源を用いることも可能である。この場合には、上式(1)に代えて、p2=p1をほぼ満たすように第1及び第2の格子21,22を構成すればよい。
また、上記各実施形態では、X線源11から射出されたX線を第1の格子21に入射させており、X線源11は単一焦点であるが、X線源11の射出側直後に、WO2006/131235号公報等に記されたマルチスリット(線源格子)を設けることにより、X焦点を分散化してもよい。マルチスリットの格子線は、Y方向に延伸しており、第1及び第2の格子21,22の少なくとも一方の格子線と平行である。この場合、マルチスリットのX方向への格子ピッチp0は、下式(18)を満たす必要がある。ここで、距離L0は、マルチスリットから第1の格子21までのZ方向への距離を表す。
その他の構成や作用については、上記各実施形態と同一である。本実施形態では、マルチスリットの位置がX線焦点の位置となるため、上式(1)において、距離L1を、距離L0で置き換えればよい。
なお、このようにマルチスリットを設ける場合においてもX線画像検出器13と第2の格子22との間に間隔L3を設けてもよい。この場合には、上式(6)または上式(13)で表されるモアレ周期Tに、上式(12)で表される倍率Rを乗じたモアレ周期T’に基づいてグループGr(x,n)またはグループGr(n,y)を設定すればよい。マルチスリット23を設けた場合においても、第2の格子22により生成されるG2像は、X線源11のX線焦点11aを原点とし、X線焦点11aからX線画像検出器20までの距離に比例して拡大されるため、G2像の倍率Rについては、上式(14)をそのまま(距離L1を距離L0で置き換えずに)用いればよい。
また、上記各実施形態では、第1の格子21が入射X線を幾何光学的に投影するように構成しているが、WO2004/058070号公報等で知られているように、第1の格子21をタルボ効果が生じる構成としてもよい。第1の格子21でタルボ効果を生じさせるためには、X線の空間干渉性を高めるように、小焦点のX線源を用いるか、上記マルチスリットを用いて小焦点化すればよい。また、第1の格子21でタルボ効果を生じさせる場合には、第1の格子21を、吸収型格子に代えて、位相型格子とすることが可能である。
第1の格子21でタルボ効果が生じる場合には、第1の格子21の自己像(G1像)が、第1の格子21からZ方向下流にタルボ距離Zmだけ離れた位置に生じるため、第1の格子21から第2の格子22までの距離L2をタルボ距離Zmに設定する必要がある。
タルボ距離Zmは、第1の格子21の構成とX線のビーム形状とに依存する。第1の格子21が吸収型格子であり、X線源11から射出されるX線がコーンビーム状である場合には、タルボ距離Zmは、下式(19)で表される。ここで、mは正の整数である。この場合には、格子ピッチp1,p2は、上式(1)をほぼ満たすように設定される(ただし、マルチスリットを用いる場合には、距離L1は距離L0に置き換えられる)。
また、第1の格子21がX線にπ/2の位相変調を与える位相型格子であり、X線源11から射出されるX線がコーンビーム状である場合には、タルボ距離Zmは、下式(20)で表される。ここで、mは0または正の整数である。この場合には、格子ピッチp1,p2は、上式(1)をほぼ満たすように設定される(ただし、マルチスリットを用いる場合には、距離L1は距離L0に置き換えられる)。
また、第1の格子21がX線にπの位相変調を与える位相型格子であり、X線源11から射出されるX線がコーンビーム状である場合には、タルボ距離Zmは、下式(21)で表される。ここで、mは0または正の整数である。この場合には、G1像のパターン周期が第1の格子21の格子周期の1/2倍となるため、格子ピッチp1,p2は、次式(22)をほぼ満たすように設定される(ただし、マルチスリットを用いる場合には、距離L1は距離L0に置き換えられる)。
また、第1の格子21が吸収型格子であり、X線源11から射出されるX線が平行ビーム状である場合には、タルボ距離Zmは、下式(23)で表される。ここで、mは正の整数である。この場合には、格子ピッチp1,p2は、p2=p1の関係をほぼ満たすように設定される。
また、第1の格子21がX線にπ/2の位相変調を与える位相型格子であり、X線源11から射出されるX線が平行ビーム状である場合には、タルボ距離Zmは、下式(24)で表される。ここで、mは0または正の整数である。この場合には、格子ピッチp1,p2は、p2=p1の関係をほぼ満たすように設定される。
そして、第1の格子21がX線にπの位相変調を与える位相型格子であり、X線源11から射出されるX線が平行ビーム状である場合には、タルボ距離Zmは、下式(25)で表される。ここで、mは0または正の整数である。この場合には、G1像のパターン周期が第1の格子21の格子周期の1/2倍となるため、格子ピッチp1,p2は、p2=p1/2の関係をほぼ満たすように設定される。
上記各実施形態は、矛盾しない範囲で相互に組み合わせてもよい。本発明は、医療診断用の放射線撮影装置の他に、工業用の放射線撮影装置等に適用することが可能である。また、放射線は、X線以外に、ガンマ線等を用いることも可能である。
10 X線撮影装置
12 格子部
13 X線画像検出器
21 第1の格子
21a X線吸収部
21b X線透過部
22 第2の格子
22a X線吸収部
22b X線透過部
30 画素部
31 画素電極
33 ゲート走査線
35 信号線
12 格子部
13 X線画像検出器
21 第1の格子
21a X線吸収部
21b X線透過部
22 第2の格子
22a X線吸収部
22b X線透過部
30 画素部
31 画素電極
33 ゲート走査線
35 信号線
Claims (14)
- 放射線源から放射された放射線を通過させて第1の周期パターン像を生成する第1の格子と、
前記第1の周期パターン像を部分的に遮蔽してモアレ縞が生じた第2の周期パターン像を生成する第2の格子と、
行方向とこれに直交する列方向に沿って2次元配列された複数の画素部により、前記第2の周期パターン像を検出して画像データを生成する放射線画像検出器と、
前記行方向及び前記列方向のうち、線状欠陥が少ない第1の方向に並ぶ所定数の前記画素部を1グループとし、このグループを前記第1の方向に所定画素ずつ変更しながら順に設定するグループ設定部と、
前記各グループ内に含まれる画素値により構成された強度変調信号の位相ズレ量を算出することにより位相微分画像を生成する位相微分画像生成部と、
を備えることを特徴とする放射線撮影装置。 - 前記モアレ縞は、前記第2の格子を、第1の格子に対して格子面内方向に相対的に傾斜して配置することにより生成されたものであり、前記第1及び第2の格子の格子方向にほぼ直交していることを特徴とする請求項1に記載の放射線撮影装置。
- 前記モアレ縞は、前記第1及び第2の格子の対向方向の位置関係、または、前記第1及び第2の格子の格子ピッチを調整することにより生成されたものであり、前記第1及び第2の格子の格子方向にほぼ平行であることを特徴とする請求項1に記載の放射線撮影装置。
- 前記第1の方向は、前記モアレ縞にほぼ直交していることを特徴とする請求項1から3いずれか1項に記載の放射線撮影装置。
- 前記モアレ縞は、前記第2の格子を、第1の格子に対して格子面内方向に相対的に傾斜して配置し、かつ、前記第1及び第2の格子の対向方向の位置関係、または、前記第1及び第2の格子の格子ピッチを調整することにより生成されたものであり、前記第1及び第2の格子の格子方向に直交せず、かつ平行でないことを特徴とする請求項1に記載の放射線撮影装置。
- 前記放射線画像検出器の画素欠陥のマップ情報を記憶する画素欠陥情報記憶部と、前記マップ情報に基づき、前記各グループ内の画素値から異常画素値を除去する異常画素値除去部と、を備え、
前記位相微分画像生成部は、前記各グループ内に含まれる前記異常画素値以外の正常画素値に基づいて前記強度変調信号の位相ズレ量を算出し、前記位相微分画像を生成することを特徴とする請求項1から5いずれか1項に記載の放射線撮影装置。 - 前記各グループ内の前記正常画素値の数をカウントする正常画素値計数部を備えることを特徴とする請求項6に記載の放射線撮影装置。
- 前記正常画素値の数が基準値より少ないグループが存在する場合にエラー報知を行うエラー報知部を備えることを特徴とする請求項7に記載の放射線撮影装置。
- 前記グループを構成する前記画素部の数は、前記モアレ縞の周期の整数倍に含まれる前記画素部の数に等しいことを特徴とする請求項1から8いずれか1項に記載の放射線撮影装置。
- 前記位相微分画像を積分処理して位相コントラスト画像を生成する位相コントラスト画像生成部を備えることを特徴とする請求項1から9いずれか1項に記載の放射線撮影装置。
- 前記第1の格子は、吸収型格子であり、入射した放射線を幾何光学的に投影することにより前記第1の周期パターン像を生成することを特徴とする請求項1から10いずれか1項に記載の放射線撮影装置。
- 前記第1の格子は、吸収型格子または位相型格子であり、入射した放射線にタルボ効果を生じさせて前記第1の周期パターン像を生成することを特徴とする請求項1から10いずれか1項に記載の放射線撮影装置。
- 前記放射線源から放射された放射線を部分的に遮蔽して焦点を分散化するマルチスリットを備えることを特徴とする請求項1から12いずれか1項に記載の放射線撮影装置。
- 放射線源から放射された放射線を第1の格子を通過させて第1の周期パターン像を生成するステップと、
前記第1の周期パターン像を第2の格子により部分的に遮蔽してモアレ縞が生じた第2の周期パターン像を生成するステップと、
行方向とこれに直交する列方向に沿って2次元配列された複数の画素部を有する放射線画像検出器により前記第2の周期パターン像を検出して画像データを生成するステップと、
前記行方向及び前記列方向のうち、線状欠陥が少ない第1の方向に並ぶ所定数の前記画素部を1グループとし、このグループを前記第1の方向に所定画素ずつ変更しながら順に設定するステップと、
前記各グループ内に含まれる画素値により構成された強度変調信号の位相ズレ量を算出することにより位相微分画像を生成するステップと、
を備えることを特徴とする放射線撮影方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011272142A JP2013123460A (ja) | 2011-12-13 | 2011-12-13 | 放射線撮影装置及び放射線撮影方法 |
PCT/JP2012/079227 WO2013088878A1 (ja) | 2011-12-13 | 2012-11-12 | 放射線撮影方法及び装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011272142A JP2013123460A (ja) | 2011-12-13 | 2011-12-13 | 放射線撮影装置及び放射線撮影方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2013123460A true JP2013123460A (ja) | 2013-06-24 |
Family
ID=48612333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011272142A Pending JP2013123460A (ja) | 2011-12-13 | 2011-12-13 | 放射線撮影装置及び放射線撮影方法 |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP2013123460A (ja) |
WO (1) | WO2013088878A1 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019239624A1 (ja) * | 2018-06-15 | 2019-12-19 | 株式会社島津製作所 | X線イメージング装置 |
US20240244346A1 (en) * | 2022-04-25 | 2024-07-18 | Beijing Xiaomi Mobile Software Co., Ltd. | Photographing device and control method thereof |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008200361A (ja) * | 2007-02-21 | 2008-09-04 | Konica Minolta Medical & Graphic Inc | X線撮影システム |
WO2012144317A1 (ja) * | 2011-04-20 | 2012-10-26 | 富士フイルム株式会社 | 放射線撮影装置及び画像処理方法 |
-
2011
- 2011-12-13 JP JP2011272142A patent/JP2013123460A/ja active Pending
-
2012
- 2012-11-12 WO PCT/JP2012/079227 patent/WO2013088878A1/ja active Application Filing
Also Published As
Publication number | Publication date |
---|---|
WO2013088878A1 (ja) | 2013-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5378335B2 (ja) | 放射線撮影システム | |
JP5731214B2 (ja) | 放射線撮影システム及びその画像処理方法 | |
JP5475737B2 (ja) | 放射線撮影装置及び画像処理方法 | |
JP5548085B2 (ja) | 回折格子の調整方法 | |
JP2011224329A (ja) | 放射線撮影システム及び方法 | |
JP2012090945A (ja) | 放射線検出装置、放射線撮影装置、放射線撮影システム | |
JP2011218147A (ja) | 放射線撮影システム | |
JP5783987B2 (ja) | 放射線撮影装置 | |
JP2011224330A (ja) | 放射線撮影システム及びそのオフセット補正方法 | |
WO2012057022A1 (ja) | 放射線撮影システム及び放射線撮影方法 | |
JP2013116313A (ja) | 放射線撮影装置及び放射線撮影方法 | |
WO2013038881A1 (ja) | 放射線撮影装置及び画像処理方法 | |
JP2011206188A (ja) | 放射線撮影システム及び方法 | |
WO2013088878A1 (ja) | 放射線撮影方法及び装置 | |
JP2013042788A (ja) | 放射線撮影装置及びアンラップ処理方法 | |
JP5635169B2 (ja) | 放射線撮影システム | |
JP5610480B2 (ja) | 放射線画像処理装置及び方法 | |
WO2013099467A1 (ja) | 放射線撮影方法及び装置 | |
JP2013106883A (ja) | 放射線撮影装置及び画像処理方法 | |
WO2013027519A1 (ja) | 放射線撮影装置及びアンラップ処理方法 | |
WO2013027536A1 (ja) | 放射線撮影装置及び放射線撮影方法 | |
WO2012147749A1 (ja) | 放射線撮影システム及び放射線撮影方法 | |
JP2012035050A (ja) | 放射線撮影システム及びその画像処理方法 | |
JP2012228301A (ja) | 放射線撮影装置 | |
JP2013063098A (ja) | 放射線撮影装置及び画像処理方法 |