JPH04139587A - 画像処理装置 - Google Patents
画像処理装置Info
- Publication number
- JPH04139587A JPH04139587A JP2263411A JP26341190A JPH04139587A JP H04139587 A JPH04139587 A JP H04139587A JP 2263411 A JP2263411 A JP 2263411A JP 26341190 A JP26341190 A JP 26341190A JP H04139587 A JPH04139587 A JP H04139587A
- Authority
- JP
- Japan
- Prior art keywords
- image
- picture
- input
- output
- motion
- 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
- 238000009499 grossing Methods 0.000 claims abstract description 23
- 238000006243 chemical reaction Methods 0.000 claims abstract description 14
- 238000012806 monitoring device Methods 0.000 abstract 1
- 238000000034 method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 238000012886 linear function Methods 0.000 description 3
- 239000008186 active pharmaceutical agent Substances 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000002583 angiography Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Studio Circuits (AREA)
- Closed-Circuit Television Systems (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
この発明は、画質の改善を図るための画像処理装置に関
し、とくにＤＳＡ（デジタルサブトラクションアンギオ
グラフィ）装置やＤＲ（デジタルラジオグラフィ）装置
などの医用ＸＩ！装置で得た画像のノイズを減少させる
ことに好適な画像処理装置に関する。
し、とくにＤＳＡ（デジタルサブトラクションアンギオ
グラフィ）装置やＤＲ（デジタルラジオグラフィ）装置
などの医用ＸＩ！装置で得た画像のノイズを減少させる
ことに好適な画像処理装置に関する。
従来より、ＤＳＡ装置やＤＲ装置などの医用Ｘ線装置で
は、画質改善のためつぎのような２つの方式が提案され
ている。 まず、第１の方式は、入力画像を時間方向にスムージン
グするというものであり、たとえばリカーシブフィルタ
回路がこれに相当する。 第２の方式は、入力画像を時間方向にスムージングする
とともに、原画像との減算を行い、その結果得られた画
像に対して非線形入出力変換を行い、スムージングされ
た画像との減算を行うというものである。
は、画質改善のためつぎのような２つの方式が提案され
ている。 まず、第１の方式は、入力画像を時間方向にスムージン
グするというものであり、たとえばリカーシブフィルタ
回路がこれに相当する。 第２の方式は、入力画像を時間方向にスムージングする
とともに、原画像との減算を行い、その結果得られた画
像に対して非線形入出力変換を行い、スムージングされ
た画像との減算を行うというものである。
しかしながら、従来の第１の方式では画像のＳＺＮ比は
向上するが動きのある被写体（たとえばカテーテル）が
ぼけてしまうという問題がある。 また、従来の第２の方式では、この第１の方式の欠点は
補えるものの、スムージングの度合や非線形人出力テー
ブルについては、被写体の動きに合わせて人間が最適な
状態となるように調整しなくてはならないという煩雑さ
がある。 この発明は、上記に鑑み、上記第２の方式を改善し、被
写体の動きを自動的に検出し予測することによりスムー
ジングの度合と非線形入出力変換特性をつねに最適なも
のとするようリアルタイムに自動調整することができる
、画像処理装置を提供することを目的とする。
向上するが動きのある被写体（たとえばカテーテル）が
ぼけてしまうという問題がある。 また、従来の第２の方式では、この第１の方式の欠点は
補えるものの、スムージングの度合や非線形人出力テー
ブルについては、被写体の動きに合わせて人間が最適な
状態となるように調整しなくてはならないという煩雑さ
がある。 この発明は、上記に鑑み、上記第２の方式を改善し、被
写体の動きを自動的に検出し予測することによりスムー
ジングの度合と非線形入出力変換特性をつねに最適なも
のとするようリアルタイムに自動調整することができる
、画像処理装置を提供することを目的とする。
上記の目的を達成するため、この発明による画像処理装
置においては、実時間でつぎつぎに入力される原画像を
時間方向にスムージングする手段と、このスムージング
後の画像と上記の原画像との間の減算を行う手段と、該
減算結果として得られた画像の各画素値を変換する手段
と、この変換後の画像と」１記スムージング後の画像と
を加算する手段と、上記減算後の画像から画像の動きを
予測し上記スムージングの度合を変化させるとともに上
記の変換特性を変化させる手段とを備えることかｑ青黴
となっている。
置においては、実時間でつぎつぎに入力される原画像を
時間方向にスムージングする手段と、このスムージング
後の画像と上記の原画像との間の減算を行う手段と、該
減算結果として得られた画像の各画素値を変換する手段
と、この変換後の画像と」１記スムージング後の画像と
を加算する手段と、上記減算後の画像から画像の動きを
予測し上記スムージングの度合を変化させるとともに上
記の変換特性を変化させる手段とを備えることかｑ青黴
となっている。
以下、この発明をＸ線ＴＶ画像の画質改善に適用した場
合の一実施例について図面を参照しながら詳細に説明す
る。第１図において、Ａ 、／ Ｄ変換器１に入力され
る信号は、Ｘ線画像のビデオ信号である。すなわち、図
示しないＸ線管からＸ線が被写体に向けて照射され、被
写体を透過したＸＩｉがイメージインテンシファイアに
入射し、イメージインテンシファイアにおいて光学的な
画像に変換され、その光学的な画像がＴＶ右カメラよっ
てビデオ信号に変換され、そのビデオ信号がたとえば３
０フレ一ム／秒でこのＡ／Ｄ変換器］に入力される。 このビデオ信号はＡ／Ｄ変換器１でデジタル信号に変換
された後、時間方向にスムージングする回路、ここでは
りカーシブフィルタ回路２に送られる。このリカーシブ
フィルタ回路２は加算器２１及び２個の入出力変換器２
２．２３で構成されている。これら入出力変換器２２．
２３はたとえばＲＯＭやＲＡＭなどで構成された画素値
の変換回路で、入出力変換器２２では画素値が１／Ｋに
、入出力変換器２３では画素値が１−１／Ｋに変換され
る。このＫはリカーシブフィルタ係数であり、このＫを
変化させるとスムージング効果が変化し、Ｋが大きいと
きスムージングの度合が大きくなる。 つぎに、このリカーシブフィルタ回路２から出力された
画像と原画像とが減算器３で減算される。 その結果、原画像の持つノイズ成分と動き成分とからな
る画像が減算器３から出力されることになる。 この減算器３の出力画像は非線形な入出力変換器４に入
力される。そして加算器５において、この入出力変換器
４の出力画像がリカーシブフィルタ回路２の出力画像と
加算され、その加算後の画像がＤ／Ａ変換器６を経て図
示しないＴＶモニター装置などに出力される。 一方、減算器３からの出力画像は分散及び最大値検出回
路７に送られ、画素値の分散と最大値が検出される。こ
こで、思考のフレーミングなどが行われている状態など
の画像全体に大きな動きがある場合について考えてみる
と、その場合には、この減算器３の出力画像における分
散、最大値とも大となる。このことから、分散が大で最
大値も大であれば、画像全体に大きな動きがあるものと
予測できる。反対に、画像に動きがなければ、分散、最
大値とも小さなものとなる。そのため、分散、最大値と
も小であれば画像全体に動きがないものと予測できる。 その中間の状態として、画像の一部分についてのみ動き
がある場合、たとえばカテーテルが挿入されている状態
の場合、分散は小さい値となるが、最大値は大きくなる
。そこで、分散が小、最大値が大であれば画像の一部に
動きがあると予測できる。なお、分散が大で、最大値が
小の場合は有り得ない。 そこで、分散・最大値を検出することは画像の動きを予
測する二ととなり、それに応じてリカーシブフィルタ回
路２のリカーシブフィルタ係数にと入出力変換器４の入
出力変換特性が変更させられる。すなわち、Ｋは第２図
に示すように分散の関数として変化させられる。入出力
変換特性は、第３図のように分散と最大値との関数とし
て定義させられ、写像され変更されるＭ Ｄ Ｓ （Ｍ
ｏｔｉｏｎｒ）ｅｔｅｃｔｉｏｎ ５ｅｎｓｉｔｉｖｉ
士ｙ ；動き補償感度）に応じて第４図、第５図、第６
図のように変化させられる。 つまり、この実施例では、ＭＤＳが大きければ入出力変
換器４の特性を第４図のように定めることにより減算器
３の出力画像をそのまま加算器５に送るので、結果的に
加算器５からの出力画像として原画像に近い画像を得る
。ＭＤＳか小さい場合は入出力変換器４の特性は第６図
のようになり、減算器３の出力画像は加算器５に送られ
ないこととなって加算器５からの出力画像はリカーシブ
フィルタ回路２の出力画像ということになる。ＭＤＳが
その中間の値のときは、第５図に示すように、入力の中
央のある幅Ｗで出力が０となり、それ以外の入力につい
ては１対１に対応する出力が得られるような入出力特性
となる。したがって、このとき、画素値の中央付近の値
についてはスムージングされた値となり、それ以外の画
素値については原画像の値となる。そして、このＭＤＳ
の値に応じて幅■′か連続的に変化させられる（Ｗが最
大値をとるとき第６図の特性となり、Ｗが最小値をとる
とき第４図の特性となる）こととなる。 辺上述Ｉ＼た、原画像とりカーシブフィルタ回路２の出
力画像との間の減算画像にお（する分散、最大値の各値
と、画像の予測される動きと、Ｋ、ＭＤＳとの関係をと
りまとめるとつぎの表のようになる 分散 最大値 予測される動き ＭＤＳ 大 大 全体の大きな動き 小 大小 大
画像の一部の動き 大 中小 小 動きなし
大 小なお、分散が大で最大値が小の
場合は上記の通り有り得ないため、省略している。 この実施例では、上記の通り分散の値の大小が画像の動
きに対応するため、分散が小さいときは画像の動きが少
ないのでノイズを低減させるためＫの値を大きくしてス
ムージング効果を大きくし、分散が大きいときは画像の
動きが大きいのでＫの値を小さくしてスムージング効果
を小さくする。 これとともに、分散と最大値との関数として定義させら
れたＭＤＳに応じてスムージングされた画素値と原画像
の画素値とを選ぶようにしている。 このような動き予測によるＫの値の変更と入出力変換器
４の特性の変更とを画像のブランキング期間に行ってリ
アルタイムでの自動調整を行うようにしている。 なお、この実施例では分散及び最大値の検出を画像の全
体について行っているが、複数個の代表点を設定し、そ
れらについて分散、最大値の検出を行うようにしてもよ
い。また、Ｋの値は分散に対し２て線形な関数とし、Ｍ
ＤＳの値は最大値に対して線形な関数としたが、これに
限らず非線形な関数とすることももちろん可能である。 さらに、動きを予測するため画像の分散と最大値とを検
出するようにしたが、他の指標を検出して動きを予測す
るようにしてもよい。 【発明の効果］ この発明の画像処理装置によれば、被写体の動きを自動
的に検出してその動きを予測し２、それに応してスムー
ジングの度合と非線形入出力変換特性をリアルタイムに
自動調整しているため、スムージングの度合と非線形入
出力変換特性がつねにＰＬ適なものとされ、被写体がい
かなる動きをしていても、つねにＳ／Ｎ比の優れた且つ
動きによるぼけのない高画質の画像を得ることができる
。
合の一実施例について図面を参照しながら詳細に説明す
る。第１図において、Ａ 、／ Ｄ変換器１に入力され
る信号は、Ｘ線画像のビデオ信号である。すなわち、図
示しないＸ線管からＸ線が被写体に向けて照射され、被
写体を透過したＸＩｉがイメージインテンシファイアに
入射し、イメージインテンシファイアにおいて光学的な
画像に変換され、その光学的な画像がＴＶ右カメラよっ
てビデオ信号に変換され、そのビデオ信号がたとえば３
０フレ一ム／秒でこのＡ／Ｄ変換器］に入力される。 このビデオ信号はＡ／Ｄ変換器１でデジタル信号に変換
された後、時間方向にスムージングする回路、ここでは
りカーシブフィルタ回路２に送られる。このリカーシブ
フィルタ回路２は加算器２１及び２個の入出力変換器２
２．２３で構成されている。これら入出力変換器２２．
２３はたとえばＲＯＭやＲＡＭなどで構成された画素値
の変換回路で、入出力変換器２２では画素値が１／Ｋに
、入出力変換器２３では画素値が１−１／Ｋに変換され
る。このＫはリカーシブフィルタ係数であり、このＫを
変化させるとスムージング効果が変化し、Ｋが大きいと
きスムージングの度合が大きくなる。 つぎに、このリカーシブフィルタ回路２から出力された
画像と原画像とが減算器３で減算される。 その結果、原画像の持つノイズ成分と動き成分とからな
る画像が減算器３から出力されることになる。 この減算器３の出力画像は非線形な入出力変換器４に入
力される。そして加算器５において、この入出力変換器
４の出力画像がリカーシブフィルタ回路２の出力画像と
加算され、その加算後の画像がＤ／Ａ変換器６を経て図
示しないＴＶモニター装置などに出力される。 一方、減算器３からの出力画像は分散及び最大値検出回
路７に送られ、画素値の分散と最大値が検出される。こ
こで、思考のフレーミングなどが行われている状態など
の画像全体に大きな動きがある場合について考えてみる
と、その場合には、この減算器３の出力画像における分
散、最大値とも大となる。このことから、分散が大で最
大値も大であれば、画像全体に大きな動きがあるものと
予測できる。反対に、画像に動きがなければ、分散、最
大値とも小さなものとなる。そのため、分散、最大値と
も小であれば画像全体に動きがないものと予測できる。 その中間の状態として、画像の一部分についてのみ動き
がある場合、たとえばカテーテルが挿入されている状態
の場合、分散は小さい値となるが、最大値は大きくなる
。そこで、分散が小、最大値が大であれば画像の一部に
動きがあると予測できる。なお、分散が大で、最大値が
小の場合は有り得ない。 そこで、分散・最大値を検出することは画像の動きを予
測する二ととなり、それに応じてリカーシブフィルタ回
路２のリカーシブフィルタ係数にと入出力変換器４の入
出力変換特性が変更させられる。すなわち、Ｋは第２図
に示すように分散の関数として変化させられる。入出力
変換特性は、第３図のように分散と最大値との関数とし
て定義させられ、写像され変更されるＭ Ｄ Ｓ （Ｍ
ｏｔｉｏｎｒ）ｅｔｅｃｔｉｏｎ ５ｅｎｓｉｔｉｖｉ
士ｙ ；動き補償感度）に応じて第４図、第５図、第６
図のように変化させられる。 つまり、この実施例では、ＭＤＳが大きければ入出力変
換器４の特性を第４図のように定めることにより減算器
３の出力画像をそのまま加算器５に送るので、結果的に
加算器５からの出力画像として原画像に近い画像を得る
。ＭＤＳか小さい場合は入出力変換器４の特性は第６図
のようになり、減算器３の出力画像は加算器５に送られ
ないこととなって加算器５からの出力画像はリカーシブ
フィルタ回路２の出力画像ということになる。ＭＤＳが
その中間の値のときは、第５図に示すように、入力の中
央のある幅Ｗで出力が０となり、それ以外の入力につい
ては１対１に対応する出力が得られるような入出力特性
となる。したがって、このとき、画素値の中央付近の値
についてはスムージングされた値となり、それ以外の画
素値については原画像の値となる。そして、このＭＤＳ
の値に応じて幅■′か連続的に変化させられる（Ｗが最
大値をとるとき第６図の特性となり、Ｗが最小値をとる
とき第４図の特性となる）こととなる。 辺上述Ｉ＼た、原画像とりカーシブフィルタ回路２の出
力画像との間の減算画像にお（する分散、最大値の各値
と、画像の予測される動きと、Ｋ、ＭＤＳとの関係をと
りまとめるとつぎの表のようになる 分散 最大値 予測される動き ＭＤＳ 大 大 全体の大きな動き 小 大小 大
画像の一部の動き 大 中小 小 動きなし
大 小なお、分散が大で最大値が小の
場合は上記の通り有り得ないため、省略している。 この実施例では、上記の通り分散の値の大小が画像の動
きに対応するため、分散が小さいときは画像の動きが少
ないのでノイズを低減させるためＫの値を大きくしてス
ムージング効果を大きくし、分散が大きいときは画像の
動きが大きいのでＫの値を小さくしてスムージング効果
を小さくする。 これとともに、分散と最大値との関数として定義させら
れたＭＤＳに応じてスムージングされた画素値と原画像
の画素値とを選ぶようにしている。 このような動き予測によるＫの値の変更と入出力変換器
４の特性の変更とを画像のブランキング期間に行ってリ
アルタイムでの自動調整を行うようにしている。 なお、この実施例では分散及び最大値の検出を画像の全
体について行っているが、複数個の代表点を設定し、そ
れらについて分散、最大値の検出を行うようにしてもよ
い。また、Ｋの値は分散に対し２て線形な関数とし、Ｍ
ＤＳの値は最大値に対して線形な関数としたが、これに
限らず非線形な関数とすることももちろん可能である。 さらに、動きを予測するため画像の分散と最大値とを検
出するようにしたが、他の指標を検出して動きを予測す
るようにしてもよい。 【発明の効果］ この発明の画像処理装置によれば、被写体の動きを自動
的に検出してその動きを予測し２、それに応してスムー
ジングの度合と非線形入出力変換特性をリアルタイムに
自動調整しているため、スムージングの度合と非線形入
出力変換特性がつねにＰＬ適なものとされ、被写体がい
かなる動きをしていても、つねにＳ／Ｎ比の優れた且つ
動きによるぼけのない高画質の画像を得ることができる
。
第１図はこの発明の一実施例のブロック図、第２図はり
カーシブフィルタ係数の変化を示すグラフ、第３図はＭ
ＤＳを表すグラフ、第４図、第５図及び第６図は入出力
変換器４の入出力変換特性をそれぞれ示すグラフである
。 〕・・・Ａ／Ｄ変換器、２・・・リカーシブフィルタ回
路、２１．５・・・加算器、２２．２３．４・・・入出
力変換器、３・・・減算器、６・・・Ｄ／Ａ変換器、７
・・・分散及び最大値検出回路。
カーシブフィルタ係数の変化を示すグラフ、第３図はＭ
ＤＳを表すグラフ、第４図、第５図及び第６図は入出力
変換器４の入出力変換特性をそれぞれ示すグラフである
。 〕・・・Ａ／Ｄ変換器、２・・・リカーシブフィルタ回
路、２１．５・・・加算器、２２．２３．４・・・入出
力変換器、３・・・減算器、６・・・Ｄ／Ａ変換器、７
・・・分散及び最大値検出回路。
Claims (1)
- （１）実時間でつぎつぎに入力される原画像を時間方向
にスムージングする手段と、このスムージング後の画像
と上記の原画像との間の減算を行う手段と、該減算結果
として得られた画像の各画素値を変換する手段と、この
変換後の画像と上記スムージング後の画像とを加算する
手段と、上記減算後の画像から画像の動きを予測し上記
スムージングの度合を変化させるとともに上記の変換特
性を変化させる手段とを備えることを特徴とする画像処
理装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2263411A JP2871054B2 (ja) | 1990-09-30 | 1990-09-30 | 画像処理装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2263411A JP2871054B2 (ja) | 1990-09-30 | 1990-09-30 | 画像処理装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH04139587A true JPH04139587A (ja) | 1992-05-13 |
| JP2871054B2 JP2871054B2 (ja) | 1999-03-17 |
Family
ID=17389125
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2263411A Expired - Lifetime JP2871054B2 (ja) | 1990-09-30 | 1990-09-30 | 画像処理装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2871054B2 (ja) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006230512A (ja) * | 2005-02-23 | 2006-09-07 | Hitachi Medical Corp | Ｘ線画像診断装置及び画像処理方法並びにプログラム |
| JP2009297077A (ja) * | 2008-06-10 | 2009-12-24 | Toshiba Corp | 画像処理装置 |
| JP2013048782A (ja) * | 2011-08-31 | 2013-03-14 | Canon Inc | 画像処理装置、画像処理方法、プログラムおよびコンピュータ記録媒体 |
| WO2015052890A1 (ja) * | 2013-10-10 | 2015-04-16 | キヤノン株式会社 | Ｘ線画像処理装置および方法、並びにｘ線撮像装置、プログラム |
-
1990
- 1990-09-30 JP JP2263411A patent/JP2871054B2/ja not_active Expired - Lifetime
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006230512A (ja) * | 2005-02-23 | 2006-09-07 | Hitachi Medical Corp | Ｘ線画像診断装置及び画像処理方法並びにプログラム |
| JP2009297077A (ja) * | 2008-06-10 | 2009-12-24 | Toshiba Corp | 画像処理装置 |
| JP2013048782A (ja) * | 2011-08-31 | 2013-03-14 | Canon Inc | 画像処理装置、画像処理方法、プログラムおよびコンピュータ記録媒体 |
| WO2015052890A1 (ja) * | 2013-10-10 | 2015-04-16 | キヤノン株式会社 | Ｘ線画像処理装置および方法、並びにｘ線撮像装置、プログラム |
| JP2015073800A (ja) * | 2013-10-10 | 2015-04-20 | キヤノン株式会社 | Ｘ線画像処理装置および方法、並びにｘ線撮像装置、プログラム |
| US10016176B2 (en) | 2013-10-10 | 2018-07-10 | Canon Kabushiki Kaisha | X-ray image processing apparatus and method, and X-ray imaging apparatus |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2871054B2 (ja) | 1999-03-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3828251B2 (ja) | 映像のダイナミックレンジ拡大装置 | |
| JPH0411466A (ja) | ノイズリデューサ | |
| JPS63271668A (ja) | Ｘ線画像処理装置 | |
| JP3159465B2 (ja) | 画像表示装置 | |
| JPH04139587A (ja) | 画像処理装置 | |
| JPH03140071A (ja) | 画像表示装置 | |
| JP3410778B2 (ja) | リカーシブフィルタ、ｘ線診断装置 | |
| JP2792133B2 (ja) | 画像処理装置 | |
| JP2885066B2 (ja) | 画像処理装置 | |
| JPH0647035A (ja) | 画像処理装置 | |
| JP3326819B2 (ja) | 画像処理装置 | |
| JPH01273487A (ja) | デジタルｘ線装置 | |
| JP2871066B2 (ja) | Ｘ線画像処理装置 | |
| JP3629875B2 (ja) | Ｘ線画像デジタル処理装置 | |
| JPS62171282A (ja) | 相関適応式雑音低減装置 | |
| JP2594987B2 (ja) | インターラインフリッカー除去回路 | |
| JPH06169906A (ja) | Ｘ線画像診断装置 | |
| JPS61125295A (ja) | テレビジヨン方式変換方式 | |
| JP3326885B2 (ja) | Ｘ線テレビジョン装置 | |
| JPH07154651A (ja) | 画像処理装置 | |
| JPS62285584A (ja) | デジタル画像処理装置 | |
| JPH01277975A (ja) | 画像ノイズリダクシヨン装置 | |
| JPH0669447B2 (ja) | Ｘ線画像処理装置 | |
| JPH03282973A (ja) | Ｘ線画像処理装置 | |
| JP3298445B2 (ja) | Ｘ線画像デジタル処理装置 |