JP2005505810A - サンプル点生成のためのリカーシブローテーション(recursiverotations)を用いた決定論的手法(strictlydeterministicmethodologies)を用いた画像レンダリング - Google Patents

サンプル点生成のためのリカーシブローテーション(recursiverotations)を用いた決定論的手法(strictlydeterministicmethodologies)を用いた画像レンダリング Download PDF

Info

Publication number
JP2005505810A
JP2005505810A JP2003502789A JP2003502789A JP2005505810A JP 2005505810 A JP2005505810 A JP 2005505810A JP 2003502789 A JP2003502789 A JP 2003502789A JP 2003502789 A JP2003502789 A JP 2003502789A JP 2005505810 A JP2005505810 A JP 2005505810A
Authority
JP
Japan
Prior art keywords
sample points
sub
ray
sample
rays
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2003502789A
Other languages
English (en)
Other versions
JP4133806B2 (ja
Inventor
アブラモブ,ジョージー
ジョンソン,クリスチャン−ヴァルール
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nvidia ARC GmbH
Original Assignee
Mental Images GmbH
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mental Images GmbH filed Critical Mental Images GmbH
Publication of JP2005505810A publication Critical patent/JP2005505810A/ja
Application granted granted Critical
Publication of JP4133806B2 publication Critical patent/JP4133806B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/50Lighting effects
    • G06T15/506Illumination models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T1/00General purpose image data processing
    • G06T1/20Processor architectures; Processor configuration, e.g. pipelining
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/50Lighting effects
    • G06T15/55Radiosity

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Image Generation (AREA)
  • Complex Calculations (AREA)
  • Image Processing (AREA)
  • Electron Beam Exposure (AREA)
  • Color Television Systems (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Image Analysis (AREA)
  • Electrophonic Musical Instruments (AREA)
  • Control Of Indicators Other Than Cathode Ray Tubes (AREA)
  • Error Detection And Correction (AREA)

Abstract

コンピュータグラフィックスシステムは、画像の中の画素の画素値を生成し、画素はシーンの中の点を表わす。コンピュータグラフィックスシステムは、画素値を選択された関数の積分の評価により生成する。コンピュータグラフィックスシステムは、サンプル点生成手段と関数評価手段とを備える。サンプル点生成手段は複数の追跡レベルを持つように構成されたレイの追跡において各々一連のレイの1つに関連付けられたサンプル点の各集合を生成するように構成される。少なくとも1つのレベルにあるレイが複数のレイに分割でき、各レイはレイインスタンス識別子に関連付けられる。サンプル点生成手段は各レベルに対して選択されたローテーション演算子が再帰的に適用される所定の決定論的な低食い違い量列としてサンプル点を生成するように構成される。関数評価手段は、サンプル点生成手段により生成されたサンプル点の1つにおける選択された関数の評価を各々代表する複数の関数値を生成するように構成され、画素値を生成するのに関数値を用いる。一実施形態では、選択されたローテーション演算子は、クランレイ・パターソンローテーション演算子である。

Description

【技術分野】
【0001】
(参照としての組み入れられるもの)
1997年6月23日にマーチン・グラベンスタイン等(Martin Grabenstein, et al.)の名前で出願され、「サンプル点を生成するために決定論的手法を用いた画像における画素の画素値を生成するシステム及び方法、」("System And Method For Generating Pixel Values For Pixels In An Image Using Strictly Deterministic Methodologies For Generating Sample Point, ")と題される本出願の譲受人に譲渡された米国特許出願第08/880,418号(以下グラベンスタインの出願として参照)は、参照として組み入れられる。
【0002】
2001年6月19日にアレキサンダー・ケラー(Alexander Keller)の名前で出願され、「サンプル点を生成するために決定論的手法を用いた画像における画素の画素値を生成するためのシステム及び方法」("System And Method For Generating Pixel Values For Pixels In An Image Using Strictly Deterministic Methodologies For Generating Sample Points,"と題される本出願の譲受人に譲渡された米国特許出願第09/884,861号(以下ケラーの出願として参照)は、参照として組み入れられる。
【0003】
(技術分野)
本発明は、広くはコンピュータグラフィックスに関し、特に、画素値を表す積分値を概算するためのサンプル点を提供するために、低食い違い量列(超一様分布列:low-discrepancy sequences)及びクランレイ・パターソンローテーション(Cranley-Patterson rotations)などのリカーシブローテーション(recursive rotations)を用いた決定論的手法を用いてレンダリングされた画像におけるピクセルの画素値を生成するシステムと方法に関する。
【背景技術】
【0004】
コンピュータグラフィックスにおいて、コンピュータは、例えばカメラ等によるシーンの記録をシミュレートするために、例えば1つ以上の光源に照明される3次元シーンにおけるオブジェクトの面の2次元画像平面への投影を表すデジタルデータを生成するために用いられる。カメラは、シーンの画像を画像平面に投影するためのレンズを備えるかもしれないが、レンズが用いられない場合には、ピンホールカメラであるかもしれない。2次元画像は、画素要素(画素(ピクセル)またはぺルと省略される)の配列から構成され、各ピクセルに対して生成されるデジタルデータは、画像平面における各々のピクセル位置での画像平面に投影されたシーンの色と輝度を表す。オブジェクトの面は、リアルに見える画像を得るために、形状、色、スペキュラー、テクスチャーなどを含めた可能な限り本物に近い形で画像に表現されるのが好ましい、いくつかのタイプの特性を含んでいる。
【0005】
一般にコンピュータグラフィックスにおいて、画像は、光源と画像平面のピクセルとの間において追跡されるレイ(rays)を評価(evaluate)するため積分を評価(evaluating)することによりレンダリングされ、レイは各々の光源と画像平面上のピクセルとの間において光子(photons)によって取られ得る径路を表す。積分の評価は、本質的には、画像平面の各ピクセルに入射する光子の流れ(光束:photon flux)の強さと色に関する情報を与える。一般に、積分は以下の形式である。
【数1】
Figure 2005505810
ここで“f”はs1+・・・+sn=s次元単位立方体[0,1)sにおいて定義される関数である。光源に対する画像平面の位置と方向によって、いくつかのレイは光源から画像平面上のピクセルへ、直接光子の軌跡(photon path)を描くかもしれない。同様に、シーンの中のオブジェクトに対する画像平面の位置と方向によって、いくつかのレイは光源からシーンの中のオブジェクトの1以上の面において反射された後、画像平面上のピクセルに達する光子の軌跡を描くかもしれない。上述したように、オブジェクトの面は、いくつものタイプの特性を持つことができるので、それらから反射される光子の色や方向などに影響を及ぼすことができる。面の色は反射された光子の色に影響を及ぼすことができる。例えば、もし面が青色であれば、反射された光子は主に(完全でないならば)青色の光子から構成される。面の他の特性はそこから反射される光子の方向に影響を及ぼすことが可能であり、したがって各々の面から追跡されるレイの軌跡に影響を及ぼし得る。様々なタイプの面の特性が光子を幾つかの方向に反射させることから、面に入射されるレイは一般に、各面の特性に依存する方向への複数の反射レイに「分割(split)」される。例えば、面が全体として鏡面的(specular)であれば、反射レイは入射角に近い似通った方向となる。一方、面が全体として拡散的(diffuse)であると、反射レイは全体として入射角とはかなり異なる多様な方向となる。
【0006】
一般的に、式(1)によって表わされる積分は、閉じた形の解(closed form solutions)をもたず、通常数値的に評価される。式(1)によって表わされる積分は次のように数値的に見積もることができる。
【数2】
Figure 2005505810
ここで、総和の中の各f(xi1 (1),・・・,xin (n))は、積分領域内のサンプル点(xi1 (n),・・・,xin (n))における関数“f”の評価値(evaluation)である。サンプル点は、幾つもの手法を用いて生成することができる;グラベンスタインの出願(Grabenstein application)は、サンプル点がハルトン列(Halton sequences)などのs次元の決定論的低食い違い量列からのサンプルを含む手法について記述し、ケラーの出願(Keller application)は、サンプル点がs次元ハマースレイ点群(Hammersley point sets)を含む手法について記述する。一般に、それらの次元s1,・・・snは、それぞれ光源と画像平面との間の光子の軌跡に関連付けられている。
【0007】
問題は、特に幾つものレイが実質的に同じ方向をもつ状況において発生する。
【発明の開示】
【発明が解決しようとする課題】
【0008】
本発明は、画素値を代表する積分値を概算するためのサンプル点を与えるために決定論的低食い違い量列とクランレイ・パターソンローテーションなどのリカーシブローテーションを用いてレンダリングされる画像のピクセルに対する画素値を生成する新規で改良されたシステムと方法を提供する。
【0009】
概略、本発明は、画像のピクセルに対する画素値を生成するコンピュータグラフィックスシステムを提供し、ピクセルはシーンにおける1つの点を代表する。このコンピュータグラフィックスシステムは、画素値を選択された関数の積分評価により生成する。コンピュータグラフィックスシステムは、サンプル点生成手段(sample point generator)と関数評価手段(function evaluator)とを備える。サンプル点生成手段は、複数の追跡レベル(trace levels)をもつように構成された光線追跡(ray trace)において、一連のレイ(a series of rays)の1つに各々関連付けられるサンプル点の集合をそれぞれ生成するように構成される。少なくとも1つのレベルにあるレイは複数のレイに分割でき、各レイはレイインスタンス識別子(ray instance identifier)に関連付けられている。サンプル点生成手段は、選択されたローテーション演算子がそれぞれのレベルに応じて再帰的に(recursively)適用される所定の決定論的低食い違い量列としてサンプル点を生成するように構成される。関数評価手段は、複数の関数値を生成するように構成され、関数値の各々はサンプル点生成手段により生成されたサンプル点の1つにおける選択された関数の評価を表わす。また、関数評価手段は、画素値を生成するのにこの関数値を用いる。
【0010】
ある実施形態では、選択されたローテーション演算子は、クランレイ・パターソンローテーション演算子である。
【0011】
更に、ある実施形態では、サンプル点生成手段はそれぞれの追跡レベルに対して生成するサンプル点をキャッシュするように構成され、次の追跡レベルに対するサンプル点を生成する間、キャッシュされたサンプル点を用いる。
【0012】
本発明は、添付された請求項における特徴により指摘される。本発明の上記及び更なる長所は、添付された図面と連携する以下の説明を参照することにより、より良く理解されるであろう。
【発明を実施するための最良の形態】
【0013】
本発明は、シーンの画像の中のピクセルに対する画素値を生成するためのコンピュータグラフィックスシステムと方法を提供し、これはシーン内の様々な点から反射された光の各画素値への寄与を表わす関数の積分を評価するためのサンプル値を生成するために用いられるサンプル点の生成に決定論的手法を用い、選択されたローテーション演算子が適用される。ここに添付された図1には、このような決定論的手法を用いる実例としてのコンピュータシステム10が描かれる。図1を参照すると、一実施形態におけるコンピュータシステム10は、プロセッサモジュール11とオペレータインターフェース部を含み、オペレータインターフェース部は、キーボード12A及び/又はマウス12Bなどのオペレータ入力部(全体としては、オペレータ入力部12として扱われる)と、ビデオディスプレイ装置13などのオペレータ出力部を備える。実例としてのコンピュータシステム10は、従来のプログラム内蔵方式のコンピュータアーキテクチャーである。プロセッサモジュール11は、例えば、1個以上のプロセッサ、メモリとディスク及び/又はテープ記録部などの大容量記憶装置(個別には示されない)を含み、これらはそれに与えられるデジタルデータに関する処理及び記録作業を実行する。プロセッサモジュール11が、複数のプロセッサデバイス含むならば、それぞれのプロセッサデバイスはシングルタスクの様々な部分を処理するように構成されてもよく、この場合タスクはそれ以外の場合よりも速く実行され得る。オペレータ入力部12は、オペレータが処理のための情報を入力することを可能とするために設けられる。ビデオディスプレイ装置13は、オペレータにプロセッサモジュール11により生成された出力情報をスクリーン14に表示するために設けられ、これには処理の間に生成される情報の他、オペレータが処理のために入力するであろうデータや、処理を制御するためにオペレータが入力するであろう情報も含まれる。プロセッサモジュール11は、いわゆる「グラフィカルユーザインターフェース」(“GUI”)を用いてビデオディスプレイ装置13での表示のための情報を生成し、これにより様々なアプリケーションプログラムの情報が個々別々の「ウインドウ」を用いて表示される。コンピュータシステム10は、入力情報をオペレータから受け取るためのキーボード12Aやマウス12B、オペレータに出力情報を表示するためのビデオディスプレイ装置13など、特定のコンポーネントを備えたかたちで示されているが、コンピュータシステム10は、様々なコンポーネントを追加的に、あるいは図1に示されたものの代わりに備えていてもよい。
【0014】
更に、プロセッサモジュール11は、1以上のネットワークポートを備え(全体として参照番号14により参照される)、これはコンピュータネットワークにコンピュータシステム10を接続するコミュニケーションリンクに接続される。ネットワークポートは、コンピュータシステム10からネットワーク上の他のコンピュータシステムや他の装置へ情報を転送し、またはそれらから情報を受信することを可能とする。例えばクライアントサーバ方式(client-server paradigm)に従って構成された一般的なネットワークにおいては、ネットワークの中のあるコンピュータシステムがサーバーとして任命され、これらには他のクライアントコンピュータシステムでの処理のためのデータやプログラム(包括的には「情報」)が蓄えられており、これによりクライアントコンピュータシステムは情報を便利に共有することができる。特定のサーバーに保持された情報にアクセスする必要があるクライアントコンピュータシステムは、サーバーがネットワークを介してそれに情報をダウンロードすることを許可する。またデータを処理した後、クライアントコンピュータシステムは、処理されたデータを保存しておくためにサーバーに戻してもよい。これらのコンピュータシステム(上述したサーバーやクライアントを含む)に加えて、ネットワークは更に、プリンターやファクシミリ装置、デジタルオーディオやビデオ記録装置と分配装置などを含んでいてもよく、これらはネットワークに接続された様々なコンピュータシステムにおいて共有されてもよい。ネットワーク上のコンピュータシステムを双方向に連結するコミュニケーションリンクは、従来周知のように、利便性のあるいかなる情報伝達媒体(information-carrying medium)を含んでいてもよく、これにはコンピュータシステム間において信号を伝達する有線、光ファイバーや他の媒体が含まれる。コンピュータシステムは、コミュニケーションリンクを介して伝達されるメッセージを用いて情報をネットワークを介して転送し、各メッセージは情報とメッセージを受け取る装置を識別するための識別子(identifier)を有する。
【0015】
本発明により扱われる問題は、図2との関係において説明される。一般的にコンピュータグラフィックスでは、光源と画像平面上のピクセルとの間において追跡されるレイの特性を決定するために積分を評価することにより画像がレンダリングされる。レイは、それぞれの光源と画像平面上のピクセル間においてシミュレートされる光子が取る経路を表す。図2には、基本的に光沢があるが略鏡面的な2つの面21、22からの反射を含み、基本的に光沢があるが略拡散的する1つの面23からの反射に引き継がれるレイの軌跡20が1例として描かれる。説明を簡略にするために、図2では反射は実際には屈折により表わされている。第1の面21に入射する光子の流れを表わすレイ24は、反射において、面21の特性に依存する方向をもつ選択された数のレイに分割される。図2に描かれたイラストでは、面21に入射するレイ24は、2つのレイ25、26に分割されるが、コンピュータグラフィックスシステム10は、面からの反射においていかなる数のレイに入射レイを分割するように制御されてもよい。図2に描かれたイラストでは、レイ25、26の双方は、面22に向けられるが、面21が略鏡面的であることから入射角は僅かにしか異ならない。
【0016】
レイ24が面21に入射すると、面22に向けられたレイ25、26は、面22の特性はもちろんのこと、それらの入射角に依存する方向をもつ選択された数のレイに分割される;図2に描かれたイラストでは、レイ25は2つのレイ27、28に分割され、レイ26はレイ29、30に分割されている。面22が略鏡面的であることから、レイ27と28は、ほとんど同じ反射角をもち、レイ29と30についても同様である。図2に描かれたイラストでは、レイ28と29も面22からほとんど同じ反射角をもつが、一方でレイ27と30の反射角はいくぶんより大きく異なる。
【0017】
レイ27乃至30は、全て面24に入射する。面21、22では、面24から反射したレイの反射角は、面24の面特性はもちろんのこと入射角に依存する。面24は略拡散的であることから、面24から反射されたレイ31乃至38の反射角は、略鏡面的な面21、22におけるものよりもより大きく異なる。しかし、略同じ入射角をもつレイ28や29などのレイでは、一対の関連する反射レイは実質的に同じ反射角をもつ。これはレイの対33、35とレイの対34、36により示される。
【0018】
レイの集合は、レイ24が図2に描かれたイラストのツリーのルートであるレイツリー(ray tree)を含む。レイツリーにおける各レイは、インスタンス番号により識別され、面21に入射するレイ24はインスタンス番号“i”により識別され、続くレイはそのインスタンス番号に基づいたものとされる。したがって、レイ24が面21からの反射において分割されるレイ25と26には、インスタンス番号“i+1”と“i+2”が付けられる。同様に、レイ25が面22からの反射において分割されるレイ27と28には、インスタンス番号“i+2”と“i+3”が付けられ、レイ26が面22からの反射において分割されるレイ29と30には、インスタンス番号“i+3”、“i+4”が付けられる。これまでに図2に描かれたイラストで示されたように、レイ27乃至30の各々は、順番に、それらが略拡散的な面24から反射されるときに選ばれた数のレイ31乃至38に分割され、各レイは入射レイのインスタンス番号に基づいたインスタンス番号をもつ。
【0019】
明らかなように、多くのレイのインスタンス番号は重複しており、重複の程度は各レイがそれぞれの面と相互作用するときに分割されるレイの総数に依存する。図2を調べると、同じインスタンス番号“i+3”をもつ3つのレイ、すなわちレイ28、29、31があることが明らかとなり、そのうちの2つ、すなわちレイ28、29は、面22から反射され面23に入射し、似通った空間的特徴をもつ。更に、同じインスタンス番号“i+4”をもつ4つのレイ、すなわちレイ30、32、33、35が存在し、そのうちの2つ、すなわちレイ33、35は面23から反射され、似通った空間的特徴をもつ。一般に、面21、22が略鏡面的であることから、レイ27乃至30は略同じ方向となり、かなり似通った空間的特徴をもつ。したがって、略拡散的な面24から反射されたレイがいくつかの方向をもつにもかかわらず、面24に入射するレイの多くが似通った角度からのものであることから、面24から反射された多くのレイも似通った方向と他の空間的特徴を備える。インスタンス番号における重複の程度は、各入射レイが面との相互作用において分割されるレイの総数に依存し、反射レイが似通った空間的特徴を備える度合いはそれらが反射される面の特性に依存するものと考えられる。
【0020】
一般に、画像は光源と画像平面のピクセルとの間において追跡されるレイを評価するため積分を評価することによりレンダリングされ、レイは各々の光源と画像平面上のピクセルとの間において光子によって取られ得る径路を表す。積分の評価は、本質的には、画像平面の各ピクセルに入射する光子の流れの強さと色に関する情報を与える。一般に、積分は以下の形式である。
【数3】
Figure 2005505810
上述のように、式(6)(式(1)に等しい)により表わされる積分は、一般に閉じた形式の解(closed form solutions)をもたず、それらは通常数値的に評価される。式(6)によって表わされる積分は、次のように数値的に評価することができる:
【数4】
Figure 2005505810
ここで、総和の中の各項f(xi1 (1),・・・,xin (n))は、積分領域内のサンプル点(xi1 (n),・・・,xin (n))における関数“f”の評価値(evaluation)である。サンプル点は、幾つもの手法を用いて生成することができる;グラベンスタインの出願(Grabenstein application)は、サンプル点がハルトン列(Halton sequences)などのs次元の決定論的低食い違い量列からのサンプルを含む手法について記述し、ケラーの出願(Keller application)は、サンプル点がs次元ハマースレイ点群(ハマースレイ点集合:Hammersley point sets)を含む手法について記述する。一般に、それらの次元s1,・・・snは、それぞれ光源と画像平面との間の光子の軌跡に関連付けられている。
【0021】
問題は特に、面22から反射され面23に入射されるレイ28、29や、面23から反射されるレイ33、35などのように、いくつものレイが実質的に同じ方向をもつ状況において発生し、それは同じサンプル点が積分を評価するのに用いられる可能性が高く、少なくともそれらのレイに関連する積分の部分が実質的に同じ数値の評価値をもつことによる。これに対処するために、式(2)を式(1)の積分の数値的評価に用いる代わりに、積分は次のように評価できる。
【数5】
Figure 2005505810
ここでSij (j)はクランレイ・パターソンローテーション演算子である。すなわち、s1+・・・+sn次元サンプル点Xij (j)までの相応する高い次元を用いて“1”を法としてシフトする(shifts modulo "one")。ここで各サンプル点Xij (j)における初めのsj座標は、式(6)の中のxij (j)に対応する。サンプル点Xij (j)の座標は、ハルトン(Halton)、スクランブルハルトン(scrambled Halton)、あるいは他のいかなる低食い違い量列から取られても良い。添え字“j”の異なる値に対しては、同一あるいは異なる低食い違い量列が用いられても良い。例えば、もしいくつかの分割された添え字の値が予め既知であれば、スクランブルあるいは非スクランブルハマースレイ点を用いてもよい。同じ次元の同じ低食い違い量列を用いた繰り返しシフトを避けるには、sjの値はj=1,・・・,nに対して0よりも大きいことが好ましい。一例として、スクランブルハルトン列が用いられた場合、s1+・・・+sl<d≦s1+・・・+sl+1の次元dに対するサンプル点は、
【数6】
Figure 2005505810
によって与えられ、ここで、各Φb(i,σb)はスクランブル根基逆数(scrambled radical inverse)であり、基数b=b,・・・、及び各基数に対する順列(置換:permutations)σ=σ1,・・・に対して、
【数7】
Figure 2005505810
によって定義される。ケラーの出願では、次のように再帰的に定義される一例としての順列σについて説明されている。基数b=2に対する順列σ2=(0,1)から始め、順列の列は以下のように定義される。
(i)もし、基数“b”が偶数であれば、順列σbは第1に2σb/2の値を取り、2σb/2+1の値を付加することにより生成される。そして、
(ii)もし、基数“b”が奇数であれば、順列σbはσb-1の値を取り、(b−1)/2以上の各値を1インクリメントさせ、b−1の値を真中に挿入することにより生成される。
この再帰的な生成手段は、以下のような順列をもたらす。
σ2=(0,1)
σ3=(0,1,2)
σ4=(0,2,1,3)
σ5=(0,3,2,1,4)
σ6=(0,2,4,1,3,5)
σ7=(0,2,5,3,1,4,6)
σ8=(0,4,2,6,1,5,3,7)…
したがって、基数“b”での値“i”に対する根基逆数が与えられ、もし根基逆数を表す「k番目」の桁が値“j”を持つならば、スクランブル根基逆数の「k番目」の桁は、上記順列σbでの“j番目”の桁の値に対応する値を持つ。
【0022】
一般に、レイトレーシングにおける再帰レベルは、トレースデプス(trace depth)を超えず、これは通常、光源と画像平面との間の経路を含むレイの数に関係する。スクランブル、非スクランブルハルトン列及びハマースレイ点群は、極めて高速に生成することができ、これらは、それらがそれ自身極めて高速である根基逆関数(radical inverse function)に基づくためである。かなり深いトレースデプスであっても、サンプル点の座標を生成するコストは、一般に各々のサンプル点に対する関数“f”の数値を生成するコストよりも安く、これは、数値的な評価に加えて、交差判定(intersection test)や、濃度(density)やテクスチャーなどの照合(look-ups)をも含むことができる。
【0023】
本発明によれば、サンプル点の座標のキャッシングを利用する手法が以下に説明され、レイトレーシングがプロセッサモジュール11(図1)により提供される複数のプロセッサに渡って分散される分散レイトレーシングにおいてどのようにキャッシュされたサンプリングが用いられるかが説明される。キャッシングの優位性は、ハルトン列の代わりに低食い違い量列を生成するためにより時間が掛かる手法が用いられた場合により明らかとなる。この新しい手法はコードセグメント(code segment)1との関連で説明される。
【0024】
コードセグメント1
Figure 2005505810
Figure 2005505810
Figure 2005505810
【0025】
コードセグメント1は2つのクラス、すなわちQMCContextとRayを含む。ライン(1)から(13)を含む1つのクラス、クラスQMCContextは、2つのメソッドを提供する。すなわち、メソッドget_sampleとメソッドsetup_child_contextである。ライン(8)から(17)を含むメソッドget_sampleは、式(8)で用いられるサンプル点を提供する。各コンテキストレベルは、本質的にレイツリー(図2参照)でのレベルと関係する。メソッドget_sampleは、再帰的(recursively)に作用される。初め、ライン(10)では再帰レベルがルートであるか否かが判定され、これは親がいない場合に当たる。もし、親レベルが存在するならば、メソッドget_sampleは、キャッシュがサンプル点をもつか否かを判定するために再帰において親レベルのために保持されるキャッシュをチェックする(ライン11)。もし、キャッシュがサンプル点をもつならば、それはそのサンプル点を用い、ライン(16)において現在(current)の再帰レベルのためにサンプル点をアップデートし、以下において説明されるクラスRayにおけるメソッドにアップデートされたサンプル点を与える。一方、もしキャッシュがサンプル点をもたなかったならば、メソッドget_sampleは、ライン(14)において、親の親レベルからサンプル点を生成し、ライン(15)において親のキャッシュにサンプル点を加える。その後、ライン(16)では、メソッドget_sampleがライン(14)において生成されたサンプル点をライン(16)における現在の再帰レベルのためにアップデートし、クラスRayにおけるメソッドにアップデートされたサンプル点を与える。最後に、もし、ライン(10)において、親が存在しない(再帰レベルがルートレベルである場合)と判定されると、メソッドget_sampleは、ライン(16)においてサンプル点を生成し、クラスRayにおけるメソッドにこのサンプル点を与える。
【0026】
クラスQMCContextのライン(18)から(23)を含むメソッドsetup_child_contextは、次の低いレベルのために(もし存在すれば)クラスQMCContextにより与えられるメソッドを初期化する。
【0027】
ライン(24)から(36)を含むクラスRayは、クラスQMCContextにより提供されるメソッドを初期化し、更にメソッドShadeを含む。ライン(27)から(36)を含み、特定のレイに対する式(8)を評価するために用いられるメソッドShadeは、クラスQMCContextにある両方のメソッドを呼び出す。メソッドShadeは、ライン(32)において、式(8)との関連で用いられるサンプルを生成するためにライン(33)で用いられるサンプル点を与えるため、クラスQMCContextのメソッドget_sampleを呼び出す。ライン(34)において、メソッドShadeは、次のコンテキストレベル(もし存在すれば)のためのクラスQMCContextを用意するためにクラスQMCContextにおけるメソッドsetup_child_contextを呼び出す。
【0028】
本発明は、画素値を表わす積分を数値的に評価するためのサンプル点を与えるために決定論的な低食い違い量列と、クランレイ・パターソンローテーションなどのリカーシブルローテーションを用いてレンダリングされた画像における画素の画素値を生成する新奇で改良されたコンピュータグラフィックスシステムを提供する。
【0029】
本発明は、リカーシブなクランレイ・パターソンローテーションとの関連において説明されたが、他の形式のローテーションが用いられてもよい。
【0030】
本発明に基づくシステムは、特定の目的のためのハードウェアや汎用目的のコンピュータシステムから全体又は一部分が構成されてもよく、これらのいかなる組み合せ、あるいは適切なプログラムにより制御されるこれらのどのような構成部分であってもよい。いかなるプログラムも、システムの全体又は一部分であってもよく、あるいは従来の方法でこのシステムに保存されていてもよく、またそれは全体あるいは一部分が、ネットワークあるいは従来の方法により情報を転送する他の仕組み(mechanism)を介してシステムの中に提供されていてもよい。更に、このシステムは操作され、かつ/または、さもなければオペレータ入力部(不図示)を用いてオペレータにより与えられた情報を利用して制御される。オペレータ入力部は、システムに直接接続されるかもしれないし、ネットワークあるいは従来の方法により情報を転送する他の仕組みを介してシステムに情報を転送してもよい。
【0031】
これまでの説明は、この発明の特定の実施形態に限られている。しかし、本発明のいくつかの、あるいは全ての優位性を達成したまま様々な変形と変更を本発明に対して行ない得る。本発明の真の精神と範囲内において、これらとそのような他の変形及び変更をカバーすることが添付された請求項の目的である。
【図面の簡単な説明】
【0032】
【図1】本発明にしたがって構成された例示的なコンピュータグラフィックスシステムを示す。
【図2】本発明を理解するのに役立つ例示的なレイの追跡を示す。

Claims (9)

  1. 画像における画素の画素値を生成するためのコンピュータグラフィックスシステムであって、前記画素はシーンの中の1点を表し、前記コンピュータグラフィックスシステムは前記画素値を選択された関数の積分の評価により生成し、
    A. 複数の追跡レベルを持つように構成されたレイの追跡において各々一連のレイの1つに関連付けられたサンプル点の各集合を生成するように構成されたサンプル点生成手段であって、少なくとも1つのレベルにある前記レイが複数のレイに分割され、各レイがレイインスタンス識別子に関連付けられ、各レベルに対して選択されたローテーション演算子が再帰的に適用される所定の決定論的な低食い違い量列として前記サンプル点を生成するように構成された前記サンプル点生成手段と、
    B. 前記サンプル点生成手段により生成された前記サンプル点の1つにおける前記選択された関数の評価を各々代表する複数の関数値を生成するように構成され、前記画素値を生成するのに前記関数値を用いる関数評価手段と
    を備えることを特徴とするレイ追跡コンピュータグラフィックスシステム。
  2. 前記選択されたローテーション演算子がクランレイ・パターソンローテーション演算子であることを特徴とする請求項1に記載のコンピュータグラフィックスシステム。
  3. 前記サンプル点生成手段が、各追跡レベルに対して生成されたサンプル点をキャッシュするように構成され、前記キャッシュされたサンプル点を続く追跡レベルのためのサンプル点の生成において用いることを特徴とした請求項1に記載のコンピュータグラフィックスシステム。
  4. 画像における画素の画素値を生成するコンピュータグラフィックスの方法であって、前記画素はシーンの中の1点を表し、前記コンピュータグラフィックスシステムは前記画素値を選択された関数の積分の評価により生成し、
    A. 複数の追跡レベルを持つように構成されたレイの追跡において各々一連のレイの1つに関連付けられたサンプル点の各集合を生成するサンプル点生成ステップであって、少なくとも1つのレベルにある前記レイが複数のレイに分割され、各レイがレイインスタンス識別子に関連付けられ、各レベルに対して選択されたローテーション演算子が再帰的に適用される所定の決定論的な低食い違い量列として前記サンプル点を生成するステップを含む前記サンプル点生成ステップと、
    B. 前記サンプル点生成ステップにおいて生成された前記サンプル点の1つにおける前記選択された関数の評価を各々代表する複数の関数値を生成し、前記画素値を生成するのに前記関数値を用いる関数評価ステップと
    を備えることを特徴とするコンピュータグラフィックスの方法。
  5. 前記選択されたローテーション演算子がクランレイ・パターソンローテーション演算子であることを特徴とする請求項4に記載のコンピュータグラフィックスの方法。
  6. 前記サンプル点生成ステップが、各追跡レベルに対して生成されたサンプル点をキャッシュするステップを含み、前記キャッシュされたサンプル点を続く追跡レベルのためのサンプル点の生成において用いることを特徴とした請求項4に記載のコンピュータグラフィックスの方法。
  7. 画像における画素の画素値を生成するコンピュータグラフィックスシステムを提供するコンピュータに関連して用いられるコンピュータプログラムであって、前記画素はシーンの中の1点を表し、前記コンピュータグラフィックスシステムは前記画素値を選択された関数の積分の評価により生成し、
    A. 複数の追跡レベルを持つように構成されたレイの追跡において前記コンピュータが各々一連のレイの1つに関連付けられたサンプル点の各集合を生成できるように構成されたサンプル点生成モジュールであって、少なくとも1つのレベルにある前記レイが複数のレイに分割され、各レイがレイインスタンス識別子に関連付けられ、各レベルに対して選択されたローテーション演算子が再帰的に適用される所定の決定論的な低食い違い量列として前記サンプル点を生成するように構成された前記サンプル点生成モジュールと、
    B. 前記サンプル点生成モジュールにより生成された前記サンプル点の1つにおける前記選択された関数の評価を各々代表する複数の関数値を前記コンピュータが生成できるように構成され、前記画素値を生成するのに前記関数値を用いる関数評価モジュールと
    を備えることを特徴とするコンピュータプログラム。
  8. 前記選択されたローテーション演算子がクランレイ・パターソンローテーション演算子であることを特徴とする請求項7に記載のコンピュータプログラム。
  9. 前記サンプル点生成モジュールが、各追跡レベルに対して生成されたサンプル点を前記コンピュータがキャッシュできるように構成され、前記キャッシュされたサンプル点を続く追跡レベルのためのサンプル点の生成において用いることを特徴とした請求項1に記載のコンピュータプログラム。
JP2003502789A 2001-06-07 2002-06-07 サンプル点生成のためのリカーシブローテーションを用いた決定論的手法を用いた画像レンダリング Expired - Lifetime JP4133806B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US29713901P 2001-06-07 2001-06-07
PCT/IB2002/003077 WO2002099754A1 (en) 2001-06-07 2002-06-07 Image rendering using strictly deterministic methodologies using recursive rotations for generating sample points

Publications (2)

Publication Number Publication Date
JP2005505810A true JP2005505810A (ja) 2005-02-24
JP4133806B2 JP4133806B2 (ja) 2008-08-13

Family

ID=23145021

Family Applications (3)

Application Number Title Priority Date Filing Date
JP2003502789A Expired - Lifetime JP4133806B2 (ja) 2001-06-07 2002-06-07 サンプル点生成のためのリカーシブローテーションを用いた決定論的手法を用いた画像レンダリング
JP2003502788A Expired - Lifetime JP4100625B2 (ja) 2001-06-07 2002-06-07 サンプル・ポイントのコースシーケンスを発生させるための確定的な方法論を用いた画像表示
JP2003504335A Expired - Lifetime JP4344606B2 (ja) 2001-06-07 2002-06-07 グローバルイルミネーションを評価するためのロシアンルーレット法を用いた描画

Family Applications After (2)

Application Number Title Priority Date Filing Date
JP2003502788A Expired - Lifetime JP4100625B2 (ja) 2001-06-07 2002-06-07 サンプル・ポイントのコースシーケンスを発生させるための確定的な方法論を用いた画像表示
JP2003504335A Expired - Lifetime JP4344606B2 (ja) 2001-06-07 2002-06-07 グローバルイルミネーションを評価するためのロシアンルーレット法を用いた描画

Country Status (9)

Country Link
US (3) US7046244B2 (ja)
EP (3) EP1397782B1 (ja)
JP (3) JP4133806B2 (ja)
AT (3) ATE373288T1 (ja)
AU (3) AU2002329516B2 (ja)
CA (3) CA2449684A1 (ja)
DE (3) DE60229129D1 (ja)
ES (3) ES2311067T3 (ja)
WO (3) WO2002101657A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010501100A (ja) * 2006-08-15 2010-01-14 メンタル イメージズ ゲーエムベーハー 準モンテカルロ法を使用するマルコフ連鎖の同時シミュレーション
CN114646290A (zh) * 2022-03-02 2022-06-21 中国地质调查局西安矿产资源调查中心 一种地球物理勘探野外点位放样方法

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7502027B1 (en) * 1999-09-13 2009-03-10 Solidworks Corporation Electronic drawing viewer
US7952583B2 (en) * 2000-06-19 2011-05-31 Mental Images Gmbh Quasi-monte carlo light transport simulation by efficient ray tracing
US7659894B2 (en) * 2000-06-19 2010-02-09 Mental Images Gmbh Terminating spatial partition hierarchies by a priori bounding memory
US7499053B2 (en) * 2000-06-19 2009-03-03 Mental Images Gmbh Real-time precision ray tracing
US8248416B2 (en) * 2000-06-19 2012-08-21 Mental Images Gmbh Efficient ray tracing without acceleration data structure
US7184042B2 (en) * 2000-06-19 2007-02-27 Mental Images Gmbh Computer graphic system and computer-implemented method for generating images using a ray tracing methodology that makes use of a ray tree generated using low-discrepancy sequences and ray tracer for use therewith
US7773088B2 (en) * 2000-06-19 2010-08-10 Mental Images Gmbh Simultaneous simulation of markov chains using quasi-monte carlo techniques
US8188997B2 (en) * 2000-06-19 2012-05-29 Mental Images Gmbh Accelerated ray tracing using shallow bounding volume hierarchies
US20050278406A1 (en) * 2002-05-15 2005-12-15 Alexander Keller System and computer-implemented method for evaluating integrals using stratification by rank-1 lattices
US7589729B2 (en) * 2002-05-15 2009-09-15 Mental Images Gmbh Image synthesis by rank-1 lattices
US7567248B1 (en) * 2004-04-28 2009-07-28 Mark William R System and method for computing intersections between rays and surfaces
US7298370B1 (en) * 2005-04-16 2007-11-20 Apple Inc. Depth ordering of planes and displaying interconnects having an appearance indicating data characteristics
US8698844B1 (en) * 2005-04-16 2014-04-15 Apple Inc. Processing cursor movements in a graphical user interface of a multimedia application
WO2007002494A2 (en) * 2005-06-23 2007-01-04 Mental Images Gmbh Real-time precision ray tracing
WO2008076852A1 (en) * 2006-12-14 2008-06-26 Mental Images Gmbh Computer graphics using meshless finite elements for light transport
US7755628B2 (en) 2006-12-29 2010-07-13 Intel Corporation Method and apparatus for multi-level ray tracing
US8623019B2 (en) * 2007-07-03 2014-01-07 Pioneer Surgical Technology, Inc. Bone plate system
US9619917B2 (en) 2008-10-03 2017-04-11 Apple Inc. Depth of field for a camera in a media-editing application
US8131770B2 (en) * 2009-01-30 2012-03-06 Nvidia Corporation System, method, and computer program product for importance sampling of partitioned domains
US8266623B2 (en) 2009-04-29 2012-09-11 Nvidia Corporation System, method, and computer program product for decomposing a sampling task into a plurality of jobs
KR101661166B1 (ko) * 2010-06-14 2016-09-29 연세대학교 산학협력단 3차원 영상 시스템에서 광선 추적 방법 및 장치
US8860725B2 (en) 2010-08-13 2014-10-14 Nvidia Corporation System, method, and computer program product for deterministically simulating light transport
JP2012181825A (ja) * 2011-02-09 2012-09-20 Canon Inc 画像処理装置およびその方法
JP5839907B2 (ja) * 2011-09-15 2016-01-06 キヤノン株式会社 画像処理装置および画像処理方法
US8847957B1 (en) * 2011-10-05 2014-09-30 Nvidia Corporation Divide-and-conquer system, method, and computer program product for providing photon mapping
CN102496170A (zh) * 2011-12-06 2012-06-13 阳赛 一种用于采样任务分解的方法
CN102496173A (zh) * 2011-12-06 2012-06-13 阳赛 一种用于采样点分离的数字偏移方法
US9202139B2 (en) * 2012-07-18 2015-12-01 Nvidia Corporation System, method, and computer program product for generating a subset of a low discrepancy sequence
US9367955B2 (en) * 2012-11-26 2016-06-14 Nvidia Corporation System, method, and computer program product for tiled screen space sample scrambling for parallel deterministic consistent light transport simulation
KR102223064B1 (ko) * 2014-03-18 2021-03-04 삼성전자주식회사 영상 처리 장치 및 방법
US9679398B2 (en) * 2015-10-19 2017-06-13 Chaos Software Ltd. Rendering images using color contribution values of render elements

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6529193B1 (en) * 1996-06-25 2003-03-04 Mental Images Gmbh & Co. Kg System and method for generating pixel values for pixels in an image using strictly deterministic methodologies for generating sample points
CA2258928A1 (en) 1996-06-25 1997-12-31 Mental Images G.M.B.H. & Co., Kg. System and method for generating pixel values for pixels in an image using strictly deterministic methodologies for generating sample points
GB9713186D0 (en) * 1997-06-24 1997-08-27 Univ Sheffield Artificial joints
US6664961B2 (en) * 2000-12-20 2003-12-16 Rutgers, The State University Of Nj Resample and composite engine for real-time volume rendering

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010501100A (ja) * 2006-08-15 2010-01-14 メンタル イメージズ ゲーエムベーハー 準モンテカルロ法を使用するマルコフ連鎖の同時シミュレーション
CN114646290A (zh) * 2022-03-02 2022-06-21 中国地质调查局西安矿产资源调查中心 一种地球物理勘探野外点位放样方法
CN114646290B (zh) * 2022-03-02 2023-08-25 中国地质调查局西安矿产资源调查中心 一种地球物理勘探野外点位放样方法

Also Published As

Publication number Publication date
US6911976B2 (en) 2005-06-28
EP1419485A1 (en) 2004-05-19
JP2005506607A (ja) 2005-03-03
JP4344606B2 (ja) 2009-10-14
ATE402458T1 (de) 2008-08-15
DE60227811D1 (de) 2008-09-04
AU2002329516B2 (en) 2008-10-02
WO2002099751A1 (en) 2002-12-12
US20030034968A1 (en) 2003-02-20
CA2449690C (en) 2008-03-25
JP2005503607A (ja) 2005-02-03
US7046244B2 (en) 2006-05-16
ES2311067T3 (es) 2009-02-01
AU2002319869B2 (en) 2008-10-02
CA2449684A1 (en) 2002-12-12
US6885370B2 (en) 2005-04-26
AU2002348894B2 (en) 2008-10-16
JP4133806B2 (ja) 2008-08-13
US20030052874A1 (en) 2003-03-20
CA2449282A1 (en) 2002-12-19
US20030063082A1 (en) 2003-04-03
CA2449690A1 (en) 2002-12-12
ES2290320T3 (es) 2008-02-16
EP1405271B1 (en) 2008-07-23
DE60222437T2 (de) 2008-06-12
ES2315396T3 (es) 2009-04-01
EP1397782A1 (en) 2004-03-17
ATE373288T1 (de) 2007-09-15
CA2449282C (en) 2010-01-05
EP1397782B1 (en) 2008-10-01
EP1419485B1 (en) 2007-09-12
WO2002099754A1 (en) 2002-12-12
EP1405271A1 (en) 2004-04-07
WO2002101657A1 (en) 2002-12-19
JP4100625B2 (ja) 2008-06-11
DE60222437D1 (de) 2007-10-25
DE60229129D1 (de) 2008-11-13
ATE409928T1 (de) 2008-10-15

Similar Documents

Publication Publication Date Title
JP4133806B2 (ja) サンプル点生成のためのリカーシブローテーションを用いた決定論的手法を用いた画像レンダリング
AU2002319869A1 (en) Image rendering using strictly deterministic methodologies using recursive rotations for generating sample points
JP4072209B2 (ja) 画素値を発生する装置および方法
US7227547B2 (en) Method for generating pixel values for pixels in an image using strictly deterministic methodologies for generating sample points
US7453461B2 (en) Image generation using low-discrepancy sequences
US7425957B2 (en) Computer graphics methods and systems using quasi-monte carlo methodology
US20090051688A1 (en) Generating images using ray tracing and ray tree generated using low-discrepancy sequences
US7499054B2 (en) Generating images using multiple photon maps
US20070211051A1 (en) Computer Graphics Systems, Methods and Computer Program Products Using Sample Points Determined Using Low-Discrepancy Sequences
US20050278406A1 (en) System and computer-implemented method for evaluating integrals using stratification by rank-1 lattices
JP4291395B2 (ja) 標本点を発生するための厳密に決定論的な方法を用いて画像中の画素の画素値を発生する装置及び方法
EP1628263B1 (en) Generating pixel values using strictly deterministic methodologies for generating sample points

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20050224

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070821

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20071119

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20071121

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20071203

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20071219

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20080107

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20080118

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20080125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080318

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080502

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20080527

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20080602

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110606

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4133806

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110606

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120606

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120606

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130606

Year of fee payment: 5

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term