JPH0991423A

JPH0991423A - 画像処理方法および装置

Info

Publication number: JPH0991423A
Application number: JP7250857A
Authority: JP
Inventors: Wataru Ito; 渡伊藤
Original assignee: Fuji Photo Film Co Ltd
Current assignee: Fujifilm Holdings Corp
Priority date: 1995-09-28
Filing date: 1995-09-28
Publication date: 1997-04-04

Abstract

(57)【要約】【課題】画像信号を強調処理するにあたり、高周波ノ
イズやすでにある程度明瞭なコントラストを有する画像
部分についての強調を抑制しつつ、かつ所望の画像部分
については必要な強調処理を行う。【解決手段】第１演算子13から出力された差信号（Ｓ
org −Ｓus）を、差信号（Ｓorg −Ｓus）と強調信号ｆ
（Ｓorg −Ｓus）とが対応付けされた、差信号が放射線
ノイズなどその大きさが極めて小さい範囲における強調
信号は略ゼロに近い小さい値に押さえられたものとな
り、所望の構造物の空間周波数に対応する信号成分に応
じた範囲における強調信号はある程度大きな値を出力
し、所望の構造物よりもコントラストが明瞭な画像部分
に対応する範囲では強調信号は略ゼロに近い小さい値に
押さえられるように設定された変換テーブル12により強
調信号ｆ（Ｓorg −Ｓus）に変換することにより、高周
波ノイズやすでにある程度明瞭なコントラストを有する
画像部分についての強調を抑制しつつ、所望の画像部分
については必要な強調処理を行う。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は画像処理方法および
装置に関するものである。

【０００２】

【従来の技術】従来より、種々の画像取得方法により得
られた画像を表す画像信号に対して、階調処理や周波数
処理等の画像処理を施し、画像の観察読影性能を向上さ
せることが行われている。特に人体を被写体とした放射
線画像のような医用画像の分野においては、医師等の専
門家が、得られた画像に基づいて患者の疾病や傷害の有
無を的確に診断する必要があり、その画像の読影性能を
向上させる画像処理は不可欠なものとなっている。

【０００３】この画像処理のうち、いわゆる周波数強調
処理としては、例えば特開昭53−163571号に示されるよ
うに、原画像の濃度値等の画像信号（オリジナル画像信
号という）Ｓorg を、Ｓproc＝Ｓorg ＋β×（Ｓorg −Ｓus）（４）なる画像信号Ｓprocに変換するものが知られている。こ
こでβはオリジナル画像信号Ｓorg （一例としては濃度
値など）に依存する周波数強調係数、Ｓusは非鮮鋭マス
ク（いわゆるボケマスク）信号である。このボケマスク
信号Ｓusは、２次元に配置された画素に対してオリジナ
ル画像信号Ｓorg を中心画素とするＮ列×Ｎ行（Ｎは奇
数）の画素マトリクスからなるマスク、すなわちボケマ
スクを設定し、Ｓus＝（ΣＳorg ）／Ｎ² （５）等として求められる超低空間周波数成分である。

【０００４】式（４）の第２項括弧内の値（Ｓorg −Ｓ
us）は、オリジナル画像信号から超低空間周波数成分で
あるボケマスク信号を減算したものであるから、オリジ
ナル画像信号のうちの、超低空間周波数成分を除去した
超低空間周波数よりも高い周波数成分Ｓspを意味する。
この比較的高い周波数成分Ｓspに周波数強調係数βを乗
じたうえで、オリジナル画像信号Ｓorg を加算すること
により、オリジナル画像信号Ｓorg のうちこの比較的高
い周波数成分Ｓsp、例えば画像信号の変化が急峻なエッ
ジ部などが相対的に強調された処理済画像信号Ｓprocを
得ることができる。

【０００５】さらに特開昭56−104645号には、下記式
（６）で表される演算をおこなって上記周波数成分Ｓsp
を強調する技術が開示されている。

【０００６】Ｓproc＝Ｓorg ＋Ｆ（ｘ）（６）（ただし、ｘ＝Ｓorg −Ｓus、Ｆ（ｘ）は｜ｘ₁｜＜｜
ｘ₂｜のときＦ′（ｘ₁）≧Ｆ′（ｘ₂）≧０であり、
少なくともｘのある値ｘ₀（｜ｘ₁｜＜｜ｘ₀｜＜｜ｘ
₂｜）を境にしてＦ′（ｘ₁）＞Ｆ′（ｘ₂）となる単
調増加関数）この技術は、Ｆ′（ｘ）すなわちＳorg とＳusの点の関
数の１次微分が｜ｘ｜が大きい程小さくなるようにして
いるので｜ｘ｜が大きいところでの周波数強調の増加率
が抑制され、これによって差信号が大きい程周波数強調
の程度を増大する一方、差信号が大きいところではその
強調の程度の増大を押えて、偽輪郭の発生を防止すると
いうものである。

【０００７】すなわち、放射線画像において差信号の小
さい部分ではボケマスク処理による周波数強調が普通に
行なわれ、差信号の大きい部分（たとえば骨と筋肉の境
界、軟部とガス部の境界、胃のＢa （バリウム）充填部
とその周囲との境界血管造影における血管陰影等）では
周波数強調の程度の増大が押えられて偽輪郭の発生を防
止することができる。

【０００８】

【発明が解決しようとする課題】ところで上述の特開昭
56−104645号に開示された技術によると、Ｆ（ｘ）は単
調増加関数であるため、オリジナル画像信号Ｓorg と非
鮮鋭マスク信号Ｓusとの差信号の増大にしたがってＦ
（ｘ）も増大する。ここで、原画像においてすでにある
程度鮮鋭度の高くコントラストの明瞭な画像部分につい
ては原画像そのままでも観察読影性は十分であるため本
来周波数強調を行う必要はないが、当該画像部分に対応
する差信号は比較的大きい値を示すため、上記式（６）
にしたがって強調処理したのでは当該コントラストの明
瞭な部分も強調され、その強調度合いは、原画像の鮮鋭
度が低いために強調しようとする所望の画像部分の強調
度合いよりも強いものとなり、その結果当該コントラス
トの明瞭な画像部分は過度の強調によりオーバーシュー
トやアンダーシュートを生じて偽輪郭（アーティファク
ト）を生じ、所望の画像部分の観察読影性能を却って低
下させることとなる。

【０００９】一方、オリジナル画像信号Ｓorg と非鮮鋭
マスク信号Ｓusとの差信号である高周波成分は、高周波
領域の放射線ノイズ等の高周波ノイズも含まれており、
この高周波ノイズは強調処理しないようにすることが画
像の観察読影性能を向上させるうえで重要である。

【００１０】本発明は上記事情に鑑みなされたものであ
って、所望の画像部分の強調を行いつつ、高周波の放射
線ノイズやすでにある程度明瞭なコントラストを有する
画像部分の強調を抑制して観察読影性能を向上した処理
済画像を表す処理済画像信号を得る画像処理方法および
装置を提供することを目的とするものである。

【００１１】

【課題を解決するための手段】本発明の画像処理方法
は、画像を表すオリジナル画像信号Ｓorg の非鮮鋭マス
ク信号Ｓusを求め、前記オリジナル画像信号Ｓorg と前
記非鮮鋭マスク信号Ｓusとの差信号を求め、該差信号に
依存する関数ｆ（ｘ）を用いて、下記式（１）にしたが
った強調処理を施すことによって所望の画像部分を周波
数強調する画像処理方法において、Ｓproc＝Ｓorg ＋ｆ（ｘ）（１）（ただし、ｆ（ｘ）は差信号に依存した関数を表す；ｘ
＝Ｓorg −Ｓus）前記所望の画像部分に対応する前記差信号の絶対値が第
１の閾値と該第１の閾値より大きい第２の閾値との間の
範囲にあるときに、前記差信号の絶対値が第１の閾値よ
り小さい範囲における前記関数ｆ（ｘ）の値の絶対値
を、前記所望の画像部分に対応する前記関数ｆ（ｘ）の
値より小さく、かつ前記差信号の絶対値が第２の閾値よ
り大きい範囲における前記関数ｆ（ｘ）の値の絶対値
を、前記所望の画像部分に対応する前記関数ｆ（ｘ）の
値より小さく設定することを特徴とするものである。

