JPH1075395A - 画像処理方法および装置 - Google Patents
画像処理方法および装置Info
- Publication number
- JPH1075395A JPH1075395A JP8182155A JP18215596A JPH1075395A JP H1075395 A JPH1075395 A JP H1075395A JP 8182155 A JP8182155 A JP 8182155A JP 18215596 A JP18215596 A JP 18215596A JP H1075395 A JPH1075395 A JP H1075395A
- Authority
- JP
- Japan
- Prior art keywords
- image signal
- band
- sorg
- signal
- limited
- 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
- 238000003672 processing method Methods 0.000 title claims description 33
- 238000006243 chemical reaction Methods 0.000 claims abstract description 84
- 230000004044 response Effects 0.000 claims abstract description 29
- 101000920618 Homo sapiens Transcription and mRNA export factor ENY2 Proteins 0.000 claims description 88
- 102100031954 Transcription and mRNA export factor ENY2 Human genes 0.000 claims description 88
- 238000000034 method Methods 0.000 claims description 87
- 230000008569 process Effects 0.000 claims description 63
- 230000001629 suppression Effects 0.000 claims description 26
- 239000002131 composite material Substances 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000003325 tomography Methods 0.000 claims description 8
- 238000012886 linear function Methods 0.000 claims description 6
- 230000002708 enhancing effect Effects 0.000 claims description 2
- 238000001914 filtration Methods 0.000 abstract description 44
- 230000010354 integration Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 51
- 230000000694 effects Effects 0.000 description 17
- 230000008859 change Effects 0.000 description 10
- 238000003384 imaging method Methods 0.000 description 9
- 230000001965 increasing effect Effects 0.000 description 8
- 230000000630 rising effect Effects 0.000 description 8
- 210000002784 stomach Anatomy 0.000 description 8
- 230000014509 gene expression Effects 0.000 description 7
- 210000000988 bone and bone Anatomy 0.000 description 6
- 230000005855 radiation Effects 0.000 description 4
- 210000004072 lung Anatomy 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000009607 mammography Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000002542 deteriorative effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
- G06T5/75—Unsharp masking
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Picture Signal Circuits (AREA)
- Transforming Light Signals Into Electric Signals (AREA)
Abstract
(57)【要約】
【課題】 ボケマスク処理を行う画像処理方法および装
置において、エッジ部近傍に発生するアーチファクトを
防止する。 【解決手段】 ボケ画像信号作成手段１において、入力
された原画像信号Ｓorg に基づいて互いに周波数応答特
性が異なる多重解像度のボケ画像信号Ｓusk （ｋ＝１〜
ｎ）を作成し、帯域制限画像信号作成手段２においてボ
ケ画像信号Ｓusk に基づいて複数の帯域制限画像信号を
作成する。変換手段３において、帯域制限画像信号のう
ち少なくとも１つの信号について、少なくとも一部を小
さくするような変換処理を施し、積算手段４において、
変換された帯域制限画像信号を積算して積算信号を作成
する。そして、周波数強調処理手段５において、積算信
号に所定の強調係数を乗じて原画像信号Ｓorg と加算す
ることにより、原画像の高周波成分が強調された処理済
画像Ｓprocを得る。
置において、エッジ部近傍に発生するアーチファクトを
防止する。 【解決手段】 ボケ画像信号作成手段１において、入力
された原画像信号Ｓorg に基づいて互いに周波数応答特
性が異なる多重解像度のボケ画像信号Ｓusk （ｋ＝１〜
ｎ）を作成し、帯域制限画像信号作成手段２においてボ
ケ画像信号Ｓusk に基づいて複数の帯域制限画像信号を
作成する。変換手段３において、帯域制限画像信号のう
ち少なくとも１つの信号について、少なくとも一部を小
さくするような変換処理を施し、積算手段４において、
変換された帯域制限画像信号を積算して積算信号を作成
する。そして、周波数強調処理手段５において、積算信
号に所定の強調係数を乗じて原画像信号Ｓorg と加算す
ることにより、原画像の高周波成分が強調された処理済
画像Ｓprocを得る。
Description
【０００１】
【発明の属する技術分野】本発明は画像信号に対して、
この画像信号の高周波成分を強調する非鮮鋭マスク処理
（いわゆるボケマスク処理）を行う画像処理方法および
装置に関するものである。
この画像信号の高周波成分を強調する非鮮鋭マスク処理
（いわゆるボケマスク処理）を行う画像処理方法および
装置に関するものである。
【０００２】
【従来の技術】画像を表す画像信号を得、この画像信号
に適切な画像処理を施した後、画像を再生表示すること
が種々の分野で行われている。例えば、放射線画像の診
断性能を向上させるために、画像信号に対して非鮮鋭マ
スク（以下ボケマスクという）を用いて周波数強調処理
を施す方法が本出願人により提案されている（特開昭55
-163472 号公報、特開昭55-87953号公報等）。この周波
数強調処理は、読み取った原画像信号Ｓorg に、この原
画像信号Ｓorg から非鮮鋭マスク画像信号（以下ボケ画
像信号という）Ｓusを引いたものに強調係数βを乗じた
ものを加える処理を施すもので、これにより、画像にお
いて所定の空間周波数成分を強調するようにしたもので
ある。これを式で表すと下記の式（１）のようになる。
に適切な画像処理を施した後、画像を再生表示すること
が種々の分野で行われている。例えば、放射線画像の診
断性能を向上させるために、画像信号に対して非鮮鋭マ
スク（以下ボケマスクという）を用いて周波数強調処理
を施す方法が本出願人により提案されている（特開昭55
-163472 号公報、特開昭55-87953号公報等）。この周波
数強調処理は、読み取った原画像信号Ｓorg に、この原
画像信号Ｓorg から非鮮鋭マスク画像信号（以下ボケ画
像信号という）Ｓusを引いたものに強調係数βを乗じた
ものを加える処理を施すもので、これにより、画像にお
いて所定の空間周波数成分を強調するようにしたもので
ある。これを式で表すと下記の式（１）のようになる。
【０００３】 Ｓproc＝Ｓorg ＋β×（Ｓorg −Ｓus） …（１） （Ｓproc：周波数強調処理された信号、Ｓorg ：原画像
信号、Ｓus：ボケ画像信号、β：強調係数） ここで、ボケ画像信号Ｓusは、例えば画像を構成する画
素の１画素おきに、各画素の周囲のＭ×Ｎの範囲内の原
画像信号Ｓorg について、 Ｓus＝ΣＳorg ／（Ｍ×Ｎ） …（２） なる演算を行うことにより求められる。
信号、Ｓus：ボケ画像信号、β：強調係数） ここで、ボケ画像信号Ｓusは、例えば画像を構成する画
素の１画素おきに、各画素の周囲のＭ×Ｎの範囲内の原
画像信号Ｓorg について、 Ｓus＝ΣＳorg ／（Ｍ×Ｎ） …（２） なる演算を行うことにより求められる。
【０００４】また、いわゆる断層撮影法（たとえば、特
開昭58-67245号公報参照）により得られた放射線画像
（断層撮影画像）中の観察の対象となる中央領域におい
て、撮影対象とされた断層面以外に存在する放射線透過
量の大きく変化した部分の像が記録シートの移動方向に
沿って生ずる障害陰影（以下これを流れ像と称する）を
除去するための画像処理方法も行われている（たとえば
特開平3-276265号等）。この方法は、断層撮影により得
られた放射線画像の画像信号のうち、流れ像に対応する
低空間周波数成分を除去することにより、流れ像が除去
された画像を生成するようにしたものである。
開昭58-67245号公報参照）により得られた放射線画像
（断層撮影画像）中の観察の対象となる中央領域におい
て、撮影対象とされた断層面以外に存在する放射線透過
量の大きく変化した部分の像が記録シートの移動方向に
沿って生ずる障害陰影（以下これを流れ像と称する）を
除去するための画像処理方法も行われている（たとえば
特開平3-276265号等）。この方法は、断層撮影により得
られた放射線画像の画像信号のうち、流れ像に対応する
低空間周波数成分を除去することにより、流れ像が除去
された画像を生成するようにしたものである。
【０００５】
【発明が解決しようとする課題】上記周波数強調処理で
は、画像の濃度が急激に変化しているエッジ部の近傍に
おいてボケ画像信号を作成する場合、図42に示すように
Ｍ×Ｎのマスク内にエッジ部が含まれてしまうため、そ
のボケ画像信号がエッジ部の濃度に影響されてしまい、
このボケ画像信号により処理が行われた画像に、オーバ
ーシュート、アンダーシュートあるいは偽輪郭などのア
ーチファクトが発生し、これにより画質が低下してしま
うという問題がある。これは、上述した断層撮影により
画像に発生する流れ像を除去するためにボケマスク処理
を行う場合にも同様に発生する問題である。
は、画像の濃度が急激に変化しているエッジ部の近傍に
おいてボケ画像信号を作成する場合、図42に示すように
Ｍ×Ｎのマスク内にエッジ部が含まれてしまうため、そ
のボケ画像信号がエッジ部の濃度に影響されてしまい、
このボケ画像信号により処理が行われた画像に、オーバ
ーシュート、アンダーシュートあるいは偽輪郭などのア
ーチファクトが発生し、これにより画質が低下してしま
うという問題がある。これは、上述した断層撮影により
画像に発生する流れ像を除去するためにボケマスク処理
を行う場合にも同様に発生する問題である。
【０００６】本発明は上記事情に鑑み、エッジ部近傍に
アーチファクトが発生することなく、良好な処理済画像
を得るボケマスク処理を行うことができる画像処理方法
および装置を提供することを目的とするものである。
アーチファクトが発生することなく、良好な処理済画像
を得るボケマスク処理を行うことができる画像処理方法
および装置を提供することを目的とするものである。
【０００７】
【課題を解決するための手段】本発明の画像処理方法お
よび装置は、原画像を表す原画像信号に、該原画像の高
周波成分に関する信号を加算することにより、該原画像
の高周波成分を強調する画像処理方法において、前記原
画像信号に基づいて、互いに周波数応答特性が異なる複
数の非鮮鋭マスク画像信号を作成し、前記原画像信号お
よび前記複数の非鮮鋭マスク画像信号、または前記複数
の非鮮鋭マスク画像信号に基づいて、前記原画像信号の
複数の周波数帯域ごとの信号を表す複数の帯域制限画像
信号を作成し、該各帯域制限画像信号のうち少なくとも
１つの信号について、該帯域制限画像信号の少なくとも
一部を小さくするような変換処理を施して複数の変換画
像信号を作成し、該各変換画像信号を積算することによ
り前記原画像信号に加算する前記高周波成分に関する信
号を得ることを特徴とするものである。
よび装置は、原画像を表す原画像信号に、該原画像の高
周波成分に関する信号を加算することにより、該原画像
の高周波成分を強調する画像処理方法において、前記原
画像信号に基づいて、互いに周波数応答特性が異なる複
数の非鮮鋭マスク画像信号を作成し、前記原画像信号お
よび前記複数の非鮮鋭マスク画像信号、または前記複数
の非鮮鋭マスク画像信号に基づいて、前記原画像信号の
複数の周波数帯域ごとの信号を表す複数の帯域制限画像
信号を作成し、該各帯域制限画像信号のうち少なくとも
１つの信号について、該帯域制限画像信号の少なくとも
一部を小さくするような変換処理を施して複数の変換画
像信号を作成し、該各変換画像信号を積算することによ
り前記原画像信号に加算する前記高周波成分に関する信
号を得ることを特徴とするものである。
【０００８】前記複数の変換画像信号の作成は、前記帯
域制限画像信号の絶対値が所定の閾値より大きい帯域制
限画像信号について、該帯域制限画像信号の絶対値を小
さくするように、該各帯域制限画像信号を変換すること
により行うことが望ましい。この際、前記各帯域制限画
像信号が所定の閾値より大きいほど該帯域制限画像信号
の絶対値を小さくするとともに、前記各帯域制限画像信
号の絶対値が該所定の閾値よりも小さい他の閾値よりも
小さいほど該帯域制限画像信号の絶対値を小さくするよ
うに、該各帯域制限画像信号を変換することにより行っ
てもよい。また、前記各帯域制限画像信号の周波数帯域
に応じて該各帯域制限画像信号の絶対値の大きさを変化
させるように、該各帯域制限画像信号を変換することに
より行ってもよい。
域制限画像信号の絶対値が所定の閾値より大きい帯域制
限画像信号について、該帯域制限画像信号の絶対値を小
さくするように、該各帯域制限画像信号を変換すること
により行うことが望ましい。この際、前記各帯域制限画
像信号が所定の閾値より大きいほど該帯域制限画像信号
の絶対値を小さくするとともに、前記各帯域制限画像信
号の絶対値が該所定の閾値よりも小さい他の閾値よりも
小さいほど該帯域制限画像信号の絶対値を小さくするよ
うに、該各帯域制限画像信号を変換することにより行っ
てもよい。また、前記各帯域制限画像信号の周波数帯域
に応じて該各帯域制限画像信号の絶対値の大きさを変化
させるように、該各帯域制限画像信号を変換することに
より行ってもよい。
【０００９】また、前記複数の変換画像信号の作成は、
前記各帯域制限画像信号を、前記周波数帯域ごとに異な
る複数の関数に基づいて、該帯域制限画像信号の絶対値
に基づいて決まる該絶対値以下の値となるように変換す
ることにより行ってもよい。この際、前記変換画像信号
の作成に使用される各関数は、前記帯域制限画像信号の
絶対値が所定の値より大きい帯域制限画像信号について
は、該変換画像信号の値がほぼ一定の値となるように該
帯域制限画像信号を変換する関数であり、該関数が高周
波帯域を処理する関数であるほど、前記所定の値が小さ
い値であることが望ましい。さらに、前記変換画像信号
の作成に使用される各関数は、該関数が低周波帯域を処
理する関数であるほど、前記帯域制限画像信号の絶対値
が０近傍の所定の範囲内の値である帯域制限画像信号を
変換した際に得られる該変換画像信号の絶対値が小さい
値であることが望ましい。
前記各帯域制限画像信号を、前記周波数帯域ごとに異な
る複数の関数に基づいて、該帯域制限画像信号の絶対値
に基づいて決まる該絶対値以下の値となるように変換す
ることにより行ってもよい。この際、前記変換画像信号
の作成に使用される各関数は、前記帯域制限画像信号の
絶対値が所定の値より大きい帯域制限画像信号について
は、該変換画像信号の値がほぼ一定の値となるように該
帯域制限画像信号を変換する関数であり、該関数が高周
波帯域を処理する関数であるほど、前記所定の値が小さ
い値であることが望ましい。さらに、前記変換画像信号
の作成に使用される各関数は、該関数が低周波帯域を処
理する関数であるほど、前記帯域制限画像信号の絶対値
が０近傍の所定の範囲内の値である帯域制限画像信号を
変換した際に得られる該変換画像信号の絶対値が小さい
値であることが望ましい。
【００１０】なお、前記帯域制限画像信号の作成、前記
変換画像信号の作成、前記高周波成分に関する信号の作
成、および該高周波成分に関する信号の前記原画像信号
への加算は、具体的には、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)＋ｆ2(Ｓus1 −Ｓus2)＋… ＋ｆk(Ｓusk-1−Ｓusk )＋…＋ｆN(ＳusN-1−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことが望ましい。
変換画像信号の作成、前記高周波成分に関する信号の作
成、および該高周波成分に関する信号の前記原画像信号
への加算は、具体的には、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)＋ｆ2(Ｓus1 −Ｓus2)＋… ＋ｆk(Ｓusk-1−Ｓusk )＋…＋ｆN(ＳusN-1−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことが望ましい。
【００１１】あるいは、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN ） ＝（1/N)・｛ｆ1（Ｓorg −Ｓus1)＋ｆ2（Ｓorg −Ｓus2)＋… ＋fk(Ｓorg −Ｓusk)＋…＋ｆN（Ｓorg −ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行ってもよい。
数） にしたがって行ってもよい。
【００１２】また、前記複数の変換画像信号を、前記帯
域制限画像信号の絶対値が所定の閾値より大きい帯域制
限画像信号について、該帯域制限画像信号の絶対値を小
さくするように、該各帯域制限画像信号を変換すること
により作成する場合には、前記帯域制限画像信号の作
成、前記変換画像信号の作成、前記高周波成分に関する
信号の作成、および該高周波成分に関する信号の前記原
画像信号への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓus1,Ｓus2,…ＳusN ） ＝｛ｆ2（Ｓus1 −Ｓus2)＋ｆ3（Ｓus2 −Ｓus3)＋… ＋fk(Ｓusk-1 −Ｓusk)＋…＋ｆN（ＳusN-1 −ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=2〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行ってもよい。
域制限画像信号の絶対値が所定の閾値より大きい帯域制
限画像信号について、該帯域制限画像信号の絶対値を小
さくするように、該各帯域制限画像信号を変換すること
により作成する場合には、前記帯域制限画像信号の作
成、前記変換画像信号の作成、前記高周波成分に関する
信号の作成、および該高周波成分に関する信号の前記原
画像信号への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓus1,Ｓus2,…ＳusN ） ＝｛ｆ2（Ｓus1 −Ｓus2)＋ｆ3（Ｓus2 −Ｓus3)＋… ＋fk(Ｓusk-1 −Ｓusk)＋…＋ｆN（ＳusN-1 −ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=2〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行ってもよい。
【００１３】あるいは、同様に、前記複数の変換画像信
号を、前記帯域制限画像信号の絶対値が所定の閾値より
大きい帯域制限画像信号について、該帯域制限画像信号
の絶対値を小さくするように、該各帯域制限画像信号を
変換することにより作成する場合には、上記処理を下記
の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓus1,Ｓus2,…ＳusN ） ＝（1/N)・｛ｆ2（Ｓus1 −Ｓus2)＋ｆ3（Ｓus1 −Ｓus3)＋… ＋fk(Ｓus1 −Ｓusk)＋…＋ｆN（Ｓus1 −ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=2〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行ってもよい。
号を、前記帯域制限画像信号の絶対値が所定の閾値より
大きい帯域制限画像信号について、該帯域制限画像信号
の絶対値を小さくするように、該各帯域制限画像信号を
変換することにより作成する場合には、上記処理を下記
の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓus1,Ｓus2,…ＳusN ） ＝（1/N)・｛ｆ2（Ｓus1 −Ｓus2)＋ｆ3（Ｓus1 −Ｓus3)＋… ＋fk(Ｓus1 −Ｓusk)＋…＋ｆN（Ｓus1 −ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=2〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行ってもよい。
【００１４】さらには、上記の変換画像信号の作成は、
全て、前記各帯域制限画像信号に対して、前記強調係数
に応じた変換処理を施すことにより行ってもよい。この
際、この変換処理は、前記各周波数帯域ごとに異なる変
換処理であって、前記帯域制限画像信号の周波数帯域が
低周波帯域であるほど、該帯域制限画像信号を強く抑制
するような変換処理であり、前記強調係数が大きいとき
ほど、高周波帯域の帯域制限画像信号に対する抑制の度
合いと、低周波帯域の帯域制限画像信号に対する抑制の
度合いの差が大きいことが望ましい。
全て、前記各帯域制限画像信号に対して、前記強調係数
に応じた変換処理を施すことにより行ってもよい。この
際、この変換処理は、前記各周波数帯域ごとに異なる変
換処理であって、前記帯域制限画像信号の周波数帯域が
低周波帯域であるほど、該帯域制限画像信号を強く抑制
するような変換処理であり、前記強調係数が大きいとき
ほど、高周波帯域の帯域制限画像信号に対する抑制の度
合いと、低周波帯域の帯域制限画像信号に対する抑制の
度合いの差が大きいことが望ましい。
【００１５】さらに、前記複数の変換画像信号の作成
は、前記帯域制限画像信号を、該帯域制限画像信号の絶
対値に基づいて決まる該絶対値以下の値となるように変
換することにより抑制画像信号を作成し、前記原画像信
号および前記複数の非鮮鋭マスク画像信号に基づいて、
前記抑制画像信号の作成に使用された帯域制限画像信号
よりも低い周波数帯域の信号を含む補助画像信号を作成
し、該補助画像信号を、該補助画像信号の絶対値が小さ
いほど１に近く、大きいほど０に近い値となるように変
換することにより前記抑制画像信号にそれぞれ対応する
倍率信号を作成し、前記抑制画像信号に、該抑制画像信
号に対応する前記倍率信号を乗ずることにより行っても
よい。
は、前記帯域制限画像信号を、該帯域制限画像信号の絶
対値に基づいて決まる該絶対値以下の値となるように変
換することにより抑制画像信号を作成し、前記原画像信
号および前記複数の非鮮鋭マスク画像信号に基づいて、
前記抑制画像信号の作成に使用された帯域制限画像信号
よりも低い周波数帯域の信号を含む補助画像信号を作成
し、該補助画像信号を、該補助画像信号の絶対値が小さ
いほど１に近く、大きいほど０に近い値となるように変
換することにより前記抑制画像信号にそれぞれ対応する
倍率信号を作成し、前記抑制画像信号に、該抑制画像信
号に対応する前記倍率信号を乗ずることにより行っても
よい。
【００１６】これは、具体的には、前記帯域制限画像信
号の作成、前記変換画像信号の作成、前記積算信号の作
成、および該積算信号の前記原画像信号への加算を、下
記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)・ｇ（Ｓus1−Ｓus2） ＋ｆ2(Ｓus1 −Ｓus2)・ｇ（Ｓus2−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk )・ｇ（Ｓusk−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（ＳusN−ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換して前記
抑制信号を作成するための関数 ｇ：前記各補助画像信号を変換して前記倍率信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことが望ましい。
号の作成、前記変換画像信号の作成、前記積算信号の作
成、および該積算信号の前記原画像信号への加算を、下
記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)・ｇ（Ｓus1−Ｓus2） ＋ｆ2(Ｓus1 −Ｓus2)・ｇ（Ｓus2−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk )・ｇ（Ｓusk−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（ＳusN−ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換して前記
抑制信号を作成するための関数 ｇ：前記各補助画像信号を変換して前記倍率信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことが望ましい。
【００１７】あるいは、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)・ｇ（Ｓorg−Ｓus2） ＋ｆ2(Ｓus1 −Ｓus2)・ｇ（Ｓorg−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk )・ｇ（Ｓorg−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（Ｓorg−ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換して前記
抑制信号を作成するための関数 ｇ：前記各補助画像信号を変換して前記倍率信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行ってもよい。
抑制信号を作成するための関数 ｇ：前記各補助画像信号を変換して前記倍率信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行ってもよい。
【００１８】さらに、前記複数の変換画像信号の作成
は、変換される帯域制限画像信号である被変換帯域制限
画像信号よりも低い周波数帯域の帯域制限画像信号であ
る低周波側帯域制限画像信号を、原点をとおり該原点に
おける傾きがほぼ０で、処理される値が大きくなるにつ
れて該傾きが漸増する非線形関数に基づいて変換するこ
とにより前記被変換帯域制限画像信号の補助画像信号を
作成し、該補助画像信号を前記被変換帯域制限画像信号
に加算することにより複合帯域制限画像信号を作成し、
該複合帯域制限画像信号を、該複合帯域制限画像信号の
絶対値に基づいて決まる該絶対値以下の値となるように
変換することにより行ってもよい。
は、変換される帯域制限画像信号である被変換帯域制限
画像信号よりも低い周波数帯域の帯域制限画像信号であ
る低周波側帯域制限画像信号を、原点をとおり該原点に
おける傾きがほぼ０で、処理される値が大きくなるにつ
れて該傾きが漸増する非線形関数に基づいて変換するこ
とにより前記被変換帯域制限画像信号の補助画像信号を
作成し、該補助画像信号を前記被変換帯域制限画像信号
に加算することにより複合帯域制限画像信号を作成し、
該複合帯域制限画像信号を、該複合帯域制限画像信号の
絶対値に基づいて決まる該絶対値以下の値となるように
変換することにより行ってもよい。
【００１９】具体的には、前記帯域制限画像信号の作
成、前記変換画像信号の作成、前記積算信号の作成、お
よび該積算信号の前記原画像信号への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝［ｆ1｛（Ｓorg −Ｓus1)＋ｇ（Ｓus1−Ｓus2）｝ ＋ｆ2｛（Ｓus1 −Ｓus2)＋ｇ（Ｓus2−Ｓus3）｝＋… ＋ｆk｛（Ｓusk-1−Ｓusk )＋ｇ（Ｓusk−Ｓusk+1）｝＋… ＋ｆN｛（ＳusN-1−ＳusN)＋ｇ（ＳusN−ＳusN+1）｝］ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各複合帯域制限画像信号を変換する
ための関数 ｇ：前記各帯域制限画像信号を変換して前記補助画像信
号を作成するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことが望ましい。
成、前記変換画像信号の作成、前記積算信号の作成、お
よび該積算信号の前記原画像信号への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝［ｆ1｛（Ｓorg −Ｓus1)＋ｇ（Ｓus1−Ｓus2）｝ ＋ｆ2｛（Ｓus1 −Ｓus2)＋ｇ（Ｓus2−Ｓus3）｝＋… ＋ｆk｛（Ｓusk-1−Ｓusk )＋ｇ（Ｓusk−Ｓusk+1）｝＋… ＋ｆN｛（ＳusN-1−ＳusN)＋ｇ（ＳusN−ＳusN+1）｝］ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各複合帯域制限画像信号を変換する
ための関数 ｇ：前記各帯域制限画像信号を変換して前記補助画像信
号を作成するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことが望ましい。
【００２０】また、上記変換画像信号の作成は、全て、
前記原画像を撮影した際の撮影部位に応じて、前記各帯
域制限画像信号を変換することにより行うことが望まし
い。
前記原画像を撮影した際の撮影部位に応じて、前記各帯
域制限画像信号を変換することにより行うことが望まし
い。
【００２１】さらに、前記原画像が断層撮影により得ら
れた流れ像を含むものである場合には、前記複数の非鮮
鋭マスク画像信号の作成を、前記原画像信号に対して該
流れ像の方向に沿う１次元非鮮鋭マスクによりマスク処
理を施すことにより行うことが望ましい。
れた流れ像を含むものである場合には、前記複数の非鮮
鋭マスク画像信号の作成を、前記原画像信号に対して該
流れ像の方向に沿う１次元非鮮鋭マスクによりマスク処
理を施すことにより行うことが望ましい。
【００２２】
【発明の効果】本発明による画像処理方法および装置
は、互いに周波数応答特性が異なる複数の非鮮鋭マスク
画像信号を作成し、この非鮮鋭マスク画像信号から原画
像信号の複数の周波数帯域ごとの信号を表す複数の帯域
制限画像信号を作成する。この各帯域制限画像信号は、
複数の非鮮鋭マスク画像信号の周波数特性に応じて、原
画像をこの周波数特性に対応した周波数帯域の信号とし
て表すものとなる。ここで、この帯域制限画像信号は、
原画像の濃度変化が比較的小さいいわゆる平坦部におい
ては、各周波数帯域において値は小さくなるものであ
る。これに対して、濃度が急激に変化するエッジ部近傍
においては、帯域制限画像信号が比較的低周波数帯域で
ある場合、すなわち非鮮鋭マスク画像信号を得る際の非
鮮鋭マスクのサイズが比較的大きい場合は、そのエッジ
部近傍の画素については、非鮮鋭マスクにエッジ部が含
まれてしまうため、帯域制限画像信号がエッジ部の影響
を受けて信号値の絶対値が比較的大きなものとなる。こ
のように、本来エッジ部でない部分がエッジ部の濃度値
の影響を受けることにより、画像処理を行うことにより
得られた画像のエッジ部にはオーバーシュート、アンダ
ーシュートなどのアーチファクトが発生してしまう。
は、互いに周波数応答特性が異なる複数の非鮮鋭マスク
画像信号を作成し、この非鮮鋭マスク画像信号から原画
像信号の複数の周波数帯域ごとの信号を表す複数の帯域
制限画像信号を作成する。この各帯域制限画像信号は、
複数の非鮮鋭マスク画像信号の周波数特性に応じて、原
画像をこの周波数特性に対応した周波数帯域の信号とし
て表すものとなる。ここで、この帯域制限画像信号は、
原画像の濃度変化が比較的小さいいわゆる平坦部におい
ては、各周波数帯域において値は小さくなるものであ
る。これに対して、濃度が急激に変化するエッジ部近傍
においては、帯域制限画像信号が比較的低周波数帯域で
ある場合、すなわち非鮮鋭マスク画像信号を得る際の非
鮮鋭マスクのサイズが比較的大きい場合は、そのエッジ
部近傍の画素については、非鮮鋭マスクにエッジ部が含
まれてしまうため、帯域制限画像信号がエッジ部の影響
を受けて信号値の絶対値が比較的大きなものとなる。こ
のように、本来エッジ部でない部分がエッジ部の濃度値
の影響を受けることにより、画像処理を行うことにより
得られた画像のエッジ部にはオーバーシュート、アンダ
ーシュートなどのアーチファクトが発生してしまう。
【００２３】本発明はこの点に鑑みてなされたものであ
り、各帯域制限画像信号のうち少なくとも１つの信号に
ついて、該帯域制限画像信号の少なくとも一部を小さく
するような変換処理を施して、この変換された各帯域制
限画像信号を積算することにより、原画像信号に加算す
るための高周波成分に関する信号を得るようにしたもの
である。このため、比較的信号値の絶対値の大きい帯域
制限画像信号は、原画像信号に加算するための高周波成
分に関する信号に対する影響力が小さくなり実質的の非
鮮鋭マスクのサイズを小さくしたものと同様の信号とな
る。そしてこれにより、濃度が急激に変化するエッジ部
近傍においても、アーチファクトの原因となる信号は影
響力が弱められるため、処理を施すことにより得られる
画像をアーチファクトのない良好なものとすることがで
きる。
り、各帯域制限画像信号のうち少なくとも１つの信号に
ついて、該帯域制限画像信号の少なくとも一部を小さく
するような変換処理を施して、この変換された各帯域制
限画像信号を積算することにより、原画像信号に加算す
るための高周波成分に関する信号を得るようにしたもの
である。このため、比較的信号値の絶対値の大きい帯域
制限画像信号は、原画像信号に加算するための高周波成
分に関する信号に対する影響力が小さくなり実質的の非
鮮鋭マスクのサイズを小さくしたものと同様の信号とな
る。そしてこれにより、濃度が急激に変化するエッジ部
近傍においても、アーチファクトの原因となる信号は影
響力が弱められるため、処理を施すことにより得られる
画像をアーチファクトのない良好なものとすることがで
きる。
【００２４】この際、帯域制限画像信号の変換方法とし
て、帯域制限画像信号の絶対値が所定の閾値よりも大き
い場合に、この絶対値が小さくなるように帯域制限画像
信号を変換することにより、特に大きな信号の影響力を
弱めることができる。あるいは帯域制限画像信号の絶対
値が所定の閾値よりも小さい他の閾値よりも小さいほど
この帯域制限画像信号の絶対値を小さくするように補正
する方法を用いれば、画像中においてノイズと見なせる
ような信号値の絶対値の小さい成分のレスポンスを小さ
くすることができ、これにより得られる画像のノイズを
低減することができる。また、帯域制限画像信号の周波
数帯域に応じて帯域制限画像信号の絶対値を変化させる
ように変換することにより、周波数帯域に応じた強調処
理を行うことも可能となる。
て、帯域制限画像信号の絶対値が所定の閾値よりも大き
い場合に、この絶対値が小さくなるように帯域制限画像
信号を変換することにより、特に大きな信号の影響力を
弱めることができる。あるいは帯域制限画像信号の絶対
値が所定の閾値よりも小さい他の閾値よりも小さいほど
この帯域制限画像信号の絶対値を小さくするように補正
する方法を用いれば、画像中においてノイズと見なせる
ような信号値の絶対値の小さい成分のレスポンスを小さ
くすることができ、これにより得られる画像のノイズを
低減することができる。また、帯域制限画像信号の周波
数帯域に応じて帯域制限画像信号の絶対値を変化させる
ように変換することにより、周波数帯域に応じた強調処
理を行うことも可能となる。
【００２５】また、閾値よりも大きいか否かによらず、
周波数帯域に応じて異なる非線形関数によって変換処理
を施せば、より適切な処理を施すことが可能となり、こ
れにより、画像処理により生成される処理済画像信号の
全体の周波数特性を自在にコントロールすることもでき
る。これは、上記アーチファクトのみならず、周波数帯
域の境界に発生する階段状のアーチファクトをも抑制す
るという効果を有する。
周波数帯域に応じて異なる非線形関数によって変換処理
を施せば、より適切な処理を施すことが可能となり、こ
れにより、画像処理により生成される処理済画像信号の
全体の周波数特性を自在にコントロールすることもでき
る。これは、上記アーチファクトのみならず、周波数帯
域の境界に発生する階段状のアーチファクトをも抑制す
るという効果を有する。
【００２６】あるいは、変換する帯域制限画像信号より
も低い周波数帯域の画像信号を用いて、２種類の信号を
加算あるいは乗算するなどして変換画像信号を作成する
ことにより、上記階段状のアーチファクトをさらに抑制
し、より滑らかな画像信号を作成して良好な処理済画像
を得ることができる。さらに、原画像を得た際の撮影部
位に応じてこのような変換を行うことにより、各撮影部
位に適した高周波成分の強調を行うことができる。
も低い周波数帯域の画像信号を用いて、２種類の信号を
加算あるいは乗算するなどして変換画像信号を作成する
ことにより、上記階段状のアーチファクトをさらに抑制
し、より滑らかな画像信号を作成して良好な処理済画像
を得ることができる。さらに、原画像を得た際の撮影部
位に応じてこのような変換を行うことにより、各撮影部
位に適した高周波成分の強調を行うことができる。
【００２７】また、上述したように断層撮影を行うこと
により得られる流れ像を除去するために、非鮮鋭マスク
処理を行う方法が提案されているが（特開平3-276265
号）、この場合においても、急激に濃度が変化するエッ
ジ部の近傍においてアーチファクトが発生するという問
題がある。この場合、流れ像を除去するためには非鮮鋭
マスクを１次元として非鮮鋭マスク画像信号を得て非鮮
鋭マスク処理を行うものであるが、この非鮮鋭マスク画
像信号を得る際に、互いに周波数応答特性が異なる複数
の非鮮鋭マスク画像信号を得、これに基づいて上述した
ような帯域制限画像信号を得、この帯域制限画像信号の
絶対値が比較的大きい場合は、その絶対値が小さくなる
ように変換してその後の処理を行うことにより、流れ像
を除去するとともにアーチファクトのない良好な画像を
得ることができる。
により得られる流れ像を除去するために、非鮮鋭マスク
処理を行う方法が提案されているが（特開平3-276265
号）、この場合においても、急激に濃度が変化するエッ
ジ部の近傍においてアーチファクトが発生するという問
題がある。この場合、流れ像を除去するためには非鮮鋭
マスクを１次元として非鮮鋭マスク画像信号を得て非鮮
鋭マスク処理を行うものであるが、この非鮮鋭マスク画
像信号を得る際に、互いに周波数応答特性が異なる複数
の非鮮鋭マスク画像信号を得、これに基づいて上述した
ような帯域制限画像信号を得、この帯域制限画像信号の
絶対値が比較的大きい場合は、その絶対値が小さくなる
ように変換してその後の処理を行うことにより、流れ像
を除去するとともにアーチファクトのない良好な画像を
得ることができる。
【００２８】
【発明の実施の形態】以下、図面を参照して本発明の実
施の形態について説明する。図１は本発明による画像処
理装置の概念を表す図である。図１に示すように本発明
による画像処理装置は、入力された原画像信号Ｓorg に
基づいて互いに周波数応答特性が異なる多重解像度のボ
ケ画像信号Ｓusk （ｋ＝１〜ｎ）を作成するボケ画像信
号作成手段１と、ボケ画像信号作成手段１において作成
されたボケ画像信号Ｓusk に基づいて複数の帯域制限画
像信号を作成する帯域制限画像信号作成手段２と、帯域
制限画像信号作成手段２において作成された複数の帯域
制限画像信号のうち少なくとも１つの信号について、該
帯域制限画像信号の少なくとも一部を小さくするような
変換処理を施す変換手段３と、変換された帯域制限画像
信号を積算して積算信号を作成する積算手段４と、積算
信号に所定の強調係数を乗じて原画像信号Ｓorg と加算
することにより、原画像の高周波成分が強調された処理
済画像Ｓprocを得る周波数強調処理手段５とからなる。
施の形態について説明する。図１は本発明による画像処
理装置の概念を表す図である。図１に示すように本発明
による画像処理装置は、入力された原画像信号Ｓorg に
基づいて互いに周波数応答特性が異なる多重解像度のボ
ケ画像信号Ｓusk （ｋ＝１〜ｎ）を作成するボケ画像信
号作成手段１と、ボケ画像信号作成手段１において作成
されたボケ画像信号Ｓusk に基づいて複数の帯域制限画
像信号を作成する帯域制限画像信号作成手段２と、帯域
制限画像信号作成手段２において作成された複数の帯域
制限画像信号のうち少なくとも１つの信号について、該
帯域制限画像信号の少なくとも一部を小さくするような
変換処理を施す変換手段３と、変換された帯域制限画像
信号を積算して積算信号を作成する積算手段４と、積算
信号に所定の強調係数を乗じて原画像信号Ｓorg と加算
することにより、原画像の高周波成分が強調された処理
済画像Ｓprocを得る周波数強調処理手段５とからなる。
【００２９】はじめに、ボケ画像信号作成手段１におい
て行われる処理について説明する。図２はボケ画像信号
作成手段１の処理を示すブロック図である。図２に示す
ように原画像を表すデジタルの原画像信号Ｓorgがフィ
ルタリング処理手段10においてローパスフィルタにより
フィルタリングされる。このローパスフィルタとして
は、例えば図３に示すような５×１のグリッド状の１次
元ガウス分布に略対応したフィルタＦを用いることがで
きる。このフィルタＦは下記の式（３）
て行われる処理について説明する。図２はボケ画像信号
作成手段１の処理を示すブロック図である。図２に示す
ように原画像を表すデジタルの原画像信号Ｓorgがフィ
ルタリング処理手段10においてローパスフィルタにより
フィルタリングされる。このローパスフィルタとして
は、例えば図３に示すような５×１のグリッド状の１次
元ガウス分布に略対応したフィルタＦを用いることがで
きる。このフィルタＦは下記の式（３）
【００３０】
【数１】
【００３１】において、σ＝１としたものである。ここ
で、フィルタＦとしてガウス信号を用いるのは、ガウス
信号は周波数空間および実空間の双方において、局在性
がよいためである。
で、フィルタＦとしてガウス信号を用いるのは、ガウス
信号は周波数空間および実空間の双方において、局在性
がよいためである。
【００３２】そして、このフィルタＦにより原画像の画
素のｘ方向およびｙ方向に対してフィルタリング処理を
施すことにより、原画像信号Ｓorg 全体に対してフィル
タリング処理が施される。
素のｘ方向およびｙ方向に対してフィルタリング処理を
施すことにより、原画像信号Ｓorg 全体に対してフィル
タリング処理が施される。
【００３３】フィルタリング処理手段10においては、こ
のようなフィルタＦにより以下のようにしてフィルタリ
ング処理が施される。図４はフィルタリング処理の詳細
を表す図である。図４に示すように、原画像信号Ｓorg
に対し、図３に示すフィルタＦにより１画素おきにフィ
ルタリング処理が施される。そしてこのフィルタリング
処理により、フィルタリング処理画像信号Ｂ1が得られ
る。このフィルタリング処理画像信号Ｂ1は、原画像に
対するサイズが１／４（ｘ方向、ｙ方向にそれぞれ１／
２）となっている。次いで、再度フィルタＦによりフィ
ルタリング処理画像信号Ｂ1に対して１画素おきにフィ
ルタリング処理が施される。そしてこのようなフィルタ
Ｆによるフィルタリング処理が繰り返し施されることに
より、ｎ個のフィルタリング処理画像信号Ｂk （k=1〜
n）が得られる。このフィルタリング処理画像信号Ｂk
は、原画像に対してサイズが１／２2kとなっている。こ
の際、フィルタリング処理画像信号Ｂk の周波数特性は
図５に示されるようなものとなる。図５に示されるよう
にフィルタリング処理画像信号Ｂk のレスポンスはｋが
大きいほど高周波成分が除去されたものとなっている
（但し、図５においてはk=1〜3としている）。
のようなフィルタＦにより以下のようにしてフィルタリ
ング処理が施される。図４はフィルタリング処理の詳細
を表す図である。図４に示すように、原画像信号Ｓorg
に対し、図３に示すフィルタＦにより１画素おきにフィ
ルタリング処理が施される。そしてこのフィルタリング
処理により、フィルタリング処理画像信号Ｂ1が得られ
る。このフィルタリング処理画像信号Ｂ1は、原画像に
対するサイズが１／４（ｘ方向、ｙ方向にそれぞれ１／
２）となっている。次いで、再度フィルタＦによりフィ
ルタリング処理画像信号Ｂ1に対して１画素おきにフィ
ルタリング処理が施される。そしてこのようなフィルタ
Ｆによるフィルタリング処理が繰り返し施されることに
より、ｎ個のフィルタリング処理画像信号Ｂk （k=1〜
n）が得られる。このフィルタリング処理画像信号Ｂk
は、原画像に対してサイズが１／２2kとなっている。こ
の際、フィルタリング処理画像信号Ｂk の周波数特性は
図５に示されるようなものとなる。図５に示されるよう
にフィルタリング処理画像信号Ｂk のレスポンスはｋが
大きいほど高周波成分が除去されたものとなっている
（但し、図５においてはk=1〜3としている）。
【００３４】なお、上記実施の形態においては、図３に
示す１次元フィルタＦにより、原画像のｘ方向およびｙ
方向に対してフィルタリング処理が施されるようになっ
ているが、図６に示すような５×５の２次元フィルタに
より原画像信号Ｓorg およびフィルタリング処理画像信
号に対して一度にフィルタリング処理が施されるように
してもよい。
示す１次元フィルタＦにより、原画像のｘ方向およびｙ
方向に対してフィルタリング処理が施されるようになっ
ているが、図６に示すような５×５の２次元フィルタに
より原画像信号Ｓorg およびフィルタリング処理画像信
号に対して一度にフィルタリング処理が施されるように
してもよい。
【００３５】次いで、このようにして得られたフィルタ
リング処理画像信号Ｂk に対して図２に示される補間演
算処理手段11において補間演算処理が施され、これによ
り原画像と同一の大きさの多重解像度のボケ画像が得ら
れる。以下この補間演算処理について説明する。
リング処理画像信号Ｂk に対して図２に示される補間演
算処理手段11において補間演算処理が施され、これによ
り原画像と同一の大きさの多重解像度のボケ画像が得ら
れる。以下この補間演算処理について説明する。
【００３６】補間演算の方法としては、Ｂスプラインに
よる方法など種々の方法が挙げられるが、本発明による
実施の形態においては、ローパスフィルタとしてガウス
信号に基づくフィルタＦを用いているため、補間演算を
行うための補間係数としてもガウス信号を用いるものと
する。ここで、ガウス信号を用いた補間係数とは、下記
の式（４）
よる方法など種々の方法が挙げられるが、本発明による
実施の形態においては、ローパスフィルタとしてガウス
信号に基づくフィルタＦを用いているため、補間演算を
行うための補間係数としてもガウス信号を用いるものと
する。ここで、ガウス信号を用いた補間係数とは、下記
の式（４）
【００３７】
【数２】
【００３８】において、σ＝２k-1と近似したものを用
いる。
いる。
【００３９】フィルタリング処理画像信号Ｂ1 を補間す
る際には、ｋ＝１であるためσ＝１となる。上記式
（４）において、σ＝１としたときの補間を行うための
フィルタは、図７に示すように５×１の１次元フィルタ
となる。まずフィルタリング処理画像信号Ｂ1 に対して
１画素おきに値が０の画素を１つずつ補間することによ
りフィルタリング処理画像信号Ｂ1 が原画像と同一のサ
イズに拡大される。このように値が０の画素が補間され
たフィルタリング処理画像信号Ｂ1 を１次元的に図８に
示す。そして、この補間されたフィルタリング処理画像
信号Ｂ1 に対して上述した図７に示すフィルタＦ1 によ
りフィルタリング処理が施される。
る際には、ｋ＝１であるためσ＝１となる。上記式
（４）において、σ＝１としたときの補間を行うための
フィルタは、図７に示すように５×１の１次元フィルタ
となる。まずフィルタリング処理画像信号Ｂ1 に対して
１画素おきに値が０の画素を１つずつ補間することによ
りフィルタリング処理画像信号Ｂ1 が原画像と同一のサ
イズに拡大される。このように値が０の画素が補間され
たフィルタリング処理画像信号Ｂ1 を１次元的に図８に
示す。そして、この補間されたフィルタリング処理画像
信号Ｂ1 に対して上述した図７に示すフィルタＦ1 によ
りフィルタリング処理が施される。
【００４０】ここで、図７に示すフィルタＦ1 は５×１
のフィルタであるが、図８に示すようにフィルタリング
処理画像信号Ｂ1 は１画素おきに値が０の画素が補間さ
れている。したがって、フィルタリング処理画像Ｂ1 に
対してフィルタＦ1 により施されるフィルタリング処理
は、実質的には２×１のフィルタ（0.5, 0.5）および３
×１のフィルタ（0.1, 0.8, 0.1）の２種類のフィルタ
により、施されるフィルタリング処理と等価なものとな
る。そしてこのフィルタリング処理により、原画像信号
Ｓorg と同一データ数、すなわち原画像と同一サイズの
ボケ画像の信号Ｓus1が得られる。
のフィルタであるが、図８に示すようにフィルタリング
処理画像信号Ｂ1 は１画素おきに値が０の画素が補間さ
れている。したがって、フィルタリング処理画像Ｂ1 に
対してフィルタＦ1 により施されるフィルタリング処理
は、実質的には２×１のフィルタ（0.5, 0.5）および３
×１のフィルタ（0.1, 0.8, 0.1）の２種類のフィルタ
により、施されるフィルタリング処理と等価なものとな
る。そしてこのフィルタリング処理により、原画像信号
Ｓorg と同一データ数、すなわち原画像と同一サイズの
ボケ画像の信号Ｓus1が得られる。
【００４１】次いで、フィルタリング処理画像信号Ｂ2
に対してフィルタリング処理が施される。フィルタリン
グ処理画像信号Ｂ2 を補間する際には、ｋ＝２であるた
め、σ＝２となる。上記式（４）において、σ＝２とし
たときの補間を行うためのフィルタは、図９に示すよう
に１１×１の１次元フィルタとなる。そして、まずフィ
ルタリング処理画像信号Ｂ2 に対して図11のように１画
素おきに値が０の画素を３つずつ補間することによりフ
ィルタリング処理画像信号Ｂ2 が原画像と同一のサイズ
に拡大される。値が０の画素が補間されたこのフィルタ
リング処理画像信号Ｂ2 に対して上述した図９に示すフ
ィルタＦ2 によるフィルタリング処理が施される。
に対してフィルタリング処理が施される。フィルタリン
グ処理画像信号Ｂ2 を補間する際には、ｋ＝２であるた
め、σ＝２となる。上記式（４）において、σ＝２とし
たときの補間を行うためのフィルタは、図９に示すよう
に１１×１の１次元フィルタとなる。そして、まずフィ
ルタリング処理画像信号Ｂ2 に対して図11のように１画
素おきに値が０の画素を３つずつ補間することによりフ
ィルタリング処理画像信号Ｂ2 が原画像と同一のサイズ
に拡大される。値が０の画素が補間されたこのフィルタ
リング処理画像信号Ｂ2 に対して上述した図９に示すフ
ィルタＦ2 によるフィルタリング処理が施される。
【００４２】ここで、図９に示すフィルタＦ2 は１１×
１のフィルタであるが、図10に示すようにフィルタリン
グ処理画像信号Ｂ2 は１画素おきに値が０の画素が３つ
ずつ補間されている。したがって、フィルタリング処理
画像信号Ｂ2 に対してフィルタＦ2 により施されるフィ
ルタリング処理は、実質的には２×１のフィルタ（0.5,
0.5）および３×１のフィルタ（0.3, 0.65, 0.05）、
（0.13, 0.74, 0.13）および（0.05, 0.65, 0.3）の４
種類のフィルタにより施されるフィルタリング処理と等
価なものとなる。そしてこのフィルタリング処理によ
り、原画像信号Ｓorg と同一データ数のボケ画像信号Ｓ
us2が得られる。
１のフィルタであるが、図10に示すようにフィルタリン
グ処理画像信号Ｂ2 は１画素おきに値が０の画素が３つ
ずつ補間されている。したがって、フィルタリング処理
画像信号Ｂ2 に対してフィルタＦ2 により施されるフィ
ルタリング処理は、実質的には２×１のフィルタ（0.5,
0.5）および３×１のフィルタ（0.3, 0.65, 0.05）、
（0.13, 0.74, 0.13）および（0.05, 0.65, 0.3）の４
種類のフィルタにより施されるフィルタリング処理と等
価なものとなる。そしてこのフィルタリング処理によ
り、原画像信号Ｓorg と同一データ数のボケ画像信号Ｓ
us2が得られる。
【００４３】そしてこのようなフィルタリング処理が全
てのフィルタリング処理画像信号Ｂk に対して行われ
る。フィルタリング処理画像信号Ｂk を補間する際に
は、上記式（４）に基づいて、３×２k −１の長さのフ
ィルタを作成し、フィルタリング処理画像信号Ｂk の各
画素の間に値が０の画素を２k −１個ずつ補間すること
により、原画像と同一サイズに拡大する。この値が０の
画素が補間されたフィルタリング処理画像信号Ｂk に対
して３×２k −１の長さのフィルタにより、フィルタリ
ング処理が行われる。
てのフィルタリング処理画像信号Ｂk に対して行われ
る。フィルタリング処理画像信号Ｂk を補間する際に
は、上記式（４）に基づいて、３×２k −１の長さのフ
ィルタを作成し、フィルタリング処理画像信号Ｂk の各
画素の間に値が０の画素を２k −１個ずつ補間すること
により、原画像と同一サイズに拡大する。この値が０の
画素が補間されたフィルタリング処理画像信号Ｂk に対
して３×２k −１の長さのフィルタにより、フィルタリ
ング処理が行われる。
【００４４】ここで、この３×２k −１の長さのフィル
タにより施されるフィルタリング処理は、２k 周期で長
さが２または３のフィルタにより施されるフィルタリン
グ処理と等価なものとなる。そしてこのフィルタリング
処理により、ｎ個のボケ画像信号Ｓusk が得られる。こ
のボケ画像信号Ｓusk を可視像として表すと、結果とし
てそれぞれ解像度が異なる、すなわち周波数応答特性が
異なる多重解像度のボケ画像となる。このように、フィ
ルタは長くなるものの、実質的には長さが２または３の
フィルタによりフィルタリング処理を施していることと
同様であるため、演算量は、フィルタが長くなろうとも
それほど多くはならないものである。このため、演算量
を減らして、多重解像度のボケ画像信号Ｓusk を高速に
作成することが可能となる。
タにより施されるフィルタリング処理は、２k 周期で長
さが２または３のフィルタにより施されるフィルタリン
グ処理と等価なものとなる。そしてこのフィルタリング
処理により、ｎ個のボケ画像信号Ｓusk が得られる。こ
のボケ画像信号Ｓusk を可視像として表すと、結果とし
てそれぞれ解像度が異なる、すなわち周波数応答特性が
異なる多重解像度のボケ画像となる。このように、フィ
ルタは長くなるものの、実質的には長さが２または３の
フィルタによりフィルタリング処理を施していることと
同様であるため、演算量は、フィルタが長くなろうとも
それほど多くはならないものである。このため、演算量
を減らして、多重解像度のボケ画像信号Ｓusk を高速に
作成することが可能となる。
【００４５】なお、本実施の形態においては、長さが３
×２k −１の１次元のフィルタにより画像のｘ方向およ
びｙ方向にフィルタリング処理が施されるようになって
いるが、予め２次元状のフィルタを作成しておき、この
フィルタによりフィルタリング処理画像に対するフィル
タリング処理を施すことにより、ボケ画像信号Ｓuskが
得られるようにしてもよい。この場合、フィルタリング
処理画像に対して、補間演算を行うためのフィルタリン
グ処理に使用するフィルタは、（３×２k −１）×（３
×２k −１）のフィルタとなるが、上述した１次元のフ
ィルタを用いる場合と同様に、これは２k 周期で２×２
または３×３のフィルタにより施されるフィルタリング
処理と等価なものとなり、上述した１次元フィルタを用
いる場合と同様に、フィルタのサイズが大きくなって
も、フィルタリング処理を行うための演算量は実質的に
はそれほど大きくならないものである。
×２k −１の１次元のフィルタにより画像のｘ方向およ
びｙ方向にフィルタリング処理が施されるようになって
いるが、予め２次元状のフィルタを作成しておき、この
フィルタによりフィルタリング処理画像に対するフィル
タリング処理を施すことにより、ボケ画像信号Ｓuskが
得られるようにしてもよい。この場合、フィルタリング
処理画像に対して、補間演算を行うためのフィルタリン
グ処理に使用するフィルタは、（３×２k −１）×（３
×２k −１）のフィルタとなるが、上述した１次元のフ
ィルタを用いる場合と同様に、これは２k 周期で２×２
または３×３のフィルタにより施されるフィルタリング
処理と等価なものとなり、上述した１次元フィルタを用
いる場合と同様に、フィルタのサイズが大きくなって
も、フィルタリング処理を行うための演算量は実質的に
はそれほど大きくならないものである。
【００４６】このようにして得られたボケ画像信号Ｓus
k の周波数特性を図11に示す。図11に示すように、ボケ
画像信号Ｓusk のｋの値が大きくなるほど、原画像信号
Ｓorg の高周波成分が除去された信号となる。
k の周波数特性を図11に示す。図11に示すように、ボケ
画像信号Ｓusk のｋの値が大きくなるほど、原画像信号
Ｓorg の高周波成分が除去された信号となる。
【００４７】図12は、上記図２のボケ画像信号作成手段
を含む本発明の画像処理装置の一実施の形態について、
その全体の構成を示す図であるが、この図に示されるよ
うに、ボケ画像信号作成手段１において生成された各ボ
ケ画像信号は、次に、帯域制限画像信号作成手段２およ
び変換手段３において処理される。図12に示すように、
まず、原画像信号Ｓorgおよびボケ画像信号作成手段１
において作成された複数のボケ画像信号Ｓusk に基づい
て帯域制限画像信号が作成されるが、この帯域制限画像
信号は減算器21により、互いに隣接する周波数帯域同士
のボケ画像信号Ｓusk の減算を行うことにより得られ
る。すなわち、Ｓorg −Ｓus1 、Ｓus1 −Ｓus2 、…Ｓ
usN-1 −ＳusN を順次計算することにより、複数の帯域
制限画像信号が求められる。この帯域制限画像信号の周
波数特性を図13に示す。図13に示すように、帯域制限画
像信号はボケ画像信号Ｓusk のｋの値が大きくなるほ
ど、原画像信号Ｓorg の低周波数成分の帯域を表す信号
となる。
を含む本発明の画像処理装置の一実施の形態について、
その全体の構成を示す図であるが、この図に示されるよ
うに、ボケ画像信号作成手段１において生成された各ボ
ケ画像信号は、次に、帯域制限画像信号作成手段２およ
び変換手段３において処理される。図12に示すように、
まず、原画像信号Ｓorgおよびボケ画像信号作成手段１
において作成された複数のボケ画像信号Ｓusk に基づい
て帯域制限画像信号が作成されるが、この帯域制限画像
信号は減算器21により、互いに隣接する周波数帯域同士
のボケ画像信号Ｓusk の減算を行うことにより得られ
る。すなわち、Ｓorg −Ｓus1 、Ｓus1 −Ｓus2 、…Ｓ
usN-1 −ＳusN を順次計算することにより、複数の帯域
制限画像信号が求められる。この帯域制限画像信号の周
波数特性を図13に示す。図13に示すように、帯域制限画
像信号はボケ画像信号Ｓusk のｋの値が大きくなるほ
ど、原画像信号Ｓorg の低周波数成分の帯域を表す信号
となる。
【００４８】次いで、変換手段３において、このように
して求められた帯域制限画像信号をこの帯域制限画像信
号の信号に応じて変換する。この変換は変換器22におい
て、例えば図14に示すような関数ｆにより行う。この関
数ｆは、帯域制限画像信号の絶対値が閾値Ｔｈ１よりも
小さい場合は傾きが１であり、閾値Ｔｈ１よりも大きい
場合は傾きが１よりも小さくなるような関数である。こ
の関数は、各帯域制限画像信号において同一のものであ
ってもよいが、各信号ごとに異なるものであってもよ
い。
して求められた帯域制限画像信号をこの帯域制限画像信
号の信号に応じて変換する。この変換は変換器22におい
て、例えば図14に示すような関数ｆにより行う。この関
数ｆは、帯域制限画像信号の絶対値が閾値Ｔｈ１よりも
小さい場合は傾きが１であり、閾値Ｔｈ１よりも大きい
場合は傾きが１よりも小さくなるような関数である。こ
の関数は、各帯域制限画像信号において同一のものであ
ってもよいが、各信号ごとに異なるものであってもよ
い。
【００４９】このような関数ｆにより変換された帯域制
限画像信号は、上述した積算手段４および周波数強調処
理手段５を内包する演算器23に入力される。この演算器
23においては下記のような処理が行われる。まず、上述
したように関数ｆにより変換された帯域制限画像信号は
積算される。そしてこの積算信号が求められると、周波
数強調処理手段５において原画像信号Ｓorg の値に応じ
た強調度βが乗じられ、さらにこの強調度βが乗じられ
た積算信号が原画像信号Ｓorg と加算されて処理済画像
信号Ｓprocが得られる。
限画像信号は、上述した積算手段４および周波数強調処
理手段５を内包する演算器23に入力される。この演算器
23においては下記のような処理が行われる。まず、上述
したように関数ｆにより変換された帯域制限画像信号は
積算される。そしてこの積算信号が求められると、周波
数強調処理手段５において原画像信号Ｓorg の値に応じ
た強調度βが乗じられ、さらにこの強調度βが乗じられ
た積算信号が原画像信号Ｓorg と加算されて処理済画像
信号Ｓprocが得られる。
【００５０】以上の帯域制限画像信号作成手段２、変換
手段３、積算手段４および周波数強調処理手段５におい
て行われる処理を下記の式（５）に示す。
手段３、積算手段４および周波数強調処理手段５におい
て行われる処理を下記の式（５）に示す。
【００５１】 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)＋ｆ2(Ｓus1 −Ｓus2)＋… ＋ｆk(Ｓusk-1−Ｓusk )＋…＋ｆN(ＳusN-1−ＳusN)｝…（５ ） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） このようにして得られた処理済画像信号Ｓprocは、例え
ば図15に示すような周波数応答特性を有するものとな
る。すなわち、上述した帯域制限画像信号は、原画像の
濃度変化が比較的小さいいわゆる平坦部においては、各
周波数帯域において信号値の絶対値は小さくなるもので
ある。これに対して、濃度が急激に変化するエッジ部近
傍においては、帯域制限画像信号が比較的低周波数帯域
である場合、すなわちボケ画像信号Ｓusk を得る際のマ
スクのサイズが比較的大きい場合は、図42に示すよう
に、そのエッジ部近傍の画素についておかれたマスクに
エッジ部が含まれてしまうため、帯域制限画像信号がエ
ッジ部の影響を受けて信号値の絶対値が比較的大きなも
のとなる。このように、本来エッジ部でない部分がエッ
ジ部の濃度値の影響を受けることにより、画像処理を行
うことにより得られた画像のエッジ部にはオーバーシュ
ート、アンダーシュートなどのアーチファクトが発生し
てしまう。
数） このようにして得られた処理済画像信号Ｓprocは、例え
ば図15に示すような周波数応答特性を有するものとな
る。すなわち、上述した帯域制限画像信号は、原画像の
濃度変化が比較的小さいいわゆる平坦部においては、各
周波数帯域において信号値の絶対値は小さくなるもので
ある。これに対して、濃度が急激に変化するエッジ部近
傍においては、帯域制限画像信号が比較的低周波数帯域
である場合、すなわちボケ画像信号Ｓusk を得る際のマ
スクのサイズが比較的大きい場合は、図42に示すよう
に、そのエッジ部近傍の画素についておかれたマスクに
エッジ部が含まれてしまうため、帯域制限画像信号がエ
ッジ部の影響を受けて信号値の絶対値が比較的大きなも
のとなる。このように、本来エッジ部でない部分がエッ
ジ部の濃度値の影響を受けることにより、画像処理を行
うことにより得られた画像のエッジ部にはオーバーシュ
ート、アンダーシュートなどのアーチファクトが発生し
てしまう。
【００５２】そこで、帯域制限画像信号の絶対値が閾値
Ｔｈ１よりも大きい場合に、上述した関数ｆによりこの
絶対値が小さくなるように帯域制限画像信号を変換し、
この変換された各帯域制限画像信号を積算し、さらに所
定の強調度で強調することにより、原画像信号Ｓorg に
加算するための高周波成分に関する信号を得るようにし
たものである。
Ｔｈ１よりも大きい場合に、上述した関数ｆによりこの
絶対値が小さくなるように帯域制限画像信号を変換し、
この変換された各帯域制限画像信号を積算し、さらに所
定の強調度で強調することにより、原画像信号Ｓorg に
加算するための高周波成分に関する信号を得るようにし
たものである。
【００５３】このため、図15に示すように、エッジ部が
存在しない平坦部においては処理済画像信号Ｓprocの周
波数特性は実線で示すようなものとなるが、エッジ部近
傍の領域においては、処理済画像信号Ｓprocは図15の破
線に示すように比較的低い周波数帯域のレスポンスが低
下されたような特性を有するものとなる。これは、エッ
ジ部近傍の領域においては、ボケ画像信号（式（１）に
おけるＳus）を得る際のマスクが実際のマスクよりも小
さくされたことと同一の効果を奏するものである。
存在しない平坦部においては処理済画像信号Ｓprocの周
波数特性は実線で示すようなものとなるが、エッジ部近
傍の領域においては、処理済画像信号Ｓprocは図15の破
線に示すように比較的低い周波数帯域のレスポンスが低
下されたような特性を有するものとなる。これは、エッ
ジ部近傍の領域においては、ボケ画像信号（式（１）に
おけるＳus）を得る際のマスクが実際のマスクよりも小
さくされたことと同一の効果を奏するものである。
【００５４】したがって、エッジ部近傍の領域に対応す
る比較的信号値の絶対値の大きい帯域制限画像信号は、
原画像信号Ｓorg に加算するための高周波成分に関する
信号に対する影響力が小さくなる。このため、濃度が急
激に変化するエッジ部近傍においても、アーチファクト
の原因となる信号は影響力が弱められるため、処理を施
すことにより得られる画像をアーチファクトのない良好
な画像とすることができる。
る比較的信号値の絶対値の大きい帯域制限画像信号は、
原画像信号Ｓorg に加算するための高周波成分に関する
信号に対する影響力が小さくなる。このため、濃度が急
激に変化するエッジ部近傍においても、アーチファクト
の原因となる信号は影響力が弱められるため、処理を施
すことにより得られる画像をアーチファクトのない良好
な画像とすることができる。
【００５５】なお、上記実施の形態においては、上記式
（５）により処理済画像信号Ｓｐｒｏｃを得るようにし
ているが、下記式（６）により処理済画像信号Ｓｐｒｏ
ｃを得るようにしてもよい。式（５）と式（６）とで異
なるのは、帯域制限画像信号を得る際に、式（５）にお
いては隣接する周波数帯域同士で減算を行っているが、
式（６）においては、全ての周波数帯域のボケ画像信号
Ｓusk と、原画像信号Ｓorg とで減算処理を行っている
点が異なるものである。式（６）により得られる処理済
画像信号Ｓprocの周波数特性を図16に示す。図16に示す
ようにエッジ部が存在しない平坦部においては、処理済
画像信号Ｓprocの周波数特性は実線で示すようなものと
なるが、エッジ部近傍の領域においては、処理済画像信
号Ｓprocは図16の破線に示すように比較的低い周波数帯
域のレスポンスが低下されたような特性を有するものと
なる。これは、エッジ部近傍の領域においては、ボケ画
像信号Ｓusk を得る際のマスクが実際のマスクよりも小
さくされたことと同一の効果を奏するものである。な
お、図15と比較すると、図16の場合は全周波数帯域に亘
ってレスポンスが低下したものとなっている。このた
め、エッジ部の近傍でない平坦部においてもレスポンス
が低下したものとなっている。このため、式（６）にし
たがって、処理済画像信号Ｓprocを求める方が平坦部の
レスポンスは低下されることなくエッジ部近傍のみのレ
スポンスが低下されるため、より好ましいものである。
（５）により処理済画像信号Ｓｐｒｏｃを得るようにし
ているが、下記式（６）により処理済画像信号Ｓｐｒｏ
ｃを得るようにしてもよい。式（５）と式（６）とで異
なるのは、帯域制限画像信号を得る際に、式（５）にお
いては隣接する周波数帯域同士で減算を行っているが、
式（６）においては、全ての周波数帯域のボケ画像信号
Ｓusk と、原画像信号Ｓorg とで減算処理を行っている
点が異なるものである。式（６）により得られる処理済
画像信号Ｓprocの周波数特性を図16に示す。図16に示す
ようにエッジ部が存在しない平坦部においては、処理済
画像信号Ｓprocの周波数特性は実線で示すようなものと
なるが、エッジ部近傍の領域においては、処理済画像信
号Ｓprocは図16の破線に示すように比較的低い周波数帯
域のレスポンスが低下されたような特性を有するものと
なる。これは、エッジ部近傍の領域においては、ボケ画
像信号Ｓusk を得る際のマスクが実際のマスクよりも小
さくされたことと同一の効果を奏するものである。な
お、図15と比較すると、図16の場合は全周波数帯域に亘
ってレスポンスが低下したものとなっている。このた
め、エッジ部の近傍でない平坦部においてもレスポンス
が低下したものとなっている。このため、式（６）にし
たがって、処理済画像信号Ｓprocを求める方が平坦部の
レスポンスは低下されることなくエッジ部近傍のみのレ
スポンスが低下されるため、より好ましいものである。
【００５６】 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN ） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN ） ＝（1/N)・｛ｆ1（Ｓorg −Ｓus1)＋ｆ2（Ｓorg −Ｓus2)＋… ＋fk(Ｓorg −Ｓusk)＋…＋ｆN（Ｓorg −ＳusN)｝ …（６） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：ボケ画像信号 ｆk(k=1〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） さらに、上述した実施の形態においては、帯域制限画像
信号を変換するための関数を図14に示すように、閾値Ｔ
ｈ１よりも絶対値が大きい場合に、この絶対値を小さく
するような関数を用いているが、例えば図17に示すよう
に、帯域制限画像信号の絶対値が閾値Ｔｈ１よりも大き
い場合に、この絶対値を小さくするとともに、閾値Ｔｈ
２よりも小さい場合に、絶対値を小さくするような関数
を用いるようにしてもよい。
数） さらに、上述した実施の形態においては、帯域制限画像
信号を変換するための関数を図14に示すように、閾値Ｔ
ｈ１よりも絶対値が大きい場合に、この絶対値を小さく
するような関数を用いているが、例えば図17に示すよう
に、帯域制限画像信号の絶対値が閾値Ｔｈ１よりも大き
い場合に、この絶対値を小さくするとともに、閾値Ｔｈ
２よりも小さい場合に、絶対値を小さくするような関数
を用いるようにしてもよい。
【００５７】このように、帯域制限画像信号の絶対値が
閾値Ｔｈ１よりも小さい閾値Ｔｈ２よりも小さいほど、
この帯域制限画像信号の絶対値を小さくするように補正
することにより、画像中においてノイズと見なせるよう
な信号値の絶対値の小さい成分のレスポンスを小さくす
ることができ、これにより処理済画像のノイズを低減す
ることができる。
閾値Ｔｈ１よりも小さい閾値Ｔｈ２よりも小さいほど、
この帯域制限画像信号の絶対値を小さくするように補正
することにより、画像中においてノイズと見なせるよう
な信号値の絶対値の小さい成分のレスポンスを小さくす
ることができ、これにより処理済画像のノイズを低減す
ることができる。
【００５８】また、上述した実施の形態においては、上
記式（５）または式（６）により処理済画像信号Ｓproc
を得るようにしているが、下記式（７）により処理済画
像信号Ｓprocを得るようにしてもよい。上述した式
（６）においては、帯域制限画像信号を得るために、原
画像信号Ｓorg からボケ画像信号Ｓusk を減算するよう
にしているが、式（７）においては、ボケ画像信号Ｓus
1 からボケ画像信号Ｓusk（k=2〜N）を減算するように
している点が異なるものである。式（７）により得られ
る処理済画像信号Ｓprocの周波数特性を図18に示す。図
18に示すように式（７）により得られる処理済画像信号
Ｓprocの周波数特性は高周波成分が除去された特性を有
するものとなっている。
記式（５）または式（６）により処理済画像信号Ｓproc
を得るようにしているが、下記式（７）により処理済画
像信号Ｓprocを得るようにしてもよい。上述した式
（６）においては、帯域制限画像信号を得るために、原
画像信号Ｓorg からボケ画像信号Ｓusk を減算するよう
にしているが、式（７）においては、ボケ画像信号Ｓus
1 からボケ画像信号Ｓusk（k=2〜N）を減算するように
している点が異なるものである。式（７）により得られ
る処理済画像信号Ｓprocの周波数特性を図18に示す。図
18に示すように式（７）により得られる処理済画像信号
Ｓprocの周波数特性は高周波成分が除去された特性を有
するものとなっている。
【００５９】式（６）により処理を行うものにおいて
は、画像中のノイズと見なせる高周波成分をも強調して
しまうため、得られる処理済画像にノイズが目立つもの
となってしまうことがある。これに対して式（７）によ
り処理を行うものにおいては、図18に示すように処理済
画像信号Ｓprocの高周波成分が除去されているため、ノ
イズが強調されることがなくなりより高画質の処理済画
像を得ることができる。
は、画像中のノイズと見なせる高周波成分をも強調して
しまうため、得られる処理済画像にノイズが目立つもの
となってしまうことがある。これに対して式（７）によ
り処理を行うものにおいては、図18に示すように処理済
画像信号Ｓprocの高周波成分が除去されているため、ノ
イズが強調されることがなくなりより高画質の処理済画
像を得ることができる。
【００６０】 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓus1,Ｓus2,…ＳusN ） ＝（1/N)・｛ｆ2（Ｓus1 −Ｓus2)＋ｆ3（Ｓus1 −Ｓus3)＋… ＋fk(Ｓus1 −Ｓusk)＋…＋ｆN（Ｓus1 −ＳusN)｝ …（７） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=2〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） さらに、下記の式（８）により処理を行うようにしても
よい。上述した式（５）においては、帯域制限画像信号
（Ｓorg −Ｓus1 ）を用いているのに対して式（８）に
おいては、帯域制限画像信号（Ｓorg −Ｓus1 ）を用い
ないものである。これにより図13における最高周波数成
分が除去されることとなるため、式（７）により処理を
行うものと同様に、処理済画像信号Ｓprocの高周波成分
が除去されてノイズが強調されることがなくなりより高
画質の処理済画像を得ることができる。
数） さらに、下記の式（８）により処理を行うようにしても
よい。上述した式（５）においては、帯域制限画像信号
（Ｓorg −Ｓus1 ）を用いているのに対して式（８）に
おいては、帯域制限画像信号（Ｓorg −Ｓus1 ）を用い
ないものである。これにより図13における最高周波数成
分が除去されることとなるため、式（７）により処理を
行うものと同様に、処理済画像信号Ｓprocの高周波成分
が除去されてノイズが強調されることがなくなりより高
画質の処理済画像を得ることができる。
【００６１】 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓus1,Ｓus2,…ＳusN ） ＝｛ｆ2（Ｓus1 −Ｓus2)＋ｆ3（Ｓus2 −Ｓus3)＋… ＋fk(Ｓusk-1 −Ｓusk)＋…＋ｆN（ＳusN-1 −ＳusN)｝ …（８） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=2〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） なお、原画像を得る際の撮影部位によって、診断に必要
な周波数帯域は異なるものである。例えば、肺と胃の二
重造影の画像とにおいて、肺は比較的低周波成分を強調
するのが好ましく、胃は胃壁のヒダの観察を行うため比
較的高周波成分を強調することが好ましい。一方、骨や
人工骨等の金属を含む画像においては、これらのエッジ
部分が強調すぎることによるアーチファクトを防止する
必要があるが、例えばマンモグラフィのように骨などの
エッジ部分を含まない画像に対しては、アーチファクト
は発生しにくく、骨を含む画像と同様に周波数成分の強
調を行うと、観察に必要な部分が逆に見にくくなるとい
う問題もある。したがって、原画像を得る際の撮影部
位、あるいは帯域制限画像信号の周波数帯域に応じて図
14に示す関数ｆの形状を変化させることが望ましい。
数） なお、原画像を得る際の撮影部位によって、診断に必要
な周波数帯域は異なるものである。例えば、肺と胃の二
重造影の画像とにおいて、肺は比較的低周波成分を強調
するのが好ましく、胃は胃壁のヒダの観察を行うため比
較的高周波成分を強調することが好ましい。一方、骨や
人工骨等の金属を含む画像においては、これらのエッジ
部分が強調すぎることによるアーチファクトを防止する
必要があるが、例えばマンモグラフィのように骨などの
エッジ部分を含まない画像に対しては、アーチファクト
は発生しにくく、骨を含む画像と同様に周波数成分の強
調を行うと、観察に必要な部分が逆に見にくくなるとい
う問題もある。したがって、原画像を得る際の撮影部
位、あるいは帯域制限画像信号の周波数帯域に応じて図
14に示す関数ｆの形状を変化させることが望ましい。
【００６２】例えば、図19に示すように、骨を含む画像
においては関数ｆをＢのようにして帯域制限画像信号の
高周波成分を抑制してエッジ部分にアーチファクトを発
生しにくくすることが好ましい。逆に、マンモグラフィ
のように骨を含まない画像においては関数ｆをＡのよう
にして高周波成分をも強調し、略全周波数帯域に亘って
帯域制限画像信号の絶対値が大きくなるように変換処理
を行うことが好ましい。
においては関数ｆをＢのようにして帯域制限画像信号の
高周波成分を抑制してエッジ部分にアーチファクトを発
生しにくくすることが好ましい。逆に、マンモグラフィ
のように骨を含まない画像においては関数ｆをＡのよう
にして高周波成分をも強調し、略全周波数帯域に亘って
帯域制限画像信号の絶対値が大きくなるように変換処理
を行うことが好ましい。
【００６３】また、胸部画像の場合、全周波数帯域に亘
って帯域制限画像信号が強調されることが好ましいた
め、帯域制限画像信号の周波数帯域に応じて関数ｆを図
20に示すように変化させることが好ましい。一方、胃の
二重造影像の場合は帯域制限画像信号の高周波成分を強
調し、低周波成分を抑制した方が好ましいため、帯域制
限画像信号の周波数帯域に応じて関数ｆを図21に示すよ
うに変化させることが好ましい。帯域制限画像信号を図
20に示す関数ｆにしたがって変換することにより得られ
た処理済画像信号Ｓprocの周波数特性を図22に、帯域制
限画像信号を図21に示す関数ｆにしたがって変換するこ
とにより得られた処理済画像信号Ｓprocの周波数特性を
図23にそれぞれ示す。図22の胸部の撮影の場合は、略全
周波数帯域に亘ってレスポンスが強調されているのに対
して、図23の胃の二重造影像の場合は、高周波帯域のレ
スポンスが他の周波数帯域よりも大きく強調されている
ことが分かる。
って帯域制限画像信号が強調されることが好ましいた
め、帯域制限画像信号の周波数帯域に応じて関数ｆを図
20に示すように変化させることが好ましい。一方、胃の
二重造影像の場合は帯域制限画像信号の高周波成分を強
調し、低周波成分を抑制した方が好ましいため、帯域制
限画像信号の周波数帯域に応じて関数ｆを図21に示すよ
うに変化させることが好ましい。帯域制限画像信号を図
20に示す関数ｆにしたがって変換することにより得られ
た処理済画像信号Ｓprocの周波数特性を図22に、帯域制
限画像信号を図21に示す関数ｆにしたがって変換するこ
とにより得られた処理済画像信号Ｓprocの周波数特性を
図23にそれぞれ示す。図22の胸部の撮影の場合は、略全
周波数帯域に亘ってレスポンスが強調されているのに対
して、図23の胃の二重造影像の場合は、高周波帯域のレ
スポンスが他の周波数帯域よりも大きく強調されている
ことが分かる。
【００６４】このように、胸部画像の場合に全周波数帯
域に亘ってレスポンスが強調されるように処理を施すこ
とにより、肺野がより観察し易いものとなり、胃の二重
造影像において高周波帯域のレスポンスが強調されるよ
うに処理を施すことにより、胃壁のヒダがより観察し易
いものとなる。そしてこのように、帯域制限画像信号の
周波数帯域に応じてあるいは撮影部位に応じて帯域制限
画像信号の絶対値を変化させるように処理を行うことに
より、撮影部位あるいは周波数帯域に応じたより観察に
適した画像を得ることができることとなる。
域に亘ってレスポンスが強調されるように処理を施すこ
とにより、肺野がより観察し易いものとなり、胃の二重
造影像において高周波帯域のレスポンスが強調されるよ
うに処理を施すことにより、胃壁のヒダがより観察し易
いものとなる。そしてこのように、帯域制限画像信号の
周波数帯域に応じてあるいは撮影部位に応じて帯域制限
画像信号の絶対値を変化させるように処理を行うことに
より、撮影部位あるいは周波数帯域に応じたより観察に
適した画像を得ることができることとなる。
【００６５】ここで、上記変換手段３における帯域制限
画像信号の変換は、上述の閾値Ｔｈ１に制限されること
なく、周波数帯域ごとに異なる複数の関数f1〜fNに基づ
いて、帯域制限画像信号の絶対値に基づいて決まる該絶
対値以下の値となるように変換してもよい。この関数f1
〜fNは、画像処理の目的に応じて適切な関数の組み合わ
せを設定することが望ましい。
画像信号の変換は、上述の閾値Ｔｈ１に制限されること
なく、周波数帯域ごとに異なる複数の関数f1〜fNに基づ
いて、帯域制限画像信号の絶対値に基づいて決まる該絶
対値以下の値となるように変換してもよい。この関数f1
〜fNは、画像処理の目的に応じて適切な関数の組み合わ
せを設定することが望ましい。
【００６６】この閾値にとらわれない実施の形態におい
ては、得られる処理済画像信号Ｓproc の周波数特性
は、変換手段３において使用される関数を周波数帯域ご
とに異ならしめることにより任意の周波数特性になるよ
うに調整可能である。したがって、前記各装置において
それぞれ処理対象である画像に求められる条件に応じ
て、処理済画像信号Ｓproc の周波数特性を調整するこ
とが可能となる。
ては、得られる処理済画像信号Ｓproc の周波数特性
は、変換手段３において使用される関数を周波数帯域ご
とに異ならしめることにより任意の周波数特性になるよ
うに調整可能である。したがって、前記各装置において
それぞれ処理対象である画像に求められる条件に応じ
て、処理済画像信号Ｓproc の周波数特性を調整するこ
とが可能となる。
【００６７】以下、前記関数を周波数帯域ごとに異なら
しめることにより得られる効果について説明する。図28
は、関数f1〜fNを全て同一の関数とする、すなわち全周
波数帯域の帯域制限画像信号を同じ変換方法で変換する
従来の方法の問題点を示す図である。これは画像の濃度
が急激に変化しているエッジ部の近傍における処理を段
階的に示したものであるが、階段状の原画像信号とその
原画像信号に基づいて作成されたボケ画像信号が(a)で
あり、これに関する帯域制限画像信号が(b)、変換画像
信号が(c)およびこの変換画像信号を積算して得た信号
が(d)としてそれぞれ示されている。
しめることにより得られる効果について説明する。図28
は、関数f1〜fNを全て同一の関数とする、すなわち全周
波数帯域の帯域制限画像信号を同じ変換方法で変換する
従来の方法の問題点を示す図である。これは画像の濃度
が急激に変化しているエッジ部の近傍における処理を段
階的に示したものであるが、階段状の原画像信号とその
原画像信号に基づいて作成されたボケ画像信号が(a)で
あり、これに関する帯域制限画像信号が(b)、変換画像
信号が(c)およびこの変換画像信号を積算して得た信号
が(d)としてそれぞれ示されている。
【００６８】この図に示されるように、通常変換画像信
号を積算して得た信号は周波数帯域の境界部において不
自然なつなぎ目ができてしまい、これがすじ状のアーチ
ファクトの原因となる。これを防止するためには、この
境界部ができるだけ自然につながるよう考慮して変換画
像信号を作成しなければならない。しかしながら、関数
が一意に定められている場合には各帯域制限画像信号に
対して境界部の影響を考慮した任意の変換を行うことが
できず、結果としてすじ状のアーチファクトを防止する
ことができなかった。これに対し、本実施の形態は、前
記関数を周波数帯域ごとに異なるものとし、これらの関
数を前記周波数帯域の境界部を考慮して設定することに
より、このすじ状のアーチファクトを防止するものであ
る。
号を積算して得た信号は周波数帯域の境界部において不
自然なつなぎ目ができてしまい、これがすじ状のアーチ
ファクトの原因となる。これを防止するためには、この
境界部ができるだけ自然につながるよう考慮して変換画
像信号を作成しなければならない。しかしながら、関数
が一意に定められている場合には各帯域制限画像信号に
対して境界部の影響を考慮した任意の変換を行うことが
できず、結果としてすじ状のアーチファクトを防止する
ことができなかった。これに対し、本実施の形態は、前
記関数を周波数帯域ごとに異なるものとし、これらの関
数を前記周波数帯域の境界部を考慮して設定することに
より、このすじ状のアーチファクトを防止するものであ
る。
【００６９】次に、このような周波数ごとに異なる関数
について例を示して説明する。図24は、上述のように変
換画像信号の作成に使用される関数の一例であり、横軸
に処理される帯域制限画像信号が、縦軸にその帯域制限
画像信号が処理された値がそれぞれ対応している。これ
らの関数は、帯域制限画像信号を、その帯域制限画像信
号の絶対値に基づいて決まる、その絶対値以下の値とな
るように変換するもので、帯域制限画像信号の絶対値が
所定の値より大きい帯域制限画像信号については、変換
画像信号の値がほぼ一定の値となるように変換を行う関
数であり、関数が高周波帯域を処理する関数であるほ
ど、前記所定の値が小さい値であることを特徴とするも
のである。
について例を示して説明する。図24は、上述のように変
換画像信号の作成に使用される関数の一例であり、横軸
に処理される帯域制限画像信号が、縦軸にその帯域制限
画像信号が処理された値がそれぞれ対応している。これ
らの関数は、帯域制限画像信号を、その帯域制限画像信
号の絶対値に基づいて決まる、その絶対値以下の値とな
るように変換するもので、帯域制限画像信号の絶対値が
所定の値より大きい帯域制限画像信号については、変換
画像信号の値がほぼ一定の値となるように変換を行う関
数であり、関数が高周波帯域を処理する関数であるほ
ど、前記所定の値が小さい値であることを特徴とするも
のである。
【００７０】言い換えれば、これらの関数はそれぞれ、
原点を通り、関数の傾きがその関数により処理される帯
域制限画像信号の値にかかわらず１以下であり、その関
数により処理される帯域制限画像信号の値の絶対値の増
加にともないその関数の傾きが０になるかまたは０に収
束するような関数であって、その傾きが０または０に近
い所定の値となるときの前記処理される値の絶対値が、
高周波帯域を処理する関数であるほど小さい値であるこ
とを特徴とする。
原点を通り、関数の傾きがその関数により処理される帯
域制限画像信号の値にかかわらず１以下であり、その関
数により処理される帯域制限画像信号の値の絶対値の増
加にともないその関数の傾きが０になるかまたは０に収
束するような関数であって、その傾きが０または０に近
い所定の値となるときの前記処理される値の絶対値が、
高周波帯域を処理する関数であるほど小さい値であるこ
とを特徴とする。
【００７１】この図24の関数は、振幅の大きな帯域制限
画像信号を抑制するような変換を行うものであり、周波
数帯域の高い帯域制限画像信号の抑制の度合いを、周波
数帯域の低い帯域制限画像信号よりも強くするものであ
るが、これは実際の放射線画像のエッジに含まれている
高周波成分が、低周波成分に比べてその振幅が小さいと
いうことを考慮したものである。すなわち図29に示され
るように、実際の放射線画像においては、かなり急峻な
エッジでさえも、図29(a)のように正確な階段状にはな
っておらず、同図(b)のように高周波成分になるほどそ
の振幅が小さくなっていることが多い。このため、各周
波数成分の振幅に合わせて、周波数の高い帯域制限画像
信号ほど小さい振幅から抑制を行うことが望ましく、本
関数によりそれを実現することができる。
画像信号を抑制するような変換を行うものであり、周波
数帯域の高い帯域制限画像信号の抑制の度合いを、周波
数帯域の低い帯域制限画像信号よりも強くするものであ
るが、これは実際の放射線画像のエッジに含まれている
高周波成分が、低周波成分に比べてその振幅が小さいと
いうことを考慮したものである。すなわち図29に示され
るように、実際の放射線画像においては、かなり急峻な
エッジでさえも、図29(a)のように正確な階段状にはな
っておらず、同図(b)のように高周波成分になるほどそ
の振幅が小さくなっていることが多い。このため、各周
波数成分の振幅に合わせて、周波数の高い帯域制限画像
信号ほど小さい振幅から抑制を行うことが望ましく、本
関数によりそれを実現することができる。
【００７２】次に図25の関数について説明する。図25の
関数は、帯域制限画像信号を、帯域制限画像信号の絶対
値に基づいて決まる、その絶対値以下の値となるように
変換を行うもので、この関数が低周波帯域を処理する関
数であるほど、帯域制限画像信号の絶対値が０近傍の所
定の範囲内の値である帯域制限画像信号を変換した際に
得られる変換画像信号の絶対値が小さい値であることを
特徴とするものである。
関数は、帯域制限画像信号を、帯域制限画像信号の絶対
値に基づいて決まる、その絶対値以下の値となるように
変換を行うもので、この関数が低周波帯域を処理する関
数であるほど、帯域制限画像信号の絶対値が０近傍の所
定の範囲内の値である帯域制限画像信号を変換した際に
得られる変換画像信号の絶対値が小さい値であることを
特徴とするものである。
【００７３】言い換えれば、これらの関数はそれぞれ、
原点を通り、関数の傾きがその関数により処理される値
にかかわらず１以下であり、その関数の０近傍における
傾きが、低周波帯域を処理する関数であるほど小さいこ
とを特徴とするものである。
原点を通り、関数の傾きがその関数により処理される値
にかかわらず１以下であり、その関数の０近傍における
傾きが、低周波帯域を処理する関数であるほど小さいこ
とを特徴とするものである。
【００７４】これらの関数は、図28に示される変換画像
信号を積算して得られる信号(d)を、原画像信号Ｓorg
に加えた場合、原画像信号Ｓorg と加算された信号との
つなぎ目、すなわち信号の立ち上がりをより自然なもの
とするという効果がある。
信号を積算して得られる信号(d)を、原画像信号Ｓorg
に加えた場合、原画像信号Ｓorg と加算された信号との
つなぎ目、すなわち信号の立ち上がりをより自然なもの
とするという効果がある。
【００７５】図26に示される関数は、周波数強調処理装
置において使用される関数で、図24および図25の関数の
両方の特徴を兼ね備えたものであり、両方の関数の効果
を得ることができるものである。
置において使用される関数で、図24および図25の関数の
両方の特徴を兼ね備えたものであり、両方の関数の効果
を得ることができるものである。
【００７６】以上のように本発明の画像処理方法および
装置により、各周波数帯域において使用される関数を目
的に応じて異ならしめて全体の周波数特性を任意のもの
とすることが可能となる。図27はその効果の一例を示す
ものであるが、この図には本発明の画像処理装置により
作成された各周波数帯域ごとの信号の周波数特性(a-
1)、(b-1)、(c-1)と、それに対応する全体の周波数特性
(a-2)、(b-2)、(c-2)が示されている。(a-1)、(a-2)は
全ての周波数帯域で傾き１の関数を設定した場合、(b-
1)、(b-2)は低周波帯域ほど傾きを小さく設定した場
合、(c-1)、(c-2)は特定の周波数帯域のみ関数の傾きを
１とし、それ以外を１より小さくすることにより全体と
してバンドパスの周波数特性を得た場合であり(c)の例
は、ある特定の周波数帯域の情報のみを得たい場合など
に有用である。
装置により、各周波数帯域において使用される関数を目
的に応じて異ならしめて全体の周波数特性を任意のもの
とすることが可能となる。図27はその効果の一例を示す
ものであるが、この図には本発明の画像処理装置により
作成された各周波数帯域ごとの信号の周波数特性(a-
1)、(b-1)、(c-1)と、それに対応する全体の周波数特性
(a-2)、(b-2)、(c-2)が示されている。(a-1)、(a-2)は
全ての周波数帯域で傾き１の関数を設定した場合、(b-
1)、(b-2)は低周波帯域ほど傾きを小さく設定した場
合、(c-1)、(c-2)は特定の周波数帯域のみ関数の傾きを
１とし、それ以外を１より小さくすることにより全体と
してバンドパスの周波数特性を得た場合であり(c)の例
は、ある特定の周波数帯域の情報のみを得たい場合など
に有用である。
【００７７】なお、上記閾値より大きい値を抑制する実
施の形態および閾値によらず全体を抑制する実施の形態
の２つの実施の形態では、変換手段３により行う変換処
理を強調係数βに応じて異なる処理としてもよい。以
下、その効果を閾値によらない関数を用いる場合を例に
説明する。
施の形態および閾値によらず全体を抑制する実施の形態
の２つの実施の形態では、変換手段３により行う変換処
理を強調係数βに応じて異なる処理としてもよい。以
下、その効果を閾値によらない関数を用いる場合を例に
説明する。
【００７８】この関数f1〜fNとして、上記強調係数βが
比較的小さい場合には図30に示されるような関数を、ま
た強調係数βが大きい場合には図31に示すような関数を
用いている。ここで、図中の「高」「低」はそれぞれ、
高周波帯域の帯域制限画像信号の変換に使用される関
数、低周波帯域の帯域制限画像信号の変換に使用される
関数を示している。
比較的小さい場合には図30に示されるような関数を、ま
た強調係数βが大きい場合には図31に示すような関数を
用いている。ここで、図中の「高」「低」はそれぞれ、
高周波帯域の帯域制限画像信号の変換に使用される関
数、低周波帯域の帯域制限画像信号の変換に使用される
関数を示している。
【００７９】これらの関数はいずれも帯域制限画像信号
をその帯域制限画像信号の絶対値よりも小さくなるよう
に抑制するものであるが、本実施の形態では、図30およ
び図31とを比較することによりわかるように、強調係数
βが大きいときほど、高周波帯域を処理する関数による
抑制の度合いと、低周波帯域を処理する関数による抑制
の度合いとの差が大きくなるように関数の組み合わせを
定義している。具体的には、高周波帯域を処理する関数
は強調係数に拘わらず一定とし、低周波帯域を処理する
関数の抑制の度合いを強調係数が大きくなるほど強くす
ることにより差を広げている。但し、これは低周波帯域
を処理する関数を一定として高周波帯域を処理する関数
の抑制の度合いを緩めてもよく、また全ての関数を抑制
の度合いの差が広がるように変化させてもよい。
をその帯域制限画像信号の絶対値よりも小さくなるよう
に抑制するものであるが、本実施の形態では、図30およ
び図31とを比較することによりわかるように、強調係数
βが大きいときほど、高周波帯域を処理する関数による
抑制の度合いと、低周波帯域を処理する関数による抑制
の度合いとの差が大きくなるように関数の組み合わせを
定義している。具体的には、高周波帯域を処理する関数
は強調係数に拘わらず一定とし、低周波帯域を処理する
関数の抑制の度合いを強調係数が大きくなるほど強くす
ることにより差を広げている。但し、これは低周波帯域
を処理する関数を一定として高周波帯域を処理する関数
の抑制の度合いを緩めてもよく、また全ての関数を抑制
の度合いの差が広がるように変化させてもよい。
【００８０】以上説明したような関数により高周波成分
を強調された画像信号を、レスポンスとして示したもの
が図32および図33である。図32は、上記図30に示される
関数の組み合わせを使用した場合であり、図33は図31に
示される関数の組み合わせを使用した場合である。これ
らの図によれば、強調係数を大きくした際に高周波成分
のみ強調の度合いが強くなり、低周波成分の強調の度合
いは強調係数が小さい場合とそれほど変わらないように
なっていることがわかる。
を強調された画像信号を、レスポンスとして示したもの
が図32および図33である。図32は、上記図30に示される
関数の組み合わせを使用した場合であり、図33は図31に
示される関数の組み合わせを使用した場合である。これ
らの図によれば、強調係数を大きくした際に高周波成分
のみ強調の度合いが強くなり、低周波成分の強調の度合
いは強調係数が小さい場合とそれほど変わらないように
なっていることがわかる。
【００８１】強調係数βによらずに関数の組み合わせを
決定する場合には、強調係数を大きくすると全ての周波
数帯域において強調の度合いが強くなる。一般に強調係
数を大きくするのは、情報量が少ない高周波成分を強調
するためであるが、これでは、低周波成分までもが強調
されてしまうことになり、その結果コントラストがつき
すぎ、アーチファクトとなることがある。したがって、
図31に示されるように、強調係数βを大きくしても低周
波成分の強調の度合いはあまり変化させないようにする
ことにより、このアーチファクトの発生を防止すること
ができる。
決定する場合には、強調係数を大きくすると全ての周波
数帯域において強調の度合いが強くなる。一般に強調係
数を大きくするのは、情報量が少ない高周波成分を強調
するためであるが、これでは、低周波成分までもが強調
されてしまうことになり、その結果コントラストがつき
すぎ、アーチファクトとなることがある。したがって、
図31に示されるように、強調係数βを大きくしても低周
波成分の強調の度合いはあまり変化させないようにする
ことにより、このアーチファクトの発生を防止すること
ができる。
【００８２】なお、上記説明においては、強調係数βが
小さい場合と大きい場合用に１組ずつ関数の組み合わせ
を示したが、これは関数の組み合わせが２組であること
を意味するものではなく、強調係数βの変化に応じて関
数の特徴をどのように変えればよいかを示すための例に
すぎない。したがって、強調係数βの変化の度合いを何
段階かのレベルに分類し、各レベルに応じて関数の組み
合わせを決定するといったことも可能であり、そのよう
な段階をより細かく設定することにより、より高画質な
画像が得られることはいうまでもない。
小さい場合と大きい場合用に１組ずつ関数の組み合わせ
を示したが、これは関数の組み合わせが２組であること
を意味するものではなく、強調係数βの変化に応じて関
数の特徴をどのように変えればよいかを示すための例に
すぎない。したがって、強調係数βの変化の度合いを何
段階かのレベルに分類し、各レベルに応じて関数の組み
合わせを決定するといったことも可能であり、そのよう
な段階をより細かく設定することにより、より高画質な
画像が得られることはいうまでもない。
【００８３】また、上記実施の形態では、強調係数を大
きくするのは高周波成分のみを強調したいためであると
の解釈に基づいて非線形関数を決定しているが、強調係
数を変化させる目的はこれに限られるものではなく、し
たがって非線形関数としては目的に応じて様々なものが
適用可能である。この際、上記例では非線形関数は各周
波数帯域ごとに異なる複数の非線形関数の組としている
が、達成したい目的によっては、全ての周波数帯域に同
一の非線形関数を用いても本発明の効果を得ることがで
きる。
きくするのは高周波成分のみを強調したいためであると
の解釈に基づいて非線形関数を決定しているが、強調係
数を変化させる目的はこれに限られるものではなく、し
たがって非線形関数としては目的に応じて様々なものが
適用可能である。この際、上記例では非線形関数は各周
波数帯域ごとに異なる複数の非線形関数の組としている
が、達成したい目的によっては、全ての周波数帯域に同
一の非線形関数を用いても本発明の効果を得ることがで
きる。
【００８４】次に、本発明の画像処理方法および装置の
他の実施の形態について説明する。但し、本実施の形態
におけるボケ画像信号作成手段１、帯域制限画像信号作
成手段２、積算手段４、周波数強調処理手段５の処理は
上記実施の形態と同じであるため、ここでは説明を省略
し、変換手段３の処理についてのみ、以下説明する。
他の実施の形態について説明する。但し、本実施の形態
におけるボケ画像信号作成手段１、帯域制限画像信号作
成手段２、積算手段４、周波数強調処理手段５の処理は
上記実施の形態と同じであるため、ここでは説明を省略
し、変換手段３の処理についてのみ、以下説明する。
【００８５】図34はこの実施の形態における本発明の画
像処理装置の全体の構成を示す図である。上記実施の形
態と同様、原画像信号Ｓorgおよびボケ画像信号作成手
段１において作成された複数のボケ画像信号Ｓusk に基
づいて帯域制限画像信号が作成される。この帯域制限画
像信号は減算器21により、互いに隣接する周波数帯域の
２つのボケ画像信号Ｓusk （但し原画像信号Ｓorgにつ
いてはＳorgとＳus1 )の減算を行うことにより得られ
る。すなわち、Ｓorg −Ｓus1 、Ｓus1 −Ｓus2、…Ｓu
sN-1 −ＳusN を順次計算することにより、複数の帯域
制限画像信号が求められる。本実施の形態では、例えば
帯域制限画像信号Ｓus1 −Ｓus2 に対し、これに対応す
る補助画像信号をＳus2 −Ｓus3 としている。したがっ
て、帯域制限画像信号作成手段と補助画像信号作成手段
は実質的に同一の手段を兼用するものとする。つまり、
帯域制限画像信号作成手段により作成された信号は、帯
域制限画像信号として処理されると同時に、隣接する帯
域制限画像信号に対応する補助画像信号としても処理さ
れるものとする。
像処理装置の全体の構成を示す図である。上記実施の形
態と同様、原画像信号Ｓorgおよびボケ画像信号作成手
段１において作成された複数のボケ画像信号Ｓusk に基
づいて帯域制限画像信号が作成される。この帯域制限画
像信号は減算器21により、互いに隣接する周波数帯域の
２つのボケ画像信号Ｓusk （但し原画像信号Ｓorgにつ
いてはＳorgとＳus1 )の減算を行うことにより得られ
る。すなわち、Ｓorg −Ｓus1 、Ｓus1 −Ｓus2、…Ｓu
sN-1 −ＳusN を順次計算することにより、複数の帯域
制限画像信号が求められる。本実施の形態では、例えば
帯域制限画像信号Ｓus1 −Ｓus2 に対し、これに対応す
る補助画像信号をＳus2 −Ｓus3 としている。したがっ
て、帯域制限画像信号作成手段と補助画像信号作成手段
は実質的に同一の手段を兼用するものとする。つまり、
帯域制限画像信号作成手段により作成された信号は、帯
域制限画像信号として処理されると同時に、隣接する帯
域制限画像信号に対応する補助画像信号としても処理さ
れるものとする。
【００８６】上記のようにして求められた帯域制限画像
信号は、変換手段３において変換される。この変換は各
帯域制限画像信号に対して、図34に示される変換器22お
よび変換器24により、周波数帯域ごとにそれぞれ、変換
器22においては関数f1〜fNにより変換処理を施して抑制
画像信号を得、変換器24においては関数ｇにより変換処
理を施して倍率信号を得て、その抑制画像信号と倍率信
号とを乗ずることにより行われる。この際、乗算は、図
34に示されるように、例えば帯域制限画像信号Ｓus1 −
Ｓus2 を変換して得た抑制画像信号と、補助画像信号Ｓ
us2 −Ｓus3 を変換して得た倍率信号というように、隣
接した周波数帯域の信号どうしで行われる。ここで、関
数ｆkにより変換を行う変換器22が、前記抑制画像信号
作成手段に相当し、関数ｇにより変換を行う変換器24
が、前記倍率信号作成手段に相当し、それらの変換器の
出力信号を乗ずる乗算器25が、前記乗算手段に相当す
る。
信号は、変換手段３において変換される。この変換は各
帯域制限画像信号に対して、図34に示される変換器22お
よび変換器24により、周波数帯域ごとにそれぞれ、変換
器22においては関数f1〜fNにより変換処理を施して抑制
画像信号を得、変換器24においては関数ｇにより変換処
理を施して倍率信号を得て、その抑制画像信号と倍率信
号とを乗ずることにより行われる。この際、乗算は、図
34に示されるように、例えば帯域制限画像信号Ｓus1 −
Ｓus2 を変換して得た抑制画像信号と、補助画像信号Ｓ
us2 −Ｓus3 を変換して得た倍率信号というように、隣
接した周波数帯域の信号どうしで行われる。ここで、関
数ｆkにより変換を行う変換器22が、前記抑制画像信号
作成手段に相当し、関数ｇにより変換を行う変換器24
が、前記倍率信号作成手段に相当し、それらの変換器の
出力信号を乗ずる乗算器25が、前記乗算手段に相当す
る。
【００８７】関数f1〜fNは、全て同じ関数であってもよ
いし、異なる関数であってもよく、画像処理の目的に応
じて任意に設定可能なものとする。本実施の形態では、
この関数f1〜fNとして、帯域制限画像信号をその帯域制
限画像信号の絶対値よりも小さくなるように抑制する関
数で、周波数帯域ごとに異なる関数を用いており、例え
ば図24に示されるような関数が使用される。
いし、異なる関数であってもよく、画像処理の目的に応
じて任意に設定可能なものとする。本実施の形態では、
この関数f1〜fNとして、帯域制限画像信号をその帯域制
限画像信号の絶対値よりも小さくなるように抑制する関
数で、周波数帯域ごとに異なる関数を用いており、例え
ば図24に示されるような関数が使用される。
【００８８】関数ｇとしては、例えば図35に示されるよ
うなものが用いられる。この図は、この関数ｇによって
補助画像信号を変換する場合、補助画像信号の絶対値が
小さいときには、変換後の値として１に近い値が得ら
れ、絶対値が大きいほど０に近い値が得られることを示
している。なお、Ｋは変換後の値が０となる値のうち最
小の値を示している。
うなものが用いられる。この図は、この関数ｇによって
補助画像信号を変換する場合、補助画像信号の絶対値が
小さいときには、変換後の値として１に近い値が得ら
れ、絶対値が大きいほど０に近い値が得られることを示
している。なお、Ｋは変換後の値が０となる値のうち最
小の値を示している。
【００８９】上記変換手段３により得られた変換画像信
号は、演算器23に入力される。演算器23は、積算手段４
および周波数強調処理手段５を内包するものである。演
算器23においては、複数の変換画像信号が積算され、こ
の積算により得られた積算値に原画像信号Ｓorgに応じ
て決定される強調係数βが乗じられ、さらにこの強調係
数βが乗じられた積算値が原画像信号Ｓorg と加算され
て処理済画像信号Ｓproc が得られる。この処理を、上
述した他の処理とともに数式として表すと、下記の式
（９） Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)・ｇ（Ｓus1−Ｓus2） ＋ｆ2(Ｓus1 −Ｓus2)・ｇ（Ｓus2−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk )・ｇ（Ｓusk−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（ＳusN−ＳusN+1）｝ …（９ ） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換して前記
抑制信号を作成するための関数 ｇ：前記各補助画像信号を変換して前記倍率信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） のようになる。
号は、演算器23に入力される。演算器23は、積算手段４
および周波数強調処理手段５を内包するものである。演
算器23においては、複数の変換画像信号が積算され、こ
の積算により得られた積算値に原画像信号Ｓorgに応じ
て決定される強調係数βが乗じられ、さらにこの強調係
数βが乗じられた積算値が原画像信号Ｓorg と加算され
て処理済画像信号Ｓproc が得られる。この処理を、上
述した他の処理とともに数式として表すと、下記の式
（９） Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)・ｇ（Ｓus1−Ｓus2） ＋ｆ2(Ｓus1 −Ｓus2)・ｇ（Ｓus2−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk )・ｇ（Ｓusk−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（ＳusN−ＳusN+1）｝ …（９ ） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換して前記
抑制信号を作成するための関数 ｇ：前記各補助画像信号を変換して前記倍率信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） のようになる。
【００９０】図36は、上記実施の形態においてエッジ付
近の画像信号を処理した際の効果を示す図である。(1)
はエッジ付近における階段状の原画像信号と、その原画
像信号に基づいて作成されたボケ画像信号を示し、(2)
は、(1)の信号の帯域制限画像信号を示し、(3)は(2)の
帯域制限画像信号より１つ低い周波数帯域の信号を含む
補助画像信号を示し、(4)は(3)の補助画像信号を関数ｇ
により処理した際に得られる倍率信号を示し、(5)は(2)
の信号を関数ｆで処理した抑制画像信号に(4)の倍率信
号を乗じたものを示している。なお(3)に示されている
値Ｋは、図35に示した値Ｋであり、(3)の補助画像信号
がＫとなる時に(4)の倍率信号が０となることが図示さ
れている。(2)の信号を単に絶対値が小さくなるように
変換した場合には、信号のピークの形状は滑らかになる
ものの、立ち上がり部分は急峻なままとなる。これに対
し、(5)に示される処理済信号は立ち上がり部分が滑ら
かになっている。このように各帯域制限画像信号の立ち
上がり部分を滑らかにすることにより、これらの信号を
積算して得られる信号の周波数帯域の境界に階段状に生
じるアーチファクトを防止することができる。
近の画像信号を処理した際の効果を示す図である。(1)
はエッジ付近における階段状の原画像信号と、その原画
像信号に基づいて作成されたボケ画像信号を示し、(2)
は、(1)の信号の帯域制限画像信号を示し、(3)は(2)の
帯域制限画像信号より１つ低い周波数帯域の信号を含む
補助画像信号を示し、(4)は(3)の補助画像信号を関数ｇ
により処理した際に得られる倍率信号を示し、(5)は(2)
の信号を関数ｆで処理した抑制画像信号に(4)の倍率信
号を乗じたものを示している。なお(3)に示されている
値Ｋは、図35に示した値Ｋであり、(3)の補助画像信号
がＫとなる時に(4)の倍率信号が０となることが図示さ
れている。(2)の信号を単に絶対値が小さくなるように
変換した場合には、信号のピークの形状は滑らかになる
ものの、立ち上がり部分は急峻なままとなる。これに対
し、(5)に示される処理済信号は立ち上がり部分が滑ら
かになっている。このように各帯域制限画像信号の立ち
上がり部分を滑らかにすることにより、これらの信号を
積算して得られる信号の周波数帯域の境界に階段状に生
じるアーチファクトを防止することができる。
【００９１】なお、関数ｇにより処理される補助画像信
号は、必ずしも上記のようなものでなくてもよく、本発
明の他の実施の形態では、周波数強調処理装置では、下
記の式（10） Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)・ｇ（Ｓorg−Ｓus2） ＋ｆ2(Ｓus1 −Ｓus2)・ｇ（Ｓorg−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk )・ｇ（Ｓorg−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（Ｓorg−ＳusN+1）｝ …（10 ） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換して前記
抑制信号を作成するための関数 ｇ：前記各補助画像信号を変換して前記倍率信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） のような式にしたがって画像処理を行っている。
号は、必ずしも上記のようなものでなくてもよく、本発
明の他の実施の形態では、周波数強調処理装置では、下
記の式（10） Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)・ｇ（Ｓorg−Ｓus2） ＋ｆ2(Ｓus1 −Ｓus2)・ｇ（Ｓorg−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk )・ｇ（Ｓorg−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（Ｓorg−ＳusN+1）｝ …（10 ） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換して前記
抑制信号を作成するための関数 ｇ：前記各補助画像信号を変換して前記倍率信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） のような式にしたがって画像処理を行っている。
【００９２】図37は、この式（10）を適用した実施の形
態においてエッジ付近の画像信号を処理した際の効果を
示す図であり、図36と同様に、(1)はエッジ付近におけ
る階段状の原画像信号と、その原画像信号に基づいて作
成されたボケ画像信号を示し、(2)は、(1)の信号の帯域
制限画像信号を示し、(3)は(2)の帯域制限画像信号より
１つ低い周波数帯域の信号を含む補助画像信号を示し、
(4)は(3)の信号を関数ｇにより処理した際に得られる倍
率信号を示し、(5)は(2)の帯域制限画像信号を関数ｆで
処理した抑制画像信号に(4)の倍率信号を乗じたものを
示している。この図に示されるように、関数ｇにより処
理する値としてＳorg−Ｓuskを用いる場合、コントラス
トの大きなエッジについては、変換により得られる信号
は小さくなり、コントラストの小さなエッジについては
変換により得られる信号はもとの帯域制限画像信号に近
い大きさとなる。
態においてエッジ付近の画像信号を処理した際の効果を
示す図であり、図36と同様に、(1)はエッジ付近におけ
る階段状の原画像信号と、その原画像信号に基づいて作
成されたボケ画像信号を示し、(2)は、(1)の信号の帯域
制限画像信号を示し、(3)は(2)の帯域制限画像信号より
１つ低い周波数帯域の信号を含む補助画像信号を示し、
(4)は(3)の信号を関数ｇにより処理した際に得られる倍
率信号を示し、(5)は(2)の帯域制限画像信号を関数ｆで
処理した抑制画像信号に(4)の倍率信号を乗じたものを
示している。この図に示されるように、関数ｇにより処
理する値としてＳorg−Ｓuskを用いる場合、コントラス
トの大きなエッジについては、変換により得られる信号
は小さくなり、コントラストの小さなエッジについては
変換により得られる信号はもとの帯域制限画像信号に近
い大きさとなる。
【００９３】この信号は積算された後に、強調のために
原画像信号に加算される。この場合、コントラストの大
きなエッジはほとんど強調されないのに対し、コントラ
ストの小さなエッジはコントラストの大きなエッジに比
べると相対的に強く強調されるという効果が得られる。
原画像信号に加算される。この場合、コントラストの大
きなエッジはほとんど強調されないのに対し、コントラ
ストの小さなエッジはコントラストの大きなエッジに比
べると相対的に強く強調されるという効果が得られる。
【００９４】以上、２種類の式について説明したが、関
数ｆおよびｇについて、あるいは関数ｇにより処理する
帯域制限画像信号の作成方法については、様々な変更が
可能である。
数ｆおよびｇについて、あるいは関数ｇにより処理する
帯域制限画像信号の作成方法については、様々な変更が
可能である。
【００９５】次に、本発明の画像処理方法および装置の
さらに他の実施の形態について説明する。但し、この実
施の形態についても、ボケ画像信号作成手段１、帯域制
限画像信号作成手段２、積算手段４、周波数強調処理手
段５の処理は上記実施の形態と同じであるため説明を省
略し、変換手段３の処理についてのみ、以下説明する。
さらに他の実施の形態について説明する。但し、この実
施の形態についても、ボケ画像信号作成手段１、帯域制
限画像信号作成手段２、積算手段４、周波数強調処理手
段５の処理は上記実施の形態と同じであるため説明を省
略し、変換手段３の処理についてのみ、以下説明する。
【００９６】図38は本実施の形態における画像処理装置
の全体の構成を示す図である。上記実施の形態と同様、
原画像信号Ｓorgおよびボケ画像信号作成手段１におい
て作成された複数のボケ画像信号Ｓusk に基づいて帯域
制限画像信号が作成される。この帯域制限画像信号は減
算器21により、互いに隣接する周波数帯域の２つのボケ
画像信号Ｓusk （但し原画像信号ＳorgについてはＳorg
とＳus1 )の減算を行うことにより得られる。すなわ
ち、Ｓorg −Ｓus1 、Ｓus1 −Ｓus2 、…ＳusN-1 −Ｓ
usN を順次計算することにより、複数の帯域制限画像信
号が求められる。
の全体の構成を示す図である。上記実施の形態と同様、
原画像信号Ｓorgおよびボケ画像信号作成手段１におい
て作成された複数のボケ画像信号Ｓusk に基づいて帯域
制限画像信号が作成される。この帯域制限画像信号は減
算器21により、互いに隣接する周波数帯域の２つのボケ
画像信号Ｓusk （但し原画像信号ＳorgについてはＳorg
とＳus1 )の減算を行うことにより得られる。すなわ
ち、Ｓorg −Ｓus1 、Ｓus1 −Ｓus2 、…ＳusN-1 −Ｓ
usN を順次計算することにより、複数の帯域制限画像信
号が求められる。
【００９７】上記のようにして求められた帯域制限画像
信号は、変換手段３において変換される。この変換で
は、図38に示されるように、ある周波数帯域の帯域制限
画像信号に、その周波数帯域よりも１つ低い周波数帯域
の帯域制限画像信号を変換器24により関数ｇを用いて変
換したものを加えて、さらにそれを変換器22により関数
ｆkを用いて変換する。ここで、関数ｇにより変換を行
う変換器24が、前記補助画像信号作成手段に相当し、加
算器26が前記複合帯域制限画像信号作成手段に相当す
る。
信号は、変換手段３において変換される。この変換で
は、図38に示されるように、ある周波数帯域の帯域制限
画像信号に、その周波数帯域よりも１つ低い周波数帯域
の帯域制限画像信号を変換器24により関数ｇを用いて変
換したものを加えて、さらにそれを変換器22により関数
ｆkを用いて変換する。ここで、関数ｇにより変換を行
う変換器24が、前記補助画像信号作成手段に相当し、加
算器26が前記複合帯域制限画像信号作成手段に相当す
る。
【００９８】関数ｇとしては、例えば図39に示されるよ
うなものが用いられる。図39に示される関数は原点をと
おり原点における傾きはほぼ０であり、処理される値が
大きくなるにつれて傾きが漸増し、最終的に傾きがほぼ
１となるような関数である。すなわち、信号の大きさが
小さいときには強く抑制し、大きくなるにつれて抑制の
度合いを緩めるというものである。上述の処理におい
て、この関数の、原点から徐々に傾きが増加する部分
は、補助画像信号の立ち上がりの部分の波形に影響を与
える。すなわち、この関数により変換を行うことで、急
峻な立ち上がり部分を滑らかにすることができる。ここ
で、実際の帯域制限画像信号の大きさには限度があるた
め、関数ｇは原点付近で傾きが０から徐々に増加するよ
うな関数であれば、どのようなものでもよく例えば、図
40に示されるようなものであってもよい。また、本実施
の形態では、信号の増幅は行わないため、傾きは最大１
としているが、本発明の効果は関数の傾きが原点付近で
漸増することにより得られるものであるため、この条件
が満たされていれば、必ずしも傾きを最大１に限定する
必要はない。
うなものが用いられる。図39に示される関数は原点をと
おり原点における傾きはほぼ０であり、処理される値が
大きくなるにつれて傾きが漸増し、最終的に傾きがほぼ
１となるような関数である。すなわち、信号の大きさが
小さいときには強く抑制し、大きくなるにつれて抑制の
度合いを緩めるというものである。上述の処理におい
て、この関数の、原点から徐々に傾きが増加する部分
は、補助画像信号の立ち上がりの部分の波形に影響を与
える。すなわち、この関数により変換を行うことで、急
峻な立ち上がり部分を滑らかにすることができる。ここ
で、実際の帯域制限画像信号の大きさには限度があるた
め、関数ｇは原点付近で傾きが０から徐々に増加するよ
うな関数であれば、どのようなものでもよく例えば、図
40に示されるようなものであってもよい。また、本実施
の形態では、信号の増幅は行わないため、傾きは最大１
としているが、本発明の効果は関数の傾きが原点付近で
漸増することにより得られるものであるため、この条件
が満たされていれば、必ずしも傾きを最大１に限定する
必要はない。
【００９９】関数f1〜fNは、全て同じ関数であってもよ
いし、異なる関数であってもよく、画像処理の目的に応
じて任意に設定可能なものとする。本実施の形態では、
この関数f1〜fNとして、例えば図24のような関数を使用
する。
いし、異なる関数であってもよく、画像処理の目的に応
じて任意に設定可能なものとする。本実施の形態では、
この関数f1〜fNとして、例えば図24のような関数を使用
する。
【０１００】上記変換手段３により得られた変換画像信
号は、演算器23に入力される。演算器23は、積算手段４
および周波数強調処理手段５を内包するものである。演
算器23においては、複数の変換画像信号が積算され、こ
の積算により得られた積算値に原画像信号Ｓorgに応じ
て決定される強調係数βが乗じられ、さらにこの強調係
数βが乗じられた積算値が原画像信号Ｓorg と加算され
て処理済画像信号Ｓproc が得られる。この処理を、上
述した他の処理とともに数式として表すと、下記の式
（11） Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝［ｆ1｛（Ｓorg −Ｓus1)＋ｇ（Ｓus1−Ｓus2）｝ ＋ｆ2｛（Ｓus1 −Ｓus2)＋ｇ（Ｓus2−Ｓus3）｝＋… ＋ｆk｛（Ｓusk-1−Ｓusk )＋ｇ（Ｓusk−Ｓusk+1）｝＋… ＋ｆN｛（ＳusN-1−ＳusN)＋ｇ（ＳusN−ＳusN+1）｝］ …（11） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各複合帯域制限画像信号を変換する
ための関数 ｇ：前記各帯域制限画像信号を変換して前記補助画像信
号を作成するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） のようになる。
号は、演算器23に入力される。演算器23は、積算手段４
および周波数強調処理手段５を内包するものである。演
算器23においては、複数の変換画像信号が積算され、こ
の積算により得られた積算値に原画像信号Ｓorgに応じ
て決定される強調係数βが乗じられ、さらにこの強調係
数βが乗じられた積算値が原画像信号Ｓorg と加算され
て処理済画像信号Ｓproc が得られる。この処理を、上
述した他の処理とともに数式として表すと、下記の式
（11） Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝［ｆ1｛（Ｓorg −Ｓus1)＋ｇ（Ｓus1−Ｓus2）｝ ＋ｆ2｛（Ｓus1 −Ｓus2)＋ｇ（Ｓus2−Ｓus3）｝＋… ＋ｆk｛（Ｓusk-1−Ｓusk )＋ｇ（Ｓusk−Ｓusk+1）｝＋… ＋ｆN｛（ＳusN-1−ＳusN)＋ｇ（ＳusN−ＳusN+1）｝］ …（11） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各複合帯域制限画像信号を変換する
ための関数 ｇ：前記各帯域制限画像信号を変換して前記補助画像信
号を作成するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） のようになる。
【０１０１】図41は、この実施の形態においてエッジ付
近の画像信号を処理した際の効果を示す図である。(1)
はエッジ付近における階段状の原画像信号と、その原画
像信号に基づいて作成されたボケ画像信号を示し、(2)
および(3)は、(1)の信号にかかる帯域制限画像信号であ
って、(2)を被変換帯域制限画像信号とした場合の低周
波側帯域制限画像信号が(3)である。(4)は(3)の低周波
側帯域制限画像信号を関数ｇにより処理した際に得られ
る補助画像信号を示し、(5)は(2)の被変換帯域制限画像
信号と(4)の補助画像信号を加算した複合帯域制限画像
信号を示し、(6)は(5)の複合帯域制限画像信号を関数ｆ
により処理して得られる変換画像信号を示している。こ
の図において、(2)の信号を単に絶対値が小さくなるよ
うに変換した場合には、信号のピークの形状は滑らかに
なるものの、立ち上がり部分は急峻なままとなる。これ
に対し、(5)に示される処理済信号は立ち上がり部分が
滑らかになっている。このように各帯域制限画像信号の
立ち上がり部分を滑らかにすることにより、これらの信
号を積算して得られる信号の周波数帯域の境界に階段状
に生じるアーチファクトを防止することができる。
近の画像信号を処理した際の効果を示す図である。(1)
はエッジ付近における階段状の原画像信号と、その原画
像信号に基づいて作成されたボケ画像信号を示し、(2)
および(3)は、(1)の信号にかかる帯域制限画像信号であ
って、(2)を被変換帯域制限画像信号とした場合の低周
波側帯域制限画像信号が(3)である。(4)は(3)の低周波
側帯域制限画像信号を関数ｇにより処理した際に得られ
る補助画像信号を示し、(5)は(2)の被変換帯域制限画像
信号と(4)の補助画像信号を加算した複合帯域制限画像
信号を示し、(6)は(5)の複合帯域制限画像信号を関数ｆ
により処理して得られる変換画像信号を示している。こ
の図において、(2)の信号を単に絶対値が小さくなるよ
うに変換した場合には、信号のピークの形状は滑らかに
なるものの、立ち上がり部分は急峻なままとなる。これ
に対し、(5)に示される処理済信号は立ち上がり部分が
滑らかになっている。このように各帯域制限画像信号の
立ち上がり部分を滑らかにすることにより、これらの信
号を積算して得られる信号の周波数帯域の境界に階段状
に生じるアーチファクトを防止することができる。
【０１０２】以上、２つの式ついて説明したが、この実
施の形態においても、例えば、関数ｆおよびｇについ
て、様々な変更が可能である。
施の形態においても、例えば、関数ｆおよびｇについ
て、様々な変更が可能である。
【０１０３】なお、本明細書において説明したいくつか
の実施の形態は全て、ガウス信号のフィルタを用いて、
補間演算処理を行うようにしているが、Ｂスプライン補
間演算によりフィルタリング処理画像に対して補間演算
処理を施すようにしてもよい。以下このＢスプライン補
間演算処理について説明する。
の実施の形態は全て、ガウス信号のフィルタを用いて、
補間演算処理を行うようにしているが、Ｂスプライン補
間演算によりフィルタリング処理画像に対して補間演算
処理を施すようにしてもよい。以下このＢスプライン補
間演算処理について説明する。
【０１０４】Ｂスプライン補間演算は、比較的鮮鋭度の
低い滑らかな２次画像を再生するための補間画像データ
を得るための補間演算方法である。このＢスプライン補
間演算は、元のサンプル点（画素）を通ることは必要と
されない代わりに、第１階微分係数および第２階微分係
数（ｆ″（Ｘ）で表す）が各区間間で連続することが必
要とされる。
低い滑らかな２次画像を再生するための補間画像データ
を得るための補間演算方法である。このＢスプライン補
間演算は、元のサンプル点（画素）を通ることは必要と
されない代わりに、第１階微分係数および第２階微分係
数（ｆ″（Ｘ）で表す）が各区間間で連続することが必
要とされる。
【０１０５】すなわち、 ｆk （ｘ）＝Ａk ｘ3 ＋Ｂk ｘ2 ＋Ｃk ｘ＋Ｄk …（12） において（式（12）においてＢk は便宜上用いられる係
数であり、フィルタリング処理画像とは異なる。）、 ｆk ′（Ｘk ）＝ｆk-1 ′（Ｘk ） …（13） ｆk ′（Ｘk+1 ）＝ｆk+1 ′（Ｘk+1 ） …（14） ｆk ″（Ｘk ）＝ｆk-1 ″（Ｘk ） …（15） ｆk ″（Ｘk+1 ）＝ｆk+1 ″（Ｘk+1 ） …（16） が条件となる。ただし、画素Ｘk における第１階微分係
数が、その画素Ｘk の前後の画素であるＸk-1 とＸk+1
とについて、これらの画像信号Ｙk-1 、Ｙk+1 の勾配
（Ｙk+1 −Ｙk-1 ）／（Ｘk+1 −Ｘk-1 ）に一致するこ
とが条件であるから、下記式（17）を満たす必要があ
る。
数であり、フィルタリング処理画像とは異なる。）、 ｆk ′（Ｘk ）＝ｆk-1 ′（Ｘk ） …（13） ｆk ′（Ｘk+1 ）＝ｆk+1 ′（Ｘk+1 ） …（14） ｆk ″（Ｘk ）＝ｆk-1 ″（Ｘk ） …（15） ｆk ″（Ｘk+1 ）＝ｆk+1 ″（Ｘk+1 ） …（16） が条件となる。ただし、画素Ｘk における第１階微分係
数が、その画素Ｘk の前後の画素であるＸk-1 とＸk+1
とについて、これらの画像信号Ｙk-1 、Ｙk+1 の勾配
（Ｙk+1 −Ｙk-1 ）／（Ｘk+1 −Ｘk-1 ）に一致するこ
とが条件であるから、下記式（17）を満たす必要があ
る。
【０１０６】 ｆk ′（Ｘk ）＝（Ｙk+1 −Ｙk-1 ）／（Ｘk+1 −Ｘk-1 ） …（17） 同様に、画素Ｘk+1 における第１階微分係数が、その画
素Ｘk+1 の前後の画素であるＸk とＸk+2 とについて、
これらの画像信号Ｙk 、Ｙk+2 の勾配（Ｙk+2−Ｙk ）
／（Ｘk+2 −Ｘk ）に一致することが条件であるから、
下記式（18）を満たす必要がある。
素Ｘk+1 の前後の画素であるＸk とＸk+2 とについて、
これらの画像信号Ｙk 、Ｙk+2 の勾配（Ｙk+2−Ｙk ）
／（Ｘk+2 −Ｘk ）に一致することが条件であるから、
下記式（18）を満たす必要がある。
【０１０７】 ｆk ′（Ｘk+1 ）＝（Ｙk+2 −Ｙk ）／（Ｘk+2 −Ｘk ） …（18） また関数ｆ（Ｘ）は一般に下記式（19）に示すもので近
似される。
似される。
【０１０８】 ｆ（Ｘ）＝ｆ(0) ＋ｆ′(0) Ｘ＋｛ｆ″(0)/２｝Ｘ2 …（19） ここで、各区間Ｘk-2 〜Ｘk-1 ，Ｘk-1 〜Ｘk ，Ｘk 〜
Ｘk+1 ，Ｘk+1 〜Ｘk+2 の間隔（格子間隔という）を１
とし、画素Ｘk からの画素Ｘk+1 方向への補間点Ｘp の
位置をｔ（０≦ｔ≦１）とすれば、式（13）〜（16）お
よび（19）より、 ｆk ′（０）＝Ｃk ＝（Ｙk+1 −Ｙk-1 ）／２ ｆk ′（１）＝３Ａk ＋２Ｂk ＋Ｃk ＝（Ｙk+2 −
Ｙk ）／２ ｆk ″（０）＝Ｙk+1 −２Ｙk ＋Ｙk-1 ＝２Ｂ したがって、 Ａk ＝（Ｙk+2 −３Ｙk+1 ＋３Ｙk −Ｙk-1 ）／６ Ｂk ＝（Ｙk+1 −２Ｙk ＋Ｙk-1 ）／２ Ｃk ＝（Ｙk+1 −Ｙk-1 ）／２ ここで、Ｄk は未知のため、 Ｄk ＝（Ｄ1 Ｙk+2 ＋Ｄ2 Ｙk+1 ＋Ｄ3 Ｙk ＋Ｄ4 Ｙ
k-1 ）／６ とおく。また、スプライン補間関数ｆk （ｘ）は上述の
通り、Ｘ＝ｔなる変数変換をしているため、 ｆk （ｘ）＝ｆk （ｔ） となる。よって、 ｆk （ｔ）＝｛（Ｙk+2 −３Ｙk+1 ＋３Ｙk −Ｙk-1 ）
／６｝ｔ3 ＋｛（Ｙk+1 −２Ｙk ＋Ｙk-1 ）／２｝ｔ2
＋｛（Ｙk+1 −Ｙk-1 ）／２｝ｔ＋（Ｄ1 Ｙk+2 ＋Ｄ2
Ｙk+1 ＋Ｄ3 Ｙk ＋Ｄ4 Ｙk-1 ）／６ となり、これを画像信号Ｙk-1 ，Ｙk ，Ｙk+1 ，Ｙk+2
について整理すると、下記式（20）で表すことができ
る。
Ｘk+1 ，Ｘk+1 〜Ｘk+2 の間隔（格子間隔という）を１
とし、画素Ｘk からの画素Ｘk+1 方向への補間点Ｘp の
位置をｔ（０≦ｔ≦１）とすれば、式（13）〜（16）お
よび（19）より、 ｆk ′（０）＝Ｃk ＝（Ｙk+1 −Ｙk-1 ）／２ ｆk ′（１）＝３Ａk ＋２Ｂk ＋Ｃk ＝（Ｙk+2 −
Ｙk ）／２ ｆk ″（０）＝Ｙk+1 −２Ｙk ＋Ｙk-1 ＝２Ｂ したがって、 Ａk ＝（Ｙk+2 −３Ｙk+1 ＋３Ｙk −Ｙk-1 ）／６ Ｂk ＝（Ｙk+1 −２Ｙk ＋Ｙk-1 ）／２ Ｃk ＝（Ｙk+1 −Ｙk-1 ）／２ ここで、Ｄk は未知のため、 Ｄk ＝（Ｄ1 Ｙk+2 ＋Ｄ2 Ｙk+1 ＋Ｄ3 Ｙk ＋Ｄ4 Ｙ
k-1 ）／６ とおく。また、スプライン補間関数ｆk （ｘ）は上述の
通り、Ｘ＝ｔなる変数変換をしているため、 ｆk （ｘ）＝ｆk （ｔ） となる。よって、 ｆk （ｔ）＝｛（Ｙk+2 −３Ｙk+1 ＋３Ｙk −Ｙk-1 ）
／６｝ｔ3 ＋｛（Ｙk+1 −２Ｙk ＋Ｙk-1 ）／２｝ｔ2
＋｛（Ｙk+1 −Ｙk-1 ）／２｝ｔ＋（Ｄ1 Ｙk+2 ＋Ｄ2
Ｙk+1 ＋Ｄ3 Ｙk ＋Ｄ4 Ｙk-1 ）／６ となり、これを画像信号Ｙk-1 ，Ｙk ，Ｙk+1 ，Ｙk+2
について整理すると、下記式（20）で表すことができ
る。
【０１０９】 ｆk （ｔ）＝｛（−ｔ3 ＋３ｔ2 −３ｔ＋Ｄ4 ）／６｝Ｙk-1 ＋｛（３ｔ3 −６ｔ2 ＋Ｄ3 ）／６｝Ｙk ＋｛（−３ｔ3 ＋３ｔ2 ＋３ｔ＋Ｄ2 ）／６｝Ｙk+1 ＋｛（ｔ3 ＋Ｄ1 ）／６｝Ｙk+2 …（20） ここで、ｔ＝１とおけば、 ｆk （１）＝｛（Ｄ4 −１）／６｝Ｙk-1 ＋｛（Ｄ3 −
３）／６｝Ｙk ＋｛（Ｄ2 ＋３）／６｝Ｙk+1 ＋｛（Ｄ
1 ＋１）／６｝Ｙk+2 次に区間Ｘk+1 〜Ｘk+2 についての式（17）は、 ｆk+1 （ｔ）＝｛（−ｔ3 ＋３ｔ2 −３ｔ＋Ｄ4 ）／６｝Ｙk ＋｛（３ｔ3 −６ｔ2 ＋Ｄ3 ）／６｝Ｙk+1 ＋｛（−３ｔ3 ＋３ｔ2 ＋３ｔ＋Ｄ2 ）／６｝Ｙk+2 ＋｛（ｔ3 ＋Ｄ1 ）／６｝Ｙk+3 …（21） ここで、ｔ＝０とおけば、 ｆk+1 （０）＝（Ｄ4 ／６）Ｙk ＋（Ｄ3 ／６）Ｙk+1
＋（Ｄ2 ／６）Ｙk+2 ＋（Ｄ1 ／６）Ｙk+3 連続性の条件（ｆk （１）＝ｆk+1 （０））、および各
フィルタリング処理画像信号に対応する係数同士が等し
いという条件により、Ｄ4 −１＝０，Ｄ3 −３＝Ｄ4 ，
Ｄ2 ＋３＝Ｄ3 ，Ｄ1 ＋１＝Ｄ2 ，Ｄ1 ＝０、となり、
したがって、 Ｄk ＝（Ｙk+1 ＋４Ｙk ＋Ｙk-1 ）／６ となる。よって、 Ｙp ＝ｆk （ｔ）＝｛（−ｔ3 ＋３ｔ2 −３ｔ＋１）／６｝Ｙk-1 ＋｛（３ｔ3 −６ｔ2 ＋４）／６｝Ｙk ＋｛（−３ｔ3 ＋３ｔ2 ＋３ｔ＋１）／６｝Ｙk+1 ＋｛ｔ3 ／６｝Ｙk+2 …（22） したがって、フィルタリング処理画像信号Ｙk-1 、
Ｙk 、Ｙk+1 、Ｙk+2 にそれぞれ対応する補間係数ｂ
k-1 、ｂk 、ｂk+1 、ｂk+2 は、 ｂk-1 ＝（−ｔ3 ＋３ｔ2 −３ｔ＋１）／６ ｂk ＝（３ｔ3 −６ｔ2 ＋４）／６ ｂk+1 ＝（−３ｔ3 ＋３ｔ2 ＋３ｔ＋１）／６ ｂk+2 ＝ｔ3 ／６ となる。
３）／６｝Ｙk ＋｛（Ｄ2 ＋３）／６｝Ｙk+1 ＋｛（Ｄ
1 ＋１）／６｝Ｙk+2 次に区間Ｘk+1 〜Ｘk+2 についての式（17）は、 ｆk+1 （ｔ）＝｛（−ｔ3 ＋３ｔ2 −３ｔ＋Ｄ4 ）／６｝Ｙk ＋｛（３ｔ3 −６ｔ2 ＋Ｄ3 ）／６｝Ｙk+1 ＋｛（−３ｔ3 ＋３ｔ2 ＋３ｔ＋Ｄ2 ）／６｝Ｙk+2 ＋｛（ｔ3 ＋Ｄ1 ）／６｝Ｙk+3 …（21） ここで、ｔ＝０とおけば、 ｆk+1 （０）＝（Ｄ4 ／６）Ｙk ＋（Ｄ3 ／６）Ｙk+1
＋（Ｄ2 ／６）Ｙk+2 ＋（Ｄ1 ／６）Ｙk+3 連続性の条件（ｆk （１）＝ｆk+1 （０））、および各
フィルタリング処理画像信号に対応する係数同士が等し
いという条件により、Ｄ4 −１＝０，Ｄ3 −３＝Ｄ4 ，
Ｄ2 ＋３＝Ｄ3 ，Ｄ1 ＋１＝Ｄ2 ，Ｄ1 ＝０、となり、
したがって、 Ｄk ＝（Ｙk+1 ＋４Ｙk ＋Ｙk-1 ）／６ となる。よって、 Ｙp ＝ｆk （ｔ）＝｛（−ｔ3 ＋３ｔ2 −３ｔ＋１）／６｝Ｙk-1 ＋｛（３ｔ3 −６ｔ2 ＋４）／６｝Ｙk ＋｛（−３ｔ3 ＋３ｔ2 ＋３ｔ＋１）／６｝Ｙk+1 ＋｛ｔ3 ／６｝Ｙk+2 …（22） したがって、フィルタリング処理画像信号Ｙk-1 、
Ｙk 、Ｙk+1 、Ｙk+2 にそれぞれ対応する補間係数ｂ
k-1 、ｂk 、ｂk+1 、ｂk+2 は、 ｂk-1 ＝（−ｔ3 ＋３ｔ2 −３ｔ＋１）／６ ｂk ＝（３ｔ3 −６ｔ2 ＋４）／６ ｂk+1 ＝（−３ｔ3 ＋３ｔ2 ＋３ｔ＋１）／６ ｂk+2 ＝ｔ3 ／６ となる。
【０１１０】以上の演算を各区間Ｘk-2 〜Ｘk-1 ，Ｘ
k-1 〜Ｘk ，Ｘk 〜Ｘk+1 ，Ｘk+1 〜Ｘk+2 について繰
り返すことにより、フィルタリング処理画像信号の全体
についてフィルタリング処理画像信号とは間隔の異なる
補間画像信号を求めることができる。
k-1 〜Ｘk ，Ｘk 〜Ｘk+1 ，Ｘk+1 〜Ｘk+2 について繰
り返すことにより、フィルタリング処理画像信号の全体
についてフィルタリング処理画像信号とは間隔の異なる
補間画像信号を求めることができる。
【０１１１】したがって、このＢスプライン補間演算処
理を各フィルタリング処理画像信号Ｂk に対して施すこ
とにより、各フィルタリング処理画像信号Ｂk に対応し
たボケ画像信号Ｓusk を得ることができる。
理を各フィルタリング処理画像信号Ｂk に対して施すこ
とにより、各フィルタリング処理画像信号Ｂk に対応し
たボケ画像信号Ｓusk を得ることができる。
【０１１２】また、上記全ての実施の形態では、原画像
信号Ｓorg に対して１画素おきにフィルタリング処理を
施し、さらに補間演算を施すことによりボケ画像信号Ｓ
uskを得るようにしているが、フィルタリング処理につ
いてもこれに限定されるものではなく、複数のサイズの
ボケマスクを用いて、原画像信号Ｓorg にフィルタリン
グ処理を施すことにより、複数の周波数応答特性が異な
るボケ画像信号Ｓuskを得るようにしてもよい。
信号Ｓorg に対して１画素おきにフィルタリング処理を
施し、さらに補間演算を施すことによりボケ画像信号Ｓ
uskを得るようにしているが、フィルタリング処理につ
いてもこれに限定されるものではなく、複数のサイズの
ボケマスクを用いて、原画像信号Ｓorg にフィルタリン
グ処理を施すことにより、複数の周波数応答特性が異な
るボケ画像信号Ｓuskを得るようにしてもよい。
【０１１３】さらに、上述した全ての実施の形態におけ
る画像処理方法および装置は、断層撮影により画像に現
れる流れ像を除去する場合にも適用できる。
る画像処理方法および装置は、断層撮影により画像に現
れる流れ像を除去する場合にも適用できる。
【０１１４】すなわち、断層撮影を行うことにより得ら
れる流れ像を除去するために、ボケマスク処理を行う方
法が提案されているが（特開平3-276265号）、この場合
においても、急激に濃度が変化するエッジ部の近傍にお
いて、アーチファクトが発生するという問題がある。こ
の場合、流れ像を除去するためにはボケマスクを１次元
としてボケ画像信号Ｓusk を得てボケマスク処理を行う
ものであるが、上述した実施の形態における画像処理に
おける互いに周波数応答特性が異なるボケ画像信号Ｓus
k を得る際に、流れ像を除去するための１次元ボケマス
クによるフィルタリング処理を行い、これに基づいて上
述したような帯域制限画像信号を得、この帯域制限画像
信号の絶対値が比較的大きい場合は、その絶対値が小さ
くなるように変換して、その後の処理を行うことによ
り、流れ像を除去するとともに、アーチファクトのない
良好な画像を得ることができる。
れる流れ像を除去するために、ボケマスク処理を行う方
法が提案されているが（特開平3-276265号）、この場合
においても、急激に濃度が変化するエッジ部の近傍にお
いて、アーチファクトが発生するという問題がある。こ
の場合、流れ像を除去するためにはボケマスクを１次元
としてボケ画像信号Ｓusk を得てボケマスク処理を行う
ものであるが、上述した実施の形態における画像処理に
おける互いに周波数応答特性が異なるボケ画像信号Ｓus
k を得る際に、流れ像を除去するための１次元ボケマス
クによるフィルタリング処理を行い、これに基づいて上
述したような帯域制限画像信号を得、この帯域制限画像
信号の絶対値が比較的大きい場合は、その絶対値が小さ
くなるように変換して、その後の処理を行うことによ
り、流れ像を除去するとともに、アーチファクトのない
良好な画像を得ることができる。
【０１１５】また、流れ像を除去するための１次元ボケ
マスクによる処理とともに、上述した２次元のボケマス
ク、あるいは１次元ボケマスクを２次元的の施す処理を
同時に行うようにしてもよい。これにより、原画像から
流れ像が除去されるとともに、エッジ部近傍のアーチフ
ァクトをより一層低減することができるため、より画質
のよい処理済画像を得ることができる。
マスクによる処理とともに、上述した２次元のボケマス
ク、あるいは１次元ボケマスクを２次元的の施す処理を
同時に行うようにしてもよい。これにより、原画像から
流れ像が除去されるとともに、エッジ部近傍のアーチフ
ァクトをより一層低減することができるため、より画質
のよい処理済画像を得ることができる。
【０１１６】また、上記一実施の形態において、帯域制
限画像信号に対して撮影部位ごとに異なる変換処理を施
すことにより得られる効果について説明したが、この撮
影部位ごとに処理を変えるという思想は上記全ての実施
の形態に適用できるものである。
限画像信号に対して撮影部位ごとに異なる変換処理を施
すことにより得られる効果について説明したが、この撮
影部位ごとに処理を変えるという思想は上記全ての実施
の形態に適用できるものである。
【図１】本発明による画像処理装置の概念を表す図
【図２】ボケ画像信号作成手段の詳細を表す図
【図３】ボケ画像信号作成手段において用いられるフィ
ルタ（１次元）を表す図
ルタ（１次元）を表す図
【図４】フィルタリング処理の詳細を表す図
【図５】フィルタリング処理画像信号の周波数特性を表
す図
す図
【図６】フィルタリング処理手段において用いられるフ
ィルタ（２次元）を表す図
ィルタ（２次元）を表す図
【図７】フィルタリング処理画像信号Ｂ1 の補間演算に
用いられるフィルタを表す図
用いられるフィルタを表す図
【図８】補間演算の詳細を表す図
【図９】フィルタリング処理画像信号Ｂ2 の補間演算に
用いられるフィルタを表す図
用いられるフィルタを表す図
【図１０】補間演算の詳細を表す図
【図１１】ボケ画像信号の周波数特性を表す図
【図１２】本発明の画像処理装置の一実施の形態におけ
る全体構成を示す図
る全体構成を示す図
【図１３】帯域制限画像信号の周波数特性を表す図
【図１４】変換手段における帯域制限画像信号の変換処
理の一例を表す図
理の一例を表す図
【図１５】処理済画像信号の周波数特性を表す図
【図１６】処理済画像信号の他の周波数特性を表す図
【図１７】変換手段における帯域制限画像信号の変換処
理の他の例を表す図
理の他の例を表す図
【図１８】帯域制限画像信号の他の周波数特性を表す図
【図１９】撮影部位に応じた帯域制限画像信号の変換処
理を表す図
理を表す図
【図２０】周波数帯域に応じた帯域制限画像信号の変換
処理を表す図
処理を表す図
【図２１】周波数帯域に応じた帯域制限画像信号の他の
変換処理を表す図
変換処理を表す図
【図２２】図２０に示す関数により得られる処理済画像
信号の周波数特性を表す図
信号の周波数特性を表す図
【図２３】図２１に示す関数により得られる処理済画像
信号の周波数特性を表す図
信号の周波数特性を表す図
【図２４】変換手段において変換処理に使用される関数
の一例を示す図
の一例を示す図
【図２５】変換手段において変換処理に使用される関数
の他の一例を示す図
の他の一例を示す図
【図２６】図２４の関数と図２５の関数の特徴を兼ね備
えた関数の一例を示す図
えた関数の一例を示す図
【図２７】周波数帯域ごとの周波数特性と全体の周波数
特性とを示す図
特性とを示す図
【図２８】全ての帯域制限画像信号を１種類の関数で変
換した際の問題点を示す図
換した際の問題点を示す図
【図２９】理想のエッジと実際のエッジを示す図
【図３０】強調係数が小さい場合に各帯域制限画像信号
を変換する関数の一例を示す図
を変換する関数の一例を示す図
【図３１】強調係数が大きい場合に各帯域制限画像信号
を変換する関数の一例を示す図
を変換する関数の一例を示す図
【図３２】図３０の関数を変換に用いることにより得ら
れる画像信号のレスポンスを示す図
れる画像信号のレスポンスを示す図
【図３３】図３１の関数を変換に用いることにより得ら
れる画像信号のレスポンスを示す図
れる画像信号のレスポンスを示す図
【図３４】本発明の画像処理装置の他の実施の形態にお
ける全体構成を示す図
ける全体構成を示す図
【図３５】倍率信号作成手段において使用される関数の
一例を示す図
一例を示す図
【図３６】図３４に示す実施の形態における効果を示す
図
図
【図３７】図３４に示す実施の形態における効果を示す
他の図
他の図
【図３８】本発明の画像処理装置のさらに他の実施の形
態における全体構成を示す図
態における全体構成を示す図
【図３９】各帯域制限画像信号を変換して補助画像信号
を作成するための関数の一例を示す図
を作成するための関数の一例を示す図
【図４０】各帯域制限画像信号を変換して補助画像信号
を作成するための関数の他の例を示す図
を作成するための関数の他の例を示す図
【図４１】図３８に示す実施の形態における効果を示す
図
図
【図４２】エッジ部近傍のボケ画像信号の作成処理を示
す図
す図
１ ボケ画像信号作成手段 ２ 帯域制限画像信号作成手段 ３ 変換手段 ４ 積算手段 ５ 周波数強調処理手段 10 フィルタリング処理手段 11 補間演算処理手段 21 減算器 22 変換器 23 演算器 24 変換器 25 乗算手段 26 加算器
───────────────────────────────────────────────────── フロントページの続き (31)優先権主張番号 特願平8−169063 (32)優先日 平８(1996)６月28日 (33)優先権主張国 日本（ＪＰ） (31)優先権主張番号 特願平8−169064 (32)優先日 平８(1996)６月28日 (33)優先権主張国 日本（ＪＰ） (31)優先権主張番号 特願平8−169065 (32)優先日 平８(1996)６月28日 (33)優先権主張国 日本（ＪＰ）
Claims (40)
- 【請求項１】 原画像を表す原画像信号に、該原画像の
高周波成分に関する信号を加算することにより、該原画
像の高周波成分を強調する画像処理方法において、 前記原画像信号に基づいて、互いに周波数応答特性が異
なる複数の非鮮鋭マスク画像信号を作成し、 前記原画像信号および前記複数の非鮮鋭マスク画像信
号、または前記複数の非鮮鋭マスク画像信号に基づい
て、前記原画像信号の複数の周波数帯域ごとの信号を表
す複数の帯域制限画像信号を作成し、 該各帯域制限画像信号のうち少なくとも１つの信号につ
いて、該帯域制限画像信号の少なくとも一部を小さくす
るような変換処理を施して複数の変換画像信号を作成
し、 該各変換画像信号を積算することにより前記原画像信号
に加算する前記高周波成分に関する信号を得ることを特
徴とする画像処理方法。 - 【請求項２】 前記複数の変換画像信号の作成を、前記
帯域制限画像信号の絶対値が所定の閾値より大きい帯域
制限画像信号について、該帯域制限画像信号の絶対値を
小さくするように、該各帯域制限画像信号を変換するこ
とにより行うことを特徴とする請求項１記載の画像処理
方法。 - 【請求項３】 前記複数の変換画像信号の作成を、前記
各帯域制限画像信号が所定の閾値より大きいほど該帯域
制限画像信号の絶対値を小さくするとともに、前記各帯
域制限画像信号の絶対値が該所定の閾値よりも小さい他
の閾値よりも小さいほど該帯域制限画像信号の絶対値を
小さくするように、該各帯域制限画像信号を変換するこ
とにより行うことを特徴とする請求項２記載の画像処理
方法。 - 【請求項４】 前記複数の変換画像信号の作成を、前記
各帯域制限画像信号の周波数帯域に応じて該各帯域制限
画像信号の絶対値の大きさを変化させるように、該各帯
域制限画像信号を変換することにより行うことを特徴と
する請求項２または３記載の画像処理方法。 - 【請求項５】 前記複数の変換画像信号の作成を、前記
各帯域制限画像信号を、前記周波数帯域ごとに異なる複
数の関数に基づいて、該帯域制限画像信号の絶対値に基
づいて決まる該絶対値以下の値となるように変換するこ
とにより行うことを特徴とする請求項１記載の画像処理
方法。 - 【請求項６】 前記変換画像信号の作成に使用される各
関数が、前記帯域制限画像信号の絶対値が所定の値より
大きい帯域制限画像信号については、該変換画像信号の
値がほぼ一定の値となるように該帯域制限画像信号を変
換する関数であり、該関数が高周波帯域を処理する関数
であるほど、前記所定の値が小さい値であることを特徴
とする請求項５記載の画像処理方法。 - 【請求項７】 前記変換画像信号の作成に使用される各
関数が、該関数が低周波帯域を処理する関数であるほ
ど、前記帯域制限画像信号の絶対値が０近傍の所定の範
囲内の値である帯域制限画像信号を変換した際に得られ
る該変換画像信号の絶対値が小さい値であることを特徴
とする請求項５または６記載の画像処理方法。 - 【請求項８】 前記帯域制限画像信号の作成、前記変換
画像信号の作成、前記高周波成分に関する信号の作成、
および該高周波成分に関する信号の前記原画像信号への
加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)＋ｆ2(Ｓus1 −Ｓus2)＋… ＋ｆk(Ｓusk-1−Ｓusk )＋…＋ｆN(ＳusN-1−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項１から７のい
ずれか１項記載の画像処理方法。 - 【請求項９】 前記帯域制限画像信号の作成、前記変換
画像信号の作成、前記高周波成分に関する信号の作成、
および該高周波成分に関する信号の前記原画像信号への
加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN ） ＝（1/N)・｛ｆ1（Ｓorg −Ｓus1)＋ｆ2（Ｓorg −Ｓus2)＋… ＋fk(Ｓorg −Ｓusk)＋…＋ｆN（Ｓorg −ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項１から７のい
ずれか１項記載の画像処理方法。 - 【請求項１０】 前記帯域制限画像信号の作成、前記変
換画像信号の作成、前記高周波成分に関する信号の作
成、および該高周波成分に関する信号の前記原画像信号
への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓus1,Ｓus2,…ＳusN ） ＝｛ｆ2（Ｓus1 −Ｓus2)＋ｆ3（Ｓus2 −Ｓus3)＋… ＋fk(Ｓusk-1 −Ｓusk)＋…＋ｆN（ＳusN-1 −ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=2〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項１から４のい
ずれか１項記載の画像処理方法。 - 【請求項１１】 前記帯域制限画像信号の作成、前記変
換画像信号の作成、前記高周波成分に関する信号の作
成、および該高周波成分に関する信号の前記原画像信号
への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓus1,Ｓus2,…ＳusN ） ＝（1/N)・｛ｆ2（Ｓus1 −Ｓus2)＋ｆ3（Ｓus1 −Ｓus3)＋… ＋fk(Ｓus1 −Ｓusk)＋…＋ｆN（Ｓus1 −ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=2〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項１から４のい
ずれか１項記載の画像処理方法。 - 【請求項１２】 前記複数の変換画像信号の作成を、前
記各帯域制限画像信号に対して、前記強調係数に応じた
変換処理を施すことにより行うことを特徴とする請求項
１から１１のいずれか１項記載の画像処理方法。 - 【請求項１３】 前記変換処理が、前記各周波数帯域ご
とに異なる変換処理であって、前記帯域制限画像信号の
周波数帯域が低周波帯域であるほど、該帯域制限画像信
号を強く抑制するような変換処理であり、 前記強調係数が大きいときほど、高周波帯域の帯域制限
画像信号に対する抑制の度合いと、低周波帯域の帯域制
限画像信号に対する抑制の度合いの差が大きいことを特
徴とする請求項１２記載の画像処理方法。 - 【請求項１４】 前記複数の変換画像信号の作成を、前
記帯域制限画像信号を、該帯域制限画像信号の絶対値に
基づいて決まる該絶対値以下の値となるように変換する
ことにより抑制画像信号を作成し、 前記原画像信号および前記複数の非鮮鋭マスク画像信号
に基づいて、前記抑制画像信号の作成に使用された帯域
制限画像信号よりも低い周波数帯域の信号を含む補助画
像信号を作成し、 該補助画像信号を、該補助画像信号の絶対値が小さいほ
ど１に近く、大きいほど０に近い値となるように変換す
ることにより前記抑制画像信号にそれぞれ対応する倍率
信号を作成し、 前記抑制画像信号に、該抑制画像信号に対応する前記倍
率信号を乗ずることにより行うことを特徴とする請求項
１記載の画像処理方法。 - 【請求項１５】 前記帯域制限画像信号の作成、前記変
換画像信号の作成、前記積算信号の作成、および該積算
信号の前記原画像信号への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)・ｇ（Ｓus1−Ｓus2） ＋ｆ2(Ｓus1 −Ｓus2)・ｇ（Ｓus2−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk )・ｇ（Ｓusk−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（ＳusN−ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換して前記
抑制信号を作成するための関数 ｇ：前記各補助画像信号を変換して前記倍率信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項１または１４
記載の画像処理方法。 - 【請求項１６】 前記帯域制限画像信号の作成、前記変
換画像信号の作成、前記積算信号の作成、および該積算
信号の前記原画像信号への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)・ｇ（Ｓorg−Ｓus2） ＋ｆ2(Ｓus1 −Ｓus2)・ｇ（Ｓorg−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk )・ｇ（Ｓorg−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（Ｓorg−ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換して前記
抑制信号を作成するための関数 ｇ：前記各補助画像信号を変換して前記倍率信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項１または１４
記載の画像処理方法。 - 【請求項１７】 前記複数の変換画像信号の作成を、変
換される帯域制限画像信号である被変換帯域制限画像信
号よりも低い周波数帯域の帯域制限画像信号である低周
波側帯域制限画像信号を、原点をとおり該原点における
傾きがほぼ０で、処理される値が大きくなるにつれて該
傾きが漸増する非線形関数に基づいて変換することによ
り前記被変換帯域制限画像信号の補助画像信号を作成
し、 該補助画像信号を前記被変換帯域制限画像信号に加算す
ることにより複合帯域制限画像信号を作成し、 該複合帯域制限画像信号を、該複合帯域制限画像信号の
絶対値に基づいて決まる該絶対値以下の値となるように
変換することにより行うことを特徴とする請求項１記載
の画像処理方法。 - 【請求項１８】 前記帯域制限画像信号の作成、前記変
換画像信号の作成、前記積算信号の作成、および該積算
信号の前記原画像信号への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝［ｆ1｛（Ｓorg −Ｓus1)＋ｇ（Ｓus1−Ｓus2）｝ ＋ｆ2｛（Ｓus1 −Ｓus2)＋ｇ（Ｓus2−Ｓus3）｝＋… ＋ｆk｛（Ｓusk-1−Ｓusk )＋ｇ（Ｓusk−Ｓusk+1）｝＋… ＋ｆN｛（ＳusN-1−ＳusN)＋ｇ（ＳusN−ＳusN+1）｝］ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各複合帯域制限画像信号を変換する
ための関数 ｇ：前記各帯域制限画像信号を変換して前記補助画像信
号を作成するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項１または１７
記載の画像処理方法。 - 【請求項１９】 前記複数の変換画像信号の作成を、前
記原画像を撮影した際の撮影部位に応じて、前記各帯域
制限画像信号を変換することにより行うことを特徴とす
る請求項１から１８のいずれか１項記載の画像処理方
法。 - 【請求項２０】 前記原画像が断層撮影により得られた
流れ像を含むものであり、前記複数の非鮮鋭マスク画像
信号の作成を、前記原画像信号に対して該流れ像の方向
に沿う１次元非鮮鋭マスクによりマスク処理を施すこと
により行うことを特徴とする請求項１から１９記載の画
像処理方法。 - 【請求項２１】 原画像を表す原画像信号に、該原画像
の高周波成分に関する信号を加算することにより、該原
画像の高周波成分を強調する画像処理装置において、 前記原画像信号に基づいて、互いに周波数応答特性が異
なる複数の非鮮鋭マスク画像信号を作成する非鮮鋭マス
ク画像信号作成手段と、 前記原画像信号および前記複数の非鮮鋭マスク画像信
号、または前記複数の非鮮鋭マスク画像信号に基づい
て、前記原画像信号の複数の周波数帯域ごとの信号を表
す複数の帯域制限画像信号を作成する帯域制限画像信号
作成手段と、 該各帯域制限画像信号のうち少なくとも１つの信号につ
いて、該帯域制限画像信号の少なくとも一部を小さくす
るような変換処理を施して複数の変換画像信号を作成す
る変換手段と、 該各変換画像信号を積算することにより前記原画像信号
に加算する前記高周波成分に関する信号を得る周波数強
調処理手段とからなることを特徴とする画像処理装置。 - 【請求項２２】 前記変換手段が、前記複数の変換画像
信号を、前記帯域制限画像信号の絶対値が所定の閾値よ
り大きい帯域制限画像信号について、該帯域制限画像信
号の絶対値を小さくするように、該各帯域制限画像信号
を変換することにより作成する手段であることを特徴と
する請求項２１記載の画像処理装置。 - 【請求項２３】 前記変換手段が、前記複数の変換画像
信号を、前記各帯域制限画像信号が所定の閾値より大き
いほど該帯域制限画像信号の絶対値を小さくするととも
に、前記各帯域制限画像信号の絶対値が該所定の閾値よ
りも小さい他の閾値よりも小さいほど該帯域制限画像信
号の絶対値を小さくするように、該各帯域制限画像信号
を変換することにより作成する手段であることを特徴と
する請求項２２記載の画像処理装置。 - 【請求項２４】 前記変換手段が、前記複数の変換画像
信号を、前記各帯域制限画像信号の周波数帯域に応じて
該各帯域制限画像信号の絶対値の大きさを変化させるよ
うに、該各帯域制限画像信号を変換することにより作成
する手段であることを特徴とする請求項２２または２３
記載の画像処理装置。 - 【請求項２５】 前記変換手段が、前記複数の変換画像
信号を、前記各帯域制限画像信号を、前記周波数帯域ご
とに異なる複数の関数に基づいて、該帯域制限画像信号
の絶対値に基づいて決まる該絶対値以下の値となるよう
に変換することにより作成する手段であることを特徴と
する請求項２１記載の画像処理装置。 - 【請求項２６】 前記変換画像信号の作成に使用される
各関数が、前記帯域制限画像信号の絶対値が所定の値よ
り大きい帯域制限画像信号については、該変換画像信号
の値がほぼ一定の値となるように該帯域制限画像信号を
変換する関数であり、該関数が高周波帯域を処理する関
数であるほど、前記所定の値が小さい値であることを特
徴とする請求項２５記載の画像処理装置。 - 【請求項２７】 前記変換画像信号の作成に使用される
各関数が、該関数が低周波帯域を処理する関数であるほ
ど、前記帯域制限画像信号の絶対値が０近傍の所定の範
囲内の値である帯域制限画像信号を変換した際に得られ
る該変換画像信号の絶対値が小さい値であることを特徴
とする請求項２５または２６記載の画像処理装置。 - 【請求項２８】 前記帯域制限画像信号の作成、前記変
換画像信号の作成、前記高周波成分に関する信号の作
成、および該高周波成分に関する信号の前記原画像信号
への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)＋ｆ2(Ｓus1 −Ｓus2)＋… ＋ｆk(Ｓusk-1−Ｓusk )＋…＋ｆN(ＳusN-1−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項２１から２７
のいずれか１項記載の画像処理装置。 - 【請求項２９】 前記帯域制限画像信号の作成、前記変
換画像信号の作成、前記高周波成分に関する信号の作
成、および該高周波成分に関する信号の前記原画像信号
への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN ） ＝（1/N)・｛ｆ1（Ｓorg −Ｓus1)＋ｆ2（Ｓorg −Ｓus2)＋… ＋fk(Ｓorg −Ｓusk)＋…＋ｆN（Ｓorg −ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項２１から２７
のいずれか１項記載の画像処理装置。 - 【請求項３０】 前記帯域制限画像信号の作成、前記変
換画像信号の作成、前記高周波成分に関する信号の作
成、および該高周波成分に関する信号の前記原画像信号
への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓus1,Ｓus2,…ＳusN ） ＝｛ｆ2（Ｓus1 −Ｓus2)＋ｆ3（Ｓus2 −Ｓus3)＋… ＋fk(Ｓusk-1 −Ｓusk)＋…＋ｆN（ＳusN-1 −ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=2〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項２１から２４
のいずれか１項記載の画像処理装置。 - 【請求項３１】 前記帯域制限画像信号の作成、前記変
換画像信号の作成、前記高周波成分に関する信号の作
成、および該高周波成分に関する信号の前記原画像信号
への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓus1,Ｓus2,…ＳusN ） ＝（1/N)・｛ｆ2（Ｓus1 −Ｓus2)＋ｆ3（Ｓus1 −Ｓus3)＋… ＋fk(Ｓus1 −Ｓusk)＋…＋ｆN（Ｓus1 −ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=2〜N)：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項２１から２４
のいずれか１項記載の画像処理装置。 - 【請求項３２】 前記変換手段が、前記複数の変換画像
信号を、前記各帯域制限画像信号に対して、前記強調係
数に応じた変換処理を施すことにより作成する手段であ
ることを特徴とする請求項２１から３１のいずれか１項
記載の画像処理装置。 - 【請求項３３】 前記変換処理が、前記各周波数帯域ご
とに異なる変換処理であって、前記帯域制限画像信号の
周波数帯域が低周波帯域であるほど、該帯域制限画像信
号を強く抑制するような変換処理であり、 前記強調係数が大きいときほど、高周波帯域の帯域制限
画像信号に対する抑制の度合いと、低周波帯域の帯域制
限画像信号に対する抑制の度合いの差が大きいことを特
徴とする請求項３２記載の画像処理装置。 - 【請求項３４】 前記変換手段が、前記帯域制限画像信
号を、該帯域制限画像信号の絶対値に基づいて決まる該
絶対値以下の値となるように変換することにより抑制画
像信号を作成する抑制画像信号作成手段と、 前記原画像信号および前記複数の非鮮鋭マスク画像信号
に基づいて、前記抑制画像信号の作成に使用された帯域
制限画像信号よりも低い周波数帯域の信号を含む補助画
像信号を作成する補助画像信号作成手段と、 該補助画像信号を、該補助画像信号の絶対値が小さいほ
ど１に近く、大きいほど０に近い値となるように変換す
ることにより前記抑制画像信号にそれぞれ対応する倍率
信号を作成する倍率信号作成手段と、 前記抑制画像信号に、該抑制画像信号に対応する前記倍
率信号を乗ずる乗算手段とからなることを特徴とする請
求項２１記載の画像処理装置。 - 【請求項３５】 前記帯域制限画像信号の作成、前記変
換画像信号の作成、前記積算信号の作成、および該積算
信号の前記原画像信号への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)・ｇ（Ｓus1−Ｓus2） ＋ｆ2(Ｓus1 −Ｓus2)・ｇ（Ｓus2−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk )・ｇ（Ｓusk−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（ＳusN−ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換して前記
抑制信号を作成するための関数 ｇ：前記各補助画像信号を変換して前記倍率信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項２１または３
４記載の画像処理装置。 - 【請求項３６】 前記帯域制限画像信号の作成、前記変
換画像信号の作成、前記積算信号の作成、および該積算
信号の前記原画像信号への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg −Ｓus1)・ｇ（Ｓorg−Ｓus2） ＋ｆ2(Ｓus1 −Ｓus2)・ｇ（Ｓorg−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk )・ｇ（Ｓorg−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（Ｓorg−ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換して前記
抑制信号を作成するための関数 ｇ：前記各補助画像信号を変換して前記倍率信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項２１または３
４記載の画像処理装置。 - 【請求項３７】 前記変換手段が、変換される帯域制限
画像信号である被変換帯域制限画像信号よりも低い周波
数帯域の帯域制限画像信号である低周波側帯域制限画像
信号を、原点をとおり該原点における傾きがほぼ０で、
処理される値が大きくなるにつれて該傾きが漸増する非
線形関数に基づいて変換することにより前記被変換帯域
制限画像信号の補助画像信号を作成する補助画像信号作
成手段と、 該補助画像信号を前記被変換帯域制限画像信号に加算す
ることにより複合帯域制限画像信号を作成する復号帯域
制限画像信号作成手段と、 該複合帯域制限画像信号を、該複合帯域制限画像信号の
絶対値に基づいて決まる該絶対値以下の値となるように
変換することにより前記変換画像信号を作成する変換画
像信号作成手段とからなることを特徴とする請求項２１
記載の画像処理装置。 - 【請求項３８】 前記帯域制限画像信号の作成、前記変
換画像信号の作成、前記積算信号の作成、および該積算
信号の前記原画像信号への加算を、下記の式 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝［ｆ1｛（Ｓorg −Ｓus1)＋ｇ（Ｓus1−Ｓus2）｝ ＋ｆ2｛（Ｓus1 −Ｓus2)＋ｇ（Ｓus2−Ｓus3）｝＋… ＋ｆk｛（Ｓusk-1−Ｓusk )＋ｇ（Ｓusk−Ｓusk+1）｝＋… ＋ｆN｛（ＳusN-1−ＳusN)＋ｇ（ＳusN−ＳusN+1）｝］ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各複合帯域制限画像信号を変換する
ための関数 ｇ：前記各帯域制限画像信号を変換して前記補助画像信
号を作成するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行うことを特徴とする請求項２１または３
７記載の画像処理装置。 - 【請求項３９】 前記変換手段が、前記複数の変換画像
信号を、前記原画像を撮影した際の撮影部位に応じて、
前記各帯域制限画像信号を変換することにより作成する
手段であることを特徴とする請求項２１から３８のいず
れか１項記載の画像処理装置。 - 【請求項４０】 前記原画像が断層撮影により得られた
流れ像を含むものであり、前記非鮮鋭マスク画像信号作
成手段が、前記複数の非鮮鋭マスク画像信号を、前記原
画像信号に対して該流れ像の方向に沿う１次元非鮮鋭マ
スクによりマスク処理を施すことにより作成する手段で
あることを特徴とする請求項２１から３９記載の画像処
理装置。
Priority Applications (8)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP18215596A JP3816151B2 (ja) | 1995-09-29 | 1996-07-11 | 画像処理方法および装置 |
| EP01118156A EP1156451B1 (en) | 1995-09-29 | 1996-09-27 | Image processing method and apparatus |
| DE69632651T DE69632651T2 (de) | 1995-09-29 | 1996-09-27 | Verfahren und Vorrichtung zur Bildverarbeitung |
| EP01118152A EP1168243B1 (en) | 1995-09-29 | 1996-09-27 | Image processing method and apparatus |
| DE69625433T DE69625433T2 (de) | 1995-09-29 | 1996-09-27 | Verfahren und Vorrichtung zur Bildverarbeitung |
| DE69632700T DE69632700T2 (de) | 1995-09-29 | 1996-09-27 | Verfahren und Vorrichtung zur Bildverarbeitung |
| EP96115572A EP0766202B1 (en) | 1995-09-29 | 1996-09-27 | Image processing method and apparatus |
| US08/723,313 US5991457A (en) | 1995-09-29 | 1996-09-30 | Multi-resolution transforming, image processing and dynamic range compressing method and apparatus |
Applications Claiming Priority (13)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7-252088 | 1995-09-29 | ||
| JP25208895 | 1995-09-29 | ||
| JP32312195 | 1995-12-12 | ||
| JP7-323121 | 1995-12-12 | ||
| JP16906496 | 1996-06-28 | ||
| JP8-169065 | 1996-06-28 | ||
| JP16906596 | 1996-06-28 | ||
| JP16906296 | 1996-06-28 | ||
| JP16906396 | 1996-06-28 | ||
| JP8-169062 | 1996-06-28 | ||
| JP8-169063 | 1996-06-28 | ||
| JP8-169064 | 1996-06-28 | ||
| JP18215596A JP3816151B2 (ja) | 1995-09-29 | 1996-07-11 | 画像処理方法および装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH1075395A true JPH1075395A (ja) | 1998-03-17 |
| JP3816151B2 JP3816151B2 (ja) | 2006-08-30 |
Family
ID=27566282
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP18215596A Expired - Fee Related JP3816151B2 (ja) | 1995-09-29 | 1996-07-11 | 画像処理方法および装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3816151B2 (ja) |
Cited By (36)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000268171A (ja) * | 1999-01-14 | 2000-09-29 | Fuji Photo Film Co Ltd | 画像処理方法および装置並びに記録媒体 |
| JP2001292308A (ja) * | 2000-02-02 | 2001-10-19 | Fuji Photo Film Co Ltd | 画像保存転送方法および装置、並びに画像処理方法、画像処理システム、および画像処理装置 |
| JP2002170114A (ja) * | 2000-12-04 | 2002-06-14 | Fuji Photo Film Co Ltd | 画像処理方法および装置並びに記録媒体 |
| JP2002525724A (ja) * | 1998-09-16 | 2002-08-13 | インテル・コーポレーション | 汎用画像強調フレームワーク |
| US6856706B2 (en) | 2000-02-08 | 2005-02-15 | Fuji Photo Film Co., Ltd. | Image processing method and system, and storage medium |
| WO2005027043A1 (ja) | 2003-09-11 | 2005-03-24 | Matsushita Electric Industrial Co., Ltd. | 視覚処理装置、視覚処理方法、視覚処理プログラム、集積回路、表示装置、撮影装置および携帯情報端末 |
| EP1217826A3 (en) * | 2000-12-19 | 2005-05-25 | Konica Corporation | Image processing apparatus |
| JP2005190468A (ja) * | 2003-12-24 | 2005-07-14 | Sharp Corp | 量子化され、デコードされた画像の画質をエンハンスする方法 |
| JP2005309570A (ja) * | 2004-04-19 | 2005-11-04 | Megachips Lsi Solutions Inc | 画像信号強調装置 |
| JP2006053853A (ja) * | 2004-08-16 | 2006-02-23 | Fuji Photo Film Co Ltd | 認証システム |
| JP2006053852A (ja) * | 2004-08-16 | 2006-02-23 | Fuji Photo Film Co Ltd | 認証システム |
| JP2006127420A (ja) * | 2004-11-01 | 2006-05-18 | Sony Corp | 画像処理装置、画像処理方法及びプログラム |
| US7301670B2 (en) | 2002-02-26 | 2007-11-27 | Konica Corporation | Medical image recording system |
| JP2008033692A (ja) * | 2006-07-28 | 2008-02-14 | Canon Inc | 画像処理装置及びその制御方法、並びに、コンピュータプログラム及びコンピュータ可読記憶媒体 |
| JP2008211840A (ja) * | 2008-05-09 | 2008-09-11 | Konica Minolta Holdings Inc | 画像処理装置 |
| US7456867B2 (en) | 2000-02-02 | 2008-11-25 | Fujifilm Corporation | Image storing-transferring method and apparatus |
| JP2009516882A (ja) * | 2005-11-23 | 2009-04-23 | セダラ ソフトウェア コーポレイション | ディジタル画像を強調する方法及びシステム |
| JP2009199104A (ja) * | 2008-02-19 | 2009-09-03 | Keyence Corp | 画像生成装置、画像生成方法及びコンピュータプログラム |
| JP2010028524A (ja) * | 2008-07-22 | 2010-02-04 | Sony Corp | 画像処理装置、画像処理方法、及びプログラム |
| JP2010136779A (ja) * | 2008-12-10 | 2010-06-24 | Univ Of Miyazaki | 医用画像の処理装置及び方法 |
| US7783126B2 (en) | 2003-09-11 | 2010-08-24 | Panasonic Corporation | Visual processing device, visual processing method, visual processing program, and semiconductor device |
| JP2010244360A (ja) * | 2009-04-07 | 2010-10-28 | Canon Inc | 画像処理装置、画像処理方法、及びコンピュータプログラム |
| JP2011040098A (ja) * | 2004-10-07 | 2011-02-24 | Mitcham Global Investments Ltd | 密度依存シャープニング |
| JP4767345B2 (ja) * | 2008-03-03 | 2011-09-07 | 三菱電機株式会社 | 画像処理装置及び方法、並びに画像表示装置及び方法 |
| JP2011239833A (ja) * | 2010-05-14 | 2011-12-01 | Shimadzu Corp | 画像処理方法およびそれを用いた放射線撮影装置 |
| CN102274037A (zh) * | 2010-06-10 | 2011-12-14 | 株式会社岛津制作所 | 图像处理方法和使用了该图像处理方法的放射线摄影装置 |
| JP2013030191A (ja) * | 2006-03-31 | 2013-02-07 | Nikon Corp | 画像処理方法、画像処理プログラム、および画像処理装置 |
| JP2013031641A (ja) * | 2011-07-05 | 2013-02-14 | Fujifilm Corp | 放射線撮影装置 |
| JP2013215593A (ja) * | 2013-06-10 | 2013-10-24 | Canon Inc | 画像処理装置、画像処理システム、画像処理方法、及び記憶媒体 |
| WO2014133010A1 (ja) * | 2013-02-28 | 2014-09-04 | 日本電気株式会社 | 画像処理方法、及び画像処理装置 |
| JP2015100539A (ja) * | 2013-11-26 | 2015-06-04 | 日立アロカメディカル株式会社 | 超音波診断装置 |
| JP2015104463A (ja) * | 2013-11-29 | 2015-06-08 | 日立アロカメディカル株式会社 | 超音波診断装置 |
| US9286702B2 (en) | 2011-06-15 | 2016-03-15 | Fujifilm Corporation | Radiographic imaging system |
| US9311695B2 (en) | 2010-03-26 | 2016-04-12 | Shimadzu Corporation | Image processing method and radiographic apparatus using the same |
| US10572981B2 (en) | 2017-03-30 | 2020-02-25 | Konica Minolta, Inc. | Radiation image processing apparatus |
| KR20200096398A (ko) * | 2019-02-01 | 2020-08-12 | 삼성전자주식회사 | 에지 보존 필터를 이용하여 이미지를 선명화하는 방법 및 장치 |
-
1996
- 1996-07-11 JP JP18215596A patent/JP3816151B2/ja not_active Expired - Fee Related
Cited By (47)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002525724A (ja) * | 1998-09-16 | 2002-08-13 | インテル・コーポレーション | 汎用画像強調フレームワーク |
| JP2000268171A (ja) * | 1999-01-14 | 2000-09-29 | Fuji Photo Film Co Ltd | 画像処理方法および装置並びに記録媒体 |
| JP2001292308A (ja) * | 2000-02-02 | 2001-10-19 | Fuji Photo Film Co Ltd | 画像保存転送方法および装置、並びに画像処理方法、画像処理システム、および画像処理装置 |
| US7460157B2 (en) | 2000-02-02 | 2008-12-02 | Fujifilm Corporation | Image storing-transferring method and apparatus, image processing method, image processing system, and image processor |
| US7456867B2 (en) | 2000-02-02 | 2008-11-25 | Fujifilm Corporation | Image storing-transferring method and apparatus |
| US8223222B2 (en) | 2000-02-02 | 2012-07-17 | Fujifilm Corporation | Image storing-transferring method and apparatus, image processing method, image processing system, and image processor |
| US6856706B2 (en) | 2000-02-08 | 2005-02-15 | Fuji Photo Film Co., Ltd. | Image processing method and system, and storage medium |
| JP2002170114A (ja) * | 2000-12-04 | 2002-06-14 | Fuji Photo Film Co Ltd | 画像処理方法および装置並びに記録媒体 |
| EP1217826A3 (en) * | 2000-12-19 | 2005-05-25 | Konica Corporation | Image processing apparatus |
| US7301670B2 (en) | 2002-02-26 | 2007-11-27 | Konica Corporation | Medical image recording system |
| US7945115B2 (en) | 2003-09-11 | 2011-05-17 | Panasonic Corporation | Visual processing device, visual processing method, visual processing program, and semiconductor device |
| US7860339B2 (en) | 2003-09-11 | 2010-12-28 | Panasonic Corporation | Visual processing device, visual processing method, visual processing program, intergrated circuit, display device, image-capturing device, and portable information terminal |
| US8165417B2 (en) | 2003-09-11 | 2012-04-24 | Panasonic Corporation | Visual processing device, visual processing method, visual processing program, integrated circuit, display device, image-capturing device, and portable information terminal |
| US7783126B2 (en) | 2003-09-11 | 2010-08-24 | Panasonic Corporation | Visual processing device, visual processing method, visual processing program, and semiconductor device |
| WO2005027043A1 (ja) | 2003-09-11 | 2005-03-24 | Matsushita Electric Industrial Co., Ltd. | 視覚処理装置、視覚処理方法、視覚処理プログラム、集積回路、表示装置、撮影装置および携帯情報端末 |
| JP4544974B2 (ja) * | 2003-12-24 | 2010-09-15 | シャープ株式会社 | 量子化され、デコードされた画像の画質をエンハンスする方法 |
| JP2005190468A (ja) * | 2003-12-24 | 2005-07-14 | Sharp Corp | 量子化され、デコードされた画像の画質をエンハンスする方法 |
| JP2005309570A (ja) * | 2004-04-19 | 2005-11-04 | Megachips Lsi Solutions Inc | 画像信号強調装置 |
| JP2006053852A (ja) * | 2004-08-16 | 2006-02-23 | Fuji Photo Film Co Ltd | 認証システム |
| JP2006053853A (ja) * | 2004-08-16 | 2006-02-23 | Fuji Photo Film Co Ltd | 認証システム |
| JP2011040098A (ja) * | 2004-10-07 | 2011-02-24 | Mitcham Global Investments Ltd | 密度依存シャープニング |
| JP2006127420A (ja) * | 2004-11-01 | 2006-05-18 | Sony Corp | 画像処理装置、画像処理方法及びプログラム |
| JP4656383B2 (ja) * | 2004-11-01 | 2011-03-23 | ソニー株式会社 | 画像処理装置、画像処理方法及びプログラム |
| US8594444B2 (en) | 2005-11-23 | 2013-11-26 | Merge Healthcare Incorporated | Method and system for enhancing digital images |
| JP2009516882A (ja) * | 2005-11-23 | 2009-04-23 | セダラ ソフトウェア コーポレイション | ディジタル画像を強調する方法及びシステム |
| JP2013030191A (ja) * | 2006-03-31 | 2013-02-07 | Nikon Corp | 画像処理方法、画像処理プログラム、および画像処理装置 |
| JP2008033692A (ja) * | 2006-07-28 | 2008-02-14 | Canon Inc | 画像処理装置及びその制御方法、並びに、コンピュータプログラム及びコンピュータ可読記憶媒体 |
| US8195009B2 (en) | 2006-07-28 | 2012-06-05 | Canon Kabushiki Kaisha | Image processing apparatus and control method thereof |
| JP2009199104A (ja) * | 2008-02-19 | 2009-09-03 | Keyence Corp | 画像生成装置、画像生成方法及びコンピュータプログラム |
| JP4767345B2 (ja) * | 2008-03-03 | 2011-09-07 | 三菱電機株式会社 | 画像処理装置及び方法、並びに画像表示装置及び方法 |
| JP2008211840A (ja) * | 2008-05-09 | 2008-09-11 | Konica Minolta Holdings Inc | 画像処理装置 |
| JP2010028524A (ja) * | 2008-07-22 | 2010-02-04 | Sony Corp | 画像処理装置、画像処理方法、及びプログラム |
| JP2010136779A (ja) * | 2008-12-10 | 2010-06-24 | Univ Of Miyazaki | 医用画像の処理装置及び方法 |
| JP2010244360A (ja) * | 2009-04-07 | 2010-10-28 | Canon Inc | 画像処理装置、画像処理方法、及びコンピュータプログラム |
| US9311695B2 (en) | 2010-03-26 | 2016-04-12 | Shimadzu Corporation | Image processing method and radiographic apparatus using the same |
| JP2011239833A (ja) * | 2010-05-14 | 2011-12-01 | Shimadzu Corp | 画像処理方法およびそれを用いた放射線撮影装置 |
| CN102274037A (zh) * | 2010-06-10 | 2011-12-14 | 株式会社岛津制作所 | 图像处理方法和使用了该图像处理方法的放射线摄影装置 |
| US9286702B2 (en) | 2011-06-15 | 2016-03-15 | Fujifilm Corporation | Radiographic imaging system |
| JP2013031641A (ja) * | 2011-07-05 | 2013-02-14 | Fujifilm Corp | 放射線撮影装置 |
| US9542617B2 (en) | 2013-02-28 | 2017-01-10 | Nec Corporation | Image processing device and image processing method for correcting a pixel using a corrected pixel statistical value |
| WO2014133010A1 (ja) * | 2013-02-28 | 2014-09-04 | 日本電気株式会社 | 画像処理方法、及び画像処理装置 |
| JPWO2014133010A1 (ja) * | 2013-02-28 | 2017-02-02 | 日本電気株式会社 | 画像処理方法、及び画像処理装置 |
| JP2013215593A (ja) * | 2013-06-10 | 2013-10-24 | Canon Inc | 画像処理装置、画像処理システム、画像処理方法、及び記憶媒体 |
| JP2015100539A (ja) * | 2013-11-26 | 2015-06-04 | 日立アロカメディカル株式会社 | 超音波診断装置 |
| JP2015104463A (ja) * | 2013-11-29 | 2015-06-08 | 日立アロカメディカル株式会社 | 超音波診断装置 |
| US10572981B2 (en) | 2017-03-30 | 2020-02-25 | Konica Minolta, Inc. | Radiation image processing apparatus |
| KR20200096398A (ko) * | 2019-02-01 | 2020-08-12 | 삼성전자주식회사 | 에지 보존 필터를 이용하여 이미지를 선명화하는 방법 및 장치 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3816151B2 (ja) | 2006-08-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3816151B2 (ja) | 画像処理方法および装置 | |
| JP3700804B2 (ja) | 画像処理方法および装置 | |
| JP3738788B2 (ja) | 画像のダイナミックレンジ圧縮処理方法および装置 | |
| JP3833177B2 (ja) | 画像処理装置、画象処理方法、記憶媒体及びプログラム | |
| EP1550981B1 (en) | System and method for image noise reduction using a minimal error spatio-temporal recursive filter | |
| US6813335B2 (en) | Image processing apparatus, image processing system, image processing method, program, and storage medium | |
| JP3690882B2 (ja) | 画像の強調処理方法および装置 | |
| US20070086671A1 (en) | Image processing apparatus | |
| JP2001057677A (ja) | 画像処理方法および装置並びに記録媒体 | |
| US7149358B2 (en) | Method and system for improving contrast using multi-resolution contrast based dynamic range management | |
| WO2015080006A1 (ja) | 超音波診断装置 | |
| JP2004242285A (ja) | ノイズ抑制処理方法および装置並びにプログラム | |
| JP3700798B2 (ja) | 画像処理方法および装置 | |
| JP3738791B2 (ja) | 画像処理方法および装置 | |
| JP3783116B2 (ja) | 放射線画像強調処理方法および装置 | |
| JP2002074356A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP3696339B2 (ja) | 画像処理方法および装置 | |
| JP4456096B2 (ja) | ノイズ除去方法および装置並びにプログラム | |
| JP3669452B2 (ja) | 画像処理方法および装置 | |
| JP2002183727A (ja) | 画像処理装置 | |
| JP3814421B2 (ja) | 画像処理装置、方法及びコンピュータ読み取り可能な記憶媒体 | |
| JP4014671B2 (ja) | 多重解像度変換方法および装置 | |
| JPH09107479A (ja) | 画像処理方法および装置 | |
| JP4274426B2 (ja) | 画像処理方法および装置並びにプログラム | |
| JP4213337B2 (ja) | 画像処理方法および装置並びに記録媒体 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060307 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060508 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20060606 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060607 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090616 Year of fee payment: 3 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090616 Year of fee payment: 3 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090616 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100616 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100616 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110616 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110616 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120616 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120616 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130616 Year of fee payment: 7 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |