JP2015208601A - X線ct装置、画像処理装置、及び投影データ生成方法 - Google Patents
X線ct装置、画像処理装置、及び投影データ生成方法 Download PDFInfo
- Publication number
- JP2015208601A JP2015208601A JP2014093627A JP2014093627A JP2015208601A JP 2015208601 A JP2015208601 A JP 2015208601A JP 2014093627 A JP2014093627 A JP 2014093627A JP 2014093627 A JP2014093627 A JP 2014093627A JP 2015208601 A JP2015208601 A JP 2015208601A
- Authority
- JP
- Japan
- Prior art keywords
- projection data
- focus
- virtual
- ray
- view
- 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 claims description 36
- 238000012545 processing Methods 0.000 claims description 62
- 238000004364 calculation method Methods 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 12
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 238000003786 synthesis reaction Methods 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 9
- 230000001678 irradiating effect Effects 0.000 claims description 4
- 230000007547 defect Effects 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims 1
- 238000003384 imaging method Methods 0.000 description 20
- 238000010586 diagram Methods 0.000 description 14
- 230000006870 function Effects 0.000 description 12
- 230000005540 biological transmission Effects 0.000 description 9
- 230000000007 visual effect Effects 0.000 description 7
- 230000014509 gene expression Effects 0.000 description 5
- 230000006872 improvement Effects 0.000 description 4
- 238000013480 data collection Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000010894 electron beam technology Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
Images
Abstract
【解決手段】 X線CT装置1の画像処理装置122は、X線源101における複数の焦点位置から被検体に対してX線を照射して撮影した焦点シフト投影データ31を焦点位置毎に分割する際に生じる欠損ビューに仮想ビューを挿入する。仮想ビューデータ生成部23は、仮想ビューのデータを他方の焦点位置についての焦点別投影データ32B(32A)を用い、かつ焦点位置とX線検出器106との位置関係に基づいて算出する。仮想ビューデータが挿入された各焦点別投影データ32A,32B(仮想投影データ)を組み合わせることにより、空間分解能を向上させた投影データを得る。
【選択図】図8
Description
まず、図1を参照して、X線CT装置1の全体構成について説明する。
図1に示すように、X線CT装置1は、スキャンガントリ部100、寝台105、及び操作卓120を備える。
図5に示すように、焦点Aは基準焦点位置(焦点オフセットがない場合の各ビューの基準位置)から距離offset_aだけ回転方向負方向へシフトし、焦点Bは上記基準位置から距離offset_bだけ回転方向正方向へシフトするものとする。
また、ビュー間の角度はΦ=2π/View数である。X線源101(焦点A,焦点B)と回転中心Oとの間の距離をSODとすると、焦点Aと回転中心Oと基準焦点位置との間の角度αと、焦点Bと回転中心Oと基準焦点位置との間の角度βは以下の式(1)、(2)のように表せる。
View1・・・焦点A(Φ−α)、
View2・・・焦点B(2Φ+β)、
View3・・・焦点A(3Φ−α)、
View4・・・焦点B(4Φ+β)、
View5・・・焦点A(5Φ−α)、
View6・・・焦点B(6Φ+β)、
View7・・・焦点A(7Φ−α)、
View8・・・焦点B(8Φ+β)、…。
View1・・・f(Φ−α)、
View3・・・f(3Φ−α)、
View5・・・f(5Φ−α)、
View7・・・f(7Φ−α)、となる。
View2、View4、View6、View8、…は欠損ビューとなる。
View2・・・g(2Φ+β)、
View4・・・g(4Φ+β)、
View6・・・g(6Φ+β)、
View8・・・g(8Φ+β)、…となる。
View1、View3、View5、View7、…は欠損ビューとなる。
生成された各焦点別投影データ32A,32Bは、仮想ビューデータ生成部23に入力される。
View2・・・f(2Φ−α)、
View4・・・f(4Φ−α)、
View6・・・f(6Φ−α)、
View8・・・f(8Φ−α)、…とする。
View1・・・g(Φ+β)、
View3・・・g(3Φ+β)、
View5・・・g(5Φ+β)、
View7・・・g(7Φ+β)、…とする。
図7(b)に示すように、焦点Aの焦点別投影データ32Aにおける仮想ビューView2のビューデータを生成する場合は、仮想ビューデータ生成部23は、他方の焦点Bの投影データ32Bにおける実ビューg(2Φ+β)のデータを参照する。ただし、仮想ビューの角度(2Φ−α)と実ビューの角度(2Φ+β)は角度θ(以下の式(3))だけ異なっているので、例えば隣接する実ビューのデータ等を利用した補間演算を行って、正確に補うものとする。
投影データ合成部24は、焦点A仮想投影データ33Aと焦点B仮想投影データ33Bを取得し、実ビュー及び仮想ビューの配置が幾何学的に正しくなるように組み合わせる。本発明適用後の理想的なビュー配置を図8(b)に示す。
図14は、本発明に係るX線CT装置1が実行する撮影・再構成処理全体の流れを説明するフローチャートである。
次に、図15を参照して本発明の第2の実施の形態について説明する。
なお、以下の説明において、第1の実施の形態と同一の構成及び動作については説明を省略し、第1の実施の形態と同一の符号を付して説明する。
次に、図16〜図22を参照して本発明の第3の実施の形態について説明する。
第1及び第2の実施の形態では、実際にX線焦点を2点にシフトさせて撮影を行い、ビュー数を2倍に増加させた仮想投影データを生成した。第3の実施の形態では、任意の位置に仮想焦点を設定し、仮想焦点から投影されたと同等のジオメトリを有する仮想ビューデータを実在の焦点から得た投影データ(実投影データ)を用いて作成する。これにより、FFS機能を持たないX線源においても投影データのビュー密度を任意の目標数まで増加させ、空間分解能を所望の目標値まで向上させる。
図16(a)は、FFS機能を用いず、焦点Aのみでスキャンする場合の焦点位置を示している。このスキャンで得られた投影データの空間分解能が2倍相当となるように仮想投影データを作成する例を説明する。
画像処理装置122は、ステップS203で求めた仮想焦点数だけ、オフセット量の異なる仮想焦点を設定する。そして上述したように、各仮想焦点についての仮想投影データ(仮想ビューデータ)を、ステップS202で入力した実投影データ及び仮想焦点に関するジオメトリから求める。そして、実投影データと各仮想投影データとを組み合わせ、X線密度の高い投影データを生成する(ステップS204)。その後、画像処理装置122はステップS204で作成した投影データを用いて画像を再構成する(ステップS205)。
これにより、FFS機能を用いずにスキャンした場合であっても、所望の空間分解能に応じた投影データを生成できる。
次に、本発明のX線CT装置1の第4の実施の形態について説明する。
FFS焦点制御において、なんらかの原因で実際の焦点位置が本来の焦点位置ではない位置にずれてしまうことがある。
第2の仮想焦点は、例えば、前のビューを参考にする。現在操作しているビューをview_nとすると、図23(e)、(f)に示すように、前のビューview_n-1で得たデータを参考データdataBとする。そして図23(g)に示すように、dataAとdataBとの補間演算により、view_nにおける理想焦点位置B(offsetB)のデータdataCを作成することができる。このときdataBの焦点がエラーであるかどうかは問わない。例えば、線形補間によりview_nにおけるoffsetBのデータdataC(理想焦点位置Bのデータ)を作成する場合、dataAの傾きをγ、dataBの傾きをθとすると、理想焦点位置BのデータdataCは以下の式(14)で表現される。
次に、図24〜図27を参照して、本発明の第5の実施の形態について説明する。
第1〜第4の実施の形態では、焦点位置のシフト方向を回転盤102の回転方向(X線検出器106のチャンネル方向)としたが、第5の実施の形態では、焦点位置のシフト方向を回転盤102の回転方向と直交する方向(X線検出器106のスライス方向(Z方向))とする場合について説明する。
View1・・・f(Φ,s1+offset_s)、
View2・・・g(2Φ,t1+offset_t)、
View3・・・f(3Φ,s1+offset_s)、
View4・・・g(4Φ,t1+offset_t)、
View5・・・f(5Φ,s1+offset_s)、
View6・・・g(6Φ,t1+offset_t)、
View7・・・f(7Φ,s1+offset_s)、
View8・・・g(8Φ,t1+offset_t)、…である。
View1・・・h(Φ,t1+offset_y)、
View2・・・i(2Φ,s1+offset_x)、
View3・・・h(3Φ,t1+offset_y)、
View4・・・i(4Φ,s1+offset_x)、
View5・・・h(5Φ,t1+offset_y)、
View6・・・i(6Φ,s1+offset_x)、
View7・・・h(7Φ,t1+offset_y)、
View8・・・i(8Φ,s1+offset_x)。
実際の撮影系(焦点位置offset_tと検出器106)の位置関係が図25(a)に示すような配置となっていると仮定する。このとき回転中心部のスライス厚dは三角形の相似の関係により、以下の式(15)の関係となっている。
なお、図25(a)と同じ検出器位置で、図25(b)に示すようにスライスs1を通るような焦点位置がoffset_sである。
次に、図28〜図31を参照して、本発明の第6の実施の形態について説明する。
第6の実施の形態では、焦点位置が理想位置からオフセットし、エラー焦点位置にある場合(第4の実施の形態)の投影データの補正方法の別の例について説明する。
まず焦点がerror_offsetにある場合のビュー角度のずれ量fを求める。図30に示すように、ビュー角度のずれ量fは、理想的な焦点ideal_offsetと、距離dだけずれた焦点error_offsetとの関係に基づいて求めることができる。すなわち、
100・・・・・・・スキャンガントリ部
101・・・・・・・X線源
102・・・・・・・回転盤
106・・・・・・・X線検出器
110・・・・・・・X線制御装置
120・・・・・・・操作卓
121・・・・・・・入力装置
122・・・・・・・画像処理装置
123・・・・・・・記憶装置
124・・・・・・・システム制御装置
125・・・・・・・表示装置
21・・・・・・・・焦点シフト投影データ取得部
22・・・・・・・・焦点別投影データ生成部
23・・・・・・・・仮想ビューデータ生成部
24・・・・・・・・投影データ合成部
25・・・・・・・・再構成演算部
31・・・・・・・・焦点シフト投影データ
32A,32B・・・焦点別投影データ
33A,33B・・・仮想投影データ
34・・・・・・・・空間分解能向上投影データ
4・・・・・・・・・仮想検出器
41・・・・・・・・焦点数と空間分解能との関係を示すデータ
51・・・・・・・・入力部
52・・・・・・・・仮想焦点数算出部52
53・・・・・・・・仮想焦点設定部
54・・・・・・・・仮想投影データ生成部54
55・・・・・・・・投影データ合成部
A,B・・・・・・・焦点
P・・・・・・・・・仮想焦点
Claims (12)
- 複数の焦点位置から被検体にX線を照射するX線源と、
前記X線源に対向配置され、前記被検体を透過したX線である透過X線を検出するX線検出器と、
前記X線源及び前記X線検出器を搭載し、前記被検体の周囲を回転する回転盤と、
前記X線源における前記焦点位置を複数箇所にシフトさせてX線を照射するX線制御部と、
前記X線制御部により前記焦点位置を複数箇所にシフトさせて照射した各X線により得た前記透過X線を収集し、焦点位置毎に分割した焦点別投影データを生成する焦点別投影データ生成部と、
一方の焦点位置についての焦点別投影データにおける欠損ビュー位置に仮想ビューを挿入し、当該仮想ビューデータを他方の焦点位置についての焦点別投影データを用い、かつ前記焦点位置と前記X線検出器との位置関係に基づいて算出する仮想ビューデータ生成部と、
前記仮想ビューデータ生成部により仮想ビューが挿入された各焦点別投影データを用いて画像を再構成する再構成演算部と、
を備えることを特徴とするX線CT装置。 - 前記仮想ビューデータ生成部は、
前記一方の焦点位置についての焦点別投影データと前記他方の焦点位置についての焦点別投影データとを、前記仮想ビューとの間の角度に応じた重みをつけて補間演算することにより前記仮想ビューデータを算出することを特徴とする請求項2に記載のX線CT装置。 - 前記仮想ビューデータ生成部により仮想ビューが挿入された各焦点位置の焦点別投影データを組み合わせる投影データ合成部を更に備え、
前記再構成演算部は、前記投影データ合成部により生成された投影データを画像を再構成することを特徴とする請求項2または請求項3に記載のX線CT装置。 - 前記再構成演算部は、前記仮想ビューデータ生成部により生成された仮想投影データと前記焦点別投影データとを用いて再構成面の画像を作成する3次元再構成処理を適用することを特徴とする請求項6に記載のX線CT装置。
- X線CT装置により撮影された実投影データを取得する実投影データ取得部と、
前記撮影時の焦点位置と異なる位置に仮想焦点を設定し、設定した仮想焦点からX線を照射した場合に得られる仮想投影データを、前記実投影データを用い、かつ前記焦点と前記仮想焦点と前記X線検出器との位置関係に基づいて算出する仮想投影データ生成部と、
前記実投影データと前記仮想投影データ生成部により生成された仮想投影データとを組み合わせてビュー数を向上させた投影データを生成する投影データ合成部と、
を備えることを特徴とする画像処理装置。 - 所望の空間分解能の目標値を入力する入力部と、
前記入力部により入力された空間分解能の目標値に応じて必要な仮想焦点数を算出する仮想焦点数算出部と、を更に備え、
前記仮想投影データ生成部は、前記仮想焦点を前記仮想焦点数算出部により算出された数だけ設定し、設定した仮想焦点からX線を照射した場合に得られる仮想投影データを、実在の焦点から得られた前記焦点別投影データを用い、かつ前記焦点と前記X線検出器と前記回転盤との位置関係に基づいて算出することを特徴とする請求項8に記載の画像処理装置。 - 前記投影データ取得部により取得した投影データの前記撮影時の実際の焦点位置が理想の焦点位置にない場合において、前記理想の焦点位置を前記仮想焦点とみなし、実際の焦点位置と前記理想の焦点位置とX線検出器位置との位置関係に基づいて、前記実際の焦点位置で得た投影データを前記理想の焦点位置において得た場合の投影データに補正するエラー焦点補正部を更に備えることを特徴とする請求項8に記載の画像処理装置。
- 画像処理装置が、
複数の焦点位置から被検体に対してX線を照射するX線CT装置により撮影された投影データを取得するステップと、
前記投影データを焦点位置毎に分割し、焦点別投影データを生成するステップと、
一方の焦点位置についての焦点別投影データにおける欠損ビュー位置に仮想ビューを挿入し、当該仮想ビューのデータを他方の焦点位置についての焦点別投影データを用い、かつ前記焦点位置とX線検出器との位置関係に基づいて算出するステップと、
を含むことを特徴とする投影データ生成方法。 - 画像処理装置が、
X線CT装置により撮影された実投影データを取得するステップと、
前記撮影時の焦点位置と異なる位置に仮想焦点を設定し、設定した仮想焦点からX線を照射した場合に得られる仮想投影データを、前記実投影データを用い、かつ前記焦点と前記仮想焦点と前記X線検出器との位置関係に基づいて算出するステップと、
前記実投影データと前記仮想投影データとを組み合わせビュー数を向上させた投影データを生成するステップと、
を含むことを特徴とする投影データ生成方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014093627A JP2015208601A (ja) | 2014-04-30 | 2014-04-30 | X線ct装置、画像処理装置、及び投影データ生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014093627A JP2015208601A (ja) | 2014-04-30 | 2014-04-30 | X線ct装置、画像処理装置、及び投影データ生成方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2015208601A true JP2015208601A (ja) | 2015-11-24 |
Family
ID=54611298
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2014093627A Pending JP2015208601A (ja) | 2014-04-30 | 2014-04-30 | X線ct装置、画像処理装置、及び投影データ生成方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2015208601A (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018004085A1 (ko) * | 2016-06-29 | 2018-01-04 | 한국생산기술연구원 | 비회전오블리크타입 ct장치 및 이를 이용한 피검사체의 3차원내부형상이미지 복원방법 |
CN112472109A (zh) * | 2019-09-12 | 2021-03-12 | 通用电气精准医疗有限责任公司 | 校准方法及x射线成像系统 |
US10945683B2 (en) | 2016-11-15 | 2021-03-16 | Canon Medical Systems Corporation | X-ray computed tomography apparatus |
WO2023132017A1 (ja) * | 2022-01-05 | 2023-07-13 | 株式会社島津製作所 | X線撮影装置およびx線撮影方法 |
CN116483025A (zh) * | 2023-04-23 | 2023-07-25 | 赛诺威盛科技(北京)股份有限公司 | 飞焦点模式下的数据采集系统、方法、电子设备及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS60179042A (ja) * | 1984-02-27 | 1985-09-12 | 株式会社日立メデイコ | Ct装置 |
JPH0810251A (ja) * | 1994-06-28 | 1996-01-16 | Hitachi Medical Corp | X線断層撮影方法および装置 |
US20120014503A1 (en) * | 2010-07-14 | 2012-01-19 | Christer Ullberg | Computed Tomography Scanning System and Method |
-
2014
- 2014-04-30 JP JP2014093627A patent/JP2015208601A/ja active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS60179042A (ja) * | 1984-02-27 | 1985-09-12 | 株式会社日立メデイコ | Ct装置 |
JPH0810251A (ja) * | 1994-06-28 | 1996-01-16 | Hitachi Medical Corp | X線断層撮影方法および装置 |
US20120014503A1 (en) * | 2010-07-14 | 2012-01-19 | Christer Ullberg | Computed Tomography Scanning System and Method |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018004085A1 (ko) * | 2016-06-29 | 2018-01-04 | 한국생산기술연구원 | 비회전오블리크타입 ct장치 및 이를 이용한 피검사체의 3차원내부형상이미지 복원방법 |
US10945683B2 (en) | 2016-11-15 | 2021-03-16 | Canon Medical Systems Corporation | X-ray computed tomography apparatus |
CN112472109A (zh) * | 2019-09-12 | 2021-03-12 | 通用电气精准医疗有限责任公司 | 校准方法及x射线成像系统 |
WO2023132017A1 (ja) * | 2022-01-05 | 2023-07-13 | 株式会社島津製作所 | X線撮影装置およびx線撮影方法 |
CN116483025A (zh) * | 2023-04-23 | 2023-07-25 | 赛诺威盛科技(北京)股份有限公司 | 飞焦点模式下的数据采集系统、方法、电子设备及介质 |
CN116483025B (zh) * | 2023-04-23 | 2024-03-22 | 赛诺威盛科技(北京)股份有限公司 | 飞焦点模式下的数据采集系统、方法、电子设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6246936B2 (ja) | X線撮像装置および画像再構成方法 | |
JPH08299322A (ja) | Ct装置 | |
JP5535532B2 (ja) | X線ct装置 | |
JP5687618B2 (ja) | コンピュータ断層撮像用スキャナ及びスキャン方法 | |
US7583782B2 (en) | X-ray CT apparatus and an image controlling method thereof | |
JP2006068523A (ja) | 検査対象の断層撮影によるスライス画像の作成方法およびコンピュータ断層撮影装置 | |
WO2012049940A1 (ja) | 医用画像処理装置、x線コンピュータ断層撮影装置および医用画像処理方法 | |
WO2013081062A1 (ja) | X線コンピュータ断層撮影装置 | |
JP2015208601A (ja) | X線ct装置、画像処理装置、及び投影データ生成方法 | |
JP2013000479A (ja) | X線ct装置及び画像再構成方法 | |
CN102770077A (zh) | X射线摄像装置 | |
JP2008012206A (ja) | X線断層撮影装置 | |
JP2021010727A (ja) | X線ctシステム及び医用処理装置 | |
JP2002172112A (ja) | 高速撮影用x線ctスキャナ | |
JP4041040B2 (ja) | 放射線断層撮影装置 | |
JP2006187453A (ja) | X線ct装置 | |
JP6377615B2 (ja) | X線ct装置及び画像再構成方法 | |
JP4025530B2 (ja) | X線ct装置 | |
JP2015116408A (ja) | X線ct装置及び不良素子補正方法 | |
JPS62284250A (ja) | 産業用ctスキヤナ | |
JP4027954B2 (ja) | コンピュータ断層撮影装置 | |
JP5595724B2 (ja) | X線ct装置 | |
JP3607285B2 (ja) | コンピュータ断層撮影装置 | |
JP2013150780A (ja) | 放射線撮影装置及び放射線撮影方法 | |
JP6220599B2 (ja) | X線ct装置及び投影データのアップサンプリング方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20160509 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20170410 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20171030 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20171030 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20180125 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20180130 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20180807 |