【００１２】なお関数ｆは具体的には、図２に示すよう
に、（Ｓorg −Ｓus）≧０においてはｆ≧０であって、
差信号（Ｓorg −Ｓus）が高周波ノイズに対応する第１
の所定値α₁（α₁＞０）までの範囲内においては関数
ｆの値が小さく、差信号（Ｓorg −Ｓus）が所望の画像
部分に対応する第２の所定値α₂（α₂＞α₁）までは
単調に増加し、第２の所定値α₂以上の範囲においては
単調減少し、（Ｓorg−Ｓus）がすでにある程度明瞭な
コントラストを有する画像部分に対応する第３の所定値
α₃（α₃＞α₂）以上の範囲は関数ｆの値が小さくな
るように設定されたものである。

【００１３】また関数ｆを、画像の高周波ノイズを表す
指標値（例えば分散値σ²等）に応じた補正用指標関数
ｇを用いて、このノイズに対応した画像部分の強調を一
層抑制するように、下記式（２）にしたがって補正する
ようにしてもよい。

【００１４】ｆ′（ｘ）＝ｇ（ｘ）×ｆ（ｘ）（２）（ただし、ｇ（ｘ）は０≦ｇ（ｘ）≦１の範囲で単調に
増加する補正用指標関数、ｆ′（ｘ）は式（１）の関数
ｆ（ｘ）に代えて適用される補正後の関数を表す）このように補正することにより式（１）は下記式
（１′）と表すことができる。

【００１５】Ｓproc＝Ｓorg ＋ｇ（ｘ）×ｆ（ｘ）（１′）なお、上記分散値σ²は所定の大きさ（例えばｍ行×ｎ
列）の画素マトリックスからなるマスクを設定し、この
マスク内の全画素のオリジナル画像信号Ｓｏｒｇの平均
値Ｓｍｅａｎを求め、下記演算（式（７）参照）により
求めるものであるが、このマスクを非鮮鋭マスクと同一
のものとすれば、Ｓmean＝Ｓusとなり演算処理を下記式
（８）に示すように簡単化することができる。

【００１６】 σ²＝｛Σ（Ｓorg −Ｓmean）²｝／（ｍ・ｎ）（７） σ²＝（Σ｜Ｓorg −Ｓus｜²）／Ｎ² （８）分散値σ²を画像のノイズを表す指標値とした補正用の
指標関数ｇの具体例を図４に示す。図示の補正用指標関
数ｇは、分散値σ²が高周波ノイズ領域に対応する範囲
内においては関数ｇの出力が小さくなるように、分散値
σ²が信号領域に対応する範囲内においては関数ｇの出
力が１に近くなるように設定されたものである。

【００１７】具体的には、補正用の関数ｇは下記式
（３）で表すことができる。

【００１８】ｇ（ｘ）＝Ｇ（Σ｜Ｓorg −Ｓus｜）（３）本発明の画像処理装置は、上記本発明の画像処理方法を
実施するための具体的な装置であって、画像を表すオリ
ジナル画像信号Ｓorg の非鮮鋭マスク信号Ｓusを求める
非鮮鋭マスク信号演算手段と、前記オリジナル画像信号
Ｓorg と前記非鮮鋭マスク信号Ｓusとの差信号を求める
差信号演算手段と、該差信号を該差信号に依存する関数
ｆ（ｘ）に変換する変換テーブルと、下記式（１）にし
たがった強調処理を施す強調処理手段とを備え、前記画
像のうち所望の画像部分を周波数強調する画像処理装置
において、Ｓproc＝Ｓorg ＋ｆ（ｘ）（１）（ただし、ｆ（ｘ）は差信号に依存した関数を表す；ｘ
＝Ｓorg −Ｓus）前記所望の画像部分に対応する前記差信号の絶対値が第
１の閾値と該第１の閾値より大きい第２の閾値との間の
範囲にあるときに、前記差信号の絶対値が第１の閾値よ
り小さい範囲における前記関数ｆ（ｘ）の値の絶対値
が、前記所望の画像部分に対応する前記関数ｆ（ｘ）の
値より小さく、かつ前記差信号の絶対値が第２の閾値よ
り大きい範囲における前記関数ｆ（ｘ）の値の絶対値
が、前記所望の画像部分に対応する前記関数ｆ（ｘ）の
値より小さく設定されていることを特徴とするものであ
る。

