JP2017130188A - 三次元モデル化オブジェクトの圧縮 - Google Patents

三次元モデル化オブジェクトの圧縮 Download PDF

Info

Publication number
JP2017130188A
JP2017130188A JP2016210954A JP2016210954A JP2017130188A JP 2017130188 A JP2017130188 A JP 2017130188A JP 2016210954 A JP2016210954 A JP 2016210954A JP 2016210954 A JP2016210954 A JP 2016210954A JP 2017130188 A JP2017130188 A JP 2017130188A
Authority
JP
Japan
Prior art keywords
image
cell
grid
significance
computer
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
JP2016210954A
Other languages
English (en)
Other versions
JP7084686B2 (ja
Inventor
ファブリス ミシェル
Michel Fabrice
ファブリス ミシェル
マリカ ブルケナフェ
Boulkenafed Malika
マリカ ブルケナフェ
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.)
Dassault Systemes SE
Original Assignee
Dassault Systemes SE
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 Dassault Systemes SE filed Critical Dassault Systemes SE
Publication of JP2017130188A publication Critical patent/JP2017130188A/ja
Application granted granted Critical
Publication of JP7084686B2 publication Critical patent/JP7084686B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N19/00Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
    • H04N19/20Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using video object coding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T9/00Image coding
    • 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/20Finite element generation, e.g. wire-frame surface description, tesselation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/20Linear translation of whole images or parts thereof, e.g. panning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N13/00Stereoscopic video systems; Multi-view video systems; Details thereof
    • H04N13/20Image signal generators
    • H04N13/275Image signal generators from 3D object models, e.g. computer-generated stereoscopic image signals
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20021Dividing image into blocks, subimages or windows

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Image Processing (AREA)
  • Image Generation (AREA)

Abstract

【課題】コンピュータで実施される三次元モデル化オブジェクトの圧縮方法を提供する。
【解決手段】三次元モデル化オブジェクトのメッシュを提供するステップと、二次元平面上でメッシュをパラメータ(u,v)化するステップと、パラメータ化したメッシュを画像Iに変換するステップと、画像Iから、それぞれが制御点の集合によって定義される複数のセルからなるグリッドを定義するステップと、グリッドの各セルの関連性を算出するステップと、関連性が所定の閾値よりも低い少なくとも一つのセルを決定するステップと、少なくとも一つの決定されたセルをリサイズし、結果として得られるグリッドを算出するステップと、結果として得られるグリッドから、画像の各画素についての変換式を算出するステップと、算出された変換式を画像Iに適用するステップと、を備える。
【選択図】図1

Description

本発明はコンピュータプログラムおよびシステムの分野に関連し、特に、三次元モデル化オブジェクトの圧縮方法、システム、およびプログラムに関連する。
メッシュ単純化の分野では、メッシュの頂点の再配置および剪定を行ってメッシュ中の弧の数を減少させ、それによりメッシュに対する演算の高速化、圧縮率の向上、およびアルゴリズム全体の効率化を目指す。再メッシュ化アルゴリズムについては、Alliez他の“Recent Advances in remeshing of Surfaces, Shape Analysis and Structuring”(2008)に包括的に概説されている。
Wang他の“Content−Aware Geometry Image Resizing”、Computer Graphics International 2010(CGI 2010) Short Papers、Singapore(2010)には、ジオメトリ画像(Geometry Image)に適用するアルゴリズムであって、当該画像のサイズと内容を不均一に縮小させ、メッシュにおける、大部分のディテールを有する箇所を強調しつつ、あまり目立たない頂点を取り除く、コンテンツアウェア型画像リサイズアルゴリズムが開示されている。ジオメトリ画像の概念については、Gu他の“Geometry Images”、SIGGRAPH(2002)において、Gu他により述べられており、ここでは、ジオメトリ画像と名付けられた二次元構造体を作成するために、メッシュのパラメータ化が行われる。ジオメトリ画像において、モデルの各頂点のxyzは、xyz値を符号化するRGB値と関連付けられている。この非常に効率的な手法を実行することにより、画像のグリッド構造によって与えられる正規化メッシュが作成される。
しかしながら、Wang他の方法には、少なくとも以下の二つの欠点がある。第一に、ここで用いられるコンテンツアウェア型の画像再ターゲット化アルゴリズムは、画像を滑らかに変換するものではないため、縮小率が大きい場合には、メッシュに、明らかなアーチファクト、例えば、偽の頂点を作る。滑らかでない変形により、ジオメトリ画像から得られる導関数に極大値が作られ、さらなる圧縮を行うのが非常に困難になる。例えば、画像が縮小されるにつれ、画像の圧縮率が低くなる。これは、連続する画像およびその派生物中にアーチファクトが作られるためである。
第二に、ここで用いられるアルゴリズムは、特に画像用に作成されたものであるため、最終的な成果物が単純化されたメッシュであるという事実は完全に無視される。例えば、ジオメトリ画像の境界の規則性を保つことは、メッシュの単純化を行ううえでは無関係であるかもしれない。
こうした背景により、依然として、三次元メッシュを単純化するための方法の改良が求められている。特に、その方法は、ジオメトリ画像を滑らかかつ規則的に、選択的に縮小するものでなければならない。
本発明は、コンピュータで実施される三次元モデル化オブジェクトの圧縮方法に関する。当該方法は、三次元モデル化オブジェクトのメッシュを提供するステップと、二次元平面上で前記メッシュをパラメータ(u,v)化するステップと、前記パラメータ化したメッシュを画像Iに変換するステップと、前記画像Iから、それぞれが制御点の集合によって定義される複数のセルからなるグリッドを定義するステップと、前記グリッドの各セルの関連性を算出するステップと、関連性が所定の閾値よりも低い少なくとも一つのセルを決定するステップと、前記少なくとも一つの決定されたセルをリサイズし、結果として得られるグリッドを算出するステップと、前記結果として得られるグリッドから、前記画像の各画素についての変換式を算出するステップと、前記算出された変換式を前記画像Iに適用するステップと、を備える。
前記方法はさらに、
・関連性が所定の閾値よりも低い前記少なくとも一つのセルの前記決定は、関連性が最も低いセルの集合の決定を含み、
・前記画像I中の各画素の有意性を算出することによって有意性画像を提供するステップと、
・前記有意性画像の元の空間領域Ωを抽出するステップと、
・制御点の集合によって補間スプラインとしてパラメータ化された、元の空間領域Ωからリサイズされた空間領域Ω’への変換式Tθを提供するステップと、を備えてもよく、
・前記画像Iから、複数のセルからなるグリッドを定義することは、前記元の空間領域Ωを、それぞれが前記制御点の集合の部分集合によって定義される複数のセルに細分化することを含み、
・前記グリッドの各セルの関連性を算出することは、各セルについて、当該セル内の前記画素の前記有意性の加重平均を算出することを含み、前記加重平均は、前記補間スプラインを用いて算出され、また、前記結果として得られるグリッドから前記画像Iの各画素について変換式を算出することは、前記結果として得られるグリッドから、前記パラメータ化されたメッシュの前記変換によって、画像Iの空間領域にわたって前記変換式Tθを算出することを含み、
・前記変換式Tθは、フリーフォーム変形補間スプラインとしてパラメータ化され、また、変位関数μθ(x)は、式

によって定義され、
ここでB、B、B、およびBは三次Bスプライン関数であり、θは前記制御点を表す関数であり、sは二つの制御点間の水平方向の間隔であり、sは二つの制御点間の垂直方向の間隔であり、xは、前記元の空間領域Ω内の画素の水平軸上の位置であり、xは、前記元の空間領域Ω内の画素の垂直軸上の位置であり、また、前記制御点の集合は、グリッドを形成し、θは、制御点のグリッドを表す関数であり、
・各画素の有意性の算出は、関数

によって実行され、
ここで、

であり、

であって、geom(c)は画像Iのc番目の色チャネルであり、Gσはガウスカーネルであり、また、

