JP4408664B2 - コーンビームx線ct装置及びそれに用いるファントム - Google Patents
コーンビームx線ct装置及びそれに用いるファントム Download PDFInfo
- Publication number
- JP4408664B2 JP4408664B2 JP2003289681A JP2003289681A JP4408664B2 JP 4408664 B2 JP4408664 B2 JP 4408664B2 JP 2003289681 A JP2003289681 A JP 2003289681A JP 2003289681 A JP2003289681 A JP 2003289681A JP 4408664 B2 JP4408664 B2 JP 4408664B2
- Authority
- JP
- Japan
- Prior art keywords
- ray
- image
- rotation
- dimensional
- angle
- 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.)
- Expired - Fee Related
Links
- 238000004364 calculation method Methods 0.000 claims description 62
- 238000012937 correction Methods 0.000 claims description 39
- 230000005855 radiation Effects 0.000 claims 1
- 238000000034 method Methods 0.000 description 42
- 238000012545 processing Methods 0.000 description 23
- 238000004519 manufacturing process Methods 0.000 description 12
- 230000005540 biological transmission Effects 0.000 description 10
- 238000003384 imaging method Methods 0.000 description 10
- 238000005259 measurement Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 238000001914 filtration Methods 0.000 description 8
- 238000007781 pre-processing Methods 0.000 description 7
- 238000010521 absorption reaction Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 238000003860 storage Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 230000005484 gravity Effects 0.000 description 4
- 239000002184 metal Substances 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
- A61B6/584—Calibration using calibration phantoms determining position of components of the apparatus or device using images of the phantom
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
このような2次元X線検出器を用いたコーンビームX線CT装置の画像再構成演算アルゴリズムとして、フェルドカンプの方法が代表的である(例えば、非特許文献1参照。)。
なお、発明の実施の形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
座標ねじれ角補正手段203は、2次元X線検出器の画像上での座標軸と実際の回転中心軸との角度の差(座標ねじれ角)の補正を行う画像回転手段である。
逆投影手段205は、フィルタリングの後の画像データに対して、文献1に記載のフェルトカンプの方法を用いて逆投影処理を行う公知の手段である。ただし本実施の形態では、[0007]で述べたように、回転前半と回転後半で幾何学パラメータが異なると考え、それぞれの幾何学パラメータを用いて、それぞれのX線透過像を用いて逆投影処理を行う。
ここで、回転軌道面(ミッドプレーン)7とは、[0016]ではX線源の焦点の軌道を含む平面であると定義したが、次の定義でも特徴づけられる平面である。被検体上の点(x,y,z)に着目する。回転によりX線源と点(x,y,z)との距離は変化するため、2次元検出器面への投影点は遠近法の原理から分かるように一般に上下方向にも移動する。しかしながら、X線源を含み回転中心軸に垂直な面では、回転によりX線源と点(x,y,z)との距離が変化しても2次元検出器面への投影点は上下方向には移動しない。このような面と定義しても回転軌道面(ミッドプレーン)7は特徴づけられる。
(ステップ20)先ず、ホールチャート20をX線イメージインテンシファイア5の前面に固定して回転撮影を行い、各角度におけるホールチャート歪み投影像、例えば150枚の収集を行う。
(ステップ301)次に画像歪み補正テーブル生成手段301により、(ステップ20)で得られた150枚のホールチャート画像歪み投影像を基にホールチャート画像歪み補正テーブルの計算を行う。
(ステップ302)次に画像歪み補正テーブル保管手段302に、(ステップ301)で計算されたホールチャート画像歪み補正テーブルを保管する。ここで保管されたホールチャート画像歪み補正テーブルは、画像歪み補正手段202でX線透視像の画像歪みを補正する際に用いられる。
(ステップ304)回転中心軸位置及び回転角度ピッチ計算手段304により、回転中心軸位置及び回転角度ピッチの実際の値(設計値との差)を計算する。計算の具体的な処理の詳細は、図6を用い後述する。ここで求められた実際の回転中心軸位置及び回転角度ピッチは、逆投影手段205で実際のX線透視像を逆投影処理をする際に用いられる。
(ステップ306)座標ねじれ角計算手段306により、[0021]で説明した座標ねじれ角を計算する。計算の具体的な処理の詳細は、図15を用い後述する。ここで求められた座標ねじれ角は、座標ねじれ角補正手段203で座標ねじれ角を補正する際に用いられる。
(ステップ201)回転軌道面(ミッドプレーン)7に限定した投影データを用いて、前処理手段201で前処理を行う。
(ステップ204)(ステップ202)で回転軌道面(ミッドプレーン)7に限定して画像歪み補正処理を行ったデータについて、フィルタリング手段204でフィルタリング処理を行う。
(ステップ205)(ステップ204)で回転軌道面(ミッドプレーン)7に限定してフィルタリング処理を行ったデータについて、逆投影手段205で逆投影処理を行う。
(ステップ402)撮影系ジオメトリパラメータの1つである回転角度ピッチの初期値を、例えば1.33°(=200°/150枚)に設定する。
(ステップ403)回転中心軸位置(以下、centerと称す)を初期位置(設計位置)に設定する。
(ステップ303)ワイヤー状ファントムの回転前半の投影データによる回転軌道面画像再構成を行う。
(ステップ410)再構成画像には複数のワイヤー断面が縦横に配置されているが、回転中心軸位置計算処理においては、最も画像***付近に配置されているワイヤーに注目する。そして、この中心ワイヤーの半径を算出する。
(ステップ411)(ステップ410)で算出した弧の半径がゼロと見なせるかを判定する。見なせれば(ステップ420)へ、見なせなければ(ステップ412)へ移行する。
(ステップ412)(ステップ411)で弧の半径がゼロと見なせなかった場合、図7の場合か図8の場合かを判断し、それに応じてcenterが正しくなるよう加減算して補正する。そして、(ステップ303)からの過程を再度行う。(ステップ411)で弧の半径がゼロと見なせた場合、回転角度ピッチの計算に移行する。
(ステップ421)(ステップ420)で算出した周辺ワイヤーの弧の半径がゼロと見なせるかを判定する。見なせれば再構成画像は図11のようになり、回転前半の回転中心軸位置および回転角度ピッチの計算を終了し、見なせなければ(ステップ422)へ移行する。(ただし、図11において、ワイヤーの再構成画像の断面が斜めの棒状で、点状に結像しないのは、左上から右下方向への投影データが欠落しているためである。)
(ステップ422)(ステップ421)で弧の半径がゼロと見なせなかった場合、図9の場合か図10の場合かを判断し、それに応じて回転角度ピッチが正しくなるよう加減算して補正する。そして、(ステップ303)からの過程を再度行う。
(ステップ501)先ず、再構成演算における回転開始角度方向を設計値−100°とする。
(ステップ502)回転前半の投影データから図11に示すような回転軌道面再構成画像を生成する。
(ステップ504)(ステップ503)で抽出されたワイヤー領域より各々のワイヤー断面の重心座標を算出する。
(ステップ505)ワイヤーの再構成断面が正確に水平に配置されているか、すなわち横方向に並ぶワイヤー断面の重心の上下座標が等しいかを判定し、等しければ(ステップ510)へ、等しくなければ(ステップ506)へ移行する。
(ステップ506)(ステップ505)で横方向に並ぶワイヤー断面の重心の上下座標が等しくないと判定された場合、再構成演算における前半回転開始角度方向をその等しくない程度に応じて加減算し、(ステップ502)からの過程を再度行う。
(ステップ511)更に、回転後半の開始角度を天井方向0°に設定する。
(ステップ512)回転後半の投影データから図12に示すようなミッドプレーン再構成画像を生成する。
(ステップ513)閾値処理により再構成画像からワイヤー領域を抽出する。
(ステップ515)ワイヤーの再構成断面が正確に水平に配置されているか、すなわち横方向に並ぶワイヤー断面の重心の上下座標が等しいかを判定し、等しければ(ステップ520)へ、等しくなければ(ステップ516)へ移行する。
(ステップ516)(ステップ515)で横方向に並ぶワイヤー断面の重心の上下座標が等しくないと判定された場合、再構成演算における後半回転開始角度方向を加減算し、(ステップ512)からの過程を再度行う。
(ステップ520)(ステップ515)で回転後半のワイヤーの再構成断面が正確に水平に配置されていると判断された場合、回転前半のワイヤーの再構成断面と回転後半のワイヤーの再構成断面を重ね合わせる。
(ステップ522)(ステップ521)で回転前半の再構成断面のワイヤー重心座標と回転後半の再構成断面のワイヤー重心座標が一致していないと判定された場合、重心座標のずれの分だけ重ね合わせ基準位置を補正して、(ステップ512)からの過程を再度行う。
(ステップ601)初めに、座標ねじれ角の初期値を0°に設定する。
(ステップ602)上記各手段で求めた幾何学パラメータを使用して、回転軌道面7に限らず、全ワイヤー投影データより3次元画像再構成処理を行う。
(ステップ604)閾値処理により再構成画像からワイヤー領域を抽出する。
(ステップ605)ワイヤーの弧の半径を算出する。
(ステップ607)(ステップ606)でワイヤーの弧の半径をゼロと見なせない場合、弧の向きに応じて座標ねじれ角の値を増減し、(ステップ602)からの過程を再度行う。
305…回転開始角度及び重ね合わせ基準位置計算手段
306…座標ねじれ角計算手段
Claims (3)
- 被検体にX線を照射するX線源と、そのX線源に対向して配置され被検体のX線透視像を得る2次元X線検出器と、これらの2次元X線検出器とX線源を同一回転中心の円軌道面上で回転移動させながら、回転中心軸上に配置された被検体のX線透視像を前記2次元X線検出器により得、得られたX線透視像を基に3次元的X線CT像を画像再構成する画像再構成手段と、画像再構成された3次元的X線CT像を表示する表示手段を備えたコーンビームX線CT装置において、
前記回転移動の回転中心軸、前記回転移動の回転開始角度及び回転角度ピッチに含まれる誤差を、前記回転中心軸に略平行に配置した複数本のワイヤーを含むファントムの再構成画像を用いて計算する誤差計算手段とを備え、前記画像再構成手段は、該計算された誤差を用いて、前記X線透視像を逆投影し、3次元的X線CT像を画像再構成することを特徴とするコーンビームX線CT装置。 - 前記2次元X線検出器の前面に予めホールチャートファントムを設定し、設置されたホールチャートファントムについてX線透視像を、前記2次元X線検出器により撮影し、その撮影されたホールチャートの前記被検体の体軸方向に沿った軸と、前記回転中心軸の前記2次元X線検出器の投影像の座標ねじれ角を計算する座標ねじれ角計算手段を備え、
前記画像再構成手段は、該計算された座標ねじれ角を用い、前記X線透視像を補正することを特徴とする請求項1に記載のコーンビームX線CT装置。 - 前記複数本のワイヤーを含むファントムの再構成画像に基づいて、前記回転中心軸の方向と、前記2次元X線検出器により撮影された前記X線透視像との座標ねじれ角を計算する座標ねじれ角計算手段を備え、
前記画像再構成手段は、該計算された座標ねじれ角を用い、前記X線透視像を回転補正する座標ねじれ角補正手段を備えることを特徴とする請求項1に記載のコーンビームX線CT装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003289681A JP4408664B2 (ja) | 2003-08-08 | 2003-08-08 | コーンビームx線ct装置及びそれに用いるファントム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003289681A JP4408664B2 (ja) | 2003-08-08 | 2003-08-08 | コーンビームx線ct装置及びそれに用いるファントム |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2005058309A JP2005058309A (ja) | 2005-03-10 |
JP2005058309A5 JP2005058309A5 (ja) | 2006-08-17 |
JP4408664B2 true JP4408664B2 (ja) | 2010-02-03 |
Family
ID=34367928
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003289681A Expired - Fee Related JP4408664B2 (ja) | 2003-08-08 | 2003-08-08 | コーンビームx線ct装置及びそれに用いるファントム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4408664B2 (ja) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101370430A (zh) | 2006-01-24 | 2009-02-18 | 株式会社岛津制作所 | X射线摄像装置 |
JP4891623B2 (ja) * | 2006-02-02 | 2012-03-07 | 株式会社東芝 | X線ct装置 |
EP2074383B1 (en) * | 2006-09-25 | 2016-05-11 | Mazor Robotics Ltd. | C-arm computerized tomography |
JP5022820B2 (ja) * | 2007-08-13 | 2012-09-12 | 日立Geニュークリア・エナジー株式会社 | 放射線検査装置及びそれを用いた配管検査方法 |
JP5011085B2 (ja) * | 2007-12-19 | 2012-08-29 | 日立Geニュークリア・エナジー株式会社 | 検査装置,検査装置のアライメント方法及び校正用ファントムの製作方法 |
WO2009128129A1 (ja) * | 2008-04-14 | 2009-10-22 | 株式会社 島津製作所 | X線診断装置 |
JP5210726B2 (ja) * | 2008-06-24 | 2013-06-12 | 株式会社東芝 | X線ct装置 |
CN102123664B (zh) * | 2008-08-13 | 2014-05-07 | 皇家飞利浦电子股份有限公司 | 使用旋转中心寻找算法进行环形伪影校正的校准方法 |
JP5415885B2 (ja) * | 2009-09-30 | 2014-02-12 | 富士フイルム株式会社 | 放射線ct装置および画像処理装置 |
KR101141045B1 (ko) * | 2010-03-16 | 2012-05-15 | 한국전기연구원 | X-선 유방암진단 장치의 좌표측정 시스템 및 그 방법 |
JP5537226B2 (ja) * | 2010-03-31 | 2014-07-02 | 株式会社日立メディコ | 放射線撮像装置 |
KR101214377B1 (ko) | 2011-09-23 | 2012-12-21 | 한국전기연구원 | 교정 팬텀을 이용한 엑스선 초점 위치 교정 방법 |
CN103519833B (zh) * | 2013-06-05 | 2015-10-07 | 东南大学 | 一种旋转c型臂x射线机的三维校正重建方法 |
JP6179394B2 (ja) * | 2013-12-27 | 2017-08-16 | 株式会社島津製作所 | 放射線撮影装置 |
CN106821405B (zh) * | 2017-01-23 | 2019-12-06 | 深圳先进技术研究院 | 一种x光机的参数标定方法、装置及*** |
EP3946059B1 (en) * | 2019-03-25 | 2022-10-19 | Koninklijke Philips N.V. | X-ray ring markers for x-ray calibration |
CN113749681A (zh) * | 2021-11-10 | 2021-12-07 | 极限人工智能(北京)有限公司 | 评估cbct机图像精度的工装 |
-
2003
- 2003-08-08 JP JP2003289681A patent/JP4408664B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2005058309A (ja) | 2005-03-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4408664B2 (ja) | コーンビームx線ct装置及びそれに用いるファントム | |
JP5019879B2 (ja) | X線ct装置、画像処理プログラム、及び画像処理方法 | |
KR101523422B1 (ko) | 방사선 촬상 장치 및 동장치에 이용하는 팬텀 | |
JP4537129B2 (ja) | トモシンセシス用途における対象物を走査するためのシステム | |
JP3548339B2 (ja) | X線撮影装置 | |
US7244063B2 (en) | Method and system for three dimensional tomosynthesis imaging | |
JP5830753B2 (ja) | X線ct撮影装置及びx線ct画像の表示方法 | |
US9532755B2 (en) | X-ray computed tomographic imaging apparatus and method for same | |
US20110075798A1 (en) | Method for correcting truncated projection data | |
JP5618292B2 (ja) | X線ct撮影装置及びx線ct画像の表示方法 | |
JP3540916B2 (ja) | 3次元x線ct装置 | |
JP2006204329A (ja) | X線断層撮影装置 | |
JP2008510504A (ja) | 傾斜した構成を用いた3次元再現 | |
JP4573593B2 (ja) | X線画像補正方法及び装置 | |
JP2009036660A (ja) | 断層撮影装置 | |
JP5016231B2 (ja) | 撮像の幾何学的パラメータを決定する方法及び装置 | |
US9066687B2 (en) | Panoramic dental x-ray unit | |
US20050133708A1 (en) | Method and system for three dimensional tomosynthesis imaging | |
JP4610304B2 (ja) | X線ct装置 | |
JP3548306B2 (ja) | X線断層撮影装置 | |
JP5572521B2 (ja) | X線ct装置、x線ct装置の画像再構成方法 | |
JP4479503B2 (ja) | 断層撮影装置 | |
JP4165319B2 (ja) | コンピュータ断層撮像方法及び装置 | |
JP4421917B2 (ja) | コーンビームx線ct装置 | |
JPH09234192A (ja) | X線撮影装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060703 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060703 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090622 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090707 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090806 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090910 |
|
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: 20091020 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20091110 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121120 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131120 Year of fee payment: 4 |
|
LAPS | Cancellation because of no payment of annual fees |