JPS59149555A - Ｘ線ｃｔ装置 - Google Patents
Ｘ線ｃｔ装置Info
- Publication number
- JPS59149555A JPS59149555A JP58022142A JP2214283A JPS59149555A JP S59149555 A JPS59149555 A JP S59149555A JP 58022142 A JP58022142 A JP 58022142A JP 2214283 A JP2214283 A JP 2214283A JP S59149555 A JPS59149555 A JP S59149555A
- Authority
- JP
- Japan
- Prior art keywords
- image
- data
- projection data
- picture
- dimensional
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔発明の技術分野〕
この発明は、画像たとえばＸ線ＣＴ装置により得られる
画像全補正する画像再処理装置に関する。
画像全補正する画像再処理装置に関する。
従来、画像診断装置たとえばＸ線ＣＴ装置で用いられる
Ｘａは、通常、多色Ｘ線であるため、Ｘ線エネルギによ
り被検体中にある骨や軟組織等の谷部位のＸ線吸収効率
が相違するので、線質硬化が生じる。線質硬化は、Ｘｌ
ｆＭＣＴ装置においてプロジェクションデータにより再
構成された画像にアーチ７アクトを生じさせる。このア
ーチ７アクトを除去するために、従来では、セカンドパ
ス線質硬化補正処理（２ｎｔｔ ｐａｓｓ ｂｅａｍ
ｈＬ：ＬｒｃｔｅｎｉｎｇＣＯγｒｔｃｔｉｏｎ ）が
行なわれている。
Ｘａは、通常、多色Ｘ線であるため、Ｘ線エネルギによ
り被検体中にある骨や軟組織等の谷部位のＸ線吸収効率
が相違するので、線質硬化が生じる。線質硬化は、Ｘｌ
ｆＭＣＴ装置においてプロジェクションデータにより再
構成された画像にアーチ７アクトを生じさせる。このア
ーチ７アクトを除去するために、従来では、セカンドパ
ス線質硬化補正処理（２ｎｔｔ ｐａｓｓ ｂｅａｍ
ｈＬ：ＬｒｃｔｅｎｉｎｇＣＯγｒｔｃｔｉｏｎ ）が
行なわれている。
Ｘ線ＣＴ装置におけるセカンドバス線質硬化補正処理の
一例を次に述べる。
一例を次に述べる。
■被検体の体軸を中心として被検体の周囲にＸ線を回転
させながら、Ｘ線を各回転角度毎に曝射して得られる全
てのプロジェクションデータを基に被検体の断層像を先
ず再構成する。
させながら、Ｘ線を各回転角度毎に曝射して得られる全
てのプロジェクションデータを基に被検体の断層像を先
ず再構成する。
■次いで、前記断層像からＢ（θ、ｔ）ｂおよびＳ（θ
Ｇ’）請求める。ここで、Ｂ（θ、１）は第１図に示す
ように、たとえば骨２と軟組織４との断層像において、
ｘｙ座標におけるｙ軸と角度θをもって交差すると共に
原点Ｏより距離ｔｆ有スるビーム（θ９１）上における
骨の総量（ｙ／龍）ヲ示すパラメータであり、Ｓ（θ＋
’）は、前記ビーム方向における軟組織のａ量（ｙ／ｍ
）を示すパラメータである。
Ｇ’）請求める。ここで、Ｂ（θ、１）は第１図に示す
ように、たとえば骨２と軟組織４との断層像において、
ｘｙ座標におけるｙ軸と角度θをもって交差すると共に
原点Ｏより距離ｔｆ有スるビーム（θ９１）上における
骨の総量（ｙ／龍）ヲ示すパラメータであり、Ｓ（θ＋
’）は、前記ビーム方向における軟組織のａ量（ｙ／ｍ
）を示すパラメータである。
■次いでビーム（θ、１）において、線質硬化の大きさ
Ｃ（θ、ｔ）’ｔ、表あるいは関数にょジ求める。
Ｃ（θ、ｔ）’ｔ、表あるいは関数にょジ求める。
■次いで線質硬化の大きさＣ（θ、１）をパラレル投影
データとして、これを再構成することにより補正用画像
を得、前記■で求めた画像からこの補正用画像を減算す
ることにより、線質硬化による影響を除去した画像を得
る。
データとして、これを再構成することにより補正用画像
を得、前記■で求めた画像からこの補正用画像を減算す
ることにより、線質硬化による影響を除去した画像を得
る。
前記Ｑ）の処理は、第（１）式により行なうことができ
る。
る。
Ｃ（θ＊ｔ）Ｆ！ＴＣＢ＜θ９ｔ）、、５（θ、 ｔ
） ）・（１）ここで、関数Ｔは、第（２）式に示すよ
うに近似することができる。
） ）・（１）ここで、関数Ｔは、第（２）式に示すよ
うに近似することができる。
Ｔ（Ｂ 、 Ｓ ） −Ｒ１（Ｊ’ ） ＋ Ｒ２（Ｊ
’す＋・−＋ Ｒｎ （Ｖ”） ・（２）ただし、Ｒ
ｈ（ｒ−は１からルまでの整数）は関数であり、またｒ
は第（３）式により表わすことができる。
’す＋・−＋ Ｒｎ （Ｖ”） ・（２）ただし、Ｒ
ｈ（ｒ−は１からルまでの整数）は関数であり、またｒ
は第（３）式により表わすことができる。
Ｖ＝ｋＢＢ＋ｋｓＳ 、・曲（３）
ただし、ｋＢ、ｋＢは定数である。
したがって、第（１）式のＣ（θ９ｔ）は、第（４）式
に置きかえることができる。
に置きかえることができる。
Ｃ（θ、 ｔ ）−！Ｒ１（Ｖ （θ、 ｔ ） ）十
Ｒ２Ｃ（Ｖ（θ、す）り十・・・・・・十Ｒｎ （（Ｖ
（θ、υ）ｎ）・・・（４）第（４）式における〆（
θ、１）は、前記■の処理により求めることもできるが
、次の処理を行なってもよい。つまハ断層像上の各ビク
セルｐ（、ｒ：。
Ｒ２Ｃ（Ｖ（θ、す）り十・・・・・・十Ｒｎ （（Ｖ
（θ、υ）ｎ）・・・（４）第（４）式における〆（
θ、１）は、前記■の処理により求めることもできるが
、次の処理を行なってもよい。つまハ断層像上の各ビク
セルｐ（、ｒ：。
ｙ）のＣＴ値μψンについて
もしμ（Ｐ）〉骨についての基準値、であれば。
Ｘ＝μψ）ＸＡＢで表わされるＸは骨を示すビクセルで
あり１ もしμ艶）〉軟組織につりでの基準値、であれば、 Ｘ＝μ（ｐ）ｘｔｃＢ で表わされるＸは軟組織を示す
ビクセルであり、 Ｘが前記以外のときはＸ＝Ｑとし、そのビクセルは空気
を表示する、 との一連のアルゴリズムによりγ（ｘ、ｙ）ｋ求める。
あり１ もしμ艶）〉軟組織につりでの基準値、であれば、 Ｘ＝μ（ｐ）ｘｔｃＢ で表わされるＸは軟組織を示す
ビクセルであり、 Ｘが前記以外のときはＸ＝Ｑとし、そのビクセルは空気
を表示する、 との一連のアルゴリズムによりγ（ｘ、ｙ）ｋ求める。
このγＣｘｅｙ）は、一枚の画像であり、しかも、γ（
ｘ、ｙ）は画像μ（” ｅ ｙ ） Ｋ対して画像処理
を行なって得たものである。したがって、Ｖ（θｖ’）
は、ｒ （” ＊ ｙ ） ＶＣより第（５）式のよう
に表示することができる。
ｘ、ｙ）は画像μ（” ｅ ｙ ） Ｋ対して画像処理
を行なって得たものである。したがって、Ｖ（θｖ’）
は、ｒ （” ＊ ｙ ） ＶＣより第（５）式のよう
に表示することができる。
×ｄｈ ・・・・・・・・・（５）
ここで、ｈはビーム（θＤｔ）におけるビーム上での位
置を示す変数である。この第（５）式は、ｒ（θｓ’）
が画像ｒ （” ｏ ｙ ）のパラレル投影データであ
ることを示している。結局、Ｖ（θ。
置を示す変数である。この第（５）式は、ｒ（θｓ’）
が画像ｒ （” ｏ ｙ ）のパラレル投影データであ
ることを示している。結局、Ｖ（θ。
ｔ）を求めるためには、まずγ（” ｏ ｙ ）を求め
、次いでγ（” ９ ｙ ）の投影データを求めればよ
い。
、次いでγ（” ９ ｙ ）の投影データを求めればよ
い。
しかしながらｔ γ（ｘ、ｙ）の投影データＶ（θ９ｔ
）を第（５）式にしたがって求めようとすると、計算量
が膨大となり、計算に時間がかかるとの問題点がある。
）を第（５）式にしたがって求めようとすると、計算量
が膨大となり、計算に時間がかかるとの問題点がある。
この発明は前記事情に鑑みてなされたものであり、画像
Ｔから投影データＶ’ｚ迅速に算出することのできる画
像再処理装置を提供することを目的とするものである。
Ｔから投影データＶ’ｚ迅速に算出することのできる画
像再処理装置を提供することを目的とするものである。
前記目的を達成するためのこの発明の概要は、画像処理
する画像再処理装置において、既に処理されて得られる
ところの画像を構成する各画像データを２次元フーリエ
変換し、次いで極座標変換した後、１次元逆フーリエ変
換するプロセッサを有することを特徴とするものである
。
する画像再処理装置において、既に処理されて得られる
ところの画像を構成する各画像データを２次元フーリエ
変換し、次いで極座標変換した後、１次元逆フーリエ変
換するプロセッサを有することを特徴とするものである
。
先ずこの発明についての原理的説明を行なう。
この発明における画像処理装置においては、得られる原
画像から抽出したところの、誤差を生じさせる部位の抽
出画像から以下の手順により抽出画像についての投影デ
ータを求めることをその特長とする。
画像から抽出したところの、誤差を生じさせる部位の抽
出画像から以下の手順により抽出画像についての投影デ
ータを求めることをその特長とする。
一般に、投影データＶから画像γを求める処理は画像再
構成処理と称され、投影データＶが平行ビームにより求
められるのであれば、第（６）式に示す２次元フーリエ
変換法により投影データＶから画像γ全再構成すること
ができる。
構成処理と称され、投影データＶが平行ビームにより求
められるのであれば、第（６）式に示す２次元フーリエ
変換法により投影データＶから画像γ全再構成すること
ができる。
ｒ −Ｆ−２（Ｐ−１（ＦＣＶ））） ・・・・・・
（６）ここで、Ｆは、ｒ〔以下、ｒ（θ、１）と示すこ
とがある。〕’ｒｔＶ？一ついて７−リエ変換してＶ（
θ、ω）を求める処理であり、υ（θ、ω〕は極座標で
示す画像である。ｐ−１は、ＩＩ＋（θ９ω）を直交座
標で示す画像Ｆ（ξ９η）１求める処理であり ｐ’−
２は、画像Ｆ（ξ９η）をξとηとについて２次元逆フ
ーリエ変換して画像γ（、ｒ、ｙ）を求める処理である
。
（６）ここで、Ｆは、ｒ〔以下、ｒ（θ、１）と示すこ
とがある。〕’ｒｔＶ？一ついて７−リエ変換してＶ（
θ、ω）を求める処理であり、υ（θ、ω〕は極座標で
示す画像である。ｐ−１は、ＩＩ＋（θ９ω）を直交座
標で示す画像Ｆ（ξ９η）１求める処理であり ｐ’−
２は、画像Ｆ（ξ９η）をξとηとについて２次元逆フ
ーリエ変換して画像γ（、ｒ、ｙ）を求める処理である
。
したがって、抽出画像をγとし、その投影データをＶと
すると、抽出画像γから第（力式にしたがってその投影
データを求めることができる。
すると、抽出画像γから第（力式にしたがってその投影
データを求めることができる。
Ｖ、＝Ｆ″’Ｉ（７’（Ｆ２（γ））） ・・・・・
・ （力第（７）式に示す処理手順は、抽出画像γ（ｘ
、ｙ）を２次元フーリエ変換処理Ｆ２シた後、極座標変
換処理Ｐｋし、次いで、１次元逆フーリエ変換処理Ｆ＝
１ｙするものである。２次元フーリエ変換処理Ｆ２およ
び１次元逆フーリエ変換処理Ｆ−１は、高速フーリエ変
換処理（ＦＦＴ）ＶＣより迅速に行なうことができる。
・ （力第（７）式に示す処理手順は、抽出画像γ（ｘ
、ｙ）を２次元フーリエ変換処理Ｆ２シた後、極座標変
換処理Ｐｋし、次いで、１次元逆フーリエ変換処理Ｆ＝
１ｙするものである。２次元フーリエ変換処理Ｆ２およ
び１次元逆フーリエ変換処理Ｆ−１は、高速フーリエ変
換処理（ＦＦＴ）ＶＣより迅速に行なうことができる。
抽出画像から投影データが求まれば、投影データ上の誤
差量も簡単に求−！り、原画像の誤差を除去することが
できる。
差量も簡単に求−！り、原画像の誤差を除去することが
できる。
次に、この発明の前記原理全適用するこの発明の一実施
例である画像再処理装置を利用してＸ線ＣＴ装置により
得られる断層像における線質硬化を補正する画像補正装
置について説明する。
例である画像再処理装置を利用してＸ線ＣＴ装置により
得られる断層像における線質硬化を補正する画像補正装
置について説明する。
第２図はこの発明の一実施例である画像再処理装置を有
する画像補正装置である。
する画像補正装置である。
第２図に示すように５画像補正装置１０は、Ｘ線ＣＴ装
置に、おいて投影データを再構成して得た原画像につい
てのデジタル画像データμＣ，Ｔ、／）を補正して線質
硬化による誤差を除去した画像についてのデジタル画像
データμ’（ｘ、ｙ）ｋ出力するものであｐ１第１の画
像メモリ１２９画像処理ユニット１４．第２の画像メモ
リ１６、画像再処理装置１８、投影データメモリ２０、
線質硬化量（ＢＥ ）推定ユニット２２、テーブル２４
、再構成ユニット２６、および演算ユニット２８を有す
る。
置に、おいて投影データを再構成して得た原画像につい
てのデジタル画像データμＣ，Ｔ、／）を補正して線質
硬化による誤差を除去した画像についてのデジタル画像
データμ’（ｘ、ｙ）ｋ出力するものであｐ１第１の画
像メモリ１２９画像処理ユニット１４．第２の画像メモ
リ１６、画像再処理装置１８、投影データメモリ２０、
線質硬化量（ＢＥ ）推定ユニット２２、テーブル２４
、再構成ユニット２６、および演算ユニット２８を有す
る。
第１の画像メモリ１２は、デジタル画像データ像処理ユ
ニット１４および演算ユニット２８に出し、また、演算
ユニッ）２８により得られた線質硬化による誤差全除去
したデジタル画像データμ′ＣＬ−，３／）を入力し、
一旦これμ’（ｘ、 ｙ）全格納した後、デジタル画像
データμ′（ｘ、ｙ）ｋたとえば画像表示装置に出力す
る記憶装置である。
ニット１４および演算ユニット２８に出し、また、演算
ユニッ）２８により得られた線質硬化による誤差全除去
したデジタル画像データμ′ＣＬ−，３／）を入力し、
一旦これμ’（ｘ、 ｙ）全格納した後、デジタル画像
データμ′（ｘ、ｙ）ｋたとえば画像表示装置に出力す
る記憶装置である。
画像処理ユニット１４は、第１の画像メモリ１２に格納
したデジタル画像データμ（ｘ、ｙ）を読み出し、デジ
タル画像データμ（−ｒ、ｙ）Ｋより示される画像中の
一定値以上のＣＴ値を与える部位の画像全抽出し、抽出
した画像についてのデジタル画像データγ（ｘ、ｙ）ｆ
出力する画像プロセッサである。
したデジタル画像データμ（ｘ、ｙ）を読み出し、デジ
タル画像データμ（−ｒ、ｙ）Ｋより示される画像中の
一定値以上のＣＴ値を与える部位の画像全抽出し、抽出
した画像についてのデジタル画像データγ（ｘ、ｙ）ｆ
出力する画像プロセッサである。
第２の画像メモリ１６は、画像処理ユニット１４より出
力されるデジタル画像データγ（ａ−、ｙ）を一旦格納
した後、これ全画像再処理装置１８に出力し、また、再
構成ユニット２６より出力されるデジタル画像データＥ
Ｃｘ、／）を一旦格納した後、これを演算ユニット２８
に出力する記憶装置である。
力されるデジタル画像データγ（ａ−、ｙ）を一旦格納
した後、これ全画像再処理装置１８に出力し、また、再
構成ユニット２６より出力されるデジタル画像データＥ
Ｃｘ、／）を一旦格納した後、これを演算ユニット２８
に出力する記憶装置である。
画像再処理装置１８は、第２の画像メモリ１６より読み
出したデジタル画像データγ（ｒ、ｙ）を高速フーリエ
変換器（ＦＦＴ）Ｋより２次元フーリエ変換処理１２μ
得られるデータＦ（ξ９η）を極座標変換処理Ｐし、次
いで得られるデータＶ（θ、ω）を高速フーリエ変換器
（ＦＦＴ）により１次元逆フーリエ変換処理Ｆ−’ Ｘ
＼薦１して投影データＶ（θ、ｔ）ＶＣするデータプロ
セッサである。
出したデジタル画像データγ（ｒ、ｙ）を高速フーリエ
変換器（ＦＦＴ）Ｋより２次元フーリエ変換処理１２μ
得られるデータＦ（ξ９η）を極座標変換処理Ｐし、次
いで得られるデータＶ（θ、ω）を高速フーリエ変換器
（ＦＦＴ）により１次元逆フーリエ変換処理Ｆ−’ Ｘ
＼薦１して投影データＶ（θ、ｔ）ＶＣするデータプロ
セッサである。
投影データメモリ２０は、投影データＶ（θｔ’）を一
旦格納した後、これをＢＥ推定ユニット２２に出力し、
また、Ｂｌｌ推定ユニット２２より出力されるパラレル
投影データＣ（θｓ’）を一旦格納した後、これを再構
成ユニット２６に出力する記憶装置である。
旦格納した後、これをＢＥ推定ユニット２２に出力し、
また、Ｂｌｌ推定ユニット２２より出力されるパラレル
投影データＣ（θｓ’）を一旦格納した後、これを再構
成ユニット２６に出力する記憶装置である。
ＢＨ推定ユニット２２は、人力する投影データＶ（θ、
１）を基に、テーブル２４に格納されている表あるいは
関数ＲＯｔ Ｒ１、・・・９Ｒｎを参照して線質硬化に
よって生ずる誤差の大きさＢＨを求め、これをパラレル
投影データＣ（θｓ ”）として出力するデータプロセ
ッサである。
１）を基に、テーブル２４に格納されている表あるいは
関数ＲＯｔ Ｒ１、・・・９Ｒｎを参照して線質硬化に
よって生ずる誤差の大きさＢＨを求め、これをパラレル
投影データＣ（θｓ ”）として出力するデータプロセ
ッサである。
再構成ユニット２６は、全てのθについてのパラレル投
影データＣ（θｏ’）を基に画像再構成し、補正用画像
についてのデジタル画像データＥ（：ｔ 、 ｙ）を第
２の画像メモ！Ｊ １６ Ｋ出力するデータプロセッサ
である。
影データＣ（θｏ’）を基に画像再構成し、補正用画像
についてのデジタル画像データＥ（：ｔ 、 ｙ）を第
２の画像メモ！Ｊ １６ Ｋ出力するデータプロセッサ
である。
演算ユニット２８は、第１の画像メモＩＪ １２ Ｋ格
納されているデジタル画像データμＣｘ 、ｙ）から、
第２の画像メモリ１６より読み出した補正用画像につい
てのデジタル画像データＥ （ｘ 、 ｙ ）を減算す
ることによＶ線質硬化の補正がさ扛た画像についてのデ
ジタル画像データμ′（！、ｙ）を出力するデータプロ
セッサである。
納されているデジタル画像データμＣｘ 、ｙ）から、
第２の画像メモリ１６より読み出した補正用画像につい
てのデジタル画像データＥ （ｘ 、 ｙ ）を減算す
ることによＶ線質硬化の補正がさ扛た画像についてのデ
ジタル画像データμ′（！、ｙ）を出力するデータプロ
セッサである。
次に、以上構成の画像処理装置１０の作用について述べ
る。
る。
図示しない）（ｌｌＪＣＴ装置において得られるゾロジ
エクションデータを再構成して得た原画像たとえば骨と
軟組織を表示する断層像についてのデジタル画像データ
μ（ｘ ？ ｙ ）を一旦第１の画像メモリ１２ＶＣ格
納した後、前記デジタル画像データμ（” ｃ ｙ ）
を画像処理ユニット１４および演算ユニット２８に出力
する。
エクションデータを再構成して得た原画像たとえば骨と
軟組織を表示する断層像についてのデジタル画像データ
μ（ｘ ？ ｙ ）を一旦第１の画像メモリ１２ＶＣ格
納した後、前記デジタル画像データμ（” ｃ ｙ ）
を画像処理ユニット１４および演算ユニット２８に出力
する。
画像処理ユニット１４では、入力するデジタル画像デー
タμ（ｘｏ ｙ ）により示される画像中の一定値以上
のＣＴ値を与える部位たとえば骨の断層像を抽出し、抽
出した断層像についてのデジタル画像データγ（：ｔ。
タμ（ｘｏ ｙ ）により示される画像中の一定値以上
のＣＴ値を与える部位たとえば骨の断層像を抽出し、抽
出した断層像についてのデジタル画像データγ（：ｔ。
ｙ）を出力する。デジタル画像データγＣｘ、ｙ）は、
第２の画像メモリ１６を介して、画像再処理装置１８に
出力され、画像再、処理装置１８において、高速７−リ
エ変換アルゴリズム（ＦＦＴ）ＶＣより２次元フーリエ
変換処理Ｆ２されてデジタル画像データｒ（４′、η）
となり、次いでデジタル画像データＦ（ξ、η〕は極座
標変換処理Ｐされてデジタル画像データν（θ、ω）と
なり、次いでデジタル画像データＶ（θ、ω）は高速フ
ーリエ変換アルゴリズム（ＦＦＴ）によ！ｌ１１次元逆
フーリエ変換されて投影データＶ（θｏ’）となる。投
影データｒ（θ、ｔ）は、投影データメモリ２０を介し
てＢＥ推足ユニツ）２．２に出力される。ＢＪＩ推定ユ
ニット２２は、テーブル２４に格納されているたとえば
関数Ｒ８＋Ｒ１ｏ・・・、 Ｒｎを参照し、前記第４式
の演算を行なって、投影データ〆（θｏ’）からパラレ
ル投影データＣ′（θ、１）を求め、得られるパラレル
投影データＣ（θ、１）を、投影データメモリ２０を介
して再構成ユニット２６に出力する。
第２の画像メモリ１６を介して、画像再処理装置１８に
出力され、画像再、処理装置１８において、高速７−リ
エ変換アルゴリズム（ＦＦＴ）ＶＣより２次元フーリエ
変換処理Ｆ２されてデジタル画像データｒ（４′、η）
となり、次いでデジタル画像データＦ（ξ、η〕は極座
標変換処理Ｐされてデジタル画像データν（θ、ω）と
なり、次いでデジタル画像データＶ（θ、ω）は高速フ
ーリエ変換アルゴリズム（ＦＦＴ）によ！ｌ１１次元逆
フーリエ変換されて投影データＶ（θｏ’）となる。投
影データｒ（θ、ｔ）は、投影データメモリ２０を介し
てＢＥ推足ユニツ）２．２に出力される。ＢＪＩ推定ユ
ニット２２は、テーブル２４に格納されているたとえば
関数Ｒ８＋Ｒ１ｏ・・・、 Ｒｎを参照し、前記第４式
の演算を行なって、投影データ〆（θｏ’）からパラレ
ル投影データＣ′（θ、１）を求め、得られるパラレル
投影データＣ（θ、１）を、投影データメモリ２０を介
して再構成ユニット２６に出力する。
再構成ユニット２６は、人力する全てのθについてのパ
ラレル投影データＣ（θ＋’）を再構成して、補正用画
像についてのデジタル画像データＥ（ｒ、、９）に変換
し、デジタル画像データＥ （ｘ。
ラレル投影データＣ（θ＋’）を再構成して、補正用画
像についてのデジタル画像データＥ（ｒ、、９）に変換
し、デジタル画像データＥ （ｘ。
ｙ）を、第２の画像メモリ１６を介して演算ユニット２
８に出力する。演算ユニット２８は、デジタル画像デー
タμ（ｘ、ｙ）からデジタル画像データＥ＜ｘ、ｙ）を
減算することにより線質硬化の補正をし、線質硬化のな
い断層像についてのデジタル画像データμ′（ｘ、ｙ）
を出力する。
８に出力する。演算ユニット２８は、デジタル画像デー
タμ（ｘ、ｙ）からデジタル画像データＥ＜ｘ、ｙ）を
減算することにより線質硬化の補正をし、線質硬化のな
い断層像についてのデジタル画像データμ′（ｘ、ｙ）
を出力する。
なお、前記デジタル画像データμ′（ｘ、ｙ）は、Ｘ線
ＣＴ装置における画像表示部に出力され、ビデオ信号に
変換された後、線質硬化の除去されたアーチファクトの
ない画像が表示されることになる。
ＣＴ装置における画像表示部に出力され、ビデオ信号に
変換された後、線質硬化の除去されたアーチファクトの
ない画像が表示されることになる。
以上詳述した画像補正装置１０は、デジタル画像データ
γ（ｘ、ｙ）をデジタル画像データＶ（θ。ｔ）に高速
変換する画像再処理装置１８を有するので、高速で短時
間のうちに原画像についての線質硬化を補正することが
できる。
γ（ｘ、ｙ）をデジタル画像データＶ（θ。ｔ）に高速
変換する画像再処理装置１８を有するので、高速で短時
間のうちに原画像についての線質硬化を補正することが
できる。
以上、この発明の一実施例について詳述したが、この発
明は前記笑施例に限定さ扛るものではなく、この発明の
要旨を変更しない範囲内で適宜に変形して笑施すること
ができるのはいうまでもない。
明は前記笑施例に限定さ扛るものではなく、この発明の
要旨を変更しない範囲内で適宜に変形して笑施すること
ができるのはいうまでもない。
この発明に係る画像再処理装置は、前記のような線質硬
化補正をする画像補正装＃に限らず、散乱線による誤差
を補正する画像補正装置にも利用することができる。
化補正をする画像補正装＃に限らず、散乱線による誤差
を補正する画像補正装置にも利用することができる。
捷た、画像再処理装置を組み込んだ画像補正装置は、第
２図に示すほかに、原画像を示すデジタル画像データμ
（、ｒ、ｙ＞から複数の画像それぞれについてのデジタ
ル画像データγ、（、：？７９３＋ ）、γ２（ｘ＋
ｙ ）、・・・、γ、Ｌ（、ｚ： 、 、９ ）を抽出
し、次いで前記デジタル画像データγＩＣＺ９ｙ）−γ
２（ｘ。
２図に示すほかに、原画像を示すデジタル画像データμ
（、ｒ、ｙ＞から複数の画像それぞれについてのデジタ
ル画像データγ、（、：？７９３＋ ）、γ２（ｘ＋
ｙ ）、・・・、γ、Ｌ（、ｚ： 、 、９ ）を抽出
し、次いで前記デジタル画像データγＩＣＺ９ｙ）−γ
２（ｘ。
ｙ）ｓ・・・、γｎ＜ｘ、ｙ）を画像再処理装置で複数
の投影データＶ、 （θ＋’）、Ｖ２（θ、１）、・・
・、Ｖｎ（θ、１）に再投影し、これら投影データＶ、
（θｅｔ）ｍ Ｖ２（θ、１）１．、、 Ｖｎ、（θＧ
’）を非線形合成して一つのパラレル投影データＣ（θ
、１）を得、この投影データＣ（θｏｔ）より画像再構
成して得た補正用画像についてのデジタル画像データＥ
（、ｔ： 、 ｙ ） ＶＣより原画像についてのデ
ジタル画像データμ（ｘ、ｙ）を補正するように、構成
してもよい。
の投影データＶ、 （θ＋’）、Ｖ２（θ、１）、・・
・、Ｖｎ（θ、１）に再投影し、これら投影データＶ、
（θｅｔ）ｍ Ｖ２（θ、１）１．、、 Ｖｎ、（θＧ
’）を非線形合成して一つのパラレル投影データＣ（θ
、１）を得、この投影データＣ（θｏｔ）より画像再構
成して得た補正用画像についてのデジタル画像データＥ
（、ｔ： 、 ｙ ） ＶＣより原画像についてのデ
ジタル画像データμ（ｘ、ｙ）を補正するように、構成
してもよい。
さらに、画像再処理装置を組み込んだ画像補正装置は、
第２図に示すほかに、原画像を示すデジタル画像データ
μ（”ｗｙ）から複数の画像それぞれについてのデジタ
ル画像データγ１ （” ｅ ’／ ） ｓγ２ （”
＋ ３’ハ・・・、ｒｒＬ（” ｗ ｙ ）を抽出し
、次いで前記デジタル画像データγｔ（”ｗｙ）、γ２
（ｘ −ｙ ）、・・・、γｎ（ｚ−、ｙ）を画像再
処理装置で複数の投影データ７（θ−ｔ）％ ’２（θ
、１）、・・・、ＶｒＬ（θ９ｔ）に再投影し、これら
投影データＶ、（θ９’）％ Ｖ２（θ。
第２図に示すほかに、原画像を示すデジタル画像データ
μ（”ｗｙ）から複数の画像それぞれについてのデジタ
ル画像データγ１ （” ｅ ’／ ） ｓγ２ （”
＋ ３’ハ・・・、ｒｒＬ（” ｗ ｙ ）を抽出し
、次いで前記デジタル画像データγｔ（”ｗｙ）、γ２
（ｘ −ｙ ）、・・・、γｎ（ｚ−、ｙ）を画像再
処理装置で複数の投影データ７（θ−ｔ）％ ’２（θ
、１）、・・・、ＶｒＬ（θ９ｔ）に再投影し、これら
投影データＶ、（θ９’）％ Ｖ２（θ。
ｔ）、・・・、Ｖｎ（θ、１）それぞれからパラレル投
影データ（？１（θ、１）、Ｃ２（θ。ｔ）、・・・、
Ｃｎ（θ９ｔ）を得、これら投影データＣ，（θ、１）
、（１’２（θｏｔ）、・・・、Ｃユ（θ９ｔ）それぞ
れにつき画像再構成して得た画像についてのデジタル画
像データＥｔ （ｘ 、ｙ ）、Ｅ２（ｘ９ｙ）、・・
・、ＥｒＬ（ｘ、ｙ）を非線形合成して補正用画像につ
いてのデジタル画像データＥ （ｘ 。
影データ（？１（θ、１）、Ｃ２（θ。ｔ）、・・・、
Ｃｎ（θ９ｔ）を得、これら投影データＣ，（θ、１）
、（１’２（θｏｔ）、・・・、Ｃユ（θ９ｔ）それぞ
れにつき画像再構成して得た画像についてのデジタル画
像データＥｔ （ｘ 、ｙ ）、Ｅ２（ｘ９ｙ）、・・
・、ＥｒＬ（ｘ、ｙ）を非線形合成して補正用画像につ
いてのデジタル画像データＥ （ｘ 。
ｙ）を得、このデジタル画像データＥ（ｘ、ｙ）を直接
ＶｃＸ線ＣＴ装置の画像表示部に出力し、あるいは、前
記デジタル画像データＥ （ｘ 、 ｙ ）により原像
についてのデジタル画像データμＣｘ、ｙ）を補正する
ように構成してもよい。
ＶｃＸ線ＣＴ装置の画像表示部に出力し、あるいは、前
記デジタル画像データＥ （ｘ 、 ｙ ）により原像
についてのデジタル画像データμＣｘ、ｙ）を補正する
ように構成してもよい。
また、第６図に示すよ５ ｖｃ ｓ第２図における第１
の画像メモリ１２５画像処理ユニット１４、第２の画像
メモリ１６、投影データメモリ２０、再構成ユニット２
６および演算ユニット２８を、Ｘ線ＣＴ装置における汎
用および再構成プロセッサ６０内に含ませ、画像補正装
置を、画像再処理装置１Ｂ、ＢＨ推定ユニット２２、お
よびテーフ゛ル２４により構成してもよい。なお、第６
図において、６１で示すのはＸ線を曝射することにより
得られる仮構体についての投影データを出力するガント
リであり、６２で示すのは画像データを記憶□する画像
メモリであり、６６で示すのは画像表示部である。
の画像メモリ１２５画像処理ユニット１４、第２の画像
メモリ１６、投影データメモリ２０、再構成ユニット２
６および演算ユニット２８を、Ｘ線ＣＴ装置における汎
用および再構成プロセッサ６０内に含ませ、画像補正装
置を、画像再処理装置１Ｂ、ＢＨ推定ユニット２２、お
よびテーフ゛ル２４により構成してもよい。なお、第６
図において、６１で示すのはＸ線を曝射することにより
得られる仮構体についての投影データを出力するガント
リであり、６２で示すのは画像データを記憶□する画像
メモリであり、６６で示すのは画像表示部である。
この発明によると、高速フーリエ変換アルコ゛リズムＦ
ＦＴを用いて高速で画像から投影データを求めること
のできる画像再処理装置を提供することができ、この発
明に係る画像再処理装置を、各種の画像を補正するため
の画像補正装置に組み込むと、高速で画像の補正をする
ことができる。
ＦＴを用いて高速で画像から投影データを求めること
のできる画像再処理装置を提供することができ、この発
明に係る画像再処理装置を、各種の画像を補正するため
の画像補正装置に組み込むと、高速で画像の補正をする
ことができる。
第１図は断層像中の特定ビーム方向（θｏ’）における
骨の総量についてのノくラメータＢ（θ。 ｔ）と軟組織の総量についてのノくラメータＳ（θ。 ｔ）とを示す説明図、第２図はこの発明の一実施例であ
る画像再処理装置を有する画像補正装置を示すブロック
図、および第６図はこの発明の一実施例である画像再処
理装置を有する他の画像補正装置を結合するＸ線ＣＴ装
置を示すブロック図である。 １８・・・画像再処理装置。
骨の総量についてのノくラメータＢ（θ。 ｔ）と軟組織の総量についてのノくラメータＳ（θ。 ｔ）とを示す説明図、第２図はこの発明の一実施例であ
る画像再処理装置を有する画像補正装置を示すブロック
図、および第６図はこの発明の一実施例である画像再処
理装置を有する他の画像補正装置を結合するＸ線ＣＴ装
置を示すブロック図である。 １８・・・画像再処理装置。
Claims (1)
- 画像処理する画像再処理装置において、既に処理されて
得られるところの画像を構成する各画像データを２次元
フーリエ変換し、次いで極座標変換した後、１次元逆フ
ーリエ変換するプロセッサを有することを特徴とする画
像再処理装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58022142A JPS59149555A (ja) | 1983-02-15 | 1983-02-15 | Ｘ線ｃｔ装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58022142A JPS59149555A (ja) | 1983-02-15 | 1983-02-15 | Ｘ線ｃｔ装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS59149555A true JPS59149555A (ja) | 1984-08-27 |
| JPH0230066B2 JPH0230066B2 (ja) | 1990-07-04 |
Family
ID=12074620
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP58022142A Granted JPS59149555A (ja) | 1983-02-15 | 1983-02-15 | Ｘ線ｃｔ装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS59149555A (ja) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6041170A (ja) * | 1983-06-07 | 1985-03-04 | エルシント リミテツド | 映像再投映システム |
| JPS6429981A (en) * | 1987-07-24 | 1989-01-31 | Matsushita Electric Works Ltd | Image processor |
| JPS6491278A (en) * | 1987-10-01 | 1989-04-10 | Toshiba Corp | System for reconstructing image |
| JPS6491279A (en) * | 1987-10-01 | 1989-04-10 | Toshiba Corp | System for reconstructing image |
| JPH01291388A (ja) * | 1988-05-18 | 1989-11-22 | Yokogawa Medical Syst Ltd | 物理量の２次元分布からストリークを識別する方法及び装置 |
| JP2009545409A (ja) * | 2006-08-03 | 2009-12-24 | ザ リージェンツ オブ ザ ユニバーシティ オブ カリフォルニア | トモグラフィーにおけるドーズ低減および画像改良のための反復方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5773461A (en) * | 1980-10-27 | 1982-05-08 | Toshiba Corp | Picture processing unit of computed tomography |
-
1983
- 1983-02-15 JP JP58022142A patent/JPS59149555A/ja active Granted
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5773461A (en) * | 1980-10-27 | 1982-05-08 | Toshiba Corp | Picture processing unit of computed tomography |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6041170A (ja) * | 1983-06-07 | 1985-03-04 | エルシント リミテツド | 映像再投映システム |
| JPS6429981A (en) * | 1987-07-24 | 1989-01-31 | Matsushita Electric Works Ltd | Image processor |
| JPS6491278A (en) * | 1987-10-01 | 1989-04-10 | Toshiba Corp | System for reconstructing image |
| JPS6491279A (en) * | 1987-10-01 | 1989-04-10 | Toshiba Corp | System for reconstructing image |
| JPH01291388A (ja) * | 1988-05-18 | 1989-11-22 | Yokogawa Medical Syst Ltd | 物理量の２次元分布からストリークを識別する方法及び装置 |
| JP2009545409A (ja) * | 2006-08-03 | 2009-12-24 | ザ リージェンツ オブ ザ ユニバーシティ オブ カリフォルニア | トモグラフィーにおけるドーズ低減および画像改良のための反復方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0230066B2 (ja) | 1990-07-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4841874B2 (ja) | 断層合成造影における直接再生方法及び装置 | |
| JP5749735B2 (ja) | Ｘ線放射線写真における骨抑制 | |
| US20050046644A1 (en) | 3D image processing apparatus | |
| JP6442243B2 (ja) | コンピュータ断層撮影画像再構成のアーチファクト低減の方法および装置 | |
| Dunn et al. | A comparison of two registration techniques for digital subtraction radiography | |
| JP2002191591A (ja) | Ｃｔ装置により取得した出力像に対する放射硬化の補正方法 | |
| JP4828920B2 (ja) | ３次元画像処理装置 | |
| BRPI1105671A2 (pt) | Aparelho e método de processamento de imagem | |
| JP2006150080A (ja) | 回転血管撮影のための血管撮影ｘ線診断装置 | |
| JPS63256843A (ja) | 角度制限ｘ線計算機式断層写真法において物体の境界情報を求める方法及び角度制限ｘ線計算機式断層写真法システムにおいて高品質の像を再生する方法 | |
| JP3987024B2 (ja) | 横方向のフィルタリング処理を用いたトモシンセシス画像を強調する方法及びシステム | |
| JPS59149555A (ja) | Ｘ線ｃｔ装置 | |
| CN111462267A (zh) | 一种x射线投影数据的获取方法和系统 | |
| JP2536024B2 (ja) | Ｘ線画像処理装置 | |
| JP2005040506A (ja) | Ｘ線画像処理方法及び装置 | |
| Montúfar et al. | Perspective and orthogonal CBCT/CT digitally reconstructed radiographs compared to conventional cephalograms | |
| KR102453897B1 (ko) | 컴퓨터 단층촬영 영상을 활용한 3차원 구강 스캔 데이터 복원 장치 및 방법 | |
| JPS62139630A (ja) | Ｘ線ｃｔ装置 | |
| Peters | Spatial filtering to improve transverse tomography | |
| US11786193B2 (en) | Metal artifacts reduction in cone beam reconstruction | |
| KR102248191B1 (ko) | 콘빔 컴퓨터 단층촬영 영상에서 콘빔아티팩트를 보정하는 방법 및 장치, 그것을 포함하는 콘빔 컴퓨터 단층촬영 장치 | |
| CN114255296B (zh) | 基于单幅x射线影像的ct影像重建方法及装置 | |
| JP2014064606A (ja) | 画像処理装置及び経時差分画像生成方法 | |
| CN114255296A (zh) | 基于单幅x射线影像的ct影像重建方法及装置 | |
| JP7471277B2 (ja) | データ整合性に基づいたコーンビームコンピュータ断層撮影における効率的な動き補償 |