任意のシーンについてのコンピュータグラフィックス画像を、光線追跡法によって生成する画像処理装置であって、デバイスの分光特性を示すデバイス特性データを入力するデバイス特性データ入力手段と、前記シーン内におけるオブジェクトの情報を示すオブジェクトデータを入力するオブジェクトデータ入力手段と、前記デバイス特性データに基づき、前記分光特性を特徴づける波長を追跡対象波長として決定する波長決定手段と、前記オブジェクトデータに基づいて前記追跡対象波長の光線に対する光線追跡を行うことによって、前記シーンにおける画素ごとの分光放射輝度を算出する光線追跡手段と、前記画素ごとの分光放射輝度を色信号に変換して前記コンピュータグラフィックス画像を生成する色変換手段と、を有することを特徴とする画像処理装置。前記デバイス特性データは、前記デバイスが備えるカラーフィルタの分光特性を示し、前記波長決定手段は、前記カラーフィルタの分光特性において、フィルタ透過率がピークとなる波長を中心とした所定範囲の第1の波長域から、前記追跡対象波長を抽出することを特徴とする請求項1に記載の画像処理装置。前記波長決定手段は、前記カラーフィルタの分光特性において、フィルタ透過率の変化が所定の閾値以上である第1の波長域から、前記追跡対象波長を抽出することを特徴とする請求項1に記載の画像処理装置。前記波長決定手段は、前記第1の波長域において、前記追跡対象波長を等間隔で抽出することを特徴とする請求項2または3に記載の画像処理装置。前記波長決定手段は、該第1の波長域とは異なる第2の波長域において、該第1の波長域における前記追跡対象波長の抽出間隔よりも大きな間隔で、前記追跡対象波長を抽出することを特徴とする請求項4に記載の画像処理装置。前記デバイスは、前記シーンと合成する実写映像を撮影する撮像装置であることを特徴とする請求項1乃至5のいずれか1項に記載の画像処理装置。前記デバイスは、前記コンピュータグラフィックス画像を表示するディスプレイであることを特徴とする請求項1乃至5のいずれか1項に記載の画像処理装置。前記オブジェクトデータは、前記シーン内におけるオブジェクトの位置および大きさ、分光反射特性および分光透過特性の情報を含むことを特徴とする請求項1乃至7のいずれか1項に記載の画像処理装置。デバイス特性データ入力手段、オブジェクトデータ入力手段、波長決定手段、光線追跡手段、および色変換手段、を有し、任意のシーンについてのコンピュータグラフィックス画像を、光線追跡法によって生成する画像処理装置における画像処理方法であって、前記デバイス特性データ入力手段が、デバイスの分光特性を示すデバイス特性データを入力し、前記オブジェクトデータ入力手段が、前記シーン内におけるオブジェクトの情報を示すオブジェクトデータを入力し、前記波長決定手段が、前記デバイス特性データに基づき、前記分光特性を特徴づける波長を追跡対象波長として決定し、前記光線追跡手段が、前記オブジェクトデータに基づいて前記追跡対象波長の光線に対する光線追跡を行うことによって、前記シーンにおける画素ごとの分光放射輝度を算出し、前記色変換手段が、前記画素ごとの分光放射輝度を色信号に変換して前記コンピュータグラフィックス画像を生成することを特徴とする画像処理方法。コンピュータ装置で実行されることにより、該コンピュータ装置を請求項1乃至8のいずれか1項に記載の画像処理装置の各手段として機能させるためのプログラム。実物に忠実なシーンをCGで再現し、実写画像と合成するためには、該シーンが実物に忠実なことだけでなく、実写画像を撮影したデバイスの特性も加味した上で、画像を合成することが重要となる。しかしながら上記特許文献1に記載の方法では、デバイスの特性については特に考慮されないため、実写画像とCGを合成して一つのシーンとして適切に構築するには、さらなる調整が必要となる。また、特許文献1に記載の方法も含め、一般的な光線追跡法では、光を波長に依存しない一様な光線として扱う。したがって、例えばオブジェクトを透過する光の屈折のように、波長によって異なる振る舞いをする物理現象を正しく再現することができない、という問題があった。本発明は上述した問題を解決するために、デバイスの特性を考慮した光線追跡を行うことで、シーンを忠実に再現したコンピュータグラフィックス画像を効率良く生成することを目的とする。上記の課題を解決するための一手段として、本発明の画像処理装置は以下の構成を備える。すなわち、任意のシーンについてのコンピュータグラフィックス画像を、光線追跡法によって生成する画像処理装置であって、デバイスの分光特性を示すデバイス特性データを入力するデバイス特性データ入力手段と、前記シーン内におけるオブジェクトの情報を示すオブジェクトデータを入力するオブジェクトデータ入力手段と、前記デバイス特性データに基づき、前記分光特性を特徴づける波長を追跡対象波長として決定する波長決定手段と、前記オブジェクトデータに基づいて前記追跡対象波長の光線に対する光線追跡を行うことによって、前記シーンにおける画素ごとの分光放射輝度を算出する光線追跡手段と、前記画素ごとの分光放射輝度を色信号に変換して前記コンピュータグラフィックス画像を生成する色変換手段と、を有することを特徴とする。本発明によれば、デバイスの特性を考慮した光線追跡を行うことで、シーンを忠実に再現したコンピュータグラフィックス画像を効率良く生成することができる。以下、本発明に係る実施形態について、図面を参照して説明する。なお、以下の実施の形態は特許請求の範囲に関る本発明を限定するものではなく、また、本実施の形態で説明されている特徴の組み合わせの全てが本発明の解決手段に必須のものとは限らない。<第1実施形態>本実施形態においては、デジタルカメラで撮影された実写映像との合成対象となるシーンのコンピュータグラフィックス画像を、光線追跡によって生成する処理について説明する。このとき、該デジタルカメラにおけるカラーフィルタの分光特性に基づいて、光線追跡を行う波長を決定することを特徴とする。まず、波長によって異なる光の振る舞いについて、図1を用いて説明する。光はオブジェクトにあたると、その表面で反射する光と透過する光に分離する。オブジェクトを透過する光は、その波長λに応じた分光屈折率によって屈折角が変化する。例えば図1に示すように、視点を開始点として、分光屈折率n1の空間を通過した光線が分光屈折率n2のオブジェクト101を透過する場合を考える。このとき、オブジェクトに対する光線の入射角をφとすると、屈折角θは以下の式で示される。θ=sin1n2sinφこのように屈折角は光の波長λの関数になるため、光線追跡法によって追跡する光の経路が波長λによって変化することになる。従って、図1でオブジェクト101を透過した光がさらに分光屈折率n3のオブジェクト102を透過する等、あるシーンにおける光の透過を忠実に再現するためには、光を分光として扱い、複数の光線から色を再現する必要がある。また、正確に色を再現するためには、光源から発せられる光の全ての波長に対して追跡計算を行う必要がある。しかし、全ての波長に対して光線追跡計算を行うことは、計算量が余りにも膨大となってしまうため現実的には不可能である。そこで、全波長域から所定間隔おきに抽出した波長について光線追跡を行って、色再現を行うことが一般的である。●装置構成本実施形態における画像処理装置の構成例について、図2を用いて説明する。図2に示すように本実施形態の画像処理装置は、CPU201、RAM202、外部記憶装置203、汎用インターフェイス204、ROM205、モニタ206、メインバス207を備える。IF204は、不図示の光ディスク装置やネットワーク装置等を、メインバス207に接続する。以下、CPU201がROM205に格納された各種ソフトウェアを動作させることで実現する各種処理について説明する。まずCPU201は、ROM205に格納されている画像処理アプリケーションを起動し、RAM202に展開する。続いて、外部記憶装置203に格納されている各種データや、IF204に接続された不図示の分光測定器や外部記憶装置から入力されるデータ等が、CPU201からの指令に基づき、メインバス207経由でRAM202に転送される。さらに、画像処理アプリケーションにしたがって、RAM202に格納されているデータにCPU201からの指令に基づく各種演算が施され、該演算結果をメインバス207経由でモニタ206に表示したり、外部記憶装置203に格納したりする。以下、上記構成において、CPU201からの指令に基づき画像処理アプリケーションが画像を生成してモニタ206に表示するまでの画像処理について、詳細に説明する。図3は、本実施形態における画像処理装置301の論理構成を示すブロック図であるが、同図に示す各構成は、上述したように画像処理アプリケーションとして実現される。図3において、デバイス特性データ入力部302は、対象デバイス、すなわちシーンへの合成対象となる実写映像を撮影した撮像デバイスの分光特性データを、デバイス特性データ309として入力する。オブジェクトデータ入力部303は、シーンを構成する各オブジェクトの情報を、オブジェクトデータ308として入力する。尚、デバイス特性データ入力部302およびオブジェクトデータ入力部303で入力されるデータは、ROM205や外部記憶装置203、IF204を介して接続される他の手段によって入力されても良い。分光放射輝度算出部304は、デバイス特性データ入力部302及びオブジェクトデータ入力部303から入力された各種データに基づき、CPU201からの指示に従って、シーンのCG画像における各画素の分光放射輝度307を算出する。色変換部305は、分光放射輝度算出部304で算出された分光放射輝度307に基づき、CPU201からの指示に従ってシーンのCG画像における各画素の色信号値を算出し、画像データ306として出力する。画像データ306は、外部記憶装置203等の記憶領域に格納されるとともに、モニタ206に表示される。●デバイス特性データここで、本実施形態におけるデバイス特性データ309の詳細について、図5を用いて説明する。図5は、対象デバイスであるデジタルカメラにおけるRGBの各カラーフィルタの分光透過率を示すデータであり、この波長域において抽出された特定の波長が、デバイス特性データ309として予めテーブル形式で保持されている。図5の横軸は波長nm、縦軸は分光透過率を示し、RGBそれぞれの波長帯の光を透過するフィルタが存在することを示している。本実施形態では、光線追跡法を適用する光線の波長追跡対象波長を、このデバイス特性データとして示された分光特性に基づいて決定する。具体的には図5に点線で示すように、RGBそれぞれの透過波長域において、フィルタ透過率が最も高いピーク部分から所定範囲で、追跡対象波長をより細かく抽出する。したがって相対的に、フィルタ感度が低い範囲では追跡対象波長の抽出数が少なくなる。ここで従来の光源追跡法によれば、追跡対象波長は全波長域から等間隔で抽出しており、例えば380~730nmの可視域において、10nm間隔で36波長を追跡対象波長として抽出し、分光計算を行っていた。対して本実施形態では、例えばRフィルタの透過ピークである600~730nmの第1の波長域については5nm間隔で追跡対象波長を抽出し、他の第2の波長域では30nm間隔で抽出する。このように、分光特性に応じて追跡対象波長の抽出間隔を変更することによって、従来の等間隔による抽出を行う場合と比べて、追跡対象波長の全体数は減らしながらも、フィルタのピーク近傍については逆に増やすことができる。したがって、全体としての演算量を削減しつつ、よりデバイス特性を反映したレンダリングを行うことが可能となる。本実施形態では、光線追跡に用いる追跡対象波長を以上のように予め決定し、デバイス特性データ309としてそれぞれの透過率とともにルックアップテーブル形式で保持している。実際には、デバイス特性データ309は例えば外部記憶装置203に保持されている。なお、デバイス特性データ309としては、複数種類のデバイスに対応したものがそれぞれ保持されており、後述するように、シーンのCGに反映させるべき対象デバイスについてのデバイス特性データ309が選択的に参照される。以降で説明する光線追跡計算では、デバイス特性データ309として保持された波長ごとに、演算を行う。●シーンCG画像の生成処理以下、本実施形態の画像処理装置301において、実写映像との合成対象となるシーンのCG画像を作成する画像生成処理について、図4のフローチャートを用いて説明する。まずS401で、シーンのCGがその特性を再現すべき対象デバイスを、例えば不図示のUIより指定する。ここで指定する対象デバイスとは、該シーンに合成される実写映像を撮影するデジタルカメラの機種名を指し、該機種名に応じたデバイス特性データ309が予め保持されているものとする。なお、デバイス特性データ309が特定できれば、機種以外の名称、たとえばカラーフィルタ名称等を指定しても良い。次にS402では、デバイス特性データ入力部302から、S401で指定された対象デバイスのデバイス特性データ309を読み込み、RAM202に格納する。そしてS403で、オブジェクトデータ入力部303からオブジェクトデータ308を入力する。ここでオブジェクトデータ308とは、CGとして生成する当該シーンに含まれるオブジェクトの位置や大きさ、反射特性等を示すデータである。次にS404で、当該シーンのCGとしての出力画像を順次走査して、分光放射輝度の算出対象となる画素を決定する。なお、この対象画素の選択順は任意である。以下、S405~S407で分光放射輝度算出部304において、上記S404で選択された対象画素についての分光放射輝度を算出する。まずS405で、上記S402で読み込まれたデバイス特性データ309に保持された波長群から、光線追跡の対象となるひとつの波長を追跡対象波長として設定する。なお、この追跡対象波長の選択順は任意である。そしてS406で、対象画素における追跡対象波長の放射輝度を算出する。この放射輝度の算出の際にオブジェクトデータ308を用いた光線追跡を行うが、その詳細については後述する。そしてS407では、デバイス特性データ309に保持された全波長について、放射輝度の算出が終了したか否か、すなわち、当該対象画素についての分光放射輝度の算出が終了したか否かを判定する。未終了であればS405に戻って次の追跡対象波長を決定し、処理を継続する。一方、全波長についての処理が終了していればS408に進み、シーンの出力画像の全画素について、分光放射輝度の算出が終了したか否かを判定する。終了していなければS404に戻って次の対象画素を決定し、処理を継続する。全画素についての処理が終了していればS407に進み、色変換部305において、出力画像の各画素について分光放射輝度から色信号値を算出することで、画像データ306が生成される。この画像データ306がすなわち、実写合成の対象となるシーンのCG画像である。なお、色信号値の算出方法の詳細については後述する。●放射輝度算出処理以下、上記S406における追跡対象波長の放射輝度の算出処理について、図6を用いて詳細に説明する。上述したように、追跡対象波長の放射輝度の算出には光線追跡法を用いる。この手法によれば、本実施形態であればデジタルカメラに相当する視点が画角内の各方向から受け取るはずの光線を、仮想的に光源まで追跡することでCG画像の描画を行う。詳細には、視点から描画する各画素の方向へ直線を伸ばし、オブジェクトとの交差の有無を数学的に判定する。そして、反射や透過により発生した光線について、再帰的に探索を繰り返し、光源と交差した場合に計算を終了する。すなわち本実施形態では、図6のフローチャートに示す処理を再帰的に実行することで複数の光線を追跡し、光線ごとに算出した放射輝度の平均により、対象画素における追跡対象波長の放射輝度を算出する。なお、図6では波長λの単色光についての放射輝度算出処理を説明するが、実際には図4のS405~S407のループ処理で示されるように、同様の処理をデバイス特性データ309として入力された波長毎に行う。これにより、デバイス特性データ309として予め設定された各追跡対象波長についての放射輝度が算出され、当該対象画素についての分光放射輝度が得られる。まずS601で、追跡の開始点から飛ばした波長λの光線がシーンにおいて交差する点を探索する。ここで追跡の開始点は、当該フローチャートが初回の実行であれば対象画素であり、再帰実行であれば光線とオブジェクトとの交点である。なお、この探索はオブジェクトデータ308に基づいて行われるが、探索方法は任意であり、例えば、シーン内の全てのオブジェクトについて順次交差判定を行って、追跡の開始点に最も近い交点を求めれば良い。次にS602で、S601で探索された光線と交差する点が、光源との交差であるか、またはその他のオブジェクトとの交差であるかを判定する。光線が光源と交差している場合にはS603に移行し、光線の強度と光源の光源特性データから、追跡対象波長λについての放射輝度を算出して、当該波長についての処理を終了する。ここで光源特性データとは、光源の配光特性や分光放射輝度等であり、RAM202等に予め保持されているとする。詳細には、光線の強度をP0、光線と光源のなす角をθ、光源からθ方向の配光特性をQ、光源の分光放射輝度をSとすると、光線の分光放射輝度E0は以下の式によって示される。E0=P0QS...一方、光線が光源以外のオブジェクトと交差している場合にはS604に移行し、まず、光線が交差したオブジェクトの分光特性を、オブジェクトデータ308から取得する。次にS605で、光線とオブジェクトとの交点から飛ばす光線数を決定する。この発射光線数は任意に決定して構わないが、数が多いほどノイズの少ない画像が得られる。以下、発射光線数に応じた回数分、S606及びS607の処理をループする。S606では、発射光線の方向および強度を算出する。光線の方向は発射光線数に応じて均等に分布すれば良く、該交点に対する反射光の強度Pr、透過光の強度Ptはそれぞれ、以下の式によって示される。なお、下式において、P0はオブジェクトと交差する前の光線の強度、frはオブジェクトの反射特性、ftはオブジェクトの透過特性である。Pr=P0fr...Pt=P0ft...ここで、オブジェクトの反射特性および透過特性について説明する。一般に、オブジェクトの反射特性は双方向反射率分布関数frで表される。オブジェクト上の位置xにおいて、ある方向ωiから入射する光の波長をλi、強度をLiとし、ある方向ωrに反射する光の波長をλr、強度をdLrとすると、反射特性frは以下の式で示される。fr=dLr/Lidω...一般に、反射光の強度分布はオブジェクトの表面形状に依存し、オブジェクト表面が平滑なほど、鏡面反射方向すなわちωi=ωr方向の反射強度が高くなり、それ以外の拡散反射方向への反射強度は小さくなる。尚、オブジェクトの反射特性frの取得方法としては、オブジェクトデータ308以外の、外部記憶装置203等に記憶されているデータを入力しても良いし、IF204に接続された分光測定器で測定したデータを直接入力しても良い。また、反射特性frのデータ形式としては、オブジェクトへの光の入射角と反射角に対応した反射率が得られれば、任意であって良い。例えば、上記式のように一意の関数として入力しても良いし、一般に拡散反射は等方的に分布するため、ωi=45°ωr=0°の反射特性を拡散反射に適用し、鏡面反射は別途LUTにする等、分割して入力しても良い。また、オブジェクトの透過特性は双方向透過率分布関数ftで表される。オブジェクト上の位置xにおいて、ある方向ωiから入射する光の波長をλi、強度をLiとし、ある方向ωtに透過する光の波長をλt、強度をdLtとすると、透過特性ftは以下の式で示される。ft=dLt/Lidω...一般に、透過光の強度分布は分光屈折率に依存するため、ωtが上記式を満たす場合に、透過強度が最大になる。尚、オブジェクトの透過特性ftの取得方法としても、オブジェクトデータ308以外の、例えば外部記憶装置203等に記憶されているデータを入力しても良い。また、透過特性ftのデータ形式も、オブジェクトへの光の入射角と屈折角に対応した透過率が得られれば、任意であって良い。例えば、上記式のように一意の関数として入力しても良いし、ωi=ωr=0°の透過特性をωtが式の屈折角θに等しい場合に適用し、それ以外の透過は別途LUTにする等、分割して入力しても良い。以上のようにS606で発射光線の方向と反射光および透過光の強度PrPtが算出されると、次にS607で図6のフローチャートに示す処理を再帰的に呼び出すことによって、当該発射光線に対する光線追跡処理が行われる。このS607の処理はすなわち、再帰処理によって最終的に交差するオブジェクトがなくなり、当該発射光線についての放射輝度が光源特性データに基づいてS603で算出された時点で、終了する。このように、発射光線についてS607の再帰処理によって最終的に上記式で算出された放射輝度E0が、当該発射光線についての放射輝度Enとなる。そして、S608で全ての発射光線について上記S606及びS607による光源追跡が終了したと判定されると、次にS609において、各発射光線の放射輝度Enの平均を算出する。この平均放射輝度値が、対象画素における追跡対象波長λについての放射輝度Exである。すなわち、光線の総数をN、各発射光線の分光放射輝度をEnとすると、対象画素における追跡対象波長λの放射輝度Exは、以下の式で算出される。なお、式においてΣはn0からN1までの整数に関する総和を示す。Ex=ΣEn...以上、S406では図6に示す処理を再帰的に実行することによって、対象画素における追跡対象波長λの単色光線を追跡し、その放射輝度Exを算出する。本実施形態では、上述したようにこの放射輝度算出処理をデバイス特性データ309が保持する各波長について行うことで、波長λごとの分光放射輝度Exを取得する。●色信号値算出処理以下、上記S407における画素ごとの色信号値の算出処理について、図7を用いて詳細に説明する。尚、本実施形態では最終的にsRGB形式の色信号値を出力するとして説明するが、AdobeRGB等の他形式による色信号値であっても構わない。まずS701で、視覚特性として、例えばCIEXYZ表色系の等色関数xyzを取得する。次にS702で、上記S405で算出した、画素ごとの分光放射輝度Exを取得する。そしてS703で、取得した視覚特性xyzと分光放射輝度Exから、以下の~式によりCIEXYZ値を算出する。そしてS704において、出力画像の色信号値を算出する。ここでは出力画像をsRGB形式とするるため、出力信号RoutGoutBoutを以下の~式によって算出する。以上、図4に示す画像生成処理によって、合成対象シーンのCGである出力画像が形成される。以上説明したように本実施形態によれば、シーンにおけるオブジェクトデータを利用した光線追跡を、デバイスの分光特性に応じて設定された光線波長について行う。これにより、一般的な等間隔で抽出された波長を用いた光線追跡を行う場合に比べて、より少ない波長数による演算で、シーン内におけるオブジェクトの色を、あたかも該デバイスで撮影したように再現することができる。したがって、該シーンを、該デバイスで撮影された実写映像との合成に適したCGとして生成することができる。。なお本実施形態では、対象デバイスとしてデジタルカメラを例として説明したが、もちろん他のデバイスであっても構わない。例えば、対象デバイスをディスプレイとした場合、デバイス特性データはディスプレイのカラーフィルタ、およびバックライトの特性となる。このとき、本実施形態に基づいてシーン画像の色決定を行うと、ディスプレイ上で最も広色域となるようにレンダリングさせることができる。<第2実施形態>以下、本発明に係る第2実施形態について、特に第1実施形態と異なる個所について説明する。第2実施形態は、上述した第1実施形態に対し、図3に示したデバイス特性データ入力部302の動作、および、デバイス特性データ309の内容が異なる。図8に、第2実施形態における画像処理装置の構成例を示すが、デバイス特性データ生成部901を除き、他の構成については上述した図3と同様であるため、説明を省略する。デバイス特性データ生成部901は、予め保持された詳細なデバイス特性データ309に基づき、分光放射輝度算出部304に入力すべきデバイス特性データを生成する。以下、第2実施形態におけるデバイス特性データの生成方法について、図9のフローチャートを用いて説明する。まず、第2実施形態におけるデバイス特性データ309としては、上記図5に示すようなフィルタの分光透過率を示す詳細なデータを有している。例えば可視光域となる380~730nmにおいて01nmおきのフィルタの分光透過率を有しているものとする。第2実施形態ではこのような詳細なデバイス特性データ309に基づき、光線追跡法の演算に適用する追跡対象波長を決定する。なお、第2実施形態におけるデバイス特性データ309としては、上記可視光域内の01nm刻みの例に限らず、もちろんその他の波長域刻み幅であっても構わない。デバイス特性データ生成部901はまずS1001において、デバイス特性データ309を読み込んでRAM202へ格納する。そしてS1002で、該読み込んだデバイス特性データ309の中から、透過率がピークを示す波長を求める。このピーク検出には、SavitzkyGolayアルゴリズム等の一般的なピーク検出アルゴリズムを用いることができるが、他の方法を用いても構わない。次にS1003では、S1002で得られたピーク波長に基づき、例えばピーク波長を中心に1nm刻みで前後5波長分を、追跡対象波長として決定する。なお、追跡対象となる波長決定の方法としては上述した方法に限定されず、デバイス特性における特徴的な部分を示す波長を決定できれば、他の方法を適用することももちろん可能である。まず、追跡対象波長の抽出範囲をフィルタ透過率が0でない範囲に限定することも有効である。また、ピーク波長周辺以外にも、透過率の変化の大きい波長域から追跡対象波長を抽出することも有効である。例えば、フィルタ透過率が0でない波長範囲を所定の波長間隔で参照し、隣り合う参照波長の透過率変化が所定の閾値以上であれば、両波長とも追跡対象波長に設定すれば良い。この場合、隣り合う参照波長の透過率変化が01以下である部分については、該参照波長の間隔よりも大きい間隔で、追跡対象波長を抽出する。以上説明したように第2実施形態によれば、詳細なデバイス特性データに基づき、光線追跡に適用する波長を柔軟に決定可能とすることで、実写合成の対象となるシーンのCGをより柔軟に生成することが可能となる。<その他の実施形態>また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェアプログラムを、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータまたはCPUやMPU等がプログラムを読み出して実行する処理である。