JPH04349458A - 画像ノイズシミュレーション方法 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は画像ノイズシミュレーシ
ョン方法に関し、詳しくは、写真感光材料などの画像形
成媒体上に形成される画像のノイズ特性を精度良くシミ
ュレーションする方法に関する。

【０００２】

【従来の技術】写真感光材料などの画像形成媒体の設計
と開発に当たって、画像のノイズ（写真感光材料の粒状
性）をシミュレーションすることは、ノイズの目標設定
や評価を行う上で有効である。かかるノイズのシミュレ
ーション方法としては、ある画素（位置）におけるノイ
ズを自己回帰モデルを用いて予測発生させるようしたも
のがある（犬井正男、日本写真学会誌　　４９、５３１
　　１９８６年発行　　等参照）。ここで、ノイズの振
幅だけが画像濃度に依存する場合には、予め発生又は測
定した均一濃度部の２次元ノイズに、画像の濃度に対応
した振幅（例えばＲＭＳ粒状度）を乗じて得ることがで
き、また、任意の同周波数成分をもつノイズは、その周
波数成分に対応する２次元ウィナースペクトルを設定し
、これにランダムな位相を与えた後、逆フーリエ変換す
ることによって得ることができる。

【０００３】

【発明が解決しようとする課題】ところで、写真感光材
料としてのカラーネガフィルムでは、粒状性（ノイズ）
は振幅，周波数特性共に各色（例えばＲ，Ｇ，Ｂ）毎の
濃度に依存して変化するという特質を有するため、ノイ
ズの振幅と周波数特性との両方を共に濃度に依存させて
シミュレーションすることが望まれるが、従来では、い
ずれか一方のみを濃度に依存させてノイズをシミュレー
ションする方法のみが提案されており、かかる従来のシ
ミュレーション方法を用いたのでは、実際のカラーネガ
フィルムの画像に対する対応性が不充分であるという問
題があった。

【０００４】本発明は上記問題点に鑑みなされたもので
あり、ノイズの振幅と周波数特性との両方を濃度に依存
させてシミュレーションする構成の画像ノイズシミュレ
ーション方法を提供することにより、特にカラーネガフ
ィルムのようにノイズの振幅と周波数特性との両方が濃
度に依存して変化するような画像形成媒体において、実
際の画像に対する対応性の良いシミュレーションが行え
るようにすることを目的とする。

【０００５】

【課題を解決するための手段】そのため本発明にかかる
画像ノイズシミュレーション方法では、画像形成媒体上
に形成される画像のノイズの平均的な振幅及び周波数特
性をそれぞれ反映するパラメータと、画像形成媒体にお
ける画像濃度との関係を求めた後、前記パラメータを用
いてノイズシミュレーションの予測式の係数を決定し、
この係数を用いた前記予測式に基づいて前記画像形成媒
体上に形成された画像のノイズをシミュレーションする
よう構成した。

【０００６】ここで、前記パラメータがノイズの平均的
な振幅を反映するパラメータと周波数特性を反映するパ
ラメータとの少なくとも２つからなり、前記画像形成媒
体の平均濃度が異なる複数試料の濃度実測値から平均濃
度別に実測値による自己相関関数を設定する一方、前記
少なくとも２つのパラメータを含めて設定されたノイズ
のモデルから定まる理論的自己相関関数を算出し、実測
値による自己相関関数と理論的自己相関関数との間の誤
差の総和を最小化することで前記少なくとも２つのパラ
メータの値を前記各平均濃度別に求めることにより、ノ
イズの平均的な振幅及び周波数特性をそれぞれ反映する
前記少なくとも２つのパラメータと画像濃度との関係を
求めるよう構成することができる。

【０００７】

【作用】かかる画像ノイズシミュレーション方法による
と、画像形成媒体上に形成される画像のノイズの平均的
な振幅及び周波数特性をそれぞれ反映するパラメータを
、画像形成媒体における画像濃度に依存させることによ
り、前記パラメータを用いて決定されるノイズシミュレ
ーションの予測式の係数を、画像濃度に依存させること
ができる。従って、画像形成媒体上に形成される画像の
もつノイズにおける振幅と周波数特性とが、共に画像濃
度に依存して変化するような場合に、ノイズシミュレー
ションの予測式を、実際の画像に精度良く対応させるこ
とができるようになる。

【０００８】また、前記パラメータがノイズの平均的な
振幅を反映するパラメータと周波数特性を反映するパラ
メータとの少なくとも２つとし、濃度実測値による自己
相関関数と理論的自己相関関数とに基づいて前記少なく
とも２つのパラメータと画像濃度との関係を求めるよう
にすれば、比較的容易にパラメータの設定が行える。

【０００９】

【実施例】以下に本発明の実施例を説明する。スペクト
ル解析の一方法であるＭＥＭ（最大エントロピー法）と
数学的に等価である自己回帰モデルによって、１次元の
ノイズを容易にシミュレーションすることができ、本実
施例では、画像形成媒体上に形成された画像のノイズシ
ミュレーションにおいて、離散的一様等方確率場を表す
以下の差分方程式に基づいた２次元自己回帰モデルを用
いる。

【００１０】 　　Ｉｍ＋１，ｎ　＋Ｉｍ−１，ｎ　＋Ｉｍ，ｎ＋１　
＋Ｉｍ，ｎ−１　−（４＋Ｋ２）Ｉｍ，ｎ　＝ＣＷｍ，
ｎ　ここで、Ｉｍ，ｎ　は画素（ｍ，ｎ）におけるノイ
ズ、Ｋはノイズの相関距離に関する（主にノイズの平均
的周波数特性を反映する）パラメータ（但し、Ｋは振幅
にも影響を与える。）、Ｃはノイズの平均的振幅を反映
するパラメータである。Ｗはホワイトノイズで、このＷ
ｍ，ｎ　の平均値は０であり、〈　　〉が期待値を表す
ものとしたときに、〈Ｗｍ，ｎ　Ｗｍ’，ｎ’　〉＝δ
ｍ，ｍ’δｎ，ｎ’である規格化されたホワイトノイズ
である。

【００１１】前記差分方程式に基づいてノイズの２次元
自己回帰モデルが、以下のように導かれる。 　　Ｉｍ，ｎ　＝ｈ０，０　（Ｃ，Ｋ）Ｗｍ，ｎ　＋Σ
ｈｐ，ｑ　（Ｃ，Ｋ）Ｉｍ＋ｐ，ｎ＋ｑ　・・・（ａ）
　　　　　　　　　　　（ｐ，ｑ）∈Ｓ ここで、（ｐ，ｑ）は座標であり、Ｓは係数ｈを用いる
範囲である。

【００１２】このように、前記２つのパラメータＣ，Ｋ
に基づいて各画素に対応する係数ｈｐ，ｑ　を求め、自
己回帰モデルを用いてこれらの係数ｈｐ，ｑ　から画素
（ｍ，ｎ）におけるノイズＩｍ，ｎ　を得ることができ
る。前記２つのパラメータＣ，Ｋの値は、前記差分方程
式から導かれる下式の自己相関関数Ｒｍ，ｎ　を用いて
濃度変動の実測データから求めることができる。

【００１３】

【数１】 　　前記濃度変動の実測データに基づく前記２つのパラ
メータの決定は、以下のようにして行われる。

【００１４】尚、本実施例では、スチル写真のカラーネ
ガフィルムを画像形成媒体として用い、かかるカラーネ
ガフィルムのノイズ（粒状ノイズ）をシミュレーション
するものとする。前記カラーネガフィルムとしては超高
感度（ＩＳＯ＝３２００程度）のものを使用した。濃度
変動の実測データを得るために、まず、前記カラーネガ
フィルムを昼光で均一露光させて現像し、３原色（Ｒ，
Ｇ，Ｂ）毎の平均濃度が、赤Ｒ；１．１８〜２．９９，
緑Ｇ；０．８３〜２．３７，青Ｂ；０．６０〜１．７４
の範囲でそれぞれ変化する１０種類の濃度レベルの試料
を作成した。

【００１５】そして、前記試料の１次元微細濃度分布（
ノイズ）のデータを、マイクロデンシトメータ（微小濃
度計）と３原色のフィルタとを用いて、各色毎及び各濃
度毎にサンプリングして、濃度変動の実測データを得る
。尚、前記マイクロデンシトメータによる濃度測定にお
いては、アパーチャサイズを幅２μｍ，高さ４００　μ
ｍ程度とし、サンプリング間隔を１μｍ程度とする。

【００１６】次に、前記濃度変動の実測データに基づい
て、各色及び各濃度毎にノイズのパワースペクトルを計
算し、空間周波数領域上でマイクロデンシトメーター測
定時のアパーチャサイズの補正と電子回路のノイズ除去
とを行った後、２次元画像領域内の画素位置に関する濃
度変動の自己相関関数Ａｃｆを計算させ、ラグＬａｇ（
相関距離；隔たり画素数）に対する実験的自己相関関数
Ａｃｆの値の変化を各色毎及び各濃度毎に得る。その後
、下式で表される理論的及び実測的自己相関関数間の誤
差の２乗和Ｅを有限個のラグＬａｇの範囲で求める関数
式を定め、Ｅの値が最小化するようにパラメータＫ，Ｃ
の値を決定すれば良い。

【００１７】

【数２】 　　尚、図１は、前記ラグＬａｇと実験的な自己相関関
数Ａｃｆの値との関係を、３原色の中の緑Ｇに関して示
したものであり、ここでは、前記緑Ｇの濃度を、０．８
３，１．０３，１．２１，１．４１，１．６０，２．１
６の６種類に変化させた場合の傾向を示してある。

【００１８】本実施例においては、前記Ｋ，Ｃの両パラ
メータの濃度に対する関係を規定する上で、前述した数
２を用いる方法と実用的に大きな誤差を生じない、やや
簡便化した方法を用いるようにしており、これを以下に
詳細に説明する。まず、前記数１の理論的自己相関関数
Ｒｍ，ｎ　（＝Ｆ（Ｋ，Ｃ，ｍ，ｎ））について、カラ
ーネガフィルムでは自己相関関数Ａｃｆが等方的である
ので、ｎを０に固定し、更に、ｍ＝０で正規化（ノーマ
ライズ）することで、ｎ，Ｃの因子を除いた自己相関関
数Ａｃｆ＝Ｒｍ，０　／Ｒ０，０　を設定し、かかる正
規化された自己相関関数から、図２に示すように、自己
相関関数Ａｃｆ＝Ｒｍ，０　／Ｒ０，０　の値と画素位
置を示すｍとの関係が、パラメータＫに対応してどのよ
うに変化するかを得る。

【００１９】ここで、図２に示す理論的自己相関関数を
解析した結果から、自己相関関数Ａｃｆが所定値（例え
ば０．５　程度であり、相関距離変化に対して自己相関
関数値の変化が大きくなる点が好ましい。）になる値と
して、ｍ（相関距離）とＫとの関係（図３参照）が定ま
り、図１に示す濃度変動の実測データに基づく解析結果
からは、やはり自己相関関数Ａｃｆが前記所定値になる
値として、ラグ数Ｌａｇ（隔たり画素数）と濃度Ｄとの
関係（図４参照）が定まることになる。

【００２０】図３におけるｍは、相関距離（ラグ数Ｌａ
ｇ）に相当するから、上記２つの関係を組み合わせて、
濃度ＤとパラメータＫとの関係を求めることができる。 ここで、２つのパラメータＫ，Ｃと濃度実測値の標準偏
差（σ）とは、以下のような関係を満たす。 σ２　＝Ｒ０，０　＝Ｃ２　／（４πＫ２）Ｃ＝２√π
・Ｋ・σ・・・・（ｂ） 従って、各画像濃度Ｄ毎の濃度実測データに基づいて、
各画像濃度Ｄに対する濃度変化の標準偏差σを求めてお
けば、所定濃度Ｄに対応するパラメータＫの値と、これ
に対応する濃度実測値の標準偏差（σ）を、上記（ｂ）
式に代入することで、各色毎に所定濃度Ｄに対応するパ
ラメータＫ，Ｃを共に決定することができる。図５は、
画像濃度Ｄに対応するパラメータＫ，Ｃを上記のように
して求めた結果を、３原色の中の緑Ｇについて示したも
のである。

【００２１】上記のようにして決定されたパラメータＫ
は、３原色いずれの場合も画像濃度Ｄに対して略単調に
増加する。しかしながら、緑Ｇの場合は濃度Ｄ変化０．
８３〜２．３７に対してＫは０．２５〜０．６４の変化
幅を示すのに対し、青Ｂの場合は濃度Ｄ変化１．１８〜
２．９９に対してＫは０．２６〜０．３７の変化幅を示
し、赤Ｒの場合は濃度Ｄ変化０．６０〜１．７４に対し
てＫは０．２７〜０．４３の変化幅を示し、緑Ｇが濃度
Ｄ変化に対してパラメータＫの変化幅が最も大きく、他
の２色Ｂ，Ｒについては濃度Ｄ変化に対してパラメータ
Ｋの変化幅は比較的小さい。

【００２２】また、パラメータＣ（ノイズの平均的振幅
に関するパラメータ）は、濃度変化に対する変化幅が、
前記パラメータＫ（ノイズの平均的周波数特性に関する
パラメータ）よりも小さく、濃度依存性はパラメータＫ
よりも小さい。尚、前記パラメータＫ，Ｃは、試料数に
応じて濃度Ｄに対する離散的なデータとなるから、かか
るデータを補間することで、濃度Ｄに対して連続するパ
ラメータＫ，Ｃの関係を求めれば良い。

【００２３】上記のようにして、パラメータＫ，Ｃが各
色及び画像濃度Ｄに関連付けられて決定されれば、濃度
Ｄの関数として前記ノイズシミュレーションの予測式（
ａ）における予測係数ｈｐ，ｑ　を各色毎に求めること
ができることになる（電子通信学会論文誌‘８２／５　
　Ｖｏｌ．Ｊ６５−Ａ　ＮＯ．５　　吉田靖夫，小倉久
直，佐倉康夫　　、「スペクトル解析」朝倉書店　　１
９８９年発行　　日野幹雄　　著　　等参照）。従って
、この係数ｈｐ，ｑ　を用いて前記ノイズシミュレーシ
ョンの予測式（ａ）に基づいて、濃度Ｄに依存するノイ
ズ（粒状ノイズ）を各色毎にシミュレーションする（発
生させる）ことができる。

【００２４】上記に説明したパラメータＫ，Ｃを画像濃
度Ｄに依存させて決定する過程を、図６のフローチャー
トに概略示してある。即ち、まず均一露光によって平均
濃度の異なる複数のネガ試料を作成し（Ｓ１）、かかる
試料の濃度測定を行う（Ｓ２）。そして、濃度測定結果
に基づいて各試料（平均濃度）毎に自己相関関数を計算
する（Ｓ３）。一方、ノイズの理論的自己相関関数を設
定する（Ｓ４）。

【００２５】ここで、実測値に基づく自己相関関数から
は、自己相関関数が所定値（０．５　）になるときの値
として相関距離（ラグ数ｍ０．５　）と濃度Ｄとの関係
を求め（Ｓ５）、理論的自己相関関数からは、ノイズの
周波数特性を反映するパラメータＫと相関距離（ラグ数
ｍ０．５　）との関係を求める（Ｓ６）。上記の２つの
関係から、濃度ＤとパラメータＫとの関係が求められ（
Ｓ７）、かかる関係と、濃度測定値に基づいて求められ
る濃度Ｄと濃度変化の標準偏差σとの関係（Ｓ８）とか
ら、前記パラメータＣと濃度Ｄとの関係を求めることが
できる（Ｓ９）。

【００２６】ここで、以下のようにして、前記シミュレ
ーション結果が実際の画像に精度良く対応しているかの
確認を行った。まず、計算機上で、全面均一濃度（試料
の各平均濃度Ｄに対応する濃度）の元画像に対して、前
記シミュレーション結果のノイズＩｍ，ｎ　を各画素毎
に乗せる一方、試料の濃度計測時におけるアパーチャサ
イズに対応する画素数のデータを平均化処理することで
、濃度測定時の状態、即ち、試料の濃度分布状態を疑似
的に再現し、この疑似的に再現した濃度分布状態と実際
の濃度実測データと比較する。

【００２７】上記比較結果によると、シミュレーション
データは実測値に対して振幅が小さくなる傾向があった
。これは、試料の濃度測定時のアパーチャサイズを、シ
ミュレーションの画素サイズ（５×５μｍ）に対してそ
のまま適用したためと推定される。そこで、Ｓｅｌｗｙ
ｎの粒状則（σ＝１／√Ｓ；Ｓ＝面積）を適用すると、
前記（ｂ）式からσとパラメータＣとは比例関係にある
から、濃度測定のアパーチャサイズとシミュレーション
における画素サイズとの面積比に準じて前記パラメータ
Ｃを補正することで、前記シミュレーションデータが実
測値の振幅に近い特性となるようにすることが好ましく
、かかる補正によって振幅，周波数特性共に実測データ
に精度良く類似したシミュレーションデータを得ること
ができる。

【００２８】ここで、パラメータＫ，Ｃの濃度に対する
関係を得てからノイズを含んだ画像を疑似的に作成する
までの作業の流れを、図７に示すフローチャートに従っ
て概略説明する。まず、前記図６に示す流れに沿って、
パラメータＫ，Ｃそれぞれと画像濃度Ｄとの関係が特定
されると（Ｓ１１）、予め設定されている前記パラメー
タＫ，Ｃから前記ノイズシミュレーションの予測式（ａ
）における予測係数ｈｐ，ｑ　を算出するプログラムに
従って演算処理して（Ｓ１２）、２次元自己回帰表現を
用いたノイズシミュレーションの予測式で用いる予測係
数ｈｐ，ｑ　を３原色毎に画像濃度Ｄに対応させて設定
する（Ｓ１３）。

【００２９】一方、シミュレーションされたノイズを乗
せる濃度均一の元画像を得るために、例えば回転走査型
濃度計によってノイズ無し画像の濃度読取りを行い（Ｓ
１４）、かかる読取り結果を、ノイズ無しの元画像信号
として磁気記憶装置に一旦記憶させる（Ｓ１５）。そし
て、前記２次元自己回帰表現を用いたノイズシミュレー
ションの予測式を用いて前記元画像の濃度に対応する各
色毎のノイズを発生させ、前述のようにして予め記憶さ
せておいた均一濃度の元画像に、前記発生させたノイズ
を乗せ、ノイズがシミュレーションされた画像信号を得
る（Ｓ１６）。

【００３０】前述のシミュレーション評価では、データ
上の比較によって行ったが、更に、上記図７におけるＳ
１６で得られたノイズを乗せたシミュレーション画像信
号に基づいて実際にプリントさせて（Ｓ１７）、カラー
ネガフィルムに基づいてプリントさせた結果と目視で対
比して、ネガフィルムの粒状ノイズのシミュレーション
データと実際のネガフィルムの画像との対応性を確認し
た。

【００３１】前記目視観察は、前述のように、全面均一
濃度の元画像に対して、前記シミュレーション結果のノ
イズＩｍ，ｎ　を各画素毎に乗せて得られるシミュレー
ションデータを、レーザープリンタを用いて、印画紙上
に出力（露光）する一方、前記試料として作成したネガ
フィルムを、通常の焼付用のプリンタ（引伸し機）を用
いて、前記レーザープリンタと同じ印画紙上に同倍率で
拡大プリントした。そして、かかる目視観察によって、
上記のノイズシミュレーション方法による結果が、画像
の濃度変化及び色の違いに対して、実際の画像に良好に
対応することが確認された。

【００３２】尚、前記レーザープリンタは、例えば図８
に示すような構成となっている。３原色の光源１〜３か
らの照射される３原色光（赤，緑，青）は、それぞれに
備えられた光量変調素子４〜６で個別に変調された後、
各ミラー７〜１０（８，９は分光透過式ミラー）を介し
て合成され、該合成光が結像用レンズ１１によって回転
式ドラム１２に装着された印画紙１３上に結像される。

【００３３】前記光量変調素子４〜６は、磁気記憶装置
１４にＲ，Ｇ，Ｂの３原色毎の信号として記憶されたデ
ジタル画像信号をＤ／Ａ変換器１５でアナログに変換し
た結果に基づいて、それぞれの光量変調レベルが制御さ
れるようになっており、前記磁気記憶装置１４に、計算
機上でノイズが乗せられた画像信号を記憶させることで
、シミュレーション画像をプリントさせることができる
ようにしてある。

【００３４】更に、上記のようにして、パラメータＫ，
Ｃを共に画像濃度に依存させることによる効果の確認の
ために、 Ａ．パラメータＫ，Ｃを共に画像濃度に対応させて変化
Ｂ．パラメータＫのみを画像濃度に対応させて変化（但
し、σ（∝Ｃ／Ｋ）は一定） Ｃ．パラメータＣのみを画像濃度に対応させて変化の３
つのバリエーションでノイズシミュレーションを試み、
それぞれにシミュレーション結果を前記レーザープリン
タを用いてプリントさせ、ネガフィルムを用いたプリン
トの結果と比較した。この結果、ネガフィルムからのプ
リントに対する対応性が、Ａ＞Ｃ＞Ｂと順となり、やは
りＫ（ノイズの周波数特性を反映したパラメータ）及び
Ｃ（ノイズの振幅を反映したパラメータ）を共に画像濃
度に依存させて変化させることが、シミュレーション精
度を上げることが確認された。但し、パラメータＫの違
いは、周波数があるレベル以上になると目視では観察さ
れなくなり、寧ろ目視に与える影響は、パラメータＣの
方が大きくなる。

【００３５】尚、本実施例では画像形成媒体としてカラ
ーネガフィルムを例としたが、このカラーネガフィルム
の他、印画紙やプリンタにおける感光ドラムなどの他の
画像形成媒体上に形成された画像のノイズをシミュレー
ションするのにも有効であり、画像形成媒体をカラーネ
ガフィルムに限定するものではない。

【００３６】

【発明の効果】以上説明したように本発明にかかる画像
ノイズシミュレーション方法によると、ノイズの平均的
な振幅及び周波数特性をそれぞれ反映するパラメータを
画像濃度に関係付けて設定し、前記パラメータに基づい
てノイズシミュレーションの予測式の係数を決定するよ
うにした。これにより、ノイズを２次元で等方的に発生
させることが可能になり、また、ノイズの振幅及び周波
数特性が共に画像濃度に依存して変化するような画像形
成媒体上に形成される画像のノイズを、精度良くシミュ
レーションすることができるようになるという効果があ
る。更に、このシミュレーションを利用することにより
、上記画像形成媒体を設計・開発する上での目標設定，
改良効果の確認等がより容易に、かつ、効率的に行える
ようになる。

【図面の簡単な説明】

【図１】ラグ数（相関距離）と自己相関関数Ａｃｆの値
との関係を実測濃度データから求めた例を示す線図。

【図２】理論的自己相関関数における相関距離ｍと自己
相関関数値とパラメータＫとの関係を示す線図。

【図３】理論的自己相関関数から求めたパラメータＫと
相関距離との関係を示す線図。

【図４】実測濃度データに基づいて得られる画像濃度と
ラグ数（相関距離）との関係を示す線図。

【図５】画像濃度とパラメータＫ，Ｃとの関係を示す線
図。

【図６】パラメータＫ，Ｃと画像濃度Ｄとの関係を特定
する処理内容を示すフローチャート。

【図７】予測式に基づきノイズの発生させノイズを含ん
だ画像を疑似的に発生させる処理内容を示すフローチャ
ート。

【図８】レーザープリンタの構成例を示すブロック図。

【符号の説明】

Ｋ　　　　ノイズの平均的周波数特性を反映するパラメ
ータＣ　　　　ノイズの平均的振幅を反映するパラメー
タＡｃｆ　　　　自己相関関数 ｍ，ｎ　　　　画素座標位置 Ｄ　　　　画像濃度 σ　　　　標準偏差

Claims (2)

【特許請求の範囲】
1. 【請求項１】画像形成媒体上に形成される画像のノイズ
   の平均的な振幅及び周波数特性をそれぞれ反映するパラ
   メータと、画像形成媒体における画像濃度との関係を求
   めた後、前記パラメータを用いてノイズシミュレーショ
   ンの予測式の係数を決定し、該係数を用いた前記予測式
   に基づいて前記画像形成媒体上に形成された画像のノイ
   ズをシミュレーションする画像ノイズシミュレーション
   方法。
2. 【請求項２】前記パラメータがノイズの平均的な振幅を
   反映するパラメータと周波数特性を反映するパラメータ
   との少なくとも２つからなり、前記画像形成媒体の平均
   濃度が異なる複数試料の濃度実測値から平均濃度別に実
   測値による自己相関関数を設定する一方、前記少なくと
   も２つのパラメータを含めて設定されたノイズのモデル
   から定まる理論的自己相関関数を算出し、実測値による
   自己相関関数と理論的自己相関関数との間の誤差の総和
   を最小化することで前記少なくとも２つのパラメータの
   値を前記各平均濃度別に求めることにより、ノイズの平
   均的な振幅及び周波数特性をそれぞれ反映する前記少な
   くとも２つのパラメータと画像濃度との関係を求めるよ
   う構成したことを特徴とする請求項１記載の画像ノイズ
   シミュレーション方法。