JP2002314816A - 画像処理装置及び方法、並びに、コンピュータプログラム及び記憶媒体 - Google Patents

画像処理装置及び方法、並びに、コンピュータプログラム及び記憶媒体

Info

Publication number
JP2002314816A
JP2002314816A JP2001117548A JP2001117548A JP2002314816A JP 2002314816 A JP2002314816 A JP 2002314816A JP 2001117548 A JP2001117548 A JP 2001117548A JP 2001117548 A JP2001117548 A JP 2001117548A JP 2002314816 A JP2002314816 A JP 2002314816A
Authority
JP
Japan
Prior art keywords
image
component
frequency
absolute value
frequency component
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
Application number
JP2001117548A
Other languages
English (en)
Other versions
JP3793036B2 (ja
Inventor
Hiroyuki Arahata
弘之 新畠
Tomohiko Matsuura
友彦 松浦
Makoto Sato
眞 佐藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Inc
Original Assignee
Canon Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Canon Inc filed Critical Canon Inc
Priority to JP2001117548A priority Critical patent/JP3793036B2/ja
Publication of JP2002314816A publication Critical patent/JP2002314816A/ja
Application granted granted Critical
Publication of JP3793036B2 publication Critical patent/JP3793036B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)

Abstract

(57)【要約】 【課題】 本発明は、対象画像のうち目的とする周波数
帯又は対象物の画像成分を効果的又は選択的に変換する
画像処理装置及び方法、並びに、コンピュータプログラ
ム及び記憶媒体を提供する。 【解決手段】 本発明の一側面としての画像処理装置
は、対象画像から当該対象画像の低周波成分と少なくと
も1つの高周波成分とを作成する作成部と、前記少なく
とも1つの高周波成分のうち所定の高周波成分における
所定の絶対値範囲の絶対値を有する要素を強調又は抑制
する成分変換部と、前記所定の絶対値範囲を決める少な
くとも1つの閾値を設定する設定部と、前記成分変換部
により変換された後の前記少なくとも1つの高周波成分
と、前記対象画像又は前記低周波成分とを用いて画像を
生成する画像生成部とを有することを特徴とする。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、画像処理装置及び
方法に関し、特に、画像に対し鮮鋭化等の周波数強調又
は抑制処理を施す画像処理装置及び方法に関する。
【0002】
【従来の技術】近年のデジタル技術の進歩により、X線
画像等の放射線画像をデジタル信号に変換し、かかるデ
ジタル画像に画像処理を施して表示装置(例えば、CR
T、液晶ディスプレイ等)に表示、又は記録装置(プリ
ンタ等)によりフィルム等の記録媒体に記録することが
行われている。かかる画像処理は、撮像装置から得られ
た画像を、撮像装置の特性等に依存して補正する前処理
と、前処理を経た画像(原画像)を診断に適した画質の
画像に変換する画質保証(QA)処理とに分類され、こ
のうちQA処理は原画像の高周波成分を強調する鮮鋭化
処理等の周波数処理を含む。
【0003】鮮鋭化処理は原画像の高周波成分を調節し
たうえで原画像又は平滑化画像に加算することで高周波
成分を強調し、原画像の鮮鋭度を上げることができる。
ここで、平滑化画像とは原画像の低周波成分を抽出する
ことで得られる画像のことであり、低周波画像又はボケ
画像とも呼ばれる(本願では交換可能に使用する)。鮮
鋭化処理では強調する周波数成分(例えば周波数帯域)
を平滑化画像を作成するための画像範囲(例えば平滑化
フィルタのマスクサイズ)に基づき調整している。
【0004】また、近年では、これと同様の効果を奏す
るものとして、原画像を表す画像信号を複数の周波数帯
の画像成分に分離し、周波数帯毎の画像成分を増加又は
減弱する多重周波数処理も用いられ始めている。一般
に、原画像中には鮮鋭度を増したい対象(例えば肺内の
血管)と、鮮鋭度を落としたい対象(例えばノイズ)が
存在する。多重周波数処理では目的とする対象の鮮鋭度
をあげるように処理されるのが普通である。例えば、鮮
鋭度を上げたい対象部を構成する主要な周波数帯が中周
波数の帯域にある場合には、中周波数帯域の画像成分を
強調することで目的とする対象部の鮮鋭度を一律に増加
することが行われている。
【0005】
【発明が解決しようとする課題】しかしながら、鮮鋭化
処理では平滑化画像を原画像から減じることにより高周
波成分を作成するため、強調したい(目的とする)周波
数帯の成分のみを選択することが困難である。一般的
に、高周波成分には広範囲(目的とする周波数帯より広
い周波数帯)の周波数成分が含まれている。このため、
鮮鋭化処理ではこの広範囲の周波数成分を含む高周波成
分を一律に強調するため、目的外の周波数成分も強調し
てしまうことがある。また、画像のエッジ部分位置に対
応する画像成分は、低周波数から高周波数までの広範囲
にわたり振幅の大きな周波数成分から構成される。この
ため、広範囲の高周波成分を一律に強調する鮮鋭化処理
では、エッジ部分にオーバーシュート等のアーティファ
クトを生じる問題がある。
【0006】一方、目的とする対象物が特定の周波数成
分(例えば、中周波数成分)をその構成要素の主とする
場合、上述したように多重周波数処理は特定の周波数の
画像成分のみを強調することができる。しかし、特定の
周波数を示す画像成分のなかには非対象物(目的と異な
る対象物)の画像成分も存在する。従って、一律に特定
周波数の画像成分を強調すると非対象物までも強調され
てしまう。また、多重周波数処理で特定周波数の周波数
成分を一律に強調する場合であっても、上述の鮮鋭化処
理と同様、エッジ部分にアーティファクトが生じる問題
がある。
【0007】このため、上述のような処理を施した画像
(処理後画像)では、非対象物が過度に強調され、又は
エッジ部分のプロファイルがオーバーシュートにより崩
れ、不自然な画像を得ることになってしまう。かかる処
理後画像に、特に診断に供される放射線画像等の医療画
像に、このようなアーティファクトが生じることは好ま
しくない。
【0008】そこで、対象画像のうち目的とする周波数
帯又は対象物の画像成分を効果的又は選択的に変換する
画像処理装置及び方法、並びに、コンピュータプログラ
ム及び記憶媒体を提供することを本発明の例示的目的と
する。
【0009】
【課題を解決するための手段】上記目的を達成するため
に、本発明の一側面としての画像処理装置は、対象画像
から当該対象画像の低周波成分と少なくとも1つの高周
波成分とを作成する作成部と、前記少なくとも1つの高
周波成分のうち所定の高周波成分における所定の絶対値
範囲の絶対値を有する要素を強調又は抑制する成分変換
部と、前記所定の絶対値範囲を決める少なくとも1つの
閾値を設定する設定部と、前記成分変換部により変換さ
れた後の前記少なくとも1つの高周波成分と、前記対象
画像又は前記低周波成分とを用いて画像を生成する画像
生成部とを有することを特徴とする。かかる画像処理装
置は、前記所定の絶対値範囲を決める前記閾値を設定可
能である。よって、かかる前記閾値を変更することで所
望の高周波成分のみを効果的に強調することができ、画
像の鮮鋭化に寄与する。その結果、本発明の画像処理装
置はより自然な感じの処理後画像を得ることができる。
また、かかる画像処理装置は前記成分変換部が前記所定
の高周波成分を変換する際に用いる変換関数を表示する
表示部を有しても良い。かかる表示部は変換関数を表示
することで、オペレータが視覚的に確認しながら処理後
画像との関係を把握することができる。また、かかる画
像処理装置は前記成分変換部が前記所定の高周波成分を
変換する際に用いる変換関数を変更するための入力部を
有してもよい。かかる入力部は前記変換関数を変更する
ことができるので、適切な処理後画像が得られるように
前記変換関数を変更することができる。
【0010】かかる画像処理装置は、強調する対象の画
像成分が強調する対象でない画像成分に比べて相対的に
高い周波数側に分布する場合、前記成分変換部は前記所
定の高周波成分の所定絶対値以上の絶対値を有する要素
を相対的に減ずるように前記所定の高周波成分を変換す
ることを特徴とする。また、かかる画像処理装置は、前
記成分変換部は前記所定の高周波成分の所定絶対値以上
の絶対値を有する要素を相対的に減ずるように前記所定
の高周波成分を変換することを特徴とする。かかる画像
処理装置は、所定の閾値を境に高周波成分を減じする変
換を行うことで、強調する対象でない画像成分が強調さ
れることを抑制し、画像の鮮鋭化に寄与する。前記少な
くとも1つの閾値は予め決められた固定値であってもよ
い。また、前記成分変換部は前記少なくとも1つの閾値
を少なくとも前記対象画像に基づいて決定してもよい。
【0011】更に、前記画像生成部は前記成分変換部に
より変換された後の前記少なくとも1つの高周波成分と
前記対象画像とを用いて画像を生成すると共に、前記成
分変換部は前記所定の高周波成分の前記絶対値範囲外の
絶対値を有する要素が実質的に0となるように前記所定
の高周波成分を変換することを特徴とする。また、前記
画像生成部は前記成分変換部により変換された後の前記
少なくとも1つの高周波成分と前記低周波成分とを用い
て画像を生成すると共に、前記成分変換部は前記所定の
高周波成分の前記絶対値範囲外の絶対値を有する要素が
実質的に不変となるように前記所定の高周波成分を変換
してもよい。
【0012】更に、前記画像処理装置は前記成分変換部
が前記所定の高周波成分を変換する際に用いる変換関数
はその傾きが連続であることを特徴とする。前記成分変
換部は前記関数の微分値が一定となるように作成するこ
とで、処理後画像の偽輪郭を抑制することができる。ま
た、上述の画像処理装置は前記設定部は前記閾値を、前
記所定の高周波成分がより低い周波数の画像成分を含む
ほど大きくすることを特徴とする。かかる画像処理装置
は、最適な閾値を決定することができ、周波数の漏れを
防ぐ効果がある。
【0013】また、本発明の画像処理装置の前記成分変
換部は、前記作成部により作成される複数の高周波成分
の各々の所定の絶対値範囲の絶対値を有する要素を強調
又は抑制することを特徴とする。例えば、強調する対象
の画像成分が強調する対象でない画像成分に比して相対
的に高い周波数側に分布する場合、前記成分変換部は前
記複数の高周波成分の各々の所定絶対値以上の絶対値を
有する要素を相対的に減ずるように前記複数の高周波成
分を変換する。また、前記成分変換部は前記複数の高周
波成分の各々の所定絶対値以上の絶対値を有する要素を
相対的に減ずるように前記複数の高周波成分を変換す
る。前記複数の高周波成分の各々の前記所定の絶対値範
囲を決める少なくとも1つの閾値は予め決められた固定
値であってもよい。また、前記成分変換部は前記複数の
高周波成分の各々の前記所定の絶対値範囲を決める少な
くとも1つの閾値を少なくとも前記対象画像に基づいて
決定してもよい。また、前記成分変換部は前記複数の高
周波成分の各々の前記所定の絶対値範囲を決める少なく
とも1つの閾値を当該高周波成分の画像成分の平均値よ
り算出してもよい。これにより、有効に他の周波数帯の
成分を除去でき、周波数の漏れを防ぐ効果がある。
【0014】本発明の別の一側面としての画像処理方法
は、対象画像から当該対象画像の低周波成分と少なくと
も1つの高周波成分とを作成するステップと、前記少な
くとも1つの高周波成分のうち所定の高周波成分におけ
る所定の絶対値範囲の絶対値を有する要素を強調又は抑
制するステップと、前記所定の絶対値範囲を決める少な
くとも1つの閾値を設定するステップと、前記強調又は
抑制ステップにより変換された後の前記少なくとも1つ
の高周波成分と、前記対象画像又は前記低周波成分とを
用いて画像を生成するステップとを有することを特徴と
する。かかる画像処理方法は、上述の画像処理装置の機
能であって、同様の作用を奏する。また、かかる処理方
法の処理ステップをコンピュータに実行させるためのプ
ログラム又は上述の画像処理装置の機能をコンピュータ
に実現させるためのプログラム、及び当該プログラムを
格納する記録媒体も本発明の一側面として機能する。
【0015】
【発明の実施の形態】以下、添付図面を参照し、本発明
の一側面としての画像処理装置100を説明する。な
お、各図において同一の参照符号は同一部材を示し、重
複説明は省略する。また、アルファベットを付した参照
符号はアルファベットのない参照符号の変形例を表し、
特に断らない限りアルファベットのない参照符号はアル
ファベットを付した参照符号を総括するものとする。本
発明の画像処理装置100は、例えば、X線撮影を行う
ときに用いられるX線撮影装置、又は当該X線撮影装置
の一部(例えば、X線撮影装置の画像処理部)として実
現される。
【0016】図1を参照するに、画像処理装置100
は、CPU110と、メモリ120と、操作パネル13
0と、表示部140と、前処理部150と、データ収集
部160と、画像処理部170とを有する。画像処理装
置100はCPUバス105を有し、各要素がCPUバ
ス105を介し互いにデータ授受可能に構成されてい
る。ここで、図1は、本発明の画像処理装置100を示
すブロック図である。なお、本実施例に示すように、画
像処理装置100はデータ収集部160に接続された撮
像部190を有し、当該撮像部190によって撮影され
たX線画像の画像処理を一装置において可能な構成とし
ている。しかし、本実施例の画像処理装置100は他の
撮像機構に任意に接続可能に構成されても良く、また、
単に後述する画像処理(鮮鋭化等の周波数処理)を達成
可能な構成のみであってもよい。
【0017】かかる構成において、画像処理装置100
は撮像部190から得られる画像(生画像)を前処理部
150によって前処理(例えば、オフセット補正、ゲイ
ン補正、Log補正)を施すことで原画像を作成する。
原画像とは生画像に対して前処理が施され、例えば撮像
部190の特性に依存した補正が行われ体裁の整えられ
た画像を意味する。しかしながら、原画像は医師が診断
を行うには不十分であるため、更に画像処理部170に
よって鮮鋭化処理等を施すことにより最も診断に適し
た、又は所望の画像を得ることができる。なお、この最
適な画像又は所望の画像をQA画像(画質保証画像)と
表現する場合もある。
【0018】CPU110はMPUなど名称の如何を問
わずいかなるプロセッサであってもよく、CPUバス1
05に接続されている。CPU110はCPUバス10
5を介し各部の動作を制御する。また、CPU110は
メモリ120に格納されたプログラムを用いて、オペレ
ータの操作パネル110への操作に従った画像処理装置
100全体の動作制御等を行う。
【0019】メモリ120は、例えば、画像処理装置1
00の動作プログラムや処理に必要な各種のデータなど
を格納するROMなどの不揮発性メモリと、画像及び必
要な制御プログラムを一時的に格納するRAMなどの揮
発性メモリとを含む。
【0020】操作パネル130は、例えば、キーボー
ド、スイッチ、マウス、タッチパネル等より構成され、
オペレータの画像処理装置100の操作を可能とする。
また、表示部140は、例えば、CRT、液晶ディスプ
レイを含む表示装置及び/又は印刷装置を含み、QA画
像を出力可能に構成されている。
【0021】前処理部150は生画像に対する前処理を
実行可能に構成された回路を有し、CPUバス105及
びデータ収集部160に接続されている。なお、本実施
例における前処理は、撮像部190の後述する2次元X
線センサ194の各画素(図示しない)の特性(暗電
流、検出感度等)の差によって生じる画像の誤差を補正
する処理を含む。より特定的には、前処理とはオフセッ
ト補正、ゲイン補正、対数変換(Log変換)等を指
し、かかる技術には当業界で周知のいかなる技術をも適
用可能である。前処理部150はデータ収集部160よ
り授受された生画像(又は、メモリ120を介しデータ
収集部160より授受された生画像)をメモリ120に
格納されたデータを基に、CPU110の制御に基づき
補正を行う。なお、前処理部150は画像処理を実行す
る部分の一つであり、後述する画像処理部170の一部
として構成されてもよい。前処理部150はかかる処理
を施した画像(原画像)を画像処理部170及び/又は
メモリ120に出力する。
【0022】データ収集部160は撮像部190、前処
理部150及びCPUバス105に接続され、撮像部1
90より出力された生画像を所定の電気信号に変換し前
処理部150及び/又はメモリ120に供給する。デー
タ収集部160は、例えば、14ビットA/D変換機を
有し、撮像部190の出力に比例したデジタル信号を前
処理部150及び/又はメモリ120に供給する。これ
により前処理部150は上述した前処理、例えば、デジ
タル信号を対数変換し、X線の線量の対数に比例したデ
ジタル信号に変換することができる。
【0023】画像処理部170は、低周波成分作成部1
72と、高周波成分作成部174と、変換曲線作成部1
76と、成分変換部178と、高周波成分加算部180
とを有し、CPUバス105を介し各々データを授受可
能に構成されている。低周波成分作成部172は原画像
から平滑化画像を作成する。高周波成分作成部174は
原画像から平滑化画像を減じることで高周波成分(高周
波画像)を作成する。変換曲線作成部176は高周波成
分作成部で得られた高周波成分を成分変換部178で変
換する際の所定の閾値を決定すると共に、変換曲線(変
換関数)を作成する。なお、かかる変換曲線について
は、後述する動作において詳細に説明するものとし、こ
こでの説明を省略する。成分変換部178は高周波成分
作成部174で作成された高周波成分を、変換曲線作成
部176で作成された変換曲線に基づき変換する。高周
波成分加算部180は、成分変換部174によって変換
された高周波成分を原画像又は平滑化画像に加算する。
【0024】本実施例の画像処理部170は成分変換部
178において高周波成分を所定の関数で変換すること
を特徴としている。画像処理部170は所定の絶対値
(閾値)以上(又は以下)の絶対値を有する高周波成
分、即ち目的とする範囲外の絶対値を有する高周波成分
の絶対値を減じて原画像(又は平滑化画像)に加算する
ことができる。これにより、処理後画像のノイズを含む
非対象物の高周波成分を相対的に抑制し、対象物の高周
波成分を相対的に強調することができる。また、画像処
理部170はオーバーシュートを抑制する効果がある。
【0025】また、画像処理部170は後述する実施例
で述べるように、表示回路182及び入力回路184を
有してもよい。本実施例において、表示回路182は処
理後画像(即ち、QA画像)、QA画像のプロファイ
ル、所定の関数形を含む画像処理に関する情報を表示部
140に表示可能とする。また、入力回路184は後述
する閾値及び所定の関数を変更するパラメータ、並びに
画像のプロファイル表示の対象となる水平ライン位置を
指定するパラメータ等を含む情報を操作パネル130を
介し入力することで、画像処理部170の変換曲線作成
部176、成分変換部178に反映する。また、かかる
情報は、表示回路182を介し表示部140に表示可能
に構成される。なお、表示回路182及び入力回路18
4には当業界で周知のいかなる技術をも適用可能であっ
て、ここでの詳細な説明は省略する。また、かかる表示
回路182及び入力回路184の機能は後述する実施例
にてより詳細に理解されるであろう。
【0026】なお、図17に示すように、画像処理部1
70は画像処理部170aに置換されてもよい。ここ
で、図17は、本発明の別の画像処理装置100aを示
すブロック図である。画像処理部170aは離散ウェー
ブレット変換部(DWT変換部)171と、変換曲線作
成部173と、成分変換部175と、逆DWT変換部1
77と、特徴抽出部179と、階調変換部181とを有
し、CPUバス105を介し各々データを授受可能に構
成されている。
【0027】DWT変換部171は原画像に対して離散
ウェーブレット変換(DWT変換と称する場合もある)
を行い、原画像を複数の周波数帯に分解したウェーブレ
ット変換係数(以下、画像成分)を得る。DWT変換部
171は図18に示すような回路より構成可能である
が、これに限定されず当業界周知のいかなる技術の適用
も制限するものではない。ここで、図18は図17に示
すDWT変換部171を構成可能な回路を例示的に示し
た図である。変換曲線作成部173はDWT変換部17
1で得られた各周波数帯の画像成分を成分変換部175
で変換する際の所定の閾値を設定すると共に、変換曲線
(変換関数)を作成する。なお、かかる変換曲線につい
ては、後述する動作において詳細に説明するものとし、
ここでの説明を省略する。成分変換部175はDWT変
換部171で得られた各周波数帯の画像成分を、変換曲
線作成部173で作成された変換曲線に基づき変換す
る。逆DWT変換部177は成分変換部175で変換さ
れた後の周波数係数に基づき逆離散ウェーブレット変換
(逆DWT変換)を行う。逆DWT変換部177は、例
えば、図19に示すような回路で構成可能であるが、当
業界周知のいかなる技術の適用も制限するものではな
い。ここで、図19は逆離散ウェーブレット変換処理を
施す例示的な回路構成を示した図である。かかる構成に
おいて、逆変換部177は変換された後の周波数係数を
逆変換することで画像を合成する。特徴抽出部179は
原画像から階調変換のための特徴量を算出する回路を備
える。また、階調変換回路181は特徴抽出部179で
算出した特徴量に基づき、逆DWT変換部177で合成
した画像を階調変換する回路を備える。
【0028】画像処理部170aは所定の絶対値(閾
値)以上(又は以下)の画像成分、即ち目的とする範囲
外の絶対値を有する画像成分の絶対値を減じることがで
きる。これにより、処理後画像はノイズを含む非対象物
の画像成分を相対的に抑制し、対象物の画像成分を相対
的に強調することができる。
【0029】また、画像処理部170aは後述する実施
例で述べるように、表示回路182及び入力回路184
を有してもよい。表示回路182及び入力回路184
は、上述した画像処理部170と同様の構成であって、
ここでの詳細な説明は省略する。なお、画像処理部17
0aにおいて、表示回路182は変換曲線作成部173
で作成した変換曲線を表示部140に表示可能に構成さ
れる。また、入力回路184は操作パネル130を介し
変換曲線のパラメータを入力することで、かかる曲線形
を変更可能に構成される。
【0030】なお、画像処理部170と画像処理部17
0aは画像の処理方法において異なるが、各処理方法は
後述する動作の説明において容易に理解されるであろ
う。
【0031】撮像部190は、X線発生部192と2次
元X線センサ194とを有し、データ収集部160に接
続される。図1に良く示されるように、撮像部190は
被検査体Pを介在させてX線発生部192と2次元X線
センサ194とが対向する位置に配置されて構成され
る。かかる構成において、撮像部190は被検査体Pに
おいて吸収及び散乱等の相互作用を経たX線を受け、X
線画像(生画像)をデータ収集部160に供給する。X
線発生部192はオペレータによって操作可能であって
(オペレータが操作パネル130を操作することによ
り、メモリ120に格納された動作プログラムに従って
CPU110が制御することの意味も含む)、例えばX
線管球を含む回路より構成される。一方、2次元センサ
194は例えば、X線の入射側から順に蛍光体とアモル
ファス・シリコン光センサとから構成される積層構造を
有し、データ収集部160に接続されている。なお、撮
像部190が上述の構成に限定されず、当業界周知の技
術を適用可能であることは言うまでもない。また、撮像
部190は画像処理装置100から独立した構成要素で
あってもよく、本発明の画像処理装置100が必ずしも
撮像部190を必要とすることを意味するものではな
い。例えば、撮像部190が独立の装置として画像処理
装置100に接続されても良い。
【0032】以下、上述の画像処理装置100の動作を
説明する。なお、以下の説明において、画像処理装置1
00は、例えば、医療用X線撮影装置として具体化され
ている。
【0033】先ず、操作パネル130にてオペレータか
ら撮影の指示がなされると、X線発生部192は被検査
体Pに対してX線ビームを放射する。X線発生部192
から放射されたX線ビームは、被検査体Pを減衰しなが
ら透過して、2次元X線センサ194に到達する。かか
るX線ビームは2次元X線センサ194によりX線画像
として出力される。ここでは、2次元X線センサ194
から出力されるX線画像を、例えば、人体部画像等とす
る。
【0034】データ収集部160は、2次元X線センサ
194から出力されたX線画像をデジタル信号に変換し
て前処理部150に供給する。前処理部150は、デー
タ収集部160から出力されたデジタル信号に対して、
オフセット補正処理やゲイン補正処理等の前処理を行
う。この前処理部150で前処理が行われた信号は原画
像として、CPU110の制御により、CPUバス10
5を介しメモリ120及び/又は画像処理部170に転
送される。
【0035】画像処理部170(又は、画像処理部17
0a)に転送された原画像は後述する実施例で説明され
る画像処理部170(又は、画像処理部170a)にお
いて所定の処理を経て、処理後画像として表示部140
に出力される。オペレータ、例えば、撮影技師又は医師
等は表示部140に出力される画像によって、撮影され
た画像の確認又は診断を行うことができる。
【0036】
【実施例】実施例1 以下、図1及び図2を参照するに、画像処理部170の
動作、即ち、本発明の画像処理装置100の画像処理方
法について説明する。ここで、図2は、本発明の一側面
としての画像処理方法を示すフローチャートである。
【0037】上述したように、画像処理部170は前処
理部150より出力された画像(原画像)を、CPUバ
ス105を介し授受する。より特定的には、前処理部1
50より出力された原画像は低周波成分作成部172に
転送される。まず、低周波成分作成部171は原画像か
ら以下に示す数式1に従い平滑化画像を作成する(ステ
ップ1000)。ここで、原画像をf(x、y)、平滑
化画像をf1(x、y)、d1、d2をマスクサイズと
する。
【0038】
【数1】
【0039】なお、ステップ1000によって作成され
る平滑化画像f1(x、y)の作成方法は数式1に示す
方法に限定されない。例えば、モルフォロジカルフィル
タリング等どのような手法を用いてもよい。
【0040】次に、低周波成分作成部172で得られた
平滑化画像f1(x、y)は高周波成分作成部174に
転送される。高周波成分作成部174は数式2に従い高
周波成分(高周波画像)を作成する。即ち、高周波成分
作成部174は低周波成分作成部172で得られた画像
f1(x、y)と原画像f(x、y)との差分を計算し
高周波画像を抽出する(ステップ1005)。ここで、
f2(x、y)を高周波成分作成部174で得られる高
周波成分とする。
【0041】
【数2】
【0042】高周波成分作成部174で得られた高周波
成分f2(x、y)は成分変換部178に転送される。
【0043】変換曲線(変換関数)G()は、例えば、
図3乃至図5及び図8乃至図11に示す曲線形で表さ
れ、変換曲線作成部176にて作成される(ステップ1
010)。ここで、図3乃至図5及び図8乃至図11
は、高周波成分f2(x、y)を変換する変換曲線
G()の例を示した図であるが、実際にはいずれも奇関
数であり、便宜上、各変換関数の第1象限のみを示して
いる。なお、図3乃至図5は、所定の高周波成分中の所
定の絶対値(閾値SL1)以下の絶対値を有する要素の
絶対値を減じる変換曲線の例を示す図である。図8乃至
図11は、所定の高周波成分中の(少なくとも)所定の
絶対値(閾値SL2)以上の絶対値を有する要素の絶対
値を減じる変換曲線の例を示す図である。
【0044】図3乃至図5に示した変換曲線G()は、
所定の高周波成分の所定の絶対値(閾値SL1)より小
さな絶対値を有する要素の絶対値を相対的に減少させ
(例えば“0”)、閾値SL1以上の絶対値を有する要
素の絶対値を相対的に増加させる(例えば不変とする)
ものである。
【0045】例えば、図6に示す様に、鮮鋭度を上げた
い対象物(例えば、血管等)の主な周波数成分が中周波
数成分であり、鮮鋭度を落としたい非対象物(例えば、
ノイズ等)の主な周波数成分が高周波数成分である場
合、中周波数に対応する高周波成分を増幅することで血
管などの鮮鋭度を選択的に上げることが可能となる。こ
こで、図6は、対象物と非対象物の周波数分布を示す図
であり、横軸が周波数、縦軸が出現頻度を示す図であ
る。なお、本出願において鮮鋭度をあげるとは目的とす
る対象物の周波数成分の振幅を相対的に増加することを
いう。
【0046】しかしながら、非対象物と対象物の周波数
分布には重なり合う領域が存在する。従って、中周波数
帯の高周波成分を一律に増加すると非対象物の強調も同
時になされてしまう問題が生じる。ただし、図7の実線
で示す様に中周波数帯の高周波成分の値としては、対象
物の高周波成分は絶対値の大きい方に偏在し、非対象物
の高周波成分は絶対値の小さい方に偏在する。ここで、
図7は、中周波数帯の高周波成分における対象成分と非
対象成分の出現頻度を示す図であり、横軸が高周波成分
の絶対値、縦軸が成分の出現頻度である。このため、中
周波数成分の成分(要素)を相対的に増加する場合に、
ノイズ成分が多く含有される範囲(例えば、成分の絶対
値が閾値SL1より小さい範囲)では当該成分の絶対値
を相対的に減ずるように変換することで非対象部を高周
波成分を相対的に抑制すると共に、対象物の成分を多く
含む領域(例えば、成分の絶対値が閾値SL1以上の範
囲)の高周波成分を相対的に増加させることで対象物の
鮮鋭度を選択的に上げることができる。
【0047】さらに、図7の点線に示す様に、より低周
波帯の高周波成分では高周波成分中に含まれる非対象物
の成分の絶対値が大きくなる傾向がある。従って、より
低周波になるに従い閾値SL1を大きくすることにより
非対象物の成分をより効果的に抑制することができる。
本実施例1では高周波成分に主に含まれる周波数帯は平
滑化フィルタのマスクサイズで決まるものであり、マス
クサイズが大きくなるほど、より低周波の成分が高周波
成分に含有されることとなる。
【0048】以上より、閾値SL1としては図14の実
線に示す様に、平滑化フィルタのマスクサイズが大きく
なる(高周波成分がより低周波数の帯域を含む)に従い
閾値SL1を大きくすることが、対象物を選択的に強調
するには有効となる。ここで、図14は閾値とマスクサ
イズの関係を示す図であり、横軸がマスクサイズであ
り、縦軸が閾値を示す。
【0049】また、例えば、高周波成分には、中周波数
帯、高周波数帯の成分も含まれる。このように高周波成
分中の特定周波数帯(所望の周波数帯)に対し、高周波
成分中に含まれる他周波数帯の成分を他周波数帯成分の
漏れと呼ぶ。この場合に、目的とする周波数帯(中周波
帯)の高周波成分に含まれる他の周波数帯(高周波数
帯、低周波数帯)の成分の絶対値は比較的に小さな値と
して現れる。そこで、高周波成分中の値の小さい成分
(要素)(例えば、閾値SL1以下の絶対値を有する成
分)を抑制することで、周波数成分の漏れの影響を抑制
できる効果があり、目的とする周波数帯の成分を選択的
に強調できる効果がある。また、特定周波数帯に含まれ
る他周波数帯の成分の絶対値は低周波帯に行くに従い大
きくなる傾向がある。よって、より低周波帯になるに従
い閾値SL1を大きくすることは、より効果的に他の周
波数帯の漏れ成分の影響を除去する効果を有する。
【0050】再び図4を参照するに、変換曲線は高周波
成分の絶対値が所定閾値(SL1)以下で0に変換され
るようにしたものである。つまり、原画像には閾値SL
1以下の高周波成分は加算されない。これにより上述の
効果を達成するものである。また、図4に示す変換曲線
は図3のそれと同様の効果を奏するものであるが、閾値
SL1以下の成分は0とすると共に、変換曲線の立ち上
がりがSL1で出力“0”から始まるようになっている
ものである。更に、図5に示す変換曲線も図3のそれと
同様な効果を有するが、閾値SL1以下の高周波成分も
出力を“0”とせず、抑制しつつ加算することで、加算
される周波数成分の入力値範囲に途切れがない。これに
より処理後画像はより自然な画像となる。
【0051】次に、図8乃至11に示す変換曲線G()
は、高周波成分が所定の閾値(SL2)より大きな絶対
値有する成分を相対的に減少するとともに、それ以外の
高周波成分を相対的に増加するものである。
【0052】例えば、図12に示す様に鮮鋭度を上げた
い対象物(血管等)の主な周波数成分が中周波数成分で
あり、鮮鋭度を落としたい非対象物(臓器等)の主な周
波数成分が低周波数成分である場合、中周波数に対応す
る高周波成分を増幅することで血管などの鮮鋭度を上げ
ることが可能となる。ここで、図12は、対象物と非対
象物の周波数分布を示す図であり、横軸が周波数、縦軸
が出現頻度である。
【0053】しかしながら、対象物と非対象物の周波数
分布には重なり合う領域が存在し、中周波数帯の高周波
成分を一律に増幅すると非対象物の高周波成分の強調が
起きてしまう問題が生じる。但し、図13の実線に示す
様に、中周波数帯の高周波成分の絶対値として対象物の
高周波成分は絶対値の小さい方に偏在し、非対象物の高
周波成分は絶対値の大きい方に偏在する。ここで、図1
3は、中周波数帯の高周波成分における対象成分と非対
象成分の絶対値分布を示す図であり、横軸が高周波成分
の絶対値、縦軸が出現頻度である。このため、中周波数
成分を相対的に強調する場合に、非対象物の成分が多く
含有される絶対値範囲(例えば、閾値SL2を超える範
囲)を相対的に減少させることで非対象物の成分を相対
的に抑制するとともに、対象物の成分を多く含む絶対値
範囲(例えば、閾値SL2以下の範囲)の高周波成分を
相対的に増加させることで対象物の鮮鋭度を選択的に相
対的に上げることができる。
【0054】更に、図13の破線で示す様に、より低周
波帯の高周波成分では高周波成分中に含まれる対象物の
成分の絶対値が大きくなる傾向があるため、より低周波
になるに従い閾値SL2を大きくすることにより、対象
物の成分をより効果的に選択強調することができる。本
実施例1では高周波成分に主に含まれる周波数帯は平滑
化フィルタのマスクサイズで決まるものであり、マスク
サイズが大きくなるほどより低周波の成分が高周波成分
の中に含有されることとなる。
【0055】以上より、閾値SL2としては図14の破
線で示す様に、平滑化フィルタのマスクサイズが大きく
なる(高周波成分がより低周波数の帯域を含む)に従い
閾値SL2を大きくすることが、対象物を選択的に強調
するには有効となる。
【0056】また、画像中のエッジ部分(例えば、人体
では骨の縁、金属片(ペースメーカ等)の縁)の位置に
対応する周波数成分は低周波から高周波まで広い範囲の
周波数成分を含有する。また、その周波数成分の振幅が
大きいため、対応する高周波成分も比較的大きな絶対値
を有する。そこで、所定値以上の高周波成分を抑制する
ことで画像中のエッジ部分にオーバーシュートなどのア
ーティファクトが生じるのを抑制することもできる。例
えば、所定閾値SL2以上の高周波成分を0とする場合
には、エッジ成分に対応する高周波成分が原画像に加算
されずエッジ構造がそのまま保存される効果がある。
【0057】図8に示す変換曲線は所定閾値SL2未満
の高周波成分を相対的に強調し、所定閾値SL2以上の
高周波成分は0としている。これにより対象物の高周波
成分は強調され、非対象物の成分は強調されず、上述の
効果が達成されるものである。また、エッジ部分に相当
する高周波成分を変更しないのでオーバーシュートなど
のアーティファクトが生じるのを抑制することもでき
る。また、図9に示す変換曲線は図3及び図8に示す関
数G()と同一の効果を同時に達成することを目的と
し、所定絶対値範囲の成分(閾値SL1とSL2の間の
絶対値を有する成分)以外の成分を0とするものであ
る。また、図10に示す変換曲線は一定閾値SL2以下
の成分を相対的に強調し、それ以外の成分を相対的に抑
制するようにしたものである。閾値SL2以上の絶対値
を有する成分を0としないことで加算される周波数成分
の絶対値範囲が連続に存在するので処理後の画像がより
自然となる。なお、かかる関数形の変換曲線であって
も、定性的には上述した効果と同様の効果を奏するもの
である。更に、図11に示す変換曲線は図8に示す変換
曲線と同様の効果をねらったものである。しかし、図1
1に示す変換曲線は、関数G()を、その両端を除き微
分連続(傾きが連続)としたことで、偽輪郭の発生等も
抑制しより自然な処理後画像が得られるものである。
【0058】また、図14の破線で示すように、目的周
波数帯に含まれる他周波数帯の漏れ成分の絶対値は低周
波帯に行くに従い大きくなる傾向があることから、平滑
化フィルタのマスクサイズが大きくなる(高周波成分が
より低い周波数の成分を含むようになる)に従い、閾値
SL2を大きくすることで、より有効に他の周波数帯の
漏れ成分の影響を除去することができる。
【0059】再び図2を参照するに、成分変換部178
は変換曲線作成部176で作成される変換曲線G()に
より高周波成分f2(x、y)を変換する(ステップ1
015)。そして、変換された変換後の画像は高周波成
分加算部180に転送される。高周波成分加算部180
はメモリ120から原画像f(x、y)を呼び出すと共
に、変換された画像を原画像f(x、y)に加算する
(ステップ1020)。かかる一連の処理は数式3で表
される。
【0060】
【数3】
【0061】ここで、f3(x、y)を高周波成分加算
部180で得られた画像とする。
【0062】次いで、画像f3(x、y)は画像処理部
170より出力されメモリ120に格納される。よっ
て、オペレータはかかる画像(即ち、処理後画像)を表
示部140により表示又は印刷することで処理後画像を
確認できる。
【0063】なお、変換曲線G()、経験的に定めら
れ、予めメモリ120にデータとして格納されていても
良い。または、高周波成分の絶対値の平均値等、原画像
又はその周波数成分等の画像データを解析した結果に基
づいて、変換曲線G()の形状(閾値や傾き等)を決め
るようにプログラムされていてもよい。また、例えば、
オペレータがメモリ120に格納された関数G()を決
めるパラメータから所望のパラメータを選択すること
で、変換曲線G()が変換曲線作成部176によって作
成されてもよい。
【0064】以上の様に実施例1では所定の高周波成分
をその絶対値に基づいて抑制又は増加するように変換曲
線を作成している。よって、目的とする対象物を選択的
に強調できる効果がある。また、周波数帯毎に変換曲線
の閾値を変更することで他の周波数帯成分の漏れの影響
を有効に抑制することができる。さらに、目的とする周
波数帯の高周波成分を増加すると共に、他の周波数帯の
高周波成分を抑制するように変換曲線を作成すること
で、目的とする周波数帯の画像成分をより選択的に強調
することができる。また、所定絶対値(閾値)以上の絶
対値を有する高周波成分を抑制することで、エッジ部分
におけるオーバーシュート等のアーティファクトを抑制
することもできる。また、平滑化フィルタのマスクサイ
ズに応じて閾値を変更することで、効果的に他の周波数
帯成分の漏れの影響を防ぐことができる。更に、鮮鋭化
処理において強調したい所望の周波数帯域の画像成分を
画像から分離する際の周波数分離性(画像から分離した
画像成分に不所望の周波数帯域の画像成分が漏れ込むこ
と)に起因する問題を改善でき、目的とする対象物(所
望の周波数帯域の画像成分)を的確に強調することがで
きる。
【0065】実施例2 以下、図1、図15及び図16を参照するに、画像処理
部170の動作における実施例1に対する変形例である
画像処理方法ついて説明する。ここで、図15は、本発
明の一側面としての画像処理方法を示すフローチャート
である。図16は、図1に示す画像処理装置100の表
示部140で出力される表示内容を例示的に示す図であ
る。図16において、SLは閾値を示し、G1()及び
G2()は変換曲線を示す。領域141は閾値SLの値
を変更するためのパラメータが入力(表示)される領域
であり、領域142は変換曲線G()の曲線形(傾き)
を変更するパラメータが入力(表示)される領域であ
る。なお、図16では閾値SLは1つのみ表示されてい
るが、閾値を2つ入力(表示)するようにしてもよい。
領域143は処理後の画像を表示する領域であり周波数
処理が行われた後の画像を表示する。領域144は領域
143に表示された画像の所定水平方向ラインのプロフ
ァイルを示す領域である。マーク(矢印)145は領域
143に表示された画像において、領域144に表示さ
れるプロファイルに対応する水平ラインの位置を示す。
領域146はプロファイルを表示するラインの位置を入
力(表示)する領域である。これらの領域141、14
2及び146への入力は操作パネル130を操作し、入
力回路184を経由することで可能となる。例えば、領
域141、142及び146は、グラフィカル・ユーザ
・インターフェイスを構成し、領域141、142及び
146に示される左向き、右向き、上向き又は下向きの
矢印を、操作パネル130を構成するマウス又はタッチ
パネル等により操作(クリック又はタッチ)することに
より、各パラメータを入力又は変更可能に構成されてい
る。その際、入力又は変更されたパラメータは表示され
る変換曲線G()の閾値及び曲線形(傾き)やマーク
(矢印)145の位置に反映されるが、それとは別に領
域141、142及び146内の矢印の近傍に数値、記
号、図等により表示されるようにしてもよい。また、操
作パネル130から入力された情報(変換曲線G()の
閾値、曲線形(傾き)、又はプロファイルの位置)は入
力回路184を介して変換曲線作成部176を含む画像
処理装置170に反映される。また、表示部140への
画像及びその他情報の表示は表示回路182を介してC
PU108により行われる。
【0066】画像処理部170に転送された原画像は、
図2に示すステップと同様な工程を経て表示部140に
出力される(ステップ1100乃至1120)。なお、
ステップ1100乃至ステップ1115の工程は、図2
のステップ1000乃至1015で説明した工程と同じ
であって、ここでの詳細な説明は省略する。また、変換
曲線作成部178は予め定められているパラメータ(例
えば、経験上妥当であると考えられるパラメータがメモ
リ120に格納されている)を用い変換曲線G()を決
定し、成分変換部178にて高周波成分の変換処理を行
うものとする。
【0067】次に、オペレータは領域143に表示され
た画像をもとに、かかる画像の所望の水平ラインの位置
を(操作パネル130を介し)領域146を利用して入
力する。そして、CPU110の制御のもと、入力され
た位置に対応するプロファイルが領域144に表示され
る(ステップ1125)。次に、オペレータはマスクサ
イズを変更したい場合には新たなマスクサイズを(操作
パネル130を介し)選択する(ステップ1130)。
オペレータにより指定されたマスクサイズ、又はメモリ
120に格納されている代表的なマスクサイズをもとに
閾値SLが決定され、変換曲線作成部176で作成され
た変換曲線を変換曲線G1()として表示部140に表
示する(ステップ1135乃至1140)。
【0068】CPU110は変換曲線G1()を表示す
ると共に、成分変換部178を介し高周波成分を変換曲
線G1()で変換する(ステップ1145)。変換され
た後の高周波成分(高周波画像)は高周波成分加算部1
80で原画像に加算される(ステップ1150)。次い
で、処理後画像は領域143に表示される(ステップ1
155)。そして、処理後画像の予め指示された位置の
プロファイルが領域144にも表示される(ステップ1
160)。
【0069】以上の工程により、領域143に表示され
た処理後画像、表示された変換曲線G1()、領域14
4に表示されたプロファイルを観察することによって、
オペレータは処理後画像の適正を判断することができ
る。なお、オペレータが変換曲線G()を変更したい場
合(ステップ1165)は、領域141及び142を利
用して操作パネル130を介し所望のパラメータを入力
する(ステップ1170)。ここで、領域141に表示
されたパラメータを変更すると、変換曲線作成部176
によって閾値SLが変更される。また、領域142に表
示されたパラメータを変更すると変換曲線作成部176
によって変換曲線G()の曲線形状(傾き等)が変更さ
れ、例えば傾きを減じた場合には曲線G1()が曲線G
2()となる。
【0070】そして、パラメータを変更した場合にはス
テップ1135乃至1165までの処理を繰り返す。か
かるステップを繰り返すことで、オペレータは所望の変
換曲線を得ることができる。また、マスクサイズを変更
したい場合には(ステップ1175)、ステップ113
0からの処理を繰り返す。マスクサイズを変更する必要
がなければ処理を終了する。
【0071】実施例2では処理後画像、画像のプロファ
イル、変換曲線形を表示することで処理の結果と変換曲
線との関係を視覚的又は直観的に把握することができ
る。さらに、パラメータを変更することができるため、
適切な又は所望の処理結果を得るようにパラメータを変
更することができる。また、実施例1で述べた効果を奏
することは言うまでもない。
【0072】実施例3 以下、図17、図20及び図21を参照するに、画像処
理部170aの動作、即ち、本発明の画像処理装置10
0aの画像処理方法について説明する。ここで、図20
は、本発明の一側面としての画像処理方法を示すフロー
チャートである。図21は、図20に示す周波数帯毎の
画像成分を変換する処理(ステップ1205)を示すフ
ローチャートである。
【0073】上述したように、画像処理部170aは前
処理部150より出力された画像(原画像)をCPUバ
ス105を介し授受する。より特定的には、前処理部1
50より出力された原画像は特徴抽出部179に転送さ
れる。まず、画像処理部170aにおいて、特徴抽出部
179が原画像から階調変換のための特徴量を算出する
(ステップ1200)。このような方法として、例えば
人体の特定領域の画素値平均を特徴量とするなどの方法
があり、このような方法として特開平11−09638
0、特開平2000−099708等が本出願人により
出願されている。
【0074】次に、原画像の周波数帯毎の画像成分を変
換する(ステップ1205)。そして、階調変換部18
1は特徴抽出部179で得られた特徴量に基づき、画像
成分が変換された後の画像の階調変換を行う(ステップ
1210)。そして階調変換された画像は表示装置14
0により表示又は印刷される。
【0075】次にステップ1205の処理を図21に従
い説明する。
【0076】図20に示すステップ1205に移行する
と、原画像はDWT変換部171に転送される。DWT
変換部171は転送された原画像を表す画像信号に対し
て2次元の離散ウェーブレット変換処理を行い、複数の
周波数帯毎の画像成分(変換係数)を出力する。メモリ
120に記憶された画像データはDWT変換部171に
より順次読み出されて変換処理され、再びメモリ120
に書きこまれる。本実施の形態におけるDWT変換部1
71において、入力された画像信号は遅延素子およびダ
ウンサンプラの組み合わせにより、偶数アドレスおよび
奇数アドレスの信号に分離され、2つのフィルタpおよ
びuによりフィルタ処理が施される。図18に示すs及
びdは、各々1次元の画像信号に対して1レベルの分解
を行った際のローパス係数およびハイパス係数を表して
おり、次に示す数式4及び数式5に従って計算される。
【0077】
【数4】
【0078】
【数5】
【0079】ただし、x(n)は変換対象となる画像信
号である。
【0080】以上の処理により、画像信号に対する1次
元の離散ウェーブレット変換処理が行われる。2次元の
離散ウェーブレット変換は、1次元の変換を画像の水平
・垂直方向に対して順次行うものであり、その詳細は公
知であるのでここでは説明を省略する。図22に示すよ
うに、画像信号は異なる周波数帯域の画像成分HH1、
HL1、LH1、...、LLに分解される(ステップ
1206)。ここで、図22は、2次元の離散ウェーブ
レット変換処理により得られる2レベルの変換係数群の
構成例を示す図である。図22においてHH1、HL
1、LH1、...、LL等(以下サブバンドと呼ぶ)
が周波数帯毎の画像成分を示す。
【0081】次に、変換曲線作成部173は、例えば、
図23乃至図25に示す変換曲線を作成する(ステップ
1207)。ここで、図23乃至図25は、変換曲線作
成部173で作成された変換曲線の例を示す図であり、
横軸が入力画像成分、縦軸が出力画像成分を示す。図2
3において、閾値SL3及び閾値SL4は各々異なる周
波数帯域に対して設定した閾値を示している。閾値は、
例えば、SL3<SL4の関係があり、より低い周波数
帯域の画像成分を変換する場合にはより閾値が大きくな
るように構成されている。
【0082】このような変換曲線では、画像成分が小さ
な領域の成分を0とすると共に、それ以外の画像成分を
増加するものである。例えば、図26に示す様に、鮮鋭
度を上げたい対象物(血管等)の主な周波数成分が中周
波数成分であり、鮮鋭度を落としたい非対象物(ノイズ
等)の主な周波数成分が高周波成分である場合、中周波
数に対応する画像成分を増幅することで血管などの鮮鋭
度を選択的に上げることが可能となる。ここで、図26
は、対象物と非対象物の各周波数分布を示した図であ
り、横軸が周波数、縦軸が対象物及び非対象物の各周波
数成分の出現頻度を示す。
【0083】しかし、非対象物と対象物の各周波数分布
に重なりあう領域が存在し、中周波数帯の画像成分を一
律に増ずると非対象物の強調も起きてしまう。但し、図
27の実線に示す様に、中周波数帯の画像成分の絶対値
として対象物の画像成分は絶対値の大きい方に偏在し、
非対象物の画像成分は絶対値の小さい方に偏在する。こ
こで、図27は、中周波数帯の画像成分中の対象物と非
対象物の各成分の出現頻度の関係を示す図であり、横軸
が画像成分の絶対値、縦軸が出現頻度である。このた
め、対象物の成分を増加する場合に、ノイズ成分が多く
含有される絶対値の範囲(例えば、閾値SL1未満の範
囲)に対し対象物の成分を多く含む絶対値の範囲の画像
成分を相対的に増加することで、非対象物の成分を強調
することなく対象物の鮮鋭度を選択的に上げることがで
きる。
【0084】さらに、図27の破線に示す様に、より低
い周波帯の画像成分では画像成分中に含まれる非対象物
の成分の絶対値が大きくなる傾向があるため、より低周
波になるに従い閾値SL1を大きくすることにより、非
対象物の成分をより効果的に抑えることができる。この
ことより、閾値としては図28に示す様に、より低い周
波数帯になるに従い閾値を大きくすることが、対象物を
選択的に強調するには効果的である。ここで、図28
は、周波数帯の周波数と閾値との関係を示し、横軸が周
波数、縦軸が閾値の値を示している。なお、図28では
線形的に閾値を変更しているが、非線形であってももち
ろんよい。
【0085】また、例えば中周波数帯の画像成分には、
低周波数帯、高周波数帯の成分も含まれる(即ち、他の
周波数帯成分の漏れ)。この場合に、目的とする周波数
帯(中周波帯)の画像成分に含まれる他の周波数帯(低
周波数帯、高周波数帯)の成分の絶対値は比較的小さな
値として現れる。そこで、絶対値の小さい範囲(例えば
閾値SL3以下の範囲)の画像成分の絶対値を相対的に
抑制することで、他周波数帯成分の漏れの影響を抑制で
き、目的とする周波数帯の成分を強調することができ
る。また、目的周波数帯に含まれる他周波数帯の成分の
絶対値はより低い周波帯に行くに従い大きくなる傾向が
ある。そこで、より低い周波帯に行くに従い、閾値Th
を大きくすることで、より有効に他の周波数帯の成分の
影響、すなわち他周波数帯成分の漏れの影響を除去又は
抑制することができる。
【0086】図24に示す変換曲線は図23に示すそれ
と同様の効果を狙ったものであり、閾値SL3、SL4
以下の成分(出力係数)を0としている点は同様である
が、かかる変換曲線は変換曲線の立ち上がりが出力係数
0から始まるようになっている点において異なる。
【0087】図25に示す変換曲線は画像成分の周波数
帯域に応じて変換曲線の傾きが異なるものであるが、図
23と同様な効果を有する。かかる変換曲線は目的とす
る周波数帯の成分は増加するように変換する曲線形(傾
きが相対的に大)とし、他の周波数帯の画像成分は減じ
る様に変換する曲線形(傾きが相対的に小)としてい
る。このことにより、目的とする周波数帯の画像成分は
相対的に強調され、目的外の周波数帯の画像成分は相対
的に抑制される。従って、目的周波数帯の画像成分のみ
を増加する場合より有効に目的とする周波数帯の画像成
分が強調される効果を有するものである。目的とする周
波数帯の成分の振幅は増加し、目的外の周波数帯の振幅
は減少するので相対的に目的とする周波数帯の成分が強
調されるためである。
【0088】なお、DWT変換部171で原画像を各周
波数帯毎の画像成分に分解した(ステップ1206)
後、変換曲線作成部173は各周波数帯における画像成
分の絶対値の平均値を算出し、この平均値に基づき閾値
を決定するステップを設けても良い。例えば、各周波数
帯における画像成分の絶対値の平均値の10%を閾値と
定めても良い。かかるステップを設けることで、周波数
帯毎の画像成分に基づき閾値を決定することができ、入
力画像の成分の大きさの変動の影響を受けず、安定した
処理を行うことができる。
【0089】一方、変換曲線作成部173は例えば図2
9乃至図31に示す変換曲線を作成する(ステップ12
07)ことも考えられる。ここで、図29乃至図31
は、変換曲線作成部173で作成された変換曲線の例を
示す図であり、横軸が入力画像成分(入力係数)、縦軸
が出力画像成分(出力係数)を示す。図29において、
閾値SL3又は閾値SL4は周波数帯毎に異なる閾値に
構成したものである。閾値には、例えば、SL3<SL
4の関係があり、より低い周波数の画像成分を変換する
場合に、より閾値が大きくなるように構成されている。
【0090】このような変換曲線は、大きな絶対値を有
する画像成分の絶対値を0とすると共に、それ以外の画
像成分の絶対値を増加するものである。
【0091】例えば、図32に示すように、鮮鋭度を上
げたい対象物(血管等)の主な周波数成分が中周波数成
分であり、鮮鋭度を落としたい非対象物(臓器等)の主
な周波数成分が低周波成分である場合に、中周波数に対
応する画像成分を増幅することで血管などの鮮鋭度を上
げることが可能となる。ここで、図32は、対象物と非
対象物の各画像成分の周波数分布を示した図であり、横
軸が周波数、縦軸が対象物及び非対象物の各周波数成分
の出現頻度を示す。
【0092】しかし、図32に示すように、非対象物と
対象物の周波数成分に重なりあう領域が存在し、中周波
数帯の画像成分を一律に増ずると非対象物の強調も起き
てしまう。但し、図33の実線に示すように、中周波数
帯の画像成分の絶対値として対象物の画像成分は絶対値
の小さい方に偏在し、非対象物の画像成分は絶対値の大
きい方に偏在する。ここで、図33は、中周波帯の画像
成分中の対象物と非対象物の出現頻度の関係を示す図で
あり、横軸が画像成分の値(絶対値)、縦軸が出現頻度
である。このため、対象物の中周波数帯の成分を増加す
る場合に、非対象物の成分が多く含有される絶対値の範
囲(例えば閾値SL3以上の範囲)の絶対値を有する画
像成分の絶対値を相対的に減ずるように変換することで
非対象物の画像成分の強調を抑制することができる。ま
た、対象物の成分を多く含む絶対値範囲の絶対値を有す
る画像成分を相対的に増加させることで対象物の鮮鋭度
を選択的に上げることができる。
【0093】さらに、図33の破線で示すように、より
低い周波数帯の画像成分では画像成分中に含まれ非対称
物の成分の絶対値が大きくなる傾向がある。従って、よ
り低周波になるに従い閾値SL3を大きくすることによ
り、非対象物の成分を相対的に抑制し、対象物の成分を
相対的に強調することができる。このことより、閾値と
しては図34に示す様に、より低い周波数帯になるに従
い閾値を大きくすることが、対象物を選択的に強調する
には効果的である。ここで、図34は、周波数帯の周波
数と閾値との関係を示し、横軸が周波数、縦軸が閾値の
値を示している。なお、図34では線形的に閾値を変更
しているが、非線形であってもよい。
【0094】また、画像中のエッジ部分(例えば人体で
は骨の縁、金属片(ペースメーカー等)の縁)を構成す
る周波数成分は低周波から高周波まで広い範囲の周波数
成分を含む。また、その周波数成分の振幅が大きいた
め、対応する画像成分(周波数係数)も比較的大きな値
を有する。そこで、図29乃至図31に示すような変換
曲線を用いて、所定の周波数帯域の画像成分の所定絶対
値(閾値)以上の絶対値を有する画像成分を相対的に抑
制する又は不変とすることは画像中のエッジ部分にオー
バーシュートなどのアーティファクトが生じるのを抑制
することにもなる。
【0095】図30に示す変換曲線は一定閾値SL3以
下の画像成分を強調している。また、一定閾値SL3以
上の画像成分は一定閾値SL3以下の画像成分に比し相
対的に抑制するか又は値を変換しないような曲線形にし
たものである。これにより対象物の画像成分は強調され
るものの、非対象物の画像成分は相対的に抑制され、又
は変更されないようにすることができる。また、上述し
たようにエッジ部分に相当する画像成分を変更しないの
で、オーバーシュートなどのアーティファクトが生じる
のを抑制することもできる。さらに、図30の変換曲線
を、その両端を除き微分連続(傾きが連続)とすること
により、偽輪郭の発生を確実に抑えるなど、より自然な
処理後画像を得ることができる。
【0096】図31に示す変換曲線は図29と同様の効
果を奏するものであるが、変換曲線が0に減ずるまでに
一定幅を持たせたものである。これにより、閾値以上の
絶対値を有する画像成分をなだらかに減少することが可
能であり、処理後の画像がより自然となる効果がある。
また、変換曲線を、その両端を除き微分連続(傾きが連
続)としており、偽輪郭の発生等を抑制し、より自然な
処理後画像を得ることができる。
【0097】以上のことにより、本来目的とする周波数
帯の画像成分は相対的に強調され、補以来目的外の周波
数帯の画像成分は相対的に抑制される。よって、所定の
周波数帯の画像成分を一律に増加する場合より、目的と
する対象物の画像成分を効果的に強調することができ
る。目的とする周波数帯の画像成分の振幅は相対的に増
加し、目的外の周波数帯の画像成分の振幅は相対的に減
少するので、相対的に目的とする周波数帯すなわち目的
とする対象物の画像成分が強調されるためである。
【0098】なお、DWT変換部171が原画像を周波
数帯毎の画像成分に分解(ステップ1206)した後、
変換曲線作成部173は各周波数帯における画像成分の
絶対値の平均値を算出し、この平均値に基づき閾値を決
定するステップを設けても良い。例えば、各周波数帯に
おける画像成分の絶対値の平均値の70%を閾値と定め
ても良い。かかるステップを設けることで、周波数帯毎
の画像成分に基づき閾値を決定することができるため、
入力画像の成分の大きさの変動の影響を受けず、安定し
た処理を行える効果がある。
【0099】次に、成分変換部175は変換曲線作成部
173で作成された変換曲線に従い、所定の周波数帯の
画像成分を変換する(ステップ1208)。なお、変換
された各周波数帯の画像成分はメモリ120に格納され
る。そして、逆DWT変換部177は逆離散ウェーブレ
ット変換を以下のように行う(ステップ1209)。メ
モリ120に記憶された変換された画像成分は逆DWT
変換部177により順次読み出されて逆ウェーブレット
変換処理が行われ、再びメモリ120に書きこまれる。
図19に示すように、入力された画像成分はuおよびp
の2つのフィルタ処理を施され、アップサンプリングさ
れた後に重ね合わされて画像信号x’が出力される。こ
れらの処理は次に示す数式6及び数式7により行われ
る。
【0100】
【数6】
【0101】
【数7】
【0102】以上の処理により、変換後の係数に対する
1次元の逆離散ウェーブレット変換処理が行われる。2
次元の逆離散ウェーブレット変換は、1次元の逆変換を
画像の水平・垂直方向に対して順次行うものであり、そ
の詳細は公知であるのでここでは説明を省略する。
【0103】以上説明したように、実施例3では周波数
帯毎に画像成分を、その絶対値に応じて抑制又は強調す
るように変換曲線を作成した。これによって、目的とす
る対象物を選択的に強調することができる。また、周波
数帯毎に変換曲線の閾値を変更することで、他周波数帯
の画像成分の漏れの影響を効果的に抑制することができ
る。さらに、目的とする周波数帯における目的とする対
象物の画像成分を相対的に強調し、他の周波数帯の漏れ
画像成分を相対的に抑制するように変換曲線を作成する
ことにより、目的とする周波数帯の画像成分における目
的とする対象物の画像成分を選択的に強調することがで
きる。
【0104】実施例4 以下、図17、図35及び図36を参照するに、画像処
理部170aの動作における実施例3に対する変形例で
ある画像処理方法ついて説明する。ここで、図35は、
本発明の一側面としての画像処理方法を示すフローチャ
ートである。図36は、図17に示す画像処理装置10
0aの表示部140により出力される表示内容を例示的
に示す図である。なお、図36に示す表示部140によ
り出力される表示内容の構成は図16に示すそれと類似
した構成であり、重複する説明は省略する。図36にお
いて、G3()及びG4()は変換曲線を示し、SLは
閾値を示す。
【0105】画像処理部170aに転送された原画像
は、まず特徴量抽出部179において解析され、その結
果原画像から階調変換のための特徴量が算出される(ス
テップ1300)。次に、原画像はDWT変換部171
で離散ウェーブレット変換処理され、複数の周波数帯の
画像成分に分解される(ステップ1305)。なお、ス
テップ1305は図21に示すステップ1206と同様
であるため、説明を省略する。原画像が各周波数帯毎の
画像成分に分解された後、変換曲線作成部173は予め
決められたパラメータに従い周波数帯毎の変換曲線を作
成する(ステップ1310)。次に、成分変換部175
はかかる変換曲線で画像成分を変換する。変換された後
の各周波数帯の画像成分は逆DWT変換部177で逆離
散ウェーブレット変換され、画像が生成される。次に、
生成された画像を階調変換部181で階調変換し、かか
る変換された画像を領域143に表示する(ステップ1
315)。
【0106】次に、領域143に表示された画像の中の
プロファイルを見たい水平ラインの位置を特定するパラ
メータを、領域146を利用して(操作パネルを介し
て)入力し、プロファイルの表示される水平ラインの位
置を決定する。入力されたパラメータに基づいて所望の
プロファイルを領域144に表示する(ステップ132
0)。そして、オペレータは変換曲線を変更したいサブ
バンドを選択する(ステップ1325)。変換曲線作成
部173は選択されたサブバンドに対応する変換曲線を
作成し、作成された変換曲線を変換曲線G3()として
表示する(ステップ1330乃至1335)。次に、か
かるサブバンドに対応する画像成分を変換曲線で変換す
る(1340)。変換された各周波数帯毎の画像成分は
逆DWT変換部177で逆離散ウェーブレット変換さ
れ、画像が合成される(ステップ1345)と共に、か
かる画像は階調変換部181で階調変換される。そし
て、階調変換後の画像は処理後画像として領域143に
表示される(ステップ1350)。そして、処理後画像
の予め指示された位置のプロファイルが領域144に表
示される(ステップ1355)。
【0107】この領域143に表示される処理後画像、
表示される変換曲線G3()、領域144に表示される
プロファイルをオペレータは観察することができる。な
お、オペレータが変換曲線G3()を変更したい場合
(ステップ1360)は、領域141及び142を利用
してパラメータを入力する(ステップ1365)。例え
ば、領域141の値を変更すると閾値SLが変更され
る。また、領域142の値を変更すると、変換曲線G3
()の曲線形(傾き等)が変更される。例えば、傾きを
減じた場合には変換曲線G3()がG4()に変更され
る。パラメータを変更した場合にはステップ1330乃
至ステップ1355までの処理が繰り返される。
【0108】そして、他のサブバンドのパラメータを変
更したい場合(ステップ1370)には、ステップ13
25からの処理を繰り返し、全ての変換曲線形を変更す
る必要がなければ処理を終了する。
【0109】実施例4では処理後画像、画像のプロファ
イル、変換曲線形を表示することで処理結果と変換曲線
(変換関数)との関係を視覚的又は直観的に把握するこ
とができる。さらに、変換曲線(変換関数)を変更する
ことができるため、適切な又は所望の処理結果を得るこ
とができる。
【0110】
【発明の効果】以上説明したように本発明では、対象画
像のうち目的とする周波数帯又は対象物の画像成分を効
果的又は選択的に変換する画像処理装置及び方法、並び
に、コンピュータプログラム及び記憶媒体を提供するこ
とができる。
【0111】本発明の画像処理装置は、変換曲線作成部
を有し、かかる作成部は所定の閾値を変更可能に作成す
ることができる。よって、かかる画像処理装置は高周波
成分を所定の閾値を境に異なる関数で変換することがで
き、目的とする周波数帯以外の高周波成分を抑制するこ
とができる。また、かかる画像処理装置は前記所定の閾
値を設定することで、所望の周波数帯のみを選択的に変
換することができる。よって、かかる画像処理装置は所
望の高周波成分のみを効果的に強調することができ、画
像の鮮鋭化に寄与する。その結果、本発明の画像処理装
置はより自然な感じの処理後画像を得ることができる。
【0112】更に、かかる画像処理装置は前記作成部に
より作成される所定の閾値を表示可能な表示部を有して
も良い。かかる表示部は処理後画像及び所定の関数を表
示することで、オペレータが視覚的に確認しながら画像
処理を行うことができる。また、かかる画像処理装置は
前記作成部により作成される前記所定の閾値を変更する
入力部を有してもよい。かかる入力部は前記閾値を変更
することができるので、適切な処理後画像が得られるよ
うに前記関数を変更することができる。これにより、か
かる処理後画像はオペレータの所望とする通りの画像を
得ることができる。例えば、かかる画像を診断に供する
場合、医師の診断の正確性を向上させることができる。
【図面の簡単な説明】
【図1】 本発明の画像処理装置を示すブロック図であ
る。
【図2】 本発明の一側面としての画像処理方法を示す
フローチャートである。
【図3】 高周波成分を変換する変換曲線の例を示した
図である。
【図4】 高周波成分を変換する変換曲線の例を示した
図である。
【図5】 高周波成分を変換する変換曲線の例を示した
図である。
【図6】 対象物と非対象物の周波数成分を示す図であ
る。
【図7】 中周波数帯の高周波成分における対象成分と
非対象成分の出現頻度を示す図である。
【図8】 高周波成分を変換する変換曲線の例を示した
図である。
【図9】 高周波成分を変換する変換曲線の例を示した
図である。
【図10】 高周波成分を変換する変換曲線の例を示し
た図である。
【図11】 高周波成分を変換する変換曲線の例を示し
た図である。
【図12】 対象物と非対象物の周波数成分を示す図で
ある。
【図13】 中周波数帯の高周波成分における対象成分
と非対象成分の出現頻度を示す図である。
【図14】 閾値とマスクサイズの関係を示す図であ
る。
【図15】 本発明の別の画像処理方法を示すフローチ
ャートである。
【図16】 図1に示す画像処理装置の表示部に出力さ
れる表示内容を例示的に示す図である。
【図17】 本発明の別の画像処理装置を示すブロック
図である。
【図18】 図17に示すDWT変換部を構成可能な回
路を例示的に示した図である。
【図19】 図17に示す逆DWT変換部を構成可能な
回路を例示的示した図である。
【図20】 本発明の別の画像処理方法を示すフローチ
ャートである。
【図21】 図20に示す周波数帯毎の画像成分を変換
する処理を示すフローチャートである。
【図22】 2次元の離散ウェーブレット変換処理によ
り得られる2レベルの変換係数群の構成例を示す図であ
る。
【図23】 変換曲線作成部で作成された変換曲線及び
閾値の例を示す図である。
【図24】 変換曲線作成部で作成された変換曲線及び
閾値の例を示す図である。
【図25】 変換曲線作成部で作成された変換曲線及び
閾値の例を示す図である。
【図26】 対象物と非対象物の周波数分布を示した図
である。
【図27】 中周波数帯の画像成分中の対象物と非対象
物の出現頻度の関係を示す図である。
【図28】 周波数帯と閾値の関係を示した図である。
【図29】 変換曲線作成部で作成された変換曲線及び
閾値の例を示す図である。
【図30】 変換曲線作成部で作成された変換曲線及び
閾値の例を示す図である。
【図31】 変換曲線作成部で作成された変換曲線及び
閾値の例を示す図である。
【図32】 対象物と非対象物の周波数分布を示した図
である。
【図33】 中周波帯の画像成分中の対象物と非対象物
の出現頻度の関係を示す図である。
【図34】 周波数帯と閾値の関係を示す図である。
【図35】 本発明の別の画像処理方法を示すフローチ
ャートである。
【図36】 図17に示す画像処理装置の表示部に出力
される表示内容を例示的に示す図である。
【符号の説明】
100 画像処理装置 100a 画像処理装置 110 CPU 120 メモリ 130 操作パネル 140 表示部 150 前処理部 160 データ収集部 170 画像処理部 170a 画像処理部 190 撮像部
───────────────────────────────────────────────────── フロントページの続き (72)発明者 佐藤 眞 東京都大田区下丸子3丁目30番2号 キヤ ノン株式会社内 Fターム(参考) 5B057 AA07 BA03 CA02 CA08 CA16 CB02 CB08 CB16 CC01 CD14 CE03 CE05 CE06 5C077 LL02 LL19 PP01 PP17

Claims (22)

    【特許請求の範囲】
  1. 【請求項1】 対象画像から当該対象画像の低周波成分
    と少なくとも1つの高周波成分とを作成する作成部と、 前記少なくとも1つの高周波成分のうち所定の高周波成
    分における所定の絶対値範囲の絶対値を有する要素を強
    調又は抑制する成分変換部と、 前記所定の絶対値範囲を決める少なくとも1つの閾値を
    設定する設定部と、 前記成分変換部により変換された後の前記少なくとも1
    つの高周波成分と、前記対象画像又は前記低周波成分と
    を用いて画像を生成する画像生成部とを有することを特
    徴とする画像処理装置。
  2. 【請求項2】 前記成分変換部が前記所定の高周波成分
    を変換する際に用いる変換関数を表示する表示部を有す
    る請求項1記載の画像処理装置。
  3. 【請求項3】 前記成分変換部が前記所定の高周波成分
    を変換する際に用いる変換関数を変更するための入力部
    を有する請求項1記載の画像処理装置。
  4. 【請求項4】 強調する対象の画像成分が強調する対象
    でない画像成分に比べて相対的に高い周波数側に分布す
    る場合、前記成分変換部は前記所定の高周波成分の所定
    絶対値以上の絶対値を有する要素を相対的に減ずるよう
    に前記所定の高周波成分を変換する請求項1記載の画像
    処理装置。
  5. 【請求項5】 前記成分変換部は前記所定の高周波成分
    の所定絶対値以上の絶対値を有する要素を相対的に減ず
    るように前記所定の高周波成分を変換する請求項1記載
    の画像処理装置。
  6. 【請求項6】 前記少なくとも1つの閾値は予め決めら
    れた固定値である請求項1乃至5記載の画像処理装置。
  7. 【請求項7】 前記成分変換部は前記少なくとも1つの
    閾値を少なくとも前記対象画像に基づいて決定する請求
    項1乃至5記載の画像処理装置。
  8. 【請求項8】 前記画像生成部は前記成分変換部により
    変換された後の前記少なくとも1つの高周波成分と前記
    対象画像とを用いて画像を生成すると共に、前記成分変
    換部は前記所定の高周波成分の前記絶対値範囲外の絶対
    値を有する要素が実質的に0となるように前記所定の高
    周波成分を変換する請求項1記載の画像処理装置。
  9. 【請求項9】 前記画像生成部は前記成分変換部により
    変換された後の前記少なくとも1つの高周波成分と前記
    低周波成分とを用いて画像を生成すると共に、前記成分
    変換部は前記所定の高周波成分の前記絶対値範囲外の絶
    対値を有する要素が実質的に不変となるように前記所定
    の高周波成分を変換する請求項1記載の画像処理装置。
  10. 【請求項10】 前記成分変換部が前記所定の高周波成
    分を変換する際に用いる変換関数はその傾きが連続であ
    る請求項1記載の画像処理装置。
  11. 【請求項11】 前記設定部は前記閾値を、前記所定の
    高周波成分がより低い周波数の画像成分を含むほど大き
    くする請求項1記載の画像処理装置。
  12. 【請求項12】 前記成分変換部は、前記作成部により
    作成される複数の高周波成分の各々の所定の絶対値範囲
    の絶対値を有する要素を強調又は抑制する請求項1記載
    の画像処理装置。
  13. 【請求項13】 強調する対象の画像成分が強調する対
    象でない画像成分に比して相対的に高い周波数側に分布
    する場合、前記成分変換部は前記複数の高周波成分の各
    々の所定絶対値以上の絶対値を有する要素を相対的に減
    ずるように前記複数の高周波成分を変換する請求項12
    記載の画像処理装置。
  14. 【請求項14】 前記成分変換部は前記複数の高周波成
    分の各々の所定絶対値以上の絶対値を有する要素を相対
    的に減ずるように前記複数の高周波成分を変換する請求
    項12記載の画像処理装置。
  15. 【請求項15】 前記複数の高周波成分の各々の前記所
    定の絶対値範囲を決める少なくとも1つの閾値は予め決
    められた固定値である請求項12記載の画像処理装置。
  16. 【請求項16】 前記成分変換部は前記複数の高周波成
    分の各々の前記所定の絶対値範囲を決める少なくとも1
    つの閾値を少なくとも前記対象画像に基づいて決定する
    請求項12記載の画像処理装置。
  17. 【請求項17】 前記成分変換部は前記複数の高周波成
    分の各々の前記所定の絶対値範囲を決める少なくとも1
    つの閾値を当該高周波成分の画像成分の平均値より算出
    する請求項12記載の画像処理装置。
  18. 【請求項18】 対象画像から当該対象画像の低周波成
    分と少なくとも1つの高周波成分とを作成するステップ
    と、 前記少なくとも1つの高周波成分のうち所定の高周波成
    分における所定の絶対値範囲の絶対値を有する要素を強
    調又は抑制するステップと、 前記所定の絶対値範囲を決める少なくとも1つの閾値を
    設定するステップと、 前記強調又は抑制ステップにより変換された後の前記少
    なくとも1つの高周波成分と、前記対象画像又は前記低
    周波成分とを用いて画像を生成するステップとを有する
    ことを特徴とする画像処理方法。
  19. 【請求項19】 請求項1乃至17のうちいずれか一項
    記載の画像処理装置の機能をコンピュータに実現させる
    ためのプログラムを記憶したコンピュータ読取可能な記
    憶媒体。
  20. 【請求項20】 請求項18記載の画像処理方法の処理
    ステップをコンピュータに実行させるためのプログラム
    を記憶したコンピュータ読取可能な記憶媒体。
  21. 【請求項21】 請求項1乃至17のうちいずれか一項
    記載の画像処理装置の機能をコンピュータに実現させる
    ためのプログラム。
  22. 【請求項22】 請求項18記載の画像処理方法の処理
    ステップをコンピュータに実行させるためのプログラ
    ム。
JP2001117548A 2001-04-16 2001-04-16 画像処理装置及び方法、並びに、コンピュータプログラム及び記憶媒体 Expired - Fee Related JP3793036B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001117548A JP3793036B2 (ja) 2001-04-16 2001-04-16 画像処理装置及び方法、並びに、コンピュータプログラム及び記憶媒体

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2001117548A JP3793036B2 (ja) 2001-04-16 2001-04-16 画像処理装置及び方法、並びに、コンピュータプログラム及び記憶媒体

Publications (2)

Publication Number Publication Date
JP2002314816A true JP2002314816A (ja) 2002-10-25
JP3793036B2 JP3793036B2 (ja) 2006-07-05

Family

ID=18968075

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001117548A Expired - Fee Related JP3793036B2 (ja) 2001-04-16 2001-04-16 画像処理装置及び方法、並びに、コンピュータプログラム及び記憶媒体

Country Status (1)

Country Link
JP (1) JP3793036B2 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004127064A (ja) * 2002-10-04 2004-04-22 Konica Minolta Holdings Inc 画像処理方法、画像処理装置、画像処理プログラム及び画像記録装置
US7995148B2 (en) 2006-12-22 2011-08-09 Kabushiki Kaisha Toshiba Image processing apparatus and image processing method
US20220130053A1 (en) * 2019-05-30 2022-04-28 Panasonic Intellectual Property Management Co., Ltd. Image processing method, image processing device, and recording medium

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004127064A (ja) * 2002-10-04 2004-04-22 Konica Minolta Holdings Inc 画像処理方法、画像処理装置、画像処理プログラム及び画像記録装置
US7995148B2 (en) 2006-12-22 2011-08-09 Kabushiki Kaisha Toshiba Image processing apparatus and image processing method
US20220130053A1 (en) * 2019-05-30 2022-04-28 Panasonic Intellectual Property Management Co., Ltd. Image processing method, image processing device, and recording medium

Also Published As

Publication number Publication date
JP3793036B2 (ja) 2006-07-05

Similar Documents

Publication Publication Date Title
JP3833177B2 (ja) 画像処理装置、画象処理方法、記憶媒体及びプログラム
US6813335B2 (en) Image processing apparatus, image processing system, image processing method, program, and storage medium
EP1501049B1 (en) Method for rendering digital radiographic images for display based on independent control of fundamental image quality parameters
JP4413504B2 (ja) 医用画像処理装置、医用画像処理方法および医用画像処理プログラム
JP3690882B2 (ja) 画像の強調処理方法および装置
EP1208536B1 (en) Contrast adjustment by computing a local mean density
JP4679710B2 (ja) ノイズ抑制処理装置並びに記録媒体
US8428324B2 (en) Method and apparatus for processing digital mammographic images with parallel multi-scale decomposition and recombining in an optimized multi-scale reconstruction
JP3700804B2 (ja) 画像処理方法および装置
JP2009516882A (ja) ディジタル画像を強調する方法及びシステム
JPH1075395A (ja) 画像処理方法および装置
WO2005110232A1 (ja) 画像処理装置及びその方法
JP2001057677A (ja) 画像処理方法および装置並びに記録媒体
JPH1075364A (ja) 画像のダイナミックレンジ圧縮処理方法および装置
JP4486784B2 (ja) ノイズ低減方法
US7024036B2 (en) Image processing apparatus, image processing method, storage medium, and program
US8374417B2 (en) Image processing method and radiographic apparatus using the same
JP3730872B2 (ja) 画像処理装置及び画像処理プログラム
JP4197392B2 (ja) ノイズ抑制処理装置並びに記録媒体
JPH1063839A (ja) 画像処理方法および装置
JP3793036B2 (ja) 画像処理装置及び方法、並びに、コンピュータプログラム及び記憶媒体
JP3783116B2 (ja) 放射線画像強調処理方法および装置
JP4235622B2 (ja) 画像処理装置、画像処理システム、画像処理方法、プログラム及び記憶媒体
JP2004086598A (ja) 診断支援装置
JP2004173078A (ja) 画像処理装置

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20051122

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20051129

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060120

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: 20060328

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060406

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: 20090414

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100414

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110414

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120414

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130414

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130414

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140414

Year of fee payment: 8

LAPS Cancellation because of no payment of annual fees