JP2000163570A - 画像処理装置、方法及びコンピュータ読み取り可能な記憶媒体 - Google Patents

Abstract

(57)【要約】 【課題】 原画像の高周波成分を保持しながらダイナミ
ックレンジ圧縮を行う画像処理装置を得る。 【解決手段】 微係数計算回路１１２ａは、階調変換回
路１１２ｅで階調変換に用いる階調変換曲線の微係数に
基づく係数を計算する。平滑化画像作成回路１１２ｂ
は、原画像の平滑化画像（低周波画像）を作成し、高周
波成分作成回路１１２ｃは、上記作成された平滑化画像
と上記原画像との差分を高周波成分として計算する。高
周波成分足し込み回路１１２ｄは、上記計算された高周
波成分を、上記計算された係数に基づいて上記原画像又
は平滑化画像に足し込み、階調変換回路１１２ａは、上
記高周波成分を足し込まれた画像を階調変換する。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、原画像の高周波成
分を保持した状態で階調変換を行うための画像処理装
置、方法及びそれらに用いられるコンピュータ読み取り
可能な記憶媒体に関するものである。

【０００２】

【従来の技術】例えば、Ｘ線胸部画像は、Ｘ線が透過し
やすい肺野の画像、及びＸ線が非常に透過しにくい縦隔
部の画像より構成されるため、画素値の存在するレンジ
が非常に広い。このため、肺野及び縦隔部の両方を同時
に観察することが可能なＸ線胸部画像を得ることは困難
であるとされてきた。

【０００３】この問題を回避する方法として、従来より
ＳＰＩＥＶｏ１．６２６Ｍｅｄｉｃｉｎｅ ＸＩＶ／Ｐ
ＡＣＳＩＶ（１９８６）に記載される方法がある。この
方法は、処理後の画素値Ｓ_D 、オリジナル画素値（入力
画素値）Ｓ_org 、オリジナル画像（入力画像）の低周波
画像の画素値Ｓ_US、定数Ａ，Ｂ，Ｃを以て（例えばＡ＝
３、Ｂ＝０．７）、 Ｓ_D ＝Ａ〔Ｓ_org −Ｓ_US〕＋Ｂ〔Ｓ_US〕＋Ｃ・・・・（１） なる式（１）で表わされるものである。

【０００４】この方法は、高周波成分（第１項）、低周
波成分（第２項）の重み付けを変えることが可能で、例
えばＡ＝３、Ｂ＝０．７では高周波成分を強調し、かつ
全体のダイナミックレンジを圧縮する効果が得られるも
のである。この方法は５人の放射線医により、処理なし
画像と比較して診断に有効であるという評価が得られて
いる。

【０００５】また、特許第２，５０９，５０３号公報に
は、処理後の画素値Ｓ_D オリジナル画素値（入力画素
値）Ｓ_org 、オリジナル画像（入力画像）のＹ方向プロ
ファイルの平均プロファイルＰｙとＸ方向プロファイル
の平均プロファイルＰｘを以て、 Ｓ_D ＝Ｓ_org ＋Ｆ〔Ｇ（Ｐｘ，Ｐｙ）〕・・・・（２） なる式（２）で表わされる方法が記載されている。

【０００６】ここで、関数ｆ（ｘ）が有する特性につい
て説明すると、まず、「ｘ＞Ｄｔｈ」ではｆ（０）が
「０」となり、「０≦ｘ≦Ｄｔｈ」ではｆ（ｘ）が切片
を「Ｅ」、傾き「Ｅｍｔｈ」として単調減少するもので
あり、（３）式で示される。 Ｆ〔ｘ〕＝Ｅ−（Ｅ／Ｄｔｈ）Ｘ・・・・（３） Ｐｙ＝（ΣＰｙｉ）／ｎ・・・・・・・・（４） Ｐｘ＝（ΣＰｘｉ）／ｎ・・・・・・・・（５） 但し、（ｉ＝１〜ｎ）、Ｐｙｉ，Ｐｘｉはプロファイ
ル。そして例えば Ｇ（Ｐｘ，Ｐｙ）＝ｍａｘ（Ｐｘ，Ｐｙ）・・・（６） で示されるものである。

【０００７】この方法によれば、低周波画像の画素値で
Ｄｔｈ以下の濃度レンジが圧縮される。

【０００８】また、上記特許第２，５０９，５０３号許
公報の方法と同様な方法が、「日本放射線技術学会雑誌
第４５巻第８号１９８９年８月１０３０頁、阿南ほか」
に記載されている。この方法は、処理後の画素値Ｓ_D 、
オリジナル画素値（入力画素値）Ｓ_org 、オリジナル画
像（入力画像）をマスクサイズＭ×Ｍ画素で移動平均を
とった時の平均画素値Ｓ_US、単調減少関数ｆ（ｘ）を以
て、 Ｓ_D ＝Ｓ_org ＋ｆ（Ｓ_US）・・・・（７） Ｓ_US＝ΣＳ_org ／Ｍ² ・・・・・・（８） なる式（７）、（８）式で表わされるものである。

【０００９】この方法は、（２）式と低周波画像の作成
方法が異なるものである。即ち、（２）式では、１次元
データで低周波画像を作成していたのに対し、２次元デ
ータで低周波画像を作成するものである。この方法も低
周波画像の画素値でＤｔｈ以下の濃度値を圧縮するもの
である。

【００１０】また、特許第２，６６３，１８９号公報に
記載された方法は、処理後の画素値Ｓ_D 、オリジナル画
素値（入力画素値）Ｓ_org 、オリジナル画像（入力画
像）をマスクサイズＭ×Ｍ画素で移動平均をとった時の
平均画素値Ｓ_US、単調増加関数ｆ１（ｘ）を以て、 Ｓ_D ＝Ｓ_org ＋ｆ１（Ｓ_US）・・・・（９） Ｓ_US＝ΣＳ_org ／Ｍ² ・・・・・・・（１０） なる（９）、（１０）式で表わされるものである。

【００１１】ここで、関数ｆ１（ｘ）が有する特性につ
いて説明すると、まず、「ｘ＜Ｄｔｈ」ではｆ１（ｘ）
が「０」となり、「Ｄｔｈ≦ｘ」ではｆ１（ｘ）が切片
を「Ｅ」、傾き「Ｅ／Ｄｔｈ」として単調減少するもの
であり、（１１）式で示される。 ｆ１〔ｘ〕＝Ｅ−（Ｅ／Ｄｔｈ）Ｘ・・・・（１１）

【００１２】この方法は、低周波画像の画素値Ｄｔｈ以
上の濃度値を圧縮するものであり、そのアルゴリズム
は、上記の日本放射線技術学会雑誌に記載された方法と
変わらない。

【００１３】また、特許第１，５３０，８３２号公報に
は、一定濃度値以上の画像の高周波成分を強調する鮮鋭
化方法が記載されている。この方法は、超低周波数成分
を強調し、コントラストを強くすることにより、診断性
能を上げること、超低周波数成分を強調すると同時に、
雑音の占める割合が大きい高周波数成分を相対的に低減
し、視覚的に見やすい画像が得られるようにすること、
偽画像を防止し、また雑音の増大を防止して診断性能を
向上させることを目的とするものである。

【００１４】即ち、この鮮鋭化方法は、補償後（処理
後）の画素値Ｓ_D 、オリジナル画素値（入力画素値）Ｓ
_org 、オリジナル画像（入力画像）をマスクサイズＭ×
Ｍ画素で移動平均をとった時の平均画素値Ｓ_USとする。
また、ＢをＳ_org 又はＳ_USの値の増大に応じて単調増加
する変数とすると、 Ｓ_D ＝Ｓ_org ＋Ｂ×（Ｓ_org −Ｓ_US）・・・・（１２） Ｓ_US＝ΣＳ_org ／Ｍ² ・・・・（１３） なる式（１２）、（１３）で表わされるものである。

【００１５】上記（１２）式を実行した場合は、特定濃
度値以上の画像の高周波成分を強調できるという効果が
ある。

【００１６】

【発明が解決しようとする課題】しかしながら、上記Ｓ
ＰＩＥＶｏｌ．６２６ＭｅｄｉｃｉｎｅＸＩＶ／ＡＣＳ
ＩＶ（１９８６）に記載される方法では、一定濃度範囲
のダイナミックレンジを圧縮する思想がないので、画像
全体のダイナミックレンジを均等に圧縮し、そのため、
一定濃度範囲だけを圧縮することができない。従って、
例えば肺正面画像に本手法を用いた場合、縦隔部のみな
らず診断に有効な師部の濃度レンジをも圧縮してしま
い、縦隔部のみを圧縮した場合よりも診断性能が下がる
という問題がある。

【００１７】一般にダイナミックレンジ圧縮された画像
は、ＣＲＴ表示、フィルム出力する場合に、再度階調変
換する。上記特許第２，５０９，５０３号公報の方法及
び上記日本放射線技術学会雑誌等の方法では、階調変換
後の画像の高周波成分の振幅を調整する思想が無いた
め、高周波成分の振幅が階調変換曲線の傾きに依存して
変換される問題がある。そのため、診断に有効な濃度領
域の高周波成分の振幅が小さくなり、有益情報が消失す
るという問題がある。つまり、一度階調変換された画像
の高周波成分を変更することができない。さらに、エッ
ジ部分でオーバーシュートやアンダーシュートが起きる
という問題がある。

【００１８】一方、特許第１，５３０，８３２号公報の
鮮鋭化方法では、高周波成分の足し込みの強さを自由に
調整できるが、ダイナミックレンジを圧縮する思想がな
く、濃度分布の広い画像を一枚のフィルムなどで表示で
きないという問題がある。

【００１９】本発明は、上記の問題を解決するために成
されたもので、画像の高周波成分の振幅を保持したま
ま、画像の濃度分布の幅を自在に調節できるようにする
ことを目的とする。

【００２０】

【課題を解決するための手段】上記の目的を達成するた
めに、本発明による画像処理装置においては、画像の階
調変換に用いる階調変換曲線の微係数に基づいて係数を
計算する微係数計算手段と、原画像の平滑化画像を作成
する平滑化画像作成手段と、上記平滑化画像作成手段で
作成された平滑化画像と上記原画像との差分を高周波成
分として計算する高周波成分作成手段と、上記高周波成
分作成手段で計算された高周波成分を、上記微係数計算
手段で計算した係数に基づいて上記原画像又は平滑化画
像に足し込む高周波成分足し込み手段と、上記高周波成
分足し込み手段で高周波成分を足し込まれた画像を階調
変換する階調変換手段とを設けている。

【００２１】また、本発明による画像処理方法において
は、画像の階調変換に用いる階調変換曲線の微係数に基
づいて係数を計算する手順と、原画像の平滑化画像を作
成する手順と、上記平滑化画像作成手段で作成された平
滑化画像と上記原画像との差分を高周波成分として計算
する手順と、上記高周波成分作成手段で計算された高周
波成分を、上記微係数計算手段で計算した係数に基づい
て上記原画像又は平滑化画像に足し込む手順と、上記高
周波成分足し込み手段で高周波成分を足し込まれた画像
を階調変換する手順とを設けている。

【００２２】また、本発明による記憶媒体においては、
画像の階調変換に用いる階調変換曲線の微係数に基づい
て係数を計算する処理と、原画像の平滑化画像を作成す
る処理と、上記平滑化画像作成手段で作成された平滑化
画像と上記原画像との差分を高周波成分として計算する
処理と、上記高周波成分作成手段で計算された高周波成
分を、上記微係数計算手段で計算した係数に基づいて上
記原画像又は平滑化画像に足し込む処理と、上記高周波
成分足し込み手段で高周波成分を足し込まれた画像を階
調変換する処理とを実行するためのプログラムを記憶し
ている。

【００２３】

【発明の実施の形態】以下、本発明の実施の形態を図面
と共に説明する。図１は、本発明の実施の形態による画
像処理装置を含むＸ線デジタル撮影装置１００を示す。
図１において、このＸ線デジタル撮影装置１００は、ダ
イナミックレンジ圧縮（ＤＲＣ）機能を有するＸ線画像
の画像処理装置であり、前処理回路１０６、ＣＰＵ１０
８、メインメモリ１０９、操作パネル１１０、画像表示
器１１１、ＤＲＣ回路１１２を備えており、ＣＰＵバス
１０７を介して互いにデータ授受されるようになされて
いる。

【００２４】また、前処理回路１０６に接続されたデー
タ収集回路１０５と、データ収集回路１０５に接続され
た２次元Ｘ線センサ１０４及びＸ線発生回路１０１とを
備えており、これらの各回路はＣＰＵバス１０７にも接
続されている。

【００２５】図２は本実施の形態の処理の流れを示すフ
ローチャートである。図３、図４はＤＲＣ回路１１２内
の階調変換回路１１２ｅで用いる階調変換曲線を示し、
横軸が入力画像の画素値、縦軸が出力画像の画素値を示
す。図３はＳ字型の階調変換曲線、図４は入力画素値Ｂ
以下の傾きがＡ／Ｂ、入力画素値Ｂより上の傾きが１と
なっている。

【００２６】図１において、メインメモリ１０９は、Ｃ
ＰＵ１０８の処理に必要な各種のデータ等が記憶される
と共に、ＣＰＵ１０８の作業用としてのワークメモリを
含む。ＣＰＵ１０８は、メインメモリ１０９を用いて、
操作パネル１１０からの操作に従って装置全体の動作制
御等を行う。これによりＸ線デジタル撮影装置１００
は、以下のように動作する。

【００２７】まず、Ｘ線発生回路１０１は、被検査体１
０２に対してＸ線ビーム１０２を放射する。放射された
Ｘ線ビーム１０２は、被検査体１０３を減衰しながら透
過して、２次元Ｘ線センサ１０４に到達し、２次元Ｘ線
センサ１０４によりＸ線画像として出力される。ここで
は、２次元Ｘ線センサ１０４から出力されるＸ線画像
を、例えば人体胸部画像等とする。

【００２８】データ収集回路１０５は、２次元Ｘ線セン
サ１０４から出力されたＸ線画像を電気信号に変換して
前処理回路１０６に供給する。前処理回路１０６は、デ
ータ収集回路１０５からの信号（Ｘ線画像信号）に対し
て、オフセット補正処理やゲイン補正処理等の前処理を
行う。前処理されたＸ線画像信号は入力画像として、Ｃ
ＰＵ１０８の制御により、ＣＰＵバス１０７を介してメ
インメモリ１０９、ＤＲＣ回路１１２に転送される。

【００２９】また、ＤＲＣ回路１１２図において、１１
２ａは、階調変換回路１１２ｅで行う階調変換の階調変
換曲線の微係数を計算する微係数計算回路であり、メイ
ンメモリ１０９は、上記計算された微係数をＣＰＵバス
１０７を介して記憶する。１１２ｂは、前処理回路１０
６から出力された原画像１１３の平滑化画像（低周波画
像）１１４を作成する平滑化画像作成回路、１１２ｃ
は、平滑化画像１１４と原画像１１３との差分を計算す
る高周波成分作成回路、１１２ｄは高周波成分作成回路
１１２ｃで作成された高周波成分を、メインメモリ１０
９に記憶された微係数の値に応じて原画像１１３に足し
込む高周波成分足し込み回路、１１２ｅは、上記高周波
成分を足し込まれた画像を、微係数計算回路１１２ａで
微係数を計算するための元になった階調曲線で階調変換
する階調変換回路である。

【００３０】次にＤＲＣ回路１１２の動作について図２
の処理の流れに従い説明する。ＣＰＵバス１０７を介し
て前処理回路１０６で処理された原画像１１３、ｆ０
（ｘ，ｙ）を、ＣＰＵ１０８の制御を介して受信したＤ
ＲＣ回路１１２は、微係数計算回路１１２ａで階調変換
回路１１２ｅで用いる階調変換曲線の微係数から係数ｃ
（ｘ）を（１４）式に従い計算する。ここで係数ｃ
（ｘ）は、階調変換曲線の傾きを１から引いた値であ
り、Ｆ（ｘ）はｘが入力画素値を示す階調変換曲線であ
る。 ｃ（ｘ）＝１／｛∂Ｆ（ｘ）／∂ｘ｝−１・・・・（１４）

【００３１】そして、この微係数ｃ（ｘ）をＣＰＵバス
１０７を介してメインメモリ１０９に記憶する（ステッ
プＳ２０１）。

【００３２】次に、平滑化画像作成回路１１２ｂは、
（１５）式に従い平滑化画像ｆｕｓ（ｘ，ｙ）（１１
４）を計算する（ステップＳ２０２）。ここでｘ，ｙは
画像上の座標を示し、ｄはマスクサイズの大きさを示す
定数である。

【００３３】

【数１】

【００３４】次に、高周波成分作成回路１１２ｂは、
（１６）式で示すように原画像ｆ０（ｘ，ｙ）と上記平
滑化画像ｆｕｓ（ｘ，ｙ）から高周波画像ｆｈ（ｘ，
ｙ）を計算する（ステップＳ２０３）。 ｆｈ（ｘ，ｙ）＝ｆ０（ｘ，ｙ）−ｆｕｓ（ｘ，ｙ）・・・・（１６）

【００３５】次に、高周波成分足し込み回路１１２ｄ
は、（１７）式で示すように原画像ｆ０（ｘ，ｙ）に上
記計算した高周波成分ｆｈ（ｘ，ｙ）を、メインメモリ
１０９に記憶した微係数ｃ（ｘ）に基づいて足し込み、
足し込み処理後画像ｆ１（ｘ，ｙ）を得る（ステップＳ
２０４）。 ｆ１（ｘ，ｙ）＝ｆ０（ｘ，ｙ）＋ａ×ｃ（ｆ０（ｘ，ｙ））×ｆｈ（ｘ，ｙ ）・・・・（１７） ここで、ａは定数で、例えば１である。

【００３６】また、（１８）式で示すように、ｆｕｓ
（ｘ，ｙ）を用いても、ほぼ同様の効果が得られること
が実験で確認されている。 ｆ１（ｘ，ｙ）＝ｆ０（ｘ，ｙ）＋ａ×ｃ（ｆｕｓ（ｘ，ｙ））×ｆｈ（ｘ， ｙ）・・・・（１８）

【００３７】さらに、（１９）、（２０）式で示すよう
にしても、同様の効果が得られることが実験で確認され
ている。 ｃ１（ｘ）＝１／｛∂Ｆ（ｘ）／∂ｘ｝・・・・（１９） ｆ１（ｘ，ｙ）＝ｆｕｓ（ｘ，ｙ）＋ａ×ｃ１（ｆ０（ｘ，ｙ））×ｆｈ（ｘ ，ｙ）・・・・（２０）

【００３８】そして、階調変換回路１１２ｅは、上記足
し込み処理後画像ｆ１（ｘ，ｙ）を（２１）式に従って
階調変換を行い、最終処理済み画像ｆｄ（ｘ，ｙ）を得
る。 ｆｄ（ｘ，ｙ）＝Ｆ（ｆ１（ｘ，ｙ））・・・・（２
１）

【００３９】また、平滑化画像ｆｕｓ（ｘ，ｙ）を（２
２）〜（２５）式で示すように、モルフォジ演算を用い
て計算してもよい。

【００４０】 ｆ２（ｘ，ｙ）＝ｍｉｎ｛ｆ０（ｘ＋ｘ１，ｙ＋ｙ１）−Ｄ（ｘ１，ｙ１）｜ ｘ１×ｘ１＋ｙ１×ｙ１≦ｒ１×ｒ１｝・・・・（２２） ｆ３（ｘ，ｙ）＝ｍａｘ｛ｆ２（ｘ＋ｘ１，ｙ＋ｙ１）＋Ｄ（ｘ１，ｙ１）｜ ｘ１×ｘ１＋ｙ１×ｙ１≦ｒ１×ｒ１｝・・・・（２３） ｆ４（ｘ，ｙ）＝ｍａｘ｛ｆ３（ｘ＋ｘ１，ｙ＋ｙ１）＋Ｄ（ｘ１，ｙ１）｜ ｘ１×ｘ１＋ｙ１×ｙ１≦ｒ１×ｒ１｝・・・・（２４） ｆｕｓ（ｘ，ｙ）＝ｍｉｎ｛ｆ４（ｘ＋ｘ１，ｙ＋ｙ１）−Ｄ（ｘ１，ｙ１） ｜ｘ１×ｘ１＋ｙ１×ｙ１≦ｒ１×ｒ１｝・・・（２５）

【００４１】ここで、Ｄ（ｘ，ｙ）を円盤状フィルタ、
ｒ１を任意の定数とし、入力画像に応じて選択される。 Ｄ（ｘ，ｙ）＝０、 ｘ×ｘ＋ｙ×ｙ≦ｒ１×ｒ１・・・・（２６） ＝−∞、 その他（２６）

【００４２】モルフォジ演算は、モルフォロジカルフィ
ルタを用いて行われ、ここで得られたｆｕｓ（ｘ，ｙ）
のプロファイルは、エッジ構造を保存しているものであ
り、従来のダイナミックレンジ圧縮の欠点であるオーバ
ーシュートやアンダーシュートが起きないものである。

【００４３】本実施の形態によれば、入力画像の任意の
階調領域の濃度分布幅を圧縮、伸張することができ、か
つ階調変換後の高周波成分の振幅を、階調変換前の画像
の高周波成分の振幅に保持できる効果がある。さらに、
オーバーシュートやアンダーシュートが起きない効果が
ある。また、平滑化画像に濃度平均を用いた場合は、計
算時間を短縮できる効果がある。

【００４４】次に、本発明の他の実施の形態としての記
憶媒体について説明する。本発明はハードウェア構成に
より実現することもできるが、ＣＰＵやメモリからなる
コンピュータシステムに構成することもできる。上記実
施の形態による図１の各機能ブロックによるシステム
を、ＣＰＵ１０８やメインメモリ１０９等からなるコン
ピュータシステムに構成する場合、上記メモリは本発明
による記憶媒体を構成する。この記憶媒体には、図２の
フローチャートを含む前述した動作を制御するための処
理手順を実行するためのプログラムが記憶される。

【００４５】また、この記憶媒体としては、半導体メモ
リ、光ディスク、光磁気ディスク、磁気媒体等を用いて
よく、これらをＲＯＭ、ＲＡＭ、ＣＤ−ＲＯＭ、フロッ
ピィディスク、磁気テープ、磁気カード、不揮発性メモ
リカード等に構成して用いてよい。

【００４６】従って、この記憶媒体を図１に示したシス
テム以外の他のシステムあるいは装置で用い、そのシス
テムあるいはコンピュータがこの記憶媒体に格納された
プログラムコードを読み出し、実行することによって
も、上記実施の形態と同等の機能を実現できると共に、
同等の効果を得ることができ、本発明の目的を達成する
ことができる。

【００４７】また、コンピュータ上で稼働しているＯＳ
等が処理の一部又は全部を行う場合、あるいは記憶媒体
から読み出されたプログラムコードが、コンピュータに
挿入された拡張機能ボードやコンピュータに接続された
拡張機能ユニットに備わるメモリに書き込まれた後、そ
のプログラムコードの指示に基づいて、上記拡張機能ボ
ードや拡張機能ユニットに備わるＣＰＵ等が処理の一部
又は全部を行う場合にも、上記実施の形態と同等の機能
を実現できると共に、同等の効果を得ることができ、本
発明の目的を達成することができる。

【００４８】

【発明の効果】以上説明したように、請求項１、４、７
記載の発明によれば、入力画像の任意の階調領域の濃度
分布幅を圧縮、伸張することができ、かつ階調変換後の
高周波成分の振幅を、階調変換前の画像の高周波成分の
振幅に保持できる効果がある。

【００４９】請求項２、５、８記載の発明によれば、平
滑化画像をモルフォロジカルフィルタで作成するため、
元画像のエッジ構造を保存することができ、アンダーシ
ュートやオーバーシュートが起きない効果がある。

【００５０】請求項３、６、９記載の発明によれば、平
滑化画像を平均濃度値を用いて作成するため、計算時間
を短縮できる効果がある。

【図面の簡単な説明】

【図１】本発明の実施の形態による画像処理装置を用い
たＸ線デジタル撮影装置示す構成図である。

【図２】本発明の実施の形態の処理の流れを示すフロー
チャートである。

【図３】本発明実施の形態による階調変換関数Ｆ１（）
を示す図である。

【図４】本発明実施の形態による階調変換関数Ｆ１（）
を示す図である。

【符号の説明】

１０８ ＣＰＵ １０９ メインメモリ １１２ ＤＲＣ回路 １１２ａ 微係数計算回路 １１２ｂ 平滑化画像作成回路 １１２ｃ 高周波成分作成回路 １１２ｄ 高周波成分足し込み回路 １１２ｅ 階調変換回路 １１３ 原画像 １１４ 平滑化画像

───────────────────────────────────────────────────── フロントページの続き Ｆターム(参考） 4C093 AA16 AA26 CA13 CA36 DA03 EB17 FC18 FC19 FD01 FD05 FD08 FD11 FD12 FF06 FF08 FF09 FF34 5B057 AA08 BA03 CA02 CA08 CA12 CA16 CB02 CB08 CB12 CB16 CC01 CE03 CE05 CE06 CE11 CH01 CH09 5C021 PA32 PA34 PA66 PA75 PA76 RB00 RB04 XA35 XB03 ZA00

Claims (9)

【特許請求の範囲】
1. 【請求項１】 画像の階調変換に用いる階調変換曲線の
   微係数に基づいて係数を計算する微係数計算手段と、 原画像の平滑化画像を作成する平滑化画像作成手段と、 上記平滑化画像作成手段で作成された平滑化画像と上記
   原画像との差分を高周波成分として計算する高周波成分
   作成手段と、 上記高周波成分作成手段で計算された高周波成分を、上
   記微係数計算手段で計算した係数に基づいて上記原画像
   又は平滑化画像に足し込む高周波成分足し込み手段と、 上記高周波成分足し込み手段で高周波成分を足し込まれ
   た画像を階調変換する階調変換手段とを備えた画像処理
   装置。
2. 【請求項２】 上記平滑化画像作成手段は、モルフォロ
   ジカルフィルタを用いて上記平滑化画像を作成すること
   を特徴とする請求項１記載の画像処理装置。
3. 【請求項３】 上記平滑化画像作成手段は、上記原画像
   の濃度平均値を用いて上記平滑化画像を作成することを
   特徴とする請求項１記載の画像処理装置。
4. 【請求項４】 画像の階調変換に用いる階調変換曲線の
   微係数に基づいて係数を計算する手順と、 原画像の平滑化画像を作成する手順と、 上記作成された平滑化画像と上記原画像との差分を高周
   波成分として計算する手順と、 上記計算された高周波成分を、上記計算した係数に基づ
   いて上記原画像又は平滑化画像に足し込む手順と、 上記高周波成分を足し込まれた画像を階調変換する手順
   とを備えた画像処理方法。
5. 【請求項５】 上記平滑化画像を作成する手順は、モル
   フォロジカルフィルタを用いて上記平滑化画像を作成す
   ることを特徴とする請求項４記載の画像処理方法。
6. 【請求項６】 上記平滑化画像を作成する手順は、上記
   原画像の濃度平均値を用いて上記平滑化画像を作成する
   ことを特徴とする請求項４記載の画像処理方法。
7. 【請求項７】 画像の階調変換に用いる階調変換曲線の
   微係数に基づいて係数を計算する処理と、 原画像の平滑化画像を作成する処理と、 上記作成された平滑化画像と上記原画像との差分を高周
   波成分として計算する処理と、 上記計算された高周波成分を、上記計算された係数に基
   づいて上記原画像又は平滑化画像に足し込む処理と、 上記高周波成分を足し込まれた画像を階調変換する処理
   とを実行するためのプログラムを記憶したコンピュータ
   読み取り可能な記憶媒体。
8. 【請求項８】 上記平滑化画像を作成する処理は、モル
   フォロジカルフィルタを用いて上記平滑化画像を作成す
   ることを特徴とする請求項７記載のコンピュータ読み取
   り可能な記憶媒体。
9. 【請求項９】 上記平滑化画像を作成する処理は、上記
   原画像の濃度平均値を用いて上記平滑化画像を作成する
   ことを特徴とする請求項７記載のコンピュータ読み取り
   可能な記憶媒体。