JPH0991420A - 画像処理方法および装置 - Google Patents
画像処理方法および装置Info
- Publication number
- JPH0991420A JPH0991420A JP7250854A JP25085495A JPH0991420A JP H0991420 A JPH0991420 A JP H0991420A JP 7250854 A JP7250854 A JP 7250854A JP 25085495 A JP25085495 A JP 25085495A JP H0991420 A JPH0991420 A JP H0991420A
- Authority
- JP
- Japan
- Prior art keywords
- ssub
- mask
- sorg
- signal
- original image
- 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.)
- Withdrawn
Links
- 238000003672 processing method Methods 0.000 title claims description 22
- 238000006243 chemical reaction Methods 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 description 13
- 230000005855 radiation Effects 0.000 description 9
- 208000004434 Calcinosis Diseases 0.000 description 7
- 101000920618 Homo sapiens Transcription and mRNA export factor ENY2 Proteins 0.000 description 5
- 102100031954 Transcription and mRNA export factor ENY2 Human genes 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000002308 calcification Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 206010006187 Breast cancer Diseases 0.000 description 1
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- 241001503485 Mammuthus Species 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Landscapes
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
(57)【要約】
【課題】 画像のアンダーシュート、オーバーシュート
を抑制し、ノイズ成分等の画像読影に不要な成分を強調
することなく、注目する特定の画像部分だけを効率よく
強調処理する。 【解決手段】 ローパスフィルタ11および演算素子14a
によりオリジナル画像信号Ｓorg のうち高周波数成分Ｓ
spを抽出し、一方、マスクサイズの１／２未満の大きさ
以下の画像部分を差信号Ｓsub としてメディアンフィル
タ12および演算子14b により抽出し、構造物については
ある程度大きい値を示すが高周波ノイズについては極小
さな値を示す、記憶手段13に記憶された変換テーブルに
より変換される強調係数β（Ｓsub ）に変換し、演算素
子14c ，14d により下記式（１）にしたがった強調処理
を行うことにより、高周波ノイズを強調することなく、
注目する特定の画像部分だけを強調処理する。 Ｓproc＝Ｓorg ＋β（Ｓsub ）×Ｓsp
（１）
を抑制し、ノイズ成分等の画像読影に不要な成分を強調
することなく、注目する特定の画像部分だけを効率よく
強調処理する。 【解決手段】 ローパスフィルタ11および演算素子14a
によりオリジナル画像信号Ｓorg のうち高周波数成分Ｓ
spを抽出し、一方、マスクサイズの１／２未満の大きさ
以下の画像部分を差信号Ｓsub としてメディアンフィル
タ12および演算子14b により抽出し、構造物については
ある程度大きい値を示すが高周波ノイズについては極小
さな値を示す、記憶手段13に記憶された変換テーブルに
より変換される強調係数β（Ｓsub ）に変換し、演算素
子14c ，14d により下記式（１）にしたがった強調処理
を行うことにより、高周波ノイズを強調することなく、
注目する特定の画像部分だけを強調処理する。 Ｓproc＝Ｓorg ＋β（Ｓsub ）×Ｓsp
（１）
Description
【０００１】
【発明の属する技術分野】本発明は画像処理方法および
装置に関し、詳細には画像中の構造物を大きさに応じて
抽出することにより、その構造物を選択的に強調処理す
る画像処理方法および装置に関するものである。
装置に関し、詳細には画像中の構造物を大きさに応じて
抽出することにより、その構造物を選択的に強調処理す
る画像処理方法および装置に関するものである。
【０００２】
【従来の技術】従来より、種々の画像取得方法により得
られた画像を表す画像信号に対して、階調処理や周波数
処理等の画像処理を施し、画像の観察読影性能を向上さ
せることが行われている。特に人体を被写体とした放射
線画像のような医用画像の分野においては、医師等の専
門家が、得られた画像に基づいて患者の疾病や傷害の有
無を的確に診断する必要があり、その画像の読影性能を
向上させる画像処理は不可欠なものとなっている。
られた画像を表す画像信号に対して、階調処理や周波数
処理等の画像処理を施し、画像の観察読影性能を向上さ
せることが行われている。特に人体を被写体とした放射
線画像のような医用画像の分野においては、医師等の専
門家が、得られた画像に基づいて患者の疾病や傷害の有
無を的確に診断する必要があり、その画像の読影性能を
向上させる画像処理は不可欠なものとなっている。
【０００３】この画像処理のうち、いわゆる周波数強調
処理としては、例えば特開昭61-169971 号に示されるよ
うに、原画像の濃度値等の画像信号（オリジナル画像信
号という）Ｓorg を、 Ｓproc＝Ｓorg ＋β×（Ｓorg −Ｓus） （４） なる画像信号Ｓprocに変換するものが知られている。こ
こでβは周波数強調係数、Ｓusは非鮮鋭マスク（いわゆ
るボケマスク）信号である。このボケマスク信号Ｓus
は、２次元に配置された画素に対してオリジナル画像信
号Ｓorg を中心画素とするＮ列×Ｎ行（Ｎは奇数）の画
素マトリクスからなるマスク、すなわちボケマスクを設
定し、 Ｓus＝（ΣＳorg ）／Ｎ2 （５） 等として求められる超低空間周波数成分である。
処理としては、例えば特開昭61-169971 号に示されるよ
うに、原画像の濃度値等の画像信号（オリジナル画像信
号という）Ｓorg を、 Ｓproc＝Ｓorg ＋β×（Ｓorg −Ｓus） （４） なる画像信号Ｓprocに変換するものが知られている。こ
こでβは周波数強調係数、Ｓusは非鮮鋭マスク（いわゆ
るボケマスク）信号である。このボケマスク信号Ｓus
は、２次元に配置された画素に対してオリジナル画像信
号Ｓorg を中心画素とするＮ列×Ｎ行（Ｎは奇数）の画
素マトリクスからなるマスク、すなわちボケマスクを設
定し、 Ｓus＝（ΣＳorg ）／Ｎ2 （５） 等として求められる超低空間周波数成分である。
【０００４】式（４）の第２項括弧内の値（Ｓorg −Ｓ
us）は、オリジナル画像信号から超低空間周波数成分で
あるボケマスク信号を減算したものであるから、オリジ
ナル画像信号のうちの、超低空間周波数成分を除去した
超低空間周波数よりも高い周波数成分Ｓspを意味する。
この比較的高い周波数成分Ｓspに周波数強調係数β（＝
一定）を乗じたうえで、オリジナル画像信号Ｓorg を加
算することにより、オリジナル画像信号Ｓorg のうちこ
の比較的高い周波数成分Ｓspが相対的に強調された処理
済画像信号Ｓprocを得ることができる。
us）は、オリジナル画像信号から超低空間周波数成分で
あるボケマスク信号を減算したものであるから、オリジ
ナル画像信号のうちの、超低空間周波数成分を除去した
超低空間周波数よりも高い周波数成分Ｓspを意味する。
この比較的高い周波数成分Ｓspに周波数強調係数β（＝
一定）を乗じたうえで、オリジナル画像信号Ｓorg を加
算することにより、オリジナル画像信号Ｓorg のうちこ
の比較的高い周波数成分Ｓspが相対的に強調された処理
済画像信号Ｓprocを得ることができる。
【０００５】しかし上記強調処理には、以下のような欠
点があることが見出された。すなわち、上記周波数強調
によると、対象物の大きさに関係なくそのエッジが強調
されてしまうので、大きな構造物と小さな構造物が重な
った部分では、かえって画像が見にくくなるという難点
がある。また、上記重なった部分以外においても、対象
物のエッジ部で余分なオーバーシュートあるいはアンダ
ーシュートが発生して、偽輪郭（アーティファクト）を
発生せしめて画像読影性能上好ましくない結果をもたら
す。
点があることが見出された。すなわち、上記周波数強調
によると、対象物の大きさに関係なくそのエッジが強調
されてしまうので、大きな構造物と小さな構造物が重な
った部分では、かえって画像が見にくくなるという難点
がある。また、上記重なった部分以外においても、対象
物のエッジ部で余分なオーバーシュートあるいはアンダ
ーシュートが発生して、偽輪郭（アーティファクト）を
発生せしめて画像読影性能上好ましくない結果をもたら
す。
【０００６】そこでこの問題を解決するものとして、上
記式（４）において、ボケマスク信号Ｓusに代えて所望
のサイズのマスク内の画像信号Ｓorg のメディアン値Ｓ
medを用いる技術が特開昭61−119250号に開示されてい
る。
記式（４）において、ボケマスク信号Ｓusに代えて所望
のサイズのマスク内の画像信号Ｓorg のメディアン値Ｓ
medを用いる技術が特開昭61−119250号に開示されてい
る。
【０００７】この技術によれば、メディアン値Ｓmed を
求めるマスクサイズの１／２未満の大きさの構造物に対
応する画像信号だけを選択的に抽出することができ、こ
の選択的に抽出された画像信号（Ｓorg −Ｓmed ）（＝
Ｓsub ）に対して下記式（６）に示すように強調係数β
を乗じることによって、観察対象たる構造物を有効に強
調処理することができる。
求めるマスクサイズの１／２未満の大きさの構造物に対
応する画像信号だけを選択的に抽出することができ、こ
の選択的に抽出された画像信号（Ｓorg −Ｓmed ）（＝
Ｓsub ）に対して下記式（６）に示すように強調係数β
を乗じることによって、観察対象たる構造物を有効に強
調処理することができる。
【０００８】 Ｓproc＝Ｓorg ＋β×Ｓsub （６）
【０００９】
【発明が解決しようとする課題】しかしながら、上記技
術は信号成分と放射線の量子ノイズ等各種の高周波ノイ
ズ成分を分離するものではないため、信号成分を強調処
理すると同時にノイズ成分も強調処理され、画像の読影
性能を向上させるためにノイズ成分の有効な分離が望ま
れている。
術は信号成分と放射線の量子ノイズ等各種の高周波ノイ
ズ成分を分離するものではないため、信号成分を強調処
理すると同時にノイズ成分も強調処理され、画像の読影
性能を向上させるためにノイズ成分の有効な分離が望ま
れている。
【００１０】すなわち上記技術によれば、選択的に抽出
された画像信号Ｓsub はメディアン値Ｓmed を求めるマ
スクサイズの１／２以上の構造物を有効に排除すること
ができるものの、１／２未満の大きさの構造物の中に
は、通常は極めて小さいサイズの高周波ノイズ等も含ま
れるため、このノイズ等を分離してより読影性能の高い
画像を再生するための処理済画像信号Ｓprocを得ること
が望まれている。
された画像信号Ｓsub はメディアン値Ｓmed を求めるマ
スクサイズの１／２以上の構造物を有効に排除すること
ができるものの、１／２未満の大きさの構造物の中に
は、通常は極めて小さいサイズの高周波ノイズ等も含ま
れるため、このノイズ等を分離してより読影性能の高い
画像を再生するための処理済画像信号Ｓprocを得ること
が望まれている。
【００１１】本発明は上記事情に鑑みなされたものであ
って、ノイズ成分等の画像読影に不要な成分を強調する
ことなく、注目する特定の画像部分だけを効率よく強調
処理し、アーティファクトの発生を抑制した画像処理方
法および装置を提供することを目的とするものである。
って、ノイズ成分等の画像読影に不要な成分を強調する
ことなく、注目する特定の画像部分だけを効率よく強調
処理し、アーティファクトの発生を抑制した画像処理方
法および装置を提供することを目的とするものである。
【００１２】
【課題を解決するための手段】本発明の第１の画像処理
方法は、非鮮鋭マスク信号Ｓusとオリジナル画像信号Ｓ
org との差信号Ｓspに対して、メディアン値Ｓmed とオ
リジナル画像信号Ｓorg との差信号Ｓsub に応じて変化
する強調係数β（Ｓsub ）を用いた下記式（１）にした
がった演算 Ｓproc＝Ｓorg ＋β（Ｓsub ）×Ｓsp （１） を施すことによって、オリジナル画像信号Ｓorg のうち
の比較的高い周波数成分Ｓspの中に含まれる放射線の量
子ノイズ等各種の高周波ノイズ成分に対する強調を抑制
するものである。
方法は、非鮮鋭マスク信号Ｓusとオリジナル画像信号Ｓ
org との差信号Ｓspに対して、メディアン値Ｓmed とオ
リジナル画像信号Ｓorg との差信号Ｓsub に応じて変化
する強調係数β（Ｓsub ）を用いた下記式（１）にした
がった演算 Ｓproc＝Ｓorg ＋β（Ｓsub ）×Ｓsp （１） を施すことによって、オリジナル画像信号Ｓorg のうち
の比較的高い周波数成分Ｓspの中に含まれる放射線の量
子ノイズ等各種の高周波ノイズ成分に対する強調を抑制
するものである。
【００１３】すなわち本願発明の画像処理方法は、第１
のマスクを用いて、画像を表すオリジナル画像信号Ｓor
g の非鮮鋭マスク信号Ｓusを求め、第２のマスクを用い
て、第２のマスク内に含まれるオリジナル画像信号Ｓor
gのメディアン値Ｓmed を求め、非鮮鋭マスク信号Ｓus
とオリジナル画像信号Ｓorg との差信号Ｓspに対して、
メディアン値Ｓmed とオリジナル画像信号Ｓorg との差
信号Ｓsub に応じて変化する強調係数β（Ｓsub ）を用
いた上記式（１）にしたがった演算を施すことを特徴と
するものである。
のマスクを用いて、画像を表すオリジナル画像信号Ｓor
g の非鮮鋭マスク信号Ｓusを求め、第２のマスクを用い
て、第２のマスク内に含まれるオリジナル画像信号Ｓor
gのメディアン値Ｓmed を求め、非鮮鋭マスク信号Ｓus
とオリジナル画像信号Ｓorg との差信号Ｓspに対して、
メディアン値Ｓmed とオリジナル画像信号Ｓorg との差
信号Ｓsub に応じて変化する強調係数β（Ｓsub ）を用
いた上記式（１）にしたがった演算を施すことを特徴と
するものである。
【００１４】ここで、第２のマスクのサイズは第１のマ
スクのサイズよりも大きく、特に第１のマスクのサイズ
の２倍以上に設定するのが望ましい。例えばマンモ画像
の石灰化陰影について強調処理を施す場合には、第１の
マスク（ボケマスク）の画素マトリクスのサイズは３列
×３行、第２のマスク（メディアン値を求めるマスク）
の画素マトリクスのサイズは９列×９行（＝81画素、し
たがって40画素以下の石灰化陰影が抽出される）程度が
望ましい。以下の発明においても同様である。
スクのサイズよりも大きく、特に第１のマスクのサイズ
の２倍以上に設定するのが望ましい。例えばマンモ画像
の石灰化陰影について強調処理を施す場合には、第１の
マスク（ボケマスク）の画素マトリクスのサイズは３列
×３行、第２のマスク（メディアン値を求めるマスク）
の画素マトリクスのサイズは９列×９行（＝81画素、し
たがって40画素以下の石灰化陰影が抽出される）程度が
望ましい。以下の発明においても同様である。
【００１５】また、強調係数β（Ｓsub ）は、例えば図
２（１）〜（３）に示すように、差信号Ｓsub の絶対値
｜Ｓsub ｜が所定の範囲内（０≦｜Ｓsub ｜≦Ｓ0 ）に
おいては、絶対値｜Ｓsub ｜が増大するにしたがって単
調に増大する関数として設定されたものが望ましい。
２（１）〜（３）に示すように、差信号Ｓsub の絶対値
｜Ｓsub ｜が所定の範囲内（０≦｜Ｓsub ｜≦Ｓ0 ）に
おいては、絶対値｜Ｓsub ｜が増大するにしたがって単
調に増大する関数として設定されたものが望ましい。
【００１６】この所定の範囲を規定する閾値Ｓ0 は、高
周波ノイズを表す信号値についての強調係数β（Ｓsub
）と強調しようとする所望の構造物を表す信号値につ
いての強調係数β（Ｓsub ）とを分離して認識しうるの
に十分な差信号Ｓsub の値であって、予め実験等により
求められる値である。
周波ノイズを表す信号値についての強調係数β（Ｓsub
）と強調しようとする所望の構造物を表す信号値につ
いての強調係数β（Ｓsub ）とを分離して認識しうるの
に十分な差信号Ｓsub の値であって、予め実験等により
求められる値である。
【００１７】本発明の第２の画像処理方法は、上記第１
の画像処理方法の第１のマスクと第２のマスクとをそれ
ぞれ互いにサイズの異なるものとして複数ずつ設定し
て、大きさの異なる複数種類の構造物に対して各別の強
調度合いの強調処理を施すことにより、各構造物の強調
度合いを画像読影の目的に応じて適切化するものであ
る。
の画像処理方法の第１のマスクと第２のマスクとをそれ
ぞれ互いにサイズの異なるものとして複数ずつ設定し
て、大きさの異なる複数種類の構造物に対して各別の強
調度合いの強調処理を施すことにより、各構造物の強調
度合いを画像読影の目的に応じて適切化するものであ
る。
【００１８】すなわち本発明の第２の画像処理方法は、
互いに異なる周波数帯域にそれぞれ対応した互いに大き
さの異なるメディアン値算出用の複数のマスクを用い
て、画像を表すオリジナル画像信号Ｓorg についての複
数のメディアン値Ｓmed n （ｎ＝１，２，…）を求め、
前記オリジナル画像信号Ｓorg を、前記複数のマスクの
大きさにそれぞれ対応した互いに異なる複数の周波数成
分Ｓn に分割し、前記複数の周波数成分Ｓn に対して、
該周波数成分に対応する前記複数のメディアン値Ｓmed
n と前記オリジナル画像信号Ｓorg との差信号Ｓsub n
に応じて変化する複数の強調係数βn （Ｓsub n ）を用
いた下記式（２）または（３）にしたがった演算を施す
ことを特徴とするものである。
互いに異なる周波数帯域にそれぞれ対応した互いに大き
さの異なるメディアン値算出用の複数のマスクを用い
て、画像を表すオリジナル画像信号Ｓorg についての複
数のメディアン値Ｓmed n （ｎ＝１，２，…）を求め、
前記オリジナル画像信号Ｓorg を、前記複数のマスクの
大きさにそれぞれ対応した互いに異なる複数の周波数成
分Ｓn に分割し、前記複数の周波数成分Ｓn に対して、
該周波数成分に対応する前記複数のメディアン値Ｓmed
n と前記オリジナル画像信号Ｓorg との差信号Ｓsub n
に応じて変化する複数の強調係数βn （Ｓsub n ）を用
いた下記式（２）または（３）にしたがった演算を施す
ことを特徴とするものである。
【００１９】 Ｓproc＝Σ｛Ｓn ×βn （Ｓsub n ）｝ （２） Ｓproc＝Ｓorg ＋Σ｛Ｓn ×βn （Ｓsub n ）｝ （３） 本発明の第１の画像処理装置は、上記本発明の第１の画
像処理方法を実施するための具体的な装置であって、第
１のマスクを用いて、画像を表すオリジナル画像信号Ｓ
org の非鮮鋭マスク信号Ｓusを求める非鮮鋭マスク信号
演算手段と、第２のマスクを用いて、第２のマスク内に
含まれるオリジナル画像信号Ｓorgのメディアン値Ｓmed
を求めるメディアン値演算手段と、メディアン値Ｓmed
とオリジナル画像信号Ｓorg との差信号Ｓsub に対応
付けされた強調係数β（Ｓsub ）に変換する変換テーブ
ルを記憶してなる記憶手段と、非鮮鋭マスク信号Ｓusと
オリジナル画像信号Ｓorg との差信号Ｓspに対して、強
調係数β（Ｓsub ）を用いた上記式（１）にしたがった
強調処理を施す強調手段とを備えてなることを特徴とす
るものである。
像処理方法を実施するための具体的な装置であって、第
１のマスクを用いて、画像を表すオリジナル画像信号Ｓ
org の非鮮鋭マスク信号Ｓusを求める非鮮鋭マスク信号
演算手段と、第２のマスクを用いて、第２のマスク内に
含まれるオリジナル画像信号Ｓorgのメディアン値Ｓmed
を求めるメディアン値演算手段と、メディアン値Ｓmed
とオリジナル画像信号Ｓorg との差信号Ｓsub に対応
付けされた強調係数β（Ｓsub ）に変換する変換テーブ
ルを記憶してなる記憶手段と、非鮮鋭マスク信号Ｓusと
オリジナル画像信号Ｓorg との差信号Ｓspに対して、強
調係数β（Ｓsub ）を用いた上記式（１）にしたがった
強調処理を施す強調手段とを備えてなることを特徴とす
るものである。
【００２０】また、記憶手段に記憶された変換テーブル
としては、強調係数β（Ｓsub ）が、例えば図２（１）
〜（３）に示すように、差信号Ｓsub の絶対値｜Ｓsub
｜が所定の範囲内（０≦｜Ｓsub ｜≦Ｓ0 ）において
は、絶対値｜Ｓsub ｜が増大するにしたがって単調に増
大するものとして設定されたものが望ましい。
としては、強調係数β（Ｓsub ）が、例えば図２（１）
〜（３）に示すように、差信号Ｓsub の絶対値｜Ｓsub
｜が所定の範囲内（０≦｜Ｓsub ｜≦Ｓ0 ）において
は、絶対値｜Ｓsub ｜が増大するにしたがって単調に増
大するものとして設定されたものが望ましい。
【００２１】本発明の第２の画像処理装置は、上記第１
の画像処理装置の第１のマスクと第２のマスクとをそれ
ぞれ互いにサイズの異なるものとして複数ずつ設定し
て、大きさの異なる複数種類の構造物に対して各別の強
調処理を施すものである。
の画像処理装置の第１のマスクと第２のマスクとをそれ
ぞれ互いにサイズの異なるものとして複数ずつ設定し
て、大きさの異なる複数種類の構造物に対して各別の強
調処理を施すものである。
【００２２】すなわち本発明の第２の画像処理装置は、
互いに異なる周波数帯域にそれぞれ対応した互いに大き
さの異なるメディアン値算出用の複数のマスクを用い
て、画像を表すオリジナル画像信号Ｓorg についての複
数のメディアン値Ｓmed n （ｎ＝１，２，…）を求める
メディアン値演算手段と、前記オリジナル画像信号Ｓor
g を、前記複数のマスクの大きさにそれぞれ対応した互
いに異なる複数の周波数成分Ｓn に分割する周波数帯域
分割手段と、前記複数のメディアン値Ｓmed n と前記オ
リジナル画像信号Ｓorg との差信号Ｓsub n を該差信号
に対応付された強調係数βn （Ｓsub n ）に変換する１
以上の変換テーブルを記憶してなる記憶手段と、前記複
数の周波数成分Ｓn に対して、該周波数成分に対応する
前記強調係数βn （Ｓsub n ）を用いた上記式（２）ま
たは（３）にしたがった強調処理を施す強調処理手段と
を備えてなることを特徴とするものである。
互いに異なる周波数帯域にそれぞれ対応した互いに大き
さの異なるメディアン値算出用の複数のマスクを用い
て、画像を表すオリジナル画像信号Ｓorg についての複
数のメディアン値Ｓmed n （ｎ＝１，２，…）を求める
メディアン値演算手段と、前記オリジナル画像信号Ｓor
g を、前記複数のマスクの大きさにそれぞれ対応した互
いに異なる複数の周波数成分Ｓn に分割する周波数帯域
分割手段と、前記複数のメディアン値Ｓmed n と前記オ
リジナル画像信号Ｓorg との差信号Ｓsub n を該差信号
に対応付された強調係数βn （Ｓsub n ）に変換する１
以上の変換テーブルを記憶してなる記憶手段と、前記複
数の周波数成分Ｓn に対して、該周波数成分に対応する
前記強調係数βn （Ｓsub n ）を用いた上記式（２）ま
たは（３）にしたがった強調処理を施す強調処理手段と
を備えてなることを特徴とするものである。
【００２３】
【発明の効果】本発明の画像処理方法，装置は、画像の
うち強調処理しようとする構造物の大きさに対応する空
間周波数に応じたマスク（第２のマスク）内に含まれる
オリジナル画像信号Ｓorg のメディアン値Ｓmed を求
め、画像を構成する画素ごとに、オリジナル画像信号Ｓ
org と求められたメディアン値Ｓmed との差信号Ｓsub
を求める。この差信号Ｓsub は、上記空間周波数に対応
する大きさの１／２以上の大きさの構造物については０
（ゼロ）を示し、１／２未満の大きさの構造物について
は０以外の値を採るため、所望とする構造物の大きさの
２倍を超える大きさのマスクを予め設定すればよい。
うち強調処理しようとする構造物の大きさに対応する空
間周波数に応じたマスク（第２のマスク）内に含まれる
オリジナル画像信号Ｓorg のメディアン値Ｓmed を求
め、画像を構成する画素ごとに、オリジナル画像信号Ｓ
org と求められたメディアン値Ｓmed との差信号Ｓsub
を求める。この差信号Ｓsub は、上記空間周波数に対応
する大きさの１／２以上の大きさの構造物については０
（ゼロ）を示し、１／２未満の大きさの構造物について
は０以外の値を採るため、所望とする構造物の大きさの
２倍を超える大きさのマスクを予め設定すればよい。
【００２４】なお、このマスクサイズの１／２未満の大
きさの構造物の中には高周波ノイズも含まれるが、この
高周波ノイズについての差信号Ｓsub は抽出した構造物
を表す差信号Ｓsub に対して極めて小さい信号値であ
る。
きさの構造物の中には高周波ノイズも含まれるが、この
高周波ノイズについての差信号Ｓsub は抽出した構造物
を表す差信号Ｓsub に対して極めて小さい信号値であ
る。
【００２５】次いでこの得られた差信号Ｓsub に応じて
変化する強調係数β（Ｓsub ）を設定する。この強調係
数β（Ｓsub ）としては、差信号Ｓsub の絶対値｜Ｓsu
b ｜が小さい範囲、具体的にはその値が高周波ノイズを
示すものである範囲においては出力（β（Ｓsub ））が
略ゼロとなり、高周波ノイズ以外の構造物を示すもので
ある範囲においてはある程度大きな値を出力するもので
あれば線形であると非線形であるとを問わない。
変化する強調係数β（Ｓsub ）を設定する。この強調係
数β（Ｓsub ）としては、差信号Ｓsub の絶対値｜Ｓsu
b ｜が小さい範囲、具体的にはその値が高周波ノイズを
示すものである範囲においては出力（β（Ｓsub ））が
略ゼロとなり、高周波ノイズ以外の構造物を示すもので
ある範囲においてはある程度大きな値を出力するもので
あれば線形であると非線形であるとを問わない。
【００２６】一方、オリジナル画像信号Ｓorg からボケ
マスク信号Ｓusを減算して比較的高い周波数成分Ｓspを
求める。この周波数成分Ｓspには高周波成分であるいわ
ゆる放射線ノイズ成分等も含まれる。
マスク信号Ｓusを減算して比較的高い周波数成分Ｓspを
求める。この周波数成分Ｓspには高周波成分であるいわ
ゆる放射線ノイズ成分等も含まれる。
【００２７】このようにして得られた比較的高い周波数
成分Ｓspに対して、画素ごとに対応する強調係数β（Ｓ
sub ）を乗じることにより、高周波ノイズ成分等につい
ては、周波数成分Ｓspはある程度の大きさの値を示すも
のの、強調係数β（Ｓsub ）が略ゼロであるため、これ
らの積の値β（Ｓsub ）×Ｓspは極めてゼロに近い小さ
な値となり、一方、この高周波ノイズ以外の所定の大き
さまでの構造物については、周波数成分Ｓspがある程度
の大きな値を示すとともに、強調係数β（Ｓsub ）もあ
る程度の大きな値を有するため、これらの積の値β（Ｓ
sub ）×Ｓspもある程度大きな値を示すこととなり、こ
の構造物は高周波ノイズ成分等に対して相対的に強調さ
れることとなる。
成分Ｓspに対して、画素ごとに対応する強調係数β（Ｓ
sub ）を乗じることにより、高周波ノイズ成分等につい
ては、周波数成分Ｓspはある程度の大きさの値を示すも
のの、強調係数β（Ｓsub ）が略ゼロであるため、これ
らの積の値β（Ｓsub ）×Ｓspは極めてゼロに近い小さ
な値となり、一方、この高周波ノイズ以外の所定の大き
さまでの構造物については、周波数成分Ｓspがある程度
の大きな値を示すとともに、強調係数β（Ｓsub ）もあ
る程度の大きな値を有するため、これらの積の値β（Ｓ
sub ）×Ｓspもある程度大きな値を示すこととなり、こ
の構造物は高周波ノイズ成分等に対して相対的に強調さ
れることとなる。
【００２８】このように本発明の画像処理方法および装
置によれば、画像中の高周波ノイズ成分等を強調するこ
となく、所望の構造物を選択的に強調処理することがで
きる。
置によれば、画像中の高周波ノイズ成分等を強調するこ
となく、所望の構造物を選択的に強調処理することがで
きる。
【００２９】また本発明の画像処理方法および装置で
は、メディアン値Ｓmed に依存したボケマスク処理を施
すことによりオーバーシュートやアンダーシュートを低
減することができ、その結果アーティファクトの発生を
抑制することもできる。
は、メディアン値Ｓmed に依存したボケマスク処理を施
すことによりオーバーシュートやアンダーシュートを低
減することができ、その結果アーティファクトの発生を
抑制することもできる。
【００３０】
【発明の実施の形態】以下、本発明の画像処理方法を具
体的に実施するための画像処理装置の実施の形態につい
て図面を用いて説明する。
体的に実施するための画像処理装置の実施の形態につい
て図面を用いて説明する。
【００３１】図１は本発明の画像処理装置の第１の実施
形態を示す概略ブロック図である。図示の画像処理装置
は、画像を表すオリジナル画像信号Ｓorg の、超低空間
周波数に対応する非鮮鋭マスク信号Ｓusを求めるローパ
スフィルタ11と、オリジナル画像信号Ｓorg から非鮮鋭
マスク信号Ｓusを減算処理してこれらの差信号Ｓsp（＝
Ｓorg −Ｓus）を求める演算素子14a と、所定の大きさ
のマスク内に含まれるオリジナル画像信号Ｓorg のメデ
ィアン値Ｓmed を求めるメディアンフィルタ12と、メデ
ィアン値Ｓmed からオリジナル画像信号Ｓorg を減算処
理してこれらの差信号Ｓsub （＝Ｓmed −Ｓorg ）を求
める演算素子14b と、この差信号Ｓsubを、差信号Ｓsub
に対応付けされた強調係数β（Ｓsub ）に変換する変
換テーブルを記憶してなる記憶手段13と、オリジナル画
像信号Ｓorg に対して、演算素子14a により得られたオ
リジナル画像信号Ｓorg の高周波成分Ｓspおよび変換テ
ーブル13により出力された強調係数β（Ｓsub ）を用い
て、下記式（１）にしたがった演算処理を施して処理済
画像信号Ｓprocを得る演算素子14c ，14d とを備えた構
成である。
形態を示す概略ブロック図である。図示の画像処理装置
は、画像を表すオリジナル画像信号Ｓorg の、超低空間
周波数に対応する非鮮鋭マスク信号Ｓusを求めるローパ
スフィルタ11と、オリジナル画像信号Ｓorg から非鮮鋭
マスク信号Ｓusを減算処理してこれらの差信号Ｓsp（＝
Ｓorg −Ｓus）を求める演算素子14a と、所定の大きさ
のマスク内に含まれるオリジナル画像信号Ｓorg のメデ
ィアン値Ｓmed を求めるメディアンフィルタ12と、メデ
ィアン値Ｓmed からオリジナル画像信号Ｓorg を減算処
理してこれらの差信号Ｓsub （＝Ｓmed −Ｓorg ）を求
める演算素子14b と、この差信号Ｓsubを、差信号Ｓsub
に対応付けされた強調係数β（Ｓsub ）に変換する変
換テーブルを記憶してなる記憶手段13と、オリジナル画
像信号Ｓorg に対して、演算素子14a により得られたオ
リジナル画像信号Ｓorg の高周波成分Ｓspおよび変換テ
ーブル13により出力された強調係数β（Ｓsub ）を用い
て、下記式（１）にしたがった演算処理を施して処理済
画像信号Ｓprocを得る演算素子14c ，14d とを備えた構
成である。
【００３２】 Ｓproc＝Ｓorg ＋β（Ｓsub ）×Ｓsp （１） ここで画像を表すオリジナル画像信号Ｓorg は、予め所
定の画像読取装置により放射線画像から読み取られて所
定の記憶手段に記憶されたものであってもよいし、画像
読取装置から直接入力されたものであってもよい。ま
た、本実施例における放射線画像としてはマンモグラム
とし、強調処理しようとする構造物としては乳がんの一
態様である微小石灰化陰影とする。
定の画像読取装置により放射線画像から読み取られて所
定の記憶手段に記憶されたものであってもよいし、画像
読取装置から直接入力されたものであってもよい。ま
た、本実施例における放射線画像としてはマンモグラム
とし、強調処理しようとする構造物としては乳がんの一
態様である微小石灰化陰影とする。
【００３３】ローパスフィルタ11はオリジナル画像信号
Ｓorg に対して、例えば３列×３行の画素マトリクスか
らなるボケマスクを設定し、下記式（５）（Ｎ＝３に設
定）にしたがって得られたボケマスク信号Ｓusを出力す
る。
Ｓorg に対して、例えば３列×３行の画素マトリクスか
らなるボケマスクを設定し、下記式（５）（Ｎ＝３に設
定）にしたがって得られたボケマスク信号Ｓusを出力す
る。
【００３４】 Ｓus＝（ΣＳorg ）／Ｎ2 （５） なお、ボケマスクとしては上記式（５）に示すようにマ
スク内の画素値の単純平均を用いるものの他、例えば図
３のマトリクスに示したような、中心画素からの距離に
応じてマスク内の画素値の重み付けを変化させたものを
用いることもできる。
スク内の画素値の単純平均を用いるものの他、例えば図
３のマトリクスに示したような、中心画素からの距離に
応じてマスク内の画素値の重み付けを変化させたものを
用いることもできる。
【００３５】メディアンフィルタ12はオリジナル画像信
号Ｓorg に対して、例えば９列×９行（＝81画素）の画
素マトリクスからなるマスクを設定し、このマスク内に
含まれる81画素分のオリジナル画像信号Ｓorg のうちの
メディアン値（中央値；累積ヒストグラムで大きい方か
ら、または小さい方から41番目の画像信号値）Ｓmedを
出力する。したがってメディアン値Ｓmed は、オリジナ
ル画像信号Ｓorg のうちメディアンフィルタ12のマスク
サイズの１／２未満の大きさ（40画素以下）の構造物に
ついてはこの構造物の存在を示す特徴的な値を出力する
ことがなく、一方、メディアンフィルタ12のサイズの１
／２以上の大きさ（41画素以上）の構造物についてはこ
の構造物の存在を示す特徴的な値を出力することとな
る。
号Ｓorg に対して、例えば９列×９行（＝81画素）の画
素マトリクスからなるマスクを設定し、このマスク内に
含まれる81画素分のオリジナル画像信号Ｓorg のうちの
メディアン値（中央値；累積ヒストグラムで大きい方か
ら、または小さい方から41番目の画像信号値）Ｓmedを
出力する。したがってメディアン値Ｓmed は、オリジナ
ル画像信号Ｓorg のうちメディアンフィルタ12のマスク
サイズの１／２未満の大きさ（40画素以下）の構造物に
ついてはこの構造物の存在を示す特徴的な値を出力する
ことがなく、一方、メディアンフィルタ12のサイズの１
／２以上の大きさ（41画素以上）の構造物についてはこ
の構造物の存在を示す特徴的な値を出力することとな
る。
【００３６】なおローパスフィルタ11によるボケマスク
のサイズ、メディアンフィルタ12によるマスクのサイズ
は、40画素以下の微小石灰化陰影を抽出するのに適正で
あることが予め実験的に明らかになっていることより上
記の通り定めたものであるが、画像全体の拡大率が異な
る場合や他の大きさの構造物を抽出する場合には、抽出
しようとする構造物に適したマスクサイズを実験的に予
め設定するのが望ましい。本実施例におけるメディアン
フィルタ12は、そのマスクサイズの１／２未満の大きさ
である40画素以下の石灰化陰影を抽出することを目的と
している。
のサイズ、メディアンフィルタ12によるマスクのサイズ
は、40画素以下の微小石灰化陰影を抽出するのに適正で
あることが予め実験的に明らかになっていることより上
記の通り定めたものであるが、画像全体の拡大率が異な
る場合や他の大きさの構造物を抽出する場合には、抽出
しようとする構造物に適したマスクサイズを実験的に予
め設定するのが望ましい。本実施例におけるメディアン
フィルタ12は、そのマスクサイズの１／２未満の大きさ
である40画素以下の石灰化陰影を抽出することを目的と
している。
【００３７】また記憶手段13に記憶された変換テーブル
は、例えば図２（１）の関数曲線に示すように、差信号
Ｓsub の絶対値｜Ｓsub ｜がゼロから所定の閾値Ｓ0 ま
での範囲内においては、この絶対値｜Ｓsub ｜が増大す
るにしたがって強調係数β（Ｓsub ）が単調に増大し、
その所定の閾値Ｓ0 より大きい範囲においては絶対値｜
Ｓsub ｜に拘らず最大値βmax （一定）を採るものであ
る。
は、例えば図２（１）の関数曲線に示すように、差信号
Ｓsub の絶対値｜Ｓsub ｜がゼロから所定の閾値Ｓ0 ま
での範囲内においては、この絶対値｜Ｓsub ｜が増大す
るにしたがって強調係数β（Ｓsub ）が単調に増大し、
その所定の閾値Ｓ0 より大きい範囲においては絶対値｜
Ｓsub ｜に拘らず最大値βmax （一定）を採るものであ
る。
【００３８】次に本実施形態の画像処理装置の作用につ
いて説明する。
いて説明する。
【００３９】画像処理装置にオリジナル画像信号Ｓorg
が入力されると、まずローパスフィルタ11はオリジナル
画像信号Ｓorg に対して、３列×３行の画素マトリクス
からなるボケマスクを設定し、上記式（５）（Ｎ＝３）
にしたがった演算を施してボケマスク信号Ｓusを出力す
る。このボケマスク信号Ｓusは３列×３行の画素マトリ
クスからなるボケマスクであるため、比較的高い周波数
成分を含む信号となっている。
が入力されると、まずローパスフィルタ11はオリジナル
画像信号Ｓorg に対して、３列×３行の画素マトリクス
からなるボケマスクを設定し、上記式（５）（Ｎ＝３）
にしたがった演算を施してボケマスク信号Ｓusを出力す
る。このボケマスク信号Ｓusは３列×３行の画素マトリ
クスからなるボケマスクであるため、比較的高い周波数
成分を含む信号となっている。
【００４０】演算素子14a は、オリジナル画像信号Ｓor
g からローパスフィルタ11から出力された比較的高い周
波数を含む信号であるボケマスク信号Ｓusを減算するこ
とにより、オリジナル画像信号Ｓorg のうち高周波成分
Ｓsp（＝Ｓorg −Ｓus）だけを出力する。
g からローパスフィルタ11から出力された比較的高い周
波数を含む信号であるボケマスク信号Ｓusを減算するこ
とにより、オリジナル画像信号Ｓorg のうち高周波成分
Ｓsp（＝Ｓorg −Ｓus）だけを出力する。
【００４１】一方、メディアンフィルタ12はオリジナル
画像信号Ｓorg に対して、９列×９行の画素マトリクス
からなるマスクを設定し、このマスク内に含まれる81画
素分のオリジナル画像信号Ｓorg のうちのメディアン値
Ｓmed を出力する。
画像信号Ｓorg に対して、９列×９行の画素マトリクス
からなるマスクを設定し、このマスク内に含まれる81画
素分のオリジナル画像信号Ｓorg のうちのメディアン値
Ｓmed を出力する。
【００４２】次いで演算素子14b が、メディアン値Ｓme
d からオリジナル画像信号Ｓorg を減算することによ
り、オリジナル画像信号Ｓorg のうち40画素以下の大き
さの構造物についてはその構造物の存在を示す特徴的な
値、すなわちゼロ以外の値を、41画素以上の大きさの構
造物についてはゼロを採る差信号Ｓsub を出力する。
d からオリジナル画像信号Ｓorg を減算することによ
り、オリジナル画像信号Ｓorg のうち40画素以下の大き
さの構造物についてはその構造物の存在を示す特徴的な
値、すなわちゼロ以外の値を、41画素以上の大きさの構
造物についてはゼロを採る差信号Ｓsub を出力する。
【００４３】メディアンフィルタ12から出力された差信
号Ｓsub は記憶手段13に入力され、この差信号Ｓsub は
記憶手段13に記憶された変換テーブルにより、その差信
号Ｓsub の大きさに応じた強調係数β（Ｓsub ）に変換
されて記憶手段13から出力される。
号Ｓsub は記憶手段13に入力され、この差信号Ｓsub は
記憶手段13に記憶された変換テーブルにより、その差信
号Ｓsub の大きさに応じた強調係数β（Ｓsub ）に変換
されて記憶手段13から出力される。
【００４４】ここで差信号Ｓsub は上述したように、40
画素以下の大きさの構造物についてはゼロ以外の値を、
41画素以上の大きさの構造物についてはゼロを採るた
め、記憶手段13に記憶された変換テーブルにより出力さ
れる強調係数β（Ｓsub ）の値としては、41画素以上の
大きさの構造物についてはゼロを示し、40画素以下の大
きさの構造物についてはゼロ以外の値を示す。また構造
物についてのオリジナル画像信号Ｓorg はその構造物の
周囲の画像信号に対してある程度の信号値差があるため
上記差信号Ｓsub もある程度大きな値を示すが、一方、
高周波ノイズについてのオリジナル画像信号Ｓorg はそ
のノイズの周囲の画像信号に対して大きな信号値差がな
いため上記差信号Ｓsub も小さい値、具体的には上記閾
値Ｓ0 よりもかなり小さい値となる。
画素以下の大きさの構造物についてはゼロ以外の値を、
41画素以上の大きさの構造物についてはゼロを採るた
め、記憶手段13に記憶された変換テーブルにより出力さ
れる強調係数β（Ｓsub ）の値としては、41画素以上の
大きさの構造物についてはゼロを示し、40画素以下の大
きさの構造物についてはゼロ以外の値を示す。また構造
物についてのオリジナル画像信号Ｓorg はその構造物の
周囲の画像信号に対してある程度の信号値差があるため
上記差信号Ｓsub もある程度大きな値を示すが、一方、
高周波ノイズについてのオリジナル画像信号Ｓorg はそ
のノイズの周囲の画像信号に対して大きな信号値差がな
いため上記差信号Ｓsub も小さい値、具体的には上記閾
値Ｓ0 よりもかなり小さい値となる。
【００４５】演算素子14c は、演算素子14a からの出力
Ｓspと記憶手段13からの出力β（Ｓsub ）とを乗算処理
して高周波成分Ｓspに強調係数β（Ｓsub ）で重み付け
を行う。ここで、強調係数β（Ｓsub ）は前述したよう
に、所定の大きさ以上の構造物についてはゼロを示すか
らそのような大きさの構造物は全く強調処理されること
がなく、また所定の大きさ未満の構造物であっても、そ
れが放射線ノイズである場合には極めて小さい値を示す
ためこの放射線ノイズはほとんど強調処理されず、放射
線ノイズ以外の所定の大きさ以下の構造物だけが選択的
に強調処理される。
Ｓspと記憶手段13からの出力β（Ｓsub ）とを乗算処理
して高周波成分Ｓspに強調係数β（Ｓsub ）で重み付け
を行う。ここで、強調係数β（Ｓsub ）は前述したよう
に、所定の大きさ以上の構造物についてはゼロを示すか
らそのような大きさの構造物は全く強調処理されること
がなく、また所定の大きさ未満の構造物であっても、そ
れが放射線ノイズである場合には極めて小さい値を示す
ためこの放射線ノイズはほとんど強調処理されず、放射
線ノイズ以外の所定の大きさ以下の構造物だけが選択的
に強調処理される。
【００４６】次いで演算素子14d は、オリジナル画像信
号Ｓorg と演算素子14c からの出力とを加算処理して処
理済画像信号Ｓprocを出力する。この出力された処理済
画像信号Ｓprocは、上述の各処理によりオリジナル画像
信号の高周波成分Ｓspのうちの微小石灰化陰影を含む所
定の大きさ以下の構造物だけが選択的に強調処理された
信号となる。
号Ｓorg と演算素子14c からの出力とを加算処理して処
理済画像信号Ｓprocを出力する。この出力された処理済
画像信号Ｓprocは、上述の各処理によりオリジナル画像
信号の高周波成分Ｓspのうちの微小石灰化陰影を含む所
定の大きさ以下の構造物だけが選択的に強調処理された
信号となる。
【００４７】図４に、従来のボケマスク処理に基づく強
調処理（式（４）参照）、メディアン値処理に基づく強
調処理（式（６）参照）、および本実施形態の式（１）
に基づく強調処理によってそれぞれ得られる処理済画像
信号Ｓprocを示す。これによれば、本実施形態の画像処
理装置では、高周波成分Ｓspの周波数帯域内にある放射
線ノイズの強調を抑制しつつ、所定の大きさ以下の微小
石灰化部分等の特徴的な画像部分を選択的に効率よく強
調処理することができ、計算機支援画像診断等において
も非常に有用である。
調処理（式（４）参照）、メディアン値処理に基づく強
調処理（式（６）参照）、および本実施形態の式（１）
に基づく強調処理によってそれぞれ得られる処理済画像
信号Ｓprocを示す。これによれば、本実施形態の画像処
理装置では、高周波成分Ｓspの周波数帯域内にある放射
線ノイズの強調を抑制しつつ、所定の大きさ以下の微小
石灰化部分等の特徴的な画像部分を選択的に効率よく強
調処理することができ、計算機支援画像診断等において
も非常に有用である。
【００４８】なお、記憶手段13に記憶された変換テーブ
ルは図２（１）に示すものの他、同図（２）、同図
（３）に示す関数形状のものとすることもできる。すな
わち図２（１）に示すものが周囲より低濃度な構造物を
抽出するのに適合した変換テーブルであるのに対し、同
図（２）に示すものは周囲より低濃度な構造物および高
濃度な構造物の双方を抽出するのに適した変換テーブル
であり、同図（３）に示すものは、所定の閾値Ｓ0 より
大きい範囲においては絶対値｜Ｓsub ｜の値が増大する
のにしたがって強調係数β（Ｓsub ）の値が減少するよ
うに設定されたものであり、この（３）に示す変換テー
ブルによれば差信号Ｓsub の値が極めて大きい領域部
分、例えばコントラストの大きい画像部分については強
調係数β（Ｓsub ）を抑制してオーバーシュートやアン
ダーシュートの発生を抑止し、その結果偽輪郭の発生を
抑制する効果をより向上させることができる。
ルは図２（１）に示すものの他、同図（２）、同図
（３）に示す関数形状のものとすることもできる。すな
わち図２（１）に示すものが周囲より低濃度な構造物を
抽出するのに適合した変換テーブルであるのに対し、同
図（２）に示すものは周囲より低濃度な構造物および高
濃度な構造物の双方を抽出するのに適した変換テーブル
であり、同図（３）に示すものは、所定の閾値Ｓ0 より
大きい範囲においては絶対値｜Ｓsub ｜の値が増大する
のにしたがって強調係数β（Ｓsub ）の値が減少するよ
うに設定されたものであり、この（３）に示す変換テー
ブルによれば差信号Ｓsub の値が極めて大きい領域部
分、例えばコントラストの大きい画像部分については強
調係数β（Ｓsub ）を抑制してオーバーシュートやアン
ダーシュートの発生を抑止し、その結果偽輪郭の発生を
抑制する効果をより向上させることができる。
【００４９】図５は本発明の画像処理装置の第２の実施
形態を示す概略ブロック図である。図示の画像処理装置
は、画像を表すオリジナル画像信号Ｓorg の超低空間周
波数に対応する非鮮鋭マスク信号Ｓusを求める第１のロ
ーパスフィルタ11a ，第２のローパスフィルタ11b ，第
３のローパスフィルタ11c と、これらの各ローパスフィ
ルタ11a ，11b ，11c の入出力信号間で差の演算を行い
オリジナル画像信号Ｓorg の高周波帯域ｆH に対応する
高周波成分ＳH 、中周波帯域ｆM に対応する中周波成分
ＳM 、低周波帯域ｆL に対応する低周波成分ＳL を求め
る演算素子14ａ1 ，14ａ2 ，14ａ3 と、画素マトリクス
として最も小さいサイズのマスク内に含まれるオリジナ
ル画像信号Ｓorg の第１のメディアン値Ｓmed 1 を求め
る第１のメディアンフィルタ12a ，画素マトリクスとし
て最も大きいサイズのマスク内に含まれるオリジナル画
像信号Ｓorg の第３のメディアン値Ｓmed 3 を求める第
３のメディアンフィルタ12c ，および第１のメディアン
フィルタ12a と第３のメディアンフィルタ12c との中間
サイズのマスク内に含まれるオリジナル画像信号Ｓorg
の第２のメディアン値Ｓmed 2 を求める第２のメディア
ンフィルタ12b と、各メディアンフィルタ12a ，12b ，
12c の出力である第１のメディアン値Ｓmed1 ，第２の
メディアン値Ｓmed 2 ，第３のメディアン値Ｓmed 3 か
らそれぞれオリジナル画像信号Ｓorg を減算処理してこ
れらの差信号Ｓsub 1 ，Ｓsub 2 ，Ｓsub 3 を求める演
算素子14ｂ1 ，14ｂ2 ，14ｂ3 と、この各差信号Ｓsub
1 ，Ｓsub 2 ，Ｓsub 3 を、各差信号Ｓsub 1 ，Ｓsub
2 ，Ｓsub 3 にそれぞれ対応付けされた強調係数β
1 （Ｓsub 1 ），β2 （Ｓsub 2 ），β3 （Ｓsub 3 ）
に変換する変換テーブルをそれぞれ各別に記憶してなる
記憶手段13a ，13b ，13c と、オリジナル画像信号Ｓor
g に対して、ローパスフィルタ11a ，11b ，11c の出力
および記憶手段13a ，13b ，13c の出力を用いて下記式
（３）にしたがった演算処理を施して処理済画像信号Ｓ
procを得る演算素子14ｃ1 ，14ｃ2 ，14ｃ3 ，14ｄ1 ，
14ｄ2 ，14ｄ3 とを備えた構成である。
形態を示す概略ブロック図である。図示の画像処理装置
は、画像を表すオリジナル画像信号Ｓorg の超低空間周
波数に対応する非鮮鋭マスク信号Ｓusを求める第１のロ
ーパスフィルタ11a ，第２のローパスフィルタ11b ，第
３のローパスフィルタ11c と、これらの各ローパスフィ
ルタ11a ，11b ，11c の入出力信号間で差の演算を行い
オリジナル画像信号Ｓorg の高周波帯域ｆH に対応する
高周波成分ＳH 、中周波帯域ｆM に対応する中周波成分
ＳM 、低周波帯域ｆL に対応する低周波成分ＳL を求め
る演算素子14ａ1 ，14ａ2 ，14ａ3 と、画素マトリクス
として最も小さいサイズのマスク内に含まれるオリジナ
ル画像信号Ｓorg の第１のメディアン値Ｓmed 1 を求め
る第１のメディアンフィルタ12a ，画素マトリクスとし
て最も大きいサイズのマスク内に含まれるオリジナル画
像信号Ｓorg の第３のメディアン値Ｓmed 3 を求める第
３のメディアンフィルタ12c ，および第１のメディアン
フィルタ12a と第３のメディアンフィルタ12c との中間
サイズのマスク内に含まれるオリジナル画像信号Ｓorg
の第２のメディアン値Ｓmed 2 を求める第２のメディア
ンフィルタ12b と、各メディアンフィルタ12a ，12b ，
12c の出力である第１のメディアン値Ｓmed1 ，第２の
メディアン値Ｓmed 2 ，第３のメディアン値Ｓmed 3 か
らそれぞれオリジナル画像信号Ｓorg を減算処理してこ
れらの差信号Ｓsub 1 ，Ｓsub 2 ，Ｓsub 3 を求める演
算素子14ｂ1 ，14ｂ2 ，14ｂ3 と、この各差信号Ｓsub
1 ，Ｓsub 2 ，Ｓsub 3 を、各差信号Ｓsub 1 ，Ｓsub
2 ，Ｓsub 3 にそれぞれ対応付けされた強調係数β
1 （Ｓsub 1 ），β2 （Ｓsub 2 ），β3 （Ｓsub 3 ）
に変換する変換テーブルをそれぞれ各別に記憶してなる
記憶手段13a ，13b ，13c と、オリジナル画像信号Ｓor
g に対して、ローパスフィルタ11a ，11b ，11c の出力
および記憶手段13a ，13b ，13c の出力を用いて下記式
（３）にしたがった演算処理を施して処理済画像信号Ｓ
procを得る演算素子14ｃ1 ，14ｃ2 ，14ｃ3 ，14ｄ1 ，
14ｄ2 ，14ｄ3 とを備えた構成である。
【００５０】 Ｓproc＝Ｓorg ＋Σ｛Ｓn ×βn （Ｓsub n ）｝ （３） ここで、各ローパスフィルタ11a ，11b ，11c のマスク
サイズと各メディアンフィルタ12a ，12b ，12c のマス
クサイズとは次の関係を有するように設定されている。
サイズと各メディアンフィルタ12a ，12b ，12c のマス
クサイズとは次の関係を有するように設定されている。
【００５１】すなわち、第１のメディアンフィルタ12a
のマスクサイズは第１のローパスフィルタ11a のマスク
サイズの２倍以上、第２のメディアンフィルタ12b のマ
スクサイズは第１のローパスフィルタ11a と第２のロー
パスフィルタ11b との各マスクサイズの積の２倍以上、
第３のメディアンフィルタ12c のマスクサイズは第１の
ローパスフィルタ11a と第２のローパスフィルタ11b と
第３のローパスフィルタ11c との各マスクサイズの積の
２倍以上である。
のマスクサイズは第１のローパスフィルタ11a のマスク
サイズの２倍以上、第２のメディアンフィルタ12b のマ
スクサイズは第１のローパスフィルタ11a と第２のロー
パスフィルタ11b との各マスクサイズの積の２倍以上、
第３のメディアンフィルタ12c のマスクサイズは第１の
ローパスフィルタ11a と第２のローパスフィルタ11b と
第３のローパスフィルタ11c との各マスクサイズの積の
２倍以上である。
【００５２】なお、第１のローパスフィルタ11a と第２
のローパスフィルタ11b と第３のローパスフィルタ11c
とは必ずしも互いに異なるマスクサイズである必要はな
く、例えばすべて第１のローパスフィルタ11a と同じも
のを用いてもよい。これはローパスフィルタを直列に接
続した場合、各ローパスフィルタを通過する毎に、より
低周波成分だけが出力されるからである。
のローパスフィルタ11b と第３のローパスフィルタ11c
とは必ずしも互いに異なるマスクサイズである必要はな
く、例えばすべて第１のローパスフィルタ11a と同じも
のを用いてもよい。これはローパスフィルタを直列に接
続した場合、各ローパスフィルタを通過する毎に、より
低周波成分だけが出力されるからである。
【００５３】また上記各ローパスフィルタ11a ，11b ，
11c を通過した後の出力信号をＳus1 ，Ｓus2 ，Ｓus3
とする。
11c を通過した後の出力信号をＳus1 ，Ｓus2 ，Ｓus3
とする。
【００５４】さらに上記記憶手段13a ，13b ，13c にそ
れぞれ記憶された各変換テーブルは前述の第１の実施の
形態における変換テーブルと同様に、例えば図２
（１）、同図（２）、同図（３）に示す関数形状を有す
るものである。
れぞれ記憶された各変換テーブルは前述の第１の実施の
形態における変換テーブルと同様に、例えば図２
（１）、同図（２）、同図（３）に示す関数形状を有す
るものである。
【００５５】次に本実施の形態の画像処理装置の作用に
ついて説明する。
ついて説明する。
【００５６】画像処理装置にオリジナル画像信号Ｓorg
が入力されると、オリジナル画像信号Ｓorg は第１のロ
ーパスフィルタ11a から順次第２のローパスフィルタ11
b 、第３のローパスフィルタ11c に入力される。第１の
ローパスフィルタ11a は入力されたオリジナル画像信号
Ｓorg に対して、例えば３列×３行の画素マトリクスか
らなるボケマスクを設定し、上記式（５）にしたがった
演算を施して第１のボケマスク信号Ｓus１を出力する。
この第１のボケマスク信号Ｓus１は３列×３行の画素マ
トリクスからなるボケマスクによるものであるため、所
定の高周波成分ＳH が除去された信号成分であり、演算
素子14ａ1 がオリジナル画像信号Ｓorgからこの第１の
ボケマスク信号Ｓus１を差し引く演算をなすことによ
り、高周波成分ＳH だけが抽出される。
が入力されると、オリジナル画像信号Ｓorg は第１のロ
ーパスフィルタ11a から順次第２のローパスフィルタ11
b 、第３のローパスフィルタ11c に入力される。第１の
ローパスフィルタ11a は入力されたオリジナル画像信号
Ｓorg に対して、例えば３列×３行の画素マトリクスか
らなるボケマスクを設定し、上記式（５）にしたがった
演算を施して第１のボケマスク信号Ｓus１を出力する。
この第１のボケマスク信号Ｓus１は３列×３行の画素マ
トリクスからなるボケマスクによるものであるため、所
定の高周波成分ＳH が除去された信号成分であり、演算
素子14ａ1 がオリジナル画像信号Ｓorgからこの第１の
ボケマスク信号Ｓus１を差し引く演算をなすことによ
り、高周波成分ＳH だけが抽出される。
【００５７】第１のローパスフィルタ11a から出力され
た第１のボケマスク信号Ｓus１は第２のローパスフィル
タ11b にも入力され、第２のローパスフィルタ11b は第
１のローパスフィルタ11a と同様の作用により、第１の
ボケマスク信号Ｓus１のうちの高い周波数成分である中
周波数成分ＳM が除去された第２のボケマスク信号Ｓus
２を出力し、演算素子14ａ2 が第１のボケマスク信号Ｓ
us１からこの第２のボケマスク信号Ｓus２を差し引く演
算をなすことにより、中周波成分ＳM だけが抽出され
る。
た第１のボケマスク信号Ｓus１は第２のローパスフィル
タ11b にも入力され、第２のローパスフィルタ11b は第
１のローパスフィルタ11a と同様の作用により、第１の
ボケマスク信号Ｓus１のうちの高い周波数成分である中
周波数成分ＳM が除去された第２のボケマスク信号Ｓus
２を出力し、演算素子14ａ2 が第１のボケマスク信号Ｓ
us１からこの第２のボケマスク信号Ｓus２を差し引く演
算をなすことにより、中周波成分ＳM だけが抽出され
る。
【００５８】第３のローパスフィルタ11c および演算素
子14ａ3 も同様にして低周波成分ＳL だけを抽出する。
なお、第３のローパスフィルタ11c からの出力である第
３のボケマスク信号Ｓus３は、低周波成分ＳL に含まれ
ない最も低い周波数成分であり、これを超低周波成分Ｓ
LL（＝Ｓus３）とする。
子14ａ3 も同様にして低周波成分ＳL だけを抽出する。
なお、第３のローパスフィルタ11c からの出力である第
３のボケマスク信号Ｓus３は、低周波成分ＳL に含まれ
ない最も低い周波数成分であり、これを超低周波成分Ｓ
LL（＝Ｓus３）とする。
【００５９】以上の作用によりオリジナル画像信号Ｓor
g は４つの周波数帯域の周波数成分ＳH ，ＳM ，ＳL ，
ＳLLに分割される。
g は４つの周波数帯域の周波数成分ＳH ，ＳM ，ＳL ，
ＳLLに分割される。
【００６０】なお、本実施例においては各周波数成分Ｓ
H ，ＳM ，ＳL ，ＳLLは周波数帯域が互いに重複しない
ように分割されるが、本発明においては必ずしも各周波
数帯域間で重複帯域がないように分割する必要はなく、
一部に重複帯域を含むように分割するものであってもよ
い。
H ，ＳM ，ＳL ，ＳLLは周波数帯域が互いに重複しない
ように分割されるが、本発明においては必ずしも各周波
数帯域間で重複帯域がないように分割する必要はなく、
一部に重複帯域を含むように分割するものであってもよ
い。
【００６１】一方、オリジナル画像信号Ｓorg は第１，
第２，第３の各メディアンフィルタ12a ，12b ，12c に
も入力される。各メディアンフィルタ12a ，12b ，12c
はこの入力されたオリジナル画像信号Ｓorg に対してそ
れぞれ前述した大きさのマスクを設定したうえで、各マ
スク内のメディアン値Ｓmed 1 ，Ｓmed 2 ，Ｓmed 3を
算出する。
第２，第３の各メディアンフィルタ12a ，12b ，12c に
も入力される。各メディアンフィルタ12a ，12b ，12c
はこの入力されたオリジナル画像信号Ｓorg に対してそ
れぞれ前述した大きさのマスクを設定したうえで、各マ
スク内のメディアン値Ｓmed 1 ，Ｓmed 2 ，Ｓmed 3を
算出する。
【００６２】これらの各メディアンフィルタ12a ，12b
，12c から出力された各メディアン値Ｓmed 1 ，Ｓmed
2 ，Ｓmed 3 は、それぞれ演算素子14ｂ1 ，14ｂ2 ，1
4ｂ3に入力されて、各メディアン値Ｓmed 1 ，Ｓmed
2 ，Ｓmed 3 からオリジナル画像信号Ｓorg がそれぞ
れ減算処理されて差信号Ｓsub 1 ，Ｓsub 2 ，Ｓsub 3
が出力される。
，12c から出力された各メディアン値Ｓmed 1 ，Ｓmed
2 ，Ｓmed 3 は、それぞれ演算素子14ｂ1 ，14ｂ2 ，1
4ｂ3に入力されて、各メディアン値Ｓmed 1 ，Ｓmed
2 ，Ｓmed 3 からオリジナル画像信号Ｓorg がそれぞ
れ減算処理されて差信号Ｓsub 1 ，Ｓsub 2 ，Ｓsub 3
が出力される。
【００６３】ここで、第１の差信号Ｓsub 1 は、高周波
成分ＳH に対応する大きさ以下の構造物についての信号
値に対して特徴的な値を示し、第２の差信号Ｓsub
2 は、中周波成分ＳM に対応する大きさ以下の構造物に
ついての信号値に対して特徴的な値を示し、第３の差信
号Ｓsub 3 は、低周波成分ＳL に対応する大きさ以下の
構造物についての信号値に対して特徴的な値を示す。
成分ＳH に対応する大きさ以下の構造物についての信号
値に対して特徴的な値を示し、第２の差信号Ｓsub
2 は、中周波成分ＳM に対応する大きさ以下の構造物に
ついての信号値に対して特徴的な値を示し、第３の差信
号Ｓsub 3 は、低周波成分ＳL に対応する大きさ以下の
構造物についての信号値に対して特徴的な値を示す。
【００６４】演算素子14ｂ1 ，14ｂ2 ，14ｂ3 から出力
された差信号Ｓsub 2 は、第１〜第３の記憶手段13a ，
13b ，13c にそれぞれ入力される。
された差信号Ｓsub 2 は、第１〜第３の記憶手段13a ，
13b ，13c にそれぞれ入力される。
【００６５】第１の記憶手段13a に記憶された第１の変
換テーブル13a は入力された第１の差信号Ｓsub 1 を第
１の強調係数β1 （Ｓsub 1 ）に変換して出力し、第２
の記憶手段13b に記憶された第１の変換テーブル13b は
入力された第２の差信号Ｓsub 2 を第２の強調係数β2
（Ｓsub 2 ）に変換して出力し、第３の記憶手段13cに
記憶された第３の変換テーブル13c は入力された第３の
差信号Ｓsub 3 を第３の強調係数β3 （Ｓsub 3 ）に変
換して出力する。
換テーブル13a は入力された第１の差信号Ｓsub 1 を第
１の強調係数β1 （Ｓsub 1 ）に変換して出力し、第２
の記憶手段13b に記憶された第１の変換テーブル13b は
入力された第２の差信号Ｓsub 2 を第２の強調係数β2
（Ｓsub 2 ）に変換して出力し、第３の記憶手段13cに
記憶された第３の変換テーブル13c は入力された第３の
差信号Ｓsub 3 を第３の強調係数β3 （Ｓsub 3 ）に変
換して出力する。
【００６６】なお、超低周波成分ＳLLに対応するメディ
アン値を仮にＳmed 4 、メディアン値Ｓmed 4 に対応す
る強調係数をβ4 （Ｓsub 4 ）とする第４の変換テーブ
ルを設けてもよいが、本実施形態においてはこの超低周
波成分ＳLLを個別に強調することを要しないため、β4
（Ｓsub 4 ）＝１（一定）とみなして省略している。す
なわち、ｎ個の周波数成分に分割された信号のうち（ｎ
−１）個の周波数成分の強調度を低くして抑制度とする
ことによって、相対的に残りの１周波数成分の強調度を
も制御することができるからである。
アン値を仮にＳmed 4 、メディアン値Ｓmed 4 に対応す
る強調係数をβ4 （Ｓsub 4 ）とする第４の変換テーブ
ルを設けてもよいが、本実施形態においてはこの超低周
波成分ＳLLを個別に強調することを要しないため、β4
（Ｓsub 4 ）＝１（一定）とみなして省略している。す
なわち、ｎ個の周波数成分に分割された信号のうち（ｎ
−１）個の周波数成分の強調度を低くして抑制度とする
ことによって、相対的に残りの１周波数成分の強調度を
も制御することができるからである。
【００６７】各記憶手段13a ，13b ，13c からの出力で
ある強調係数β1 （Ｓsub 1 ），β2 （Ｓsub 2 ），β
3 （Ｓsub 3 ）はそれぞれ対応する演算素子14ｃ1 ，14
ｃ2，14ｃ3 に入力される。一方、この演算素子14
ｃ1 ，14ｃ2 ，14ｃ3 には対応する演算素子14ａ1 ，14
ａ2 ，14ａ3 からそれぞれ前述の高周波成分ＳH 、中周
波成分ＳM 、低周波成分ＳL も入力され、演算素子14ｃ
1 ，14ｃ2 ，14ｃ3 はこれらを対応する組ごと（β
1 （Ｓsub 1 ）とＳH ，β2 （Ｓsub 2 ）とＳM ，β3
（Ｓsub 3 ）とＳL ）に乗算し、それぞれ演算素子14ｄ
1 ，14ｄ2 ，14ｄ3 に出力する。この作用は各周波数成
分を各周波数成分毎のメディアン値に応じた重み付けを
行うことを意味している。
ある強調係数β1 （Ｓsub 1 ），β2 （Ｓsub 2 ），β
3 （Ｓsub 3 ）はそれぞれ対応する演算素子14ｃ1 ，14
ｃ2，14ｃ3 に入力される。一方、この演算素子14
ｃ1 ，14ｃ2 ，14ｃ3 には対応する演算素子14ａ1 ，14
ａ2 ，14ａ3 からそれぞれ前述の高周波成分ＳH 、中周
波成分ＳM 、低周波成分ＳL も入力され、演算素子14ｃ
1 ，14ｃ2 ，14ｃ3 はこれらを対応する組ごと（β
1 （Ｓsub 1 ）とＳH ，β2 （Ｓsub 2 ）とＳM ，β3
（Ｓsub 3 ）とＳL ）に乗算し、それぞれ演算素子14ｄ
1 ，14ｄ2 ，14ｄ3 に出力する。この作用は各周波数成
分を各周波数成分毎のメディアン値に応じた重み付けを
行うことを意味している。
【００６８】演算素子14ｄ1 ，14ｄ2 ，14ｄ3 は演算素
子14ｃ1 ，14ｃ2 ，14ｃ3 からの出力を総和するととも
に、第３のローパスフィルタ11c からの出力である超低
周波成分ＳLLを加算し、上記式（３）に示す処理済画像
信号Ｓprocを出力する。
子14ｃ1 ，14ｃ2 ，14ｃ3 からの出力を総和するととも
に、第３のローパスフィルタ11c からの出力である超低
周波成分ＳLLを加算し、上記式（３）に示す処理済画像
信号Ｓprocを出力する。
【００６９】この処理済画像信号Ｓprocは、オリジナル
画像信号Ｓorg に対して、各周波数成分ごとにメディア
ン値に応じた強調処理を施した信号であり、上述のマス
クサイズ、分割する周波数帯域の数などを種々変化させ
ることにより、また変換テーブルによる強調関数の形状
を種々変えることにより、各周波数帯域に含まれる所望
の信号成分だけを選択的に所望の強調度合いで強調処理
を行うことができ、従来の強調処理と比べてよりきめ細
かい強調の調節が可能となり、画像診断性能を向上させ
ることができる。
画像信号Ｓorg に対して、各周波数成分ごとにメディア
ン値に応じた強調処理を施した信号であり、上述のマス
クサイズ、分割する周波数帯域の数などを種々変化させ
ることにより、また変換テーブルによる強調関数の形状
を種々変えることにより、各周波数帯域に含まれる所望
の信号成分だけを選択的に所望の強調度合いで強調処理
を行うことができ、従来の強調処理と比べてよりきめ細
かい強調の調節が可能となり、画像診断性能を向上させ
ることができる。
【図１】本発明の画像処理装置の第１の実施の形態を示
すブロック図
すブロック図
【図２】（１）〜（３）変換テーブルの具体例を示すグ
ラフ
ラフ
【図３】図１に示したローパスフィルタのマスクの重み
付け係数を示す図
付け係数を示す図
【図４】従来のボケマスク処理に基づく強調処理（式
（４）参照）、メディアン値処理に基づく強調処理（式
（６）参照）、および本実施形態の式（１）に基づく強
調処理によってそれぞれ得られる処理済画像信号Ｓproc
を示す図
（４）参照）、メディアン値処理に基づく強調処理（式
（６）参照）、および本実施形態の式（１）に基づく強
調処理によってそれぞれ得られる処理済画像信号Ｓproc
を示す図
【図５】本発明の画像処理装置の第２の実施の形態を示
すブロック図
すブロック図
【符号の説明】 11 ローパスフィルタ 12 メディアンフィルタ 13 記憶手段 14ａ，14ｂ，14ｃ，14ｄ 演算素子 Ｓorg オリジナル画像信号 Ｓmed メディアン信号 Ｓus ボケマスク信号 Ｓsp 高周波数成分 Ｓsub 差信号 β（Ｓsub ） 強調係数 Ｓproc 処理済画像信号
Claims (10)
- 【請求項１】 第１のマスクを用いて、画像を表すオリ
ジナル画像信号Ｓorg の非鮮鋭マスク信号Ｓusを求め、 第２のマスクを用いて、該第２のマスク内に含まれる前
記オリジナル画像信号Ｓorg のメディアン値Ｓmed を求
め、 前記非鮮鋭マスク信号Ｓusと前記オリジナル画像信号Ｓ
org との差信号Ｓspに対して、前記メディアン値Ｓmed
と前記オリジナル画像信号Ｓorg との差信号Ｓsub に応
じて変化する強調係数β（Ｓsub ）を用いた下記式
（１）にしたがった演算を施すことを特徴とする画像処
理方法。 Ｓproc＝Ｓorg ＋β（Ｓsub ）×Ｓsp （１） - 【請求項２】 前記第２のマスクのサイズを、前記第１
のマスクのサイズよりも大きく設定することを特徴とす
る請求項１記載の画像処理方法。 - 【請求項３】 前記強調係数β（Ｓsub ）は、前記差信
号Ｓsub の絶対値｜Ｓsub ｜が所定の範囲内において
は、該絶対値｜Ｓsub ｜が増大するにしたがって単調に
増大する関数として設定されたことを特徴とする請求項
１または２記載の画像処理方法。 - 【請求項４】 互いに異なる周波数帯域にそれぞれ対応
した互いに大きさの異なるメディアン値算出用の複数の
マスクを用いて、画像を表すオリジナル画像信号Ｓorg
についての複数のメディアン値Ｓmed n （ｎ＝１，２，
…）を求め、 前記オリジナル画像信号Ｓorg を、前記複数のマスクの
大きさにそれぞれ対応した互いに異なる複数の周波数成
分Ｓn に分割し、 前記複数の周波数成分Ｓn に対して、該周波数成分に対
応する前記複数のメディアン値Ｓmed n と前記オリジナ
ル画像信号Ｓorg との差信号Ｓsub n に応じて変化する
複数の強調係数βn （Ｓsub n ）を用いた下記式（２）
または（３）にしたがった演算を施すことを特徴とする
画像処理方法。 Ｓproc＝Σ｛Ｓn ×βn （Ｓsub n ）｝ （２） Ｓproc＝Ｓorg ＋Σ｛Ｓn ×βn （Ｓsub n ）｝ （３） - 【請求項５】 前記強調係数βn （Ｓsub n ）は、前記
差信号Ｓsub n の絶対値｜Ｓsub n ｜が所定の範囲内に
おいては、該絶対値｜Ｓsub n ｜が増大するにしたがっ
て単調に増大する関数として設定されたことを特徴とす
る請求項４記載の画像処理方法。 - 【請求項６】 第１のマスクを用いて、画像を表すオリ
ジナル画像信号Ｓorg の非鮮鋭マスク信号Ｓusを求める
非鮮鋭マスク信号演算手段と、 第２のマスクを用いて、該第２のマスク内に含まれる前
記オリジナル画像信号Ｓorg のメディアン値Ｓmed を求
めるメディアン値演算手段と、 前記メディアン値Ｓmed と前記オリジナル画像信号Ｓor
g との差信号Ｓsub を該差信号に対応付された強調係数
β（Ｓsub ）に変換する変換テーブルを記憶してなる記
憶手段と、 前記非鮮鋭マスク信号Ｓusと前記オリジナル画像信号Ｓ
org との差信号Ｓspに対して、前記強調係数β（Ｓsub
）を用いた下記式（１）にしたがった強調処理を施す
強調手段とを備えてなることを特徴とする画像処理装
置。 Ｓproc＝Ｓorg ＋β（Ｓsub ）×Ｓsp （１） - 【請求項７】 前記第２のマスクのサイズは、前記第１
のマスクのサイズよりも大きく設定されていることを特
徴とする請求項６記載の画像処理装置。 - 【請求項８】 前記記憶手段に記憶された前記変換テー
ブルは、前記差信号Ｓsub の絶対値｜Ｓsub ｜が所定の
範囲内においては、該絶対値｜Ｓsub ｜が増大するにし
たがって前記強調係数β（Ｓsub ）が単調に増大するよ
うに設定されたものであることを特徴とする請求項６ま
たは７記載の画像処理装置。 - 【請求項９】 互いに異なる周波数帯域にそれぞれ対応
した互いに大きさの異なるメディアン値算出用の複数の
マスクを用いて、画像を表すオリジナル画像信号Ｓorg
についての複数のメディアン値Ｓmed n （ｎ＝１，２，
…）を求めるメディアン値演算手段と、 前記オリジナル画像信号Ｓorg を、前記複数のマスクの
大きさにそれぞれ対応した互いに異なる複数の周波数成
分Ｓn に分割する周波数帯域分割手段と、 前記複数のメディアン値Ｓmed n と前記オリジナル画像
信号Ｓorg との差信号Ｓsub n を該差信号にそれぞれ対
応付された強調係数βn （Ｓsub n ）に変換する１以上
の変換テーブルを記憶してなる記憶手段と、 前記複数の周波数成分Ｓn に対して、該周波数成分に対
応する前記強調係数βn （Ｓsub n ）を用いた下記式
（２）または（３）にしたがった強調処理を施す強調処
理手段とを備えてなることを特徴とする画像処理装置。 Ｓproc＝Σ｛Ｓn ×βn （Ｓsub n ）｝ （２） Ｓproc＝Ｓorg ＋Σ｛Ｓn ×βn （Ｓsub n ）｝ （３） - 【請求項１０】 前記記憶手段に記憶された前記変換テ
ーブルは、前記差信号Ｓsub n の絶対値｜Ｓsub n ｜が
所定の範囲内においては、該絶対値｜Ｓsubn ｜が増大
するにしたがって前記強調係数βn （Ｓsub n ）が単調
に増大するように設定されたものであることを特徴とす
る請求項９記載の画像処理装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7250854A JPH0991420A (ja) | 1995-09-28 | 1995-09-28 | 画像処理方法および装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7250854A JPH0991420A (ja) | 1995-09-28 | 1995-09-28 | 画像処理方法および装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH0991420A true JPH0991420A (ja) | 1997-04-04 |
Family
ID=17214001
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP7250854A Withdrawn JPH0991420A (ja) | 1995-09-28 | 1995-09-28 | 画像処理方法および装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0991420A (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1075364A (ja) * | 1995-09-29 | 1998-03-17 | Fuji Photo Film Co Ltd | 画像のダイナミックレンジ圧縮処理方法および装置 |
| JP2006209578A (ja) * | 2005-01-31 | 2006-08-10 | Konica Minolta Holdings Inc | 画像処理方法および画像処理装置ならびに画像処理プログラム |
| JP2007274532A (ja) * | 2006-03-31 | 2007-10-18 | Sony Corp | 撮像装置、高周波成分検出回路、高周波成分検出方法及びコンピュータプログラム |
| JP2012048756A (ja) * | 2011-12-05 | 2012-03-08 | Canon Inc | 画像処理装置及びその制御方法、並びに、コンピュータプログラム及びコンピュータ可読記憶媒体 |
| US8385674B2 (en) | 2006-07-28 | 2013-02-26 | Canon Kabushiki Kaisha | Image processing apparatus and control method thereof |
-
1995
- 1995-09-28 JP JP7250854A patent/JPH0991420A/ja not_active Withdrawn
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1075364A (ja) * | 1995-09-29 | 1998-03-17 | Fuji Photo Film Co Ltd | 画像のダイナミックレンジ圧縮処理方法および装置 |
| JP2006209578A (ja) * | 2005-01-31 | 2006-08-10 | Konica Minolta Holdings Inc | 画像処理方法および画像処理装置ならびに画像処理プログラム |
| JP4661238B2 (ja) * | 2005-01-31 | 2011-03-30 | コニカミノルタホールディングス株式会社 | 画像処理方法および画像処理装置ならびに画像処理プログラム |
| JP2007274532A (ja) * | 2006-03-31 | 2007-10-18 | Sony Corp | 撮像装置、高周波成分検出回路、高周波成分検出方法及びコンピュータプログラム |
| US7889272B2 (en) | 2006-03-31 | 2011-02-15 | Sony Corporation | Imaging apparatus, high frequency component detection circuit, high frequency component detection method and computer program product |
| US8385674B2 (en) | 2006-07-28 | 2013-02-26 | Canon Kabushiki Kaisha | Image processing apparatus and control method thereof |
| JP2012048756A (ja) * | 2011-12-05 | 2012-03-08 | Canon Inc | 画像処理装置及びその制御方法、並びに、コンピュータプログラム及びコンピュータ可読記憶媒体 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3690882B2 (ja) | 画像の強調処理方法および装置 | |
| JP3816151B2 (ja) | 画像処理方法および装置 | |
| US20110125030A1 (en) | Medical diagnostic device and method of improving image quality of medical diagnostic device | |
| Suzuki et al. | Neural filter with selection of input features and its application to image quality improvement of medical image sequences | |
| EP1952344A4 (en) | METHOD AND SYSTEM FOR ENHANCING DIGITAL IMAGES | |
| JP3669455B2 (ja) | 画像処理方法および画像処理装置 | |
| JPH10171983A (ja) | 画像処理方法および装置 | |
| JP3738788B2 (ja) | 画像のダイナミックレンジ圧縮処理方法および装置 | |
| Stahl et al. | Noise-resistant weak-structure enhancement for digital radiography | |
| JP3700798B2 (ja) | 画像処理方法および装置 | |
| JP3738791B2 (ja) | 画像処理方法および装置 | |
| Mustafa et al. | A systematic review: Contrast enhancement based on spatial and frequency domain | |
| JP3932489B2 (ja) | 画像処理方法および装置 | |
| JPH0991420A (ja) | 画像処理方法および装置 | |
| Agaian et al. | Novel medical image enhancement algorithms | |
| JPH0991421A (ja) | 画像処理方法および装置 | |
| JPH09107479A (ja) | 画像処理方法および装置 | |
| US20230133074A1 (en) | Method and Apparatus for Noise Reduction | |
| JPH10105701A (ja) | 放射線画像強調処理方法および装置 | |
| JP3731762B2 (ja) | 画像処理方法および装置 | |
| JP3690844B2 (ja) | 画像処理方法および装置 | |
| JP2003337942A (ja) | 画像処理装置 | |
| JPH0991423A (ja) | 画像処理方法および装置 | |
| JPH09163227A (ja) | 画像処理方法および装置 | |
| JP3731761B2 (ja) | 画像処理方法および装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20021203 |