JPH0650523B2 - 画像処理装置 - Google Patents
画像処理装置Info
- Publication number
- JPH0650523B2 JPH0650523B2 JP61014395A JP1439586A JPH0650523B2 JP H0650523 B2 JPH0650523 B2 JP H0650523B2 JP 61014395 A JP61014395 A JP 61014395A JP 1439586 A JP1439586 A JP 1439586A JP H0650523 B2 JPH0650523 B2 JP H0650523B2
- Authority
- JP
- Japan
- Prior art keywords
- input
- image
- output
- value
- dimensional
- 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 - Fee Related
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)
- Processing Or Creating Images (AREA)
- Closed-Circuit Television Systems (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Picture Signal Circuits (AREA)
Description
【発明の詳細な説明】 〔産業上の利用分野〕 この発明は、特に工業用，医用Ｘ線テレビに使用され、
ディジタル化された画像信号を処理し、画像の改善を行
なう画像処理装置に関するものである。
ディジタル化された画像信号を処理し、画像の改善を行
なう画像処理装置に関するものである。
工業用，医用Ｘ線テレビにおいて、放射線画像の統計雑
音を低減させるための従来の画像処理装置として第６図
に示すものであった。図において、１０はＭ×Ｎ個の要
素からなる画像メモリ、３０は加算器である。
音を低減させるための従来の画像処理装置として第６図
に示すものであった。図において、１０はＭ×Ｎ個の要
素からなる画像メモリ、３０は加算器である。
このような構成になる画像処理装置において、時刻ｔｉ
におけるディジタル化された入力画像信号Ｉｉ（Ｘｍ，
Ｙｎ）と時刻ｔi-1における画像メモリ（２次元メモ
リ）１０の値Ｃi-1（Ｘｍ，Ｙｎ）が加算器３０で加算
された後、画像メモリ１０に送られ、時刻ｔｉにおける
画像メモリの値Ｃｉ（Ｘｍ，Ｙｎ）となる。（Ｘｍ，Ｙ
ｎ）は画素の位置であり、上記演算は各画素毎に順次行
なわれ、Δｔ＝ｔｉ−ｔi-1の時間の間にＭ×Ｎ個の全
画素の画素毎の加算と画像メモリ１０の更新が行なわれ
る。Δｔは通常テレビレートの３０ｍｓであることが多
い。そして画像メモリ１０の内容は処理後の画像信号と
して出力される。
におけるディジタル化された入力画像信号Ｉｉ（Ｘｍ，
Ｙｎ）と時刻ｔi-1における画像メモリ（２次元メモ
リ）１０の値Ｃi-1（Ｘｍ，Ｙｎ）が加算器３０で加算
された後、画像メモリ１０に送られ、時刻ｔｉにおける
画像メモリの値Ｃｉ（Ｘｍ，Ｙｎ）となる。（Ｘｍ，Ｙ
ｎ）は画素の位置であり、上記演算は各画素毎に順次行
なわれ、Δｔ＝ｔｉ−ｔi-1の時間の間にＭ×Ｎ個の全
画素の画素毎の加算と画像メモリ１０の更新が行なわれ
る。Δｔは通常テレビレートの３０ｍｓであることが多
い。そして画像メモリ１０の内容は処理後の画像信号と
して出力される。
次に動作についてより詳細に説明する。画像信号は、本
来その画像化処理の過程において統計的性質に基く量子
雑音を含むため、信号強度が弱いときには量子雑音の影
響が大きくなり、コントラスト，分解能等を劣化させ
る。そこで静止画像の場合、画像信号をディジタル化
し、各画素毎に積算して量子雑音の影響を低減するよう
な処理がなされる。第６図において、時刻ｔｉにおける
画像メモリ１０の値Ｃｉ（Ｘｍ，Ｙｎ）は、明らかにそ
れまでの画像信号Ｉｊ（Ｘｍ，Ｙｎ）の積算 となっている。従ってこの画像メモリ１０の内容を読み
出し、作成した処理後の画像においては、上記加算回数
をＮとすると、量子雑音による相対統計誤差（以下、統
計誤差と称する）は 倍改善されることになる。
来その画像化処理の過程において統計的性質に基く量子
雑音を含むため、信号強度が弱いときには量子雑音の影
響が大きくなり、コントラスト，分解能等を劣化させ
る。そこで静止画像の場合、画像信号をディジタル化
し、各画素毎に積算して量子雑音の影響を低減するよう
な処理がなされる。第６図において、時刻ｔｉにおける
画像メモリ１０の値Ｃｉ（Ｘｍ，Ｙｎ）は、明らかにそ
れまでの画像信号Ｉｊ（Ｘｍ，Ｙｎ）の積算 となっている。従ってこの画像メモリ１０の内容を読み
出し、作成した処理後の画像においては、上記加算回数
をＮとすると、量子雑音による相対統計誤差（以下、統
計誤差と称する）は 倍改善されることになる。
この画像処理装置の１ライン分の入出力信号を第７図
(a)，(b)に示す。出力信号は各画素毎に入力信号のＮ倍
となっており、統計誤差は 倍改善されている。ところが、この図では一般的なシェ
ーディング特性の場合を示しているが、画像上の信号強
度は場所により異なるので、処理後においても統計誤差
の大きさは画素により異なり、画面全体で一様でない。
このため、図の場合には周辺部の統計誤差は中央部の統
計誤差より大きく、図の凸状信号の上に重畳された信号
Ｓ１及びＳ２は、凸状信号に対する相対強度（コントラ
スト）が同じであっても、中央にあるＳ１信号は判別さ
れ、周辺にあるＳ２信号は判別されないことがある。ま
た十分積算され、周囲の統計的量子雑音が小さくなりＳ
２信号も判別可能な場合にも、この積算画像をモニター
画面上に表示した場合、人は肉眼でＳ１信号とＳ２信号
を同時には視認できないことがある。これは人の視覚系
のダイナミックレンジが狭いためであり、積算された画
像の明るい部分と暗い部分で明度差が大きくなりすぎて
いることによる。
(a)，(b)に示す。出力信号は各画素毎に入力信号のＮ倍
となっており、統計誤差は 倍改善されている。ところが、この図では一般的なシェ
ーディング特性の場合を示しているが、画像上の信号強
度は場所により異なるので、処理後においても統計誤差
の大きさは画素により異なり、画面全体で一様でない。
このため、図の場合には周辺部の統計誤差は中央部の統
計誤差より大きく、図の凸状信号の上に重畳された信号
Ｓ１及びＳ２は、凸状信号に対する相対強度（コントラ
スト）が同じであっても、中央にあるＳ１信号は判別さ
れ、周辺にあるＳ２信号は判別されないことがある。ま
た十分積算され、周囲の統計的量子雑音が小さくなりＳ
２信号も判別可能な場合にも、この積算画像をモニター
画面上に表示した場合、人は肉眼でＳ１信号とＳ２信号
を同時には視認できないことがある。これは人の視覚系
のダイナミックレンジが狭いためであり、積算された画
像の明るい部分と暗い部分で明度差が大きくなりすぎて
いることによる。
従来の画像処理装置は以上のように構成されているの
で、積算された画像の量子雑音の低減効果は画素毎に異
なり、同じコントラストを持つ信号であっても、周辺信
号の強度によって判別できたりできなかったりすること
がある。また積算された画像をモニター画面上に表示し
た場合、画像の明部と暗部の明度差が大きくなっている
ので、ダイナミックレンジが限られている人の視覚にと
っては、量子雑音から判別されている信号も視認されな
いなどの欠点があった。
で、積算された画像の量子雑音の低減効果は画素毎に異
なり、同じコントラストを持つ信号であっても、周辺信
号の強度によって判別できたりできなかったりすること
がある。また積算された画像をモニター画面上に表示し
た場合、画像の明部と暗部の明度差が大きくなっている
ので、ダイナミックレンジが限られている人の視覚にと
っては、量子雑音から判別されている信号も視認されな
いなどの欠点があった。
この発明は上記のような従来のものの欠点を除去するた
めになされたもので、積算処理された画像は全画素にわ
たり同じ相対的大きさの量子雑音を持ち、またこの画像
をモニター画面上に表示した場合、入力画像における明
部の範囲を広くすることができ、小さな信号を肉眼で視
認できる画像処理装置を提供するものである。
めになされたもので、積算処理された画像は全画素にわ
たり同じ相対的大きさの量子雑音を持ち、またこの画像
をモニター画面上に表示した場合、入力画像における明
部の範囲を広くすることができ、小さな信号を肉眼で視
認できる画像処理装置を提供するものである。
本願の第１の発明に係る画像処理装置は、画像信号を画
素毎に逐次積算する２次元積算器と、この２次元積算器
の出力を指数変換する指数関数器と、この指数関数器の
出力と入力画像信号との差を対応する画素毎にとる減算
器とを設け、上記減算器の出力を上記２次元積算器の対
応する要素に入力するようにしたものである。
素毎に逐次積算する２次元積算器と、この２次元積算器
の出力を指数変換する指数関数器と、この指数関数器の
出力と入力画像信号との差を対応する画素毎にとる減算
器とを設け、上記減算器の出力を上記２次元積算器の対
応する要素に入力するようにしたものである。
また本願の第２の発明に係る画像処理装置は、ディジタ
ル画像信号を逐次記憶する２次元記憶器と、この２次元
記憶器の出力からその指数変換値を減算する関数変換器
と、この関数変換器の出力と入力画像信号との和を対応
する画素毎にとる加算器とを設け、上記加算器の出力を
上記２次元記憶器の対応する要素に入力するようにした
ものである。
ル画像信号を逐次記憶する２次元記憶器と、この２次元
記憶器の出力からその指数変換値を減算する関数変換器
と、この関数変換器の出力と入力画像信号との和を対応
する画素毎にとる加算器とを設け、上記加算器の出力を
上記２次元記憶器の対応する要素に入力するようにした
ものである。
第１の発明においては、画像信号の画素毎に２次元積算
器の出力を指数関数器を通してフィードバックし、これ
と対応する次フレーム画素の入力画像信号との差をと
り、これを上記２次元積算器の対応する要素に入力する
ループを作成し、画面のどの位置についても出力画像の
量子雑音の相対的標準偏差を一定にする。
器の出力を指数関数器を通してフィードバックし、これ
と対応する次フレーム画素の入力画像信号との差をと
り、これを上記２次元積算器の対応する要素に入力する
ループを作成し、画面のどの位置についても出力画像の
量子雑音の相対的標準偏差を一定にする。
また第２の発明においては、画像信号の画素毎に２次元
記憶器の出力を関数変換器を通してフィードバックし、
これと対応する次フレーム画素の入力画像信号との和を
とり、これを上記２次元記憶器の対応する要素に入力す
るループを作成し、画面のどの位置についても出力画像
の量子雑音の相対的標準偏差を一定にする。
記憶器の出力を関数変換器を通してフィードバックし、
これと対応する次フレーム画素の入力画像信号との和を
とり、これを上記２次元記憶器の対応する要素に入力す
るループを作成し、画面のどの位置についても出力画像
の量子雑音の相対的標準偏差を一定にする。
以下、本願の第１の発明の一実施例を図について説明す
る。第１図において、１は（Ｍ×Ｎ）個の要素を持つ２
次元積算器、２は２次元積算器１の出力Ｃｉ（Ｘｍ，Ｙ
ｎ）を入力とする指数関数器、３はディジタル化され時
系列的に送り込まれる入力画像信号Ｉｉ（Ｘｍ，Ｙｎ）
を正入力とし、入力画像信号の画素（Ｘｍ，Ｙｎ）に対
応する２次元積算器出力Ｃｉ（Ｘｍ，Ｙｎ）の指数関数
器出力であるｂ・ｅｘｐ｛ａ・Ｃｉ（Ｘｍ，Ｙｎ）｝を
負入力とし、両者の減算を行なう減算器であり、この減
算器出力が２次元積算器１に入力されその対応する要素
に加えられてループを形成している。一画面の要素は
（Ｍ×Ｎ）個でできており、これら全画素の１サイクル
分の処理がテレビレートの３０ｍｓで実行されるように
なっている。Ｍ＝Ｎ＝５１２の場合を考えると、一画素
につき３０ｍｓ／（５１２×５１２）≒１１５ｎｓ以下
で１サイクル処理を行なう必要があるが、指数関数器２
として予め指数関数のテーブルをメモリに書込んでお
き、入力に応じてテーブルを参照して出力するようにし
ておけば、スタティックＲＡＭのアクセスタイム４０ｎ
ｓ以内であり、十分可能である。
る。第１図において、１は（Ｍ×Ｎ）個の要素を持つ２
次元積算器、２は２次元積算器１の出力Ｃｉ（Ｘｍ，Ｙ
ｎ）を入力とする指数関数器、３はディジタル化され時
系列的に送り込まれる入力画像信号Ｉｉ（Ｘｍ，Ｙｎ）
を正入力とし、入力画像信号の画素（Ｘｍ，Ｙｎ）に対
応する２次元積算器出力Ｃｉ（Ｘｍ，Ｙｎ）の指数関数
器出力であるｂ・ｅｘｐ｛ａ・Ｃｉ（Ｘｍ，Ｙｎ）｝を
負入力とし、両者の減算を行なう減算器であり、この減
算器出力が２次元積算器１に入力されその対応する要素
に加えられてループを形成している。一画面の要素は
（Ｍ×Ｎ）個でできており、これら全画素の１サイクル
分の処理がテレビレートの３０ｍｓで実行されるように
なっている。Ｍ＝Ｎ＝５１２の場合を考えると、一画素
につき３０ｍｓ／（５１２×５１２）≒１１５ｎｓ以下
で１サイクル処理を行なう必要があるが、指数関数器２
として予め指数関数のテーブルをメモリに書込んでお
き、入力に応じてテーブルを参照して出力するようにし
ておけば、スタティックＲＡＭのアクセスタイム４０ｎ
ｓ以内であり、十分可能である。
次に動作について説明する。
入力画像として量子雑音の大きい同じ画面の像（例えば
放射線画像）が順次送り込まれているとする。ある数の
画面が入力され平衡に達した時、 Ｉｉ（Ｘｍ，Ｙｎ） ≒ｂ・ｅｘｐ｛ａ・Ｃｉ（Ｘｍ，Ｙｎ）｝…(1) となり、Ｃｉ（Ｘｍ，Ｙｎ）がほぼ一定値Ｃ（Ｘｍ，Ｙ
ｎ）となるので、この系の入出力特性は式(1)より、 となり、対数特性を持つ。
放射線画像）が順次送り込まれているとする。ある数の
画面が入力され平衡に達した時、 Ｉｉ（Ｘｍ，Ｙｎ） ≒ｂ・ｅｘｐ｛ａ・Ｃｉ（Ｘｍ，Ｙｎ）｝…(1) となり、Ｃｉ（Ｘｍ，Ｙｎ）がほぼ一定値Ｃ（Ｘｍ，Ｙ
ｎ）となるので、この系の入出力特性は式(1)より、 となり、対数特性を持つ。
また、この系の積算に関する等価時定数を求めるため、
第２図の積分器を考える。第２図では積算器出力ＣのＫ
倍がフィードバックされており、入力ＩからＫ・Ｃを減
じたものが積算器入力となり、ループを形成している。
この系はよく知られた一次遅れの特性を持ち、その時定
数γは１／Ｋで与えられる。入力Ｉｉとして、ある平均
値Ｉのまわりに正規分布をしたものを考えると、出力Ｃ
の相対的標準偏差σは、 で与えられるので、 となる。第１図に示した実施例による系では、上記Ｋに
相当する値が２次元積算器出力の値に依存し一定ではな
いが、上記のような同じ画面像が順次入力されている場
合、平衡時には等価的な一次遅れの積算時定数を考える
ことができる。この場合、第２図のＫに相当する値とし
て を与えればよい。
第２図の積分器を考える。第２図では積算器出力ＣのＫ
倍がフィードバックされており、入力ＩからＫ・Ｃを減
じたものが積算器入力となり、ループを形成している。
この系はよく知られた一次遅れの特性を持ち、その時定
数γは１／Ｋで与えられる。入力Ｉｉとして、ある平均
値Ｉのまわりに正規分布をしたものを考えると、出力Ｃ
の相対的標準偏差σは、 で与えられるので、 となる。第１図に示した実施例による系では、上記Ｋに
相当する値が２次元積算器出力の値に依存し一定ではな
いが、上記のような同じ画面像が順次入力されている場
合、平衡時には等価的な一次遅れの積算時定数を考える
ことができる。この場合、第２図のＫに相当する値とし
て を与えればよい。
従って、この場合 となり、積分時定数は となる。第２図の場合と同様に、入力画像信号は平均値
Ｉ（Ｘｍ，Ｙｎ）まわりに正規分布をしていると考えて
よいので、積算器１の相対的標準偏差σは となり、入力の信号レベルによらず一定の値となる。つ
まり、画面のどの位置についても出力画像の量子雑音の
相対的標準偏差は一定となる。
Ｉ（Ｘｍ，Ｙｎ）まわりに正規分布をしていると考えて
よいので、積算器１の相対的標準偏差σは となり、入力の信号レベルによらず一定の値となる。つ
まり、画面のどの位置についても出力画像の量子雑音の
相対的標準偏差は一定となる。
なお、入力画像信号が８ビットでディジタル化されてい
る場合、画像メモリの深さを１６ビットにしておくと、
ｂ＝１とし、ａ＝ｌｎ256／65536＝8.5×10-5とできる
ので、 つまり、出力画像の標準偏差を0.65％にすることができ
る。このときの入出力特性を第３図に示した。このとき
の指数関数器２の入出力特性は第３図の入出力を入れ替
えたものにしておけばよい。
る場合、画像メモリの深さを１６ビットにしておくと、
ｂ＝１とし、ａ＝ｌｎ256／65536＝8.5×10-5とできる
ので、 つまり、出力画像の標準偏差を0.65％にすることができ
る。このときの入出力特性を第３図に示した。このとき
の指数関数器２の入出力特性は第３図の入出力を入れ替
えたものにしておけばよい。
また第７図(c)に、本実施例による出力画像の１ライン
分を従来例と対比して示している。
分を従来例と対比して示している。
このように本実施例装置では、重量フレームの画像レベ
ル変化率が一定となるよう、２次元積算器１の出力を指
数変換し、これと次フレームとの減算を行なって上記２
次元積算器１に入力するようにしたので、画像の明部加
算ではその加算量が一定時刻以降増加せず平衡し、他方
暗部はさらに加算し続け明部と同じ統計精度が得られた
時点で平衡する。従って全画面にわたりノイズコントラ
ストの一定な画像が得られ、また入力画像の濃淡差はレ
ンジ圧縮されてモニター画面上に抽出できる。
ル変化率が一定となるよう、２次元積算器１の出力を指
数変換し、これと次フレームとの減算を行なって上記２
次元積算器１に入力するようにしたので、画像の明部加
算ではその加算量が一定時刻以降増加せず平衡し、他方
暗部はさらに加算し続け明部と同じ統計精度が得られた
時点で平衡する。従って全画面にわたりノイズコントラ
ストの一定な画像が得られ、また入力画像の濃淡差はレ
ンジ圧縮されてモニター画面上に抽出できる。
次に本願の第２の発明の一実施例を図について説明す
る。第４図において、１０は（Ｍ×Ｎ）個の要素を持つ
２次元記憶器、２０はこの２次元記憶器１０の出力Ｃｉ
（Ｘｍ，Ｙｎ）を入力とし、これを指数変換［ｂ・ｅｘ
ｐ｛ａ・Ｃｉ（Ｘｍ，Ｙｎ）｝するとともに、上記入力
値Ｃｉ（Ｘｍ，Ｙｎ）からこの指数変換値を減算して出
力する関数変換器、３０はディジタル化された時系列的
に送り込まれる入力画像信号Ｉｉ（Ｘｍ，Ｙｎ）を正入
力とし、入力画像信号の画素（Ｘｍ，Ｙｎ）に対応する
２次元記憶器Ｃｉ（Ｘｍ，Ｙｎ）の関数変換器出力であ
るＣｉ（Ｘｍ，Ｙｎ）−ｂ・ｅｘｐ｛ａ・Ｃｉ（Ｘｍ，
Ｙｎ）｝を正入力とし、両者の加算を行なう加算器であ
り、この加算器出力が２次元記憶器１０に入力され、そ
の対応する要素に加えられてループを形成している。
る。第４図において、１０は（Ｍ×Ｎ）個の要素を持つ
２次元記憶器、２０はこの２次元記憶器１０の出力Ｃｉ
（Ｘｍ，Ｙｎ）を入力とし、これを指数変換［ｂ・ｅｘ
ｐ｛ａ・Ｃｉ（Ｘｍ，Ｙｎ）｝するとともに、上記入力
値Ｃｉ（Ｘｍ，Ｙｎ）からこの指数変換値を減算して出
力する関数変換器、３０はディジタル化された時系列的
に送り込まれる入力画像信号Ｉｉ（Ｘｍ，Ｙｎ）を正入
力とし、入力画像信号の画素（Ｘｍ，Ｙｎ）に対応する
２次元記憶器Ｃｉ（Ｘｍ，Ｙｎ）の関数変換器出力であ
るＣｉ（Ｘｍ，Ｙｎ）−ｂ・ｅｘｐ｛ａ・Ｃｉ（Ｘｍ，
Ｙｎ）｝を正入力とし、両者の加算を行なう加算器であ
り、この加算器出力が２次元記憶器１０に入力され、そ
の対応する要素に加えられてループを形成している。
即ち本実施例の構成は、第５図に示すように、上記第１
図に示した実施例の構成と、第６図に示した従来の２次
元積算器の構成とからなるものであり、これが簡略化さ
れて第４図の構成が得られる。従って、関数変換器２０
の変換式ｆは、 ｆ＝Ｃｉ（Ｘｍ，Ｙｎ） −ｂ・ｅｘｐ｛ａ・Ｃｉ（Ｘｍ，Ｙｎ）｝ となり、本実施例全体の作用，効果は上記第１図の実施
例で示したものと全く同様である。
図に示した実施例の構成と、第６図に示した従来の２次
元積算器の構成とからなるものであり、これが簡略化さ
れて第４図の構成が得られる。従って、関数変換器２０
の変換式ｆは、 ｆ＝Ｃｉ（Ｘｍ，Ｙｎ） −ｂ・ｅｘｐ｛ａ・Ｃｉ（Ｘｍ，Ｙｎ）｝ となり、本実施例全体の作用，効果は上記第１図の実施
例で示したものと全く同様である。
また、一画面の画素は（Ｍ×Ｎ）個となっており、これ
ら全画素の一サイクル分の処理がテレビレートの３０ｍ
ｓで実行されるようになっている。Ｍ＝Ｎ＝５１２の場
合を考えると、一画素につき３０ｍｓ／（５１２×５１
２）≒１１５ｎｓ以下で１サイクル処理を行なう必要が
あるが、関数変換器２０として予め所定の関数のテーブ
ルをメモリに書込んでおき、入力に応じてテーブルを参
照し出力するようにしておけばスタティックＲＡＭのア
クセス４０ｎｓ以内で十分可能である。
ら全画素の一サイクル分の処理がテレビレートの３０ｍ
ｓで実行されるようになっている。Ｍ＝Ｎ＝５１２の場
合を考えると、一画素につき３０ｍｓ／（５１２×５１
２）≒１１５ｎｓ以下で１サイクル処理を行なう必要が
あるが、関数変換器２０として予め所定の関数のテーブ
ルをメモリに書込んでおき、入力に応じてテーブルを参
照し出力するようにしておけばスタティックＲＡＭのア
クセス４０ｎｓ以内で十分可能である。
なお、上記各実施例ではディジタル化された入力画像信
号が時系列的にテレビレートで入力されている場合につ
いて説明したが、画像の入力レートはどのような値であ
ってもよい。また、１画面の画素数を５１２×５１２と
したが、この数もどのような値であってもよい。さら
に、上記実施例では(2)式のｂ＝１とした例を示した
が、入力信号レベルがある値以上である場合にはｂとし
てそのレベルに対応する値を選ぶこともできる。
号が時系列的にテレビレートで入力されている場合につ
いて説明したが、画像の入力レートはどのような値であ
ってもよい。また、１画面の画素数を５１２×５１２と
したが、この数もどのような値であってもよい。さら
に、上記実施例では(2)式のｂ＝１とした例を示した
が、入力信号レベルがある値以上である場合にはｂとし
てそのレベルに対応する値を選ぶこともできる。
以上のように、本願の第１の発明によれば、２次元積算
器の出力を指数関数器を通してフィードバックし、入力
画像信号とこの指数関数器出力を入力とする減算器の出
力を上記２次元積算器の入力とするようループを構成
し、対応する画素毎に演算処理できるようにしたので、
また第２の発明では、２次元記憶器の出力を関数変換器
を通してフィードバックし、入力画像信号とこの関数変
換器出力を入力とする加算器の出力を上記２次元記憶器
の入力とするようループを構成し、対応する画素毎に演
算処理できるようにしたので、出力画像における雑音の
標準偏差が信号レベルによらず画面全体にわたり一定と
なり、また出力は入力の対数となるので、モニター画面
に表示される出力画像をダイナミックレンジの狭い人の
肉眼で見たときにも、入力画像の信号レベルによらず広
範囲において小さな信号が視認できるという効果があ
る。
器の出力を指数関数器を通してフィードバックし、入力
画像信号とこの指数関数器出力を入力とする減算器の出
力を上記２次元積算器の入力とするようループを構成
し、対応する画素毎に演算処理できるようにしたので、
また第２の発明では、２次元記憶器の出力を関数変換器
を通してフィードバックし、入力画像信号とこの関数変
換器出力を入力とする加算器の出力を上記２次元記憶器
の入力とするようループを構成し、対応する画素毎に演
算処理できるようにしたので、出力画像における雑音の
標準偏差が信号レベルによらず画面全体にわたり一定と
なり、また出力は入力の対数となるので、モニター画面
に表示される出力画像をダイナミックレンジの狭い人の
肉眼で見たときにも、入力画像の信号レベルによらず広
範囲において小さな信号が視認できるという効果があ
る。
第１図は本願の第１の発明の一実施例による画像処理装
置を示すブロック図、第２図は一次遅れ積分器を説明す
るためのブロック図、第３図は上記画像処理装置の入出
力特性を示す図、第４図は本願の第２の発明の一実施例
による画像処理装置を示すブロック図、第５図は該装置
を等価回路的に示すブロック図、第６図は従来の画像処
理装置を示すブロック図、第７図(a)，(b)，(c)は入力
画像，従来の画像処理装置による出力画像，本発明の一
実施例による画像処理装置による出力画像のそれぞれ１
ライン分を示す図である。 １…２次元積算器、２…指数変換器、３…減算器、１０
…２次元記憶器、２０…関数変換器、３０…加算器。 なお図中同一符号は同一又は相当部分を示す。
置を示すブロック図、第２図は一次遅れ積分器を説明す
るためのブロック図、第３図は上記画像処理装置の入出
力特性を示す図、第４図は本願の第２の発明の一実施例
による画像処理装置を示すブロック図、第５図は該装置
を等価回路的に示すブロック図、第６図は従来の画像処
理装置を示すブロック図、第７図(a)，(b)，(c)は入力
画像，従来の画像処理装置による出力画像，本発明の一
実施例による画像処理装置による出力画像のそれぞれ１
ライン分を示す図である。 １…２次元積算器、２…指数変換器、３…減算器、１０
…２次元記憶器、２０…関数変換器、３０…加算器。 なお図中同一符号は同一又は相当部分を示す。
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 JPS61267875A (ja) | 1986-11-27 |
| JPH0650523B2 true 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) |
Families Citing this family (4)
| 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 | 株式会社ヤマトヤ商会 | カラーテレビジヨンの画像処理方法及びその装置 |
| JPH07109624B2 (ja) * | 1990-07-30 | 1995-11-22 | 株式会社島津製作所 | 超音波診断装置 |
| WO2013048436A1 (en) * | 2011-09-30 | 2013-04-04 | Analogic Corporation | Photon count correction |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4458267A (en) | 1981-10-06 | 1984-07-03 | Analogic Corporation | Digital x-ray system |
Family Cites Families (20)
| 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 |
| 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
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4458267A (en) | 1981-10-06 | 1984-07-03 | Analogic Corporation | Digital x-ray system |
Also Published As
| Publication number | Publication date |
|---|---|
| GB2170375A (en) | 1986-07-30 |
| US4785409A (en) | 1988-11-15 |
| GB2170375B (en) | 1988-06-22 |
| JPS61267875A (ja) | 1986-11-27 |
| GB8602122D0 (en) | 1986-03-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0289152B1 (en) | Apparatus for processing an x-ray image | |
| US4636850A (en) | Apparatus and method for enhancement of video images | |
| US9424632B2 (en) | System and method for generating high dynamic range images | |
| EP4083903B1 (en) | Method and apparatus for processing image, and storage medium | |
| CN110458827A (zh) | 医学影像坏点的检测方法、装置、设备及介质 | |
| CN108510557B (zh) | 图像色调映射方法及装置 | |
| JPH07325920A (ja) | ディジタルｘ線画像を可視画像として表示する方法及び装置 | |
| KR100206267B1 (ko) | 화상처리장치의 화상데이타 경계성분 강조회로 | |
| CN102223480A (zh) | 图像处理装置及图像处理方法 | |
| JPH0650523B2 (ja) | 画像処理装置 | |
| US7760264B2 (en) | Method of obtaining an image | |
| TWI482147B (zh) | 影像處理方法及影像顯示裝置 | |
| US12568314B2 (en) | Image signal processor, method of operating the image signal processor, and application processor including the image signal processor | |
| JPS60214181A (ja) | 符号化方式 | |
| EP1411472A2 (en) | Gradation conversion of an image | |
| US11394893B2 (en) | Image generation apparatus, image generation method, and storage media | |
| JP4206562B2 (ja) | 動体表示用画像処理装置及び方法 | |
| CN110120015A (zh) | 降低图像的噪声的方法和设备 | |
| KR20200094940A (ko) | 히트맵을 생성하는 장치 및 방법 | |
| Garrett et al. | Real-time digital image processing for radiological imaging | |
| JP2696997B2 (ja) | デジタル画像データ変換装置 | |
| JPH09237345A (ja) | クラス分類適応処理における学習方法 | |
| JPH03198836A (ja) | Ｘ線画像処理装置 | |
| JPH0833922B2 (ja) | 画像ノイズリダクシヨン装置 | |
| JP2837969B2 (ja) | 画像処理方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |