JPS61267875A - 画像処理装置 - Google Patents
画像処理装置Info
- Publication number
- JPS61267875A JPS61267875A JP61014395A JP1439586A JPS61267875A JP S61267875 A JPS61267875 A JP S61267875A JP 61014395 A JP61014395 A JP 61014395A JP 1439586 A JP1439586 A JP 1439586A JP S61267875 A JPS61267875 A JP S61267875A
- Authority
- JP
- Japan
- Prior art keywords
- input
- output
- image
- dimensional
- value
- 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
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/14—Picture signal circuitry for video frequency region
- H04N5/21—Circuitry for suppressing or minimising disturbance, e.g. moiré or halo
- H04N5/213—Circuitry for suppressing or minimising impulsive noise
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Processing Or Creating Images (AREA)
- Closed-Circuit Television Systems (AREA)
- Picture Signal Circuits (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔産業上の利用分野〕
この発明は、特に工業用、医用Ｘ線テレビに使用され、
ディジタル化された画像信号を処理し、画像の改善を行
なう画像処理装置に関するものである。
ディジタル化された画像信号を処理し、画像の改善を行
なう画像処理装置に関するものである。
工業用、医用Ｘ線テレビにおいて、放射線画像の統計雑
音を低減させるための従来の画像処理装置として第６図
に示すものがあった。図において、１０はＭＸＮ個の要
素からなる画像メモリ、３０は加算器である。
音を低減させるための従来の画像処理装置として第６図
に示すものがあった。図において、１０はＭＸＮ個の要
素からなる画像メモリ、３０は加算器である。
このような構成になる画像処理装置において、時刻ｔｉ
におけるディジタル化された入力画像信号Ｉ ｉ （
Ｘｍ、Ｙｎ）と時刻ｔ２−１における画像メモリ （２
次元メモリ）１０の値ＣＪ、−，（Ｘｍ、 Ｙｎ）が加
算器３０で加算された後、画像メモリ１０に送られ、時
刻ｔｉにおける画像メモリの値Ｃｔ（Ｘｍ、Ｙｎ）とな
る。（Ｘｍ、Ｙｎ）は画素の位置であり、上記演算は各
画素毎に順次行なわれ、Δｔ＝ｔｉ−ｔｌ−１の時間の
間にＭＸＮ個の全画素の画素毎の加算と画像メモリ１０
の更新が行なわれる。Δｔは通常テレビレートの３ Ｑ
ｍ ｓであることが多い。そして画像メモリ１０の内
容は処理後の画像信号として出力される。
におけるディジタル化された入力画像信号Ｉ ｉ （
Ｘｍ、Ｙｎ）と時刻ｔ２−１における画像メモリ （２
次元メモリ）１０の値ＣＪ、−，（Ｘｍ、 Ｙｎ）が加
算器３０で加算された後、画像メモリ１０に送られ、時
刻ｔｉにおける画像メモリの値Ｃｔ（Ｘｍ、Ｙｎ）とな
る。（Ｘｍ、Ｙｎ）は画素の位置であり、上記演算は各
画素毎に順次行なわれ、Δｔ＝ｔｉ−ｔｌ−１の時間の
間にＭＸＮ個の全画素の画素毎の加算と画像メモリ１０
の更新が行なわれる。Δｔは通常テレビレートの３ Ｑ
ｍ ｓであることが多い。そして画像メモリ１０の内
容は処理後の画像信号として出力される。
次に動作についてより詳細に説明する。画像信号は、本
来その映像化処理の過程において統計的性質に基く量子
雑音を含むため、信号強度が弱いときには量子雑音の影
響が大きくなり・、コントラスト、分解能等を劣化させ
る。そこで静止画像の場合、画像信号をディジタル化し
、各画素毎に積算して量子雑音の影響を低減するような
処理がなされる。第６図において、時刻ｔｉにおける画
像メモリｌＯの値Ｃｉ （Ｘｍ、Ｙｎ）は、明らかに
それまでの画像信号Ｉ ｊ （Ｘｍ、 Ｙｎ）の積算
Σ／＝ｌ Ｉ ｊ （Ｘｍ、Ｙｎ）となっている。従ってこの画
像メモリ１０の内容を読み出し、作成した処理後の画像
においては、上記加算回数をＮとすると、量子雑音によ
る相対統計誤差は５倍改善されることになる。
来その映像化処理の過程において統計的性質に基く量子
雑音を含むため、信号強度が弱いときには量子雑音の影
響が大きくなり・、コントラスト、分解能等を劣化させ
る。そこで静止画像の場合、画像信号をディジタル化し
、各画素毎に積算して量子雑音の影響を低減するような
処理がなされる。第６図において、時刻ｔｉにおける画
像メモリｌＯの値Ｃｉ （Ｘｍ、Ｙｎ）は、明らかに
それまでの画像信号Ｉ ｊ （Ｘｍ、 Ｙｎ）の積算
Σ／＝ｌ Ｉ ｊ （Ｘｍ、Ｙｎ）となっている。従ってこの画
像メモリ１０の内容を読み出し、作成した処理後の画像
においては、上記加算回数をＮとすると、量子雑音によ
る相対統計誤差は５倍改善されることになる。
この画像処理装置の１ライン分の入出力信号を第７図（
ａ）、 （ｂ）に示す。出力信号は各画素毎に入力信号
のＮ倍となっており、統計誤差は５倍改善されている。
ａ）、 （ｂ）に示す。出力信号は各画素毎に入力信号
のＮ倍となっており、統計誤差は５倍改善されている。
ところが、この図では一般的なシェーディング特性の場
合を示しているが、画像上の信号強度は場所により異な
るので、処理後においても統計誤差の大きさは画素によ
り異なり、画面全体で一様でない。このため、図の場合
には周辺部の統計誤差は中央部の統計誤差より大きく、
図の凸状信号の上に重畳された信号Ｓ１及びＳ２は、凸
状信号に対する相対強度（コントラスト）が同じであっ
ても、中央にある８１信号は判別され、周辺にあるＳ２
信号は判別されないことがある。
合を示しているが、画像上の信号強度は場所により異な
るので、処理後においても統計誤差の大きさは画素によ
り異なり、画面全体で一様でない。このため、図の場合
には周辺部の統計誤差は中央部の統計誤差より大きく、
図の凸状信号の上に重畳された信号Ｓ１及びＳ２は、凸
状信号に対する相対強度（コントラスト）が同じであっ
ても、中央にある８１信号は判別され、周辺にあるＳ２
信号は判別されないことがある。
また十分積算され、周囲の統計的量子雑電が小さくなり
８２信号も判別可能な場合にも、この積算画像をモニタ
ー画面上に表示した場合、人は肉眼で８１信号と８２信
号を同時には視認できないことがある。これは人の視覚
系のダイナミックレンジが狭いためであり、積算された
画像の明るい部分と暗い部分で明度差が太き（なりすぎ
ていることによる。
８２信号も判別可能な場合にも、この積算画像をモニタ
ー画面上に表示した場合、人は肉眼で８１信号と８２信
号を同時には視認できないことがある。これは人の視覚
系のダイナミックレンジが狭いためであり、積算された
画像の明るい部分と暗い部分で明度差が太き（なりすぎ
ていることによる。
従来の画像処理装置は以上のように構成されているので
、積算された画像の量子雑音の低減効果は画素毎に異な
り、同じコントラストを持つ信号であっても、周辺信号
の強度によって判別できたりできなかったりすることが
ある。また積算された画像をモニター画面上に表示した
場合、画像の明部と暗部の明度差が大きくなっているの
で、ダイナミックレンジが限られている人の視覚にとっ
ては、量子雑音から判別されている信号も視認されない
などの欠点があった。
、積算された画像の量子雑音の低減効果は画素毎に異な
り、同じコントラストを持つ信号であっても、周辺信号
の強度によって判別できたりできなかったりすることが
ある。また積算された画像をモニター画面上に表示した
場合、画像の明部と暗部の明度差が大きくなっているの
で、ダイナミックレンジが限られている人の視覚にとっ
ては、量子雑音から判別されている信号も視認されない
などの欠点があった。
この発明は上記のような従来のものの欠点を除去するた
めになされたもので、積算処理された画像は全画素にわ
たり同じ相対的大きさの量子雑音を持ち、またこの画像
をモニター画面上に表示した場合、入力画像における明
部の範囲を広くすることができ、小さな信号を肉眼で視
認できる画像処理装置を提供するものである。
めになされたもので、積算処理された画像は全画素にわ
たり同じ相対的大きさの量子雑音を持ち、またこの画像
をモニター画面上に表示した場合、入力画像における明
部の範囲を広くすることができ、小さな信号を肉眼で視
認できる画像処理装置を提供するものである。
本願の第１の発明に係る画像処理装置は、画像信号を画
素毎に逐次積算する２次元積算器と、この２次元積算器
の出力を指数変換する指数関数器と、この指数関数器の
出力と入力画像信号との差を対応する画素毎にとる減算
器とを設け、上記減算器の出力を上記２次元積算器の対
応する要素に入力するようにしたものである。
素毎に逐次積算する２次元積算器と、この２次元積算器
の出力を指数変換する指数関数器と、この指数関数器の
出力と入力画像信号との差を対応する画素毎にとる減算
器とを設け、上記減算器の出力を上記２次元積算器の対
応する要素に入力するようにしたものである。
また本願の第２の発明に係る画像処理装置は、ディジタ
ル画像信号を逐次記憶する２次元記憶器と、この２次元
記憶器の出力からその指数変換値を減算する関数変換器
と、この関数変換器の出力と入力画像信号との和を対応
する画素毎にとる加算器とを設け、上記加算器の出力を
上記２次元記憶器の対応する要素に入力するようにした
ものである。
ル画像信号を逐次記憶する２次元記憶器と、この２次元
記憶器の出力からその指数変換値を減算する関数変換器
と、この関数変換器の出力と入力画像信号との和を対応
する画素毎にとる加算器とを設け、上記加算器の出力を
上記２次元記憶器の対応する要素に入力するようにした
ものである。
第１の発明においては、画像信号の画素毎に２次元積算
器の出力を指数関数器を通してフィードバックし、これ
と対応する次フレーム画素の入力画像信号との差をとり
、これを上記２次元積算器の対応する要素に入力するル
ープを作成し、画面のどの位置についても出力画像の量
子雑音の相対的標準偏差を一定にする。
器の出力を指数関数器を通してフィードバックし、これ
と対応する次フレーム画素の入力画像信号との差をとり
、これを上記２次元積算器の対応する要素に入力するル
ープを作成し、画面のどの位置についても出力画像の量
子雑音の相対的標準偏差を一定にする。
また第２の発明においては、画像信号の画素毎に２次元
記憶器の出力を関数変換器を通してフィードバックし、
これと対応する次フレーム画素の入力画像信号との和を
とり、これを上記２次元記憶器の対応する要素に入力す
るループを作成し、画面のどの位置についても出力画像
の量子雑音の相対的標準偏差を一定にする。
記憶器の出力を関数変換器を通してフィードバックし、
これと対応する次フレーム画素の入力画像信号との和を
とり、これを上記２次元記憶器の対応する要素に入力す
るループを作成し、画面のどの位置についても出力画像
の量子雑音の相対的標準偏差を一定にする。
以下、本願の第１の発明の一実施例を図について説明す
る。第１図において、１は（ＭＸＮ）個の要素を持つ２
次元積算器、２は２次元積算器１の出力Ｃｔ （Ｘｍ
、Ｙｎ）を入力とする指数関数器、３はディジタル化さ
れ時系列的に送り込まれる入力画像信号１ ｉ （Ｘ
ｍ、Ｙｎ）を正入力とし、入力画像信号の画素（Ｘｍ、
Ｙｎ）に対応する２次元積算器出力Ｃ１（Ｘｍ、Ｙｎ）
の指数関数器出力であるｂ−ｅｘｐ （ａ−ｃｉ
（Ｘｍ、Ｙｎ））を負入力とし、両者の減算を行なう減
算器であり、この減算器出力が２次元積算器１に入力さ
れその対応する要素に加えられてループを形成している
。
る。第１図において、１は（ＭＸＮ）個の要素を持つ２
次元積算器、２は２次元積算器１の出力Ｃｔ （Ｘｍ
、Ｙｎ）を入力とする指数関数器、３はディジタル化さ
れ時系列的に送り込まれる入力画像信号１ ｉ （Ｘ
ｍ、Ｙｎ）を正入力とし、入力画像信号の画素（Ｘｍ、
Ｙｎ）に対応する２次元積算器出力Ｃ１（Ｘｍ、Ｙｎ）
の指数関数器出力であるｂ−ｅｘｐ （ａ−ｃｉ
（Ｘｍ、Ｙｎ））を負入力とし、両者の減算を行なう減
算器であり、この減算器出力が２次元積算器１に入力さ
れその対応する要素に加えられてループを形成している
。
一画面の要素は（Ｍ Ｘ Ｎ）個でできており、これら
全画素の１サイクル分の処理がテレビレートの３０ ｍ
ｓで実行されるようになっている。 Ｍ−Ｎ＝５１２
の場合を考えると、一画素につき３０ ｍ ｓ／ （５
１２Ｘ５１２） ′、１１５ｎｓ以下で１サイクル処理
を行なう必要があるが、指数関数器２として予め指数関
数のテーブルをメモリに書込んでおき、入力に応じてテ
ーブルを参照し出力するようにしておけば、スタティッ
クＲＡＭのアクセスタイム４Ｏｎｓ以内であり、十分可
能である。
全画素の１サイクル分の処理がテレビレートの３０ ｍ
ｓで実行されるようになっている。 Ｍ−Ｎ＝５１２
の場合を考えると、一画素につき３０ ｍ ｓ／ （５
１２Ｘ５１２） ′、１１５ｎｓ以下で１サイクル処理
を行なう必要があるが、指数関数器２として予め指数関
数のテーブルをメモリに書込んでおき、入力に応じてテ
ーブルを参照し出力するようにしておけば、スタティッ
クＲＡＭのアクセスタイム４Ｏｎｓ以内であり、十分可
能である。
次に動作について説明する。
入力画像として量子雑音の大きい同じ画面の像（例えば
放射線画像）が順次送り込まれているとする。ある数の
画面が入力され平衡に達した時、Ｉ ｉ （Ｘｍ、Ｙ
ｎ） ＝ｂ −ｅｘｐ （ａ−Ｃｉ （Ｘｍ、Ｙｎ））
・−（１）となり、Ｃｉ （Ｘｍ、Ｙｎ）がほぼ一定
値Ｃ（Ｘｍ、 Ｙｎ）となるので、この系の入出力特性
は式％式％ となり、対数特性を持つ。
放射線画像）が順次送り込まれているとする。ある数の
画面が入力され平衡に達した時、Ｉ ｉ （Ｘｍ、Ｙ
ｎ） ＝ｂ −ｅｘｐ （ａ−Ｃｉ （Ｘｍ、Ｙｎ））
・−（１）となり、Ｃｉ （Ｘｍ、Ｙｎ）がほぼ一定
値Ｃ（Ｘｍ、 Ｙｎ）となるので、この系の入出力特性
は式％式％ となり、対数特性を持つ。
また、この系の積算に関する等優待定数を求めるため、
第２図の積分器を考える。第２図では積算器出力Ｃのに
倍がフィードバックされており、入カニからに−Ｃを減
じたものが積算器入力となり、ループを形成している。
第２図の積分器を考える。第２図では積算器出力Ｃのに
倍がフィードバックされており、入カニからに−Ｃを減
じたものが積算器入力となり、ループを形成している。
この系はよく知られた一次遅れの特性を持ち、その時定
数ｒは１／にで与えられる。入力Ｉｆとして、ある平均
値■のまわりに正規分布をしたものを考えると、出力Ｃ
の相対的標準偏差σは、 σ＝１１５〒７ で与えられるので、 となる。第１図に示した実施例による系では、上記Ｋに
相当する値が２次元積算器出力の値に依存し一定ではな
いが、上記のような同じ画面像が順次入力されている場
合、平衡時には等価的な一次遅れの積算時定数を考える
ことができる。この場合、第２図のＫに相当する値とし
て を与えればよい。
数ｒは１／にで与えられる。入力Ｉｆとして、ある平均
値■のまわりに正規分布をしたものを考えると、出力Ｃ
の相対的標準偏差σは、 σ＝１１５〒７ で与えられるので、 となる。第１図に示した実施例による系では、上記Ｋに
相当する値が２次元積算器出力の値に依存し一定ではな
いが、上記のような同じ画面像が順次入力されている場
合、平衡時には等価的な一次遅れの積算時定数を考える
ことができる。この場合、第２図のＫに相当する値とし
て を与えればよい。
従って、この場合
ｄ ［Ｃ（Ｘｍ、Ｙｎ）］
＝ａ −ｂ −ｅｘｐ （ａ −Ｃ（Ｘｍ、 Ｙｎ）
）＝ａ ・Ｉ （Ｘｍ、Ｙｎ） となり、積分時定数は Ｋ ａ −１（Ｘｍ、 Ｙｎ） −（３）と
なる。第２図の場合と同様に、入力画像信号は平均値１
（Ｘｍ、Ｙｎ）まわりに正規分布をしていると考え
てよいので、積算器１の相対的標準偏差σは となり、入力の信号レベルによらず一定の値となる。つ
まり、画面のどの位置についても出力画像の量子雑音の
相対的標準偏差は一定となる。
）＝ａ ・Ｉ （Ｘｍ、Ｙｎ） となり、積分時定数は Ｋ ａ −１（Ｘｍ、 Ｙｎ） −（３）と
なる。第２図の場合と同様に、入力画像信号は平均値１
（Ｘｍ、Ｙｎ）まわりに正規分布をしていると考え
てよいので、積算器１の相対的標準偏差σは となり、入力の信号レベルによらず一定の値となる。つ
まり、画面のどの位置についても出力画像の量子雑音の
相対的標準偏差は一定となる。
なお、入力画像信号が８ビツトでディジタル化されてい
る場合、画像メモリの深さを１６ビツトにしておくと、
ｂ−１とし、ａ ＝ ｌ ｎ ２５６／６５５３６＝
８．５ｘｌＯとできるので、 つまり、出力画像の標準偏差を０．６５％にすることが
できる。このときの入出力特性を第３図に示した。この
ときの指数関数器２の入出力特性は第３図の入出力を入
れ替えたものにしておけばよい。
る場合、画像メモリの深さを１６ビツトにしておくと、
ｂ−１とし、ａ ＝ ｌ ｎ ２５６／６５５３６＝
８．５ｘｌＯとできるので、 つまり、出力画像の標準偏差を０．６５％にすることが
できる。このときの入出力特性を第３図に示した。この
ときの指数関数器２の入出力特性は第３図の入出力を入
れ替えたものにしておけばよい。
また第７図（Ｃ）に、本実施例による出力画像の１ライ
ン分を従来例と対比して示している。
ン分を従来例と対比して示している。
このように本実施例装置では、重畳フレームの画像レベ
ル変化率が一定となるよう、２次元積算器１の出力を指
数変換し、これと次フレームとの減算を行なって上記２
次元積算器１に入力するようにしたので、画像の明部加
算ではその加算量が一定時刻以降増加せず平衡し、他方
暗部はさらに加算し続は明部と同じ統計精度が得られた
時点で平衡する。従って全画面にわたりノイズコントラ
ストの一定な画像が得られ、また入力画像の濃淡差はレ
ンジ圧縮されてモニター画面上に抽出できる。
ル変化率が一定となるよう、２次元積算器１の出力を指
数変換し、これと次フレームとの減算を行なって上記２
次元積算器１に入力するようにしたので、画像の明部加
算ではその加算量が一定時刻以降増加せず平衡し、他方
暗部はさらに加算し続は明部と同じ統計精度が得られた
時点で平衡する。従って全画面にわたりノイズコントラ
ストの一定な画像が得られ、また入力画像の濃淡差はレ
ンジ圧縮されてモニター画面上に抽出できる。
次に本願の第２の発明の一実施例を図について説明する
。第４図において、１０は（ＭＸＮ）個の要素を持つ２
次元記憶器、２０はこの２次元記憶器ｌＯの出力Ｑｉ
（Ｘｍ、Ｙｎ）を入力とし、これを指数変換［ｂ−ｅｘ
ｐ （ａ−Ｃｉ （Ｘｍ、 Ｙｎ））するとともに
、上記入力値Ｃｉ（Ｘｍ、Ｙｎ）からこの指数変換値を
減算して出力する関数変換器、３０はディジタル化され
た時系列的に送り込まれる入力画像信号Ｉ ｉ （Ｘ
ｍ、Ｙｎ）を正入力とし、入力画像信号の画素（Ｘｍ、
Ｙｎ）に。
。第４図において、１０は（ＭＸＮ）個の要素を持つ２
次元記憶器、２０はこの２次元記憶器ｌＯの出力Ｑｉ
（Ｘｍ、Ｙｎ）を入力とし、これを指数変換［ｂ−ｅｘ
ｐ （ａ−Ｃｉ （Ｘｍ、 Ｙｎ））するとともに
、上記入力値Ｃｉ（Ｘｍ、Ｙｎ）からこの指数変換値を
減算して出力する関数変換器、３０はディジタル化され
た時系列的に送り込まれる入力画像信号Ｉ ｉ （Ｘ
ｍ、Ｙｎ）を正入力とし、入力画像信号の画素（Ｘｍ、
Ｙｎ）に。
対応する２次元記憶器Ｃｉ （Ｘｍ、Ｙｎ）の関数変
換器出力であるＣ ｉ （Ｘｍ、 Ｙｎ） −ｂ −
ｅｘｐ（ａ−Ｃｉ （Ｘｍ、Ｙｎ））を正入力とし、
両者の加算を行なう加算器であり、この加算器出力が２
次元記憶器１０に入力され、その対応する要素に加えら
れてループを形成している。
換器出力であるＣ ｉ （Ｘｍ、 Ｙｎ） −ｂ −
ｅｘｐ（ａ−Ｃｉ （Ｘｍ、Ｙｎ））を正入力とし、
両者の加算を行なう加算器であり、この加算器出力が２
次元記憶器１０に入力され、その対応する要素に加えら
れてループを形成している。
即ち本実施例の構成は、第５図に示すように、上記第１
図に示した実施例の構成と、第６図に示した従来の２次
元積算器の構成とからなるものであり、これが簡略化さ
れて第４図の構成が得られる。従って、関数変換器２０
の変換式ｆは、ｆ＝ｃｉ （Ｘｍ、Ｙｎ） −ｂ −ｅｘｐ （ａ−Ｃｉ （Ｘｍ、 Ｙｎ）
）となり、本実施例全体の作用、効果は上記第１図の実
施例で示したものと全く同様である。
図に示した実施例の構成と、第６図に示した従来の２次
元積算器の構成とからなるものであり、これが簡略化さ
れて第４図の構成が得られる。従って、関数変換器２０
の変換式ｆは、ｆ＝ｃｉ （Ｘｍ、Ｙｎ） −ｂ −ｅｘｐ （ａ−Ｃｉ （Ｘｍ、 Ｙｎ）
）となり、本実施例全体の作用、効果は上記第１図の実
施例で示したものと全く同様である。
また、一画面の画素は（ＭＸＮ）個となっており、これ
ら全画素の一すイクル分の処理がテレビレートの３０
ｍ ｓで実行されるようになっている。
ら全画素の一すイクル分の処理がテレビレートの３０
ｍ ｓで実行されるようになっている。
Ｍ＝Ｎ＝５１２の場合を考えると、一画素につき３０ｍ
５／’（５１２Ｘ５１２）＃１１５ｎｓ以下で１サイク
ル処理を行なう必要があるが、関数変換器２０として予
め所定の関数のテーブルをメモリに書込んでおき、入力
に応じてテーブルを参照し出力するようにしておけばス
タティックＲＡＭのアクセスｄＱｎｓ以内で十分可能で
ある。
５／’（５１２Ｘ５１２）＃１１５ｎｓ以下で１サイク
ル処理を行なう必要があるが、関数変換器２０として予
め所定の関数のテーブルをメモリに書込んでおき、入力
に応じてテーブルを参照し出力するようにしておけばス
タティックＲＡＭのアクセスｄＱｎｓ以内で十分可能で
ある。
なお、上記各実施例ではディジタル化された入−万両像
信号が時系列的にテレビレートで入力されている場合に
ついて説明したが、画像の入力レートはどのような値で
あってもよい。また、１画面の画素数を５１２Ｘ５１２
としたが、この数もどのような値であってもよい。さら
に、上記実施例では（２）式のｂ＝１とした例を示した
が、入力信号レベルがある値以上である場合にはｂとし
てそのレベルに対応する値を選ぶこともできる。
信号が時系列的にテレビレートで入力されている場合に
ついて説明したが、画像の入力レートはどのような値で
あってもよい。また、１画面の画素数を５１２Ｘ５１２
としたが、この数もどのような値であってもよい。さら
に、上記実施例では（２）式のｂ＝１とした例を示した
が、入力信号レベルがある値以上である場合にはｂとし
てそのレベルに対応する値を選ぶこともできる。
以上のように、本願の第１の発明によれば、２次元積算
器の出力を指数関数器を通してフィードバックし、入力
画像信号とこの指数関数器出力を入力とする減算器の出
力を上記２次元積算器の入力とするようループを構成し
、対応する画素毎に演算処理できるようにしたので、ま
た第２の発明では、２次元記憶器の出力を関数変換器を
通してフィードバックし、入力画像信号とこの関数変換
器出力を入力とする加算器の出力を上記２次元記憶器の
入力とするようループを構成し、対応する画素毎に演算
処理できるようにしたので、出力画像における雑音の標
準偏差が信号レベルによらず画面全体にわたり一定とな
り、また出力は入力の対数となるので、モニター画面に
表示される出力画像をダイナミックレンジの狭い人の肉
眼で見たときにも、入力画像の信号レベルによらず広範
囲において小さな信号が視認できるという効果がある。
器の出力を指数関数器を通してフィードバックし、入力
画像信号とこの指数関数器出力を入力とする減算器の出
力を上記２次元積算器の入力とするようループを構成し
、対応する画素毎に演算処理できるようにしたので、ま
た第２の発明では、２次元記憶器の出力を関数変換器を
通してフィードバックし、入力画像信号とこの関数変換
器出力を入力とする加算器の出力を上記２次元記憶器の
入力とするようループを構成し、対応する画素毎に演算
処理できるようにしたので、出力画像における雑音の標
準偏差が信号レベルによらず画面全体にわたり一定とな
り、また出力は入力の対数となるので、モニター画面に
表示される出力画像をダイナミックレンジの狭い人の肉
眼で見たときにも、入力画像の信号レベルによらず広範
囲において小さな信号が視認できるという効果がある。
第１図は本願の第１の発明の一実施例による画像処理装
置を示すブロック図、第２図は一次遅れ積分器を説明す
るためのブロック図、第３図は上記画像処理装置の入出
力特性を示す図、第４図は本願の第２の発明の一実施例
による画像処理装置を示すブロック図、第５図は該装置
を等価回路的に示すブロック図、第６図は従来の画像処
理装置を示すブロック図、第７図（ａ）、 （ｂ）、
（Ｃ１は入力画像。 従来の画像処理装置による出力画像９本発明の一実施例
による画像処理装置による出力画像のそれぞれ１ライン
分を示す図である。 １・・・２次元積算器、２・・・指数変換器、３・・・
減算器、１０・・・２次元記憶器、２０・・・関数変換
器、３０・・・加算器。 なお図中同一符号は同−又は相当部分を示す。
置を示すブロック図、第２図は一次遅れ積分器を説明す
るためのブロック図、第３図は上記画像処理装置の入出
力特性を示す図、第４図は本願の第２の発明の一実施例
による画像処理装置を示すブロック図、第５図は該装置
を等価回路的に示すブロック図、第６図は従来の画像処
理装置を示すブロック図、第７図（ａ）、 （ｂ）、
（Ｃ１は入力画像。 従来の画像処理装置による出力画像９本発明の一実施例
による画像処理装置による出力画像のそれぞれ１ライン
分を示す図である。 １・・・２次元積算器、２・・・指数変換器、３・・・
減算器、１０・・・２次元記憶器、２０・・・関数変換
器、３０・・・加算器。 なお図中同一符号は同−又は相当部分を示す。
Claims (4)
- （１）ディジタル化された画像信号を逐次積算する２次
元積算器と、 該２次元積算器の積算出力の時系列信号を指数変換する
指数関数器と、 該指数関数器からの出力信号を負入力、入力画像信号を
正入力として対応する画素毎に演算を行なう減算器とを
備え、 該減算器の差出力を上記２次元積算器の対応する要素に
入力するようにしたことを特徴とする画像処理装置。 - （２）上記指数関数器は、指数変換テーブルを蓄えたメ
モリを有し、入力値に応じ該テーブルを参照して指数変
換するものであることを特徴とする特許請求の範囲第１
項記載の画像処理装置。 - （３）ディジタル化された画像信号を逐次記憶する２次
元記憶器と、 該２次元記憶器の出力の時系列信号を指数変換するとと
もに該指数変換値を該時系列信号から減算する関数変換
器と、 該関数変換器からの出力信号を正入力、入力画像信号を
正入力として対応する画素毎に演算を行なう加算器とを
備え、 該加算器の和出力を上記２次元記憶器の対応する要素に
入力するようにしたことを特徴とする画像処理装置。 - （４）上記関数変換器は、該変換器への入力値から該入
力値の指数値を減じた値である、各入力値に対応した変
換値を変換テーブル上に蓄えたメモリを有し、入力値に
応じ該テーブルを参照して数値変換するものであること
を特徴とする特許請求の範囲第３項記載の画像処理装置
。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1625585 | 1985-01-29 | ||
| JP60-16255 | 1985-01-29 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS61267875A true JPS61267875A (ja) | 1986-11-27 |
| JPH0650523B2 JPH0650523B2 (ja) | 1994-06-29 |
Family
ID=11911451
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP61014395A Expired - Fee Related JPH0650523B2 (ja) | 1985-01-29 | 1986-01-24 | 画像処理装置 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US4785409A (ja) |
| JP (1) | JPH0650523B2 (ja) |
| GB (1) | GB2170375B (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0486583A (ja) * | 1990-07-30 | 1992-03-19 | Shimadzu Corp | 超音波診断装置 |
| JP2014535039A (ja) * | 2011-09-30 | 2014-12-25 | アナロジック コーポレイション | 光子数の補正方法および装置｛ｐｈｏｔｏｎｃｏｕｎｔｃｏｒｒｅｃｔｉｏｎ｝ |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62274471A (ja) * | 1986-05-23 | 1987-11-28 | Fanuc Ltd | 画像処理装置 |
| JP2833652B2 (ja) * | 1988-07-15 | 1998-12-09 | 株式会社ヤマトヤ商会 | カラーテレビジヨンの画像処理方法及びその装置 |
Family Cites Families (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3872246A (en) * | 1971-10-14 | 1975-03-18 | Siemens Ag | Apparatus for elimination of interference from video signals |
| CH583910A5 (ja) * | 1974-10-25 | 1977-01-14 | Contraves Ag | |
| US3984865A (en) * | 1975-03-26 | 1976-10-05 | Rca Corporation | Transient suppression in television video systems |
| GB1515551A (en) * | 1975-04-25 | 1978-06-28 | British Broadcasting Corp | Noise reduction in electrical signals |
| US4521808A (en) * | 1979-03-22 | 1985-06-04 | University Of Texas System | Electrostatic imaging apparatus |
| US4392123A (en) * | 1980-06-02 | 1983-07-05 | The Dindima Group Pty. Ltd. | Signal-to-noise improving system |
| DE3068972D1 (en) * | 1980-11-28 | 1984-09-20 | Ibm | Raster crt flicker reducing apparatus |
| US4405920A (en) * | 1980-12-31 | 1983-09-20 | Naomi Weisstein | Enhancing the perceptibility of barely perceptible images |
| US4367490A (en) * | 1981-04-24 | 1983-01-04 | General Electric Company | Noise reduction in digital fluoroscopy systems |
| US4458267A (en) | 1981-10-06 | 1984-07-03 | Analogic Corporation | Digital x-ray system |
| US4488245A (en) * | 1982-04-06 | 1984-12-11 | Loge/Interpretation Systems Inc. | Method and means for color detection and modification |
| DE3313430C2 (de) * | 1982-04-15 | 1986-09-25 | Victor Company Of Japan, Ltd., Yokohama, Kanagawa | Rauschunterdrückungsschaltung |
| US4463375A (en) * | 1982-09-07 | 1984-07-31 | The Board Of Trustees Of The Leland Standford Junior University | Multiple-measurement noise-reducing system |
| GB8310239D0 (en) * | 1983-04-15 | 1983-05-18 | Micro Consultants Ltd | Image processing system |
| DE3342353A1 (de) * | 1983-11-23 | 1985-05-30 | Siemens AG, 1000 Berlin und 8000 München | Verfahren zum betrieb eines computertomographen |
| US4638456A (en) * | 1984-01-27 | 1987-01-20 | E. I. Du Pont De Nemours And Company | Spot quantitation |
| US4571635A (en) * | 1984-02-17 | 1986-02-18 | Minnesota Mining And Manufacturing Company | Method of image enhancement by raster scanning |
| US4573035A (en) * | 1984-07-25 | 1986-02-25 | Analogic Corporation | Digital x-ray system |
| US4631598A (en) * | 1984-10-15 | 1986-12-23 | Burkhardt Norman S | High speed, high resolution image processing system |
| US4639784A (en) * | 1985-02-19 | 1987-01-27 | Rca Corporation | Video signal recursive filter adaptively controlled responsive to the history of image motion |
| US4688095A (en) * | 1986-02-07 | 1987-08-18 | Image Technology Incorporated | Programmable image-transformation system |
-
1986
- 1986-01-24 JP JP61014395A patent/JPH0650523B2/ja not_active Expired - Fee Related
- 1986-01-28 US US06/823,278 patent/US4785409A/en not_active Expired - Lifetime
- 1986-01-29 GB GB08602122A patent/GB2170375B/en not_active Expired
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0486583A (ja) * | 1990-07-30 | 1992-03-19 | Shimadzu Corp | 超音波診断装置 |
| JP2014535039A (ja) * | 2011-09-30 | 2014-12-25 | アナロジック コーポレイション | 光子数の補正方法および装置｛ｐｈｏｔｏｎｃｏｕｎｔｃｏｒｒｅｃｔｉｏｎ｝ |
Also Published As
| Publication number | Publication date |
|---|---|
| GB2170375B (en) | 1988-06-22 |
| GB8602122D0 (en) | 1986-03-05 |
| US4785409A (en) | 1988-11-15 |
| JPH0650523B2 (ja) | 1994-06-29 |
| GB2170375A (en) | 1986-07-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0289152B1 (en) | Apparatus for processing an x-ray image | |
| US10580120B2 (en) | Image filtering based on image gradients | |
| EP2552099B1 (en) | Method and camera for providing an estimation of a mean signal to noise ratio value for an image | |
| JPH07325920A (ja) | ディジタルｘ線画像を可視画像として表示する方法及び装置 | |
| US6175657B1 (en) | Adaptive intrafield reducing of Gaussian noise by fuzzy logic processing | |
| US7206460B2 (en) | Method for contrast matching of multiple images of the same object or scene to a common reference image | |
| Ismail et al. | Adapted single scale Retinex algorithm for nighttime image enhancement | |
| US20100092044A1 (en) | Image Processing Apparatus and Method | |
| US20130308841A1 (en) | Method and apparatus for image processing | |
| JPH0650523B2 (ja) | 画像処理装置 | |
| JPH0566751A (ja) | 階調補正回路 | |
| WO2008056305A1 (en) | Noise reduction of an image signal | |
| Watanabe et al. | Improvement of color quality with modified linear multiscale retinex | |
| US11647298B2 (en) | Image processing apparatus, image capturing apparatus, image processing method, and storage medium | |
| US7319544B2 (en) | Processing of digital images | |
| EP3364368B1 (en) | Decoding a bayer-mask or like coded image | |
| JPH0669447B2 (ja) | Ｘ線画像処理装置 | |
| JP2841675B2 (ja) | Ｘ線画像処理装置 | |
| JPS61172541A (ja) | ダブルウインドウ表示装置 | |
| JP2526376B2 (ja) | コントラスト自動変換方式 | |
| JPH09237345A (ja) | クラス分類適応処理における学習方法 | |
| JPH0833922B2 (ja) | 画像ノイズリダクシヨン装置 | |
| Garrett et al. | Real-time digital image processing for radiological imaging | |
| JP3443948B2 (ja) | 放射線撮像装置 | |
| JPH0787398A (ja) | Ｘ線画像診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |