JP3937616B2

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

Info

Publication number: JP3937616B2
Application number: JP33987398A
Authority: JP
Inventors: 弘之新畠
Original assignee: Canon Inc
Current assignee: Canon Inc
Priority date: 1998-11-30
Filing date: 1998-11-30
Publication date: 2007-06-27
Anticipated expiration: 2018-11-30
Also published as: JP2000163570A

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／Ｍ²・・・・（１３）
なる式（１２）、（１３）で表わされるものである。
【００１５】
上記（１２）式を実行した場合は、特定濃度値以上の画像の高周波成分を強調できるという効果がある。
【００１６】
【発明が解決しようとする課題】
しかしながら、上記ＳＰＩＥＶｏｌ．６２６ＭｅｄｉｃｉｎｅＸＩＶ／ＡＣＳＩＶ（１９８６）に記載される方法では、一定濃度範囲のダイナミックレンジを圧縮する思想がないので、画像全体のダイナミックレンジを均等に圧縮し、そのため、一定濃度範囲だけを圧縮することができない。従って、例えば肺正面画像に本手法を用いた場合、縦隔部のみならず診断に有効な師部の濃度レンジをも圧縮してしまい、縦隔部のみを圧縮した場合よりも診断性能が下がるという問題がある。
【００１７】
一般にダイナミックレンジ圧縮された画像は、ＣＲＴ表示、フィルム出力する場合に、再度階調変換する。上記特許第２，５０９，５０３号公報の方法及び上記日本放射線技術学会雑誌等の方法では、階調変換後の画像の高周波成分の振幅を調整する思想が無いため、高周波成分の振幅が階調変換曲線の傾きに依存して変換される問題がある。そのため、診断に有効な濃度領域の高周波成分の振幅が小さくなり、有益情報が消失するという問題がある。つまり、一度階調変換された画像の高周波成分を変更することができない。さらに、エッジ部分でオーバーシュートやアンダーシュートが起きるという問題がある。
【００１８】
一方、特許第１，５３０，８３２号公報の鮮鋭化方法では、高周波成分の足し込みの強さを自由に調整できるが、ダイナミックレンジを圧縮する思想がなく、濃度分布の広い画像を一枚のフィルムなどで表示できないという問題がある。
【００１９】
本発明は、上記の問題を解決するために成されたもので、画像の高周波成分の振幅を保持したまま、画像の濃度分布の幅を自在に調節できるようにすることを目的とする。
【００２０】
【課題を解決するための手段】
本発明の画像処理方法は、原画像データに対し階調変換処理及び周波数処理を施す画像処理方法であって、上記原画像データから高周波成分を生成する工程と、上記階調変換処理に基づいて上記高周波成分の値を変更する工程と、上記変更された高周波成分に基づいて第二の画像データを生成する工程と、上記第二の画像データの画素値を階調変換処理で変更する工程とを備え、上記高周波成分の値を、階調変換曲線の微分値に基づいて変更することを特徴とする。
【００２１】
また、本発明の画像処理装置は、原画像データに対し階調変換処理及び周波数処理を施す画像処理装置であって、上記原画像データから高周波成分を生成する手段と、上記階調変換処理に基づいて上記高周波成分の値を変更する手段と、上記変更された高周波成分に基づいて第二の画像データを生成する手段と、上記第二の画像データの画素値を階調変換処理で変更する手段とを備え、上記高周波成分の値を、階調変換曲線の微分値に基づいて変更することを特徴とする。
【００２２】
また、本発明のコンピュータ読み取り可能な記憶媒体は、原画像データに対し階調変換処理及び周波数処理を施す画像処理方法をコンピュータに実行させるためのプログラムを記憶したコンピュータ読み取り可能な記憶媒体であって、上記原画像データから高周波成分を生成する処理と、上記高周波成分の値を、階調変換曲線の微分値に基づいて変更する処理と、上記変更された高周波成分に基づいて第二の画像データを生成する処理と、上記第二の画像データの画素値を階調変換処理で変更する処理とをコンピュータに実行させるためのプログラムを記憶したことを特徴とする。
【００２３】
【発明の実施の形態】
以下、本発明の実施の形態を図面と共に説明する。
図１は、本発明の実施の形態による画像処理装置を含むＸ線デジタル撮影装置１００を示す。
図１において、このＸ線デジタル撮影装置１００は、ダイナミックレンジ圧縮（ＤＲＣ）機能を有するＸ線画像の画像処理装置であり、前処理回路１０６、ＣＰＵ１０８、メインメモリ１０９、操作パネル１１０、画像表示器１１１、ＤＲＣ回路１１２を備えており、ＣＰＵバス１０７を介して互いにデータ授受されるようになされている。
【００２４】
また、前処理回路１０６に接続されたデータ収集回路１０５と、データ収集回路１０５に接続された２次元Ｘ線センサ１０４及びＸ線発生回路１０１とを備えており、これらの各回路はＣＰＵバス１０７にも接続されている。
【００２５】
図２は本実施の形態の処理の流れを示すフローチャートである。
図３、図４はＤＲＣ回路１１２内の階調変換回路１１２ｅで用いる階調変換曲線を示し、横軸が入力画像の画素値、縦軸が出力画像の画素値を示す。図３はＳ字型の階調変換曲線、図４は入力画素値Ｂ以下の傾きがＡ／Ｂ、入力画素値Ｂより上の傾きが１となっている。
【００２６】
図１において、メインメモリ１０９は、ＣＰＵ１０８の処理に必要な各種のデータ等が記憶されると共に、ＣＰＵ１０８の作業用としてのワークメモリを含む。ＣＰＵ１０８は、メインメモリ１０９を用いて、操作パネル１１０からの操作に従って装置全体の動作制御等を行う。これによりＸ線デジタル撮影装置１００は、以下のように動作する。
【００２７】
まず、Ｘ線発生回路１０１は、被検査体１０２に対してＸ線ビーム１０２を放射する。放射されたＸ線ビーム１０２は、被検査体１０３を減衰しながら透過して、２次元Ｘ線センサ１０４に到達し、２次元Ｘ線センサ１０４によりＸ線画像として出力される。ここでは、２次元Ｘ線センサ１０４から出力されるＸ線画像を、例えば人体胸部画像等とする。
【００２８】
データ収集回路１０５は、２次元Ｘ線センサ１０４から出力されたＸ線画像を電気信号に変換して前処理回路１０６に供給する。前処理回路１０６は、データ収集回路１０５からの信号（Ｘ線画像信号）に対して、オフセット補正処理やゲイン補正処理等の前処理を行う。前処理されたＸ線画像信号は入力画像として、ＣＰＵ１０８の制御により、ＣＰＵバス１０７を介してメインメモリ１０９、ＤＲＣ回路１１２に転送される。
【００２９】
また、ＤＲＣ回路１１２図において、１１２ａは、階調変換回路１１２ｅで行う階調変換の階調変換曲線の微係数を計算する微係数計算回路であり、メインメモリ１０９は、上記計算された微係数をＣＰＵバス１０７を介して記憶する。１１２ｂは、前処理回路１０６から出力された原画像１１３の平滑化画像（低周波画像）１１４を作成する平滑化画像作成回路、１１２ｃは、平滑化画像１１４と原画像１１３との差分を計算する高周波成分作成回路、１１２ｄは高周波成分作成回路１１２ｃで作成された高周波成分を、メインメモリ１０９に記憶された微係数の値に応じて原画像１１３に足し込む高周波成分足し込み回路、１１２ｅは、上記高周波成分を足し込まれた画像を、微係数計算回路１１２ａで微係数を計算するための元になった階調曲線で階調変換する階調変換回路である。
【００３０】
次にＤＲＣ回路１１２の動作について図２の処理の流れに従い説明する。
ＣＰＵバス１０７を介して前処理回路１０６で処理された原画像１１３、ｆ０（ｘ，ｙ）を、ＣＰＵ１０８の制御を介して受信したＤＲＣ回路１１２は、微係数計算回路１１２ａで階調変換回路１１２ｅで用いる階調変換曲線の微係数から係数ｃ（ｘ）を（１４）式に従い計算する。ここで係数ｃ（ｘ）は、階調変換曲線の傾きを１から引いた値であり、Ｆ（ｘ）はｘが入力画素値を示す階調変換曲線である。
ｃ（ｘ）＝１／｛∂Ｆ（ｘ）／∂ｘ｝−１・・・・（１４）
【００３１】
そして、この微係数ｃ（ｘ）をＣＰＵバス１０７を介してメインメモリ１０９に記憶する（ステップＳ２０１）。
【００３２】
次に、平滑化画像作成回路１１２ｂは、（１５）式に従い平滑化画像ｆｕｓ（ｘ，ｙ）（１１４）を計算する（ステップＳ２０２）。ここでｘ，ｙは画像上の座標を示し、ｄはマスクサイズの大きさを示す定数である。
【００３３】
【数１】

【００３４】
次に、高周波成分作成回路１１２ｂは、（１６）式で示すように原画像ｆ０（ｘ，ｙ）と上記平滑化画像ｆｕｓ（ｘ，ｙ）から高周波画像ｆｈ（ｘ，ｙ）を計算する（ステップＳ２０３）。
ｆｈ（ｘ，ｙ）＝ｆ０（ｘ，ｙ）−ｆｕｓ（ｘ，ｙ）・・・・（１６）
【００３５】
次に、高周波成分足し込み回路１１２ｄは、（１７）式で示すように原画像ｆ０（ｘ，ｙ）に上記計算した高周波成分ｆｈ（ｘ，ｙ）を、メインメモリ１０９に記憶した微係数ｃ（ｘ）に基づいて足し込み、足し込み処理後画像ｆ１（ｘ，ｙ）を得る（ステップＳ２０４）。
ｆ１（ｘ，ｙ）＝ｆ０（ｘ，ｙ）＋ａ×ｃ（ｆ０（ｘ，ｙ））×ｆｈ（ｘ，ｙ）・・・・（１７）
ここで、ａは定数で、例えば１である。
【００３６】
また、（１８）式で示すように、ｆｕｓ（ｘ，ｙ）を用いても、ほぼ同様の効果が得られることが実験で確認されている。
ｆ１（ｘ，ｙ）＝ｆ０（ｘ，ｙ）＋ａ×ｃ（ｆｕｓ（ｘ，ｙ））×ｆｈ（ｘ，ｙ）・・・・（１８）
【００３７】
さらに、（１９）、（２０）式で示すようにしても、同様の効果が得られることが実験で確認されている。
ｃ１（ｘ）＝１／｛∂Ｆ（ｘ）／∂ｘ｝・・・・（１９）
ｆ１（ｘ，ｙ）＝ｆｕｓ（ｘ，ｙ）＋ａ×ｃ１（ｆ０（ｘ，ｙ））×ｆｈ（ｘ，ｙ）・・・・（２０）
【００３８】
そして、階調変換回路１１２ｅは、上記足し込み処理後画像ｆ１（ｘ，ｙ）を（２１）式に従って階調変換を行い、最終処理済み画像ｆｄ（ｘ，ｙ）を得る。
ｆｄ（ｘ，ｙ）＝Ｆ（ｆ１（ｘ，ｙ））・・・・（２１）
【００３９】
また、平滑化画像ｆｕｓ（ｘ，ｙ）を（２２）〜（２５）式で示すように、モルフォロジー演算を用いて計算しても良い。
【００４０】
ｆ２（ｘ，ｙ）＝ｍｉｎ｛ｆ０（ｘ＋ｘ１，ｙ＋ｙ１）−Ｄ（ｘ１，ｙ１）｜ｘ１×ｘ１＋ｙ１×ｙ１≦ｒ１×ｒ１｝・・・・（２２）
ｆ３（ｘ，ｙ）＝ｍａｘ｛ｆ２（ｘ＋ｘ１，ｙ＋ｙ１）＋Ｄ（ｘ１，ｙ１）｜ｘ１×ｘ１＋ｙ１×ｙ１≦ｒ１×ｒ１｝・・・・（２３）
ｆ４（ｘ，ｙ）＝ｍａｘ｛ｆ３（ｘ＋ｘ１，ｙ＋ｙ１）＋Ｄ（ｘ１，ｙ１）｜ｘ１×ｘ１＋ｙ１×ｙ１≦ｒ１×ｒ１｝・・・・（２４）
ｆｕｓ（ｘ，ｙ）＝ｍｉｎ｛ｆ４（ｘ＋ｘ１，ｙ＋ｙ１）−Ｄ（ｘ１，ｙ１）｜ｘ１×ｘ１＋ｙ１×ｙ１≦ｒ１×ｒ１｝・・・（２５）
【００４１】
ここで、Ｄ（ｘ，ｙ）を円盤状フィルタ、ｒ１を任意の定数とし、入力画像に応じて選択される。
Ｄ（ｘ，ｙ）＝０、ｘ×ｘ＋ｙ×ｙ≦ｒ１×ｒ１・・・・（２６）
＝−∞、その他（２６）
【００４２】
モルフォロジー演算は、モルフォロジカルフィルタを用いて行われ、ここで得られたｆｕｓ（ｘ，ｙ）のプロファイルは、エッジ構造を保存しているものであり、従来のダイナミックレンジ圧縮の欠点であるオーバーシュートやアンダーシュートが起きないものである。
【００４３】
本実施の形態によれば、入力画像の任意の階調領域の濃度分布幅を圧縮、伸張することができ、かつ階調変換後の高周波成分の振幅を、階調変換前の画像の高周波成分の振幅に保持できる効果がある。さらに、オーバーシュートやアンダーシュートが起きない効果がある。また、平滑化画像に濃度平均を用いた場合は、計算時間を短縮できる効果がある。
【００４４】
次に、本発明の他の実施の形態としての記憶媒体について説明する。
本発明はハードウェア構成により実現することもできるが、ＣＰＵやメモリからなるコンピュータシステムに構成することもできる。上記実施の形態による図１の各機能ブロックによるシステムを、ＣＰＵ１０８やメインメモリ１０９等からなるコンピュータシステムに構成する場合、上記メモリは本発明による記憶媒体を構成する。この記憶媒体には、図２のフローチャートを含む前述した動作を制御するための処理手順を実行するためのプログラムが記憶される。
【００４５】
また、この記憶媒体としては、半導体メモリ、光ディスク、光磁気ディスク、磁気媒体等を用いてよく、これらをＲＯＭ、ＲＡＭ、ＣＤ−ＲＯＭ、フロッピィディスク、磁気テープ、磁気カード、不揮発性メモリカード等に構成して用いてよい。
【００４６】
従って、この記憶媒体を図１に示したシステム以外の他のシステムあるいは装置で用い、そのシステムあるいはコンピュータがこの記憶媒体に格納されたプログラムコードを読み出し、実行することによっても、上記実施の形態と同等の機能を実現できると共に、同等の効果を得ることができ、本発明の目的を達成することができる。
【００４７】
また、コンピュータ上で稼働しているＯＳ等が処理の一部又は全部を行う場合、あるいは記憶媒体から読み出されたプログラムコードが、コンピュータに挿入された拡張機能ボードやコンピュータに接続された拡張機能ユニットに備わるメモリに書き込まれた後、そのプログラムコードの指示に基づいて、上記拡張機能ボードや拡張機能ユニットに備わるＣＰＵ等が処理の一部又は全部を行う場合にも、上記実施の形態と同等の機能を実現できると共に、同等の効果を得ることができ、本発明の目的を達成することができる。
【００４８】
【発明の効果】
以上説明したように本発明によれば、入力画像の任意の階調領域の濃度分布幅を圧縮、伸張することができ、かつ階調変換後の高周波成分の振幅を、階調変換前の画像の高周波成分の振幅に保持できる効果がある。
【００４９】
また、平滑化画像をモルフォロジカルフィルタを用いて作成するようにした場合には、元画像のエッジ構造を保存することができ、アンダーシュートやオーバーシュートが起きない効果がある。
【図面の簡単な説明】
【図１】本発明の実施の形態による画像処理装置を用いたＸ線デジタル撮影装置示す構成図である。
【図２】本発明の実施の形態の処理の流れを示すフローチャートである。
【図３】本発明実施の形態による階調変換関数Ｆ１（）を示す図である。
【図４】本発明実施の形態による階調変換関数Ｆ１（）を示す図である。
【符号の説明】
１０８ＣＰＵ
１０９メインメモリ
１１２ＤＲＣ回路
１１２ａ微係数計算回路
１１２ｂ平滑化画像作成回路
１１２ｃ高周波成分作成回路
１１２ｄ高周波成分足し込み回路
１１２ｅ階調変換回路
１１３原画像
１１４平滑化画像

Claims

原画像データに対し階調変換処理及び周波数処理を施す画像処理方法であって、
上記原画像データから高周波成分を生成する工程と、
上記階調変換処理に基づいて上記高周波成分の値を変更する工程と、
上記変更された高周波成分に基づいて第二の画像データを生成する工程と、
上記第二の画像データの画素値を階調変換処理で変更する工程とを備え、
上記高周波成分の値を、階調変換曲線の微分値に基づいて変更することを特徴とする画像処理方法。
上記階調変換処理は、上記階調変換曲線に基づいて上記第二の画像データの画素値を変更することを特徴とする請求項１に記載の画像処理方法。
上記高周波成分は、上記原画像データの平滑化画像を生成し、該平滑化画像を原画像データから減算することで生成されることを特徴とする請求項１に記載の画像処理方法。
上記平滑化画像はモルフォロジカルフィルタを用いて生成されることを特徴とする請求項３に記載の画像処理方法。
原画像データに対し階調変換処理及び周波数処理を施す画像処理装置であって、
上記原画像データから高周波成分を生成する手段と、
上記階調変換処理に基づいて上記高周波成分の値を変更する手段と、
上記変更された高周波成分に基づいて第二の画像データを生成する手段と、
上記第二の画像データの画素値を階調変換処理で変更する手段とを備え、
上記高周波成分の値を、階調変換曲線の微分値に基づいて変更することを特徴とする画像処理装置。
Ｘ線発生回路と、
該Ｘ線発生回路が発するＸ線を上記原画像データに変換するための２次元Ｘ線センサとを備えることを特徴とする請求項５に記載の画像処理装置。
原画像データに対し階調変換処理及び周波数処理を施す画像処理方法をコンピュータに実行させるためのプログラムを記憶したコンピュータ読み取り可能な記憶媒体であって、
上記原画像データから高周波成分を生成する処理と、
上記高周波成分の値を、階調変換曲線の微分値に基づいて変更する処理と、
上記変更された高周波成分に基づいて第二の画像データを生成する処理と、
上記第二の画像データの画素値を階調変換処理で変更する処理とをコンピュータに実行させるためのプログラムを記憶したことを特徴とするコンピュータ読み取り可能な記憶媒体。