JPH0230066B2 - - Google Patents
Info
- Publication number
- JPH0230066B2 JPH0230066B2 JP58022142A JP2214283A JPH0230066B2 JP H0230066 B2 JPH0230066 B2 JP H0230066B2 JP 58022142 A JP58022142 A JP 58022142A JP 2214283 A JP2214283 A JP 2214283A JP H0230066 B2 JPH0230066 B2 JP H0230066B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- data
- projection data
- digital image
- image data
- 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 - Lifetime
Links
- 238000012958 reprocessing Methods 0.000 claims description 21
- 238000012545 processing Methods 0.000 claims description 20
- 238000000034 method Methods 0.000 description 18
- 238000005510 radiation hardening Methods 0.000 description 16
- 238000003702 image correction Methods 0.000 description 14
- 210000000988 bone and bone Anatomy 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 7
- 210000004872 soft tissue Anatomy 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000005855 radiation Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000001678 irradiating effect Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
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]
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Description
【発明の詳細な説明】
〔発明の技術分野〕
この発明は、画像たとえばＸ線CT装置により
得られる画像を補正する画像処理装置に関する。
得られる画像を補正する画像処理装置に関する。
従来、画像診断装置たとえばＸ線CT装置で用
いられるＸ線は、通常、多色Ｘ線であるため、Ｘ
線エネルギにより被検体中にある骨や軟組織等の
各部位のＸ線吸収効率が相違するので、線質硬化
が生じる。線質硬化は、Ｘ線CT装置においてプ
ロジエクシヨンデータにより再構成された画像に
アーチフアクトを生じさせる。このアーチフアク
トを除去するために、従来では、セカンドパス線
質硬化補正処理（2nd pass beam hardening
correction）が行なわれている。
いられるＸ線は、通常、多色Ｘ線であるため、Ｘ
線エネルギにより被検体中にある骨や軟組織等の
各部位のＸ線吸収効率が相違するので、線質硬化
が生じる。線質硬化は、Ｘ線CT装置においてプ
ロジエクシヨンデータにより再構成された画像に
アーチフアクトを生じさせる。このアーチフアク
トを除去するために、従来では、セカンドパス線
質硬化補正処理（2nd pass beam hardening
correction）が行なわれている。
Ｘ線CT装置におけるセカンドパス線質硬化捕
正処理の一例を次に述べる。
正処理の一例を次に述べる。
被検体の体軸を中心として被検体の周囲にＸ
線を回転させながら、Ｘ線を各回転角度毎に曝
射して得られる全てのプロジエクシヨンデータ
を基に被検体の断層像を先ず再構成する。
線を回転させながら、Ｘ線を各回転角度毎に曝
射して得られる全てのプロジエクシヨンデータ
を基に被検体の断層像を先ず再構成する。
次いで、前記断層像からＢ（θ、ｔ）、および
Ｓ（θ、ｔ）を求める。ここで、Ｂ（θ、ｔ）は
第１図に示すように、たとえば骨２と軟組織４
との断層像において、xy座標におけるｙ軸と
角度θをもつて交差すると共に原点Ｏより距離
ｔを有するビーム（θ、ｔ）上における骨の総
量（ｇ／mm）を示すパラメータであり、Ｓ（θ、
ｔ）は、前記ビーム方向における軟組織の総量
（ｇ／mm）を示すパラメータである。
Ｓ（θ、ｔ）を求める。ここで、Ｂ（θ、ｔ）は
第１図に示すように、たとえば骨２と軟組織４
との断層像において、xy座標におけるｙ軸と
角度θをもつて交差すると共に原点Ｏより距離
ｔを有するビーム（θ、ｔ）上における骨の総
量（ｇ／mm）を示すパラメータであり、Ｓ（θ、
ｔ）は、前記ビーム方向における軟組織の総量
（ｇ／mm）を示すパラメータである。
次いでビーム（θ、ｔ）において、線質硬化
の大きさＣ（θ、ｔ）を、表あるいは関数によ
り求める。
の大きさＣ（θ、ｔ）を、表あるいは関数によ
り求める。
次いで線質硬化の大きさＣ（θ、ｔ）をパラ
レル投影データとして、これを再構成すること
により補正用画像を得、前記で求めた画像か
らこの補正用画像を減算することにより、線質
硬化による影響を除去した画像を得る。
レル投影データとして、これを再構成すること
により補正用画像を得、前記で求めた画像か
らこの補正用画像を減算することにより、線質
硬化による影響を除去した画像を得る。
前記の処理は、第(1)式により行なうことがで
きる。
きる。
Ｃ（θ、ｔ）≒Ｔ（Ｂ（θ、ｔ）、Ｓ（θ、ｔ））
………(1) ここで、関数Ｔは、第(2)式に示すように近似す
ることができる。
………(1) ここで、関数Ｔは、第(2)式に示すように近似す
ることができる。
Ｔ（Ｂ、Ｓ）
≒R1（Ｖ）＋R2（V2）＋…＋Rn（Vn） ………(2)
ただし、Rk（ｋは１からｎまでの整数）は関数
であり、またＶは第(3)式により表わすことができ
る。
であり、またＶは第(3)式により表わすことができ
る。
Ｖ＝kBＢ＋kSＳ ………(3)
ただし、kB、kSは定数である。
したがつて、第(1)式のＣ（θ、ｔ）は、第(4)式
に置きかえることができる。
に置きかえることができる。
Ｃ（θ、ｔ）
≒R1（Ｖ（θ、ｔ））＋R2（（Ｖ（θ、ｔ））2）
＋…＋Rn（（Ｖ（θ、ｔ））n） ………(4)
第(4)式におけるＶ（θ、ｔ）は、前記の処理
により求めることもできるが、次の処理を行なつ
てもよい。つまり、断層像上の各ピクセルＰ（ｘ、
ｙ）のCT値μ（ｐ）について もしμ（ｐ）＞骨についての基準値、であれば、 Ｘ＝μ（ｐ）×kBで表わされるＸは骨を示すピク
セルであり、 もしμ（ｐ）＞軟組織についての基準値、であれ
ば、 Ｘ＝μ（ｐ）×kSで表わされるＸは軟組織を示す
ピクセルであり、 Ｘが前記以外のときはＸ＝０とし、そのピクセ
ルは空気を表示する、との一連のアルゴリズムに
よりγ（ｘ、ｙ）を求める。このγ（ｘ、ｙ）は、
１枚の画像であり、しかも、γ（ｘ、ｙ）は画像
μ（ｘ、ｙ）に対して画像処理を行なつて得たも
のである。したがつて、Ｖ（θ、ｔ）は、γ（ｘ、
ｙ）により第(5)式のように表示することができ
る。
により求めることもできるが、次の処理を行なつ
てもよい。つまり、断層像上の各ピクセルＰ（ｘ、
ｙ）のCT値μ（ｐ）について もしμ（ｐ）＞骨についての基準値、であれば、 Ｘ＝μ（ｐ）×kBで表わされるＸは骨を示すピク
セルであり、 もしμ（ｐ）＞軟組織についての基準値、であれ
ば、 Ｘ＝μ（ｐ）×kSで表わされるＸは軟組織を示す
ピクセルであり、 Ｘが前記以外のときはＸ＝０とし、そのピクセ
ルは空気を表示する、との一連のアルゴリズムに
よりγ（ｘ、ｙ）を求める。このγ（ｘ、ｙ）は、
１枚の画像であり、しかも、γ（ｘ、ｙ）は画像
μ（ｘ、ｙ）に対して画像処理を行なつて得たも
のである。したがつて、Ｖ（θ、ｔ）は、γ（ｘ、
ｙ）により第(5)式のように表示することができ
る。
Ｖ（θ、ｔ）＝∫∞ -∞γ（tcosθ
＋hsioθ、hsioθ−tsioθ）×dh ………(5)
ここで、ｈはビーム（θ、ｔ）におけるビーム
上での位置を示す変数である。この第(5)式は、Ｖ
（θ、ｔ）が画像γ（ｘ、ｙ）のパラレル投影デー
タであることを示している。結局、Ｖ（θ、ｔ）
を求めるためには、まずγ（ｘ、ｙ）を求め、次
いでγ（ｘ、ｙ）の投影データを求めればよい。
上での位置を示す変数である。この第(5)式は、Ｖ
（θ、ｔ）が画像γ（ｘ、ｙ）のパラレル投影デー
タであることを示している。結局、Ｖ（θ、ｔ）
を求めるためには、まずγ（ｘ、ｙ）を求め、次
いでγ（ｘ、ｙ）の投影データを求めればよい。
しかしながら、γ（ｘ、ｙ）の投影データＶ
（θ、ｔ）を第(5)式にしたがつて求めようとする
と、計算量が膨大となり、計算に時間がかかると
の問題点がある。
（θ、ｔ）を第(5)式にしたがつて求めようとする
と、計算量が膨大となり、計算に時間がかかると
の問題点がある。
この発明は前記事情に鑑みてなされたものであ
り、画像γから投影データＶを迅速に算出するこ
とのできる画像再処理装置を提供することを目的
とするものである。
り、画像γから投影データＶを迅速に算出するこ
とのできる画像再処理装置を提供することを目的
とするものである。
前記目的を達成するためのこの発明の概要は、
画像処理する画像再処理装置において、既に処理
されて得られるところの画像を構成する各画像デ
ータを２次元フーリエ変換し、次いで極座標変換
した後、１次元逆フーリエ変換するプロセツサを
有することを特徴とするものである。
画像処理する画像再処理装置において、既に処理
されて得られるところの画像を構成する各画像デ
ータを２次元フーリエ変換し、次いで極座標変換
した後、１次元逆フーリエ変換するプロセツサを
有することを特徴とするものである。
先ずこの発明についての原理的説明を行なう。
この発明における画像処理装置においては、得
られる原画像から押出したところの、誤差を生じ
させる部位の抽出画像から以下に手順により抽出
画像についての投影データを求めることをその特
長とする。
られる原画像から押出したところの、誤差を生じ
させる部位の抽出画像から以下に手順により抽出
画像についての投影データを求めることをその特
長とする。
一般に、投影データＶから画像γを求める処理
は画像再構成処理と称され、投影データＶが平行
ビームにより求められるのであれば、第(6)式に示
す２次元フーリエ変換法により投影データＶから
画像γを再構成することができる。
は画像再構成処理と称され、投影データＶが平行
ビームにより求められるのであれば、第(6)式に示
す２次元フーリエ変換法により投影データＶから
画像γを再構成することができる。
γ＝F-2（P-1（Ｆ（Ｖ））） ………(6)
ここで、Ｆは、Ｖ〔以下、Ｖ（θ、ｔ）と示すこ
とがある。〕をｔについてフーリエ変換してυ
（θ、ω）を求める処理であり、υ（θ、ω）は極
座標で示す画像である。P-1は、υ（θ、ω）を
直交座標で示す画像Γ（ξ、η）を求める処理で
あり、F-2は、画像Γ（ξ、η）をξとηとにつ
いて２次元フーリエ変換して画像γ（ｘ、ｙ）を
求める処理である。
とがある。〕をｔについてフーリエ変換してυ
（θ、ω）を求める処理であり、υ（θ、ω）は極
座標で示す画像である。P-1は、υ（θ、ω）を
直交座標で示す画像Γ（ξ、η）を求める処理で
あり、F-2は、画像Γ（ξ、η）をξとηとにつ
いて２次元フーリエ変換して画像γ（ｘ、ｙ）を
求める処理である。
したがつて、抽出画像をγとし、その投影デー
タをＶとすると、抽出画像γから第(7)式にしたが
つてその投影データを求めることができる。
タをＶとすると、抽出画像γから第(7)式にしたが
つてその投影データを求めることができる。
Ｖ＝F-1（Ｐ（F2（γ））） ………(7)
第(7)式に示す処理手順は、抽出画像γ（ｘ、ｙ）
を２次元フーリエ変換処理F2した後、極座標変
換処理Ｐをし、次いで、１次元逆フーリエ変換処
理F-1をするものである。２次元フーリエ変換処
理F2および１次元逆フーリエ変換処理F-1は、高
速フーリエ変換処理（FFT）により迅速に行な
うことができる。抽出画像から投影データが求ま
れば、投影データ上の誤差量も簡単に求まり、原
画像の誤差を除去することができる。
を２次元フーリエ変換処理F2した後、極座標変
換処理Ｐをし、次いで、１次元逆フーリエ変換処
理F-1をするものである。２次元フーリエ変換処
理F2および１次元逆フーリエ変換処理F-1は、高
速フーリエ変換処理（FFT）により迅速に行な
うことができる。抽出画像から投影データが求ま
れば、投影データ上の誤差量も簡単に求まり、原
画像の誤差を除去することができる。
次に、この発明の前記原理を適用するこの発明
の一実施例である画像再処理装置を利用してＸ線
CT装置により得られる断層像における線質硬化
を補正する画像補正装置について説明する。
の一実施例である画像再処理装置を利用してＸ線
CT装置により得られる断層像における線質硬化
を補正する画像補正装置について説明する。
第２図はこの発明の一実施例である画像再処理
装置を有する画像補正装置である。
装置を有する画像補正装置である。
第２図に示すように、画像補正装置１０は、Ｘ
線CT装置において投影データを再構成して得た
原画像についてのデジタル画像データμ（ｘ、ｙ）
を補正して線質硬化による誤差を除去した画像に
ついてのデジタル画像データμ′（ｘ、ｙ）を出力
するものであり、第１の画像メモリ１２、画像処
理ユニツト１４、第２の画像メモリ１６、画像再
処理装置１８、投影データメモリ２０、線質硬化
量（BH）推定ユニツト２２、テーブル２４、再
構成ユニツト２６、および演算ユニツト２８を有
する。
線CT装置において投影データを再構成して得た
原画像についてのデジタル画像データμ（ｘ、ｙ）
を補正して線質硬化による誤差を除去した画像に
ついてのデジタル画像データμ′（ｘ、ｙ）を出力
するものであり、第１の画像メモリ１２、画像処
理ユニツト１４、第２の画像メモリ１６、画像再
処理装置１８、投影データメモリ２０、線質硬化
量（BH）推定ユニツト２２、テーブル２４、再
構成ユニツト２６、および演算ユニツト２８を有
する。
第１の画像メモリ１２は、デジタル画像データ
μ（ｘ、ｙ）を格納した後、これを読み出して画
像処理ユニツト１４および演算ユニツト２８に出
力し、また、演算ユニツト２８により得られた線
質硬化による誤差を除去したデジタル画像データ
μ′（ｘ、ｙ）を入力し、一旦これμ′（ｘ、ｙ）を格
納した後、デジタル画像データμ′（ｘ、ｙ）をた
とえば画像表示装置に出力する記憶装置である。
μ（ｘ、ｙ）を格納した後、これを読み出して画
像処理ユニツト１４および演算ユニツト２８に出
力し、また、演算ユニツト２８により得られた線
質硬化による誤差を除去したデジタル画像データ
μ′（ｘ、ｙ）を入力し、一旦これμ′（ｘ、ｙ）を格
納した後、デジタル画像データμ′（ｘ、ｙ）をた
とえば画像表示装置に出力する記憶装置である。
画像処理ユニツト１４は、第１の画像メモリ１
２に格納したデジタル画像データμ（ｘ、ｙ）を
読み出し、デジタル画像データμ（ｘ、ｙ）によ
り示される画像中の一定値以上のCT値を与える
部位の画像を抽出し、抽出した画像についてのデ
ジタル画像データγ（ｘ、ｙ）を出力する画像プ
ロセツサである。
２に格納したデジタル画像データμ（ｘ、ｙ）を
読み出し、デジタル画像データμ（ｘ、ｙ）によ
り示される画像中の一定値以上のCT値を与える
部位の画像を抽出し、抽出した画像についてのデ
ジタル画像データγ（ｘ、ｙ）を出力する画像プ
ロセツサである。
第２の画像メモリ１６は、画像処理ユニツト１
４より出力されるデジタル画像データγ（ｘ、ｙ）
を一旦格納した後、これを画像再処理装置１８に
出力し、また、再構成ユニツト２６より出力され
るデジタル画像データＥ（ｘ、ｙ）を一旦格納し
た後、これを演算ユニツト２８に出力する記憶装
置である。
４より出力されるデジタル画像データγ（ｘ、ｙ）
を一旦格納した後、これを画像再処理装置１８に
出力し、また、再構成ユニツト２６より出力され
るデジタル画像データＥ（ｘ、ｙ）を一旦格納し
た後、これを演算ユニツト２８に出力する記憶装
置である。
画像再処理装置１８は、第２の画像メモリ１６
より読み出したデジタル画像データγ（ｘ、ｙ）
を高速フーリエ変換器（FFT）により２次元フ
ーリエ変換処理F2し、得られるデータΓ（ξ、
η）を極座標変換処理Ｐし、次いで得られるデー
タυ（θ、ω）を高速フーリエ変換器（FFT）に
より１次元逆フーリエ変換処理F-1して投影デー
タＶ（θ、ｔ）にするデータプロセツサである。
より読み出したデジタル画像データγ（ｘ、ｙ）
を高速フーリエ変換器（FFT）により２次元フ
ーリエ変換処理F2し、得られるデータΓ（ξ、
η）を極座標変換処理Ｐし、次いで得られるデー
タυ（θ、ω）を高速フーリエ変換器（FFT）に
より１次元逆フーリエ変換処理F-1して投影デー
タＶ（θ、ｔ）にするデータプロセツサである。
投影データメモリ２０は、投影データＶ（θ、
ｔ）を一旦格納した後、これをBH推定ユニツト
２２に出力し、また、BH推定ユニツト２２より
出力されるパラレル投影データＣ（θ、ｔ）を一
旦格納した後、これを再構成ユニツト２６に出力
する記憶装置である。
ｔ）を一旦格納した後、これをBH推定ユニツト
２２に出力し、また、BH推定ユニツト２２より
出力されるパラレル投影データＣ（θ、ｔ）を一
旦格納した後、これを再構成ユニツト２６に出力
する記憶装置である。
BH推定ユニツト２２は、入力する投影データ
Ｖ（θ、ｔ）を基に、テーブル２４に格納されて
いる表あるいは関数R0、R1、…、Rnを参照して
線質硬化によつて生ずる誤差の大きさBHを求
め、これをパラレル投影データＣ（θ、ｔ）とし
て出力するデータプロセツサである。
Ｖ（θ、ｔ）を基に、テーブル２４に格納されて
いる表あるいは関数R0、R1、…、Rnを参照して
線質硬化によつて生ずる誤差の大きさBHを求
め、これをパラレル投影データＣ（θ、ｔ）とし
て出力するデータプロセツサである。
再構成ユニツト２６は、全てのθについてのパ
ラレル投影データＣ（θ、ｔ）を基に画像再構成
し、補正用画像についてのデジタル画像データＥ
（ｘ、ｙ）を第２の画像メモリ１６に出力するデ
ータプロセツサである。
ラレル投影データＣ（θ、ｔ）を基に画像再構成
し、補正用画像についてのデジタル画像データＥ
（ｘ、ｙ）を第２の画像メモリ１６に出力するデ
ータプロセツサである。
演算ユニツト２８は、第１の画像メモリ１２に
格納されているデジタル画像データμ（ｘ、ｙ）
から、第２の画像メモリ１６より読み出した補正
用画像についてのデジタル画像データＥ（ｘ、ｙ）
を減算することにより線質硬化の補正がされた画
像についてのデジタル画像データμ′（ｘ、ｙ）を
出力するデータプロセツサである。
格納されているデジタル画像データμ（ｘ、ｙ）
から、第２の画像メモリ１６より読み出した補正
用画像についてのデジタル画像データＥ（ｘ、ｙ）
を減算することにより線質硬化の補正がされた画
像についてのデジタル画像データμ′（ｘ、ｙ）を
出力するデータプロセツサである。
次に、以上構成の画像処理装置１０の作用につ
いて述べる。
いて述べる。
図示しないＸ線CT装置において得られるプロ
ジエクシヨンデータを再構成して得た原画像たと
えば骨と軟組織を表示する断層像についてのデジ
タル画像データμ（ｘ、ｙ）を一旦第１の画像メ
モリ１２に格納した後、前記デジタル画像データ
μ（ｘ、ｙ）を画像処理ユニツト１４および演算
ユニツト２８に出力する。画像処理ユニツト１４
では、入力するデジタル画像データμ（ｘ、ｙ）
により示される画像中の一定値以上のCT値を与
える部位たとえば骨の断層像を抽出し、抽出した
断層像についてのデジタル画像データγ（ｘ、ｙ）
を出力する。デジタル画像データγ（ｘ、ｙ）は、
第２の画像メモリ１６を介して、画像再処理装置
１８に出力され、画像再処理装置１８において、
高速フーリエ変換アルゴリズム（FFT）により
２次元フーリエ変換処理F2されてデジタル画像
データΓ（ξ、η）となり、次いでデジタル画像
データΓ（ξ、η）は極座標変換処理Ｐされてデ
ジタル画像データυ（θ、ω）となり、次いでデ
ジタル画像データυ（θ、ω）は高速フーリエ変
換アルゴリズム（FFT）により１次元逆フーリ
エ変換されて投影データＶ（θ、ｔ）となる。投
影データＶ（θ、ｔ）は、投影データメモリ２０
を介してBH推定ユニツト２２に出力される。
BH推定ユニツト２２は、テーブル２４に格納さ
れているたとえば関数R0、R1、…、Rnを参照
し、前記第４式の演算を行なつて、投影データＶ
（θ、ｔ）からパラレル投影データＣ（θ、ｔ）を
求め、得られるパラレル投影データＣ（θ、ｔ）
を、投影データメモリ２０を介して再構成ユニツ
ト２６に出力する。再構成ユニツト２６は、入力
する全てのθについてのパラレル投影データＣ
（θ、ｔ）を再構成して、補正用画像についての
デジタル画像データＥ（ｘ、ｙ）に変換し、デジ
タル画像データＥ（ｘ、ｙ）を、第２の画像メモ
リ１６を介して演算ユニツト２８に出力する。演
算ユニツト２８は、デジタル画像データμ（ｘ、
ｙ）からデジタル画像データＥ（ｘ、ｙ）を減算
することにより線質硬化の補正をし、線質硬化の
ない断層像についてのデジタル画像データμ′（ｘ、
ｙ）を出力する。
ジエクシヨンデータを再構成して得た原画像たと
えば骨と軟組織を表示する断層像についてのデジ
タル画像データμ（ｘ、ｙ）を一旦第１の画像メ
モリ１２に格納した後、前記デジタル画像データ
μ（ｘ、ｙ）を画像処理ユニツト１４および演算
ユニツト２８に出力する。画像処理ユニツト１４
では、入力するデジタル画像データμ（ｘ、ｙ）
により示される画像中の一定値以上のCT値を与
える部位たとえば骨の断層像を抽出し、抽出した
断層像についてのデジタル画像データγ（ｘ、ｙ）
を出力する。デジタル画像データγ（ｘ、ｙ）は、
第２の画像メモリ１６を介して、画像再処理装置
１８に出力され、画像再処理装置１８において、
高速フーリエ変換アルゴリズム（FFT）により
２次元フーリエ変換処理F2されてデジタル画像
データΓ（ξ、η）となり、次いでデジタル画像
データΓ（ξ、η）は極座標変換処理Ｐされてデ
ジタル画像データυ（θ、ω）となり、次いでデ
ジタル画像データυ（θ、ω）は高速フーリエ変
換アルゴリズム（FFT）により１次元逆フーリ
エ変換されて投影データＶ（θ、ｔ）となる。投
影データＶ（θ、ｔ）は、投影データメモリ２０
を介してBH推定ユニツト２２に出力される。
BH推定ユニツト２２は、テーブル２４に格納さ
れているたとえば関数R0、R1、…、Rnを参照
し、前記第４式の演算を行なつて、投影データＶ
（θ、ｔ）からパラレル投影データＣ（θ、ｔ）を
求め、得られるパラレル投影データＣ（θ、ｔ）
を、投影データメモリ２０を介して再構成ユニツ
ト２６に出力する。再構成ユニツト２６は、入力
する全てのθについてのパラレル投影データＣ
（θ、ｔ）を再構成して、補正用画像についての
デジタル画像データＥ（ｘ、ｙ）に変換し、デジ
タル画像データＥ（ｘ、ｙ）を、第２の画像メモ
リ１６を介して演算ユニツト２８に出力する。演
算ユニツト２８は、デジタル画像データμ（ｘ、
ｙ）からデジタル画像データＥ（ｘ、ｙ）を減算
することにより線質硬化の補正をし、線質硬化の
ない断層像についてのデジタル画像データμ′（ｘ、
ｙ）を出力する。
なお、前記デジタル画像データμ′（ｘ、ｙ）は、
Ｘ線CT装置における画像表示部に出力され、ビ
デオ信号に変換された後、線質硬化の除去された
アーチフアクトのない画像が表示されることにな
る。
Ｘ線CT装置における画像表示部に出力され、ビ
デオ信号に変換された後、線質硬化の除去された
アーチフアクトのない画像が表示されることにな
る。
以上詳述した画像補正装置１０は、デジタル画
像データγ（ｘ、ｙ）をデジタル画像データＶ
（θ、ｔ）に高速変換する画像再処理装置１８を
有するので、高速で短時間のうちに原画像につい
ての線質硬化を補正することができる。
像データγ（ｘ、ｙ）をデジタル画像データＶ
（θ、ｔ）に高速変換する画像再処理装置１８を
有するので、高速で短時間のうちに原画像につい
ての線質硬化を補正することができる。
以上、この発明の一実施例について詳述した
が、この発明は前記実施例に限定されるものでは
なく、この発明の要旨を変更しない範囲内で適宜
に変形して実施することができるのはいうまでも
ない。
が、この発明は前記実施例に限定されるものでは
なく、この発明の要旨を変更しない範囲内で適宜
に変形して実施することができるのはいうまでも
ない。
この発明に係る画像再処理装置は、前記のよう
な線質硬化補正をする画像補正装置に限らず、散
乱線による誤差を補正する画像補正装置にも利用
することができる。
な線質硬化補正をする画像補正装置に限らず、散
乱線による誤差を補正する画像補正装置にも利用
することができる。
また、画像再処理装置を組み込んだ画像補正装
置は、第２図に示すほかに、原画像を示すデジタ
ル画像データμ（ｘ、ｙ）から複数の画像それぞ
れについてのデジタル画像データγ1（ｘ、ｙ）、γ2
（ｘ、ｙ）、…、γo（ｘ、ｙ）を抽出し、次いで前
記デジタル画像データγ1（ｘ、ｙ）、γ2（ｘ、ｙ）、
…、γo（ｘ、ｙ）を画像再処理装置で複数の投影
データV1（θ、ｔ）、V2（θ、ｔ）、…、Vo（θ、
ｔ）に再投影し、これら投影データV1（θ、ｔ）、
V2（θ、ｔ）、…Vo（θ、ｔ）を非線形合成して一
つのパラレル投影データＣ（θ、ｔ）を得、この
投影データＣ（θ、ｔ）より画像再構成して得た
補正用画像についてのデジタル画像データＥ（ｘ、
ｙ）により原画像についてのデジタル画像データ
μ（ｘ、ｙ）を補正するように、構成してもよい。
置は、第２図に示すほかに、原画像を示すデジタ
ル画像データμ（ｘ、ｙ）から複数の画像それぞ
れについてのデジタル画像データγ1（ｘ、ｙ）、γ2
（ｘ、ｙ）、…、γo（ｘ、ｙ）を抽出し、次いで前
記デジタル画像データγ1（ｘ、ｙ）、γ2（ｘ、ｙ）、
…、γo（ｘ、ｙ）を画像再処理装置で複数の投影
データV1（θ、ｔ）、V2（θ、ｔ）、…、Vo（θ、
ｔ）に再投影し、これら投影データV1（θ、ｔ）、
V2（θ、ｔ）、…Vo（θ、ｔ）を非線形合成して一
つのパラレル投影データＣ（θ、ｔ）を得、この
投影データＣ（θ、ｔ）より画像再構成して得た
補正用画像についてのデジタル画像データＥ（ｘ、
ｙ）により原画像についてのデジタル画像データ
μ（ｘ、ｙ）を補正するように、構成してもよい。
さらに、画像再処理装置を組み込んだ画像補正
装置は、第２図に示すほかに、原画像を示すデジ
タル画像データμ（ｘ、ｙ）から複数の画像それ
ぞれについてのデジタル画像データγ1（ｘ、ｙ）、
γ2（ｘ、ｙ）、…、γo（ｘ、ｙ）を抽出し、次いで
前記デジタル画像データγ1（ｘ、ｙ）、γ2（ｘ、
ｙ）、…、γo（ｘ、ｙ）を画像再処理装置で複数の
投影データV1（θ、ｔ）、V2（θ、ｔ）、…、Vo
（θ、ｔ）に再投影し、これら投影データV1（θ、
ｔ）、V2（θ、ｔ）、…、Vo（θ、ｔ）それぞれか
らパラレル投影データC1（θ、ｔ）、C2（θ、ｔ）、
…、Co（θ、ｔ）を得、これら投影データC1（θ、
ｔ）、C2（θ、ｔ）、…、Co（θ、ｔ）それぞれに
つき画像再構成して得た画像についてのデジタル
画像データE1（ｘ、ｙ）、E2（ｘ、ｙ）、…、Eo
（ｘ、ｙ）を非線形合成して補正用画像について
のデジタル画像データＥ（ｘ、ｙ）を得、このデ
ジタル画像データＥ（ｘ、ｙ）を直接にＸ線CT装
置の画像表示部に出力し、あるいは、前記デジタ
ル画像データＥ（ｘ、ｙ）により原像についての
デジタル画像データμ（ｘ、ｙ）を補正するよう
に構成してもよい。
装置は、第２図に示すほかに、原画像を示すデジ
タル画像データμ（ｘ、ｙ）から複数の画像それ
ぞれについてのデジタル画像データγ1（ｘ、ｙ）、
γ2（ｘ、ｙ）、…、γo（ｘ、ｙ）を抽出し、次いで
前記デジタル画像データγ1（ｘ、ｙ）、γ2（ｘ、
ｙ）、…、γo（ｘ、ｙ）を画像再処理装置で複数の
投影データV1（θ、ｔ）、V2（θ、ｔ）、…、Vo
（θ、ｔ）に再投影し、これら投影データV1（θ、
ｔ）、V2（θ、ｔ）、…、Vo（θ、ｔ）それぞれか
らパラレル投影データC1（θ、ｔ）、C2（θ、ｔ）、
…、Co（θ、ｔ）を得、これら投影データC1（θ、
ｔ）、C2（θ、ｔ）、…、Co（θ、ｔ）それぞれに
つき画像再構成して得た画像についてのデジタル
画像データE1（ｘ、ｙ）、E2（ｘ、ｙ）、…、Eo
（ｘ、ｙ）を非線形合成して補正用画像について
のデジタル画像データＥ（ｘ、ｙ）を得、このデ
ジタル画像データＥ（ｘ、ｙ）を直接にＸ線CT装
置の画像表示部に出力し、あるいは、前記デジタ
ル画像データＥ（ｘ、ｙ）により原像についての
デジタル画像データμ（ｘ、ｙ）を補正するよう
に構成してもよい。
また、第３図に示すように、第２図における第
１の画像メモリ１２、画像処理ユニツト１４、第
２の画像メモリ１６、投影データメモリ２０、再
構成ユニツト２６および演算ユニツト２８を、Ｘ
線CT装置における汎用および再構成プロセツサ
３０内に含ませ、画像補正装置を、画像再処理装
置１８、BH推定ユニツト２２、およびテーブル
２４により構成してもよい。なお、第３図におい
て、３１で示すのはＸ線を曝射することにより得
られる被検体についての投影データを出力するガ
ントリであり、３２で示すのは画像データを記憶
する画像メモリであり、３３で示すのは画像表示
部である。
１の画像メモリ１２、画像処理ユニツト１４、第
２の画像メモリ１６、投影データメモリ２０、再
構成ユニツト２６および演算ユニツト２８を、Ｘ
線CT装置における汎用および再構成プロセツサ
３０内に含ませ、画像補正装置を、画像再処理装
置１８、BH推定ユニツト２２、およびテーブル
２４により構成してもよい。なお、第３図におい
て、３１で示すのはＸ線を曝射することにより得
られる被検体についての投影データを出力するガ
ントリであり、３２で示すのは画像データを記憶
する画像メモリであり、３３で示すのは画像表示
部である。
この発明によると、高速フーリエ変換アルゴリ
ズムFFTを用いて高速で画像から投影データを
求めることのできる画像再処理装置を提供するこ
とができ、この発明に係る画像再処理装置を、各
種の画像を補正するための画像補正装置に組み込
むと、高速で画像の補正をすることができる。
ズムFFTを用いて高速で画像から投影データを
求めることのできる画像再処理装置を提供するこ
とができ、この発明に係る画像再処理装置を、各
種の画像を補正するための画像補正装置に組み込
むと、高速で画像の補正をすることができる。
第１図は断層像中の特定ビーム方向（θ、ｔ）
における骨の総量についてのパラメータＢ（θ、
ｔ）と軟組織の総量についてのパラメータＳ（θ、
ｔ）とを示す説明図、第２図はこの発明の一実施
例である画像再処理装置を有する画像補正装置を
示すブロツク図、および第３図はこの発明の一実
施例である画像再処理装置を有する他の画像補正
装置を結合するＸ線CT装置を示すブロツク図で
ある。 １８……画像再処理装置。
における骨の総量についてのパラメータＢ（θ、
ｔ）と軟組織の総量についてのパラメータＳ（θ、
ｔ）とを示す説明図、第２図はこの発明の一実施
例である画像再処理装置を有する画像補正装置を
示すブロツク図、および第３図はこの発明の一実
施例である画像再処理装置を有する他の画像補正
装置を結合するＸ線CT装置を示すブロツク図で
ある。 １８……画像再処理装置。
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 JPS59149555A (ja) | 1984-08-27 |
| JPH0230066B2 true 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) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4616318A (en) * | 1983-06-07 | 1986-10-07 | Elscint, Inc. | System for reprojecting images using transform techniques |
| JP2530660B2 (ja) * | 1987-07-24 | 1996-09-04 | 松下電工株式会社 | 画像処理装置 |
| JPS6491279A (en) * | 1987-10-01 | 1989-04-10 | Toshiba Corp | System for reconstructing image |
| JPS6491278A (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 | 物理量の２次元分布からストリークを識別する方法及び装置 |
| MX2009001136A (es) * | 2006-08-03 | 2009-02-10 | Univ California | Metodos iterativos para reduccion de dosis y mejoramiento de imagen en tomografia. |
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 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPS59149555A (ja) | 1984-08-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4841874B2 (ja) | 断層合成造影における直接再生方法及び装置 | |
| US7623691B2 (en) | Method for helical windmill artifact reduction with noise restoration for helical multislice CT | |
| JP5711241B2 (ja) | コーンビームｃｔ画像におけるアーチファクト低減のための方法 | |
| JP2002191591A (ja) | Ｃｔ装置により取得した出力像に対する放射硬化の補正方法 | |
| CN105469366B (zh) | 一种消减ct图像金属伪影的解析方法 | |
| JP2017143943A (ja) | 放射線画像処理装置、方法およびプログラム | |
| US6751284B1 (en) | Method and system for tomosynthesis image enhancement using transverse filtering | |
| JP2007097977A (ja) | Ｘ線ｃｔ装置 | |
| JP2007089674A (ja) | 外観形状計測装置およびｘ線ｃｔ装置 | |
| JP2006175213A (ja) | ３次元画像処理装置 | |
| US6845143B2 (en) | CT image reconstruction | |
| JP2536024B2 (ja) | Ｘ線画像処理装置 | |
| JP6185023B2 (ja) | 断層画像生成装置、方法およびプログラム | |
| JPH0230066B2 (ja) | ||
| JP4584550B2 (ja) | Ｘ線計測装置 | |
| US20060251313A1 (en) | Method of producing a cross-sectional image | |
| JP2005040506A (ja) | Ｘ線画像処理方法及び装置 | |
| JP3223195B2 (ja) | Ｘ線ｃｔ装置 | |
| US20220319072A1 (en) | Medical image processing apparatus and medical image processing method | |
| JPH07275236A (ja) | Ｘ線ｃｔ装置 | |
| KR102248191B1 (ko) | 콘빔 컴퓨터 단층촬영 영상에서 콘빔아티팩트를 보정하는 방법 및 장치, 그것을 포함하는 콘빔 컴퓨터 단층촬영 장치 | |
| JPH11113892A (ja) | 画像再構成装置 | |
| KR102492949B1 (ko) | 의료영상의 처리장치 및 처리방법 | |
| Narkbuakaew et al. | Metal artifact reduction based on 2D-projection correction for dental cone-beam CT images | |
| JPH0570462B2 (ja) |