であって、norm(c)は、法線の画像のc番目の色チャネルであり、Gσはガウスカーネルであり、
・各画素の有意性の算出は、さらに、ユーザ操作に応じて、前記メッシュの二以上の頂点を選択するステップと、前記ユーザによる選択を、頂点のチャートに変換するステップと、係数α=0を設定して、前記チャートの頂点にリンクされたジオメトリ画像Iの画素の有意性を最小化することを含み、
・各画素の有意性の算出は、さらに、ユーザ操作に応じて、前記メッシュの二以上の頂点を選択するステップと、前記ユーザによる選択を、頂点のチャートに変換するステップと、係数α=1を設定して、前記チャートの頂点にリンクされたジオメトリ画像Iの画素の有意性を最大化することを含み、
・前記画像Iから、複数のセルからなるグリッドを定義するステップと、前記グリッドの各セルの関連性を算出するステップと、関連性が所定の閾値よりも低い少なくとも一つのセルを決定するステップと、前記少なくとも一つの決定されたセルをリサイズし、結果として得られるグリッドを算出するステップと、前記結果として得られるグリッドから、前記画像の各画素についての変換式を算出するステップと、前記算出された変換式を前記画像Iに適用するステップとは、前記画像Iに対する前記算出された変換式が圧縮基準を満たすまで反復され、
・前記反復は、前記グリッドに、前記所定の閾値よりも関連性が低いセルがもう存在しないこと、前記画像Iに前記算出した変換式を適用することによって得られた新しい画像Jの圧縮率が、所定の閾値に達したこと、前記画像Iに前記算出した変換式を適用することによって得られた新しい画像Jの画素数が所定の閾値に達したことのうち少なくとも一つの圧縮基準が満たされたときに停止し、
・前記集合における制御点の数は、反復のたびに増加し、
・前記少なくとも一つの決定されたセルのリサイズは、前記少なくとも一つの決定されたセルを縮小することを含み、前記少なくとも一つの決定されたセルの縮小は、すべての頂点を、前記少なくとも一つの決定されたセルの頂点に近づくように、当該セルの左および右に移動させ、すべての頂点を、前記少なくとも一つの決定されたセルの頂点に近づくように、当該セルの上および下に移動させることによって行われ、
・前記算出した変換式を前記画像Iに適用することによって得られた新しい画像Jから三次元モデル化オブジェクトを構築することを備える。
さらに、上記方法を実行するための命令を備えたコンピュータプログラムが提供される。
さらに、上記コンピュータプログラムを記録した、コンピュータ読み取り可能な記憶媒体が提供される。
さらに、メモリに接続された処理回路を備えるシステムであって、前記メモリに上記コンピュータプログラムが記録されているシステムが提供される。
以下、非限定的な例として、本発明の実施の形態を、添付の図面を参照しつつ説明する。
本方法の一例についてのフローチャートを示す。 三次元モデル化オブジェクトのメッシュの一例を示す。 図2のメッシュ上におけるチャート作成を示す。 パラメータ化後の図2のメッシュを示す。 図4のパラメータ化後のメッシュから得られたジオメトリ画像を示す。 本システムの一例における画像有意性算出の一例を示す。 本システムの一例における画像有意性算出の一例を示す。 サマリー画像算出の一例を示す。 サマリー画像算出の一例を示す。 本発明により得られたジオメトリ画像を示す。 順方向マッピング変換の一例を示す。 逆方向マッピング変換の一例を示す。 Bスプライン補間グリッドの一例を示す。 Bスプライン補間グリッドの制御点変位の効果の一例を示す。 本発明を実施するコンピュータプログラムの演算モジュールの一例を示す。 本発明の方法を実行するコンピュータシステムの一例を示す。
図1のフローチャートを参照して、コンピュータで実施される三次元(3D)モデル化オブジェクトの圧縮方法を提案する。本方法は、3Dモデル化オブジェクトのメッシュを設定するステップを有する。そして、メッシュは、二次元(2D)平面上でパラメータ(u,v)化される。パラメータ化されたメッシュは、その後、画像Iに変換される。例えば、画像Iは、ジオメトリ画像である。そして、複数のセルからなるグリッドが、画像Iから定義される。グリッドの各セルは制御点の集合で定義される。このように、グリッドは、制御点の集合により形成される。そして、グリッドの各セルの関連性が算出され、すべてのセルがまとまって、一つのサマリー画像を形成する。次に、関連性が所定の閾値よりも低い一つ以上のセルを決定する。例えば最も関連性の低いk個のセルが選択される。そして、決定されたセルが、制御点を移動(変位)させることによってリサイズされ、その結果として得られるグリッドが算出される。この結果として得られるグリッドは、一つ以上の制御点が変位した後に得られる新たなグリッドである。そして、この結果として得られるグリッドから、画像の各画素の変換式が算出され、この変換式が画像Iに適用される。このような方法は、三次元メッシュの単純化を改善する。
特に、本発明により、ジオメトリ画像を、滑らかかつ規則的に縮小することが可能となる。実際、本発明では、ジオメトリ画像のうち、色が変化しない、あるいはわずかにのみ変化する領域を取り除く。ジオメトリ画像によって表されるモデルは、このように、ジオメトリ画像のうち、メッシュに関する冗長な(したがって不要な)情報を備える領域を取り除くことによって縮小される。さらに、本発明では、矩形のジオメトリ画像を維持しない。これにより、画像の境界に関係なく情報の最小化が行われる。
本方法は、コンピュータで実施される。すなわち、本方法のステップ(あるいは略全てのステップ)が少なくとも一つのコンピュータ、または類似の任意のシステムによって実行される。よって本方法のステップは、コンピュータにより、完全に自動的に、あるいは半自動的に実行される。例えば、本方法の少なくともいくつかのステップは、ユーザとコンピュータの対話を通じて始動されてもよい。求められるユーザとコンピュータの対話レベルは、想定される自動性のレベルに応じたものであって、ユーザの要望を実装する必要性との間でバランスをとるものとしてもよい。例えば、このレベルは、ユーザが定義し、かつ/または、あらかじめ定義するものであってもよい。
例えば、提供するステップは、ユーザの動作、例えば、圧縮する3Dモデルをユーザが選択することによって始動されてもよい。
コンピュータによる本方法の実施の典型的な例は、この目的に適したシステムを用いて本方法を実行することである。このシステムは、メモリに接続したプロセッサ(処理回路とも呼ぶ)を備える。システムは、さらに、グラフィカル・ユーザ・インターフェイス(GUI)を備えていてもよい。メモリ上には、本方法を実行するための命令を備えるコンピュータプログラムが記録されている。また、メモリは、データベースを記憶していてもよい。メモリは、そのような記憶に適した任意のハードウェアであり、場合により、物理的に区別可能ないくつかの部分(例えば、プログラム用に一つ、場合によりデータベース用に一つ)を備える。
本方法は、広く、モデル化オブジェクトを扱うものである。モデル化オブジェクトは、記憶したデータ、例えばデータベースに記憶したデータによって定義される任意の物体である。さらには、「モデル化オブジェクト」という表現は、データそのものも示す。システムの種類に応じて、異なる種類のデータによってモデル化オブジェクトが定義されてもよい。実際、システムは、CADシステム、CAEシステム、CAMシステム、PDMシステム、および/またはPLMシステムの任意の組み合わせであってもよい。これらの異なるシステムにおいて、モデル化オブジェクトは、対応するデータによって定義される。したがって、CADオブジェクト、PLMオブジェクト、PDMオブジェクト、CAEオブジェクト、CAMオブジェクト、CADデータ、PLMデータ、PDMデータ、CAMデータ、CAEデータについて言及することがある。しかしながら、モデル化オブジェクトは、これらのシステムの任意の組み合わせに対応するデータによって定義されてもよいため、これらのシステムは、互いに排他的なものではない。したがって、システムは、CAD兼PLMシステムである可能性もある。
CADシステムは、少なくとも、CATIAのようなモデル化オブジェクトのグラフィック表現に基づくモデル化オブジェクトの設計に適した任意のシステムをも意味する。この場合、モデル化オブジェクトを定義するデータは、モデル化オブジェクトの表現が可能なデータを備える。CADシステムは、例えば、辺や線を用い、また、ある場合には面や曲面を用いて、CADモデル化オブジェクトの表現を提供してもよい。線、辺、あるいは面は、例えば非一様有理Bスプライン(NURBS)など、様々な様式で表現されてもよい。具体的には、CADファイルは仕様を含み、その仕様に基づきジオメトリが生成可能であり、よって表現が生成可能となる。モデル化オブジェクトの仕様は一つまたは複数のCADファイルに記憶されていてもよい。CADシステムでモデル化オブジェクトを表現するファイルの典型的なサイズは、一部分あたり1メガバイトの範囲である。また、モデル化オブジェクトは、典型的には、数千の部分の集合体であってもよい。
CADの文脈において、モデル化オブジェクトは、典型的には、3Dモデル化オブジェクト、例えば、一つの部品や部品の集合体などの製品、あるいは製品の集合体を表現するものであってもよい。「3Dモデル化オブジェクト」は、3D表現が可能なデータによってモデル化される任意のオブジェクトを意味する。3D表現は、その部品をすべての角度から見ることを可能にする。たとえば、3Dで表現された3Dモデル化オブジェクトは、その軸のうちの任意の軸、あるいは、その表現が表示された画面中の任意の軸を中心に、処理して回転させることが可能である。これは、特に、3Dモデル化されていない2Dアイコンについては除外される。3D表現の表示は、設計を容易にする(すなわち、設計者が作業を達成するスピードを統計的に速める)。製品の設計は生産工程の一部であるから、これによって当該産業における生産工程が迅速化する。
3Dモデル化オブジェクトは、CADソフトウェア・ソリューションやCADシステム等を用いた仮想的デザインの完了後に実世界において生産される製品、すなわち(例えば機械的な)部品や、部品の集合(あるいは同様に部品の集合。なぜなら部品の集合は本方法の観点では一つの部品としてみることができ、また、本方法は集合体における各部品に独立して適用することができるからである)、より一般的には任意の剛体の集合体(例えば移動機構)などの製品のジオメトリを表現してもよい。CADソフトウェア・ソリューションは、航空宇宙、建築、建設、消費財、ハイテク機器、産業機器、輸送、海洋、および/または海洋石油/ガス生産、または交通を含む、限定されることのない様々な産業分野において製品の設計を可能にする。本方法により設計された3Dモデル化オブジェクトは、このように、地上車両の部品(例えば、自動車および軽トラック機器、レーシングカー、オートバイ、トラックおよびモーター機器、トラック、バス、電車を含む)、航空車両の部品(例えば、航空機体機器、航空宇宙機器、推進機器、防衛製品、航空路線機器、宇宙機器を含む)、海軍車両の部品(例えば、海軍用機器、商業用船舶、オフショア機器、ヨットおよび作業船、船舶用機器を含む)、一般的な機械部品(例えば、工業用製造機械、大型移動機械または機器、設置機器、産業機器製品、加工金属製品、タイヤ製造製品を含む)、電気機械部品または電子部品(例えば、家電、セキュリティおよび/または制御および/または計装製品、コンピューティングおよび通信機器、半導体、医療装置および設備を含む)、消費者製品(例えば、家具、家庭用および庭用製品、レジャー用品、ファッション用品、耐久消費財小売用品、織物類小売用品を含む)、包装(例えば、食品および飲料およびたばこ、美容およびパーソナルケア、家庭用製品包装を含む)などのように、任意の機械的部品であり得る工業製品を表してもよい。
図16は、本システムの一例を示す図であり、本システムは、クライアントコンピュータシステム、例えば、ユーザのワークステーションである。
本例のクライアントコンピュータは、内部通信バス1000に接続された中央演算処理装置(CPU)1010、および同じくバスに接続されたランダムアクセスメモリ(RAM)1070とを備える。クライアントコンピュータはさらに、バスに接続されたビデオランダムアクセスメモリ1100と関連付けられたグラフィックス処理装置(GPU)1110を備える。ビデオRAM1100は、当該技術分野において、フレームバッファとしても知られる。大容量記憶装置コントローラ1020は、ハードドライブ1030などの大容量記憶装置へのアクセスを管理する。コンピュータプログラムの命令及びデータを有形に具現化するのに適した大容量メモリ装置は、例として、EPROM、EEPROM及びフラッシュメモリ装置のような半導体メモリ装置、内蔵ハードディスクやリムーバブルディスクなどの磁気ディスク、光磁気ディスク、およびCD−ROMディスク1040を含む、全ての形式の不揮発性メモリを含む。前述のいずれも、特別に設計されたASIC(特定用途向け集積回路)によって補完されてもよいし、組み入れられてもよい。ネットワークアダプタ1050は、ネットワーク1060へのアクセスを管理する。クライアントコンピュータはまた、カーソル制御装置、キーボードなどの触覚装置1090を含んでいてもよい。カーソル制御装置は、ユーザがディスプレイ1080上の任意の所望の位置にカーソルを選択的に位置させることを可能にするために、クライアントコンピュータ内で使用される。さらに、カーソル制御デバイスは、ユーザが様々なコマンドを選択し、制御信号を入力することを可能にする。カーソル制御装置は、システムに制御信号を入力するための多数の信号生成装置を含む。典型的には、カーソル制御装置はマウスであってもよく、マウスのボタンは信号を生成するために使用される。あるいは、または追加的に、クライアントコンピュータシステムは、感知パッドおよび/または感知スクリーンを備えてもよい。
コンピュータプログラムは、コンピュータによって実行可能な命令を含んでいてもよく、命令は、上記システムに方法を実行させるための手段を含む。プログラムは、システムのメモリを含む任意のデータ記憶媒体に記録可能であってもよい。プログラムは、例えば、デジタル電子回路、またはコンピュータハードウェア、ファームウェア、ソフトウェア、またはそれらの組み合わせで実装されてもよい。プログラムは、例えばプログラマブルプロセッサによる実行のための機械読み取り可能な記憶装置に有形に具体化された製品のような装置として実装されてもよい。方法ステップは、プログラム可能なプロセッサが命令のプログラムを実行し、入力データを操作して出力を生成することによって方法の機能を実行することによって実行されてもよい。したがって、プロセッサは、データ記憶システム、少なくとも一つの入力デバイス、および少なくとも一つの出力デバイスからデータおよび命令を受信し、また、それらにデータおよび命令を送信するようにプログラム可能であってもよく、またそのように接続されていてもよい。アプリケーションプログラムは、高水準の手続き型またはオブジェクト指向のプログラミング言語で、または必要に応じてアセンブリ言語または機械語で実装されていてもよい。いずれの場合も、言語はコンパイラ型言語またはインタープリタ型言語であってもよい。プログラムは、フルインストールプログラムまたは更新プログラムであってもよい。いずれの場合も、プログラムをシステムに適用すると、本方法を実行するための指示が得られる。
図1に戻って、ステップS10において、3Dオブジェクトのメッシュが設定される。3Dオブジェクトは3Dモデル化オブジェクトであってもよい。このステップは、ユーザ操作、例えば、ユーザが、メモリ上に記憶され得る3Dモデル化オブジェクトを選択することによって実行される。これは、本方法を実行するコンピュータシステムによって自動的に達成されてもよい。ここで、設定するとは、3Dモデル化オブジェクトのメッシュがコンピュータシステムによって利用可能となる状態にすること、すなわち、システムがメッシュに関連するデータにアクセス可能となる状態にすることを意味する。典型的には、コンピュータシステムは、3Dモデル化オブジェクトのメッシュがシステムで利用可能になると、当該メッシュを、例えば、頭部をモデル化した3Dモデル化オブジェクトのメッシュを示す図2に示されるように表示してもよい。
次に、ステップS20において、ステップS10において設定されたメッシュが、2D平面においてパラメータ(u,v)化される。メッシュを2D平面においてパラメータ化することは、2D空間内で3Dモデル化オブジェクトの仕様を有するために必要なパラメータが決定されることを意味する。実際には、パラメータ化は、2D平面内のメッシュの各頂点の有効座標(u,v)の完全な集合を識別することを含む。パラメータ化は、当技術分野で知られているように実行され、例えば、ABF++パラメータ化アルゴリズムを使用することができる。これにより、例えば、メッシュ全体が2次元平面内で少ない歪みで展開されるようになり、または少数のチャートに切断されたメッシュが2次元平面で少ない歪みで展開されるようになる。
図4は、図2のメッシュを2D平面でパラメータ化したものを示す。メッシュのパラメータ化により、2D平面における2次元座標(u,v)を有する頂点の集合が得られる。興味深いことに、3Dモデル化オブジェクトの元のメッシュの各頂点は、2D平面内の一つの頂点によって表される。
パラメータ化の前に、メッシュをチャートにカットするステップを行うことができる。これは、ステップS10における入力とされた図2のメッシュ上でのチャート作成を示す図3に示されている。これにより、任意の属(形状におけるトポロジー的「穴」の数)のメッシュを扱うことが可能になり、当該技術分野で知られているように、単一のトポロジー的円板に縮小することが可能になる。
そして、ステップS30において、ステップS20でパラメータ化されたメッシュを、画像、すなわち画像Iに変換する。これは、メッシュのパラメータ化によって、ラスタ画像を形成する2D座標を有する頂点の集合が得られることを意味し、ここで、当該集合における各頂点はラスタ画像の画素であり、当該集合における各頂点の座標(x,y,z)は、ラスタ画像の赤/青/緑のカラーチャネルに格納される。
このように、ステップS30での出力はジオメトリ画像である。ジオメトリ画像は、Xianfeng Gu、Steven J. Gortler、Hugues Hoppe、“Proceedings of the 29th annual conference on Computer graphics and interactive techniques”、355−361ページ、ACM Press(2002)における“Geometry images”に開示されている。ジオメトリ画像は、ポリゴンメッシュの2D画像表現への変換を改善する。ジオメトリ画像は、当該技術分野で知られているように、3Dモデル化オブジェクトの任意の面を、3Dモデル化オブジェクトの元のジオメトリを完全に規則的に再メッシュ化したものである2D画像として、静的メッシュ上へ変換することによって得られ、逆変換をサポートする。ジオメトリ画像によって、ジオメトリを量子化された点の単純な2D配列として捉えた、完全に規則的な構造が提供される。法線や色のような表面信号は、同じ暗示的な面のパラメータ化(テクスチャ座標は存在しない)により、同様の2D配列に格納される。
図5は、図2のメッシュから得られたジオメトリ画像を示す。この画像を形成する各画素は頂点であり、画素の色は頂点の座標(x,y,z)を符号化する。
説明のみを目的として、以下に画像の数学的な定義について述べる。画像は、空間領域Ω(空間領域Ωは次元dの領域)から画像値の空間であるV(Vは次元cの空間)へのマッピングである。最も一般的な意味で画像を定義するために、ΩやVの性質について仮定はしない。
例えば、コンピュータメモリにおいて256階調にわたるm×n次元の2D画像Iは、マッピング(式1)として定義することができる。
この例では、空間領域Ωは次元d=2の領域であり、画像値Vの空間は次元c=1の空間である。2D画像は、離散空間領域Ω上で定義されるように、離散画像と考えることができる。これは、以下の形式の空間領域Ωで定義された連続画像とは対照的である。

離散変換は、通常、離散画像に作用する変換を指し、一方、連続変換は連続画像に作用する。離散画像空間のみがコンピュータのメモリにおいて表現可能であるため、当該技術分野で知られているように、補間のような数値法を用いて以下の整数位置における画像の値を算出する。
続くステップS40〜S80は、ジオメトリ画像、例えば、ステップS30で出力されたジオメトリ画像から、いわゆるサマリー画像を出力することを目的とする。サマリー画像は、有意性画像(Significance Image)を、スプラインに基づく画像補間に結びつける計算方法である。例えば、スプラインは、B−Splineベースのフリーフォーム変形であってもよい。しかしながら、ステップS40〜S80は、他のスプラインベースの補間、例えばNURBSを用いて実行することができる。有意性画像は、画像内の関連する部分を識別することを可能にする。有意性画像は、元の画像(例えば、ジオメトリ画像)と同じサイズであり、元の画像の各画素について、有意性画像の画素が、元の画像における対応する画素を変更するコストを定義する。ジオメトリ画像の有意性画像演算については、ステップ40を参照して述べる。
元の画像領域にわたって、フリーフォーム変形グリッド(または他の任意のスプラインベースの補間グラフ)が定義される。あるいは、複数のセルからなるグリッドが画像Iから定義され、各セルは制御点の集合によって定義される。グラフの制御点を移動すると、画像に適用される変換式が変更され、画像が変換される。例えば、この変換は微分同相変換であってもよく、この場合、結果として得られる画像が、滑らかな導関数によって滑らかになる。
セルの変形は、頂点の変位に明確にリンクされ、順に、画像の変形に影響を及ぼす。セルは、少なくとも三つの頂点で区切られたグラフの閉サイクルであり、これは、フリーフォーム変形の正則グリッドの場合には、グリッドのセルとまったく同じである。そして、変形グリッド内のどの頂点をどの程度移動すべきかを定義する問題が生じ、頂点の動きが、ある仮定の下で常に微分同相変換となることが示される。この問題に答えるために、サマリー画像が定義される。サマリー画像は、グリッドの各セルについて、関連性の演算結果を示す。したがって、サマリー画像は有意性画像の代わりとして機能し、有意性画像によって定義されるグラフ内の一つのセルを変形するコストを定量化する。サマリー画像は、変形グラフのセル数と同じ数の画素をカウントし、サマリー画像の各画素は、グラフ中の対応するセルを変形するコストを定量化する。サマリー画像の算出は、ステップ50〜80を参照して説明する。
図1のフローチャートを再び参照すると、ステップS40において、画像I、例えばジオメトリ画像における各画素の有意性(Significance)を計算することによって、有意性画像Sが得られる。画像の有意性は、当該画像の各画素に関連付けられたコストとして定義される。このコストは、画素が伝える情報量に直接関係している。画素によって伝達される情報の量は、画像内の画素の関連性であってもよいが、これに限定されるものではない。例えば、関連性はユーザによって決定されてもよい。画像Iの画素の有意性は、画像の整合性のために、当該画素の重要度を定量化するものである。したがって、画像Iの画像有意性は、画像Iにおける、他よりも大きな情報を伝達する画素を示す。
例えば図6では、キャラクタの左右の白いゾーンを縮小しても、シーンの視覚的理解は変わらない。実際、見る者は、図6の画像からも、白色の領域が縮小された画像からも、同じ情報を収集する。理想的には、画像の外観を大幅に変更することなく内容を調整することができない領域では、画像の有意性が高くなる。
実際には、画像の有意性自体が画像として表される。さらに、画像も有意性画像も、その空間領域は同じである。言い換えれば、どちらの画像における画素位置も、同じ空間位置を表す。
有意性画像の各画素に割り当てられた値は、当該画素によって伝達される情報の量を表す。これは、図6の画像から得られた有意性画像を示す図7に示すとおりである。すなわち、図7では、暗い領域は、内容の情報量が多い領域を表す。有意性画像の表現により、画像全体の情報の力に大きな影響を与えることなく、すべての白色の領域を縮小可能であることが容易に理解できる。
当該技術分野で知られている任意の画像有意性演算法を使用することができる。例えば、画像の有意性を算出するための最も基本的な方法の一つ、すなわち、画像の勾配ノルムの、平滑化されたバージョンを使用することができ、ここでの平滑化は2Dガウスフィルタを用いて行われる。図7は、画像の勾配ノルムが平滑化された図6の画像を示す。画像有意性演算の別の例として、画像顕著性演算を用いてもよい。画像の顕著性は、画像のセグメント化のために、画像における関連する部分を見つけることを目的とする。興味深いことに、画像有意性は、画像が、ビデオの場合のように、静的ではなく動的である場合には、動き情報を使用することによっても算出され得る。これは、ビデオのターゲットを変更する際、再ターゲット化されたビデオフレーム間の時間的な一貫性を維持するのに特に便利である。動き情報は、通常、画像のオプティカルフローとして符号化され、当技術分野で知られている任意の方法を用いて算出することができる。
好ましくは、有意性画像Sig Iは、式7の関数で計算される。

ここで、

である。
geomという項は、画像Iのc番目の色チャネルであり、Enormという項は、法線の画像のc番目の色チャネルである。
geomは、以下のように定義される(式7a)。

また、Enormは、式7bによって定義される。
式7aと7bの両方において、有意性をマルチスケール演算のためのスケールσで定義することを可能にする、ガウスカーネルGσを使用する。特に、ジオメトリ画像では、冗長な頂点を有する優先順位の高い領域を減らすことが重要であり、これにより、緩勾配のスムーズな領域に変換される。これが、ジオメトリ画像の有意性Sig Iの定義の要点である。パラメータ化されたメッシュ法線を、利用可能なときには利用することにより、メッシュ内の顕著な隆線を有する領域を、さらに保護することができる。
興味深いことに、各画素の有意性の算出(すなわち、有意性画像の算出)は、メッシュ内の圧縮する領域を保護するためのユーザ対話を含んでいてもよい。ユーザは、保護されるべきメッシュの頂点を選択する(例えば、ユーザは、メッシュの領域を選択する)。これは、保護される頂点のチャートに変換される。チャート全体のコストは、画像有意性において非常に高く設定される。すなわち、係数αは、チャートの頂点にリンクされたジオメトリ画像の画素に対してα=0に設定される。αの値は、0に近い値に設定してもよい。例えば、αの値は、値の範囲[0;0,2]に含まれ得る。
逆に、メッシュの領域は大幅に単純化する必要がある可能性がある。ユーザは、保護されるべきメッシュの頂点を選択する(例えば、ユーザは、メッシュの領域を選択する)。ユーザによるこの選択は、保護されるべき頂点のチャートに変換される。チャート全体のコストは、画像有意性において非常に低く設定されている。すなわち、係数αは、チャートの頂点にリンクされたジオメトリ画像の画素に対してα=1に設定される。αの値は1に近い値に設定してもよい。例えば、αの値は、値の範囲[0,8;1]に含まれ得る。
次に、ステップS50において、有意性画像の元の空間領域が抽出される。元の空間領域はΩとする。空間領域は、デジタル画像が、その画素の空間座標によって定義される領域(例えば平面)である(2D画像の場合、空間座標は、一般的には

で表される)。上述のように、重要画像の空間領域は、圧縮すべき画像Iの空間領域と同じであってもよい。有意性画像と画像Iの空間領域がともに同じ場合、これらの画像における画素位置は、ともに同じ空間位置を表す。
次に、ステップS40において、Tθで表される変換式が設定される。Tθは、元の空間領域(Ω)からリサイズされた空間領域(Ω’)までの制御点の集合によって補間スプラインとしてパラメータ化される。変換式という用語は、画像Iの画素位置に適用されてIの既存の画素位置を単純に置き換える演算子を意味し、これによりリサイズされた画像Jが生じる。よって、変換式は画素の値ではなく、画素の座標に作用する。リサイズする画像Iの画素座標の部分集合のみが関係する離散変換から、画像Iのすべての画素に関係する連続変換まで、あらゆる変換が可能である。この変換はTθで表され、Tは変換を指定することを意味し、θは変換のパラメータを指定することを意味する。これについては後述する。
説明のために、変換式Tθの数学的定義について述べる。変換式Tθのような画像変換は、入力空間領域Ωから出力空間領域Ω’へのマッピングである。入力空間領域Ωは、ステップS40で抽出された有意性画像の元の空間領域であってもよい。出力空間領域Ω’は、リサイズされた画像Jの空間領域であり得る。Ω’は、Ωのリサイズされたバージョンである。ΩとΩ’の密度は同じである。密度という用語は、空間内の値の種類(実数値、量子化値、積分値・・・)を特徴づけるものであるが、範囲が異なる可能性がある。変換はΩからΩ’へのマッピングのすべての特性を有する。すなわち、マッピングが線形である場合、変換は線形となり、マッピングが連続的であれば変換は連続的である。
リサイズされた画像Jは、変換式Tθの適用対象である、リサイズ対象の画像Iとして定義される。したがって、Ωの画素のすべての位置xに対して、次の式(式2)が適用される。
言い換えれば、画素の新しい位置y=Tθ(x)のそれぞれにおいて、I(x)の値がJ(y)にコピーされる。画素の値は変更されず、その位置のみが変更され得る。画素は、リサイズ対象の画像の座標フレームから、リサイズされた画像へ、「順方向」に移動される。すなわち、この場合、図11に示すように、順方向マッピングが実行される。ここで、斜線の入った四角は空の画素を示す。なぜなら、当該変換では、新しい画像内の非整数位置へのマッピングが行われるからである。
興味深いことに、変換式Tθは可逆的であってもよい。すなわち、当該変換は逆にすることができる。この場合、リサイズ対象の画像Iの空間領域Ωのすべての位置xについて、次の式(式3)が適用される。
式3の結果として、リサイズされた画像J内のすべての画素に、画像I内の位置が割り当てられる。画素は、リサイズされた画像Jの座標フレームから、リサイズ対象の画像Jへ、「逆方向」に移動される。すなわち、この場合、図12に示すように、逆方向マッピングが実行される。この逆方向マッピングが実行可能なのは、元の位置空間内の非整数位置を扱うときに補間手法が利用できるからである。実際には、任意の変換反転手法を用いることができる。例えば、Gary E. Christensen、Hans J. Johnson: Consistent Image Registration. IEEE Trans. Med. Imaging 20(7): 568−582(2001)を用いることができる。
ステップS60で設定される変換式Tθは、微分同相変換式Tθであってもよい。微分同相変換は、連続的で、全単射で、可微分で、全単射の導関数である変換であり、その逆もまた連続的で可微分である。有利なことに、微分同相変換は、リサイズ対象の画像について、リサイズ対象の画像の導関数についてもアーチファクトを発生させることのないように、滑らかな変換を行う。これは、区分的な滑らかな変換とは対照的である。区分的な滑らかな変換では、入力空間領域Ωが区画化され、その区画のエッジにアーチファクトが生じる。アーチファクトの各々は、画像を圧縮しにくくする高い周波数を生成する。微分同相変換のもう一つの利点は、変換が可逆的であることである。すなわち、不自然に見える変換はすべて破棄されるため、公知の折り重なり(fold−over)の問題は回避される。実際、変換が、ある箇所で可逆的でない場合、空間領域Ωにおける二つの異なる画素位置が、リサイズ後の空間領域Ω’の同一の画素位置にマッピングされることがあり得る。これにより、画像がそれ自体の上に折り重なるという視覚的効果が生じる。
変換式Tθはパラメータ化される。これは、空間領域Ωが非常に大きい可能性があり(無限大になる可能性がある)、計算が非実用的になるため、変換が空間領域Ω全体に対して明示されないことを意味する。パラメータ化された変換式は、Ωの濃度よりも小さい限られた数のパラメータを除き、Ωにおける各位置で評価可能である。
変換式Tθは、補間スプラインとしてパラメータ化される。当該技術分野で知られているように、スプラインは区分的に定義された多項式関数であり、スプラインの次数に応じて、区画間の接続における滑らかさの程度を設定することができる。区画の接点におけるスプラインの滑らかさの程度は、スプライン全体の滑らかさの程度でもある。スプラインは補間ノードとも呼ばれる制御点を有する。制御点の変位は、スプラインの形状を変更する。したがって、変換式Tθは、Tθをパラメータ化する補間スプラインの制御点の集合によってパラメータ化される。
変換関数Tθのパラメータ化に際し、任意の種類のスプラインが考えられる。実際には、画像処理に関して特に興味深い以下のようなスプラインがあるが、これらに限定されない。
・薄板スプライン。Bookstein、Fred L.“Principal warps: Thinplate splines and the decomposition of deformations”、IEEE Transactions on Pattern Analysis and Machine Intelligence archive、第11巻第6号、1989年6月、567−585ページ。
・フリーフォーム変形。Sederberg、Thomas W.、およびScott R. Parry。“Free−form deformation of solid geometric models”、Proceeding of SIGGRAPH‘86 Proceedings of the 13th annual conference on Computer graphics and interactive techniques、151−160ページ。
・非一様有理Bスプライン(NURBS)。Piegl、Les A.、およびWayne Tiller.“The NURBS book”。
この変換は、入力空間領域Ωから出力空間領域Ω’へのマッピングである。変換式Tθは、各画素について、画素位置xとその変位μθ(x)の和であってもよく、ここで、μθはΩからΩ’へのマッピングである。Ωにおけるすべての位置xについて、変換式Tθは次の式(式4)で表すことができる。

変換式Tθのこの表現を使用することにより、一般性を失うことなく、画素の変位の観点からTθを表すことが可能になり、有利である。この変換は、その位置に関連する成分と画素位置の実際の変位とに分解され、(大きい領域上で)変位が小さい場合でも、変換をかなり大きくすることができる。
ここで図13を参照すると、図13には、有意性画像から抽出された元の空間領域Ωにわたって定義された制御点の集合の例が示されている。なお、いくつかの制御点が元の空間領域Ωの外側にあることがわかる。この例では、元の空間領域Ωは次元d=2の領域であり、斜線付きの矩形で表されている。制御点の集合は、変換式Tθをパラメータ化する補間スプラインのパラメータを制御し、それらは黒点で表されている。説明のみを目的として、フリーフォーム変形(Free Form Deformations:FFD)補間スプラインに焦点を絞って議論する。
各制御点の位置はθ(i、j)で表すことができ、θ(i、j)は、(i、j)制御点の位置を表す2Dベクトルである。制御点は、インデックスiおよびjを用いてそれぞれ、左から右へ、上から下へと任意にインデックス化される。二つの制御点間の水平方向の間隔はsで表され、垂直方向の間隔はsで表すことができ、同様に、水平軸上の元の空間領域Ω内の位置xの投影はxで表され、垂直軸上の元の空間領域Ω内の位置xの投影はxで表される。リサイズ対象の画像Iの空間領域も、左から右、上から下に向いているものとみなされる。
水平方向および垂直方向は、任意に選択され、参照フレーム(x、y)を使用して定義することができる。例えば、図13において、参照フレームは、Ωの左上隅に配置された、共通の原点を有する二つの矢印で表されている。したがって、図13において、x軸は、水平方向を表し、y軸は、垂直方向を表す。
さらに図13を参照すると、制御点の集合は正則グリッドを形成し、制御点は元の空間領域Ωにわたって均一に広がっている。制御点のグリッドはFFDから取得されたものである。この例では、スプラインの次数に応じて、四つの制御点でサポートされている2Dベジェ補間スプラインが算出されている。正則グリッドという用語は、二つの制御点間の元の水平方向の間隔sの値が常に同じであることを伴う。同様に、二つの制御点間の元の垂直方向の間隔sの値も常に同じである。水平方向の間隔sは、図13に示されるように、垂直方向間隔sと同じ値を有していてもよい。なお、以下に説明されるように、図13の制御点の集合のうちの制御点は本発明のプロセス中に移動され得る。
興味深いことに、制御点の集合は、元の空間領域Ωにわたって定義される新しい領域を生成する。この新しい領域をΛとし、これについては後述する。ΛはΩに作用し、これは、制御点(例えば、図14の制御点140)の変位がスプラインに影響を及ぼし、今度はスプラインがΩに局所的に影響を及ぼすことを意味する。
前述のとおり、変換式Tθは、入力空間領域Ωから出力空間領域Ω’までのマッピングであり、ΩからΩ’までの遷移は、画素位置xとその変位μθ(x)との和である。どの画素位置からでも、変位関数μθ(x)は、制御点位置から取得され得る。したがって、式4の変位関数μθは、以下の式5に示されるように、制御点の位置の加重線型結合として定義され得る。
式5において、FFDの表現には、三次Bスプラインが使用される。すなわち関数BとBは、当該技術分野で知られているように、同じ次数のバーンスタイン多項式とも一致する三次Bスプラインである。異なる次数の任意のBスプラインを使用することができ、本発明は三次のBスプラインに限定されない。本発明は、FFDにリンクされるBスプラインに限定されず、前述のとおり、あらゆる種類のスプラインを使用することができる。
前述のとおり、各制御点の位置はθ(i、j)で表すことができる。式5において、θは、パラメータ化された変換式の制御点の集合を表す関数である。特に、関数θは、変換式Tθの唯一のパラメータである。実際、θについて知っていれば、変位μθを完全に記述することができ、ひいては、変換式Tθを完全に記述することができる。結果として、グリッド内の制御点の変位は、それに応じて変換式Tθを変化させる。制御点の変位は空間領域Ωの局所部分にのみ影響を及ぼす。実際、式5で説明されるように、制御点の位置の周囲のサイズ4s×4sの近傍のΩ内の画素位置xのみがその動きによって影響を受ける。B−Splinesの次数が低いと近傍のサイズが異なる可能性がある。
式5の三次Bスプラインに基づくFFDは、意図的に、三次まで滑らかとされている。この次数は、ベルンシュタイン多項式の次数が低いほど小さい。変位の可逆性と全単射性、およびその第一の導関数は、制御点の変位上の制約を介して達成される。すなわち、FFD変位の結果得られる変換は、制御点の変位が、制御点の間隔sとsの0.48倍よりも小さい場合、微分同相変換である。この結果は、当該技術分野において知られている。例えば、IEEE Transactions on Visualization and Computer Graphics archive、第2巻第4号、1996年12月、337〜354ページのLee、Seungyoung他、“Image metamorphosis with scattered feature constraints”。したがって、制御点の移動に起因するΩ内の位置xの変位は、制御点の変位に追従する。その結果、二つの制御点を接近させると、それら制御点の場所の間において画像が縮小され、制御点の間隔が広がると、それら制御点の場所の間において画像が拡大される。なお、変換の連続性を維持するため、Ω内における変位位置xの周囲の領域(4s×4s)もまた影響を受ける。
再び図1を参照すると、ステップS70において、元の空間領域Ωがセルに細分化される。各セルは、制御点の集合における部分集合によって定義される。したがって、セルという用語は、リサイズされた空間領域Ω’の下位部分を定義し、この下位部分は、制御点の部分集合によって区切られる。セルは、当該集合における制御点を接続するエッジで区切られる閉曲面である。セルの集合は、空間領域Ωにわたって空間領域を形成する。セルは、閉曲面であり、これは、セルが少なくとも三つの制御点によって形成されることを伴う。
図13において、セル130は、四つの制御点を接続する四つのエッジによって定義される。セル130の閉曲面内にある二つの対角線エッジは示されていない。これらの二つのエッジは、セルの閉曲面を区切ることに寄与しない。さらに図13を参照して、制御点は、セルの正則グリッドを形成する。例えば、四つの制御点によって定義されたセル130は、四つの制御点によって区切られたポリゴンである。なお、制御点の集合は、必ずしも正則グリッドを形成するわけではない。例えば、制御点の数が四より多い、あるいは四より少ないセルも使用され得る。
制御点の集合から得られたグリッドは、空間領域Ω内の制御点の座標の表現として用いることができる。
制御点は、空間領域Ω内には定義されないが、いくつかの制御点の位置は空間領域Ωと重複し、セルは、Ωと重複する制御点を用いて空間領域Ω内に定義される。
セルの集合は、元の空間領域Ωにわたるレイヤーとして見てもよい。各セルは、元の空間領域Ωにおける一つ以上の画素と関連付けられる。すなわち、各セルによって形成される領域は、一つ以上の画素をカバーする。ΩおよびΛは、分離されない。Λは、Ωを細分化したものである。ΛがΩにわたる四つのセルの正則グリッドである場合、セルインデックス(0,0)、(0,1)、(1,0)および(1,1)が得られ、各セルは、Ωのサブ領域であり、例えば、セル(0,0)は、Ωから、以下を含む。
これは、ポリゴンが画像上で描画され、ポリゴン内の一つ一つの画素がセルに属すると言うのに等しい。
セルのサイズを縮小すること、すなわちその表面積を縮小することは、セルを囲む制御点を互いに近づけることと同じであり、したがって、セルの位置で画像領域を縮小することと同じである。これは、セルのサイズを大きくする場合にも当てはまる。画像領域は、実際、セルが変形されるのと同じようにふるまう。そして、空間領域Ωの外に定義された制御点が、空間領域Ωの内部に位置する制御点の距離sとsに留まるように固定される場合、出力画像領域の空間領域Ω’は、内部セルの累積面積と同じ面積を有する。なお、固定される必要がある制御点の数は、スプラインの次数によって異なる。この変換式Tθから得られる結果は、式4に定義される。
ステップS70において、元の空間領域Ωを細分化するセルの集合が判定されている。このセルの集合は、Λで表される元の空間領域Ωの離散領域を定義する。各セルは、画像に対する画素配置と同じようにみなすことができ、同じようにインデックス化できる。例えば、四つのセルが式1の空間領域Ω上で定義されていれば、四つのセルは、それぞれ、(0,0)、(0,1)、(1,0)および(1,1)にインデックス化され、そして各セルは、空間領域Ωのサブ領域である。例えば、セル(0,0)は、Ωからの[0,1,・・・,(m/2)−1]×[0,1,・・・,(n/2)−1]を包含する。
次に、ステップS80において、セル内の画素の有意性の加重平均がΛの各セルについて算出される。その結果、Λにわたる有意性が定義される。サマリー画像Sは、Λを形成するセルの集合から形成される。サマリー画像の各セルは、サマリー画像の画素であると考えることができる。
ΛはΩの区画から成るので、領域Λは、

と表すことができる。ここで、KおよびLは、Ωの各方向における区画の数であり、Ω(i,j)は、Ωの区画であり、セル(i,j)によってカバーされる領域である。Ωの水平方向および垂直方向は、任意に選択可能であり、上述のように、参照フレーム(x、y)を用いて定義可能である。
加重平均は、補間スプラインを用いて算出される。実際には、Λ内の各位置y=(i,j)に対して、三次BスプラインFFDのサマリー画像を定義することが可能である。サマリー画像SIは、以下の式6で定義可能である。

ここで、S(x)は、設定された有意性画像を表す関数である。
ここで図8を参照すると、図7の画像の有意性画像が示されており、この画像上でセルの集合が元の空間領域Ωの離散領域Λを定義している。領域Λは、有意性画像の元の空間領域Ω上に重ねられている。図8において、領域Λは、図7の元の空間領域Ωを完全にカバーしている。さらに、セルの集合は、図13を参照して述べたように、元の空間領域Ωにわたって一様に広がる制御点のグリッドを形成する、制御点の集合から取得されている。
ここで図9を参照すると、Λを形成するセルの集合から算出されるサマリー画像が示されており、各セルについて、当該セル内の画素の有意性の加重平均が算出されている。サマリー画像の各セルは、セルの有意性を表すために用いられるグレースケールを使用して色付けされており、そのため、画素の集合の有意性がセルに含まれている。例えば、キャラクタの左側および右側の領域は白色で表され、したがって、これらの領域がシーンの視覚的理解にとって有意ではないことを意味する。
図1のフローチャートに戻り、ステップS90において、一つ以上のセルがサマリー画像から選択され、それによってセルの集合が形成される。この選択は、サマリー画像における各セルの関連性に応じて実行される。所定の閾値より低い関連性を有するセルが選択される。この例では、関連性は有意性であり、すなわち、所定の閾値よりも低い有意性を有する(一または複数の)セルが選択される。
典型的には、サマリー画像の一または複数のセルの選択は、一または複数のセルであって最も低い一つの関連性を有するセル、あるいはこの後者の場合には、最も低い複数の関連性を有する複数のセルを決定することを含み得る。例えば、有限個のk個のセルを同時に縮小するためには、k個のセルは、最も低いサマリー画像の関連性を有する。kは正の自然数である。別の例として、関連性が最も低いセルのパーセンテージpを用いることができる。なお、セルの選択は、上述の選択法に限定されない。興味深いことに、本発明のセルを縮小する際の最適化法は、矩形画像境界の制約を必要としないため、探索空間が広いほど恩恵が得られる。したがって、画像内のセルの経路に規則性を課すことを試みることなく、一度に一つずつセルを減少させることが可能である。
これで、サマリー画像のセルがその関連性に応じて選択された。サマリー画像は、ステップS100からS120で縮小される。すなわち、画像Iのメッシュに関する無用な情報を含むジオメトリ画像の着色領域が単純化される。この縮小は、サマリー画像の異なるセルについて、繰り返して行い、また、有意性画像における極小値を避け、最終的なメッシュ内の小さな領域同様、大きな領域を縮小できるように、マルチスケール法で行う。
図1のフローチャートに戻り、ステップS100において、ステップS90で前もって決定された一つ以上のセルが修正される。これは、スプラインの制御点が変位することを意味する。元の空間領域Ωのセルは、セルの関連性が各セルの画素の有意性の加重平均から計算されると、加重平均がより小さいセルの一つ以上の制御点を変位させることによって変形される。
画像をリサイズするために画素を除去または追加するコンテンツアウェア型アルゴリズムの分野で使用されることが知られているアルゴリズムのうちの一つが、サマリー画像に直接適用可能であるが、これらのアルゴリズムは、画素上ではなく、セル(すなわち、制御点のグループ)上で動作する。これらのアルゴリズムは、画像の有意性によって導かれ、画像サマリーを最小限変更しつつ、いくつかの最適性制約に応じて画像を変更するために、変換パラメータを調整する。上述のとおり、画像が真っ直ぐなエッジを保つ必要はなく、画像を少しも歪ませる必要はない。
したがって、セル上での操作は、(i)アルゴリズムが、より小さい加重平均を有する一つ以上のセルを識別し(ステップS90)、(ii)これらのアルゴリズムが行うように、画素を取り除いたり追加したりする代わりに、識別されたセルのサイズを変更する(ステップS100)ことを含む。前述のように、対応するセルのサイズが以下のように変更される。
・セルを取り除かなければならないとアルゴリズムが判定した場合、セルのサイズが縮小される。
・セルに画素を追加しなければならないとアルゴリズムが判定した場合、セルのサイズが拡大される。実際には、セルのリサイズは、変換の微分同相性を維持するために、制御点の間隔の0.48倍まで行うことができる。
ジオメトリ画像の縮小において最も典型的なケースであるセルのサイズの縮小は、すべての頂点を、縮小するセルの頂点に近づけるよう、当該セルの左および右に移動し、すべての頂点を、縮小するセルの頂点に近づけるよう、当該セルの上および下に移動することにより行う。これは、正則グリッドのセルを縮小する最も効率的な方法である。正則グリッドは、長方形または正方形のセルのみを含む。セルの縮小は、グリッドの他のセルの変位を伴うことが好ましい。例えば、縮小されたセルを囲むグリッドのセルは、一つのセルの縮小が画像の縮小を引き起こすように移動される。そうでなければ、セルが縮小されているときに縮小されていないセルは強制的に移動される。セルの縮小後に得られる空白フィールド(空き領域とも呼ばれる)は、変位した周囲のセルで満たされる。ここで、周囲のセルは、(縮小されたセルとは対照的に)縮小されていない他のセルである。なお、セルを縮小するための制御点の移動は、この例に限定されない。すなわち、単一のセルを縮小するための多くの発見的方法が採用され得る。
次に、ステップS110において、結果として得られたグリッド、すなわち、ステップS100で変形されたグリッドから、画像の各画素について変換式が算出される。
フローチャートの例では、変換式Tθは、変位した一つ以上の制御点を用いて画像Iの空間領域にわたって計算される。式5の関数θは、パラメータ化された変換式の制御点の集合を表す関数であり、関数θは変換式Tθの唯一のパラメータである。したがって、ステップS70で実行されるように、制御点が変位すると、変換式Tθのパラメータが設定される。したがって、式5で定義される変換式Tθを算出することが可能である。
最後に、ステップS120において、算出された変換式Tθをジオメトリ画像Iの空間領域に適用することによって、画像J(新しい画像Jとも呼ぶことができる)が算出される。画像Jは、リサイズされた空間領域Ω’上で定義される。算出された変換式Tθを適用することは、画像Iの入力空間領域Ωから画像Jの出力空間領域Ω’へのマッピングが実行されることを意味する。すなわち、各画素の位置が、式3および4にしたがって修正される。
上述のとおり、サマリー画像内の異なるセルについて反復的に縮小が行われ、ステップS40〜S120が、許容可能な結果が得られるまで繰り返される。反復を停止するのに、多くの基準が利用可能である。例えば、所定の閾値より低い関連性を有するグリッドのセルがそれ以上ない場合、反復を停止してもよい。別の例として、すべてのセルが変更されたときに停止してもよい。別の基準は、新しい画像Jの圧縮率が、特定の値に達するか、またはそれをさらに超えることであってもよい。反復を停止する基準のさらなる例は、新しい画像Jの画素数が所定の閾値に達することである。例えば、Jの画素数が画像Iの画素数の二分の一となることである。
ステップS40〜S120の反復処理は、所定回の反復の各セットにおいて、集合における制御点の数を増加させることをさらに含んでいてもよい。これは、マルチスケール変換グラフ更新とも呼ばれる。これにより、有意性画像における極小値を回避し、最終的なメッシュ内の小さな領域同様、大きな領域も縮小することができる。
新しい画像Jは、本方法の入力時に設定されるメッシュを表すジオメトリ画像Iの縮小版である。この画像Jは、JPEG・・・などの既知の画像圧縮アルゴリズムを用いてさらに圧縮することができるが、画像圧縮アルゴリズムがこれに限定されることはない。
新しい画像Jは、新しいメッシュを再構築するのに役立つ。構築は、画像Jを用いて行う。興味深いことに、連結性データを記憶する必要はない。実際、設定されたメッシュの変換は、図3を参照して説明したように、メッシュをチャートに分割するステップを含むことが好ましい。したがって、設定されたメッシュが穴に沿って切断されるので、任意の再メッシュ化アルゴリズムを圧縮メッシュ上で(すなわち、画像J上で)用いることができる。
図10は、図5のジオメトリ画像の新しい画像Jを示す。興味深いことに、黒で示された領域100a〜100dは新しい画像Jに属さず、単に、圧縮画像Jが真っ直ぐなエッジを保持する必要がないことを示すために存在する。なぜならジオメトリ画像を圧縮するのに必要ではないからである。
ここで図15を参照すると、本発明による3Dメッシュの圧縮を実行するためのコンピュータプログラムの演算モジュールの例が示されている。各モジュールは、コンピュータ、例えば図16のコンピュータシステムの処理回路によって実行可能な命令を含む。
ジオメトリ画像算出モジュール1510は、圧縮すべき3Dモデル化されたメッシュ1500を入力として受け取る(ステップS10)。モジュール1520は、モジュール1510によって算出されたジオメトリ画像Iの各画素の有意性を算出する(ステップS20〜30)。モジュール1520の出力は、ステップS40に従って画像Iの有意性画像を算出する。モジュール1520はさらに、ステップS50からS80に従って、サマリー画像を算出する。なお、二つの別個のモジュールがモジュール1520による算出を実行する可能性がある。すなわち、有意性画像を算出するための一つのモジュールと、サマリー画像を算出するための一つのモジュールである。このように、モジュール1520は、変換式Tθをパラメータ化する補間スプラインの制御点の集合を算出する。サマリー画像算出モジュール1520は、いくつかの種類の補間スプライン関数を実装することができる。一つのスプライン関数を、ユーザの操作に応じて、あるいはデフォルトで選択してもよい。補間スプラインから得られた制御点の集合は、グリッド、例えば、制御点が画像Iの元の空間領域Ω上に一様に広がる正則グリッドとして表現することができる。
最小セル縮小/最適化モジュールは、サマリー画像を入力として取り込む。オプティマイザモジュールは、サマリー画像のどのセルが所定の閾値よりも関連性(有意性)が低いかを返す(ステップS90)。このモジュールは、加重平均がより小さい一つ以上のセルを識別し、識別されたセルのサイズ変更を決定することが可能な、一つ以上の既知のアルゴリズムを実装する。したがって、このモジュールは、サマリー画像の変換パラメータを提供する。いくつかのアルゴリズムの中から、一つのアルゴリズムを、ユーザの操作に応じて、あるいはデフォルトで選択してもよい。
制御点移動モジュール1540は、ステップS100で定義されるように、一つ以上の制御点の変位を定義する。したがって、モジュール1540の出力は、制御点のグリッドを変更する変換パラメータである。このモジュールは、モジュール40の出力から制御点の位置を導き出す。
変換モジュール1550は、変換式Tθを算出し(ステップS110)、算出した変換式Tθを、圧縮するメッシュの画像Iの空間領域に適用する(ステップS120)。この変換モジュールの出力は、リサイズされた画像Jである。
次に、画像Iの新しいサイズ(すなわち、画像Jのサイズ)が基準に適合するかどうかがテストモジュール1570によってテストされる。画像Jが許容可能な結果に達していない場合、プロセスは繰り返される。反復は反復モジュール1580によって管理され、反復の制御はモジュール1570および1580に依存する。画像Jが期待される基準(例えばJの画素数)に適合する場合、変換モジュール1550の出力は画像Jであるとみなされる。リサイズされた画像が期待されるサイズに適合しない場合、新しい変換パラメータを用いてプロセス全体(ステップS40〜S120)が繰り返される。すなわち、モジュール1530および1540によって制御点の新しい変位が決定される。
興味深いことに、有意性画像モジュール1520で使用されるアルゴリズムで画像有意性を再計算する必要がある場合、現在の変位後のグリッドを用いて画像有意性をいつでも再計算することができる。最終的な画像は、得られた変換式が元の画像Iに適用されたときに、最後に一度だけ算出される。
本発明では、有利なことに、コンテンツアウェア型画像リサイズのためのシームカービング(Seam carving)に記載されているような、現在のコンテンツアウェア型画像再ターゲット化で既に用いられているいくつかの計算モジュールを再利用できる。既存の有意性算出モジュール1520、最適化セルモジュール1530、および変換モジュール1550が再利用可能である。これにより、開発コストを抑えることができ、また本発明を、現在のコンテンツアウェア型画像再ターゲット化ソフトウェアと容易に統合できて、有利である。
モジュールメッシュ再構成/ジッパリング(zippering)モジュール1590は、画像J、すなわち圧縮した画像Iから、メッシュを再構成する。このモジュールの出力1592は、ステップS10で設定されたメッシュと略同一のメッシュである。すなわち、復元したメッシュがステップS10の元のメッシュとわずかに異なるように、圧縮によって元のメッシュのいくつかの頂点が除去されている。図3を参照して前述したように、入力メッシュは任意の属であってもよく、まず、それを切り取って、単一のトポロジー的円板に縮小する。ジオメトリ画像には、任意の形状のチャートが敷き詰められている。ジッパリングアルゴリズムを使用して、チャートの境界間の不連続性を除去し、当該技術分野で知られているように水密面(watertight surface)を作成する。再メッシュ化したチャートは結びついている。このジッパリング処理は、通常、モジュール1590によって実行される。
以上、本発明の好適な実施形態について説明した。本発明の要旨および本発明の範囲から逸脱することなく、様々な変更がなされ得ることが理解されるであろう。したがって、他の実施態様は、以下の特許請求の範囲内にある。例えば、演算速度を増加させる多重スケール/多重解像度アプローチ(Lee、Seungyong、George Wolberg、およびSung Yong Shinに記載。「Scattered data interpolation with multilevel Bsplines」、IEEE Transactions on Visualization and Computer Graphics、第3巻第3号、1997年7月、228〜244ページ)も、本発明において、アルゴリズムに何ら変更を加えることなく利用できる。すなわち、画像の有意性を、新しいスケール/解像度レベルのそれぞれの後で再計算する必要があるだけである。

Claims (15)

  1. コンピュータで実施される、三次元モデル化オブジェクトの圧縮方法であって、
    三次元モデル化オブジェクトのメッシュを設定するステップ(S10)と、
    二次元平面上で前記メッシュをパラメータ(u,v)化するステップ(S20)と、
    前記パラメータ化したメッシュを画像Iに変換するステップ(S30)と、
    前記画像Iから、各々制御点の集合によって定義される複数のセルからなるグリッドを定義するステップ(S70)と、
    前記グリッドの各セルの関連性を算出するステップ(S80)と、
    関連性が所定の閾値よりも低い少なくとも一つのセルを決定するステップ(S90)と、
    前記少なくとも一つの決定されたセルをリサイズし、結果として得られるグリッドを算出するステップ(S100)と、
    前記結果として得られるグリッドから、前記画像Iの各画素についての変換式を算出するステップと(S110)と、
    前記算出された変換式を前記画像Iに適用するステップ(S120)と
    を有する方法。
  2. 前記関連性が所定の閾値よりも低い前記少なくとも一つのセルを決定するステップにおいて、前記関連性が最も低いセルの集合が決定される、
    請求項1に記載の、コンピュータで実施される方法。
  3. 前記画像I中の各画素の有意性を算出することによって有意性画像を提供すること(S40)と、
    前記有意性画像の元の空間領域(Ω)を抽出すること(S50)と、
    制御点の集合によって補間スプラインとしてパラメータ化された、元の空間領域(Ω)からリサイズされた空間領域(Ω’)への変換式Tθを提供すること(S60)と、
    をさらに備え、
    前記画像Iから、複数のセルからなるグリッドを定義することは、前記元の空間領域(Ω)を、それぞれが前記制御点の集合の部分集合によって定義される複数のセルに細分化することを含み、
    前記グリッドの各セルの関連性を算出することは、各セルについて、当該セル内の前記画素の前記有意性の加重平均を算出することを含み、前記加重平均は、前記補間スプラインを用いて算出され、また、
    前記結果として得られるグリッドから前記画像Iの各画素について変換式を算出することは、前記結果として得られるグリッドから、前記パラメータ化されたメッシュの前記変換によって、画像Iの空間領域にわたって前記変換式Tθを算出することを含む、
    請求項1〜2のいずれか一つに記載の、コンピュータで実施される方法。
  4. 前記変換式Tθは、フリーフォーム変形補間スプラインとしてパラメータ化され、また、変位関数μθ(x)は、式

    によって定義され、
    ここでB、B、B、およびBは三次Bスプライン関数であり、θは前記制御点を表す関数であり、sは二つの制御点間の水平方向の間隔であり、sは二つの制御点間の垂直方向の間隔であり、xは、前記元の空間領域(Ω)内の画素の水平軸上の位置であり、xは、前記元の空間領域(Ω)内の画素の垂直軸上の位置であり、また、前記制御点の集合は、グリッドを形成し、θは、制御点のグリッドを表す関数である、
    請求項3に記載の、コンピュータで実施される方法。
  5. 各画素の有意性の算出は、関数

    によって実行され、
    ここで、

    であり、

    であって、geom(c)は画像Iのc番目の色チャネルであり、Gσはガウスカーネルであり、また、

    であって、norm(c)は、法線の画像のc番目の色チャネルであり、Gσはガウスカーネルである、
    請求項3または4に記載の、コンピュータで実施される方法。
  6. 各画素の有意性の算出は、さらに、
    ユーザ操作に応じて、前記メッシュの二以上の頂点を選択するステップと、
    前記ユーザによる選択を、頂点のチャートに変換するステップと、
    係数α=0を設定して、前記チャートの頂点にリンクされた画像Iの画素の有意性を最小化するステップとを含む、
    請求項5に記載の、コンピュータで実施される方法。
  7. 各画素の有意性の算出は、さらに、
    ユーザ操作に応じて、前記メッシュの二以上の頂点を選択するステップと、
    前記ユーザによる選択を、頂点のチャートに変換するステップと、
    係数α=1を設定して、前記チャートの頂点にリンクされた画像Iの画素の有意性を最大化するステップとを含む、
    請求項5に記載の、コンピュータで実施される方法。
  8. 前記画像Iから、複数のセルからなるグリッドを定義するステップと、前記グリッドの各セルの関連性を算出するステップと、関連性が所定の閾値よりも低い少なくとも一つのセルを決定するステップと、前記少なくとも一つの決定されたセルをリサイズし、結果として得られるグリッドを算出するステップと、前記結果として得られるグリッドから、前記画像の各画素についての変換式を算出するステップと、前記算出された変換式を前記画像Iに適用するステップとは、前記画像Iに対する前記算出された変換式が圧縮基準を満たすまで反復される、
    請求項1〜7のいずれか一つに記載の、コンピュータで実施される方法。
  9. 前記反復は、
    前記グリッドに、前記所定の閾値よりも関連性が低いセルがもはや存在しないこと、
    前記画像Iに前記算出した変換式を適用することによって得られた新しい画像Jの圧縮率が、所定の閾値に達したこと、
    前記画像Iに前記算出した変換式を適用することによって得られた新しい画像Jの画素数が所定の閾値に達したこと
    のうち少なくとも一つの基準が満たされたときに停止する、
    請求項8に記載の、コンピュータで実施される方法。
  10. 前記集合における制御点の数は、反復のたびに増加する、
    請求項8または9に記載の、コンピュータで実施される方法。
  11. 前記少なくとも一つの決定されたセルのリサイズは、前記少なくとも一つの決定されたセルを縮小することを含み、前記少なくとも一つの決定されたセルの縮小は、
    すべての頂点を、前記少なくとも一つの決定されたセルの頂点に近づくように、当該セルの左右に移動させ、
    すべての頂点を、前記少なくとも一つの決定されたセルの頂点に近づくように、当該セルの上下に移動させる
    ことによって行われる、
    請求項1〜10のいずれか一つに記載の、コンピュータで実施される方法。
  12. さらに、前記算出した変換式を前記画像Iに適用することによって得られた新しい画像Jから三次元モデル化オブジェクトを構築することを備える、
    請求項1〜11のいずれか一つに記載の、コンピュータで実施される方法。
  13. 請求項1〜12のいずれかに記載の方法を実行するための命令を備えたコンピュータプログラム。
  14. 請求項13に記載のコンピュータプログラムを記録した、コンピュータ読み取り可能な記憶媒体。
  15. 請求項13に記載のコンピュータプログラムが記録されるメモリに接続された処理回路を備えるシステム。
JP2016210954A 2015-10-30 2016-10-27 三次元モデル化オブジェクトの圧縮 Active JP7084686B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP15306739.2A EP3163536B1 (en) 2015-10-30 2015-10-30 Compression of a three-dimensional modeled object
EP15306739.2 2015-10-30

Publications (2)

Publication Number Publication Date
JP2017130188A true JP2017130188A (ja) 2017-07-27
JP7084686B2 JP7084686B2 (ja) 2022-06-15

Family

ID=54608458

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016210954A Active JP7084686B2 (ja) 2015-10-30 2016-10-27 三次元モデル化オブジェクトの圧縮

Country Status (4)

