JP2019521417A - シーンを表す点群を符号化する方法、符号化器システム、及びプログラムを記憶した非一時的コンピューター可読記録媒体 - Google Patents

シーンを表す点群を符号化する方法、符号化器システム、及びプログラムを記憶した非一時的コンピューター可読記録媒体 Download PDF

Info

Publication number
JP2019521417A
JP2019521417A JP2018559911A JP2018559911A JP2019521417A JP 2019521417 A JP2019521417 A JP 2019521417A JP 2018559911 A JP2018559911 A JP 2018559911A JP 2018559911 A JP2018559911 A JP 2018559911A JP 2019521417 A JP2019521417 A JP 2019521417A
Authority
JP
Japan
Prior art keywords
point
points
input
residual data
model parameters
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
JP2018559911A
Other languages
English (en)
Other versions
JP6676193B2 (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.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
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 Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Publication of JP2019521417A publication Critical patent/JP2019521417A/ja
Application granted granted Critical
Publication of JP6676193B2 publication Critical patent/JP6676193B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T9/00Image coding
    • G06T9/001Model-based coding, e.g. wire frame
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/10Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T9/00Image coding
    • G06T9/004Predictors, e.g. intraframe, interframe coding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T9/00Image coding
    • G06T9/005Statistical coding, e.g. Huffman, run length coding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T9/00Image coding
    • G06T9/007Transform coding, e.g. discrete cosine transform
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00Indexing scheme for image data processing or generation, in general
    • G06T2200/04Indexing scheme for image data processing or generation, in general involving 3D image data

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Geometry (AREA)
  • Discrete Mathematics (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)
  • Image Processing (AREA)
  • Image Generation (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

メモリと通信するプロセッサを備える符号化器を用いて、シーンを表す点群を符号化する方法は、パラメーター化された表面を、3次元空間におけるロケーションを表す入力点を有する点群上にフィッティングするステップと、パラメーター化された表面からモデルパラメーターを生成するステップと、パラメーター化された表面から対応する点を計算するステップであって、この対応する点は、入力点に対応する、ステップと、対応する点と点群の入力点とに基づいて残差データを計算するステップと、モデルパラメーター及び残差データを圧縮して、それぞれ、符号化されたモデルパラメーター及び符号化された残差データを得るステップと、パラメーター化された表面の符号化されたモデルパラメーターと符号化された残差データとからビットストリームを生成するステップとを含む。

Description

本発明は、包括的には、点群のジオメトリを圧縮する方法に関し、より詳細には、この点群のジオメトリを予測的に符号化する方法に関する。
点群は、3D空間における点のセットを含むことができる。例えば、所与の点は、そのロケーション又はジオメトリを指定する具体的な(x,y,z)座標を有することができる。この点ロケーションは、センサー分解能によって決定されるか、又は点群を生成するのに実行される任意の前処理によって決定される分解能で、3D空間における任意の箇所に配置することができる。細かい分解能について、又は広大な3D空間に拡がる点群を有する粗い分解能について、点ロケーションの整数バイナリ表現又は浮動小数点バイナリ表現には、複数のビットを要し得る。高精度で全ての座標を記憶又はシグナリングすることにより、捕捉システムが、高忠実度で点群座標を保存することが可能になるが、その一方で、そのような表現は、保存する際には記憶空間を、又はシグナリングする際には帯域幅を大量に消費する可能性がある。したがって、後の記憶又は送信のために、点群ロケーション又はジオメトリを表すのに必要なビット数を低減する必要がある。圧縮を伴う場合でも、圧縮された表現のサイズは、著しく大きなものとなる可能性があり、したがって、点群ジオメトリの粗い表現を、ファイル又はビットストリームの全体を復号する必要なく、迅速に又は容易に復号及び再構成することが可能になる圧縮された表現もまた必要とされている。
本発明のいくつかの実施形態は、所与の表面をパラメーター化して、このパラメーター化された表面を、点群上にフィッティングすることによって、この点群を効果的に符号化することができるという認識に基づいている。
したがって、1つの実施形態は、メモリと通信するプロセッサを備える符号化器を用いて、シーンを表す点群を符号化する方法であって、この方法は、
パラメーター化された表面を、入力点によって形成される点群上にフィッティングするステップと、
パラメーター化された表面からモデルパラメーターを生成するステップと、
パラメーター化された表面上の対応する点を計算するステップであって、この対応する点は、入力点に対応する、ステップと、
対応する点と点群の入力点とに基づいて残差データを計算するステップと、
モデルパラメーター及び残差データを圧縮して、符号化されたモデルパラメーター及び符号化された残差データを得るステップと、
パラメーター化された表面の符号化されたモデルパラメーターと符号化された残差データとからビットストリームを生成するステップと、
を含む、方法を開示する。
さらに、本発明のいくつかの実施形態は、パラメーター化された表面のモデルパラメーターと元の点群から計算された残差データとを含むビットストリームを受信することによって、符号化された点群データを効果的に復号することができるという認識に基づいている。
したがって、1つの実施形態は、シーンを表す点群を符号化する符号化器システムであって、点群の各点は、3次元(3D)空間におけるロケーションであり、この符号化器システムは、
メモリと通信するプロセッサと、
メモリに記憶される符号化器モジュールであって、ステップを実行することによってシーンを表す点群を符号化するように構成される、符号化器モジュールと、
を備え、ステップは、
パラメーター化された表面を、入力点によって形成される点群上にフィッティングするステップと、
パラメーター化された表面からモデルパラメーターを生成するステップと、
パラメーター化された表面上の対応する点を計算するステップであって、この対応する点は、入力点に対応する、ステップと、
対応する点と点群の入力点とに基づいて残差データを計算するステップと、
モデルパラメーター及び残差データを圧縮して、それぞれ、符号化されたモデルパラメーター及び符号化された残差データを得るステップと、
パラメーター化された表面の符号化されたモデルパラメーターと符号化された残差データとからビットストリームを生成するステップと、
を含む、符号化器システムを開示する。
本発明の実施形態による符号化プロセスのブロック図である。 本発明の実施形態による、表面モデルと編成された点群についての点ロケーションとの間の対応関係を計算するプロセスのブロック図である。 本発明の実施形態による、表面モデルと編成されていない点群についての点ロケーションとの間の対応関係を計算するプロセスのブロック図である。 本発明の実施形態による、点ロケーションと表面モデルの対応する点との間の残差を計算するプロセスである。 本発明の実施形態による復号プロセスの図である。 本発明の実施形態による、階層的分割プロセスの図である。
本発明の様々な実施形態が、図面を参照して以下で説明される。図面は縮尺どおり描かれておらず、類似の構造又は機能の要素は、図面全体にわたって同様の参照符号によって表されることに留意されたい。図面は、本発明の特定の実施形態の説明を容易にすることのみを意図することにも留意されたい。図面は、本発明の網羅的な説明として意図されるものでもなければ、本発明の範囲を限定するものとして意図されるものでもない。加えて、本発明の特定の実施形態と併せて説明される態様は、必ずしもその実施形態に限定されず、本発明の任意の他の実施形態において実施することができる。
本発明の実施形態は、3次元(3D)点群を、3D点の点ロケーションの予測子として機能する点の群又はセット上にフィッティングされた表面のパラメトリックモデルを用いて圧縮し、パラメーター及び/又は予測誤差を量子化及びシグナリングすることによって圧縮を達成する、方法及びシステムを提供する。
いくつかの実施形態は、メモリと通信するプロセッサを備える符号化器を用いて、シーンを表す点群を符号化する方法であって、点群の各点は、3次元(3D)空間におけるロケーションであり、この方法は、
パラメーター化された表面を、入力点によって形成される点群上にフィッティングするステップと、
パラメーター化された表面からモデルパラメーターを生成するステップと、
パラメーター化された表面上の対応する点を計算するステップであって、この対応する点は、入力点に対応する、ステップと、
対応する点と点群の入力点とに基づいて残差データを計算するステップと、
モデルパラメーター及び残差データを圧縮して、符号化されたモデルパラメーター及び符号化された残差データを得るステップと、
パラメーター化された表面の符号化されたモデルパラメーターと符号化された残差データとからビットストリームを生成するステップと、
を含む、方法を開示する。
符号化プロセス
図1は、本発明のいくつかの実施形態による、符号化器100によって実行される符号化プロセスのブロック図である。符号化器100は、この符号化プロセスを実行するために、メモリ(図示せず)と通信するプロセッサ(図示せず)を備える。
符号化器100への入力は、N個の点p、i={1,2,...,N}のセットを含む点群101であり、ここで、各pは、3D空間における点ロケーションである。点ロケーションは、3D空間における座標によって表すことができる。例えば、デカルト座標系を用いる場合、pは、三つ組{x,y,z}によって表すことができ、ここで、x、y、zは、点の座標である。別の例では、球面座標系を用いることができ、この球面座標系において各点pは、三つ組{r,θi,φi}によって表され、ここで、rは、動径を示し、θiは、極角を示し、φiは、点の方位角を表す。
入力点p101のセット又はサブセットについて、表面モデルフィッティングプロセス102が、入力点p101のロケーション又は入力点によって表されるオブジェクトの表面を近似する表面モデルを計算する。例えば、入力点が3D空間における球面オブジェクトの表現又は境界ロケーションである場合、表面モデル104は、2つのパラメーター、すなわち、半径と球の原点のロケーションの座標とを有する球面の方程式とすることができる。別の例では、表面モデル104は、ベジェ曲面(Bezier surface:ベジェ表面)又はベジェ曲面パッチとすることができる。この場合、表面モデル104は、3D空間における位置又は座標のパラメトリック表現である、p(u,v)によって表すことができる。ベジェ曲面パッチは、以下の方程式によって表すことができる。
Figure 2019521417
ここで、p(u,v)は、2つのパラメーターを用いる、3D空間におけるロケーション又は座標のパラメトリック表現である。(u,v)が単位正方形であるパラメーターu及びvを所与とすると、p(u,v)は、3D空間における(x,y,z)座標である。ベジェ曲面モデルについて、
Figure 2019521417
及び
Figure 2019521417
は、
Figure 2019521417
の形態のバーンスタイン多項式であり、以下の二項係数、
Figure 2019521417
を用いる。
ベジェ曲面パッチの形状は、モデルパラメーターbijによって求められ、このモデルパラメーターは、「制御点」として知られる。制御点は、整数i及びjによってインデックス付けされ、ここで、iは、u次元に沿った0〜nのインデックスであり、jは、v次元に沿った0〜mのインデックスである。例えば、m=3かつn=3の場合、全部で16個の制御点が存在する。制御点bijは、モデルパラメーター103に含まれる。
後に説明される、編成された点群について、予め設定された編成情報117は、点が編成される方法を指定し、この編成情報は、表面モデルフィッティングプロセス102に入力することができる。この編成情報117は、各入力点p101と、編成行列、グリッド、ベクトルにおけるこの入力点のロケーション、又はメモリにおけるアドレスとの間のマッピングの明細とすることができる。いくつかの場合、編成情報117は、インデックス、又は各入力点p101のメモリにおけるアドレス若しくは順序に基づいて推測することができる。
表面モデルフィッティングプロセス102は、歪みを最小化するようにモデルパラメーター103を選択することができる。歪みは、各入力点pと表面モデル104上の各対応する点fmj106との間の全歪み又は平均歪みによって表すことができる。この場合、表面モデル104は、入力点p101のセット又はサブセットの予測として機能する。また、表面モデルフィッティングプロセス102は、再構成誤差も最小化することができる。再構成誤差は、各入力点pと対応する再構成された点との間の全誤差又は平均誤差である。この場合、再構成された点は、復号プロセスにおいて再構成された点と同一であり、このことは、後に説明される。
モデルパラメーター103が計算されると、点のセット又はサブセットについての表面モデル104が、これらのモデルパラメーターを用いて生成される。代替的に、表面モデル104は、モデルパラメーター103を計算する前の、表面モデルフィッティングプロセス102中に生成することができる。この場合、表面モデル104は、既に利用可能であり、再生成される必要がない。
表面モデル104は、連続的モデルf105又は離散的モデル106とすることができる。連続的モデルは、単位正方形における(u,v)について連続値を使用し、その結果、表面モデル
Figure 2019521417
をもたらし、ここで、(f,f,f)は、デカルト空間における表面モデル104のロケーションを表す。離散的モデルについて、3D空間における表面モデル104のロケーションのパラメトリック表現p(u,v)は、離散的パラメーター(u,v)を使用する(ここで、iは、一般的なインデックス項として用いられ、上述されたベジェ曲面パッチ方程式において用いられていたiとは異なるインデックスであることに留意されたい)。離散的モデルにおいて、表面モデル104は、点
Figure 2019521417
106のセットによって表され、ここで、各対応する点fm106は、連続的表面モデルfの離散点を表す。符号化プロセスのこの部分の意図は、N個の入力点のセット又はサブセットの良好な表現であるM個の表面モデル点のセットを生成することである。
表面モデル14が取得されると、符号化器100は、残差計算プロセス107を実行して、各入力点p101に対応する残差
Figure 2019521417
108のセットを計算する。各残差rは、予測誤差として用いることができるが、これは何故なら、表面モデルは、入力点の予測であるものと考えることができるためである。いくつかの場合、予測誤差は、残差に基づいて計算することができ、例えば、残差rが各座標間の差分の三つ組によって表される場合(この場合、r={x−f ,y−f ,z−f }、i={1,2,...,N}である)、予測誤差は、残差の各成分の二乗を合計することによる二乗誤差歪みを用いて計算することができる。
残差r108を計算するために、各入力点p101は、表面モデルfm106上に対応する点を有さなければならない。N個の入力点p101及びM個の表面モデル点fm106は、表面モデル104上の対応する点f110を計算するプロセス109に入力される。別の実施形態では、連続的表面f105の点を、対応関係計算プロセス109において用いることができる。プロセス109は、N個の対応する表面モデル点
Figure 2019521417
110のセットを出力し、ここにおいて、所与の{1,2,...,N}におけるiについて、f110は、入力点p101に対応する表面モデル104上の点を示す。
表面モデル104上の対応する点f110が計算されると、上述されたように、残差計算プロセス107において、対応する点f110と入力点p101との間の残差r108が計算される。
残差r、i={1,2,...,N}108を変換プロセス111に入力して、変換係数112のセットを生成することができる。例えば、変換プロセス111は、変換係数112のセットを出力するのに、残差r108に、離散コサイン変換(DCT)又は他の空間領域−周波数領域変換を適用することができる。いくつかの場合、変換プロセス111は、パススルー(pass-through:素通りの)変換プロセスとすることができ、ここにおいて、データを改変する動作は行われず、変換係数112は残差r108と同一のものとされる。
変換プロセス111の後、変換係数112は、量子化プロセス113に入力される。量子化プロセスは、変換係数112を量子化して、量子化された変換係数114のセットを出力する。量子化プロセス113の目的は、データの各要素が多数の異なる値を含む空間からの或る値を有する可能性のある入力データのセットを、データの各要素がより少数の異なる値を含む空間からの或る値を有する出力データのセットによって表すことである。例えば、量子化プロセス113は、浮動小数点の入力データを整数の出力データに量子化することができる。量子化プロセス113が可逆的なものである場合、ロスレス圧縮を達成することができる。これは、量子化プロセス113からの各々の生じ得る出力要素を、対応する量子化前の入力要素に逆マッピングすることができることを示す。
続いて、量子化された変換係数114は、エントロピー符号化器115に入力される。エントロピー符号化器115は、変換係数115のバイナリ表現を出力ビットストリーム116に出力する。出力ビットストリーム116は、続いて、送信するか、又はメモリ若しくはコンピューターファイルに記憶することができる。モデルパラメーター103もまた、ビットストリーム116に出力するために、エントロピー符号化器115に入力される。
別の実施形態では、エントロピー符号化器115は、固定長の符号化器を含むことができる。
表面モデル上の対応する点を計算するプロセス
図2及び図3は、対応する点の計算プロセス200及び300のブロック図を示している。対応する点の計算プロセスには、2つのタイプがある。対応する点の計算プロセス200は、編成された点群に適用され、対応する点の計算プロセス300は、編成されていない点群に適用される。
表面モデル上の対応する点を計算するプロセス109の目的は、後の段階で各入力点p101と表面モデル104上の各対応する点f110との間の残差r108を計算することができるように、入力点p101ごとに、表面モデル104上の対応する点f110を求めることである。プロセス109において表面モデル104上の対応する点を計算するこのプロセスは、編成された点群又は編成されていない点群上で行うことができる。これらのプロセスは、編成された点群及び編成されていない点群の以下の説明の後に、詳述される。
編成された点群において、各点p101は、メモリ又は1つ以上の要素を含む行列若しくはベクトルにおける位置又はアドレスに関連付けられる。例えば、3×4行列は、12個の要素を含む。点群における各点p101は、行列における要素若しくは位置に関連付けるか、又はこの要素若しくは位置に割り当てることができる。この関連付けにより、3D空間における点p101の実際の座標とは独立とすることができる、点間の編成が生成される。したがって、行列又はメモリにおけるこれらの点の編成に基づいて、点に演算を適用することができ、例えば、点は、編成行列又はグリッドにおけるこれらの点の近接性に基づいて演算することができる。点群における点p、i={1,2,...,N}101の数がメモリアドレス又は行列要素の総数よりも小さい場合、メモリアドレス又は行列要素のうちのいくつかは、これらに関連付けられる入力点p101を有しない。これらの要素は、空白の空間、又は編成グリッド内の穴である「ヌル点」と呼ぶことができる。
編成されていない点群において、入力点p101に関連付けられる、上述されたような編成行列又はグリッドは存在しない。3D空間における各点のロケーションは既知であるが、行列若しくはグリッドにおける要素のロケーション又はメモリにおけるアドレスロケーションとの関連付けによって規定される点間に関連性は存在しない。
編成された点群についての対応する点の計算
図2は、対応する点の計算プロセス200のブロック図を示している。N個の入力点p、i={1,2,...,N}101のセットを所与とすると、入力点pは、予め設定された編成又はマッピングに従って、メモリ内又は行列内に配置される。
符号化プロセス100において、上述されたように、予め設定された編成又はマッピングは、編成情報117によって表すことができる。編成情報117は、指定されたマッピングとすることもできるし、各入力点p101がメモリに入力又は記憶された順序に基づいて推測することもできる。N個の入力点p101は、M個の要素、ロケーション、又はアドレスを有するメモリ又は行列若しくはベクトルに配置されて記憶される(201)。M≧Nの場合、「ヌル点」のプレースホルダーがメモリ、行列、又はベクトルのロケーションに記憶される(202)。これらのヌル点のプレースホルダーは、どの入力点p101にも関連付けられていない。この場合、配置された入力点p101及びヌル点のプレースホルダーは、組み合わされてM個のロケーションに記憶され(210)、このM個のロケーションは、N個のp101及びM−N個のヌル点のプレースホルダーを含む。
編成された点群において、ステップS1に示すように、各入力点p101は、メモリ又は1つ以上の要素を含む行列若しくはベクトルにおける位置又はアドレスに関連付けられる。編成された点群における各入力点p101は、行列における予め設定された位置に割り当てられた要素に関連付けることができる。この場合、関連付けは、入力点p101と行列における要素との間でなされる。この関連付けは、3D空間における入力点p101の実際の座標とは独立して実行することができる。ステップS2において、配置オペレーションを、この関連付けに基づいて行列における要素に適用することができる。ステップS2において、行列における要素は、行列におけるこれらの要素の近接性に基づいてM個のメモリロケーションに配置することができる。
点群における入力点p、i={1,2,...,N}101の数がM個のメモリアドレス又は行列要素の総数よりも小さい場合、メモリアドレス又は行列要素のうちのいくつかは、入力点p101を有しない。これらの入力点pを有しない要素は、空白の空間、又は編成行列における穴である「ヌル点」と呼ぶことができる。ステップS1及びS2が実行された後、行列要素は、M個のメモリアドレス又はメモリにおける行列要素に記憶される。例えば、メモリに12個の要素が記憶される行列が、図2に示されている。ヌル点が
Figure 2019521417
によって示されていることに留意されたい。すなわち、M≧Nの場合、「ヌル点」のプレースホルダーが、どの入力点p101にも関連付けられていない、メモリ、行列又はベクトルのロケーションに記憶される(202)。この場合、配置された入力点p101及びヌル点のプレースホルダーは、組み合わされてM個のロケーションに記憶され(210)、このM個のロケーションは、N個の入力点p101及びM−N個のヌル点のプレースホルダーを含む。
ステップS4において、N個の入力点p101が、メモリ、M個のアドレス又は要素を有する行列又はベクトルに編成されると、表面モデル点104、例えば、離散的表面モデル点fm、j={1,2,...,M}106が、ステップS4において表面モデルフィッティングプロセス102を用いてプロセス203において生成される。
ステップS5において、組み合わされた(210)ロケーションに含まれるものと同じ数及び配置のアドレス又は要素を含む表面モデル点fmが、メモリに配置され(204)、その結果、表面モデル点211の順序付けられた又は編成された配置がもたらされる。
いくつかの場合、表面モデル点は、フィッティングプロセス102中に生成することができる。
プロセス109において表面モデル上の対応する点を計算するプロセスは、N個の対応する表面モデル点f、i={1,2,...,N}110を出力するので、モデルと編成された点群の点ロケーションとの間の対応する点を計算するプロセス200の最終ステップ(ステップ6)は、M個の表面モデル点fm、j={1,2,...,M}106からN個の点を選択する(205)方法である。この選択は、組み合わされた(210)ロケーションのコンテンツを、例えば、最初の要素から最後の要素までの順で検査することによって達成することができる。この順序における位置は、jによってインデックス付けすることができ、非ヌル点を追跡するインデックスは、iとして示すことができる。組み合わされた(210)ロケーションの第1の要素、すなわち、j=1が非ヌル点を含む場合、この非ヌル点は、pとなり、そのため、対応する点fは、表面モデル点fmに等しく設定され、iは増分される。この第1の要素がヌルである場合、このロケーションについて対応する点は存在せず、そのため、割り当ては行われない。次に、jが増分され、これらのステップは、組み合わされた(210)ロケーションの次の要素について反復される。最終結果は、表面モデル点fm、j={1,2,...,M}の要素から選択又はマッピングされる(212)N個の対応する点f、i={1,2,...,N}110の出力であり、ここで、この表面モデル点のロケーションは、組み合わされた(210)ロケーションを記憶するメモリにおける非ヌルエントリに対応する。
上述したように、復号プロセスにおいて、点群は、編成された点群とすることができる。別の場合では、点群は、編成されていない点群とすることができる。
編成されていない点群についての対応する点の計算
図3は、対応する点を計算するプロセス300のブロック図を示している。編成されていない点群の表面モデル上の対応する点を計算するプロセス300が図3に示されている。N個の入力点p、i={1,2,...,N}101のセットを所与とすると、プロセス300が編成されていない点群上で行われ、そのため、入力点p、i={1,2,...,N}101は、予め設定された編成又はマッピングに従ってメモリ又は行列内に配置される必要がない。
表面モデルを、3D空間におけるロケーション又は座標のパラメトリック表現であるp(u,v)によって表すことができることが上記で記載された。編成されていない点群及び任意選択で編成された点群について、表面モデルフィッティングプロセス102は、所与の入力点に対応する表面モデルパッチ上のロケーションを生成するのに用いられるパラメーター(u,v)を計算することができる。したがって、入力点p、i={1,2,...,N}101の各々について、パラメーターの対(u,v)301が存在する。パラメーター対(u,v)301のセットは、付加的に、モデルパラメーター103としてシグナリングすることができる。各入力点p101に関連付けられる各パラメーター対(u,v)301の選択は、表面モデルフィッティングプロセス102中、例えば、各入力点p101と表面パッチp(u,v)上のロケーションとの間の差分又は歪みを最小化することによって計算することができる。
入力点p101及び計算されたパラメーター(u,v)301を所与とすると、表面モデルfm、i={1,2,...,N}311上の対応する点(110)は、所与の(u,v)301のパラメーター対について、表面モデルfm311上の対応する点を、表面パッチp(u,v)上のロケーションに設定することによって計算される(303)。
N個の入力点p、i={1,2,...,N}101が、表面モデル点fm、i={1,2,...,N}311とともに順序付けられて配置され(310)、この順序は、パラメーター対(u,v)301の順序付けられたリスト320の順序と一致する。入力点101の数と表面モデル点fmの数とは双方ともにNであるので、入力点pとモデル点fmとの間に一対一の対応関係が存在し、そのため、表面モデル104の対応する点f、i={1,2,...,N}110は、プロセス305において、表面モデル点fm、i={1,2,...,N}311への一対一のマッピング312を用いて設定される。
点ロケーションと表面モデルの対応する点との間の残差を計算するプロセス
編成された点群の点ロケーションと表面モデル104の対応する点との間の残差を計算するプロセスは、図4に示されている。図2の配置プロセス201及び記憶プロセス202からの出力として入力点p、i={1,2,...,N}101及びヌル点のプレースホルダーを含む、組み合わされた(210)ロケーションのセットを所与とし、かつ図1の符号化プロセス100中のプロセス109において計算された表面モデルの対応する点f、i={1,2,...,N}110を所与とすると、残差r、i={1,2,...,N}108が、各入力点p101と表面モデル104の対応する点fとの間の差分として計算される(440)。この差分は、点の各座標について計算することができ、ここで、r={x ,y−f ,z−f }、i={1,2,...,N}である。編成されていない点群を処理する場合、配置された点210にヌル点のプレースホルダーは存在しないので、そのため、残差のリストは、入力点101と表面モデル104上の対応する点110とに対して一対一の対応関係を有する。編成された点群を処理する場合、残差を計算することにおいて一対一の対応関係が依然として存在するが、残差は、配置された点210と同じように配置される。図4において、配置された点210のリストにおけるヌル点によって占有されるロケーションが、配置された点210のリストにおけるヌル点ロケーションに対応するエントリを有しないように、対応する点のリスト430がより大きなメモリ内で配置される。配置された点210のリストから対応する点のリスト430のこのリストを減じることによって、残差441を記憶するメモリは、配置された点210を記憶するメモリ内のロケーションと同じロケーションにおいてヌル点エントリを含むことができる。したがって、編成された点群について、残差計算プロセス107から出力された残差108は、記憶ロケーションに基づいて明示的又は暗示的に、ヌル点エントリのロケーションに関する情報も含む。編成された点群について、ヌル点エントリは存在せず、そのため、出力残差108は、残差計算プロセス107によって計算されるものとして、差分のリストを含む。
一実施形態では、残差データは、対応する点と入力点との間の距離を表す。
さらに、いくつかの場合、入力点p101は、属性を含むことができる。この場合、対応する点は、属性を含むことができる。例えば、属性は、色情報を含むことができる。
さらに、いくつかの実施形態は、シーンを表す点群を符号化する符号化器システムを開示し、点群の各点は、3次元(3D)空間におけるロケーションであり、本符号化器システムは、メモリと通信するプロセッサと、メモリに記憶される符号化器モジュールであって、ステップを実行することによってシーンを表す点群を符号化するように構成される、符号化器モジュールとを備え、ステップは、パラメーター化された表面を、入力点によって形成される点群上にフィッティングするステップと、パラメーター化された表面からモデルパラメーターを生成するステップと、パラメーター化された表面上の対応する点を計算するステップであって、この対応する点は、入力点に対応する、ステップと、対応する点と点群の入力点とに基づいて残差データを計算するステップと、モデルパラメーター及び残差データを圧縮して、それぞれ、符号化されたモデルパラメーター及び符号化された残差データを得るステップと、パラメーター化された表面の符号化されたモデルパラメーターと符号化された残差データとからビットストリームを生成するステップとを含む。
復号プロセス
本発明のいくつかの実施形態は、パラメーター化された表面のモデルパラメーターと元の点群から計算された残差データとを含むビットストリームを受信することによって、符号化された点群データを効果的に復号することができるという認識に基づいている。いくつかの実施形態は、メモリと通信するプロセッサを備える復号器を用いて、シーンを表す点群を復号する方法を開示し、本方法は、パラメーター化された表面についてのモデルパラメーターを受信するステップと、残差データを受信するステップと、モデルパラメーターを用いてパラメーター化された表面を求めるステップと、予め設定された配置に従ってパラメーター化された表面から対応する点を計算するステップと、残差データと対応する点とを合成することによって再構成された入力点を計算するステップとを含む。
図5は、本発明の一実施形態の復号器又は復号プロセスの図である。復号器への入力はビットストリーム116である。ビットストリーム116は、エントロピー復号器515によって復号され、これにより、モデルパラメーター103のセットと量子化された変換係数114とが生成される。量子化された変換係数は、逆量子化され(513)、逆量子化された変換係数512が生成される。逆量子化された変換係数512は逆変換511に入力され、これにより、再構成された残差
Figure 2019521417
508のセットが生成され、この残差は、例えば、デカルト座標系を用いる場合、
Figure 2019521417
とすることができる。
エントロピー復号器515から復号されたモデルパラメーター103は、表面モデル生成プロセス515に入力され、これにより、符号化プロセス100の表面モデルフィッティングプロセス102の方法と同じ方法において表面モデル104が生成されるが、これには、ビットストリーム116からモデルパラメーターを復号することによってこれらが提供されるので、このプロセスのフィッティング部分を伴わない。表面モデル104は、連続的なものである
Figure 2019521417
105とすることもできるし、離散的なものである
Figure 2019521417
106とすることもできる。表面モデル104がベジェ曲面パッチである場合、モデルパラメーター103は、制御点bijを含む。表面モデル104は、プロセス109において表面モデル上の対応する点を計算するプロセスに入力され、これにより、予め設定された配置に従って表面モデル104の対応する点f、i={1,2,...,N}110が計算される。例えば、予め設定された配置は、モデルにおけるパラメーターの近接性によって実行することができる。この場合、パラメーター化された表面からの対応する点は、或る配置に従って計算することができ、モデルパラメーターは、配置の明細を含む。
再構成された残差
Figure 2019521417
と表面モデル104上の対応する点fとが入力として与えられると、合成器507が、再構成された点群又は再構成された入力点
Figure 2019521417
501を計算する。合成器507において実行されるプロセスは、例えば、加算とすることができ、この加算において、
Figure 2019521417
である。
別の例では、合成器507によるプロセスは、残差データ及び対応する点を加算することによって実行することができる。
符号化プロセス100における入力点p101と同様に、再構成された点群
Figure 2019521417
501は、編成されていないものとすることも、編成されたものとすることもできる。
編成されていない場合について、ビットストリーム116から復号された(515)モデルパラメーターエントロピーは、(u,v)301のパラメーター対のリスト320を含むことができ、このパラメーター対は、ベジェ曲面モデル104の表面パッチp(u,v)等の、表面モデル104を生成するためのパラメーターとして用いられる。表面モデル104上の対応する点は、これらの点が符号化器において計算される方法と同様に計算され(109)、表面モデル104上のN個の対応する点110のセットがもたらされる。
編成される場合について、再構成された残差
Figure 2019521417
508は、ヌル点のロケーションを含み、すなわち、これらの残差は、符号化プロセス100内で生成された残差441と同じ方法でメモリに配置される。復号プロセス500において表面モデル104上の対応する点110を計算するプロセス109は、符号化プロセス100内で行われる方法と同じ方法で行うことができる。
ヌル点を有する再構成された残差508上で行われる場合、合成器507は、入力ヌル点をその出力にコピーする。
別の実施形態では、合成器は、表面モデル104上の対応する点110をその出力にコピーすることができ、この場合、再構成された点群501におけるヌル点のプレースホルダーの全てが表面モデル104からの対応する点110に置き換えられる。この実施形態は、「穴」すなわち点群から欠落した点を、表面モデル104からの点に置き換えるように用いることができる。
階層的分割プロセス
2つ以上の表面モデル104が必要な場合、階層的分割プロセスを用いて点群を分解して、各々がそれら自体のモデルを有するより小さな点群にすることができる。
符号化の階層的分割プロセス
編成された点群について、表面モデル104は、表面パッチp(u,v)とすることができ、ここで、パラメーターの対(u,v)301が、表面モデル104のために3D空間における点ロケーションを計算するのに用いられる。連続的表面モデル104について、表面パッチp(u,v)のパラメーター(u,v)は、単位正方形上に拡がり、すなわち、0≦u、v≦1である。離散的表面モデルについて、パラメーター空間(u,v)は、パラメーター値(u,v)、j={1,2,...,J}、k={1,2,...,K}のJ×Kグリッドとして離散化することができる。この場合、表面パッチp(u,v)は、3D空間においてM×N点を含む。
図6は、本発明の実施形態による階層的分割プロセスの図である。初期パラメーター空間、この場合、2Dグリッド(u,v)、j={1,2,...,J}、k={1,2,...,K}601は、幅w=J及び高さh=Kを有する初期矩形602によって表すことができる。
いくつかの場合、近接性を、2Dグリッド上の近傍位置によって規定することができる。さらに、近傍位置は、水平方向、垂直方向、及び対角方向におけるものとすることができる。
サイズw×hを有する矩形602によって表される初期パラメーター空間(u,v)601を所与とすると、表面モデルパッチp(u,v)603は、表面モデルフィッティングプロセス102によって生成される。符号化プロセス100及び復号プロセス500は、fm、l={1,2,...,L}106として離散的表面モデルを表し、そのため、表面モデルマッピングに対するパッチ604は、p(u,v)j={1,2,...,J}、k={1,2,...,K}をfm、l={1,2,...,L}にマッピングする。したがって、L=J×Kである。
全ての入力点を始めに表すことができる、現在の矩形表現に関連付けられた入力点と、離散的表面モデル104、106との間のフィッティング誤差e603が計算される(605)。フィッティング誤差eは、例えば、入力点の各成分(例えば、x成分、y成分、及びz成分)とその対応する表面モデル点との間の全誤差又は平均二乗誤差として測定することができる。別の実施形態では、フィッティング誤差eは、表面モデル104の表面法線及び事前に計算された表面モデル104の表面法線の全偏差又は平均偏差として測定することができる。
フィッティング誤差e606が予め設定されたフィッティング誤差閾値T607よりも小さい場合、階層的分割プロセス600は、成功したものとみなされ、そのため、矩形がこれ以上分割されないことを示す分割フラグ610が出力され(608)、モデルパラメーター103とともに現在の矩形に関連付けられた点についての離散的表面モデルfm106も出力され(608)、このプロセスは終了する。表面モデル104がベジェ曲面である場合、モデルパラメーター103は、式(1)に示される制御点bijを含むことができる。これらのモデルパラメーターは、現在の矩形表現の幅w及び高さhも含むことができる。
フィッティング誤差e606が予め設定されたフィッティング誤差閾値T607以上である場合、現在の表面モデル104は、現在の矩形に関連付けられた入力点への良好なフィットではないものとみなされる。この場合、矩形は、矩形分割プロセス609によって分割されて2つ以上の矩形にされ、矩形が分割されることを示す分割フラグ610が出力される(611)。
矩形分割プロセス609は、パラメーター空間の現在の矩形表現を取り、幅w及び高さhを有する矩形を分割して、2つ以上の矩形を含む分割された矩形612にする。例えば、バイナリ分割は、矩形を、各々幅w/2及び高さhを有する2つの矩形に分割することができる。別の例では、分割後、各矩形は、幅w及び高さh/2を有することができる。w/2又はh/2が整数ではない場合、丸め技法を行うことができ、例えば、一方の矩形は、幅floor(w/2)を有することができ、他方はfloor(w/2)+1を有することができ、ここでfloor()は、端数を切り捨てて最も近い整数にする。バイナリ分割について、幅に沿って分割するか又は高さに沿って分割するかは、予め設定されたプロセスによって判断することもできるし、いずれの次元を分割するのかについての判断は、ビットストリームにおいてフラグとして明示的にシグナリングすることもできる。予め設定されたプロセスによって判断する場合、そのような予め設定されたプロセスの一例は、より長い方の次元に交差するように矩形を分割することである。例えば、w=10及びh=4である場合、矩形は、幅に交差するように分割することができ、そのため、分割された矩形612は、各々幅w=5及び高さh=4を有する。
別の実施形態では、矩形分割プロセス609は、1つの矩形を分割して3つ以上の矩形、例えば、各々矩形の1つの4分割部分を表す4つの矩形にすることができる。
別の実施形態では、矩形分割プロセス609は、矩形に関連付けられた入力点101の密度に基づいて矩形を分割することができる。例えば、その矩形についての入力点が3D空間における2つの別個のオブジェクトを表す場合、矩形分割プロセス609は、この矩形を分割して、各々が1つのオブジェクトを含む2つの部分にすることができる。矩形を最も良好に分割する場所を求めるために、例えば、矩形についての入力点の密度を測定することができ、分割は、各分割された矩形内の密度を最大化するラインに沿って行うことができ、又は、別の例では、分割は、ラインと分割された矩形内の各点との間の距離の総和を、例えば、最小二乗近似を用いてラインをフィッティングすることによって、最小化するラインに沿って行うことができる。分割のロケーションを示すパラメーターを、モデルパラメーターとして符号化してビットストリームにおいてシグナリングすることができる。
分割された矩形612が取得されると、階層的分割プロセス600は、分割された矩形612を、階層的分割プロセス600内の表面モデルフィッティングプロセス102に入力することによって、これらの分割された矩形612の各々について反復される。
矩形又は分割された矩形のフィッティング誤差606が閾値607よりも小さい場合にこの矩形又は分割された矩形について分割プロセス609を終端することに加えて、矩形についての分割プロセス609は、幅w又は高さh、又は面積、すなわち、幅と高さとの積が予め設定された値よりも小さい場合に終端することができる。例えば、16個の制御点を有するベジェ曲面パッチが表面モデル104として用いられる場合、及び、矩形の面積が10である場合、16個の制御点のシグナリングを必要とする表面モデル104を用いて矩形をフィッティングする代わりに、矩形に関連付けられた10個の入力点101を直接シグナリングするのがより効率的であり得る。
別の実施形態では、表面モデル104に用いる制御点の数は、矩形の幅と高さとに依拠することができる。例えば、矩形の面積が16よりも小さい場合、より少数の制御点を有する表面モデル104を用いることができる。ルックアップテーブルに対する値又はインデックスを、所与の矩形のためにいくつの制御点が表面モデル104において用いられるかを示すのに、ビットストリームにおいてシグナリングすることができる。
別の実施形態では、ランク、すなわち、入力点101の行列の線形独立な行若しくは列の数、又は矩形に関連付けられた入力点の行列の分解を測定することができ、表面モデル104を生成するのに用いる制御点の数は、ランク以下の値に設定することができる。
階層的分割プロセス600が完了すると、分割フラグ610のシーケンスが出力される。例えば、バイナリ分割がより長い次元に交差するように常に生じることを認識することによって、分割フラグと初期矩形表現の幅及び高さとによって完全に指定することができるという点において分割プロセスが予め設定されたものである場合、全ての分割された矩形のロケーションを復元するのに十分な情報が復号器に利用可能である。したがって、初期矩形表現602の幅及び高さ、分割フラグ610のシーケンス、並びに各矩形についての制御点等のモデルパラメーター103を復号器が用いて、各矩形についての表面モデル104を生成することができる(515)。
復号器における分割された矩形の処理
1つの実施形態では、復号器は、ビットストリームから、初期矩形表現602の幅及び高さ、すなわち、2D編成グリッドを復号することができる。次に、復号器は、ビットストリームから、分割フラグ610を復号し、ここで、1すなわち真の分割フラグは、矩形が分割されて、例えば、2つの矩形(第1の矩形及び第2の矩形)にされることを示し、そして例えば、その分割は、矩形のより大きい次元に交差するように生じることができる。矩形が分割される場合、ビットストリームから復号された次の分割フラグ610は、第1の矩形についての分割フラグ610である。
その分割フラグが0すなわち偽である場合、第1の矩形は、この後分割されず、ペイロードフラグをビットストリームから復号して、いずれの種類のデータが第1の矩形について復号されるのかを示すことができる。ペイロードフラグが1すなわち真である場合、制御点、及び第1の矩形についての量子化された変換係数114等の残差108を表すデータが、ビットストリームからエントロピー復号される。残差を表すデータを復号した後、この矩形についての表面モデル104の制御点等のモデルパラメーター103をビットストリームからエントロピー復号することができる。ペイロードフラグが0すなわち偽である場合、第1の矩形について残差は利用可能でなく、これは、例えば表面モデル104が用いられず、符号化器が入力点101をビットストリームに直接シグナリングした場合に起こり得る。この場合、復号器は、次に、ビットストリームから、第1の矩形についての入力点、量子化された入力点、又は変換された入力点の量子化された変換係数をエントロピー復号する。別の実施形態では、ペイロードフラグは用いられない。この場合、表面モデル104が常に用いられ、そのため、残差を表すデータがビットストリームから復号される。
その分割フラグが1すなわち真である場合、第1の矩形は、更に分割される。
階層又はツリーを表すデータ構造は、幅優先又は深さ優先でトラバースすることができる。幅優先の場合、ビットストリームから復号される次のデータは、第2の矩形についての分割フラグである。深さ優先の場合、ビットストリームから復号される次のデータは、現在の矩形についての分割フラグであり、この場合、この分割フラグは、第1の矩形が更に分割されるか否かを示す。
復号器における分割された矩形のこの処理は、全ての矩形が処理されるまで各矩形上で実行される。
符号化プロセス中の場合には、付加的な基準を用いて或るブロックが分割されるか又は分割されないかを判断することができる。例えば、矩形の次元(幅及び/又は高さ、又は面積)が予め設定された閾値未満である場合、その矩形についての分割プロセスは、ビットストリームから分割フラグを復号する必要なく終端することができる。各矩形の次元は、初期矩形表現602の高さ及び幅から分割プロセス中に推測することができる。
付加的な実施形態
別の実施形態では、矩形又は分割された矩形612の選択されたセットが処理され、再構成された点群501が復号プロセス600によって出力された後、例えば、スケーリング可能な復号システムにおいて、付加的な矩形を処理して、再構成された点群501のための付加的な点を生成することができる。
別の実施形態では、表面モデルが符号化プロセス100及び復号プロセス600によって用いられる場合にもたらされる、離散的表面モデル106と再構成された点群501との間の差分としてフィッティング誤差606が計算される。
別の実施形態では、差分に基づくフィッティング誤差606、例えば、入力点101と離散的表面モデル106との間の平均二乗誤差を計算する代わりに、誤差メトリックは、処理される矩形に関連付けられた全てのデータを表すのに、平均二乗誤差とビットストリームにおいて占有されるビットのスケーリングされた数との合計とすることができる。
別の実施形態では、矩形の分割は、矩形のより短い次元に交差するように生じさせることができる。
別の実施形態では、制御点、他のモデルパラメーター103、及び残差108に関連付けられたデータ等の、各矩形に関連付けられたどのデータよりも前に、全ての分割フラグ610がビットストリームから復号される。この実施形態により、残りのデータがビットストリームから復号される前に、完全な分割階層が復号器によって認識されることが可能になる。
別の実施形態では、ペイロードフラグと関連付けられたデータとが復号される前に、分割フラグのうちのいくつか又は全てがビットストリームから復号され、その後、ペイロードフラグと矩形の選択されたサブセットからのデータとをビットストリームから復号することができる。初期矩形表現602上の対象の領域を指定する、例えば、コンピューターディスプレイ上の編成グリッドの表現のアウトラインを引くことによって、所望の矩形を選択することができ、その後、対象の領域に含まれるか又は部分的に含まれる矩形に関連付けられたビットストリーム内のデータのみをビットストリームから復号することができる。
別の実施形態では、3D空間から対象の領域が選択され、例えば、3D境界ボックスが対象の入力点101を包含し、その後、3Dの対象の領域における点ロケーションについて、2D矩形表現602におけるそれらの対応するロケーションが識別され、すなわち、検索又は逆マッピングが行われ、2D境界ボックスが、全てのこれらの選択された対応するロケーションを含む2D矩形表現上で計算される。この境界ボックスは、選択された対応する入力点101によってポピュレートされた新たな最初の部分矩形を含む。3D境界ボックスの内部に対応する3D点ロケーションを有しないこの新たな部分矩形の内部の任意のロケーションについて、これらのロケーションはヌル点でポピュレートされ、ここで、ヌル点は、3D空間において対応する点がない、例えば、「穴」であることを示す値である。
別の実施形態では、矩形分割プロセス609は、各分割された矩形612に関連付けられた入力点101の数が等しくなるように又は概ね等しくなるように矩形を分割する。
別の実施形態では、各入力点101は、関連付けられた属性を有する。属性は、限定するものではないが、色値、反射、又は温度を含む追加データである。
別の実施形態では、表面モデルパッチの制御点は、符号化プロセス100中に量子化され、任意選択で変換され、復号プロセス500中に逆量子化され、任意選択で逆変換される。
別の実施形態では、入力点を含まないメモリロケーションにおけるヌル点のプレースホルダーとともに入力点がメモリ内に配置されない、編成されていない点群について、パラメーター化されたモデルに関連付けられたマニホールド上の対応する入力点ロケーションの位置は、ビットストリームにおいてシグナリングされる。例えば、ベジェ曲面又はベジェパッチを用いて、パラメーターu及びvによって表される単位正方形が3D空間における(x,y,z)座標にマッピングされ、ここで、u及びvの双方は、0.0〜1.0(0.0及び1.0を含む)である。好ましい実施形態の場合である編成された点群について、単位正方形は、3Dにおける各(x,y,z)点が(u,v)単位正方形平面上のサンプル位置に対応し、3D空間において対応する点を有しない(u,v)単位正方形平面上のサンプル位置はヌル点によってポピュレートすることができるように、均一なグリッドを用いてサンプリングすることができる。(u,v)単位正方形平面におけるサンプル位置は、予め設定された順序を用いて、入力点101によってポピュレートする、すなわち、関連付けることができる。例えば、入力点101のシーケンスは、2Dサンプル位置を、2Dサンプル位置における各行が、その行内の最初の要素から最後の要素まで埋められて、その後、次の入力点が、次の行へと向かうラスタスキャン順でポピュレートすることができる。
別の実施形態では、パラメーター空間601は、均一なサンプリングを有する2Dグリッド又はマニホールドである。
別の実施形態では、パラメーター空間601は、不均一なサンプリングを有する2Dグリッド又はマニホールドである。例えば、グリッドの中心は、エッジよりも高密度で点を有することができる。更なる一例として、均一なグリッドは、全ての整数の位置において
サンプリングすることができ、その後、パラメーター空間全体の半分の幅及び高さを有する、中心を合わせた正方形を、全ての整数の2分の1又は4分の1の位置において付加的にサンプリングすることができる。
別の実施形態では、2D変換111が、残差データ108の成分に、2D編成グリッド117上のこれらの対応するロケーションに従って、適用される。
別の実施形態では、分割された矩形112に対応する残差データ108は、分割された矩形612の次元に基づいた順序、例えば、最大面積から最小面積の順序、又は別の例では、最小面積から最大面積の順序でビットストリーム116にシグナリングされ、ここで面積は、分割された矩形612の面積すなわち幅×高さである。
別の実施形態では、再構成された入力点は、3次元マップ(3Dマップ)を含む。或る車両は、3D空間におけるその位置を、この車両上に配置されたセンサーから点群を捕捉することによって求め、その後、捕捉された点群を、再構成された入力点又は再構成された点群と比較することができる。捕捉された点群上の点を、再構成された入力点と位置合わせすること、すなわち、アライメントすることによって、捕捉された点群におけるオブジェクト又は点の位置を、再構成された点群におけるオブジェクト又は点に関連付けることができる。車両のセンサーによって捕捉されたオブジェクトの、車両又は車両のセンサーの位置に対する位置が既知であることを考慮し、及び、位置合わせ処理の後、車両のセンサーによって捕捉されたオブジェクトの、再構成された点群に対する位置が既知であることを考慮すると、再構成された点群における車両の位置を推測することができ、したがって、3Dマップにおける車両の位置を推測、ひいては認識することができる。
発明の効果
本発明のいくつかの実施形態によれば、点群を効果的に符号化及び復号することができ、実施形態は、点群におけるオブジェクトの3次元表現を圧縮するのに有用とすることができる。さらに、本発明のいくつかの実施形態による符号化及び復号の方法は、ファイル又はビットストリーム全体を復号する必要なく、点群ジオメトリの粗い表現を迅速かつ容易に復号及び再構成することを可能にする圧縮表現を生成することができる。

Claims (20)

  1. メモリと通信するプロセッサを備える符号化器を用いて、シーンを表す点群を符号化する方法であって、前記点群の各点は、3次元(3D)空間におけるロケーションであり、該方法は、
    パラメーター化された表面を、入力点によって形成される前記点群上にフィッティングするステップと、
    前記パラメーター化された表面からモデルパラメーターを生成するステップと、
    前記パラメーター化された表面上の対応する点を計算するステップであって、該対応する点は、前記入力点に対応する、ステップと、
    前記対応する点と前記点群の前記入力点とに基づいて残差データを計算するステップと、
    前記モデルパラメーター及び前記残差データを圧縮して、符号化されたモデルパラメーター及び符号化された残差データを得るステップと、
    前記パラメーター化された表面の前記符号化されたモデルパラメーターと前記符号化された残差データとからビットストリームを生成するステップと、
    を含む、方法。
  2. 前記点群は、編成された点群である、請求項1に記載の方法。
  3. 前記点群は、編成されていない点群である、請求項1に記載の方法。
  4. 前記残差データは、前記対応する点と前記入力点との間の距離を表す、請求項1に記載の方法。
  5. 前記入力点及び前記対応する点は、属性を含む、請求項1に記載の方法。
  6. 前記属性は、色情報を含む、請求項5に記載の方法。
  7. 前記圧縮するステップは、エントロピー符号化を含む、請求項1に記載の方法。
  8. 前記圧縮するステップは、変換ステップ及び量子化ステップを含む、請求項1に記載の方法。
  9. 前記モデルパラメーターを前記生成するステップは、前記フィッティングするステップ中に実行される、請求項1に記載の方法。
  10. 前記残差データを前記計算するステップは、前記対応する点の前記成分と前記点群の前記入力点の前記対応する成分との間の差分である、請求項1に記載の方法。
  11. 前記パラメーター化された表面についての前記モデルパラメーターを受信するステップと、
    前記残差データを受信するステップと、
    前記モデルパラメーターを用いて前記パラメーター化された表面を求めるステップと、
    予め設定された配置に従って前記パラメーター化された表面から前記対応する点を計算するステップと、
    前記残差データと前記対応する点とを合成することによって再構成された入力点を計算するステップと、
    を更に含む、請求項1に記載の方法。
  12. 前記パラメーター化された表面からの前記対応する点は、或る配置に従って計算され、前記モデルパラメーターは、前記配置の明細を含む、請求項11に記載の方法。
  13. 前記予め設定された配置は、予め設定された配置における近接性によって求められる、請求項11に記載の方法。
  14. 前記合成するステップは、前記残差データと前記対応する点とを加算することによって実行される、請求項11に記載の方法。
  15. 前記再構成された入力点は、3次元マップを含み、車両が、該車両から取得されたデータと前記3次元マップにおけるデータとを登録することによって、前記3次元マップ上での該車両の位置を求める、請求項11に記載の方法。
  16. 前記車両から取得された前記データは、点群であり、前記登録することは、前記車両によって取得された前記点群を、3次元マップを含む前記再構成された入力点と比較することを含む、請求項15に記載の方法。
  17. モデルパラメーターのサブセット及び残差データのサブセットが受信されて、入力点のサブセットを再構成するように用いられ、後続する付加的なモデルパラメーターのサブセット及び付加的な残差データのサブセットが受信されて、再構成された前記入力点のサブセットを精緻化するように用いられる、請求項11に記載の方法。
  18. 2次元グリッドを包含する矩形を画定するステップと、
    各入力点を、前記2次元グリッド上のインデックスに関連付けるステップと、
    前記矩形の各々について、該矩形内にインデックス付けされた前記入力点に、パラメーター化された表面をフィッティングするステップと、
    前記矩形の各々について、前記パラメーター化された表面と前記矩形内にインデックス付けされた前記入力点との間のフィッティング誤差を測定するステップと、
    前記フィッティング誤差が予め設定された閾値を超えた場合、前記矩形の各々を階層的に分割してより小さな矩形にするステップと、
    を更に含む、請求項1に記載の方法。
  19. シーンを表す点群を符号化する符号化器システムであって、前記点群の各点は、3次元(3D)空間におけるロケーションであり、該符号化器システムは、
    メモリと通信するプロセッサと、
    前記メモリに記憶される符号化器モジュールであって、ステップを実行することによってシーンを表す点群を符号化するように構成される、符号化器モジュールと、
    を備え、前記ステップは、
    パラメーター化された表面を、入力点によって形成される前記点群上にフィッティングするステップと、
    前記パラメーター化された表面からモデルパラメーターを生成するステップと、
    前記パラメーター化された表面上の対応する点を計算するステップであって、該対応する点は、前記入力点に対応する、ステップと、
    前記対応する点と前記点群の前記入力点とに基づいて残差データを計算するステップと、
    前記モデルパラメーター及び前記残差データを圧縮して、それぞれ、符号化されたモデルパラメーター及び符号化された残差データを得るステップと、
    前記パラメーター化された表面の符号化された前記モデルパラメーターと前記符号化された残差データとからビットストリームを生成するステップと、
    を含む、符号化器システム。
  20. シーンを表す点群を符号化するプログラムが記憶された非一時的コンピューター可読記録媒体であって、前記点群の各点は、3次元(3D)空間におけるロケーションであり、プロセッサによって実行されると、前記プログラムは、該プロセッサに、
    パラメーター化された表面を、入力点によって形成される前記点群上にフィッティングするステップと、
    前記パラメーター化された表面からモデルパラメーターを生成するステップと、
    前記パラメーター化された表面上の対応する点を計算するステップであって、該対応する点は、前記入力点に対応する、ステップと、
    前記対応する点と前記点群の前記入力点とに基づいて残差データを計算するステップと、
    前記モデルパラメーター及び前記残差データを圧縮して、それぞれ、符号化されたモデルパラメーター及び符号化された残差データを得るステップと、
    前記パラメーター化された表面の前記符号化されたモデルパラメーターと前記符号化された残差データとからビットストリームを生成するステップと、
    を実行させる、非一時的コンピューター可読記録媒体。
JP2018559911A 2016-08-19 2017-08-04 シーンを表す点群を符号化する方法、符号化器システム、及びプログラムを記憶した非一時的コンピューター可読記録媒体 Active JP6676193B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US15/241,112 US20180053324A1 (en) 2016-08-19 2016-08-19 Method for Predictive Coding of Point Cloud Geometries
US15/241,112 2016-08-19
PCT/JP2017/029238 WO2018034253A1 (en) 2016-08-19 2017-08-04 Method, encoder system and non-transitory computer readable recording medium storing thereon program for encoding point cloud of representing scene

Publications (2)

Publication Number Publication Date
JP2019521417A true JP2019521417A (ja) 2019-07-25
JP6676193B2 JP6676193B2 (ja) 2020-04-08

Family

ID=59829427

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018559911A Active JP6676193B2 (ja) 2016-08-19 2017-08-04 シーンを表す点群を符号化する方法、符号化器システム、及びプログラムを記憶した非一時的コンピューター可読記録媒体

Country Status (4)

Country Link
US (1) US20180053324A1 (ja)
EP (1) EP3501005B1 (ja)
JP (1) JP6676193B2 (ja)
WO (1) WO2018034253A1 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021210549A1 (ja) * 2020-04-14 2021-10-21 パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカ 三次元データ符号化方法、三次元データ復号方法、三次元データ符号化装置、及び三次元データ復号装置
JP2022512110A (ja) * 2019-03-01 2022-02-02 テンセント・アメリカ・エルエルシー 点群圧縮の方法および装置
JP2022528526A (ja) * 2020-01-09 2022-06-14 テンセント・アメリカ・エルエルシー 点群圧縮のための変換係数のアルファベットパーティションコーディングの方法及び装置並びにプログラム
JP2023510822A (ja) * 2020-01-10 2023-03-15 上海交通大学 点群データのカプセル化方法及び伝送方法
WO2024062938A1 (ja) * 2022-09-20 2024-03-28 パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカ 復号方法及び復号装置

Families Citing this family (67)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10496336B2 (en) * 2016-11-17 2019-12-03 Google Llc K-D tree encoding for point clouds using deviations
US10430975B2 (en) * 2016-11-17 2019-10-01 Google Llc Advanced k-D tree encoding for point clouds by most significant axis selection
EP3399757A1 (en) * 2017-05-04 2018-11-07 Thomson Licensing Method and apparatus to encode and decode two-dimension point clouds
US11818401B2 (en) 2017-09-14 2023-11-14 Apple Inc. Point cloud geometry compression using octrees and binary arithmetic encoding with adaptive look-up tables
US10897269B2 (en) * 2017-09-14 2021-01-19 Apple Inc. Hierarchical point cloud compression
US10861196B2 (en) 2017-09-14 2020-12-08 Apple Inc. Point cloud compression
US11113845B2 (en) 2017-09-18 2021-09-07 Apple Inc. Point cloud compression using non-cubic projections and masks
US10909725B2 (en) 2017-09-18 2021-02-02 Apple Inc. Point cloud compression
US10607373B2 (en) 2017-11-22 2020-03-31 Apple Inc. Point cloud compression with closed-loop color conversion
US11651523B2 (en) 2017-12-28 2023-05-16 Nokia Technologies Oy Apparatus, a method and a computer program for volumetric video
WO2019142164A1 (en) * 2018-01-19 2019-07-25 Interdigital Vc Holdings, Inc. Processing a point cloud
EP3515068A1 (en) * 2018-01-19 2019-07-24 Thomson Licensing A method and apparatus for encoding and decoding three-dimensional scenes in and from a data stream
CN108492329B (zh) * 2018-03-19 2022-04-12 北京航空航天大学 一种三维重建点云精度和完整度评价方法
US10909726B2 (en) 2018-04-10 2021-02-02 Apple Inc. Point cloud compression
US10939129B2 (en) 2018-04-10 2021-03-02 Apple Inc. Point cloud compression
US10867414B2 (en) 2018-04-10 2020-12-15 Apple Inc. Point cloud attribute transfer algorithm
US10909727B2 (en) 2018-04-10 2021-02-02 Apple Inc. Hierarchical point cloud compression with smoothing
US10964067B2 (en) 2018-04-12 2021-03-30 Samsung Electronics Co., Ltd. Visual quality enhancement of reconstructed point clouds via color smoothing
US10956790B1 (en) * 2018-05-29 2021-03-23 Indico Graphical user interface tool for dataset analysis
KR20210018254A (ko) * 2018-06-06 2021-02-17 파나소닉 인텔렉츄얼 프로퍼티 코포레이션 오브 아메리카 삼차원 데이터 부호화 방법, 삼차원 데이터 복호 방법, 삼차원 데이터 부호화 장치, 및 삼차원 데이터 복호 장치
WO2019244931A1 (ja) * 2018-06-19 2019-12-26 パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカ 三次元データ符号化方法、三次元データ復号方法、三次元データ符号化装置、及び三次元データ復号装置
CN112438047B (zh) * 2018-06-26 2022-08-09 华为技术有限公司 用于点云译码的高级语法设计
US11017566B1 (en) 2018-07-02 2021-05-25 Apple Inc. Point cloud compression with adaptive filtering
US11202098B2 (en) 2018-07-05 2021-12-14 Apple Inc. Point cloud compression with multi-resolution video encoding
US11095908B2 (en) 2018-07-09 2021-08-17 Samsung Electronics Co., Ltd. Point cloud compression using interpolation
CN112470480B (zh) 2018-07-11 2024-06-18 索尼公司 图像处理装置和方法
US11012713B2 (en) 2018-07-12 2021-05-18 Apple Inc. Bit stream structure for compressed point cloud data
JP2021182650A (ja) * 2018-07-20 2021-11-25 ソニーグループ株式会社 画像処理装置および方法
US10891758B2 (en) 2018-07-23 2021-01-12 Google Llc Geometry encoder
CN112514396A (zh) 2018-08-02 2021-03-16 索尼公司 图像处理设备和图像处理方法
CN109166160B (zh) * 2018-09-17 2022-07-01 华侨大学 一种采用图形预测的三维点云压缩方法
CN110958455B (zh) * 2018-09-26 2022-09-23 华为技术有限公司 点云编解码方法、编解码器、编解码装置和存储介质
JP2022003716A (ja) 2018-10-01 2022-01-11 ソニーグループ株式会社 画像処理装置および方法
US11367224B2 (en) 2018-10-02 2022-06-21 Apple Inc. Occupancy map block-to-patch information compression
US10853973B2 (en) 2018-10-03 2020-12-01 Apple Inc. Point cloud compression using fixed-point numbers
US11430155B2 (en) 2018-10-05 2022-08-30 Apple Inc. Quantized depths for projection point cloud compression
JP7434175B2 (ja) 2018-12-07 2024-02-20 パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカ 三次元データ符号化方法、三次元データ復号方法、三次元データ符号化装置、及び三次元データ復号装置
CN118781209A (zh) * 2018-12-14 2024-10-15 交互数字Vc控股公司 用于对空间数据进行程序化着色的方法和装置
US11288843B2 (en) 2019-01-04 2022-03-29 Samsung Electronics Co., Ltd. Lossy compression of point cloud occupancy maps
US11120581B2 (en) * 2019-03-01 2021-09-14 Tencent America LLC Method and apparatus for point cloud compression
US11127166B2 (en) * 2019-03-01 2021-09-21 Tencent America LLC Method and apparatus for enhanced patch boundary identification for point cloud compression
US11272158B2 (en) * 2019-03-01 2022-03-08 Tencent America LLC Method and apparatus for point cloud compression
US11057564B2 (en) 2019-03-28 2021-07-06 Apple Inc. Multiple layer flexure for supporting a moving image sensor
WO2021002657A1 (ko) * 2019-07-04 2021-01-07 엘지전자 주식회사 포인트 클라우드 데이터 송신 장치, 포인트 클라우드 데이터 송신 방법, 포인트 클라우드 데이터 수신 장치 및 포인트 클라우드 데이터 수신 방법
WO2021010200A1 (ja) * 2019-07-17 2021-01-21 ソニー株式会社 情報処理装置および方法
CN110418135B (zh) * 2019-08-05 2022-05-27 北京大学深圳研究生院 一种基于邻居的权重优化的点云帧内预测方法及设备
US11627314B2 (en) 2019-09-27 2023-04-11 Apple Inc. Video-based point cloud compression with non-normative smoothing
US11562507B2 (en) 2019-09-27 2023-01-24 Apple Inc. Point cloud compression using video encoding with time consistent patches
US11538196B2 (en) 2019-10-02 2022-12-27 Apple Inc. Predictive coding for point cloud compression
US11895307B2 (en) 2019-10-04 2024-02-06 Apple Inc. Block-based predictive coding for point cloud compression
US11798196B2 (en) 2020-01-08 2023-10-24 Apple Inc. Video-based point cloud compression with predicted patches
US11625866B2 (en) 2020-01-09 2023-04-11 Apple Inc. Geometry encoding using octrees and predictive trees
WO2021210837A1 (ko) * 2020-04-13 2021-10-21 엘지전자 주식회사 포인트 클라우드 데이터 송신 장치, 포인트 클라우드 데이터 송신 방법, 포인트 클라우드 데이터 수신 장치 및 포인트 클라우드 데이터 수신 방법
WO2021246796A1 (ko) * 2020-06-05 2021-12-09 엘지전자 주식회사 포인트 클라우드 데이터 송신 장치, 포인트 클라우드 데이터 송신 방법, 포인트 클라우드 데이터 수신 장치 및 포인트 클라우드 데이터 수신 방법
US11615557B2 (en) 2020-06-24 2023-03-28 Apple Inc. Point cloud compression using octrees with slicing
US11620768B2 (en) 2020-06-24 2023-04-04 Apple Inc. Point cloud geometry compression using octrees with multiple scan orders
CN114598891B (zh) * 2020-12-07 2023-05-26 腾讯科技(深圳)有限公司 点云数据编码方法、解码方法、点云数据处理方法及装置
US11948338B1 (en) 2021-03-29 2024-04-02 Apple Inc. 3D volumetric content encoding using 2D videos and simplified 3D meshes
EP4071718A1 (en) * 2021-04-09 2022-10-12 Beijing Xiaomi Mobile Software Co., Ltd. Method and apparatus of encoding/decoding point cloud geometry data captured by a spinning sensors head
EP4072145A1 (en) * 2021-04-09 2022-10-12 Beijing Xiaomi Mobile Software Co., Ltd. Method and apparatus of encoding/decoding point cloud geometry data captured by a spinning sensors head
CN115412717B (zh) * 2021-05-26 2024-05-17 荣耀终端有限公司 一种点云方位角信息的预测编解码方法及装置
CN115412713B (zh) * 2021-05-26 2024-06-25 荣耀终端有限公司 一种点云深度信息的预测编解码方法及装置
CN117529912A (zh) * 2021-06-04 2024-02-06 华为技术有限公司 3d地图的编解码方法及装置
CN115474046A (zh) * 2021-06-11 2022-12-13 维沃移动通信有限公司 点云属性信息编码方法、解码方法、装置及相关设备
US12079934B1 (en) 2021-07-27 2024-09-03 Apple Inc. Point cloud re-sampling using hierarchical sphere or disk distributions
US20230298216A1 (en) * 2022-03-15 2023-09-21 Tencent America LLC Predictive coding of boundary geometry information for mesh compression
CN118246797B (zh) * 2024-03-21 2024-09-06 苏州奥特兰恩自动化设备有限公司 一种基于人工智能的工厂控制方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7907784B2 (en) * 2007-07-09 2011-03-15 The United States Of America As Represented By The Secretary Of The Commerce Selectively lossy, lossless, and/or error robust data compression method

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2022512110A (ja) * 2019-03-01 2022-02-02 テンセント・アメリカ・エルエルシー 点群圧縮の方法および装置
JP7141535B2 (ja) 2019-03-01 2022-09-22 テンセント・アメリカ・エルエルシー 点群圧縮の方法および装置
JP2022528526A (ja) * 2020-01-09 2022-06-14 テンセント・アメリカ・エルエルシー 点群圧縮のための変換係数のアルファベットパーティションコーディングの方法及び装置並びにプログラム
JP7246515B2 (ja) 2020-01-09 2023-03-27 テンセント・アメリカ・エルエルシー 点群圧縮のための変換係数のアルファベットパーティションコーディングの方法及び装置並びにプログラム
JP2023510822A (ja) * 2020-01-10 2023-03-15 上海交通大学 点群データのカプセル化方法及び伝送方法
WO2021210549A1 (ja) * 2020-04-14 2021-10-21 パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカ 三次元データ符号化方法、三次元データ復号方法、三次元データ符号化装置、及び三次元データ復号装置
WO2024062938A1 (ja) * 2022-09-20 2024-03-28 パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカ 復号方法及び復号装置

Also Published As

Publication number Publication date
US20180053324A1 (en) 2018-02-22
JP6676193B2 (ja) 2020-04-08
EP3501005B1 (en) 2024-10-02
WO2018034253A1 (en) 2018-02-22
EP3501005A1 (en) 2019-06-26

Similar Documents

Publication Publication Date Title
JP6676193B2 (ja) シーンを表す点群を符号化する方法、符号化器システム、及びプログラムを記憶した非一時的コンピューター可読記録媒体
CN113678466A (zh) 用于预测点云属性编码的方法和设备
JP7330306B2 (ja) 変換方法、逆変換方法、エンコーダ、デコーダ及び記憶媒体
US20220292730A1 (en) Method and apparatus for haar-based point cloud coding
US11711535B2 (en) Video-based point cloud compression model to world signaling information
JP7520223B2 (ja) 点群階層化方法、デコーダ、エンコーダ及び記憶媒体
CN114820831A (zh) 划分方法、编码器、解码器以及计算机存储介质
JP7564348B2 (ja) 2次元正則化平面投影に基づくポイントクラウド符号化及び復号方法並びに装置
CN118614061A (zh) 编解码方法、编码器、解码器以及存储介质
KR20120085134A (ko) 3차원 메쉬 부호화 장치 및 방법, 그리고 3차원 메쉬 복호화 장치 및 방법
WO2022131948A1 (en) Devices and methods for sequential coding for point cloud compression
CN112995758B (zh) 点云数据的编码方法、解码方法、存储介质及设备
CN114009014A (zh) 颜色分量预测方法、编码器、解码器及计算机存储介质
Wei et al. Enhanced intra prediction scheme in point cloud attribute compression
CN114915793B (zh) 基于二维规则化平面投影的点云编解码方法及装置
CN114915792B (zh) 基于二维规则化平面投影的点云编解码方法及装置
KR20230119690A (ko) 대규모 포인트 클라우드-지향형 2 차원 정규화된 평면투영, 및 인코딩 및 디코딩 방법
KR20210152186A (ko) 인트라 예측을 이용한 포인트 클라우드 데이터 부호화/복호화 방법 및 장치
EP4244813B1 (en) Devices and methods for scalable coding for point cloud compression
WO2023142133A1 (zh) 编码方法、解码方法、编码器、解码器及存储介质
EP4233006A2 (en) Devices and methods for spatial quantization for point cloud compression
WO2024012381A1 (en) Method, apparatus, and medium for point cloud coding
WO2023123284A1 (zh) 一种解码方法、编码方法、解码器、编码器及存储介质
WO2022217472A1 (zh) 点云编解码方法、编码器、解码器及计算机可读存储介质
Daribo et al. Dynamic compression of curve-based point cloud

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20181112

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20181112

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20191114

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20191119

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20191220

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: 20200212

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200311

R150 Certificate of patent or registration of utility model

Ref document number: 6676193

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250