JPS60254375A - Ｘ線ｃｔ画像処理装置 - Google Patents
Ｘ線ｃｔ画像処理装置Info
- Publication number
- JPS60254375A JPS60254375A JP59111869A JP11186984A JPS60254375A JP S60254375 A JPS60254375 A JP S60254375A JP 59111869 A JP59111869 A JP 59111869A JP 11186984 A JP11186984 A JP 11186984A JP S60254375 A JPS60254375 A JP S60254375A
- Authority
- JP
- Japan
- Prior art keywords
- image
- storage device
- image storage
- bone
- adder
- 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
- 230000006870 function Effects 0.000 claims abstract description 51
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 29
- 238000012545 processing Methods 0.000 claims description 13
- 238000000034 method Methods 0.000 claims description 7
- 230000000873 masking effect Effects 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 8
- 238000007796 conventional method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 210000004872 soft tissue Anatomy 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/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Processing Or Creating Images (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)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
（技術分野）
本発明は、Ｘ線ＣＴ （ＣＴはＣＯｍｐＵｔｅｄＴ’
Ｏｍ ＯＱ ｒａ ｐ １１ Ｖの略）像に生じる骨の
ビームノ＼−１：ニングアーチファクトを再構成像のみ
を使って高速に除去するようにしたＸ線ＣＴ画像処理装
置に関する。（従来技術） 骨のビームハードニング補正は、従来ジョセフ他の文献
（Ｊ ｏｓｅｐｌ＋−８ｐｉｔａｌ−ＡＭ ｅｔｈｏｄ
ｔ’ｏｒＣｏｒｒｅｃｔｉｎｇ ３ｏｎｅ Ｉｎｄｕ
ｃｅｄ Ａｒｔｉｆａｃｔｓ ｉｎＣｏ’ｍｐｕｔｅｄ
Ｔ ｏｍｏａｒａｐｈｙ Ｓ ｃａｎｎｅｒｓ、 Ｊ
ｏｕｒｎａｌ ｏｆＣｏｍｐｕｔｅ＋’ △５ｓｉｓ
ｔｅｔｌ Ｔ−ｏｍｏｇｒａｐｈｙ Ｖｏ１２ｐ１００
〜１ｏ８． Ｊａｎ、１９７Ｂ）に見られるように、再
構成画像を画像上でスキャンしその中に含まれるソフト
ティシュ一部と骨（０丁値で分離される）の経路長から
補正用を泪算し、これをもとに実スキャンデータを補正
した上で再び再構成Ｊ−る方法で行われていた。
Ｏｍ ＯＱ ｒａ ｐ １１ Ｖの略）像に生じる骨の
ビームノ＼−１：ニングアーチファクトを再構成像のみ
を使って高速に除去するようにしたＸ線ＣＴ画像処理装
置に関する。（従来技術） 骨のビームハードニング補正は、従来ジョセフ他の文献
（Ｊ ｏｓｅｐｌ＋−８ｐｉｔａｌ−ＡＭ ｅｔｈｏｄ
ｔ’ｏｒＣｏｒｒｅｃｔｉｎｇ ３ｏｎｅ Ｉｎｄｕ
ｃｅｄ Ａｒｔｉｆａｃｔｓ ｉｎＣｏ’ｍｐｕｔｅｄ
Ｔ ｏｍｏａｒａｐｈｙ Ｓ ｃａｎｎｅｒｓ、 Ｊ
ｏｕｒｎａｌ ｏｆＣｏｍｐｕｔｅ＋’ △５ｓｉｓ
ｔｅｔｌ Ｔ−ｏｍｏｇｒａｐｈｙ Ｖｏ１２ｐ１００
〜１ｏ８． Ｊａｎ、１９７Ｂ）に見られるように、再
構成画像を画像上でスキャンしその中に含まれるソフト
ティシュ一部と骨（０丁値で分離される）の経路長から
補正用を泪算し、これをもとに実スキャンデータを補正
した上で再び再構成Ｊ−る方法で行われていた。
しかしながら、この方法では次のような問題点があった
。
。
０２度の再構成と画像上でのスキャンが必要なため、演
算間が莫人であり、高速処理が期待できない。
算間が莫人であり、高速処理が期待できない。
■実スキャンデータが必要なため事後処理が容易でない
。
。
（発明の目的）
本発明は、このような点に鑑みてなされたもので、その
目的は、アーチファクトの除去効果は損わず、高速かつ
簡単に画像データのみから骨のビームハードニング補正
を行い１ｑるＸ線ＣＴ画像処理装置を提供づ−ることに
ある。
目的は、アーチファクトの除去効果は損わず、高速かつ
簡単に画像データのみから骨のビームハードニング補正
を行い１ｑるＸ線ＣＴ画像処理装置を提供づ−ることに
ある。
（発明の構成）
このような目的を達成する第１の発明は、外部装置に対
しＸＩＩＣＴの再構成濃淡像を入出力することのできる
第１の画像記憶装置と、この第１の画像記憶装置から与
えられる画像を予め定められたＣＴ値をもとに２鎮化す
る２埴化装置と、この２値画像を記憶しておくための第
２の画像記憶装置と、重み関数を蓄えた重み関数記憶装
置と、前記２値画像をマスクとして前記重み関数を加算
器るマスクイ４加算器と、このマスク付加算器で得られ
た加算結果を蓄える第３の画像記ｍ装置と、前記第１の
画像記憶装置と前記第３の画像記憶装置の画像間の加算
を行い、加算結果を第１の画像記憶装置に与える画像加
算器とを具備し、再構成像の骨成分を表づ全画素でこの
点を中心として広がる２次元の重み関数の総和をとって
補正像とし、再構成像と補正像をＪＩＧ算づることによ
り骨のビームハードニング補正を行うようにしたことを
特徴とするもので、第２の発明は、外部装置に対してＸ
線ＣＴの再構成濃淡像を入出力することのできる第１の
画像記憶装置と、この第１の画像記憶装置から与えられ
る画像を縦横ともに縮小する画像縮小装置と、この縮小
画像を記憶する第４の画像記憶装置と、この第４の画像
記憶装置り冒ろ与えられる画像を予め定められたＣＴ値
をもとに２値化する２値化装置と、この２値画像を記憶
しておくための第２の画像記憶装置と、重み関数を蓄え
た重み関数記憶装置と、前記２値画像をマスクとして前
記重み関数を加算するマスク（−１加算器と、このマス
ク付加算器で得られｌ、二加算結果を蓄える第３の画像
記憶装置ど、この第３の画像記憶装置から与えられる画
像を補間を施し”＜縦横ともにもとの大きさに拡大する
画像拡大装置と、前記第１の画像記憶装置と前記画像拡
大装置の画像間の加算を行い、加算結果を第１の画像記
憶装置に与える画像加算器とを具備し、骨を表１２１ａ
１ｉｊｉ像及び重み関数を縦横１／Ｎに縮小した後で
２値画像をマスクとして重み関数を２次元で加算する処
理を骨を表す全画素について繰り返し、得られ／ζ加算
像を補間を使って縦横Ｎ倍に拡大した上でもとの再構成
濶淡像に加算して骨のビームハードニング補正を行うよ
うにしたことを’ｔ？ｉ徴とするもので・ある。
しＸＩＩＣＴの再構成濃淡像を入出力することのできる
第１の画像記憶装置と、この第１の画像記憶装置から与
えられる画像を予め定められたＣＴ値をもとに２鎮化す
る２埴化装置と、この２値画像を記憶しておくための第
２の画像記憶装置と、重み関数を蓄えた重み関数記憶装
置と、前記２値画像をマスクとして前記重み関数を加算
器るマスクイ４加算器と、このマスク付加算器で得られ
た加算結果を蓄える第３の画像記ｍ装置と、前記第１の
画像記憶装置と前記第３の画像記憶装置の画像間の加算
を行い、加算結果を第１の画像記憶装置に与える画像加
算器とを具備し、再構成像の骨成分を表づ全画素でこの
点を中心として広がる２次元の重み関数の総和をとって
補正像とし、再構成像と補正像をＪＩＧ算づることによ
り骨のビームハードニング補正を行うようにしたことを
特徴とするもので、第２の発明は、外部装置に対してＸ
線ＣＴの再構成濃淡像を入出力することのできる第１の
画像記憶装置と、この第１の画像記憶装置から与えられ
る画像を縦横ともに縮小する画像縮小装置と、この縮小
画像を記憶する第４の画像記憶装置と、この第４の画像
記憶装置り冒ろ与えられる画像を予め定められたＣＴ値
をもとに２値化する２値化装置と、この２値画像を記憶
しておくための第２の画像記憶装置と、重み関数を蓄え
た重み関数記憶装置と、前記２値画像をマスクとして前
記重み関数を加算するマスク（−１加算器と、このマス
ク付加算器で得られｌ、二加算結果を蓄える第３の画像
記憶装置ど、この第３の画像記憶装置から与えられる画
像を補間を施し”＜縦横ともにもとの大きさに拡大する
画像拡大装置と、前記第１の画像記憶装置と前記画像拡
大装置の画像間の加算を行い、加算結果を第１の画像記
憶装置に与える画像加算器とを具備し、骨を表１２１ａ
１ｉｊｉ像及び重み関数を縦横１／Ｎに縮小した後で
２値画像をマスクとして重み関数を２次元で加算する処
理を骨を表す全画素について繰り返し、得られ／ζ加算
像を補間を使って縦横Ｎ倍に拡大した上でもとの再構成
濶淡像に加算して骨のビームハードニング補正を行うよ
うにしたことを’ｔ？ｉ徴とするもので・ある。
（実施例）
以下、図面を用いて本発明の実施例を詳しく説明する。
第１図は第１の発明に係るＸ線Ｃ１−画像処理装置の一
実施例を示づ要部のブロック構成図である。この図にお
いて、１は再構成濃淡像を記憶づる第１の画像記憶装置
、２は２俯化装置、３は第２の画像記４＠装置、４は重
み関数記憶装置、５５はマスクイ」加算器、６は第３の
画像記憶装置、７は画像加Ｑ器である。
実施例を示づ要部のブロック構成図である。この図にお
いて、１は再構成濃淡像を記憶づる第１の画像記憶装置
、２は２俯化装置、３は第２の画像記４＠装置、４は重
み関数記憶装置、５５はマスクイ」加算器、６は第３の
画像記憶装置、７は画像加Ｑ器である。
第１の画像記ｆｆｉ装置１は外部装＠（図示せず）に対
し画像を入出力することができるようになっている。２
値化装置２は第１の画像記憶装置１の画像を予め設定さ
れたＣＴ値をしきい値として２値化づるものであり、第
２の画像記憶装置３はこの２値画像を記１ ｔｌるもの
である。マスク付加算器５は重み関数記憶装置４と第２
の画像記憶装置３との出力を受（プ、２値画像上骨を表
すＣＴ値の高い画素位置に重み関数の中心を合わせて２
次元に広がる重み関数を第３の画像記憶装置６に加算記
憶する。画像加睦器７は第３の画像記憶装置６の画像と
第１の画像記憶装置１の画像とを加算するもので、加算
されＩこ画像は第１の画像記憶装置１に記憶されるよう
になっている。
し画像を入出力することができるようになっている。２
値化装置２は第１の画像記憶装置１の画像を予め設定さ
れたＣＴ値をしきい値として２値化づるものであり、第
２の画像記憶装置３はこの２値画像を記１ ｔｌるもの
である。マスク付加算器５は重み関数記憶装置４と第２
の画像記憶装置３との出力を受（プ、２値画像上骨を表
すＣＴ値の高い画素位置に重み関数の中心を合わせて２
次元に広がる重み関数を第３の画像記憶装置６に加算記
憶する。画像加睦器７は第３の画像記憶装置６の画像と
第１の画像記憶装置１の画像とを加算するもので、加算
されＩこ画像は第１の画像記憶装置１に記憶されるよう
になっている。
このような構成において、外部装置より第１の画像記憶
装置１に入力されたＸ線ＣＴの再構成濃淡像は、２値化
装置２において予め設定されたしきい値と比較され２値
化される。この２値画像を第２の画像記憶装置３に順次
記憶してゆく。マスク（＝Ｉ加棹器５では、この第２の
画像記憶装置３より２値画像を受（プ、この２値画像上
骨を表すＣＴ値の高い画素位置に重み関数記憶装置４か
らの重み関数の中心を合わゼて２次元に広がる重み関数
を第３の画像記憶装置６に加算記憶Ｊ−る。尚、第３の
画像記憶装置６の初期記憶値は零に設定されている。こ
の加算処理をしきい値」メ−）−の０丁値を持つすべて
の画素に対して繰り返す。
装置１に入力されたＸ線ＣＴの再構成濃淡像は、２値化
装置２において予め設定されたしきい値と比較され２値
化される。この２値画像を第２の画像記憶装置３に順次
記憶してゆく。マスク（＝Ｉ加棹器５では、この第２の
画像記憶装置３より２値画像を受（プ、この２値画像上
骨を表すＣＴ値の高い画素位置に重み関数記憶装置４か
らの重み関数の中心を合わゼて２次元に広がる重み関数
を第３の画像記憶装置６に加算記憶Ｊ−る。尚、第３の
画像記憶装置６の初期記憶値は零に設定されている。こ
の加算処理をしきい値」メ−）−の０丁値を持つすべて
の画素に対して繰り返す。
このようにしてできあがった画像は、画像加算器７にＪ
−って初めに人力された再構成濃淡像（第１の画像記憶
装置１に記憶されている）に加算される。これにより、
骨のビームハードニングアーヂフｊ・り［・を除去Ｊる
ことができ、その画像は改めて第１の画像記憶装置１に
書き込まれ、必要に応じて外部装置に出力される。
−って初めに人力された再構成濃淡像（第１の画像記憶
装置１に記憶されている）に加算される。これにより、
骨のビームハードニングアーヂフｊ・り［・を除去Ｊる
ことができ、その画像は改めて第１の画像記憶装置１に
書き込まれ、必要に応じて外部装置に出力される。
第２図は本発明の他の実施例を示す説明図である。第１
図と異なるところは第３の画像記憶装置６と画像加算器
７を省略した点である。第２図においては、重み関数の
総和をとって後から２枚の画像間で加算する代りに、重
み関数を再構成像の対応Ｊる位置に直接加算し、結果を
第１の画像記憶装置 構成によっても上述と同様に骨のビームハードニングア
ーチフアクＬ−を除去覆ることができる。
図と異なるところは第３の画像記憶装置６と画像加算器
７を省略した点である。第２図においては、重み関数の
総和をとって後から２枚の画像間で加算する代りに、重
み関数を再構成像の対応Ｊる位置に直接加算し、結果を
第１の画像記憶装置 構成によっても上述と同様に骨のビームハードニングア
ーチフアクＬ−を除去覆ることができる。
尚、あるＣＴ（ａ以上の画素に対して同一の重み関数を
加算しているが、あるＣＴ値以上の画素に対してそのＣ
Ｔ値により変調した重み関数を加算するように構成する
ことも可能である。
加算しているが、あるＣＴ値以上の画素に対してそのＣ
Ｔ値により変調した重み関数を加算するように構成する
ことも可能である。
又、マスク付加算器５は第２の画像記憶装置３の２値画
像と重み関数記憶装置４の重み関数の２次元コンボリュ
ーションて冒ぎ換えることも可能で必る。
像と重み関数記憶装置４の重み関数の２次元コンボリュ
ーションて冒ぎ換えることも可能で必る。
以上）ホベた第１の発明の実施例においては、骨を表す
画素の数だけ２次元の重み関数を加算しており、場合に
よってはかなりの演算量になる。そこで、骨を表１２値
画像１重み関数を縦横１／Ｎに縮小した後で２値画像を
マスクとして重み関数を２次元で加算し、骨を表す全画
素について加算を繰り返すようにすれば、演算量を減少
させることができ、それだけ高速処理を実現することが
できる。
画素の数だけ２次元の重み関数を加算しており、場合に
よってはかなりの演算量になる。そこで、骨を表１２値
画像１重み関数を縦横１／Ｎに縮小した後で２値画像を
マスクとして重み関数を２次元で加算し、骨を表す全画
素について加算を繰り返すようにすれば、演算量を減少
させることができ、それだけ高速処理を実現することが
できる。
第３図はこのような高速処理を図った第２の発明におけ
るＸ線ＣＴ画像処理装置の一実施例を示す構成図である
。第３図において、第１図と同等部分には同一符号を付
しである。第１図と異なる構成部分は画像縮小装置３１
．第４の画像記ｆＩＡ装置３２１画像拡大装Ｒ３３，第
５の画像記憶装置３４である。
るＸ線ＣＴ画像処理装置の一実施例を示す構成図である
。第３図において、第１図と同等部分には同一符号を付
しである。第１図と異なる構成部分は画像縮小装置３１
．第４の画像記ｆＩＡ装置３２１画像拡大装Ｒ３３，第
５の画像記憶装置３４である。
画像縮小装置３１は第１の画像記憶装置１の画像を受け
、これを縦横両方向に１／Ｎに縮小するもので、結果は
第４の画像記憶装置３２に格納される。第４の画像記憶
装置３２の出力は２値化装・ 置２に導かれるようにな
っている。
、これを縦横両方向に１／Ｎに縮小するもので、結果は
第４の画像記憶装置３２に格納される。第４の画像記憶
装置３２の出力は２値化装・ 置２に導かれるようにな
っている。
画像拡大装置３３は第３の画像記憶装置６の縮小画像を
Ｎ倍し、もとの大きさの画像に拡大するものである。拡
大画像は第５の画像記憶装＠３４に一旦格納された後画
像加算器７に送られる。
Ｎ倍し、もとの大きさの画像に拡大するものである。拡
大画像は第５の画像記憶装＠３４に一旦格納された後画
像加算器７に送られる。
このような構成における動作を次に説明づ−る。
外部装置から第１の画像記憶装置１に入ツノされたＸ線
ＣＨの再構成濃淡像は、画像縮小装置３１で縦横両方向
に１／Ｎに縮小され、第４の画像記憶装＠３２に記憶さ
れる。尚、Ｎは通常整数が用いられる。
ＣＨの再構成濃淡像は、画像縮小装置３１で縦横両方向
に１／Ｎに縮小され、第４の画像記憶装＠３２に記憶さ
れる。尚、Ｎは通常整数が用いられる。
記憶装置３２の縮小画像は予め設定されたＣＴ値に従っ
て２値化装置２で２値化され、第２の画像拡大装置３に
記憶される。
て２値化装置２で２値化され、第２の画像拡大装置３に
記憶される。
重み関数記憶装置４０重み関数も前述の縮小率１／Ｎに
対応して予め縦横１／Ｎに縮小されているものとする。
対応して予め縦横１／Ｎに縮小されているものとする。
マスク付加算器５において、重み関数記憶装置４の重み
関数を２値画像（記憶装置３）の骨を族１−各画素に対
し中心を合わＵて第３の画像記憶装置６に加算記憶ブー
る。尚、第３の画像記憶装置６の初期記憶値は零に設定
されている。
関数を２値画像（記憶装置３）の骨を族１−各画素に対
し中心を合わＵて第３の画像記憶装置６に加算記憶ブー
る。尚、第３の画像記憶装置６の初期記憶値は零に設定
されている。
この加算処理を骨を表ねづ全画素について行った結果は
縮小像に対重る補正像となっているので、これを画像拡
大装置３３でＮ倍し、もとの大ぎざに復元する。この拡
大には割算時間の短縮のため線型補間を用い、１次元毎
に２段階に分けて２次元の補間を行う。拡大された画像
は第５の画像記憶装＠３４に蓄えられ、この拡大補正像
は画像加算器７において第１の画像記１８置１上の入力
画像と２次元加算され、その結果が第１の画像記憶装置
１上に出力される。
縮小像に対重る補正像となっているので、これを画像拡
大装置３３でＮ倍し、もとの大ぎざに復元する。この拡
大には割算時間の短縮のため線型補間を用い、１次元毎
に２段階に分けて２次元の補間を行う。拡大された画像
は第５の画像記憶装＠３４に蓄えられ、この拡大補正像
は画像加算器７において第１の画像記１８置１上の入力
画像と２次元加算され、その結果が第１の画像記憶装置
１上に出力される。
第４図はＮ＝４の場合の拡大補間を説明する説明図であ
る。既知の値（○印）をもとにまずＸ印での値をそれを
挟むＱ印の値から次式で算出する。
る。既知の値（○印）をもとにまずＸ印での値をそれを
挟むＱ印の値から次式で算出する。
ｆ （ｉ ＞ ＝ｆ （１）＋（ｉ ／４）・（ｆ （
２）−ｆ （１）） ここで、＝＝１．２．３ 次に＊印での値をそれを挟む左右のＸ印の値から次式で
輝出する。
２）−ｆ （１）） ここで、＝＝１．２．３ 次に＊印での値をそれを挟む左右のＸ印の値から次式で
輝出する。
・＜ｒ ＜２＞−ｒ （１））
ここで、ｒ ＝１．２．３
このようにして補間を行う。Ｎ＝４以外の場合も同様で
ある。
ある。
尚、実施例に限らず、次のようにしでもよい。
０画像縮小装置３１では実際には画像を縮小ぜず、２値
化装冒２と一体化して縦横ともＮ画素毎に１ノ°ンプル
して２値化する。
化装冒２と一体化して縦横ともＮ画素毎に１ノ°ンプル
して２値化する。
■画像拡大装跨３３は第５の画像記憶装置３４を通さず
に直接第１の画像記憶装置と加算Ｊる。
に直接第１の画像記憶装置と加算Ｊる。
■マスク付加算器５を２次元コンポリコーション装置に
向き換える。
向き換える。
０画像縮小装置３１．第４の画像記憶装置３２と第２の
画像記憶装置３を用いずに、２値化装置２が出力づるＮ
＋Ｊンブル旬の値に従って加算器５をオン・Δノする
。
画像記憶装置３を用いずに、２値化装置２が出力づるＮ
＋Ｊンブル旬の値に従って加算器５をオン・Δノする
。
（発明の効宋）
以−１説明したように、本発明によれば、重み関数を適
切に設定することにより、ジ、、１Ｆ：！フ他の従来の
方法によるのと同等の骨のじ一ムハードニング補正がで
き、脳と頭蓋骨の画＋ＩＪ！１のきれが良く、しかも骨
によるカッピングのない画像を１ｑることができる。
切に設定することにより、ジ、、１Ｆ：！フ他の従来の
方法によるのと同等の骨のじ一ムハードニング補正がで
き、脳と頭蓋骨の画＋ＩＪ！１のきれが良く、しかも骨
によるカッピングのない画像を１ｑることができる。
又、従来のような画像−［：で改めてスキャンを繰り返
し、その経路内での成分を調べ、それによる実スキャン
データの補正を行い、それを再び再構成するという複雑
かつ莫大な演算量を要する処理と比べ、再構成像をもど
に簡単な画像間の演算を繰り返すだ１）で処理ができる
など、アルゴリズムの簡単さ、演算量の少なさから極め
て実用的である。
し、その経路内での成分を調べ、それによる実スキャン
データの補正を行い、それを再び再構成するという複雑
かつ莫大な演算量を要する処理と比べ、再構成像をもど
に簡単な画像間の演算を繰り返すだ１）で処理ができる
など、アルゴリズムの簡単さ、演算量の少なさから極め
て実用的である。
例えば、使用画像３２０ｘ ３２０画素画素用重み関数
６１Ｘ ６１画素としたとき、ジ３廿フ他の従来方法で
は約１５時間ががっていたが、本発明によれば約５０分
の処理時間でづむ。
６１Ｘ ６１画素としたとき、ジ３廿フ他の従来方法で
は約１５時間ががっていたが、本発明によれば約５０分
の処理時間でづむ。
又、第３図に示すような本発明によれば、最初に画像を
１／Ｎ２に縮小し、この縮小画像に対して２値化、重み
関数の加算を行うため、匝み関数ももどの１／Ｎ２に縮
小しＩこものを使うことから加筒回数はＭ ２ ×ｆ＜
から＜Ｍ／Ｎ） ２ｘｋ ／Ｎ２と１．／ Ｎ４に減少
する。ここでＭ２は重み関数の面積を画素数で表したも
の、ｋは骨を表す画素の数である。
１／Ｎ２に縮小し、この縮小画像に対して２値化、重み
関数の加算を行うため、匝み関数ももどの１／Ｎ２に縮
小しＩこものを使うことから加筒回数はＭ ２ ×ｆ＜
から＜Ｍ／Ｎ） ２ｘｋ ／Ｎ２と１．／ Ｎ４に減少
する。ここでＭ２は重み関数の面積を画素数で表したも
の、ｋは骨を表す画素の数である。
計ＩＩＩによるシミュレーションににれば、Ｎ＝４、Ｍ
＝６１で画質の劣化は±１ＣＴ値以下であり、４算時間
は縮小しない場合に比べて １／ ２０３倍どｔＪ′つ
Ｉこ。
＝６１で画質の劣化は±１ＣＴ値以下であり、４算時間
は縮小しない場合に比べて １／ ２０３倍どｔＪ′つ
Ｉこ。
この方式は畳のビームハードニング補正像が低い空間周
波数成分しか含んでいないという現象を巧みに利用した
結果のたまものである。
波数成分しか含んでいないという現象を巧みに利用した
結果のたまものである。
第１図は第１の発明に係るＸ線ＣＴ画像処理装置の一実
施例を承り構成図、第２図は第１の発明の他の実施例を
示す説明図、第３図は第２の発明の一実施例を示す構成
図、第４図は拡大補間を説明するための説明図である。 １・・・第１の画像記憶装置 ２・・・２１ｍ化装置 ３・・・第２の画像記憶装置 ４・・・重み関数記憶装置 ５・・・マスク付加算器 ６・・・第３の画像記憶装置 ７・・・画像加算器 ３１・・・画像縮小！！、置 ３２・・・第４の画像記憶装置 ３３・・・画像拡大装置 ３４・・・第５の画像記憶装置 ｐＰ）２図 ｏｘｘｘ。 Ｘ来米来× ×僚１ｇ辛 Ｘ Ｘ米殊来Ｘ
施例を承り構成図、第２図は第１の発明の他の実施例を
示す説明図、第３図は第２の発明の一実施例を示す構成
図、第４図は拡大補間を説明するための説明図である。 １・・・第１の画像記憶装置 ２・・・２１ｍ化装置 ３・・・第２の画像記憶装置 ４・・・重み関数記憶装置 ５・・・マスク付加算器 ６・・・第３の画像記憶装置 ７・・・画像加算器 ３１・・・画像縮小！！、置 ３２・・・第４の画像記憶装置 ３３・・・画像拡大装置 ３４・・・第５の画像記憶装置 ｐＰ）２図 ｏｘｘｘ。 Ｘ来米来× ×僚１ｇ辛 Ｘ Ｘ米殊来Ｘ
Claims (2)
- （１）外部装置に対しＸ線ＣＴの再構成濃淡像を入出力
づることのできる第１の画像記憶装置と、この第１の画
像記憶装置から与えられる画像を予め定められたＣＴ値
をもとに２値化Ｊる２値化装置と、この２値画像を記憶
しており１．：めの第２の画像記憶装置と、重み関数を
蓄えた重み関数記４Ｊａ装置と、前記２値画像をマスク
として前記重み関数を加算するマスク付加算器と、この
マスク付加算器で得られＩζ加算結果を蓄える第３の画
像記憶装置と、前記第１の画像記憶装置と前記第３の画
像記憶装置の画像間の加算を行い、加算結果を第１の画
像記憶装置に与える画像加算器とを具備し、再構成像の
青成分を表す全画素でこの点を中心どして広がる２次元
の重み関数の総和をとって補正像とし、再構成像と補正
像を加算することにより骨のビームハードニング補正を
行うようにしたことを特徴とするＸ￥ＡＣＴ画像処理装
置。 - （２）外部装置に対してＸ線ＣＴの再構成濃淡像を入出
力することのできる第１の画像記憶装置と、この第１の
画像記憶装置から与えられる画像を縦横ともに縮小する
画像縮小装置と、この縮小画像を記憶する第４の画像記
憶装置と、この第４の画像記憶装置から与えられる画像
を予め定められたＣＴ値をもどに２（ｌＩｊ化する２値
化装置と、この２値画像を記憶しておくための第２の画
像記憶装置と、重み関数を蓄えた重み関数記憶装置と、
前記２値画像をマスクとして前記重み関数を加算するマ
スク付加算器と、このマスク付加算器で得られた加算結
果を蓄える第３の画像記憶装置と、この第３のｍ像記憶
装置から与えられる画像を補間を施しで縦横ともにもと
の大きさに拡大する画像拡大装置と、前記第１の画像記
憶装置と前記画像拡大装置の画像間の加算を行い、加算
結果を第１の画像記憶装置に与える画像加算器とを具備
し、骨を表す２値画像及び重み関数を縦横１／Ｎに縮小
した後で２１画像をマスクとして重み関数を２次元で加
算する処理を骨を表す全画素について繰り返し、得られ
た加算像を補間を使って縦横Ｎ倍に拡大した上でもとの
再構成Ｉ淡像に加算して骨のビームハードニング補正を
行うようにしたことを特徴とするＸ線ＣＴｉ１ｉｉ像処
理装置。
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. |
| EP85901557A EP0204844B1 (en) | 1984-03-15 | 1985-03-13 | X-ray ct image processor |
| US06/800,077 US4739481A (en) | 1984-03-15 | 1985-03-13 | X-ray CT image processing apparatus |
| PCT/JP1985/000124 WO1985004090A1 (en) | 1984-03-15 | 1985-03-13 | X-ray ct image processor |
| DE198585901557T DE204844T1 (de) | 1984-03-15 | 1985-03-13 | Bildverarbeitungsgeraet fuer roentgentomographie. |
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 true JPS60254375A (ja) | 1985-12-16 |
| JPH0454455B2 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) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH03144860A (ja) * | 1989-10-31 | 1991-06-20 | Shimadzu Corp | 画像フイルタリング装置 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS538589A (en) * | 1976-07-01 | 1978-01-26 | Emi Ltd | Method of correcting error in radiation absorbing display |
-
1984
- 1984-05-31 JP JP59111869A patent/JPS60254375A/ja active Granted
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS538589A (en) * | 1976-07-01 | 1978-01-26 | Emi Ltd | Method of correcting error in radiation absorbing display |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH03144860A (ja) * | 1989-10-31 | 1991-06-20 | Shimadzu Corp | 画像フイルタリング装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0454455B2 (ja) | 1992-08-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3995854B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| US7072524B1 (en) | Method of and system for image processing and recording medium for carrying out the method | |
| JP3096054B2 (ja) | ハーフトーンデジタル映像から連続トーンデジタル映像への逆変換方法 | |
| US4698843A (en) | Method for compensating for void-defects in images | |
| JPH04306775A (ja) | デジタル画像データの圧縮方法 | |
| JP2001509971A (ja) | フィルタされた影像の中の妨害生成物を最小にする方法及び装置 | |
| JPH0944657A (ja) | 画像処理方法および装置 | |
| JP2004242285A (ja) | ノイズ抑制処理方法および装置並びにプログラム | |
| JP4244094B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| JPS63266982A (ja) | 画像信号の補間方法及びそれを実施する画像信号処理装置 | |
| WO1985004090A1 (en) | X-ray ct image processor | |
| JP3738791B2 (ja) | 画像処理方法および装置 | |
| JP3999432B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| JPH06325165A (ja) | 放射線画像の部分の表示方法 | |
| JP3485454B2 (ja) | 画像階調変換装置、画像階調変更方法およびその方法を実行させるプログラムを記録した媒体ならびに赤外線カメラ | |
| KR101341617B1 (ko) | 단일 영상의 오차모델을 기반으로 한 고해상도 영상 복원장치 및 방법 | |
| JP5514132B2 (ja) | 画像縮小装置、画像拡大装置、及びこれらのプログラム | |
| US7551322B2 (en) | Image edge filtering | |
| JPS60254375A (ja) | Ｘ線ｃｔ画像処理装置 | |
| JP3965302B2 (ja) | 空間フィルタリング処理方法 | |
| WO2011087083A1 (ja) | データ処理方法、データ処理装置、及びデータ処理プログラム | |
| JP5498970B2 (ja) | 画像縮小装置、画像拡大装置、及びこれらのプログラム | |
| JP2002335400A (ja) | 画像処理装置、画像処理方法および記録媒体 | |
| JP2019049865A (ja) | 階調超解像装置及びプログラム | |
| JPH08154188A (ja) | 動画像処理装置 |