JPH03144860A - 画像フイルタリング装置 - Google Patents
画像フイルタリング装置Info
- Publication number
- JPH03144860A JPH03144860A JP1284007A JP28400789A JPH03144860A JP H03144860 A JPH03144860 A JP H03144860A JP 1284007 A JP1284007 A JP 1284007A JP 28400789 A JP28400789 A JP 28400789A JP H03144860 A JPH03144860 A JP H03144860A
- Authority
- JP
- Japan
- Prior art keywords
- line
- image
- picture element
- picture
- processing
- 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
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
この発明は、たとえばＸ線ＣＴ装置やＭＲＩ装置などで
得た医療の診断用画像に対してフィルタ処理を行う、画
像フィルタリング装置に関する。
得た医療の診断用画像に対してフィルタ処理を行う、画
像フィルタリング装置に関する。
画像にシャープあるいはスムーズ効果をもたらすピクチ
ャーフィルタ処理は、たとえば第２図に示すような各画
素の値を持つ画像マトリクスＭに対して第３図に示すよ
うなテンプレートを作用させるものとして説明すると、
Ｂラインの各画素については Ｂ１＝ａ−ｂ１＋β・（ａｌ＋ｃ１＋ｂｌ−１＋ｂｉ＋
ｘ）＋γ°（ａｉ−ｘ”ａｉ＋ｔ＋ｃｉ−４＋ｃｉ＋ｉ
）の計算を行って新たな画素値Ｂ１を得ることに相当す
る。したがってこの計算を全画素について行う必要があ
り、画像マトリクスＭのサイズがＮ×Ｎであるとすると
、その計算を、ラフ１−ウェアまたは専用ハードウェア
によりほぼＮ２回繰り返すことになる。
ャーフィルタ処理は、たとえば第２図に示すような各画
素の値を持つ画像マトリクスＭに対して第３図に示すよ
うなテンプレートを作用させるものとして説明すると、
Ｂラインの各画素については Ｂ１＝ａ−ｂ１＋β・（ａｌ＋ｃ１＋ｂｌ−１＋ｂｉ＋
ｘ）＋γ°（ａｉ−ｘ”ａｉ＋ｔ＋ｃｉ−４＋ｃｉ＋ｉ
）の計算を行って新たな画素値Ｂ１を得ることに相当す
る。したがってこの計算を全画素について行う必要があ
り、画像マトリクスＭのサイズがＮ×Ｎであるとすると
、その計算を、ラフ１−ウェアまたは専用ハードウェア
によりほぼＮ２回繰り返すことになる。
しかしながら、このように各画素ごとに計算を繰り返す
のでは、画像マトリクスの１辺の画素数の２乗に比例し
て処理数が増え、処理時間が長くかかるという問題があ
る。とくに医療の診断用画像については、近年、画像の
マトリクスサイズは増加する方向にあり、処理スピード
の低下が懸念される。 この発明は、処理回数を減少させて処理スピードを向上
させることができる、画像フィルタリング装置を提供す
ることを目的とする。
のでは、画像マトリクスの１辺の画素数の２乗に比例し
て処理数が増え、処理時間が長くかかるという問題があ
る。とくに医療の診断用画像については、近年、画像の
マトリクスサイズは増加する方向にあり、処理スピード
の低下が懸念される。 この発明は、処理回数を減少させて処理スピードを向上
させることができる、画像フィルタリング装置を提供す
ることを目的とする。
上記目的を達成するため、この発明による画像フィルタ
リング装置においては、画像の各ラインの画素値に対し
てそれぞれテンプレートの各係数を乗じる乗算手段と、
乗算後の各ラインのそれぞれの画素値列をプロファイル
として各画素にパックプロジェクションする手段とを有
することが特徴となっている。
リング装置においては、画像の各ラインの画素値に対し
てそれぞれテンプレートの各係数を乗じる乗算手段と、
乗算後の各ラインのそれぞれの画素値列をプロファイル
として各画素にパックプロジェクションする手段とを有
することが特徴となっている。
画像の各ラインの画素値に対して、テンプレートの各係
数をそれぞれ乗じ、得られた各ラインごとの画素値列を
並べてみると、その同一カラムの値同士を加算した加算
値は、上記テンプレートによりピクチャーフィルタ処理
して得た値に相当する。 そこで、乗算後の各ラインのそれぞれの画素値列をプロ
ファイルとして各画素にバックプロジェクションすれば
、１ライン分のピクチャーフィルタ処理が１度に行える
ことになり、処理回数を減少させ、全体の処理スピード
を大幅に向上させることができる。
数をそれぞれ乗じ、得られた各ラインごとの画素値列を
並べてみると、その同一カラムの値同士を加算した加算
値は、上記テンプレートによりピクチャーフィルタ処理
して得た値に相当する。 そこで、乗算後の各ラインのそれぞれの画素値列をプロ
ファイルとして各画素にバックプロジェクションすれば
、１ライン分のピクチャーフィルタ処理が１度に行える
ことになり、処理回数を減少させ、全体の処理スピード
を大幅に向上させることができる。
つぎにこの発明の一実施例について図面を参照しながら
説明する。第１図において、原画像の各画素の値につき
、ラインごとに乗算器１によりα、β、γの係数が乗算
される。ここで、原画像は第２図で示すような、ＮＸＮ
のマトリクスサイズの画像マトリクスＭで表現されてい
るものとする。 そして、ピクチャーフィルタとして作用させるテンプレ
ート力（第３図のようであるとする。この場合、第４図
に示すように、Ａラインの各画素値ａｌ＋ａ ２． ａ
３．・・・についてはβを乗算したものと、γを乗算
したものとを得る。Ｂラインの各画素値ｂ１゜ｂ ２．
ｂ ３．・・・についてはβを乗算しｆＳものと、α
を乗算したものとをえる。また、Ｃラインの各画素値Ｃ
Ｉ＋ Ｃ２＋ Ｃ３，・・・についてはγを乗算したも
のと、βを乗算したものとを得る。 こうしてα、入β、γを乗算した画素値の列が各ライン
ごとに得られる。その乗算後の各ラインごとの画素値の
列を、第４図のように、左右方向（カラム方向〉に１つ
ずつずらしながら並べて９つの画素値列Ｐ １＋ Ｐ２
＋ Ｐ ３＋・・・、Ｐ９を得る。そしてこのように並
べた画素値列の同一コラムの値をＰｌからＰ９まで加算
する（Ｍ１方向に加算する〉と、その加算後の値は上記
の式で示されるピクチャーフィルタ処理後の画素値Ｂ
ｉ （Ｂ２．Ｂ３．Ｂ４．・・・）に相当することが分
かる。 そこで、このような９つの画素値列Ｐ１〜Ｐ、をバッフ
ァ２に形成し、これら画素値列をバックプロジェクショ
ンすべきプロファイルとしてパックプロジェクタ３に送
る（なお、画素値列Ｐ□とＰ３、Ｐ４とＰ６、Ｐ７とＰ
８とはまったく同じものであり、単にカラム方向の位置
が違えであるに過ぎないので、第４図のように並べたこ
とと同等の関係が得られるならば、必ずしもバッファ２
において実際に並べる必要はない〉。 パックプロジェクタ３では、通常、第５図（イ）で示す
ようなバックプロジェクションが行われる。 すなわち、角度δの方向に分布しているデータであるプ
ロファイルＰを、矢印で示すように、その分布方向に直
角な方向に、画像マトリクスＭの各画素（に相当するメ
モリの記憶区分）に逆投影するのである。そして、この
パックプロジェクタは、通常、Ｘ線ＣＴ装置やＭＲＩ装
置などでプロファイルデータから原画像を再構成するた
めに用いられ、その場合は３６０°の各方向でのプロフ
ァイルデータを画像マトリクスＭの各画素につぎつぎに
逆投影することになる。 ここでは、上記のようにして得た９つの画素値列Ｐ１〜
Ｐ９を、第５図（ロ）に示すように、δ＝Ｏのプロファ
イルＰとしてＢラインの各画素にバックプロジェクショ
ンする。すると、Ｂラインのすべての画素についてのピ
クチャーフィルタ処理が終了したことになる。パックプ
ロジェクタ３が、この９つの画素値列Ｐ、〜Ｐ９のバッ
クプロジェクションを一度に行うことができるパックプ
ロジェクタである場合には、１度の処理で１ライン分の
ピクチャーフィルタ処理ができるため、マトリクスサイ
ズがＮＸＮのとき、Ｎ回の処理で全画像のピクチャーフ
ィルタ処理が行えることになる。あるいは、第４図のよ
うに並べた関係としておいて縦方向に加算したものをバ
ックブロジクションする場合も同様に１ライン分の処理
が１度の処理でできる。 パックプロジェクタは一般に画像再構成手段として用い
られており、ソフトウェアあるいは処理速度を高めるた
め専用ハードウェアとして構成され、Ｘ線ＣＴ装置やＭ
ＲＩ装置では必ず備えられている。そこで、この画像再
構成用のパックプロジェクタを画像フィルタリングのた
めに流用することとすれば、なんらの特別な構成の追加
も不要どなる。また、画像再構成用のパックプロジェク
タに画像フィルタリング機能を受は持たせることにより
、本来の画像処理ハードウェアに別の画像処理を独立に
行わせることもできるようになる。 さらに、たとえばＸ線ＣＴ装置やＭＲＩ装置において、
スキャンによりプロファイルデータを収集した後パック
プロジェクタにより画像再構成を行い、得られた画像に
対して再びパックプロジェクタを使用してピクチャーフ
ィルタ処理を行い、処理後の画像を磁気ディスク装置な
どに格納するようなシステムを構成すれば、スキャン中
は画像処理ハードウェアが切り離された格好となり、磁
気ディスク装置などから他の画像を読み出して処理する
こともできる。 なお、上記では乗算器１、バッファ２、パックプロジェ
クタ３はハードウェア的に表現しているが、ソフトウェ
アでもハードウェアでも構成可能であることはもちろん
である。
説明する。第１図において、原画像の各画素の値につき
、ラインごとに乗算器１によりα、β、γの係数が乗算
される。ここで、原画像は第２図で示すような、ＮＸＮ
のマトリクスサイズの画像マトリクスＭで表現されてい
るものとする。 そして、ピクチャーフィルタとして作用させるテンプレ
ート力（第３図のようであるとする。この場合、第４図
に示すように、Ａラインの各画素値ａｌ＋ａ ２． ａ
３．・・・についてはβを乗算したものと、γを乗算
したものとを得る。Ｂラインの各画素値ｂ１゜ｂ ２．
ｂ ３．・・・についてはβを乗算しｆＳものと、α
を乗算したものとをえる。また、Ｃラインの各画素値Ｃ
Ｉ＋ Ｃ２＋ Ｃ３，・・・についてはγを乗算したも
のと、βを乗算したものとを得る。 こうしてα、入β、γを乗算した画素値の列が各ライン
ごとに得られる。その乗算後の各ラインごとの画素値の
列を、第４図のように、左右方向（カラム方向〉に１つ
ずつずらしながら並べて９つの画素値列Ｐ １＋ Ｐ２
＋ Ｐ ３＋・・・、Ｐ９を得る。そしてこのように並
べた画素値列の同一コラムの値をＰｌからＰ９まで加算
する（Ｍ１方向に加算する〉と、その加算後の値は上記
の式で示されるピクチャーフィルタ処理後の画素値Ｂ
ｉ （Ｂ２．Ｂ３．Ｂ４．・・・）に相当することが分
かる。 そこで、このような９つの画素値列Ｐ１〜Ｐ、をバッフ
ァ２に形成し、これら画素値列をバックプロジェクショ
ンすべきプロファイルとしてパックプロジェクタ３に送
る（なお、画素値列Ｐ□とＰ３、Ｐ４とＰ６、Ｐ７とＰ
８とはまったく同じものであり、単にカラム方向の位置
が違えであるに過ぎないので、第４図のように並べたこ
とと同等の関係が得られるならば、必ずしもバッファ２
において実際に並べる必要はない〉。 パックプロジェクタ３では、通常、第５図（イ）で示す
ようなバックプロジェクションが行われる。 すなわち、角度δの方向に分布しているデータであるプ
ロファイルＰを、矢印で示すように、その分布方向に直
角な方向に、画像マトリクスＭの各画素（に相当するメ
モリの記憶区分）に逆投影するのである。そして、この
パックプロジェクタは、通常、Ｘ線ＣＴ装置やＭＲＩ装
置などでプロファイルデータから原画像を再構成するた
めに用いられ、その場合は３６０°の各方向でのプロフ
ァイルデータを画像マトリクスＭの各画素につぎつぎに
逆投影することになる。 ここでは、上記のようにして得た９つの画素値列Ｐ１〜
Ｐ９を、第５図（ロ）に示すように、δ＝Ｏのプロファ
イルＰとしてＢラインの各画素にバックプロジェクショ
ンする。すると、Ｂラインのすべての画素についてのピ
クチャーフィルタ処理が終了したことになる。パックプ
ロジェクタ３が、この９つの画素値列Ｐ、〜Ｐ９のバッ
クプロジェクションを一度に行うことができるパックプ
ロジェクタである場合には、１度の処理で１ライン分の
ピクチャーフィルタ処理ができるため、マトリクスサイ
ズがＮＸＮのとき、Ｎ回の処理で全画像のピクチャーフ
ィルタ処理が行えることになる。あるいは、第４図のよ
うに並べた関係としておいて縦方向に加算したものをバ
ックブロジクションする場合も同様に１ライン分の処理
が１度の処理でできる。 パックプロジェクタは一般に画像再構成手段として用い
られており、ソフトウェアあるいは処理速度を高めるた
め専用ハードウェアとして構成され、Ｘ線ＣＴ装置やＭ
ＲＩ装置では必ず備えられている。そこで、この画像再
構成用のパックプロジェクタを画像フィルタリングのた
めに流用することとすれば、なんらの特別な構成の追加
も不要どなる。また、画像再構成用のパックプロジェク
タに画像フィルタリング機能を受は持たせることにより
、本来の画像処理ハードウェアに別の画像処理を独立に
行わせることもできるようになる。 さらに、たとえばＸ線ＣＴ装置やＭＲＩ装置において、
スキャンによりプロファイルデータを収集した後パック
プロジェクタにより画像再構成を行い、得られた画像に
対して再びパックプロジェクタを使用してピクチャーフ
ィルタ処理を行い、処理後の画像を磁気ディスク装置な
どに格納するようなシステムを構成すれば、スキャン中
は画像処理ハードウェアが切り離された格好となり、磁
気ディスク装置などから他の画像を読み出して処理する
こともできる。 なお、上記では乗算器１、バッファ２、パックプロジェ
クタ３はハードウェア的に表現しているが、ソフトウェ
アでもハードウェアでも構成可能であることはもちろん
である。
この発明の画像フィルタリング装置によれば、フィルタ
リング処理回数が画像マトリクスの１辺の画素数に比例
する程度に減少でき、処理速度を高めることが可能とな
る。
リング処理回数が画像マトリクスの１辺の画素数に比例
する程度に減少でき、処理速度を高めることが可能とな
る。
第１図はこの発明の一実施例にかかる画像フィルタリン
グ装置を示すブロック図、第２図は原画像の各画素の値
を示す図、第３図はピクチャーフィルタとして作用させ
るテンプレートを示す図、第４図はパックプロジェクシ
ョンすべきプロファイルとして形成されるデータ列を示
す図、第５図（イ）、（ロ）はパックプロジェクション
を説明するための図である。 １・・・乗算器、２・・・バッファ、３・・・パックプ
ロジェクタ、Ｍ・・・画像マトリクス、Ｐ・・・プロフ
ァイル。
グ装置を示すブロック図、第２図は原画像の各画素の値
を示す図、第３図はピクチャーフィルタとして作用させ
るテンプレートを示す図、第４図はパックプロジェクシ
ョンすべきプロファイルとして形成されるデータ列を示
す図、第５図（イ）、（ロ）はパックプロジェクション
を説明するための図である。 １・・・乗算器、２・・・バッファ、３・・・パックプ
ロジェクタ、Ｍ・・・画像マトリクス、Ｐ・・・プロフ
ァイル。
Claims (1)
- （１）画像の各ラインの画素値に対してそれぞれテンプ
レートの各係数を乗じる乗算手段と、乗算後の各ライン
のそれぞれの画素値列をプロファイルとして各画素にバ
ックプロジェクションする手段とを有することを特徴と
する画像フィルタリング装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1284007A JPH0644280B2 (ja) | 1989-10-31 | 1989-10-31 | 画像フイルタリング装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1284007A JPH0644280B2 (ja) | 1989-10-31 | 1989-10-31 | 画像フイルタリング装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH03144860A true JPH03144860A (ja) | 1991-06-20 |
| JPH0644280B2 JPH0644280B2 (ja) | 1994-06-08 |
Family
ID=17673089
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1284007A Expired - Fee Related JPH0644280B2 (ja) | 1989-10-31 | 1989-10-31 | 画像フイルタリング装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0644280B2 (ja) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS60254375A (ja) * | 1984-05-31 | 1985-12-16 | Yokogawa Medical Syst Ltd | Ｘ線ｃｔ画像処理装置 |
-
1989
- 1989-10-31 JP JP1284007A patent/JPH0644280B2/ja not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS60254375A (ja) * | 1984-05-31 | 1985-12-16 | Yokogawa Medical Syst Ltd | Ｘ線ｃｔ画像処理装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0644280B2 (ja) | 1994-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3708434B2 (ja) | 反復円錐ビームｃｔ画像再構成 | |
| US5815421A (en) | Method for transposing a two-dimensional array | |
| DE3726595A1 (de) | Vorrichtung zur korrektur von zwischen verschiedenen ansichten stattfindender bewegung in einem abbildungssystem sowie verfahren dafuer | |
| JP2597612B2 (ja) | 画像圧縮装置 | |
| JPS63256843A (ja) | 角度制限ｘ線計算機式断層写真法において物体の境界情報を求める方法及び角度制限ｘ線計算機式断層写真法システムにおいて高品質の像を再生する方法 | |
| JPH0123136B2 (ja) | ||
| JP2008306298A (ja) | 画像処理装置、画像処理方法、及びプログラム | |
| JPH069061B2 (ja) | 画像デ−タの平滑化方法 | |
| EP1264275A2 (en) | Fast hierarchical reprojection algorithm for tomography | |
| JPH03267879A (ja) | スカラー・データ処理方式 | |
| CN110321089B (zh) | 一种全息体视图打印方法及系统 | |
| JPH03144860A (ja) | 画像フイルタリング装置 | |
| JP3006722B2 (ja) | コンピュータ断層撮影装置 | |
| JPS61283967A (ja) | 画像記録再生装置 | |
| CN114627161A (zh) | 源图片到驱动图片的运动迁移的方法和系统 | |
| JPH071310B2 (ja) | デ−タ収集回路 | |
| JPS58206726A (ja) | 画像処理装置 | |
| Galt et al. | Filtering in frequency space | |
| JPH0581951B2 (ja) | ||
| CN111476743A (zh) | 一种基于分数阶微分的数字信号滤波与图像处理方法 | |
| JPS59155237A (ja) | Ｃｔ画像処理方式 | |
| JPS6219973A (ja) | デ−タ圧縮方法 | |
| JPS6272072A (ja) | 画像処理装置 | |
| JPH0231283A (ja) | 画像再構成装置 | |
| JPH02299370A (ja) | データのコード化 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |