"任意のシーンについてのコンピュータグラフィックス画像を、光線追跡法によって生成する画像処理装置であって、光源の分光特性および配光特性を示す光源特性データを入力する光源特性入力手段と、前記シーン内におけるオブジェクトの分光特性を示すオブジェクト特性データを入力するオブジェクト特性入力手段と、前記光源特性データおよび前記オブジェクト特性データに基づいて、互いに異なる波長の光線ごとに光線追跡を行うことによって、前記シーンにおける各画素に対応する分光放射輝度を算出する光線追跡手段と、前記各画素に対応する分光放射輝度から色信号値を算出して前記コンピュータグラフィックス画像を生成する色信号値算出手段と、を有することを特徴とする画像処理装置。前記光線追跡手段は光線ごとに、前記光源特性データと、該光線が前記シーン内において最初に交差する初回交差オブジェクトの前記オブジェクト特性データに基づいて最大交差回数を設定し、該シーン内における該光線とオブジェクトとの交差回数が該最大交差回数を超えないように光線追跡を行うことを特徴とする請求項1に記載の画像処理装置。前記光線追跡手段は、前記初回交差オブジェクトとの交差後に残存する前記光線の強度が大きいほど、前記最大交差回数を大きく設定することを特徴とする請求項2に記載の画像処理装置。前記光線追跡手段は、前記初回交差オブジェクトに対する前記光線の反射成分および透過成分のそれぞれについて、前記最大交差回数を設定することを特徴とする請求項2または3に記載の画像処理装置。前記光線追跡手段は光線ごとに、前記シーン内で該光線と交差するオブジェクトごとに、前記光源特性データと当該オブジェクトの前記オブジェクト特性データに基づいて発射光線数を設定し、該光線と該オブジェクトとの交点から該発射光線数の光線を飛ばして光線追跡を行うことを特徴とする請求項1乃至4のいずれか1項に記載の画像処理装置。前記光線追跡手段は、前記オブジェクトとの交差後に残存する前記光線の強度が大きいほど、前記発射光線数を大きく設定することを特徴とする請求項5に記載の画像処理装置。前記光線追跡手段は、前記オブジェクトに対する前記光線の反射成分および透過成分のそれぞれについて、前記発射光線数を設定することを特徴とする請求項5または6に記載の画像処理装置。さらに、視覚特性を入力する視覚特性入力手段を有し、前記色信号値算出手段は、前記各画素に対応する分光放射輝度および前記視覚特性から、当該画素の色信号値を算出することを特徴とする請求項1乃至7のいずれか1項に記載の画像処理装置。前記光源特性データは、光源から任意の方向に放射される光の強度分布を示す情報を、前記光源の配光特性として含むことを特徴とする請求項1乃至8のいずれか1項に記載の画像処理装置。前記光源特性データは、光源の分光放射輝度を示す情報を前記光源の分光特性として含み、前記オブジェクト特性データは、オブジェクトの反射特性と透過特性および屈折率を示す情報を前記オブジェクトの分光特性として含むことを特徴とする請求項9に記載の画像処理装置。前記オブジェクトの反射特性の情報は、該オブジェクトに対する光線の入射角と反射角に応じた特性を示すことを特徴とする請求項10に記載の画像処理装置。前記オブジェクトの透過特性の情報は、該オブジェクトに対する光線の入射角と屈折角に応じた特性を示すことを特徴とする請求項10に記載の画像処理装置。光源特性入力手段、オブジェクト特性入力手段、光線追跡手段、および色信号値算出手段、を有し、任意のシーンについてのコンピュータグラフィックス画像を、光線追跡法によって生成する画像処理装置における画像処理方法であって、前記光源特性入力手段が、光源の分光特性および配光特性を示す光源特性データを入力し、前記オブジェクト特性入力手段が、前記シーン内におけるオブジェクトの分光特性を示すオブジェクト特性データを入力し、前記光線追跡手段が、前記光源特性データおよび前記オブジェクト特性データに基づいて、互いに異なる波長の光線ごとに光線追跡を行うことによって、前記シーンにおける各画素に対応する分光放射輝度を算出し、前記色信号値算出手段が、前記各画素に対応する分光放射輝度から色信号値を算出して前記コンピュータグラフィックス画像を生成することを特徴とする画像処理方法。コンピュータ装置で実行されることにより、該コンピュータ装置を請求項1乃至12のいずれか1項に記載の画像処理装置の各手段として機能させるためのプログラム。しかしながら、上記CG画像を生成する従来の各手法においては、光源の配光特性が考慮されていないため、例えば蛍光灯とLED照明のように、配光特性の異なる光源によってシーンに生じる陰影の変化を再現できないという問題があった。また、一般的な光線追跡法では光を波長に依存しない一様な光線として扱うため、オブジェクトを透過する光の屈折のように、波長によって異なるふるまいをする物理現象を正しく再現することができない、という問題があった。本発明は上述した問題を解決するために、光源やオブジェクトの分光特性に基づいた光線追跡を行うことで、シーンを忠実に再現したコンピュータグラフィックス画像を生成することを目的とする。上記問題を解決するための一手段として、本発明の画像処理装置は以下の構成を備える。すなわち、任意のシーンについてのコンピュータグラフィックス画像を、光線追跡法によって生成する画像処理装置であって、光源の分光特性および配光特性を示す光源特性データを入力する光源特性入力手段と、前記シーン内におけるオブジェクトの分光特性を示すオブジェクト特性データを入力するオブジェクト特性入力手段と、前記光源特性データおよび前記オブジェクト特性データに基づいて、互いに異なる波長の光線ごとに光線追跡を行うことによって、前記シーンにおける各画素に対応する分光放射輝度を算出する光線追跡手段と、前記各画素に対応する分光放射輝度から色信号値を算出して前記コンピュータグラフィックス画像を生成する色信号値算出手段と、を有することを特徴とする。本発明によれば、光源やオブジェクトの分光特性に基づいた光線追跡を行うことで、シーンを忠実に再現したコンピュータグラフィックス画像を生成することが可能となる。以下、本発明に係る実施形態について、図面を参照して説明する。なお、以下の実施の形態は特許請求の範囲に関る本発明を限定するものではなく、また、本実施の形態で説明されている特徴の組み合わせの全てが本発明の解決手段に必須のものとは限らない。<第1実施形態>本実施形態においては、光源の配光特性と分光特性、及びオブジェクトの分光特性を考慮して、光線追跡法によりシーンを忠実に再現したコンピュータグラフィックス画像を生成する処理について説明する。まず、波長によって異なる光の振る舞いについて、図1を用いて説明する。光はオブジェクトにあたると、その表面で反射する光と透過する光に分離する。オブジェクトを透過する光は、その波長λに応じた屈折率によって屈折角が変化する。例えば図1に示すように、視点を開始点として、屈折率n1の空間を通過した光線が屈折率n2のオブジェクト101を透過する場合を考える。このとき、オブジェクト101に対する光線の入射角をφとすると、屈折角θは以下の式で示される。θ=sin1n2sinφこのように屈折角は光の波長λの関数になるため、光線追跡法によって追跡する光の経路が波長λによって変化することになる。従って、図1でオブジェクト101を透過した光がさらに屈折率n3のオブジェクト102を透過する等、あるシーンにおける光の透過を忠実に再現するためには、光を分光として扱う必要がある。●装置構成本実施形態における画像処理装置の構成例について、図2を用いて説明する。図2に示すように本実施形態の画像処理装置は、CPU201、メインメモリ202、HDD203、汎用インターフェイス204、モニタ205、メインバス206を備える。IF204は、分光測定器207、外部記憶装置208等をメインバス206に接続する。以下、CPU201がHDD203に格納された各種ソフトウェアを動作させることで実現する各種処理について説明する。まずCPU201は、HDD203に格納されている画像処理アプリケーションを起動し、メインメモリ202に展開する。続いて、HDD203に格納されている各種データや、IF204に接続された分光測定器207や外部記憶装置208から入力されるデータ等が、CPU201からの指令に基づき、メインバス206経由でメインメモリ202に転送される。さらに、画像処理アプリケーション内の処理にしたがって、メインメモリ202に格納されているデータにCPU201からの指令に基づく各種演算が施され、該演算結果をメインバス206経由でモニタ205に表示したり、HDD203に格納したりする。以下、上記構成において、CPU201からの指令に基づき画像処理アプリケーションが画像を生成してモニタ205に表示するまでの画像処理について、詳細に説明する。図3は、本実施形態における画像処理装置301の論理構成を示すブロック図であるが、同図に示す各構成は、上述したように画像処理アプリケーションとして実現される。図3において、光源特性入力部302は、シーンを照明する光源の配光特性および分光放射輝度を示す光源特性データ307を入力する。オブジェクト特性入力部303は、シーンを構成する各オブジェクトの反射特性、透過特性及び屈折率を示すオブジェクト特性データ308を入力する。視覚特性入力部304は、シーンを観察する際の視覚特性データ309を入力する。尚、光源特性データ307、オブジェクト特性データ308、視覚特性データ309は、HDD203からの入力、またはIF204を介して分光測定器207、外部記憶装置208からIF204を介した入力等、所望の手段によって入力される。分光放射輝度算出部305は、光源特性入力部302及びオブジェクト特性入力部303から入力された情報に基づき、CPU201からの指示に従って、シーンのCG画像における各画素の分光放射輝度310を算出する。色信号値算出部306は、分光放射輝度算出部305で算出した各画素の分光放射輝度310と、視覚特性入力部304から入力された視覚特性データ309に基づき、CPU201からの指示に従ってシーンのCG画像における各画素の色信号値を算出する。該色信号値は例えばRGB値として算出され、画像データ311として出力される。画像データ311は、HDD203等の記憶領域に格納されるとともに、モニタ205に表示される。●シーンCG画像の生成処理以下、本実施形態の画像処理装置301において、シーンのCG画像を作成する画像生成処理について、図4のフローチャートを用いて説明する。まずS401で、光源特性入力部302が光源特性データ307を入力する。この光源特性データ307は上述したように、シーンを照明する光源の配光特性及び分光放射輝度を示す。ここで、光源の配光特性について図5を用いて説明する。配光特性は、光源から放射される光の各方向への強度分布を示すものであり、光源の指向性の違いによって一つの光源がシーンを照らす範囲が異なる。図5は、指向性の低い蛍光灯などの光源における配光特性の例であり、シーンを広く照らしていることが分かる。一方図5は、指向性の高いLEDなどの光源における配光特性の例であり、シーンの一部だけが照らされていることが分かる。このように、光源の配光特性が異なるとシーンに生じる陰影が大きく変化するので、光源の配光特性はシーンを忠実に再現するために重要なファクタとなる。尚、配光特性はHDD203等に記憶されているデータを入力しても良いし、IF204に接続された分光測定器207で光源を測定したデータを入力しても良い。また、入力するデータ形式としては、光源からの光の出力角と光の強度との対応が得られれば、どのような形式であっても構わない。例えば、出力角に対する光の強度の関数形式であって良いし、出力角と光の強度の対応を示すルックアップテーブルの形式であっても良い。次に、分光放射輝度について図6を用いて説明する。光源の分光放射輝度は、波長毎の輝度分布を示すものであり、光源によって特性が異なる。図6は蛍光灯の分光放射輝度分布の例であり、図6はLED照明の分光放射輝度分布の例である。このように、光源の種類が異なると分光放射輝度分布も変化し、その結果、オブジェクトの色の見えが変化するので、光源の分光特性はシーンを忠実に再現するために重要なファクタとなる。尚、分光特性も配光特性と同様に、HDD203等に記憶されているデータを入力しても良いし、IF204に接続された分光測定器207で光源を測定したデータを入力しても良い。また、入力するデータ形式としては、可視光域の波長と放射輝度との対応が得られれば、どのような形式であっても構わない。例えば、光の波長に対する放射輝度の関数であっても良いし、波長と放射輝度の対応を示すLUTの形式でも良い。次にS402で、オブジェクト特性入力部303がオブジェクト特性データ308を入力する。このオブジェクト特性データ308は上述したように、シーン内の各オブジェクトの反射特性、透過特性及び屈折率を示す。ここで、オブジェクトの反射特性および透過特性について、図7を用いて説明する。一般に、オブジェクトの反射特性は双方向反射率分布関数frで表される。オブジェクト上の位置xにおいて、ある方向ωiから入射する光の波長をλi、強度をLiとし、ある方向ωrに反射する光の波長をλr、強度をdLrとすると、反射特性frは以下の式で示される。fr=dLr/Lidω...一般に、反射光の強度分布はオブジェクトの表面形状に依存し、オブジェクト表面が平滑なほど、鏡面反射方向すなわちωi=ωr方向の反射強度が高くなり、それ以外の拡散反射方向への反射強度は小さくなる。尚、オブジェクトの反射特性frの取得方法としては、HDD203等に記憶されているデータを入力しても良いし、IF204に接続された分光測定器207で測定したデータを入力しても良い。また、反射特性frのデータ形式としては、オブジェクトへの光の入射角と反射角に対応した反射率が得られれば、どのような形式であっても構わない。例えば、上記式のように一意の関数として入力しても良いし、一般に拡散反射は等方的に分布するため、ωi=45°ωr=0°の反射特性を拡散反射に適用し、鏡面反射は別途LUTにする等、分割して入力しても構わない。また、オブジェクトを透過する光は双方向透過率分布関数ftで表される。オブジェクト上の位置xにおいて、ある方向ωiから入射する光の波長をλi、強度をLiとし、ある方向ωtに透過する光の波長をλt、強度をdLtとすると、透過特性ftは以下の式で示される。ft=dLt/Lidω...一般に、透過光の強度分布は屈折率に依存するため、ωtが上記式を満たす場合に、透過強度が最大になる。尚、オブジェクトの透過特性ftの取得方法としても、HDD203等に記憶されているデータを入力しても良いし、IF204に接続された分光測定器207で測定したデータを入力しても良い。また、透過特性ftのデータ形式も、オブジェクトへの光の入射角と屈折角に対応した透過率が得られれば、どのような形式であっても構わない。例えば、上記式のように一意の関数として入力しても良いし、ωi=ωr=0°の透過特性をωtが式の屈折角θに等しい場合に適用し、それ以外の透過は別途LUTにする等、分割して入力しても良い。次にS403で、視覚特性入力部304が視覚特性データ309を入力する。この視覚特性データ309は、分光放射輝度を画像の色信号値に変換する際に用いられ、ここで視覚特性データ309の例として、図8にCIEXYZ表色系で用いる等色関数を示す。等色関数は、可視光域の光に対する目の感度を示している。尚、視覚特性データ309は、HDD203等に予め記憶しておいて入力すれば良い。次にS404で分光放射輝度算出部305が、当該シーンのCGとしての出力画像を順次走査して、分光放射輝度の算出対象となる画素を決定する。なお、この対象画素の選択順は任意である。次にS405で分光放射輝度算出部305において、上記S404で選択された対象画素における分光放射輝度を算出する。この分光放射輝度の算出には光線追跡法を用いるが、処理の詳細については後述する。次にS406で色信号値算出部306において、対象画素の色信号値を算出する。この色信号値の算出方法の詳細については後述する。そしてS407で、出力画像の全画素について分光放射輝度の算出が終了したか否かを判定する。終了していれば生成した画像データを出力して処理を終了するが、未終了であればS404に戻り、次の対象画素を決定して処理を継続する。●分光放射輝度算出処理以下、上記S405における分光放射輝度の算出処理について、図9を用いて詳細に説明する。上述したように、本実施形態における分光放射輝度の算出には光線追跡法を用いる。この手法によれば、視点が画角内の各方向から受け取るはずの光線を、仮想的に光源まで追跡することでCG画像の描画を行う。詳細には、視点から描画する各画素の方向へ直線を伸ばし、オブジェクトとの交差の有無を数学的に判定する。そして、反射や透過により発生した光線について、再帰的に探索を繰り返し、光源と交差した場合に計算を終了する。すなわち本実施形態では、図9のフローチャートに示す処理を再帰的に実行することで複数の光線を追跡し、光線ごとに算出した分光放射輝度の平均により、対象画素における分光放射輝度を算出する。なお、図9では波長λの単色光についての分光放射輝度算出処理を説明するが、実際には同様の処理を所定の波長毎に行うことで、当該対象画素についての分光放射輝度が得られる。なお、処理対象となる波長は全波長域から等間隔で抽出されるものとする。例えば380~730nmの可視域において、10nm間隔で36波長を処理対象波長として抽出する。まずS901で、追跡の開始点から飛ばした波長λの光線がシーンにおいて交差する点を探索する。ここで追跡の開始点は、当該フローチャートが初回の実行であれば対象画素であり、再帰実行であれば光線とオブジェクトとの交点である。なお、この交点の探索方法は任意であり、例えば、シーン内の全てのオブジェクトについて順次交差判定を行って、追跡の開始点に最も近い交点を求めれば良い。次にS902で、S901で探索された光線と交差する点が、光源との交差であるか、またはその他のオブジェクトとの交差であるかを判定する。光線が光源と交差している場合にはS903に移行し、光線の強度と交差した光源の光源特性データ307から、波長λについての分光放射輝度を算出して、当該波長についての処理を終了する。ここで光源特性データ307は光源の配光特性および分光放射輝度を示すものであるから、光線の分光放射輝度E0は以下の式によって示される。なお下式においては、光線の強度をP0、光線と光源のなす角をθ、光源からθ方向の配光特性をQ、光源の分光放射輝度をSとする。E0=P0QS...一方、光線が光源以外のオブジェクトと交差している場合にはS904に移行し、まず、光線が交差したオブジェクトのオブジェクト特性データ308を取得する。次にS905で、光線とオブジェクトとの交点から飛ばす光線数を決定する。この発射光線数は任意に決定して構わないが、数が多いほどノイズの少ない画像が得られる。以下、発射光線数に応じた回数分、S906及びS907の処理をループする。S906では、発射光線の方向および強度を算出する。光線の方向は発射光線数に応じて均等に分布すれば良く、該交点に対する反射光の強度Pr、透過光の強度Ptはそれぞれ、以下の式によって示される。なお、下式において、P0はオブジェクトと交差する前の光線の強度、frは上記式に示したオブジェクトの反射特性、ftは上記式に示したオブジェクトの透過特性である。Pr=P0fr...Pt=P0ft...次にS907で図9のフローチャートに示す処理を再帰的に呼び出すことによって、当該発射光線に対する光線追跡処理が行われる。このS907の処理はすなわち、再帰処理によって最終的に交差するオブジェクトがなくなり、当該発射光線についての分光放射輝度がS903で算出された時点で、終了する。このように、発射光線についてS907の再帰処理によって最終的に上記式で算出された分光放射輝度E0が、当該発射光線についての分光放射輝度Enとなる。そしてS908で、全ての発射光線について上記S906及びS907による光源追跡が終了したと判定されると、次にS909において、各発射光線の分光放射輝度Enの平均を算出する。この平均分光放射輝度値が、対象画素における追跡対象波長λについての分光放射輝度Exである。すなわち、光線の総数をN、各発射光線の分光放射輝度をEnとすると、対象画素における追跡対象波長λの分光放射輝度Exは、以下の式で算出される。なお、式においてΣはn0からN1までの整数に関する総和を示す。Ex=ΣEn...以上、S405では図9に示す光線追跡処理を再帰的に実行することによって、対象画素における分光放射輝度Exを算出する。●色信号値算出処理以下、上記S406における画素ごとの色信号値の算出処理について、図10を用いて詳細に説明する。尚、本実施形態では最終的にsRGB形式の色信号値を出力するとして説明するが、AdobeRGB等の他形式による色信号値であっても構わない。まずS1001で、S403で入力した視覚特性データ309として、例えばCIEXYZ表色系の等色関数xyzを取得する。次にS1002で、上記S405で算出した、画素ごとの分光放射輝度Exを取得する。そしてS1003で、取得した視覚特性xyzと分光放射輝度Exから、以下の~式によりCIEXYZ値を算出する。そしてS1004において、出力画像の色信号値を算出する。ここでは出力画像をsRGB形式とするるため、出力信号RoutGoutBoutを以下の~式によって算出する。以上、図4に示す画像生成処理によって、シーンのCGである出力画像が形成される。以上説明したように本実施形態によれば、シーンにおける光線追跡を、光源の配光特性と分光特性、及びオブジェクトの分光特性を利用して行うことで、該シーンを忠実に再現したCG画像を生成することができる。<変形例>以下、第1実施形態の変形例について説明する。上述した第1実施形態では、光源特性データ307によって示される配光特性を、光の強度分布として扱う例を示した。しかしながら、光の強度分布は波長によって異なる場合がある。例えばLED照明の場合、青色を発光するLEDと、黄色を発光する蛍光体を組み合わせて光源を作る場合がある。この時、図5にも示したようにLEDは指向性の高い光を発光するため、蛍光体と組み合わせた際に色ムラが発生する。すなわち、図11に示すように、光源の配光特性が波長によって変化するため、分光放射輝度も光線と光源のなす角度θに依存することになる。そこで、光源情報として波長によって異なる配光特性を入力することが考えられる。この場合における分光放射輝度の算出処理は第1実施形態とほぼ同様であるが、図9におけるS903の最終的な分光放射輝度の算出処理のみが異なる。すなわち、光源が波長によって異なる配光特性を有する場合、その配光特性は、波長λと光線と光線のなす角θの関数Qとして表わされる。この場合、第1実施形態で用いた光源の分光放射輝度Sは、配光特性Qから得られるため不要になる。従って、S903における光線の分光放射輝度E0の算出は、以下の式に従う。E0=P0Q...以上のように、光源の配光特性が波長によって異なる場合であっても、上記第1実施形態と同様にシーンを忠実に再現したCG画像を生成することができる。<第2実施形態>以下、本発明に係る第2実施形態について説明する。上述した第1実施形態では、光源の配光特性や分光特性、およびオブジェクトの分光特性に基づいた光線追跡を行うことによって、該シーンを忠実に再現したCG画像を生成する例を示した。しかしながら、第1実施形態のように分光による光線追跡を行う場合、波長数に応じて演算量が増加し、処理速度が低下してしまう。一般に光線追跡法においては、1本の光線におけるオブジェクトとの交差回数や、追跡を行う光線数を削減することで演算量を抑制して処理速度を上げることが可能である反面、生成画像にノイズが発生して画質が低下してしまう。ここで、上記特許文献1に記載された空間分割法を用いた光線追跡の効率化手法を、分光による光線追跡に適用することを考える。この場合に、更なる高速化のために光の波長によらず交差回数や光線数を削減すると、やはりノイズが発生してしまう。そこで第2実施形態においては、上述した第1実施形態に加え、光源やオブジェクトの分光特性に応じて光線追跡時における演算量を制御することで、処理の高速化とノイズ低減化を実現する。なお、第2実施形態における画像処理装置の構成は、上述した第1実施形態で図2、図3に示した構成と同様であるため、同符号を参照するとして説明を省略する。●シーンCG画像の生成処理以下、第2実施形態の画像処理装置301において、シーンのCG画像を作成する画像生成処理について、図12のフローチャートを用いて説明する。まずS1401で、光源特性入力部302が光源特性データ307を入力する。この光源特性データ307は、第1実施形態で図5図6を用いて説明したような、シーンを照明する光源の配光特性及び分光放射輝度を示す。次にS1402で、オブジェクト特性入力部303がオブジェクト特性データ308を入力する。このオブジェクト特性データ308は、第1実施形態で図7を用いて説明したように、シーン内の各オブジェクトの反射特性、透過特性及び屈折率を示す。次にS1403で、視覚特性入力部304が視覚特性データ309を入力する。この視覚特性データ309は、分光放射輝度を画像の色信号値に変換する際に用いられ、例えば図8に示す等色関数である。次にS1404で分光放射輝度算出部305が、当該シーンのCGとしての出力画像を順次走査して、分光放射輝度の算出対象となる画素を決定する。なお、この対象画素の選択順は任意である。次にS1405で分光放射輝度算出部305が、対象画素から発せられた光線が経路追跡中のオブジェクトと交差した回数をカウントするための変数である、交差回数Nを0に初期化する。次にS1406で分光放射輝度算出部305が、波長λごとの交差回数の上限となる最大交差回数Mを初期化する。最大交差回数Mは、対象画素からの光線が光源に到達するまでに多数のオブジェクトと交差する場合に、光源への到達前に交差回数に応じて処理を打ち切るための閾値となる変数である。本来、光線が光源に到達するまで追跡を行うことが望ましいが、それには膨大な計算が必要となるため、この最大交差回数Mまでに光線が光源に到達しない場合には光線追跡を打ち切ることで計算量を削減し、計算時間の高速化を実現する。最大交差回数Mの初期値としては、光の波長λごとに1以上の任意の値を設定する。なお、最大交差回数Mが小さいと、計算を高速化できる半面、生成結果画像にノイズが発生してしまうため、求められる計算時間とノイズ量を考慮して、初期値を設定することが望ましい。また、最大交差回数Mの初期値としては、光の波長λによらず等しい値を設定することが望ましい。次にS1407で分光放射輝度算出部305において、上記S1404で選択された対象画素における分光放射輝度を算出する。この分光放射輝度の算出には光線追跡法を用いるが、処理の詳細については後述する。次にS1408で色信号値算出部306において、対象画素の色信号値を算出する。この色信号値の算出方法は、上述した第1実施形態で図10を用いて説明した方法と同様であるため、説明を省略する。そしてS1409で、出力画像の全画素について分光放射輝度の算出が終了したか否かを判定する。終了していれば生成した画像データを出力して処理を終了するが、未終了であればS1404に戻り、次の対象画素を決定して処理を継続する。●分光放射輝度算出処理以下、上記S1407における分光放射輝度の算出処理について、図13を用いて詳細に説明する。第2実施形態においても、図13のフローチャートに示す処理を再帰的に実行することで複数の光線を追跡し、光線ごとに算出した分光放射輝度の平均により、対象画素における分光放射輝度を算出する。なお、図13では波長λの単色光についての分光放射輝度算出処理を説明するが、実際には同様の処理を所定の波長毎に行うことで、当該対象画素についての分光放射輝度が得られる。まずS1901で、追跡の開始点から飛ばした波長λの光線がシーンにおいて交差する点を探索する。ここで追跡の開始点は、当該フローチャートが初回の実行であれば対象画素であり、再帰実行であれば光線とオブジェクトとの交点である。なお、この交点の探索方法は任意であり、例えば、シーン内の全てのオブジェクトについて順次交差判定を行って、追跡の開始点に最も近い交点を求めれば良い。次にS1902で、S1901で探索された光線と交差する点が、光源との交差であるか、またはその他のオブジェクトとの交差であるかを判定する。光線が光源と交差している場合にはS1914に移行し、光線の強度と交差した光源の光源特性データ307から、波長λについての分光放射輝度E0を上記式より算出して、当該波長についての処理を終了する。一方、S1902で光線が光源以外のオブジェクトと交差している場合にはS1903に移行し、光線に対する交差回数Nに1を加えて値を更新する。そしてS1904で、光線が交差したオブジェクトのオブジェクト特性データ308を取得する。次にS1905で、光線の交差回数Nと光線の波長における最大交差回数Mの大小を判定する。交差回数Nが最大交差回数Mを超える場合は、光線が光源に未到達であっても光線追跡処理を終えるためにS1913に移行する。一方、交差回数Nが最大交差回数M以下の場合は、光線とオブジェクトの交点から新たに複数の光線を飛ばし、各光線についての光線追跡処理を行うためにS1906に移行する。このように第2実施形態では、予め設定した最大交差回数Mまでに光線が光源に到達しない場合は、光源到達前に処理を打ち切ることで、処理の高速化を実現する。しかし、上述したように最大交差回数Mを小さくすると、計算を高速化できる半面、生成画像中にノイズが多くなってしまう。そこで第2実施形態では詳細は後述するが、光線が最初に交差したオブジェクトの反射率や透過率と光源の放射輝度、およびS1406で設定された初期値に応じて、光の波長ごとに最適な最大交差回数Mを設定する。S1906では、光線とオブジェクトとの交点から新たに飛ばす発射光線数を決定する。この発射光線数は任意に決定して構わないが、数が多いほどノイズの少ない画像が得られる反面、計算量が増大する。なお、新たに飛ばす各光線の交差回数N、および、最大交差回数Mの値は、元の光線の値を引き継ぐものとする。以下、発射光線数に応じた回数分、S1907~S1910の処理をループする。S1907では、発射光線の方向および強度を算出する。光線の方向は発射光線数に応じて均等に分布すれば良く、該交点に対する反射光の強度Pr、透過光の強度Ptはそれぞれ、第1実施形態と同様に式によって示される。次にS1908では、光線における現在の交差回数Nの値が1であるか否かの判定を行う。交差回数Nが1、すなわち初回の交差である場合はS1909に移行し、1以外、すなわち2回目以降の交差である場合はS1910へ移行する。S1909では、光源の分光特性とオブジェクトの分光特性に応じて、波長λにおける最大交差回数Mを算出し、更新する。ここで、最大交差回数Mの設定方法について詳細に説明する。最大交差回数Mは、シーン内のオブジェクトに対する光線の反射成分および透過成分のそれぞれについて設定される。例えば、オブジェクトの反射率が小さい波長の光線では、該オブジェクトとの交差によって放射輝度が大きく減衰してしまうため、該交差後の他オブジェクトとの交差による放射輝度の変化量は相対的に小さくなる。そのため、少ない交差回数の後に得られる放射輝度と、十分大きな交差回数の後に得られる放射輝度とで、差がほとんど生じない。一方、オブジェクトの反射率が大きい波長の光線では、該オブジェクトとの交差後も放射輝度があまり減衰しないため、該交差後の他オブジェクトとの交差による放射輝度の変化量は相対的に大きくなる。そのため、少ない交差回数の後に得られる放射輝度と、十分大きな交差回数の後に得られる放射輝度とで、大きな差が生じ、すなわちノイズが発生してしまう。このような特性は、オブジェクトに対する透過成分についても同様である。そこで第2実施形態では、ある波長の光線が最初に交差するオブジェクトについて、その反射率または透過率が大きければ最大交差回数Mを大きく設定し、逆に小さければ最大交差回数Mを小さく設定する。また、光源からの放射輝度の大きさはそのまま、交差後の放射輝度の大きさに反映されるため、最大交差回数Mの設定時にも、光源の放射輝度が参照される。すなわち、初回交差オブジェクトとの交差後に残存する光線の強度が大きいほど、最大交差回数Mを大きく設定する。このように、光源の放射輝度とオブジェクトの分光特性に応じて最大交差回数Mを設定することで、演算量の削減とノイズ発生の抑制を両立することが可能となる。ここで図14を用いて、特にオブジェクトからの反射成分についての最大交差回数Mの設定方法の概念について説明する。ある入射角度、出射角度における、光線と交差するオブジェクトの分光反射率が図14で示される場合、例えば480nmのように反射率が高い波長では最大交差回数Mを大きく設定する。一方、例えば680nmのように反射率が低い波長では最大交差回数Mを小さく設定する。具体的には、オブジェクトに対する反射光、透過光のそれぞれの最大交差回数Mr,Mtは、以下の式で記述される。下式によれば最大交差回数Mr,Mtは、S406で初期化した最大交差回数Mと、上記式で表される分光特性frft、および光源の分光放射輝度Sの関数gとして表わされる。Mr=gfrS...Mt=gftS...尚、関数gとしては、反射率や透過率が高い波長において最大交差回数Mが大きくなるような関数であれば良く、具体的には、以下の式や式に示すような関数を用いることができる。gftS=M×ft×S...gftS=M×√×S...次にS1910で図13のフローチャートに示す処理を再帰的に呼び出すことによって、当該発射光線に対する光線追跡処理が行われる。このS1910の処理はすなわち、再帰処理によって最終的に交差するオブジェクトがなくなるか、または発射光線とオブジェクトの交差回数Nが最大交差回数Mを超えた場合に、終了する。このように、発射光線についてS1910の再帰処理によって最終的に上記式で算出された分光放射輝度E0が、当該発射光線についての分光放射輝度Enとなる。そしてS1911で、全ての発射光線について上記S1907~S1910による光源追跡が終了したと判定されると、次にS1912において、各発射光線の分光放射輝度Enの平均を、上記式により算出する。この平均分光放射輝度値が、対象画素における追跡対象波長λについての分光放射輝度Exである。一方、S1905で、光線とオブジェクトの交差回数Nが最大交差回数Mを超える場合にはS1913に移行し、光線追跡処理を打ち切る場合の光線の方向と強度を算出する。具体的には、光線の方向としては交点から光源の方向を設定する。また光線の強度Pとしては、交点から光源までの直線上に他のオブジェクトが存在する場合には0を設定し、他のオブジェクトが存在しない場合には以下の式に従い、反射のみを考慮する。なお、下式において、P0はオブジェクトと交差する前の光線の強度、frは上記式に示したオブジェクトの反射特性である。P=P0fr...そして、S1913で光線追跡処理を打ち切る場合の光線の方向と強度が算出されると、次にS1914で、波長λについての分光放射輝度E0を上記式より算出して、当該波長についての処理を終了する。この演算には、S1913で算出された光線の強度Pと、光源の配光特性Qおよび分光放射輝度Sが適用される。以上説明したように第2実施形態によれば、分光による光線追跡を行う際に、光源の放射輝度と、オブジェクトの分光特性を考慮して光線の波長ごとに交差回数を制限する。これにより、波長ごとに必要最低限の交差による演算が可能となるため、計算時間の高速化とノイズ発生の低減を実現することが可能になる。<第3実施形態>以下、本発明に係る第3実施形態について説明する。上述した第2実施形態では、光源とオブジェクトの分光特性に応じて、光線の波長ごとに最大交差回数を設定する例を示した。第3実施形態では、これら分光特性に応じて、交点から飛ばす光線数を波長ごとに設定することで、計算時間の高速化とノイズの発生を低減する例を示す。なお、第3実施形態における画像処理装置の構成は、上述した第1実施形態で図2、図3に示した構成と同様であるため、同符号を参照するとして説明を省略する。●分光放射輝度算出処理第2実施形態の画像処理装置301においてCG画像を作成する画像生成処理は、第2実施形態で図12のフローチャートに示した処理とほぼ同様であり、S1407における分光放射輝度の算出処理が異なるのみである。したがって、以下では第3実施形態における分光放射輝度の算出処理について、図15を用いて説明する。第3実施形態においても、図15のフローチャートに示す処理を再帰的に実行することで複数の光線を追跡し、光線ごとに算出した分光放射輝度の平均により、対象画素における分光放射輝度を算出する。また、図15でも波長λの単色光についての分光放射輝度算出処理を説明するが、実際には同様の処理を所定の波長毎に行うことで、当該対象画素についての分光放射輝度が得られる。図15において、発射光線数を決定するS1201以外の各処理は、第2実施形態で図13のフローチャートに示した各処理と同様の処理を行う。なお、第3実施形態においても、第2実施形態と同様にS1406で初期設定された最大交差回数Mを用いるが、この最大交差回数Mは最初に交差するオブジェクトに関わらす固定とし、更新されない。以下、第3実施形態のS1201における、分光特性に応じた光線数の決定処理について説明する。第2実施形態でも説明したように、オブジェクトの反射率や透過率が小さい波長の光線では、該オブジェクトとの交差によって放射輝度が大きく減衰するため、その後の他オブジェクトとの交差による放射輝度の変化量は小さい。そのため、交点から改めて飛ばす光線数が少ない場合に得られる放射輝度と、十分多い光線数によって得られる放射輝度とで、差がほとんど生じない。一方、オブジェクトの反射率や透過率が大きい波長の光線では、該オブジェクトとの交差後も放射輝度があまり減衰しないため、その後の他オブジェクトとの交差による放射輝度の変化量は相対的に大きくなる。そのため、少ない発射光線数によって得られる放射輝度と、十分多い光線数によって得られる放射輝度とで、大きな差が生じ、すなわちノイズが発生してしまう。そこで第3実施形態では、ある波長の光線が交差するオブジェクトごとに、その反射率または透過率が大きければ発射光線数を大きく設定し、逆に小さければ発射光線数を小さく設定する。なお、発射光線数の設定時に光源の放射輝度にも参照されるべきであることは、第2実施形態と同様である。すなわち、オブジェクトとの交差後に残存する光線の強度が大きいほど、発射光線数を大きく設定する。このように、光源の放射輝度と、オブジェクトの分光特性に応じて光線追跡時の光線数を設定することで、演算量の削減とノイズ発生の抑制を両立することが可能となる。ここで図16を用いて、特に交点からの反射成分についての発射光線数の設定方法の概念について説明する。ある入射角度、出射角度における、光線と交差するオブジェクトの分光反射率が図16で示される場合、例えば480nmのように反射率が高い波長では光線数を多く設定する。一方、例えば680nmのように反射率が低い波長では光線数を少なく設定する。具体的には、オブジェクトに対する反射光、透過光のそれぞれの光線数Rr,Rtは、上記式で表される分光特性frft、および光源の分光放射輝度Sの関数gとして、以下のように記述される。Rr=hS...Rt=hS...尚、関数hとしては、反射率が高い波長において光線数が大きくなるような関数であれば良く、具体的には、以下の式や式に示すような関数を用いることができる。hS=Rmax×ft×S...hS=Rmax×√×S...なお、式においてRmaxは任意の値を持つ定数であり、反射率や透過率が1である場合の光線数である。Rmaxが大きいほど、ノイズの少ない画像が得られる半面、計算量が増大する。第3実施形態では、以上のようにオブジェクトの分光特性に応じて設定した光線数で、光線とオブジェクトとの交点から新たな光線を飛ばし、光線経路の追跡を繰り返し行うことで、当該対象画素における放射輝度を算出する。以上説明したように第3実施形態によれば、分光による光線追跡を行う際に、光源の放射輝度と、オブジェクトの分光特性を考慮して発射光線数を設定する。これにより、波長ごとに必要最低限の光線数による演算が可能となるため、計算時間の高速化とノイズ発生の低減を実現することが可能になる。なお、第3実施形態を上述した第2実施形態と同時に適用することも可能である。すなわち、光源とオブジェクトの分光特性に応じて、光線の波長ごとに最大交差回数を設定し、さらに、各交点から飛ばす光線数をオブジェクトごとに設定すれば良い。<変形例>以下、第2および第3実施形態の変形例について説明する。上述した第2および第3実施形態では、光源特性やオブジェクトの反射透過特性には光の波長依存性がある、すなわち分光情報を持つものとして説明を行った。しかしながら本発明は、分光情報を利用した光線追跡処理への適用に限られるものではなく、例えば濃淡情報のみを持つ1色での光線追跡や、RGBの3色での光線追跡にも適用可能である。その場合、光線の放射輝度に対するオブジェクトの反射特性または透過特性が低い場合は、光線数または最大交差数を小さく設定することで、処理の高速化が実現できる。他方、光線の放射輝度に対するオブジェクトの反射特性または透過特性が高い場合は、光線数または最大交差数を大きく設定することで、処理の高速化よりもノイズ低減化を優先することができる。<その他の実施形態>また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェアプログラムを、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータまたはCPUやMPU等がプログラムを読み出して実行する処理である。"