Country Link
US (1) US9972129B2 (ja)
EP (1) EP3163536B1 (ja)
JP (1) JP7084686B2 (ja)
CN (1) CN106713923B (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101942370B1 (ko) * 2017-10-19 2019-01-25 인천대학교 산학협력단 자유 정점의 이동 순서 할당에 기초하여 메쉬 품질의 향상을 지원하는 그래픽 처리 장치 및 방법

Families Citing this family (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10748312B2 (en) 2016-02-12 2020-08-18 Microsoft Technology Licensing, Llc Tagging utilizations for selectively preserving chart elements during visualization optimizations
US10347017B2 (en) * 2016-02-12 2019-07-09 Microsoft Technology Licensing, Llc Interactive controls that are collapsible and expandable and sequences for chart visualization optimizations
EP3416103A1 (en) * 2017-06-16 2018-12-19 Dassault Systèmes Dataset for learning a function taking images as inputs
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
US10909725B2 (en) 2017-09-18 2021-02-02 Apple Inc. Point cloud compression
US11113845B2 (en) 2017-09-18 2021-09-07 Apple Inc. Point cloud compression using non-cubic projections and masks
US10607373B2 (en) 2017-11-22 2020-03-31 Apple Inc. Point cloud compression with closed-loop color conversion
US10699444B2 (en) * 2017-11-22 2020-06-30 Apple Inc Point cloud occupancy map compression
US10939129B2 (en) 2018-04-10 2021-03-02 Apple Inc. Point cloud compression
US11010928B2 (en) 2018-04-10 2021-05-18 Apple Inc. Adaptive distance based point cloud compression
US10909726B2 (en) 2018-04-10 2021-02-02 Apple Inc. Point cloud compression
US10909727B2 (en) 2018-04-10 2021-02-02 Apple Inc. Hierarchical point cloud compression with smoothing
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
US11012713B2 (en) 2018-07-12 2021-05-18 Apple Inc. Bit stream structure for compressed point cloud data
US10825199B2 (en) * 2018-09-12 2020-11-03 Cognex Corporation Methods and apparatus for processing image data for machine vision
US10878299B2 (en) 2018-09-12 2020-12-29 Cognex Corporation Methods and apparatus for testing multiple fields for machine vision
US10846563B2 (en) 2018-09-12 2020-11-24 Cognex Corporation Methods and apparatus for generating a dense field of three dimensional data for machine vision
US11386524B2 (en) 2018-09-28 2022-07-12 Apple Inc. Point cloud compression image padding
US11367224B2 (en) 2018-10-02 2022-06-21 Apple Inc. Occupancy map block-to-patch information compression
US11430155B2 (en) 2018-10-05 2022-08-30 Apple Inc. Quantized depths for projection point cloud compression
CN111191390B (zh) * 2018-10-26 2023-09-01 中国航发商用航空发动机有限责任公司 表面带有内凹部的零件建模方法、设备以及电子设备
CN109767480B (zh) * 2018-12-17 2023-08-18 浙江工商大学 基于显著性和空间特征的三维网格动画压缩方法
EP3671492A1 (en) * 2018-12-21 2020-06-24 Dassault Systèmes Adaptive compression of simulation data for visualization
US11348284B2 (en) 2019-01-08 2022-05-31 Apple Inc. Auxiliary information signaling and reference management for projection-based point cloud compression
WO2020154448A1 (en) * 2019-01-23 2020-07-30 Eloupes, Inc. Aligning pre-operative scan images to real-time operative images for a mediated-reality view of a surgical site
US11057564B2 (en) 2019-03-28 2021-07-06 Apple Inc. Multiple layer flexure for supporting a moving image sensor
EP4000043A4 (en) * 2019-07-19 2023-04-19 Nokia Solutions and Networks Oy MECHANISM FOR REDUCING LOG ENTRIES BASED ON CONTENT
US11562507B2 (en) 2019-09-27 2023-01-24 Apple Inc. Point cloud compression using video encoding with time consistent patches
US11627314B2 (en) 2019-09-27 2023-04-11 Apple Inc. Video-based point cloud compression with non-normative smoothing
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
US11475605B2 (en) 2020-01-09 2022-10-18 Apple Inc. Geometry encoding of duplicate points
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
CN111862046B (zh) * 2020-07-21 2023-11-17 江苏省人民医院(南京医科大学第一附属医院) 一种心脏冠脉剪影中导管位置判别系统和方法
US11948338B1 (en) 2021-03-29 2024-04-02 Apple Inc. 3D volumetric content encoding using 2D videos and simplified 3D meshes
US11544419B1 (en) 2021-07-26 2023-01-03 Pointlab, Inc. Subsampling method for converting 3D scan data of an object for marine, civil, and architectural works into smaller densities for processing without CAD processing

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013178727A (ja) * 2011-06-03 2013-09-09 Nintendo Co Ltd 立体ディスプレイにおける複数のオブジェクトを制御する装置および方法
JP2015028782A (ja) * 2013-07-30 2015-02-12 ダッソー システムズDassault Systemes 三次元モデル化オブジェクトの圧縮
JP2015122069A (ja) * 2013-12-20 2015-07-02 ダッソー システムズDassault Systemes 画像のリサイジング

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6192156B1 (en) * 1998-04-03 2001-02-20 Synapix, Inc. Feature tracking using a dense feature array
US6508768B1 (en) * 2000-11-22 2003-01-21 University Of Kansas Medical Center Ultrasonic elasticity imaging
US20050168460A1 (en) * 2002-04-04 2005-08-04 Anshuman Razdan Three-dimensional digital library system
US7265752B2 (en) * 2004-01-09 2007-09-04 Microsoft Corporation Multi-chart geometry images
US7545979B2 (en) * 2005-04-12 2009-06-09 General Electric Company Method and system for automatically segmenting organs from three dimensional computed tomography images
CN101626509B (zh) * 2009-08-10 2011-02-02 北京工业大学 三维网格编码、解码方法及编码、解码装置
US8644608B2 (en) * 2009-10-30 2014-02-04 Eiffel Medtech Inc. Bone imagery segmentation method and apparatus
CN101754039A (zh) * 2009-12-22 2010-06-23 中国科学技术大学 面向移动设备的3d参数解码系统
EP2387004B1 (en) * 2010-05-11 2016-12-14 Dassault Systèmes Lossless compression of a structured set of floating point numbers, particularly for CAD systems
KR101794537B1 (ko) * 2011-01-21 2017-11-07 삼성전자주식회사 데이터 처리 장치 및 방법
JP2012257198A (ja) * 2011-05-17 2012-12-27 Canon Inc 立体画像符号化装置、その方法、および立体画像符号化装置を有する撮像装置
CN102625126B (zh) * 2012-03-22 2014-10-15 北京工业大学 一种基于预测的三维网格编码方法
US9797225B2 (en) * 2013-11-27 2017-10-24 Saudi Arabian Oil Company Data compression of hydrocarbon reservoir simulation grids
CN103581652A (zh) * 2013-11-27 2014-02-12 重庆卓美华视光电有限公司 多视点立体视频数据处理方法及装置
US9519999B1 (en) * 2013-12-10 2016-12-13 Google Inc. Methods and systems for providing a preloader animation for image viewers
US10055857B2 (en) * 2014-08-29 2018-08-21 Ati Technologies Ulc Extension of the MPEG/SC3DMC standard to polygon meshes
US9734595B2 (en) * 2014-09-24 2017-08-15 University of Maribor Method and apparatus for near-lossless compression and decompression of 3D meshes and point clouds

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013178727A (ja) * 2011-06-03 2013-09-09 Nintendo Co Ltd 立体ディスプレイにおける複数のオブジェクトを制御する装置および方法
JP2015028782A (ja) * 2013-07-30 2015-02-12 ダッソー システムズDassault Systemes 三次元モデル化オブジェクトの圧縮
JP2015122069A (ja) * 2013-12-20 2015-07-02 ダッソー システムズDassault Systemes 画像のリサイジング

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
SHU-FAN WANG ET AL.: "Content-Aware Geometry Image Resizing", [ONLINE], JPN7020004176, June 2010 (2010-06-01), ISSN: 0004765406 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101942370B1 (ko) * 2017-10-19 2019-01-25 인천대학교 산학협력단 자유 정점의 이동 순서 할당에 기초하여 메쉬 품질의 향상을 지원하는 그래픽 처리 장치 및 방법

Also Published As

Publication number Publication date
US20170124761A1 (en) 2017-05-04
CN106713923A (zh) 2017-05-24
JP7084686B2 (ja) 2022-06-15
CN106713923B (zh) 2023-05-12
EP3163536A1 (en) 2017-05-03
EP3163536B1 (en) 2021-12-08
US9972129B2 (en) 2018-05-15

Similar Documents

Publication Publication Date Title
JP7084686B2 (ja) 三次元モデル化オブジェクトの圧縮
JP7388811B2 (ja) オートエンコーダの学習
US11144679B2 (en) Engraving a 2D image on a subdivision surface
JP7491685B2 (ja) ニューラルネットワークのセット
JP7473335B2 (ja) ソリッドcadフィーチャを推測するためのニューラルネットワークの学習
JP6463625B2 (ja) 画像のリサイジング
JP7473336B2 (ja) ソリッドcadフィーチャを推測するためのニューラルネットワークの学習
EP3379495B1 (en) Seamless fracture in an animation production pipeline
JP2018109948A (ja) パラメトリックビュー関数に基づくデータベースの照会
EP3846123B1 (en) 3d reconstruction with smooth maps
JP2022036024A (ja) パラメータ化された3dモデルを出力するニューラルネットワーク
Schollmeyer et al. Direct trimming of NURBS surfaces on the GPU
JP2022036023A (ja) 3dモデルを出力する変分オートエンコーダ
KR100392516B1 (ko) 보간되지 않은 볼륨 데이터의 실시간 렌더링 방법
JP2023160791A (ja) 画像からの3d再構成
JP2012256280A (ja) 動画像処理方法、動画像処理装置及び動画像処理プログラム

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20161028

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20191007

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20201208

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210105

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20210331

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210601

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20211207

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220330

C60 Trial request (containing other claim documents, opposition documents)

Free format text: JAPANESE INTERMEDIATE CODE: C60

Effective date: 20220330

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20220408

C21 Notice of transfer of a case for reconsideration by examiners before appeal proceedings

Free format text: JAPANESE INTERMEDIATE CODE: C21

Effective date: 20220412

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220603

R150 Certificate of patent or registration of utility model

Ref document number: 7084686

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150