【００１９】また、前記差信号を、前記画像のうち該差
信号に依存する高周波ノイズを表す補正用指標関数ｇ
（ｘ）に変換する補正用変換テーブルと、下記式（２）
にしたがって前記関数ｆ（ｘ）をさらに備えた構成を採
用することもできる。

【００２０】ｆ′（ｘ）＝ｇ（ｘ）×ｆ（ｘ）（２）（ただし、ｇ（ｘ）は０≦ｇ（ｘ）≦１の範囲で単調に
増加する補正用指標関数、ｆ′（ｘ）は式（１）の関数
ｆ（ｘ）に代えて適用される補正後の関数を表す）なお上記強調関数ｆ（ｘ）としては、上記本発明の画像
処理方法における強調関数と同様に、具体的には図２に
示すものなどを適用することができる。

【００２１】さらに、この強調関数ｆ（ｘ）を、オリジ
ナル画像信号のうちノイズを表す分散値σ²等の指標値
に応じた補正用指標関数ｇ（ｘ）により、このノイズに
対応した画像部分の強調が抑制されるように、下記式
（２）にしたがって補正する補正手段をさらに備えた構
成を採ることもできる。

【００２２】ｆ′（ｘ）＝ｇ（ｘ）×ｆ（ｘ）（２）（ただし、ｇ（Ｓorg −Ｓus）は指標関数；０≦ｇ（Ｓ
org −Ｓus）≦１、ｆ′（Ｓorg −Ｓus）は式（１）の
関数ｆ（Ｓorg −Ｓus）に代えて適用される補正後の関
数を表す）具体的には、指標値に応じた関数ｇは下記式（３）で表
すことができる。

【００２３】ｇ（ｘ）＝Ｇ（Σ｜Ｓorg −Ｓus｜²）（３）

【００２４】

【発明の効果】本発明の画像処理方法，装置によれば、
オリジナル画像信号Ｓorg とその非鮮鋭マスク信号との
差信号が小さい範囲においては、その差信号の中に占め
る高周波ノイズの割合が高いため、これの強調処理を抑
制することによって画像の観察読影性能を向上せしめ、
一方、差信号が所望の画像部分に対応する差信号より大
きい範囲における原画像はある程度明瞭なコントラスト
を有するため、これの強調処理を抑制することによって
当該部分に偽輪郭が生じるのを防止し画像の観察読影性
能を向上させることができる。

【００２５】

【発明の実施の形態】以下、本発明の画像処理方法を具
体的に実施するための画像処理装置の実施の形態につい
て図面を用いて説明する。

【００２６】図１は本発明の画像処理装置の第１の実施
の形態を示す概略ブロック図である。図示の画像処理装
置は、画像を表すオリジナル画像信号（濃度値）Ｓorg
の、超低空間周波数に対応する非鮮鋭マスク信号Ｓusを
求めるローパスフィルタ11と、このローパスフィルタ11
により得られた非鮮鋭マスク信号Ｓusとオリジナル画像
信号Ｓorg との差信号（Ｓorg −Ｓus）を算出する第１
演算子13と、得られた差信号（Ｓorg −Ｓus）が入力さ
れて、その差信号（Ｓorg −Ｓus）に対応づけされた所
定の強調信号ｆ（Ｓorg −Ｓus）を出力する変換テーブ
ル12と、強調信号ｆ（Ｓorg −Ｓus）ともとのオリジナ
ル画像信号Ｓorg とを加算処理して処理済画像信号Ｓpr
ocを算出する第２演算子14とを備えた構成である。

【００２７】ローパスフィルタ11は入力されたオリジナ
ル画像信号Ｓorg に対して、例えば３列×３行の画素マ
トリクスからなるボケマスクを設定し、下記式（５）
（Ｎ＝３に設定）にしたがって得られたボケマスク信号
Ｓusを画像強調手段16に出力する。

【００２８】Ｓus＝（ΣＳorg ）／Ｎ² （５）なお、ボケマスクとしては式（５）に示すようにマスク
内の画素値の単純平均を用いるものの他、中心画素から
の距離に応じてマスク内の画素値の重み付けを変化させ
たものを用いることもできる。

【００２９】変換テーブル12は例えば図２に示すような
関数形状にしたがって、差信号（Ｓorg −Ｓus）と強調
信号ｆ（Ｓorg −Ｓus）とが対応付けされている。

【００３０】この関数形状は、（Ｓorg −Ｓus）≧０に
おいてはｆ≧０であって、（Ｓorg−Ｓus）が原画像の
高周波ノイズに対応する第１の所定値α₁（α₁＞０）
までの範囲内においては関数ｆの値が小さく、（Ｓorg
−Ｓus）が所望の画像部分に対応する第２の所定値α₂
（α₂＞α₁）までは単調に増加し、第２の所定値α₂
以上の範囲においては単調減少し、（Ｓorg −Ｓus）が
すでにある程度明瞭なコントラストを有する画像部分に
対応する第３の所定値α₃（α₃＞α₂）以上の範囲は
関数ｆの値が小さくなるように設定されたものである。
なお、（Ｓorg−Ｓus）≦０における関数形状は、（Ｓo
rg −Ｓus）≧０の範囲における関数形状を原点につい
て対称にしたものとなっている。

【００３１】次に本実施形態の画像処理装置の作用につ
いて説明する。

【００３２】まず、所定の画像情報読取装置等により得
られたオリジナル画像信号Ｓorg がローパスフィルタ1
1、第１演算子13、第２演算子14にそれぞれ入力され
る。

【００３３】ローパスフィルタ11は入力されたオリジナ
ル画像信号Ｓorg に対して、例えば３列×３行の画素マ
トリクスからなるボケマスクを設定し、上記式（５）に
したがって得られたボケマスク信号Ｓusを第１演算子13
に出力する。

【００３４】なお、ボケマスクとしては式（５）に示す
ようにマスク内の画素値の単純平均を用いるものの他、
前述したように中心画素からの距離に応じてマスク内の
画素値の重み付けを適宜変化させてもよい。

【００３５】第１演算子13は、入力されたオリジナル画
像信号Ｓorg およびボケマスク信号Ｓusに基づいて、こ
れらの差信号（Ｓorg −Ｓus）を算出し、その結果を変
換テーブル12に出力する。

【００３６】変換テーブル12は入力された差信号（Ｓor
g −Ｓus）をその差信号に応じた強調信号ｆ（Ｓorg −
Ｓus）に変換して第２演算子14に出力する。

【００３７】この強調信号ｆ（Ｓorg −Ｓus）は、前述
したような関数形状にしたがって差信号（Ｓorg −Ｓu
s）と強調信号ｆ（Ｓorg −Ｓus）とが対応付けされて
いるため、差信号が放射線ノイズなどその大きさが極め
て小さい範囲、すなわち第１の所定値の範囲（−α₁＜
Ｓorg −Ｓus＜α₁）における強調信号は略ゼロに近い
小さい値に押さえられたものとなり、この差信号が放射
線ノイズによるものより大きく、所望の構造物の空間周
波数に対応する信号成分に応じた範囲、すなわち差信号
が第１の所定値と第３の所定値の範囲（α₁≦Ｓorg −
Ｓus≦α₃または−α₃≦Ｓorg −Ｓus≦−α₁）にお
ける強調信号はある程度大きな値を出力し、所望の構造
物よりも鮮鋭度が高い画像部分に対応する範囲（α₃＜
Ｓorg −ＳusまたはＳorg −Ｓus＜−α₃）では強調信
号は略ゼロに近い小さい値に押さえられる。

