JPH0454455B2 - - Google Patents
Info
- Publication number
- JPH0454455B2 JPH0454455B2 JP59111869A JP11186984A JPH0454455B2 JP H0454455 B2 JPH0454455 B2 JP H0454455B2 JP 59111869 A JP59111869 A JP 59111869A JP 11186984 A JP11186984 A JP 11186984A JP H0454455 B2 JPH0454455 B2 JP H0454455B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- storage device
- image storage
- adder
- mask
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T12/00—Tomographic reconstruction from projections
- G06T12/30—Image post-processing, e.g. metal artefact correction
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Closed-Circuit Television Systems (AREA)
- Image Analysis (AREA)
- Medical Treatment And Welfare Office Work (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Processing Or Creating Images (AREA)
Description
【発明の詳細な説明】
（技術分野）
本発明は、Ｘ線CT（CTはComputed
Tomographyの略）像に生じる骨のビームハー
ドニングアーチフアクトを再構成像のみを使つて
高速に除去するようにしたＸ線CT画像処理装置
に関する。
Tomographyの略）像に生じる骨のビームハー
ドニングアーチフアクトを再構成像のみを使つて
高速に除去するようにしたＸ線CT画像処理装置
に関する。
（従来技術）
骨のビームハードニング補正は、従来ジヨセフ
他の文献（Joseph−Spital−Ａ Method for
Correcting Bone Induced Artifacts in
Computed Tomography Scanners，Journal of
Computer Assisted Tomography Vol2 ｐ
100〜108，Jan.1978）に見られるように、再構成
像を画像上でスキヤンしその中に含まれるソフト
テイシユー部と骨（CT値で分離される）の経路
長から補正量を計算し、これをもとに実スキヤン
データを補正した上で再び再構成する方法で行わ
れていた。
他の文献（Joseph−Spital−Ａ Method for
Correcting Bone Induced Artifacts in
Computed Tomography Scanners，Journal of
Computer Assisted Tomography Vol2 ｐ
100〜108，Jan.1978）に見られるように、再構成
像を画像上でスキヤンしその中に含まれるソフト
テイシユー部と骨（CT値で分離される）の経路
長から補正量を計算し、これをもとに実スキヤン
データを補正した上で再び再構成する方法で行わ
れていた。
しかしながら、この方法では次のような問題点
があつた。
があつた。
２度の再構成と画像上でのスキヤンが必要な
ため、演算量が莫大であり、高速処理が期待で
きない。
ため、演算量が莫大であり、高速処理が期待で
きない。
実スキヤンデータが必要なため事後処理が容
易でない。
易でない。
（発明の目的）
本発明は、このような点に鑑みてなされたもの
で、その目的は、アーチフアクトの除去効果は損
わず、高速かつ簡単に画像データのみから骨のビ
ームハードニング補正を行い得るＸ線CT画像処
理装置を提供することにある。
で、その目的は、アーチフアクトの除去効果は損
わず、高速かつ簡単に画像データのみから骨のビ
ームハードニング補正を行い得るＸ線CT画像処
理装置を提供することにある。
（発明の構成）
このような目的を達成する第１の発明は、外部
装置に対しＸ線CTの再構成像を入出力すること
のできる第１の画像記憶装置と、この第１の画像
記憶装置から与えられる画像を予め定められた
CT値をもとに２値化する２値化装置と、この２
値画像を記憶しておくための第２の画像記憶装置
と、重み関数を蓄えた重み関数記憶装置と、前記
２値画像をマスクとして前記重み関数を加算する
マスク付加算器と、このマスク付加算器で得られ
た加算結果を蓄える第３の画像記憶装置と、前記
第１の画像記憶装置と前記第３の画像記憶装置の
画像間の加算を行い、加算結果を第１の画像記憶
装置に与える画像加算器とを具備し、再構成像の
骨成分を表す全画素でこの点を中心として広がる
２次元の重み関数の総和をとつて補正像とし、こ
の補正像をもとの再構成像に加算することにより
骨のビームハードニング補正を行うようにしたこ
とを特徴とするもので、第２の発明は、外部装置
に対してＸ線CTの再構成像を入出力することの
できる第１の画像記憶装置と、この第１の画像記
憶装置から与えられる画像を縦横ともに縮小する
画像縮小装置と、この縮小画像を記憶する第４の
画像記憶装置と、この第４の画像記憶装置から与
えられる画像を予め定められたCT値をもとに２
値化する２値化装置と、この２値画像を記憶して
おくための第２の画像記憶装置と、重み関数を蓄
えた重み関数記憶装置と、前記２値画像をマスク
として前記重み関数を加算するマスク付加算器
と、このマスク付加算器で得られた加算結果を蓄
える第３の画像記憶装置と、この第３の画像記憶
装置から与えられる画像を補間を施して縦横とも
にもとの大きさに拡大する画像拡大装置と、前記
第１の画像記憶装置と前記画像拡大装置の画像間
の加算を行い、加算結果を第１の画像記憶装置に
与える画像加算器とを具備し、骨を表す２値画像
及び重み関数を縦横１／Ｎに縮小した後で２値画
像をマスクとして重み関数を２次元で加算する処
理を骨を表す全画素について繰り返し、得られた
加算像を補間を使つて縦横Ｎ倍に拡大した上でも
との再構成像に加算して骨のビームハードニング
補正を行うようにしたことを特徴とするものであ
る。
装置に対しＸ線CTの再構成像を入出力すること
のできる第１の画像記憶装置と、この第１の画像
記憶装置から与えられる画像を予め定められた
CT値をもとに２値化する２値化装置と、この２
値画像を記憶しておくための第２の画像記憶装置
と、重み関数を蓄えた重み関数記憶装置と、前記
２値画像をマスクとして前記重み関数を加算する
マスク付加算器と、このマスク付加算器で得られ
た加算結果を蓄える第３の画像記憶装置と、前記
第１の画像記憶装置と前記第３の画像記憶装置の
画像間の加算を行い、加算結果を第１の画像記憶
装置に与える画像加算器とを具備し、再構成像の
骨成分を表す全画素でこの点を中心として広がる
２次元の重み関数の総和をとつて補正像とし、こ
の補正像をもとの再構成像に加算することにより
骨のビームハードニング補正を行うようにしたこ
とを特徴とするもので、第２の発明は、外部装置
に対してＸ線CTの再構成像を入出力することの
できる第１の画像記憶装置と、この第１の画像記
憶装置から与えられる画像を縦横ともに縮小する
画像縮小装置と、この縮小画像を記憶する第４の
画像記憶装置と、この第４の画像記憶装置から与
えられる画像を予め定められたCT値をもとに２
値化する２値化装置と、この２値画像を記憶して
おくための第２の画像記憶装置と、重み関数を蓄
えた重み関数記憶装置と、前記２値画像をマスク
として前記重み関数を加算するマスク付加算器
と、このマスク付加算器で得られた加算結果を蓄
える第３の画像記憶装置と、この第３の画像記憶
装置から与えられる画像を補間を施して縦横とも
にもとの大きさに拡大する画像拡大装置と、前記
第１の画像記憶装置と前記画像拡大装置の画像間
の加算を行い、加算結果を第１の画像記憶装置に
与える画像加算器とを具備し、骨を表す２値画像
及び重み関数を縦横１／Ｎに縮小した後で２値画
像をマスクとして重み関数を２次元で加算する処
理を骨を表す全画素について繰り返し、得られた
加算像を補間を使つて縦横Ｎ倍に拡大した上でも
との再構成像に加算して骨のビームハードニング
補正を行うようにしたことを特徴とするものであ
る。
（実施例）
以下、図面を用いて本発明の実施例を詳しく説
明する。第１図は第１の発明に係るＸ線CT画像
処理装置の一実施例を示す要部のブロツク構成図
である。この図において、１は再構成像を記憶す
る第１の画像記憶装置、２は２値化装置、３は第
２の画像記憶装置、４は重み関数記憶装置、５は
マスク付加算器、６は第３の画像記憶装置、７は
画像加算器である。
明する。第１図は第１の発明に係るＸ線CT画像
処理装置の一実施例を示す要部のブロツク構成図
である。この図において、１は再構成像を記憶す
る第１の画像記憶装置、２は２値化装置、３は第
２の画像記憶装置、４は重み関数記憶装置、５は
マスク付加算器、６は第３の画像記憶装置、７は
画像加算器である。
第１の画像記憶装置１は外部装置（図示せず）
に対し画像を入出力することができるようになつ
ている。２値化装置２は第１の画像記憶装置１の
画像を予め設定されたCT値をしきい値として２
値化するものであり、第２の画像記憶装置３はこ
の２値画像を記憶するものである。マスク付加算
器５は重み関数記憶装置４と第２の画像記憶装置
３との出力を受け、２値画像上骨を表すCT値の
高い画素位置に重み関数の中心を合わせて２次元
に広がる重み関数を第３の画像記憶装置６に加算
記憶する。画像加算器７は第３の画像記憶装置６
の画像と第１の画像記憶装置１の画像とを加算す
るもので、加算された画像は第１の画像記憶装置
１に記憶されるようになつている。
に対し画像を入出力することができるようになつ
ている。２値化装置２は第１の画像記憶装置１の
画像を予め設定されたCT値をしきい値として２
値化するものであり、第２の画像記憶装置３はこ
の２値画像を記憶するものである。マスク付加算
器５は重み関数記憶装置４と第２の画像記憶装置
３との出力を受け、２値画像上骨を表すCT値の
高い画素位置に重み関数の中心を合わせて２次元
に広がる重み関数を第３の画像記憶装置６に加算
記憶する。画像加算器７は第３の画像記憶装置６
の画像と第１の画像記憶装置１の画像とを加算す
るもので、加算された画像は第１の画像記憶装置
１に記憶されるようになつている。
このような構成において、外部装置より第１の
画像記憶装置１に入力されたＸ線CTの再構成像
は、２値化装置２において予め設定されたしきい
値と比較され２値化される。この２値画像を第２
の画像記憶装置３に順次記憶してゆく。マスク付
加算器５では、この第２の画像記憶装置３より２
値画像を受け、この２値画像上骨を表すCT値の
高い画素位置に重み関数記憶装置４からの重み関
数の中心を合わせて２次元に広がる重み関数を第
３の画像記憶装置６に加算記憶する。尚、第３の
画像記憶装置６の初期記憶値は零に設定されてい
る。この加算処理をしきい値以上のCT値を持つ
すべての画素に対して繰り返す。
画像記憶装置１に入力されたＸ線CTの再構成像
は、２値化装置２において予め設定されたしきい
値と比較され２値化される。この２値画像を第２
の画像記憶装置３に順次記憶してゆく。マスク付
加算器５では、この第２の画像記憶装置３より２
値画像を受け、この２値画像上骨を表すCT値の
高い画素位置に重み関数記憶装置４からの重み関
数の中心を合わせて２次元に広がる重み関数を第
３の画像記憶装置６に加算記憶する。尚、第３の
画像記憶装置６の初期記憶値は零に設定されてい
る。この加算処理をしきい値以上のCT値を持つ
すべての画素に対して繰り返す。
このようにしてできあがつた画像は、画像加算
器７によつて初めに入力された再構成像（第１の
画像記憶装置１に記憶されている）に加算され
る。これにより、骨のビームハードニングアーチ
フアクトを除去することができ、その画像は改め
て第１の画像記憶装置１に書き込まれ、必要に応
じて外部装置に出力される。
器７によつて初めに入力された再構成像（第１の
画像記憶装置１に記憶されている）に加算され
る。これにより、骨のビームハードニングアーチ
フアクトを除去することができ、その画像は改め
て第１の画像記憶装置１に書き込まれ、必要に応
じて外部装置に出力される。
第２図は本発明の他の実施例を示す説明図であ
る。第１図と異なるところは第３の画像記憶装置
６と画像加算器７を省略した点である。第２図に
おいては、重み関数の総和をとつて後から２枚の
画像間で加算する代りに、重み関数を再構成像の
対応する位置に直接加算し、結果を第１の画像記
憶装置１に書き込むようにしている。このような
構成によつても上述と同様に骨のビームハードニ
ングアーチフアクトを除去することができる。
る。第１図と異なるところは第３の画像記憶装置
６と画像加算器７を省略した点である。第２図に
おいては、重み関数の総和をとつて後から２枚の
画像間で加算する代りに、重み関数を再構成像の
対応する位置に直接加算し、結果を第１の画像記
憶装置１に書き込むようにしている。このような
構成によつても上述と同様に骨のビームハードニ
ングアーチフアクトを除去することができる。
尚、あるCT値以上の画素に対して同一の重み
関数を加算しているが、あるCT値以上の画素に
対してそのCT値により変調した重み関数を加算
するように構成することも可能である。
関数を加算しているが、あるCT値以上の画素に
対してそのCT値により変調した重み関数を加算
するように構成することも可能である。
又、マスク付加算器５は第２の画像記憶装置３
の２値画像と重み関数記憶装置４の重み関数の２
次元コンボリユーシヨンで置き換えるとも可能で
ある。
の２値画像と重み関数記憶装置４の重み関数の２
次元コンボリユーシヨンで置き換えるとも可能で
ある。
以上述べた第１の発明の実施例においては、骨
を表す画素の数だけ２次元の重み関数を加算して
おり、場合によつてはかなりの演算量になる。そ
こで、骨を表す２値画像重み関数を縦横１／Ｎに
縮小した後で２値画像をマスクとして重み関数を
２次元で加算し、骨を表す全画素について加算を
繰り返すようにすれば、演算量を減少させること
ができ、それだけ高速処理を実現することができ
る。
を表す画素の数だけ２次元の重み関数を加算して
おり、場合によつてはかなりの演算量になる。そ
こで、骨を表す２値画像重み関数を縦横１／Ｎに
縮小した後で２値画像をマスクとして重み関数を
２次元で加算し、骨を表す全画素について加算を
繰り返すようにすれば、演算量を減少させること
ができ、それだけ高速処理を実現することができ
る。
第３図はこのような高速処理を図つた第２の発
明におけるＸ線CT画像処理装置の一実施例を示
す構成図である。第３図において、第１図と同等
部分には同一符号を付してある。第１図と異なる
構成部分は画像縮小装置３１、第４の画像記憶装
置３２、画像拡大装置３３、第５の画像記憶装置
３４である。
明におけるＸ線CT画像処理装置の一実施例を示
す構成図である。第３図において、第１図と同等
部分には同一符号を付してある。第１図と異なる
構成部分は画像縮小装置３１、第４の画像記憶装
置３２、画像拡大装置３３、第５の画像記憶装置
３４である。
画像縮小装置３１は第１の画像記憶装置１の画
像を受け、これを縦横両方向に１／Ｎに縮小する
もので、結果は第４の画像記憶装置３２に格納さ
れる。第４の画像記憶装置３２の出力は２値化装
置２に導かれるようになつている。
像を受け、これを縦横両方向に１／Ｎに縮小する
もので、結果は第４の画像記憶装置３２に格納さ
れる。第４の画像記憶装置３２の出力は２値化装
置２に導かれるようになつている。
画像拡大装置３３は第３の画像記憶装置６の縮
小画像をＮ倍し、もとの大きさの画像に拡大する
ものである。拡大画像は第５の画像記憶装置３４
に一旦格納された後画像加算器７に送られる。
小画像をＮ倍し、もとの大きさの画像に拡大する
ものである。拡大画像は第５の画像記憶装置３４
に一旦格納された後画像加算器７に送られる。
このような構成における動作を次に説明する。
外部装置から第１の画像記憶装置１に入力された
Ｘ線CTの再構成像は、画像縮小装置３１で縦横
両方向に１／Ｎに縮小され、第４の画像記憶装置
３２に記憶される。尚、Ｎは通常整数が用いられ
る。
外部装置から第１の画像記憶装置１に入力された
Ｘ線CTの再構成像は、画像縮小装置３１で縦横
両方向に１／Ｎに縮小され、第４の画像記憶装置
３２に記憶される。尚、Ｎは通常整数が用いられ
る。
記憶装置３２の縮小画像は予め設定されたCT
値に従つて２値化装置２で２値化され、第２の画
像記憶装置３に記憶される。
値に従つて２値化装置２で２値化され、第２の画
像記憶装置３に記憶される。
重み関数記憶装置４の重み関数も前述の縮小率
１／Ｎに対応して予め縦横１／Ｎに縮小されてい
るものとする。マスク付加算器５において、重み
関数記憶装置４の重み関数を２値画像（記憶装置
３）の骨を表す各画素に対し中心を合わせて第３
の画像記憶装置６に加算記憶する。尚、第３の画
像記憶装置６の初期記憶値は零に設定されてい
る。この加算処理を骨を表す全画素について行つ
た結果は縮小像に対する補正像となつているの
で、これを画像拡大装置３３でＮ倍し、もとの大
きさに復元する。この拡大には計算時間の短縮の
ため線型補間を用い、１次元毎に２段階に分けて
２次元の補間を行う。拡大された画像は第５の画
像記憶装置３４に蓄えられ、この拡大補正像は画
像加算器７において第１の画像記憶装置１上の入
力画像と２次元加算され、その結果が第１の画像
記憶装置１上に出力される。
１／Ｎに対応して予め縦横１／Ｎに縮小されてい
るものとする。マスク付加算器５において、重み
関数記憶装置４の重み関数を２値画像（記憶装置
３）の骨を表す各画素に対し中心を合わせて第３
の画像記憶装置６に加算記憶する。尚、第３の画
像記憶装置６の初期記憶値は零に設定されてい
る。この加算処理を骨を表す全画素について行つ
た結果は縮小像に対する補正像となつているの
で、これを画像拡大装置３３でＮ倍し、もとの大
きさに復元する。この拡大には計算時間の短縮の
ため線型補間を用い、１次元毎に２段階に分けて
２次元の補間を行う。拡大された画像は第５の画
像記憶装置３４に蓄えられ、この拡大補正像は画
像加算器７において第１の画像記憶装置１上の入
力画像と２次元加算され、その結果が第１の画像
記憶装置１上に出力される。
第４図はＮ＝４の場合の拡大補間を説明する説
明図である。既知の値（○印）をもとにまず×印
での値をそれを挟む○印の値から次式で算出す
る。
明図である。既知の値（○印）をもとにまず×印
での値をそれを挟む○印の値から次式で算出す
る。
f^（ｉ）＝ｆ(1)＋（ｉ／４）・｛ｆ(2)−ｆ(1)｝
ここで、ｉ＝１，２，３
次に＊印での値をそれを挟む左右の×印の値か
ら次式で算出する。
ら次式で算出する。
ｆ（ｉ）＝f^(1)＋ｉ／４）・｛f^(2)f^(1)｝
ここで、ｉ＝１，２，３
このようにして補間を行う。Ｎ＝４以外の場合
も同様である。
も同様である。
尚、実施例に限らず、次のようにしてもよい。
画像縮小装置３１では実際には画像を縮小せ
ず、２値化装置２と一体化して縦横ともＮ画素
毎にサンプルして２値化する。
ず、２値化装置２と一体化して縦横ともＮ画素
毎にサンプルして２値化する。
画像拡大装置３３は第５の画像記憶装置３４
を通さずに直接第１の画像記憶装置と加算す
る。
を通さずに直接第１の画像記憶装置と加算す
る。
マスク付加算器５を２次元コンボリユーシヨ
ン装置に置き換える。
ン装置に置き換える。
画像縮小装置３１、第４の画像記憶装置３２
と第２の画像記憶装置３を用いずに、２値化装
置２が出力するＮサンプル毎の値に従つて加算
器５をオン・オフする。
と第２の画像記憶装置３を用いずに、２値化装
置２が出力するＮサンプル毎の値に従つて加算
器５をオン・オフする。
（発明の効果）
以上説明したように、本発明によれば、重み関
数を適切に設定することにより、ジヨセフ他の従
来の方法によるのと同等の骨とビームハードニン
グ補正ができ、脳と頭蓋骨の画像上のきれが良
く、しかも骨によるカツピングのない画像を得る
ことができる。
数を適切に設定することにより、ジヨセフ他の従
来の方法によるのと同等の骨とビームハードニン
グ補正ができ、脳と頭蓋骨の画像上のきれが良
く、しかも骨によるカツピングのない画像を得る
ことができる。
又、従来のような画像上で改めてスキヤンを繰
り返し、その経路内での成分を調べ、それによる
実スキヤンデータの補正を行い、それを再び再構
成するという複雑かつ莫大な演算量を要する処理
と比べ、再構成像をもとに簡単な画像間の演算を
繰り返すだけで処理ができるなど、アルゴリズム
の簡単さ、演算量の少なさから極めて実用的であ
る。
り返し、その経路内での成分を調べ、それによる
実スキヤンデータの補正を行い、それを再び再構
成するという複雑かつ莫大な演算量を要する処理
と比べ、再構成像をもとに簡単な画像間の演算を
繰り返すだけで処理ができるなど、アルゴリズム
の簡単さ、演算量の少なさから極めて実用的であ
る。
例えば、使用画像320×320画素、使用重み関数
61×61画素としたとき、ジヨセフ他の従来方法で
は約15時間かかつていたが、本発明によれば約50
分の処理時間ですむ。
61×61画素としたとき、ジヨセフ他の従来方法で
は約15時間かかつていたが、本発明によれば約50
分の処理時間ですむ。
又、第３図に示すような本発明によれば、最初
に画像を１／N2に縮小し、この縮小画像に対し
て２値化、重み関数の加算を行うため、重み関数
ももとの１／N2に縮小したものを使うことから
加算回数はM2×ｋから（Ｍ／Ｎ）2×ｋ／N2と
１／N4に減少する。ここでM2は重み関数の面積
を画素数で表したもの、ｋは骨を表す画素の数で
ある。
に画像を１／N2に縮小し、この縮小画像に対し
て２値化、重み関数の加算を行うため、重み関数
ももとの１／N2に縮小したものを使うことから
加算回数はM2×ｋから（Ｍ／Ｎ）2×ｋ／N2と
１／N4に減少する。ここでM2は重み関数の面積
を画素数で表したもの、ｋは骨を表す画素の数で
ある。
計算機によるシミユレーシヨンによれば、Ｎ＝
４、Ｍ＝61で画質の劣化は±1CT値以下であり、
計算時間は縮小しない場合に比べて1/203倍とな
つた。
４、Ｍ＝61で画質の劣化は±1CT値以下であり、
計算時間は縮小しない場合に比べて1/203倍とな
つた。
この方式は骨のビームハードニング補正像が低
い空間周波数成分しか含んでいないという現象を
巧みに利用した結果のたまものである。
い空間周波数成分しか含んでいないという現象を
巧みに利用した結果のたまものである。
第１図は第１の発明に係るＸ線CT画像処理装
置の一実施例を示す構成図、第２図は第１の発明
の他の実施例を示す説明図、第３図は第２の発明
の一実施例を示す構成図、第４図は拡大補間を説
明するための説明図である。 １……第１の画像記憶装置、２……２値化装
置、３……第２の画像記憶装置、４……重み関数
記憶装置、５……マスク付加算器、６……第３の
画像記憶装置、７……画像加算器、３１……画像
縮小装置、３２……第４の画像記憶装置、３３…
…画像拡大装置、３４……第５の画像記憶装置。
置の一実施例を示す構成図、第２図は第１の発明
の他の実施例を示す説明図、第３図は第２の発明
の一実施例を示す構成図、第４図は拡大補間を説
明するための説明図である。 １……第１の画像記憶装置、２……２値化装
置、３……第２の画像記憶装置、４……重み関数
記憶装置、５……マスク付加算器、６……第３の
画像記憶装置、７……画像加算器、３１……画像
縮小装置、３２……第４の画像記憶装置、３３…
…画像拡大装置、３４……第５の画像記憶装置。
Claims (1)
- 【特許請求の範囲】 １ 外部装置に対しＸ線CTの再構成像を入出力
することのできる第１の画像記憶装置と、この第
１の画像記憶装置から与えられる画像を予め定め
られたCT値をもとに２値化する２値化装置と、
この２値画像を記憶しておくための第２の画像記
憶装置と、重み関数を蓄えた重み関数記憶装置
と、前記２値画像をマスクとして前記重み関数を
加算するマスク付加算器と、このマスク付加算器
で得られた加算結果を蓄える第３の画像記憶装置
と、前記第１の画像記憶装置と前記第３の画像記
憶装置の画像間の加算を行い、加算結果を第１の
画像記憶装置に与える画像加算器とを具備し、再
構成像の骨成分を表す企画素でこの点を中心とし
て広がる２次元の重み関数の総和をとつて補正像
とし、この補正像をもとの再構成像に加算するこ
とにより骨のビームハードニング補正を行うよう
にしたことを特徴とするＸ線CT画像処理装置。 ２ 外部装置に対してＸ線CTの再構成像を入出
力することのできる第１の画像記憶装置と、この
第１の画像記憶装置から与えられる画像を縦横と
もに縮小する画像縮小装置と、この縮小画像を記
憶する第４の画像記憶装置と、この第４の画像記
憶装置から与えられる画像を予め定められたCT
値をもとに２値化する２値化装置と、この２値画
像を記憶しておくための第２の画像記憶装置と、
重み関数を蓄えた重み関数記憶装置と、前記２値
画像をマスクとして前記重み関数を加算するマス
ク付加算器と、このマスク付加算器で得られた加
算結果を蓄える第３の画像記憶装置と、この第３
の画像記憶装置から与えられる画像を補間を施し
て縦横ともにもとの大きさに拡大する画像拡大装
置と、前記第１の画像記憶装置と前記画像拡大装
置の画像間の加算を行い、加算結果を第１の画像
記憶装置に与える画像加算器とを具備し、骨を表
す２値画像及び重み関数を縦横１／Ｎに縮小した
後で２値画像をマスクとして重み関数を２次元で
加算する処理を骨を表す全画素について繰り返
し、得られた加算像を補間を使つて縦横Ｎ倍に拡
大した上でもとの再構成像に加算して骨のビーム
ハードニング補正を行うようにしたことを特徴と
するＸ線CT画像処理装置。
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59111869A JPS60254375A (ja) | 1984-05-31 | 1984-05-31 | Ｘ線ｃｔ画像処理装置 |
| DE8585901557T DE3586203T2 (de) | 1984-03-15 | 1985-03-13 | Bildverarbeitungsgeraet fuer roentgentomographie. |
| PCT/JP1985/000124 WO1985004090A1 (fr) | 1984-03-15 | 1985-03-13 | Processeur d'images ct a rayons x |
| US06/800,077 US4739481A (en) | 1984-03-15 | 1985-03-13 | X-ray CT image processing apparatus |
| DE198585901557T DE204844T1 (de) | 1984-03-15 | 1985-03-13 | Bildverarbeitungsgeraet fuer roentgentomographie. |
| EP85901557A EP0204844B1 (en) | 1984-03-15 | 1985-03-13 | X-ray ct image processor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59111869A JPS60254375A (ja) | 1984-05-31 | 1984-05-31 | Ｘ線ｃｔ画像処理装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS60254375A JPS60254375A (ja) | 1985-12-16 |
| JPH0454455B2 true JPH0454455B2 (ja) | 1992-08-31 |
Family
ID=14572198
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP59111869A Granted JPS60254375A (ja) | 1984-03-15 | 1984-05-31 | Ｘ線ｃｔ画像処理装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS60254375A (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0644280B2 (ja) * | 1989-10-31 | 1994-06-08 | 株式会社島津製作所 | 画像フイルタリング装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB1578442A (en) * | 1976-07-01 | 1980-11-05 | Emi Ltd | Radiography |
-
1984
- 1984-05-31 JP JP59111869A patent/JPS60254375A/ja active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPS60254375A (ja) | 1985-12-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3882651B2 (ja) | 画像処理装置およびプログラム | |
| JP3995854B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP4598507B2 (ja) | 最小誤差時空間再帰フィルタを使用する画像ノイズ低減のためのシステム及び方法 | |
| US6289133B1 (en) | Image processing method and apparatus | |
| EP1316919A2 (en) | Method for contrast-enhancement of digital portal images | |
| JPS6041170A (ja) | 映像再投映システム | |
| JP2002133410A (ja) | ノイズ抑制処理装置並びに記録媒体 | |
| US6947176B1 (en) | Method for correcting lightness of image | |
| JPH10171983A (ja) | 画像処理方法および装置 | |
| Van Tran et al. | Flexible mask subtraction for digital angiography | |
| JP2004242285A (ja) | ノイズ抑制処理方法および装置並びにプログラム | |
| Jeong et al. | Multi-frame example-based super-resolution using locally directional self-similarity | |
| EP0351062A1 (en) | Method and apparatus for generating composite images | |
| CN113793272A (zh) | 图像降噪方法及装置、存储介质、终端 | |
| JPH08294001A (ja) | 画像処理方法および画像処理装置 | |
| WO1985004090A1 (fr) | Processeur d'images ct a rayons x | |
| JP3738791B2 (ja) | 画像処理方法および装置 | |
| Cox et al. | Automatic registration of temporal image pairs for digital subtraction angiography | |
| JP3999432B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| JPH06325165A (ja) | 放射線画像の部分の表示方法 | |
| JP3973181B2 (ja) | 目的物の一連の放射線医学画像を処理する方法 | |
| JPH1063839A (ja) | 画像処理方法および装置 | |
| JPH08140964A (ja) | Ｘ線ｃｔ装置 | |
| JPH0454455B2 (ja) | ||
| JP2002183727A (ja) | 画像処理装置 |