【００３８】第２演算子14は入力されたオリジナル画像
信号Ｓorg と強調信号ｆ（Ｓorg −Ｓus）とを加算処理
し、下記式（１）に示す処理済画像信号Ｓprocを出力す
る。

【００３９】Ｓproc＝Ｓorg ＋ｆ（Ｓorg −Ｓus）（１）この結果、処理済画像信号Ｓprocは、高周波の放射線ノ
イズおよびすでにある程度明瞭なコントラストを有する
画像部分の強調が抑制され、かつ所望の画像部分の鮮鋭
度が強調されたものとなる。

【００４０】したがって本実施形態の画像処理装置によ
れば、原画像においてすでにある程度明瞭なコントラス
トを有する画像部分については強調が抑制されるため、
当該画像部分の過度の強調によるオーバーシュート、ア
ンダーシュートを防止し、アーティファクトの形成を低
減して、観察読影性能の優れた可視画像を再生するのに
適した処理済画像信号を得ることができる。

【００４１】図３は本発明の画像処理装置の第２の実施
の形態を示す概略ブロック図である。図示の画像処理装
置は、図１に示した第１の実施形態の画像処理装置に対
して、強調信号ｆ（Ｓorg −Ｓus）を補正する補正手段
である分散値演算手段15と補正値テーブル16と第３演算
子17とをさらに備えた構成である。

【００４２】すなわち、分散値演算手段15は上述の差信
号（Ｓorg −Ｓus）に基づいて、画像の高周波ノイズ成
分を含む指標値として分散値σ²を求め、補正値テーブ
ル16はこの分散値σ²に応じた補正係数ｇ（０≦ｇ≦
１）を求めるものであり、第３演算子17は得られた補正
係数ｇを強調信号ｆ（Ｓorg −Ｓus）に乗じてその積を
求めるものである。なお分散値σ²は下記式（８）によ
り求められ、これは差信号（Ｓorg −Ｓus）に応じた関
数ｇ（Ｓorg −Ｓus）として表すことができる。

【００４３】 σ²＝（Σ｜Ｓorg −Ｓus｜²）／Ｎ² （８）ここで補正値テーブル16が出力する補正係数ｇは、図４
に示すように、分散値σ²が極小さい範囲においては１
を大幅に下回り、分散値σ²がある程度以上の大きさの
範囲においては１に略等しい値を採るような関数形状と
して分散値σ²に対応付けされている。すなわち分散値
σ²が小さい範囲においては、画像信号よりも高周波ノ
イズを含む割合が大きいため、この範囲においては強調
信号を小さく抑制する。

【００４４】次に本実施形態の画像処理装置の作用につ
いて説明する。

【００４５】まず、所定の画像情報読取装置等により得
られたオリジナル画像信号Ｓorg がローパスフィルタ1
1、第１演算子13、第２演算子14にそれぞれ入力され
る。

【００４６】ローパスフィルタ11は入力されたオリジナ
ル画像信号Ｓorg に対して、例えば３列×３行の画素マ
トリクスからなるボケマスクを設定し、上記式（５）に
したがって得られたボケマスク信号Ｓusを第１演算子13
に出力する。

【００４７】第１演算子13は、入力されたオリジナル画
像信号Ｓorg およびボケマスク信号Ｓusに基づいて、こ
れらの差信号（Ｓorg −Ｓus）を算出し、その結果を変
換テーブル12および分散値演算手段12に出力する。

【００４８】変換テーブル12は入力された差信号（Ｓor
g −Ｓus）をその差信号に応じた強調信号ｆ（Ｓorg −
Ｓus）に変換して第２演算子14に出力する。

【００４９】この強調信号ｆ（Ｓorg −Ｓus）は、差信
号が放射線ノイズなどその大きさが極小さい範囲（−α
₁＜Ｓorg −Ｓus＜α₁）における強調信号は略ゼロに
近い小さい値に押さえられたものとなり、この差信号が
放射線ノイズによるものより大きく、所望の空間周波数
の構造物による信号成分に対応する範囲（α₁≦Ｓorg
−Ｓus≦α₃または−α₃≦Ｓorg −Ｓus≦−α₁）に
おける強調信号はある程度大きな値を出力し、所望の構
造物よりもコントラストが明瞭な画像部分に対応する部
分（α₃＜Ｓorg −ＳusまたはＳorg −Ｓus＜−α₃）
では強調信号は略ゼロに近い小さい値に押さえられ、第
３演算子17に入力される。

【００５０】一方、分散値演算手段15は入力された差信
号（Ｓorg −Ｓus）に応じた分散値σ²を求める。ここ
で、分散値σ²を求めるための局所的に設定されるマス
クは、演算処理の簡単化のため非鮮鋭マスクと同一のも
のとする。これによって分散値マスク内の平均値Ｓmean
とボケマスク信号Ｓusと等しいものとなり、得られた分
散値σ²は補正値テーブル16に入力される。

【００５１】補正値テーブル16は入力された分散値σ²
に対応する補正係数ｇ（Ｓorg −Ｓus）を出力する。す
なわち、分散値σ²が極めて小さい範囲においては１を
大幅に下回る補正係数ｇ（Ｓorg −Ｓus）を、分散値σ
²がある程度以上の大きさの範囲においては１に略等し
い補正係数ｇ（Ｓorg −Ｓus）をそれぞれ出力する。

【００５２】第３演算子17は入力された補正係数ｇ（Ｓ
org −Ｓus）と強調信号ｆ（Ｓorg−Ｓus）との積を求
め、得られた積を第２演算子14に出力する。

【００５３】第２演算子14は入力されたオリジナル画像
信号Ｓorg と補正係数ｇ（Ｓorg −Ｓus）により補正さ
れた強調信号ｆ（Ｓorg −Ｓus）・ｇ（Ｓorg −Ｓus）
とを加算処理し、下記式（１′）に示す処理済画像信号
Ｓprocを出力する。

【００５４】Ｓproc＝Ｓorg ＋ｇ（Ｓorg −Ｓus）×ｆ（Ｓorg −Ｓus）（１′）この結果、処理済画像信号Ｓprocは、高周波の放射線ノ
イズおよびすでにある程度明瞭なコントラストを有する
画像部分の強調が抑制され、かつ所望の画像部分の鮮鋭
度が強調されたものとなる。

【００５５】このように本実施形態の画像処理装置によ
れば、原画像においてすでにある程度明瞭なコントラス
トを有する画像部分については強調が抑制されるため、
過度の強調によるオーバーシュート、アンダーシュート
を防止し、アーティファクトの形成を低減して、観察読
影性能の優れた可視画像を再生するのに適した画像信号
を得ることができる。

【図面の簡単な説明】

【図１】本発明の画像処理装置の第１の実施の形態を示
す概略ブロック図

【図２】差信号（Ｓorg −Ｓus）と強調信号ｆ（Ｓorg
−Ｓus）とを対応付ける変換テーブルを示すグラフ

【図３】本発明の画像処理装置の第２の実施の形態を示
す概略ブロック図

【図４】分散値σ²と補正係数ｇ（Ｓorg −Ｓus）とを
対応付ける変換テーブルを示すグラフ

【符号の説明】

11 ローパスフィルタ 12 変換テーブル 13，14，17 演算子 15 分散値演算手段 16 補正値テーブル

Claims

【特許請求の範囲】

【請求項１】画像を表すオリジナル画像信号Ｓorg の
非鮮鋭マスク信号Ｓusを求め、前記オリジナル画像信号
Ｓorg と前記非鮮鋭マスク信号Ｓusとの差信号を求め、
該差信号に依存する関数ｆ（ｘ）を用いて、下記式
（１）にしたがった強調処理を施すことによって所望の
画像部分を周波数強調する画像処理方法において、Ｓproc＝Ｓorg ＋ｆ（ｘ）（１）（ただし、ｆ（ｘ）は差信号に依存した関数を表す；ｘ
＝Ｓorg −Ｓus）前記所望の画像部分に対応する前記差信号の絶対値が第
１の閾値と該第１の閾値より大きい第２の閾値との間の
範囲にあるときに、前記差信号の絶対値が第１の閾値より小さい範囲におけ
る前記関数ｆ（ｘ）の値の絶対値を、前記所望の画像部
分に対応する前記関数ｆ（ｘ）の値より小さく、かつ前
記差信号の絶対値が第２の閾値より大きい範囲における
前記関数ｆ（ｘ）の値の絶対値を、前記所望の画像部分
に対応する前記関数ｆ（ｘ）の値より小さく設定するこ
とを特徴とする請求項１記載の画像処理方法。
【請求項２】前記差信号に基づく、前記画像のうち高
周波ノイズを表す補正用指標関数ｇ（ｘ）を用いて、下
記式（２）にしたがって前記関数ｆ（ｘ）を補正せしめ
ることを特徴とする請求項１記載の画像処理方法。ｆ′（ｘ）＝ｇ（ｘ）×ｆ（ｘ）（２）（ただし、ｇ（ｘ）は０≦ｇ（ｘ）≦１の範囲で単調に
増加する補正用指標関数、ｆ′（ｘ）は式（１）の関数
ｆ（ｘ）に代えて適用される補正後の関数を表す）
【請求項３】前記補正用指標関数ｇ（ｘ）は、下記式
（３）に示す前記差信号の分散値に応じた関数Ｇ（ｘ）
であることを特徴とする請求項２記載の画像処理方法。ｇ（ｘ）＝Ｇ（Σ｜Ｓorg −Ｓus｜²）（３）
【請求項４】画像を表すオリジナル画像信号Ｓorg の
非鮮鋭マスク信号Ｓusを求める非鮮鋭マスク信号演算手
段と、前記オリジナル画像信号Ｓorg と前記非鮮鋭マス
ク信号Ｓusとの差信号を求める差信号演算手段と、該差
信号を該差信号に依存する関数ｆ（ｘ）に変換する変換
テーブルと、下記式（１）にしたがった強調処理を施す
強調処理手段とを備え、前記画像のうち所望の画像部分
を周波数強調する画像処理装置において、Ｓproc＝Ｓorg ＋ｆ（ｘ）（１）（ただし、ｆ（ｘ）は差信号に依存した関数を表す；ｘ
＝Ｓorg −Ｓus）前記所望の画像部分に対応する前記差信号の絶対値が第
１の閾値と該第１の閾値より大きい第２の閾値との間の
範囲にあるときに、前記差信号の絶対値が第１の閾値より小さい範囲におけ
る前記関数ｆ（ｘ）の値の絶対値が、前記所望の画像部
分に対応する前記関数ｆ（ｘ）の値より小さく、かつ前
記差信号の絶対値が第２の閾値より大きい範囲における
前記関数ｆ（ｘ）の値の絶対値が、前記所望の画像部分
に対応する前記関数ｆ（ｘ）の値より小さく設定されて
いることを特徴とする画像処理装置。
【請求項５】前記差信号を、前記画像のうち該差信号
に依存する高周波ノイズを表す補正用指標関数ｇ（ｘ）
に変換する補正用変換テーブルと、下記式（２）にしたがって前記関数ｆ（ｘ）を補正せし
める補正手段とをさらに備えたことを特徴とする請求項
４記載の画像処理装置。ｆ′（ｘ）＝ｇ（ｘ）×ｆ（ｘ）（２）（ただし、ｇ（ｘ）は０≦ｇ（ｘ）≦１の範囲で単調に
増加する補正用指標関数、ｆ′（ｘ）は式（１）の関数
ｆ（ｘ）に代えて適用される補正後の関数を表す）
【請求項６】前記補正用指標関数ｇ（ｘ）は、下記式
（３）に示す前記差信号の分散値に応じた関数Ｇ（ｘ）
であることを特徴とする請求項５記載の画像処理装置。ｇ（ｘ）＝Ｇ（Σ｜Ｓorg −Ｓus｜²）（３）