(57)【特許請求の範囲】
【請求項１】
ポリゴンの集合として表現されたポリゴンモデルのデータを、立体物の造形のために削減する装置であって、
前記ポリゴンモデルを構成する各ポリゴンの各頂点に対して、同一の座標値をもつ異なるポリゴンの頂点を探索し、取得された一連の頂点に対して同一の頂点ＩＤを付与することにより、各ポリゴンに含まれる各頂点を、同一座標の頂点を同一値とする頂点ＩＤで表現し、各頂点を含むポリゴンと対応付けて記録した頂点構成テーブルと、前記頂点構成テーブルに付与された頂点ＩＤと座標値を対応付けて記録した頂点座標テーブルと、を作成するポリゴンモデル構造化手段と、
前記頂点座標テーブルに記録されている頂点と同一座標となる重複先の頂点との対応関係を記録する頂点重複テーブルを定義するテーブル定義手段と、
前記頂点構成テーブルに記録されているポリゴンに対して、面積と周長の積、面積のいずれかが所定のしきい値以下で、全ての頂点が既に削除対象に設定された全てのポリゴンの削除対象稜線の２頂点と共有しないポリゴンを削除対象ポリゴンとして順次設定し、更に各々設定されたポリゴンより削除対象稜線を選定する削除対象ポリゴン設定手段と、
前記頂点構成テーブルを参照して前記削除対象ポリゴン設定手段にて設定された全てのポリゴンの削除対象稜線の２端点を削除対象頂点とし、前記削除対象頂点を基に前記頂点座標テーブルを参照して２つの前記削除対象頂点の平均点を求め、前記削除対象頂点の１つを代表頂点とし、前記代表頂点以外の前記削除対象頂点の情報を頂点座標テーブルから削除するとともに、前記代表頂点の頂点ＩＤに、前記平均点の座標値を対応付けて記録する頂点座標テーブル更新手段と、
前記頂点重複テーブルにおいて、前記代表頂点以外の前記削除対象頂点の重複先を前記代表頂点に設定する頂点重複テーブル更新手段と、
前記頂点構成テーブルの各頂点ＩＤの値を、前記頂点重複テーブルに設定された重複先に変更し、変更された結果、同一のポリゴンに対応する頂点ＩＤが２つ以上重複する場合、当該ポリゴンを前記頂点構成テーブルより削除する頂点構成テーブル更新手段と、
前記頂点構成テーブル更新手段により更新された頂点構成テーブルの各頂点ＩＤを前記頂点座標テーブル更新手段により更新された頂点座標テーブルを参照しながら座標値に置換することにより、ポリゴンモデルを更新するポリゴンモデル更新手段と、
を有することを特徴とする立体物造形用データ削減装置。
【請求項２】
前記ポリゴンモデル更新手段により更新されたポリゴンモデルを構成するポリゴン数が所定の目標値より大きい場合、前記削除対象ポリゴン設定手段、前記頂点座標テーブル更新手段、前記頂点重複テーブル更新手段、前記頂点構成テーブル更新手段、前記ポリゴンモデル更新手段に対して繰り返して処理を実行させる制御を行う制御手段と、
を更に備えることを特徴とする請求項１に記載の立体物造形用データ削減装置。
【請求項３】
更に頂点更新テーブル更新手段を備え、
前記テーブル定義手段は、前記頂点座標テーブルに記録されている頂点が、前記削除対象頂点に設定されたことを記録する頂点更新テーブルを更に定義し、
前記削除対象ポリゴン設定手段は、前記頂点更新テーブルを参照して、前記削除対象頂点に設定されたことが記録されていないものを、前記全てのポリゴンの削除対象稜線の２頂点と共有しないポリゴンとし、
前記頂点更新テーブル更新手段は、前記頂点重複テーブルにおける頂点について更新が行われた場合に、頂点更新テーブルの当該頂点を削除対象頂点に設定することを特徴とする請求項１または請求項２に記載の立体物造形用データ削減装置。
【請求項４】
前記削除対象ポリゴン設定手段は、
前記頂点構成テーブルに記録されている全てのポリゴンに対して面積を算出し、面積の最小値に所定の１以上の値を乗算した値を前記所定のしきい値に設定するようにしていることを特徴とする請求項１から請求項３のいずれか一項に記載の立体物造形用データ削減装置。
【請求項５】
前記削除対象ポリゴン設定手段は、
前記頂点構成テーブルに記録されている全てのポリゴンに対して面積と周長を算出し、面積と周長の積の最小値に所定の１以上の値を乗算した値を前記所定のしきい値に設定するようにていることを特徴とする請求項１から請求項３のいずれか一項に記載の立体物造形用データ削減装置。
【請求項６】
前記削除対象ポリゴン設定手段は、削除対象に設定されたポリゴンより削除対象稜線を選定するにあたり、最短の稜線を選定するようにしていることを特徴とする請求項１から請求項５のいずれか一項に記載の立体物造形用データ削減装置。
【請求項７】
請求項１から請求項６のいずれか一項に記載の立体物造形用データ削減装置と、
前記立体物造形用データ削減装置から出力されたポリゴンモデルを用いて立体物を造形する立体物造形装置と、
を有することを特徴とする立体物造形システム。
【請求項８】
コンピュータを、請求項１から請求項６のいずれか一項に記載の立体物造形用データ削減装置として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、立体物を表すデータを基に、樹脂等を加工して立体物を造形する３Ｄプリンタ等の立体物造形装置の利用に際して、立体物造形装置への出力用データを削減するための技術に関する。
【背景技術】
【０００２】
近年、立体物を表すデータをもとに、樹脂や石膏等を加工して造形する３Ｄプリンタが立体物造形装置として普及してきている。３Ｄプリンタの医療応用として、医療診断で撮影されるCT/MRI画像（DICOM形式３次元ボクセルデータ）を基に３Ｄプリンタ出力用モデリングデータを作成し、手術シミュレーション・医学教育・インフォームドコンセント向けの臓器模型を出力したり、体内埋め込み用の人工臓器（血管、骨、関節など）を作成したりする試みが始まっている。
【０００３】
立体物造形用データである３Ｄプリンタ出力用モデリングデータは、ポリゴンの集合体で表現されたポリゴンモデルである。CT/MRI画像から臓器模型の３Ｄプリンタ出力用モデリングデータを作成する方法としては、特許文献１、特許文献２に示されるように、対象とする臓器領域のボクセルを抽出し、Marching Cube法（特許文献３など）により臓器領域の境界面を多面体（三角形の集合）で近似したサーフェースに変換して表面形状を表現したサーフェースデータを得る手法が一般的である。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第３６９０５０１号公報
【特許文献２】特開２０１２−２１６１０２号公報
【特許文献３】米国特許第７２０９１３６号公報
【特許文献４】特開平８−１５３２１１号公報
【特許文献５】特許第３３５０４７３号公報
【特許文献６】特許第４４６３５９７号公報
【特許文献７】特許第５０１８７２１号公報
【特許文献８】特開平１１−２３２４８９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上記のサーフェース変換手法では５１２の３乗程度のボクセルを基に、３次元的に三角形を生成する方法をとるため、三角形ポリゴンが膨大になり、ＣＧ分野においては、画面描画の際に時間がかかり対話操作の応答性が悪くなるという問題が指摘されている。これは３Ｄプリンタ出力においても問題を起こし、プリンタ出力前処理（ポリゴンから積層造形用のデータに変換する処理）に時間がかかり、場合によりワークメモリーオーバーフローを引き起こして出力不能になる。
【０００６】
これに対してＣＧ分野においては、画面の解像度に合せてポリゴン分割を粗くして曲面を表現したり、可視化に供しないポリゴンを事前に削除したりする方法がとられる。例えば、平板オブジェクトでは、単一面だけで表現し裏面を構成するポリゴンを省略したり、特許文献４のように、視点方向から隠面にあたり可視化に供しないポリゴンを事前に削除したり、特許文献５のように、他のポリゴンによって隠される内側のポリゴンを判定して事前に削除し、描画対象ポリゴンを削減する方法が提案されている。また、特許文献６のように、微細な凹凸表現に対してポリゴンを用いず、テクスチャやバンプマッピングで代用する方法が提案されている。
【０００７】
しかし、３Ｄプリンタ出力においては、画面より細かい解像度が要求され、ポリゴンにより閉空間を構成する必要があるため、平板オブジェクトにおいて、裏面を構成するポリゴンを省略することはできない。また、造形物は種々のアングルから鑑賞されるため基本的に隠面は存在せず、造形物内部にも別の形状を造形可能なため、ＣＧ分野で提案されているオクルージョンカリングという手法は意味をもたない。微細な凹凸に対しても陰影などで視覚的に擬似表現することはできず、ポリゴン形状で表現せざるを得ない。特許文献７のように、建築模型においては極端にポリゴンを削減したデフォルメされた表現が要望される場合もあるが、医療応用ではそれなりにリアリティが要求される。そのため、特許文献８のように、造形物の外観を損なわない程度に、三角形ポリゴン自体を削減する手法をとらざるを得ない。
【０００８】
三角形ポリゴンの削減にあたっては、特許文献８の手法が基本になる。これは、(1)与えられたポリゴンから削減しても全体形状にあまり影響を与えない三角形および稜線を探索する処理と、(2)対象稜線を削減するとともにそれと頂点を共有する隣接三角形を探索し補正を行う処理との２種の探索処理を所望のデータ量に削減されるまで繰り返す。特許文献８の手法では、双方の処理において削除候補とする稜線と頂点を共有する隣接三角形を探索する処理を必要とし、与えられたポリゴン数の２乗に比例して処理負荷が増大するという問題を抱えていた。
【０００９】
処理負荷を削減する方法として、削減対象の入力ポリゴンモデルのデータ量を少なくすれば良いため、入力ポリゴンモデルを分割して処理をすればよい。特許文献１では、ポリゴンに変換前のボクセル画像の段階で分割を行い、分割されたボクセルデータごとにサーフェース変換を行ってポリゴンを削減する手法を提案している。これにより確かに高速化は実現できるが、ボクセル分割された境界部にサーフェースが形成されてしまうため、ポリゴンデータを統合させると、造形物上において分割境界部に対応する不連続な亀裂が発生してしまうという問題がある（ただし、ＣＧ用途においては、レンダリングにより境界部は目立たない。）。
【００１０】
また、見かけ上は単一の閉空間にまとまって見えても、データ的には複数の閉空間に分離されている場合がある。例えば、主たる閉空間に接触した小空間が存在する場合、そのままプリンタ出力すれば支障無いが、特許文献８等によりポリゴン削減を施すと、小空間が削られて主たる閉空間からの距離が長くなり、完全に遊離した状態になることがある。この状態でプリンタ出力すると、プリンタ出力前処理でエラー停止したり、空中に浮遊するよう造形指示された小空間が出力途上で落下し出力トラブルを発生させたりする可能性がある（ただし、ＣＧ用途においては、特に問題を生じない。）。
【００１１】
そこで、本発明は、出力造形物の品質の劣化を抑えながら、ポリゴンで構成されたポリゴンモデルから、少ない処理負荷で、効率的にポリゴンを削減して立体物造形用データを得ることが可能な立体物造形用データ削減装置を提供することを課題とする。
【課題を解決するための手段】
【００１２】
上記課題を解決するため、本発明第１の態様では、
ポリゴンの集合として表現されたポリゴンモデルのデータを、立体物の造形のために削減する装置であって、
前記ポリゴンモデルを構成する各ポリゴンの各頂点に対して、同一の座標値をもつ異なるポリゴンの頂点を探索し、取得された一連の頂点に対して同一の頂点ＩＤを付与することにより、各ポリゴンに含まれる各頂点を、同一座標の頂点を同一値とする頂点ＩＤで表現し、各頂点を含むポリゴンと対応付けて記録した頂点構成テーブルと、前記頂点構成テーブルに付与された頂点ＩＤと座標値を対応付けて記録した頂点座標テーブルと、を作成するポリゴンモデル構造化手段と、
前記頂点座標テーブルに記録されている頂点と同一座標となる重複先の頂点との対応関係を記録する頂点重複テーブルを定義するテーブル定義手段と、
前記頂点構成テーブルに記録されているポリゴンに対して、面積と周長の積、面積のいずれかが所定のしきい値以下で、全ての頂点が既に削除対象に設定された全てのポリゴンの削除対象稜線の２頂点と共有しないポリゴンを削除対象ポリゴンとして順次設定し、更に各々設定されたポリゴンより削除対象稜線を選定する削除対象ポリゴン設定手段と、
前記頂点構成テーブルを参照して前記削除対象ポリゴン設定手段にて設定された全てのポリゴンの削除対象稜線の２端点を削除対象頂点とし、前記削除対象頂点を基に前記頂点座標テーブルを参照して２つの前記削除対象頂点の平均点を求め、前記削除対象頂点の１つを代表頂点とし、前記代表頂点以外の前記削除対象頂点の情報を頂点座標テーブルから削除するとともに、前記代表頂点の頂点ＩＤに、前記平均点の座標値を対応付けて記録する頂点座標テーブル更新手段と、
前記頂点重複テーブルにおいて、前記代表頂点以外の前記削除対象頂点の重複先を前記代表頂点に設定する頂点重複テーブル更新手段と、
前記頂点構成テーブルの各頂点ＩＤの値を、前記頂点重複テーブルに設定された重複先に変更し、変更された結果、同一のポリゴンに対応する頂点ＩＤが２つ以上重複する場合、当該ポリゴンを前記頂点構成テーブルより削除する頂点構成テーブル更新手段と、
前記頂点構成テーブル更新手段により更新された頂点構成テーブルの各頂点ＩＤを前記頂点座標テーブル更新手段により更新された頂点座標テーブルを参照しながら座標値に置換することにより、ポリゴンモデルを更新するポリゴンモデル更新手段と、
を有することを特徴とする立体物造形用データ削減装置を提供する。
【００１３】
本発明第１の態様によれば、ポリゴンモデルを構成する各ポリゴンの各頂点に対して、同一の座標値をもつ異なるポリゴンの頂点を探索し、取得された一連の頂点に対して同一の頂点ＩＤを付与することにより、各ポリゴンに含まれる各頂点を、同一座標の頂点を同一値とする頂点ＩＤで表現し、各頂点を含むポリゴンと対応付けて記録した頂点構成テーブルと、頂点構成テーブルに付与された頂点ＩＤと座標値を対応付けて記録した頂点座標テーブルとを作成し、頂点座標テーブルに記録されている頂点と同一座標となる重複先の頂点との対応関係を記録する頂点重複テーブルを定義し、頂点構成テーブルに記録されているポリゴンに対して、面積と周長の積、面積のいずれかが所定のしきい値以下で、全ての頂点が既に削除対象に設定された全てのポリゴンの削除対象稜線の２頂点と共有しないポリゴンを削除対象ポリゴンとして順次設定し、更に各々設定されたポリゴンより削除対象稜線を選定し、頂点構成テーブルを参照して削除対象ポリゴン設定手段にて設定された全てのポリゴンの削除対象稜線の２端点を削除対象頂点とし、削除対象頂点を基に頂点座標テーブルを参照して複数の削除対象頂点の平均点を求め、削除対象頂点の１つを代表頂点とし、代表頂点以外の削除対象頂点の情報を頂点座標テーブルから削除するとともに、代表頂点の頂点ＩＤに、平均点の座標値を対応付けて記録するようにし、頂点重複テーブルにおいて、代表頂点以外の削除対象頂点の重複先を代表頂点に設定するようにし、頂点構成テーブルの各頂点ＩＤの値を、更新された頂点重複テーブルに設定された重複先に変更し、変更された結果、同一のポリゴンに対応する頂点ＩＤが２つ以上重複する場合、そのポリゴンを頂点構成テーブルより削除するようにし、更新された頂点構成テーブルの各頂点ＩＤを、更新された頂点座標テーブルを参照しながら座標値に置換することにより、ポリゴンモデルを更新するようにしたので、設定された削除対象ポリゴンに隣接するポリゴンの探索負荷を減少させることができる。このため、出力造形物の品質の劣化を抑えながら、ポリゴンで構成されたポリゴンモデルから、少ない処理負荷で、効率的にポリゴンを削減して立体物造形用データを得ることが可能となる。特に、医療用のCT/MRI画像をもとに生成されたサーフェースデータ等のポリゴンモデルを効率良く３Ｄプリンタ等の立体物造形装置に出力することが可能となる。
【００１４】
また、本発明第２の態様では、前記ポリゴンモデル更新手段により更新されたポリゴンモデルを構成するポリゴン数が所定の目標値より大きい場合、前記削除対象ポリゴン設定手段、前記頂点座標テーブル更新手段、前記頂点重複テーブル更新手段、前記頂点構成テーブル更新手段、前記ポリゴンモデル更新手段に対して繰り返して処理を実行させる制御を行う制御手段と、を更に備えることを特徴とする。
【００１５】
本発明第２の態様によれば、更新されたポリゴンモデルを構成するポリゴン数が目標値より大きい場合、削除対象ポリゴンの設定、頂点座標テーブルの更新、頂点重複テーブルの更新、頂点構成テーブルの更新、ポリゴンモデルの更新を繰り返して行うようにしたので、ポリゴンモデルを構成するポリゴン数を目標値まで削減することが可能となる。
【００１６】
また、本発明第３の態様では、
更に頂点更新テーブル更新手段を備え、
前記テーブル定義手段は、前記頂点座標テーブルに記録されている頂点が、前記削除対象頂点に設定されたことを記録する頂点更新テーブルを更に定義し、
前記削除対象ポリゴン設定手段は、前記頂点更新テーブルを参照して、前記削除対象頂点に設定されたことが記録されていないものを、前記全てのポリゴンの削除対象稜線の２頂点と共有しないポリゴンとし、
前記頂点更新テーブル更新手段は、前記頂点重複テーブルにおける頂点について更新が行われた場合に、頂点更新テーブルの当該頂点を削除対象頂点に設定することを特徴とする。
【００１７】
本発明第３の態様によれば、頂点座標テーブルに記録されている頂点が、削除対象頂点に設定されたことを記録する頂点更新テーブルを更に定義し、削除対象ポリゴン設定の際、頂点更新テーブルを参照して、削除対象頂点に設定されたことが記録されていないものを、全てのポリゴンの削除対象稜線の２頂点と共有しないポリゴンとし、頂点重複テーブルにおいて所定の頂点について更新が行われた場合に、頂点更新テーブルの対応する頂点を削除対象頂点に設定するようにしたので、削除対象に設定しようとするポリゴンのいずれの頂点と共有する頂点をもつ他の隣接ポリゴンが削除対象に事前に設定されていないことを容易に確認でき、隣接する複数のポリゴンを同時に削除対象に設定することを回避でき、削除後のポリゴンモデルを歪ませることを防ぐことができる。
【００１８】
また、本発明第４の態様では、前記削除対象ポリゴン設定手段は、前記頂点構成テーブルに記録されている全てのポリゴンに対して面積を算出し、面積の最小値に所定の１以上の値を乗算した値を前記所定のしきい値に設定するようにしていることを特徴とする。
【００１９】
本発明第４の態様によれば、頂点構成テーブルに記録されている全てのポリゴンに対して面積を算出し、面積の最小値に所定の１以上の値を乗算した値を前記所定のしきい値に設定するようにしたので、ポリゴン面積が小さいものを優先してポリゴンの削除を行うことが可能となる。
【００２０】
また、本発明第５の態様では、前記削除対象ポリゴン設定手段は、前記頂点構成テーブルに記録されている全てのポリゴンに対して面積と周長を算出し、面積と周長の積の最小値に所定の１以上の値を乗算した値を前記所定のしきい値に設定するようにていることを特徴とする。
【００２１】
本発明第５の態様によれば、頂点構成テーブルに記録されている全てのポリゴンに対して面積と周長を算出し、面積と周長の積の最小値に所定の１以上の値を乗算した値を前記所定のしきい値に設定するようにしたので、面積と周長の積が小さいものを優先してポリゴンの削除を行うことが可能となる。
【００２２】
また、本発明第６の態様では、前記削除対象ポリゴン設定手段は、削除対象に設定されたポリゴンより削除対象稜線を選定するにあたり、最短の稜線を選定するようにしていることを特徴とする。
【００２３】
本発明第６の態様によれば、削除対象ポリゴン設定手段が、削除対象に設定されたポリゴンより削除対象稜線を選定するにあたり、最短の稜線を選定するようにしたので、ポリゴンモデルの全体形状に歪をあまり与えずにポリゴンの削除を行うことが可能となる。
【００２４】
また、本発明第７の態様では、
前記立体物造形用データ削減装置と、
前記立体物造形用データ削減装置から出力されたポリゴンモデルを用いて立体物を造形する立体物造形装置と、
を有することを特徴とする立体物造形システム
を有することを特徴とする立体物造形システムを提供する。
【００２５】
本発明第７の態様によれば、立体物造形用データ削減装置と、立体物造形用データ削減装置から出力されたポリゴンモデルを用いて立体物を造形する立体物造形装置により立体物造形システムを実現するようにしたので、ボードコンピュータを組み込んだ３Ｄプリンタ等の形態で、立体物造形システムを提供することが可能となる。
【発明の効果】
【００２６】
本発明によれば、出力造形物の品質の劣化を抑えながら、ポリゴンで構成されたポリゴンモデルから、少ない処理負荷で、効率的にポリゴンを削減して立体物造形用データを得ることが可能となる。
【図面の簡単な説明】
【００２７】
【図１】本発明の一実施形態に係る立体物造形用データ削減装置のハードウェア構成図である。
【図２】本発明の一実施形態に係る立体物造形用データ削減装置の構成を示す機能ブロック図である。
【図３】DICOM形式ボクセルデータからポリゴンモデルへの変換の様子を示す図である。
【図４】DICOM形式ボクセルデータを、複数のスライス画像として示した図である。
【図５】DICOM形式ボクセルデータを、４つの形態で表示した例である。
【図６】図５に示したDICOM形式ボクセルデータを「3D-Slicer」により変換して得られたＳＴＬデータを示す図である。
【図７】領域別ポリゴンモデル（ＳＴＬデータ）を示す図である。
【図８】領域別ＳＴＬデータを合成したポリゴンモデル（ＳＴＬデータ）のレンダリング像を示す図である。
【図９】本発明の一実施形態に係る立体物造形用データ削減装置の処理概要を示すフローチャートである。
【図１０】処理対象とするポリゴン群の一例を示す図である。
【図１１】図１０に示したポリゴン群を実現するための頂点座標配列データを示す図である。
【図１２】ステップＳ２０におけるポリゴン削減処理の詳細を示すフローチャートである。
【図１３】ステップＳ２１のポリゴンモデルの構造化の詳細を示すフローチャートである。
【図１４】ポリゴン群とハッシュ値の関係を示す図である。
【図１５】ポリゴンモデル構造化の初期化時における各テーブルの状態を示す図である。
【図１６】ポリゴンモデル構造化における各テーブルの変化の様子を示す図である。
【図１７】ポリゴンモデル構造化における各テーブルの変化の様子を示す図である。
【図１８】ポリゴンモデル構造化における各テーブルの変化の様子を示す図である。
【図１９】ポリゴンモデル構造化における各テーブルの変化の様子を示す図である。
【図２０】ポリゴンモデル構造化における各テーブルの変化の様子を示す図である。
【図２１】初期化時における頂点構成テーブル、頂点座標テーブル、頂点重複テーブル、頂点更新テーブルを示す図である。
【図２２】Ｐ６とＰ１０の２つの削除対象ポリゴンの一辺が削除対象辺として設定された状態のポリゴン群を示す図である。
【図２３】稜線ＪＫの処理時における各テーブルの様子を示す図である。
【図２４】稜線ＤＥの処理時における各テーブルの様子を示す図である。
【図２５】頂点重複テーブルによる頂点構成テーブルの更新の様子を示す図である。
【図２６】頂点構成テーブルからのポリゴンの削除の様子を示す図である。
【図２７】ポリゴン削除処理後の頂点座標配列データを示す図である。
【図２８】ポリゴン削除処理後のポリゴン群の状態を示す図である。
【図２９】ポリゴン削減処理により生じる問題を説明するための図である。
【図３０】隔離された小体積グループの削除処理の詳細を示すフローチャートである。
【図３１】図３０のステップＳ４４における最小間隔の算出処理の詳細を示すフローチャートである。
【発明を実施するための形態】
【００２８】
以下、本発明の好適な実施形態について図面を参照して詳細に説明する。
＜１．装置構成＞
図１は、本発明の一実施形態に係る立体物造形用データ削減装置２００を含む立体物造形システムのハードウェア構成図である。本実施形態に係る立体物造形用データ削減装置２００は、汎用のコンピュータで実現することができ、図１に示すように、ＣＰＵ（Central Processing Unit）１と、コンピュータのメインメモリであるＲＡＭ（Random Access Memory）２と、ＣＰＵ１が実行するプログラムやデータを記憶するためのハードディスク、フラッシュメモリ等の大容量の記憶装置３と、キーボード、マウス等のキー入力Ｉ／Ｆ（インターフェース）４と、３Ｄプリンタやデータ記憶媒体等の外部装置とデータ通信するためのデータ入出力Ｉ／Ｆ（インターフェース）５と、液晶ディスプレイ等の表示デバイスである表示部６と、を備え、互いにバスを介して接続されている。
【００２９】
３Ｄプリンタ７は、汎用の３Ｄプリンタであり、立体物の三次元形状をポリゴンの集合で表現したポリゴンモデルである立体物造形用データを基に樹脂、石膏等の素材を加工して立体物を造形する立体物造形装置である。３Ｄプリンタ７は、データ処理部７ａと出力部７ｂを有している。３Ｄプリンタ７のデータ処理部７ａは、データ入出力Ｉ／Ｆ５に接続されており、データ入出力Ｉ／Ｆ５から受け取った削減された立体物造形用データを基に、出力部７ｂが立体物を造形するようになっている。
【００３０】
図１では、立体物造形用データ削減装置２００と３Ｄプリンタ７は分離した形態で示されているが、現在市販されている殆どの３Ｄプリンタ製品には立体物造形用データ削減装置２００の構成要素である、ＣＰＵ１、ＲＡＭ２、記憶装置３、キー入力Ｉ／Ｆ４（汎用コンピュータ向けキーボード・マウスではなく、テンキーレベルの数種のボタン）、データ入出力Ｉ／Ｆ５、表示部６（数行の文字を表示可能な小型液晶パネル、タッチパネルを重畳させキー入力Ｉ／Ｆ４を兼ねることも多い）も小規模ながら重複して備えている。従って、３Ｄプリンタ７自体が外部記憶媒体経由で立体物造形用データを直接受け取り、単独で立体物を造形する運用も可能になっている（特に民生用の３Ｄプリンタではこちらの形態の方が多い）。すなわち、図１に示した立体物造形システムを１つの筐体に収めて、“３Ｄプリンタ”という製品として流通することも多い。
【００３１】
図２は、本実施形態に係る立体物造形用データ削減装置の構成を示す機能ブロック図である。図２において、１０は制御手段、２０はグループ分類手段、３０はポリゴンモデル構造化手段、４０はテーブル定義手段、５０は削除対象ポリゴン設定手段、６０は頂点座標テーブル更新手段、７０は頂点重複テーブル更新手段、７５は頂点更新テーブル更新手段、８０は頂点構成テーブル更新手段、９０はポリゴンモデル更新手段、１００はグループ合成手段、１１０はポリゴンモデル合成手段、１２０は領域別ポリゴンモデル記憶手段、１３０は合成ポリゴンモデル記憶手段である。
【００３２】
制御手段１０は、立体物造形用データ削減装置全体の制御を行う。グループ分類手段２０は、削除対象ポリゴンの選定対象のポリゴンモデルに対して、各グループに属するポリゴンがそのグループ内の他のいずれかのポリゴンと、ポリゴンを構成する辺（２頂点）を共有するようにグループ分類する。ポリゴンモデル構造化手段３０は、ポリゴンモデルから、各ポリゴンに含まれる各頂点を、同一座標の頂点を同一値とする頂点ＩＤで表現し、ポリゴンを構成する各頂点を頂点ＩＤで表現してポリゴンと対応付けて記録した頂点構成テーブルと、頂点構成テーブルに付与された頂点ＩＤと座標値を対応付けて記録した頂点座標テーブルを作成する。テーブル定義手段４０は、頂点座標テーブルに記録されている頂点ＩＤと同一座標となる重複先の頂点ＩＤとの対応関係を記録する頂点重複テーブルと、頂点ＩＤと更新フラグを対応付けた頂点更新テーブルを定義する。削除対象ポリゴン設定手段５０は、頂点構成テーブルに記録されているポリゴンに対して所定の間隔で一部のポリゴンを削除対象ポリゴンとして設定する。
【００３３】
頂点座標テーブル更新手段６０は、頂点構成テーブルを参照して削除対象ポリゴンを構成する複数の頂点ＩＤを削除対象頂点として特定し、削除対象頂点を基に頂点座標テーブルを参照して複数の削除対象頂点の平均点を求め、削除対象頂点の１つを代表頂点とし、代表頂点以外の削除対象頂点の情報を頂点座標テーブルから削除するとともに、代表頂点の頂点ＩＤに、平均点の座標値を対応付けて記録する。頂点重複テーブル更新手段７０は、頂点重複テーブルにおいて、代表頂点以外の全ての削除対象頂点の重複先を代表頂点に設定する。頂点更新テーブル更新手段７５は、頂点更新テーブルにおいて、更新された頂点の更新フラグを設定する。頂点構成テーブル更新手段８０は、頂点構成テーブルの中で、更新された頂点重複テーブルにおいて更新された頂点ＩＤの値を、更新された重複先に変更し、変更された結果、ポリゴンを構成する同一の値をもつ頂点ＩＤが２つ以上重複する場合、当該ポリゴンを頂点構成テーブルより削除する。ポリゴンモデル更新手段９０は、更新された頂点構成テーブルの各頂点ＩＤを、更新された頂点座標テーブルを参照しながら座標値に置換することにより、削減されたポリゴンの配列で表現されたポリゴンモデルを作成する。
【００３４】
グループ合成手段１００は、分類された各グループに対するポリゴンモデル構造化手段３０、テーブル定義手段４０、削除対象ポリゴン設定手段５０、頂点座標テーブル更新手段６０、頂点重複テーブル更新手段７０、頂点更新テーブル更新手段７５、頂点構成テーブル更新手段８０、ポリゴンモデル更新手段９０による一連の処理の実行により作成された複数のグループの更新されたポリゴンモデルを、領域別に単一の領域別ポリゴンモデルに合成する。ポリゴンモデル合成手段１１０は、複数セットの更新された領域別ポリゴンモデルを単一の合成ポリゴンモデルに合成する。制御手段１０、グループ分類手段２０、ポリゴンモデル構造化手段３０、テーブル定義手段４０、削除対象ポリゴン設定手段５０、頂点座標テーブル更新手段６０、頂点重複テーブル更新手段７０、頂点更新テーブル更新手段７５、頂点構成テーブル更新手段８０、ポリゴンモデル更新手段９０、グループ合成手段１００、ポリゴンモデル合成手段１１０は、ＣＰＵ１が、記憶装置３に記憶されているプログラムを実行することにより実現される。
【００３５】
領域別ポリゴンモデル記憶手段１２０は、領域別ポリゴンモデルを記憶した記憶手段であり、記憶装置３により実現される。領域別ポリゴンモデルとは、１つの領域を１セットとして構成されたポリゴンモデルである。合成ポリゴンモデル記憶手段１３０は、複数セットの領域別ポリゴンモデルを合成した単一のポリゴンモデルである合成ポリゴンモデルを記憶する記憶手段であり、記憶装置３により実現される。合成ポリゴンモデル記憶手段１３０に記憶される合成ポリゴンモデルは、３Ｄプリンタ７への出力のために削減された立体物造形用データとなる。領域別ポリゴンモデル記憶手段１２０に記憶された領域別ポリゴンモデル、合成ポリゴンモデル記憶手段１３０に記憶される合成ポリゴンモデルは、いずれもポリゴンの集合体であり、ポリゴン単位で各ポリゴンが有する頂点を記録した形式となっている。
【００３６】
図２に示した各構成手段は、現実には図１に示したように、コンピュータおよびその周辺機器等のハードウェアに専用のプログラムを搭載することにより実現される。すなわち、コンピュータが、専用のプログラムに従って各手段の内容を実行することになる。なお、本明細書において、コンピュータとは、ＣＰＵ等の演算処理部を有し、データ処理が可能な装置を意味し、パーソナルコンピュータなどの汎用コンピュータだけでなく、製品としての“３Ｄプリンタ”に組み込まれたボードコンピュータも含む。
【００３７】
図１に示した記憶装置３には、ＣＰＵ１を動作させ、コンピュータを、立体物造形用データ削減装置として機能させるための専用のプログラムが実装されている。この専用のプログラムを実行することにより、ＣＰＵ１は、制御手段１０、グループ分類手段２０、ポリゴンモデル構造化手段３０、テーブル定義手段４０、削除対象ポリゴン設定手段５０、頂点座標テーブル更新手段６０、頂点重複テーブル更新手段７０、頂点更新テーブル更新手段７５、頂点構成テーブル更新手段８０、ポリゴンモデル更新手段９０、グループ合成手段１００、ポリゴンモデル合成手段１１０としての機能を実現することになる。また、記憶装置３は、領域別ポリゴンモデル記憶手段１２０、合成ポリゴンモデル記憶手段１３０として機能するだけでなく、立体物造形用データ削減装置としての処理に必要な様々なデータを記憶する。
【００３８】
＜２．処理動作＞
＜２．１．前処理＞
次に、図１、図２に示した立体物造形用データ削減装置の処理動作について説明する。まず、医療用ボクセルデータからサーフェースデータである領域別ポリゴンモデルを作成する前処理について説明する。サーフェースデータとは、サーフェースモデリングにより立体構造を表面形状として表現したデータであり、多角形のデータであるポリゴンにより構成されるポリゴンモデルである。前処理は、コンピュータが公知のプログラムを実行することにより行われる。ここでは、公知のプログラムとして、米国ハーバード大学が中心に開発したオープンソースソフトウェアである「3D-Slicer」を用いた場合を例にとって説明する。まず、コンピュータが、全体ボクセルデータから、特定臓器領域のボクセルデータを抽出する。次に、コンピュータが、特定臓器領域のボクセルデータをポリゴンモデルであるサーフェースデータに変換する。このようにして得られた領域別ポリゴンモデルは、領域別ポリゴンモデル記憶手段１２０に格納される。前処理は、図１に示した立体物造形用データ削減装置を実現するコンピュータに、上記プログラムを実行させることにより行ってもよいし、図１に示した立体物造形用データ削減装置を実現するコンピュータとは別のコンピュータにより実行し、得られた領域別ポリゴンモデルを記憶装置３により実現される領域別ポリゴンモデル記憶手段１２０に格納するようにしてもよい。
【００３９】
本実施形態では、医療用ボクセルデータとして医用画像データ交換に関する国際標準規格であるDICOM形式のボクセルデータを用いる。また、ポリゴンモデルであるサーフェースデータとしてＳＴＬ（Standard Triangulated Language）データを用いる。ＳＴＬデータとは、元来３Ｄグラフィックス分野におけるモデリングデータ交換に関する業界標準で、最近では３Ｄプリンタでも標準的に使用されるようになったもので、三次元形状を三角形のポリゴンの集合体で表現したデータである。そして、コンピュータが上記「3D-Slicer」を実行することにより、図３（ａ）に示すようなDICOM形式ボクセルデータから、図３（ｂ）に示すような領域別のボクセルデータを抽出し、サーフェースモデル変換を行って図３（ｃ）に示すようなポリゴンモデル（ＳＴＬデータ）が得られる。領域別のボクセルデータの抽出、サーフェースモデル変換には、特許文献３等に記載されたMarching Cube法が用いられている。
【００４０】
図４〜図８に、「3D-Slicer」による処理画面の一例を示す。図４は、「3D-Slicer」の処理対象となるDICOM形式ボクセルデータを、複数のスライス画像として示した図である。図４の例では、人体頭部のサジタル像を２５６枚のスライス画像で示している。各スライス画像は、５１２×５１２画素で構成されている。図５は、DICOM形式ボクセルデータを、４つの形態で表示した例である。図５において、上部に表示されているのは、ボリュームレンダリング像である。下部に表示されているのは、左側から順にアキシャル像、サジタル像、コロナル像である。
【００４１】
図６は、図５に示したDICOM形式ボクセルデータを「3D-Slicer」により変換して得られたポリゴンモデルのＳＴＬデータを示す図である。図６において、上部に表示されているのは、ＳＴＬデータのサーフェースレンダリング像である。図６の例では、眼球・脳領域しきい値：626.75以上（MRI撮影でボクセル値10bits幅の場合）、眼球結合組織Ｌしきい値：298.36-458.74、眼球結合組織Ｒしきい値：298.36-431.32を設定することにより、領域（組織）別に分類された領域別ＳＴＬデータをサーフェースレンダリングしたものを示している。
【００４２】
図７は、「3D-Slicer」により得られた領域別ＳＴＬデータをシェーディングにより陰影表現した図である。図８は、領域別ＳＴＬデータを合成したポリゴンモデル（ＳＴＬデータ）のサーフェースレンダリング像を示す図である。図８に示すFront Viewは、図６の上部に示したものと同一である。前処理の後は、図７に示したような領域別ポリゴンモデルが領域別ポリゴンモデル記憶手段１２０に記憶される。
【００４３】
＜２．２．処理概要＞
次に、図１、図２に示した立体物造形用データ削減装置の処理動作について説明する。図９は、本発明の一実施形態に係る立体物造形用データ削減装置の処理概要を示すフローチャートである。まず、グループ分類手段２０が、領域別ポリゴンモデル記憶手段１２０から、１つの領域についての領域別ポリゴンモデルを読み込み、複数のグループに分類する（ステップＳ１０）。次に、各グループのポリゴンモデルに対して、ポリゴンモデル構造化手段３０、テーブル定義手段４０、削除対象ポリゴン設定手段５０、頂点座標テーブル更新手段６０、頂点重複テーブル更新手段７０、頂点更新テーブル更新手段７５、頂点構成テーブル更新手段８０、ポリゴンモデル更新手段９０が、ポリゴン削減処理を実行する（ステップＳ２０）。制御手段１０が、全グループについてポリゴン削減処理が終了したか否かを判定し（ステップＳ３０）、全グループについて終了していない場合は、未処理のグループに対してステップＳ２０のポリゴン削減処理を実行する。全グループに対するポリゴン削減処理が終了した場合は、グループ合成手段１００が、隔離された小体積グループの削除を行う（ステップＳ４０）。続いて、グループ合成手段１００は、全グループのポリゴンモデルの合成を行う（ステップＳ５０）。制御手段１０が、全領域について処理が終了したか否かを判定し（ステップＳ６０）、全領域について終了していない場合は、未処理の領域に対してステップＳ１０〜Ｓ５０の処理を実行する。全領域に対する処理が終了した場合は、ポリゴンモデル合成手段１１０が、全領域の領域別ポリゴンモデルの合成を行う（ステップＳ７０）。
【００４４】
＜２．３．グループ分類処理＞
まず、ステップＳ１０において、グループ分類手段２０がグループ分類処理を行う。具体的には、領域別ポリゴンモデルに対して、各グループに属するポリゴンが当該グループ内の他のいずれかのポリゴンと、ポリゴンを構成する辺を共有するようにグループ分類する。グループ分類処理の具体的な手法としては、公知の様々な手法を用いることができる。
【００４５】
＜２．４．ポリゴン削減処理＞
次に、ステップＳ２０におけるポリゴンモデル構造化手段３０、テーブル定義手段４０、削除対象ポリゴン設定手段５０、頂点座標テーブル更新手段６０、頂点重複テーブル更新手段７０、頂点更新テーブル更新手段７５、頂点構成テーブル更新手段８０、ポリゴンモデル更新手段９０によるポリゴン削減処理について説明する。ポリゴン削減処理におけるポリゴンの削除の手法として、ポリゴンを構成する稜線を特定し、その稜線を削除することによりポリゴンを削除する手法を用いる。ここで、稜線とは、１つのポリゴンの２頂点を結ぶ辺を意味する。
【００４６】
図１０は、処理対象とするポリゴン群の一例を示す図である。ポリゴンモデルは、図１１に示されるように空間における三次元の値をもつポリゴン群であるが、図１０の例では、紙面上での説明の便宜上、二次元空間に投影したポリゴン群を示している。図１０の例では、削減前、１８個の三角形のポリゴンからなるポリゴン群を示している。図１０において、大文字のアルファベットＡ〜Ｎは特定のポリゴンに属する形式でなく、座標が異なる共有頂点を識別するために示している。大文字のアルファベットＡ〜Ｎは、ポリゴン群全体において共有頂点を識別するために、便宜上付されたものであり、後述する頂点座標配列データや、頂点構成テーブル、頂点座標テーブル、頂点重複テーブル、頂点更新テーブルには記録されていない（一部記載されているものは、理解を容易にするために用いられている。）。また、ポリゴンＩＤＰ１〜Ｐ１８には、各ポリゴンＩＤに対応付けて表現されたポリゴン別頂点ＩＤＶ１ａ〜Ｖ１８ｃが記載されており、各々は共有頂点Ａ〜Ｎのいずれかに対応している。ポリゴン別頂点ＩＤＶ１ａ〜Ｖ１８ｃは、説明の便宜上設けたものであり、頂点座標配列データ、頂点構成テーブル、頂点座標テーブル、頂点重複テーブル、頂点更新テーブルには記録されていない。
【００４７】
図１１は、図１０に示したポリゴン群を実現するための頂点座標配列データを示す図である。頂点座標配列データは、ポリゴンモデルを規定するデータであり、ポリゴンの集合体であるポリゴンモデルを配列構造で定義している。すなわち、図１１は、ＳＴＬデータより法線ベクトルの情報を削除した構造を示している。図１１に示すように、頂点座標配列データは、ポリゴンを識別するポリゴン識別情報であるポリゴンＩＤＰ１〜Ｐ１８に対応付けて、ポリゴンを構成する第１頂点、第２頂点、第３頂点の３つの頂点順に各座標値が記録されている。
【００４８】
ポリゴン別頂点ＩＤは、その頂点が属するポリゴンのＩＤと当該ポリゴンにおける頂点順が特定される形式でポリゴンと頂点を識別する識別情報である。具体的には、図１０に示すように、頂点を示すＶの後に、対応するポリゴンＩＤに含まれる数字、そして第１頂点、第２頂点、第３頂点の別を示すａ、ｂ、ｃのいずれかの小文字のアルファベットで表現する。例えば、図１１の１行目に示すように、ポリゴンＰ１の第１頂点のポリゴン別頂点ＩＤは、“Ｖ１ａ”と表現される。ポリゴン別頂点ＩＤ“Ｖ１ａ”に数字“１”が入っているため、ポリゴンＰ１の頂点であることが明確になっている。実際には、このポリゴン別頂点ＩＤは単なるシアリル番号で、この番号の順序で各頂点の座標値が頂点座標配列データに記録されているが、上記のポリゴン別頂点ＩＤでは記号を付加し、図１０との対応関係、頂点がどのポリゴンに所属するかを明らかになるように表記したものである。また、図１１の頂点座標配列データにおいて、各頂点の座標は、図１０で示したＡ〜Ｎのいずれかの共有頂点に対応するため、各頂点の座標値にＡ〜Ｎを識別するアルファベットの小文字を添えて表現している。例えば、図１１の１行目に示すように、ポリゴンＰ１の第１頂点Ｖ１ａの座標値は、図１０において共有頂点“Ｄ”に対応するので、“（Ｘｄ，Ｙｄ，Ｚｄ）”と表現される。各ポリゴンごとに、構成される３つの頂点に対応付けて３つの頂点の座標値が記録されている。上記のような規則に従って作成されるため、頂点座標配列データにおいては、同一座標の頂点が複数箇所に重複して記録されることになる。例えば、図１０に示した頂点Ｄは、６つのポリゴンＰ１，Ｐ２，Ｐ４，Ｐ５，Ｐ６，Ｐ７に共有されるため、座標値“（Ｘｄ，Ｙｄ，Ｚｄ）”は、頂点Ｖ１ａ，Ｖ２ａ，Ｖ４ｂ，Ｖ５ｃ，Ｖ６ｃ，Ｖ７ｃとして６箇所に記録される。
【００４９】
図１２は、ポリゴン削減処理の詳細を示すフローチャートである。まず、ポリゴンモデル構造化手段３０が、頂点構成テーブル、頂点座標テーブルを用意する（ステップＳ２１）。ステップＳ２１においては、頂点構成テーブル、頂点座標テーブルを用意することができれば、どのような手法を用いてもよいが、ポリゴンモデルの構造化を行って作成することもできる。ポリゴンモデルの構造化とは、ポリゴンモデルを実現するデータである頂点座標配列データを、頂点構成テーブルおよび頂点座標テーブルで構成される構造に変換することにより、ポリゴンモデルにおける頂点関係（トポロジー・位相情報）が座標値（数値データ）から分離され、あるポリゴンの頂点の修正により影響を受ける周辺のポリゴンの頂点の探索が容易になり、編集し易い構造に組み替えることを意味する。
【００５０】
図１３は、ステップＳ２１のポリゴンモデルの構造化の詳細を示すフローチャートである。ステップＳ２１のポリゴンモデルの構造化は、ポリゴンモデル構造化手段３０により行われる。ポリゴンモデル構造化手段３０は、ハッシュ値算出手段、座標値照合手段、頂点構成テーブル処理手段、頂点座標テーブル処理手段としての機能を備えている。
【００５１】
まず、頂点構成テーブル、頂点座標テーブル、頂点ＩＤの初期化を行う（Ｓ１０１）。頂点構成テーブルは、各ポリゴンに含まれる各頂点を、同一座標の頂点を同一値とする共有頂点の頂点ＩＤで表現し、ポリゴンを構成する各頂点を頂点ＩＤで表現してポリゴンと対応付けて記録したテーブルである。頂点ＩＤは、共有頂点を識別するための頂点識別情報であり、共有頂点を一意に特定することができれば、どのようなものを用いてもよいが、演算処理の効率化のため、本実施形態では、０から始まるシリアルな整数値を用いている。頂点座標テーブルは、頂点構成テーブルに付与された頂点ＩＤと座標値を対応付けて記録したテーブルである。頂点構成テーブル、頂点座標テーブルの初期化は、各項目が記録されていない状態にすることを意味する。また、頂点ＩＤの初期化として現在の頂点ＩＤ＝０に設定する。さらに、ハッシュ値とリンク先頂点ＩＤを対応付けた先頭位置テーブルも初期化する。先頭位置テーブルの初期化は、リンク先頂点ＩＤを全て“−１”に設定することにより行われる。頂点座標テーブルと先頭位置テーブルの２つのテーブルによりハッシュ値別頂点座標テーブルを実現する。
【００５２】
次に、頂点座標配列データより頂点座標値を抽出するにあたり、ポリゴンＩＤと頂点順を定義し、ポリゴンＩＤは“Ｐ１”、頂点順は“０”に初期化する（Ｓ１０１）。そして、現在のポリゴンＩＤと頂点順に対応する１つの頂点座標値を抽出する（Ｓ１０２）。例えば、図１１に示した頂点座標配列データからは、ポリゴンＩＤ“Ｐ１”と頂点順“０”である第１頂点の頂点座標値（Ｘｄ，Ｙｄ，Ｚｄ）が抽出される。
【００５３】
続いて、抽出された頂点座標値（Ｘｄ，Ｙｄ，Ｚｄ）よりハッシュ値を算出する（Ｓ１０３）。具体的には、まず、ポリゴンモデル内の全てのポリゴンの３次元座標値ｘ，ｙ，ｚの最大値Ｘｍａｘ，Ｙｍａｘ，Ｚｍａｘ、最小値Ｘｍｉｎ，Ｙｍｉｎ，Ｚｍｉｎを求める。これは、図１１に示したような頂点配列データの座標ごとの最大値、最小値を求めることになる。そして、座標値ｘ，ｙ，ｚの値を均等なＳ段階に分類し、整数値ｈｘ，ｈｙ，ｈｚ（０≦ｈｘ，ｈｙ，ｈｚ≦Ｓ−１）に変換する。整数Ｓとしては、例えばＳ＝３２を設定することができる。具体的には、以下の〔数式１〕に従った処理を実行することにより、ｈｘ，ｈｙ，ｈｚを算出する。ただし、〔数式１〕に従った処理の後、小数点以下を切り捨てた値を、整数値ｈｘ，ｈｙ，ｈｚとして得るものとする。
【００５４】
〔数式１〕
ｈｘ＝（ｘ−Ｘｍｉｎ）・Ｓ／（Ｘｍａｘ−Ｘｍｉｎ）
ｈｙ＝（ｙ−Ｙｍｉｎ）・Ｓ／（Ｙｍａｘ−Ｙｍｉｎ）
ｈｚ＝（ｚ−Ｚｍｉｎ）・Ｓ／（Ｚｍａｘ−Ｚｍｉｎ）
【００５５】
そして、以下の〔数式２〕に従った処理を実行することにより、ハッシュ値ｈを算出する。
【００５６】
〔数式２〕
ｈ＝ｈｘ＋ｈｙ・Ｓ＋ｈｚ・Ｓ2（０≦ｈ≦Ｓ3−１）
【００５７】
次に、算出されたハッシュ値を用いて先頭位置テーブルを参照する（Ｓ１０４）。そして、そのハッシュ値に対応するリンク先頂点ＩＤが先頭位置テーブルに登録されているか否かを判定する（Ｓ１０５）。ハッシュ値に対応するリンク先頂点ＩＤ（図中“リンク頂点ＩＤ”と表示）が“−１”、すなわち初期状態である場合は、そのハッシュ値に対応するリンク先頂点ＩＤが未登録、すなわち既登録でないことを意味する。
【００５８】
この場合、先頭位置テーブルに頂点ＩＤの記録を行う（Ｓ１０６）。具体的には、先頭位置テーブルのハッシュ値に対応するリンク先頂点ＩＤに、現在の頂点ＩＤ（初期状態では“０”）を書き込む。続いて、頂点座標テーブルの現在の頂点ＩＤの位置にステップＳ１０２で抽出した座標値を書き込む。この時、頂点座標テーブルにもリンク先頂点ＩＤを書き込む欄があるが、初期状態の“−１”のままにする。また、頂点構成テーブルにステップＳ１０２で抽出した現在のポリゴンＩＤと頂点順に対応付けて、現在の頂点ＩＤを書き込む。ここでは、新規な頂点ＩＤとして記録されることになる。そして、現在の頂点ＩＤをインクリメントする（Ｓ１０７）。
【００５９】
一方、ステップＳ１０５において、ハッシュ値に対応するリンク先頂点ＩＤが“−１”でない場合は、そのハッシュ値に対応するリンク先頂点ＩＤが既登録であることを意味する。この場合、ステップＳ１０２で抽出した座標値と、頂点座標テーブルにおける既登録のリンク先頂点ＩＤに対応する座標値が一致するか否かを判定する（Ｓ１０８）。判定の結果、両座標値が一致する場合は、頂点構成テーブルの現在のポリゴンＩＤ、頂点順の位置に、現在の頂点ＩＤを記録する（Ｓ１０９）。ステップＳ１０８における判定の結果、両座標値が異なる場合は、頂点座標テーブルのリンク先頂点ＩＤを参照し、リンク先頂点ＩＤが０以上であるか否かを判定する（Ｓ１１０）。判定の結果、リンク先頂点ＩＤが０以上である場合は、その値を頂点ＩＤとし、ステップＳ１０８に戻って、その頂点ＩＤに対応する座標値とステップＳ１０２で抽出した座標値を照合する。ステップＳ１１０における判定の結果、リンク先頂点ＩＤが負値（例えば“−１”）である場合は、頂点座標テーブルのその頂点ＩＤに対応するリンク先頂点ＩＤに現在の頂点ＩＤを記録し、ステップＳ１０７に戻って、頂点座標テーブルの現在の頂点ＩＤの位置にステップＳ１０２で抽出した座標値を記録する。この時、頂点座標テーブルの現在の頂点ＩＤにおけるリンク先頂点ＩＤについては、初期状態の“−１”のままにする。また、頂点構成テーブルにステップＳ１０２で抽出した現在のポリゴンＩＤと頂点順に対応付けて頂点ＩＤを記録する。そして、頂点ＩＤをインクリメントする。
【００６０】
ステップＳ１０７またはＳ１０９の処理を終えたら、全ポリゴンに対して処理を終えたか否かを判定する（Ｓ１１１）。具体的には、まず、頂点順をインクリメントし、頂点順が“２”を超えたら、ポリゴンＩＤをインクリメントして頂点順を“０”にする。例えば、現在のポリゴンＩＤが“Ｐ１”の場合、ポリゴンＩＤをインクリメントすると、ポリゴンＩＤは“Ｐ２”になる。ポリゴンＩＤがポリゴン総数を超えていなければ、ステップＳ１０２に戻って処理を繰り返す。ポリゴンＩＤがポリゴン総数を超えていれば、全ポリゴンに対して処理を終えたと判定して処理を終了する。
【００６１】
図１３に示したステップＳ２１のポリゴンモデルの構造化処理の具体例について説明する。図１０に示したポリゴン群が、その頂点の座標値により図１４に示すようなハッシュ値をとる場合を例にとって説明する。図１４の例では、説明の都合上２次元で表現し、上記〔数式２〕においてＳ＝２とすると、ｈ＝ｈｘ＋ｈｙ・２となり、頂点Ａ，Ｂ，Ｄ，Ｅのハッシュ値がｈ＝０、頂点Ｃ，Ｆのハッシュ値がｈ＝１、頂点Ｇ，Ｈ，Ｊ，Ｌ，Ｍのハッシュ値がｈ＝２、頂点Ｉ，Ｋ，Ｎのハッシュ値がｈ＝３であることを示している。実際には上記〔数式２〕において、Ｓ＝３２に設定し、ｈ＝ｈｘ＋ｈｙ・３２＋ｈｚ・１０２４となり、ハッシュ値ｈは、０〜３２７６７の値をとるが、以下では、説明を簡略化するため、ハッシュ値は前述の０〜３の値をとるものとして説明していく。
【００６２】
図１５は、Ｓ１０１における初期化時の各テーブルの状態を示す図である。図１５（ａ）は頂点構成テーブル、図１５（ｂ）は先頭位置テーブル、図１５（ｃ）は頂点座標テーブルを示している。先頭位置テーブルには、ハッシュ値に頂点ＩＤが対応付けられてリンク先頂点ＩＤとして記録される。ここでは、図１４に示したようにハッシュ値が０〜３の値をとるため、初期状態の先頭位置テーブルには、全てのハッシュ値０〜３に対応付けてリンク先頂点ＩＤとして“−１”になっている。頂点座標テーブルには、頂点ＩＤに対応付けられて頂点座標が記録されるとともに、リンク先頂点ＩＤが記録される。
【００６３】
図１５（ａ）に示すように、頂点構成テーブルは、ポリゴンＩＤに対応付けて頂点順に３つの頂点ＩＤが記録されるようになっているが、頂点順は、左から０，１，２となっており、図１１に示した第１頂点、第２頂点、第３頂点に対応している。ステップＳ１０１における初期化の際、頂点順＝０に初期化される。ステップＳ１０２において、頂点座標配列データより、ポリゴンＰ１の第１頂点の座標値（Ｘｄ，Ｙｄ，Ｚｄ）を抽出して、ステップＳ１０３において、ハッシュ値“０”が算出されたものとする。ステップＳ１０４において、ハッシュ値“０”で先頭位置テーブルを参照すると、ステップＳ１０５において、リンク先の頂点ＩＤが登録されていないため、先頭位置テーブルのハッシュ値“０”に対応付けて、ステップＳ１０６において、現在の頂点ＩＤ“０”をリンク先として記録する。
【００６４】
そして、ステップＳ１０７において、頂点座標テーブルにおける現在の頂点ＩＤ“０”の位置に、ステップＳ１０２において抽出した座標値（Ｘｄ，Ｙｄ，Ｚｄ）を記録し、リンク先頂点ＩＤとして“−１”を記録する。さらに、ステップＳ１０７において、現在のポリゴンＩＤ“Ｐ１”の現在の頂点順“０”の位置に、現在の頂点ＩＤ“０”を記録する。この結果、各テーブルは、図１６に示すような状態となる。そして、頂点ＩＤをインクリメントして現在の頂点ＩＤを“０”→“１”とする。また、頂点順をインクリメントして現在の頂点順を“０”→“１” とする。現在のポリゴンＩＤは“Ｐ１”であり、“Ｐ１３”を超えていないため、ステップＳ１１１において、全ポリゴンに対して処理終了していないと判断され、ステップＳ１０２に戻る。
【００６５】
ステップＳ１０２において、頂点座標配列データより、現在のポリゴンＩＤ“Ｐ１”、現在の頂点順“１”に基づき、ポリゴンＰ１の第２頂点の座標値（Ｘｂ，Ｙｂ，Ｚｂ）を抽出する。そして、ステップＳ１０３において、座標値（Ｘｂ，Ｙｂ，Ｚｂ）に対してハッシュ値“０”が算出されたものとする。ステップＳ１０４において、ハッシュ値“０”で先頭位置テーブルを参照すると、ステップＳ１０５において、リンク先の頂点ＩＤとして“０”が既登録であるため、ステップＳ１０８において、既登録の頂点ＩＤ“０”に対応する頂点座標テーブルの座標値（Ｘｄ，Ｙｄ，Ｚｄ）と、ステップＳ１０２で抽出した座標値（Ｘｂ，Ｙｂ，Ｚｂ）を照合する。照合の結果、不一致であるので、ステップＳ１１０において、頂点座標テーブルのリンク先の頂点ＩＤが０以上であるか負の値であるかを判定する。
【００６６】
図１６に示すように、リンク先の頂点ＩＤが“−１”であるので、ステップＳ１０７において、頂点座標テーブルの現在の頂点ＩＤ“１”の位置に、ステップＳ１０２において抽出した座標値（Ｘｂ，Ｙｂ，Ｚｂ）を記録し、リンク先頂点ＩＤとして“−１”を記録する。さらに、ステップＳ１０７において、頂点構成テーブルの、現在のポリゴンＩＤ“Ｐ１”の現在の頂点順“１”の位置に、現在の頂点ＩＤ“１”を記録する。この結果、各テーブルは、図１７に示すような状態となる。そして、頂点ＩＤをインクリメントして現在の頂点ＩＤを“１”→“２”とする。また、頂点順をインクリメントして現在の頂点順を“１”→“２” とする。現在のポリゴンＩＤは“Ｐ１”であり、“Ｐ１３”を超えていないため、ステップＳ１１１において、全ポリゴンに対して処理終了していないと判断され、ステップＳ１０２に戻る。
【００６７】
ポリゴンＩＤ“Ｐ１”、頂点ＩＤ“２”、頂点順“２”の場合は、ポリゴンＩＤ“Ｐ１”、頂点ＩＤ“１”、頂点順“１”の場合と同様に処理され、各テーブルは、図１８に示すような状態となる。そして、頂点ＩＤをインクリメントして現在の頂点ＩＤを“２”→“３”とする。また、頂点順は、１つのポリゴンについて“０”“１”“２”の３つの値のみをとるため、ポリゴンＩＤをインクリメントして“Ｐ１”→“Ｐ２”とし、頂点順をリセットして現在の頂点順を“２”→“０” とする。現在のポリゴンＩＤは“Ｐ２”となり、“Ｐ１８”を超えていないため、ステップＳ１１１において、全ポリゴンに対して処理終了していないと判断され、ステップＳ１０２に戻る。
【００６８】
ステップＳ１０２において、頂点座標配列データより、現在のポリゴンＩＤ“Ｐ２”、現在の頂点順“０”に基づき、ポリゴンＰ２の第１頂点の座標値（Ｘｄ，Ｙｄ，Ｚｄ）を抽出する。そして、ステップＳ１０３において、座標値（Ｘｄ，Ｙｄ，Ｚｄ）に対してハッシュ値“０”が算出されたものとする。ステップＳ１０４において、ハッシュ値“０”で先頭位置テーブルを参照すると、ステップＳ１０５において、リンク先の頂点ＩＤとして“０”が既登録であるため、ステップＳ１０８において、既登録の頂点ＩＤ“０”に対応する頂点座標テーブルの座標値（Ｘｄ，Ｙｄ，Ｚｄ）と、ステップＳ１０２で抽出した座標値（Ｘｄ，Ｙｄ，Ｚｄ）を照合する。照合の結果、一致するので、ステップＳ１０９において、頂点座標テーブルの座標値が一致した頂点ＩＤ“０”を、頂点構成テーブルの、現在のポリゴンＩＤ“Ｐ２”の現在の頂点順“０”の位置に記録する。この結果、各テーブルは、図１９に示すような状態となる。
【００６９】
上記のようにして、ポリゴンＰ１８まで処理を行ってステップＳ１１１の判定により処理を終了する。このとき、各テーブルは、図２０に示すような状態となる。すなわち、ハッシュ値が同一の他の座標値を頂点ＩＤより最近傍に記録した単一の頂点ＩＤがリンク先として対応付けられているものとなる。上記説明のように、ポリゴンモデルの構造化処理においては、ステップＳ１０８において座標値の照合を行い、不一致であれば頂点座標テーブルに追加登録するが、一致している場合は、次々に照合を繰り返していく。この処理は、頂点座標テーブルが大きくなり、既登録の座標値が多くなる程、負荷が高くなることになる。本実施形態では、ハッシュ値を用いてハッシュ値別頂点座標テーブルを参照することにより、頂点座標テーブルの全ての座標値の照合を行うことなく、ハッシュ値が異なる座標値とは一致しないことが自明なため照合対象から除外することができ、同一ハッシュ値の頂点座標テーブルのみを参照すればよいことになる。頂点座標テーブルの中で同一ハッシュ値をもつ座標値は、本実施形態では高々１／４程度であるが、前述の通り実際には上記〔数式２〕において、Ｓ＝３２に設定することが多く、そうすると頂点座標テーブルの中で同一ハッシュ値をもつ座標値は１／３２７６７程度に削減され、全件と照合する場合に比べ照合時間も１／３２７６７程度に短縮されることになる。このため、ポリゴンモデルの構造化処理を桁違いに高速に行うことが可能となる。
【００７０】
図２０の例で、ハッシュ値が“０”の場合、図２０（ｃ）の頂点座標テーブルでは、リンク先頂点ＩＤに従って、頂点ＩＤ“０”“１”“２”“６”のものだけ照合すれば、新たな頂点座標を追加登録するべきか否かの判断が可能となる。図２０（ｂ）に示した先頭位置テーブルを削除し、図２０（ｃ）に示した頂点座標テーブルよりリンク先頂点ＩＤの欄を削除することにより、図２１に示したような頂点構成テーブルと頂点座標テーブルを得ることができ、ポリゴンモデルの構造化は完了する。
【００７１】
図２１は、ステップＳ２１において用意された頂点構成テーブルと頂点座標テーブルを示す図である。図２１（ａ）に示す頂点構成テーブルは、ポリゴンを構成する各頂点をポリゴンモデル内で一意に識別する共有頂点と対応付け、共有頂点の識別情報である頂点ＩＤを、ポリゴンを識別するポリゴン識別情報であるポリゴンＩＤに対応付けて記録されたテーブルである。本実施形態ではポリゴンが三角形であるため、図２１（ａ）の例では、１つのポリゴンＩＤに３つの頂点ＩＤが対応付けられている。図２１（ｂ）に示す頂点座標テーブルは、共有頂点の識別情報である頂点ＩＤと、ポリゴンモデル内における頂点のｘｙｚ座標値である頂点座標が対応付けて記録されたテーブルである。
【００７２】
次に、テーブル定義手段４０が、頂点重複テーブルを初期化する（ステップＳ２２）。頂点重複テーブルは、頂点ＩＤと重複先の頂点ＩＤ（図中「重複先」と示す）を対応付けたテーブルである。図２１（ｃ）に頂点重複テーブルの一例を示す。頂点重複テーブルには、初期値として、重複先がないことを示す“−１”が設定される。重複先がないことを示すことができれば、初期値として、“−１”以外を設定してもよい。
【００７３】
次に、テーブル定義手段４０が、頂点更新テーブルを初期化する（ステップＳ２３）。頂点更新テーブルは、頂点ＩＤと更新フラグ（図中「更新」と示す）を対応付けたテーブルである。図２１（ｄ）に頂点更新テーブルの一例を示す。頂点更新テーブルには、初期値として、更新されていないことを示す“０”が設定される。更新されていないことを示すことができれば、初期値として、“０”以外を設定してもよい。
【００７４】
次に、削除対象ポリゴン設定手段５０が、ポリゴンモデルを構成する全ポリゴンを対象として、各ポリゴンの面積および周長を算出し、算出された面積と周長の積の最小値を算出する（ステップＳ２４）。すなわち、全てのポリゴンについて、面積と周長を乗じて積を算出し、面積と周長の積が各ポリゴンの中で最小となる値を特定する。
【００７５】
次に、削除対象ポリゴン設定手段５０は、ステップＳ２４で算出された最小値に所定の係数αを乗算した値を所定のしきい値とし、面積と周長の積が所定のしきい値未満であって、３頂点のうち１つの共有頂点も更新されていないポリゴンを削除対象ポリゴンとし、削除対象ポリゴンの最短稜線（最短辺）の２端点をその中点に統合する処理を行う（ステップＳ２５）。所定の係数αとしては、１以上の値をとる実数値に設定し、１．５≦α≦２．５の範囲であることが好ましい。本実施形態では、α＝２．０としている。３頂点のうち１つの共有頂点も更新されていないとは、削除対象に設定しようとするポリゴンのいずれの頂点と共有する頂点をもつ他の隣接ポリゴンが削除対象に事前に設定されていないことを意味し、実際に削除される最短稜線の２端点である２つの頂点だけでなく残りのもう１つの頂点も更新されていないことを確認する必要があり、隣接する複数のポリゴンを同時に削除対象に設定することにより、削除後のポリゴンモデルを歪ませることを防いでいる。具体的には、頂点更新テーブルを参照して、ポリゴンＰ６については、頂点ＩＤ“６”“９”“１０”に対応する更新フラグが全て“０”であることを確認して、削除対象ポリゴンとする。このように、頂点更新テーブルを参照することにより、ポリゴンの全ての頂点が、既に削除対象に設定された全てのポリゴンの削除対象稜線の２頂点と共有しないことを確認することができる。削除対象稜線の２頂点は、削除対象頂点となる。
【００７６】
ステップＳ２５による処理において、図１０に示したポリゴン群のうち、面積と周長の積がしきい値未満であるポリゴンが、ポリゴンＰ６とポリゴンＰ１０の２つであったとする。そして、図２２に示すように、ポリゴンＰ６の最短稜線は稜線ＤＥ、ポリゴンＰ１０の最短稜線は稜線ＪＫであったとする。この場合、２つの最短稜線のうち、短い方の最短稜線である稜線ＪＫから処理を行う。そして、稜線ＪＫを有するポリゴンＰ１０の３頂点が更新されていないか否かを、頂点更新テーブルを参照して確認する。図２１に示すように、初期状態では、ポリゴンＰ１０の３頂点の更新フラグは、いずれも“０”であるため、ポリゴンＰ１０は、削除対象ポリゴンとなる。この場合、２頂点Ｊ，Ｋをそれぞれ中点Ｏに統合する。
【００７７】
次に、頂点座標テーブル更新手段６０、頂点重複テーブル更新手段７０、頂点更新テーブル更新手段７５が、頂点座標テーブル、頂点重複テーブル、頂点更新テーブルを更新する（ステップＳ２６）。具体的には、稜線ＪＫの中点をＯ（Ｘｏ，Ｙｏ，Ｚｏ）とし、頂点座標テーブルにおいて、頂点Ｊに対応していた頂点ＩＤ“９”の頂点の座標値（Ｘｊ，Ｙｊ，Ｚｊ）を中点Ｏの座標値（Ｘｏ，Ｙｏ，Ｚｏ）に変更する。また、頂点座標テーブルにおいて、頂点Ｋに対応していた頂点ＩＤ“１０”の頂点の座標値（Ｘｋ，Ｙｋ，Ｚｋ）を削除する。
【００７８】
頂点ＩＤ“１０”の頂点は、頂点ＩＤ“９”の頂点の座標（Ｘｏ，Ｙｏ，Ｚｏ）と同一であり、重複することになる。このため、頂点重複テーブルにおいて、頂点座標テーブルから削除された頂点ＩＤ“１０”の重複先を頂点ＩＤ“９”に設定する。また、頂点更新テーブルにおいて、所定の頂点について、頂点重複テーブルにおいて更新が行われた旨を記録するため、頂点重複テーブルが更新された頂点ＩＤ“９”、“１０”の更新フラグを“１”に設定する。これは、その頂点ＩＤで特定される頂点が、削除対象頂点に設定されたことを意味している。この結果、各テーブルは、図２１に示した初期状態から図２３に示したような状態に変更される。
【００７９】
そして、ステップＳ２５に戻り、２つの最短稜線のうち、長い方の最短稜線である稜線ＤＥの処理を行う。具体的には、稜線ＤＥを有するポリゴンＰ６の３頂点が更新されていないか否かを、頂点更新テーブルを参照して確認する。図２３に示すように、ポリゴンＰ１０の稜線ＪＫの処理後においても、ポリゴンＰ６の３頂点の更新フラグは、いずれも“０”であるため、ポリゴンＰ６は、削除対象ポリゴンとなる。この場合、２頂点Ｄ，Ｅをそれぞれ中点Ｐに統合する。
【００８０】
そして、ステップＳ２６において、頂点座標テーブル、頂点更新テーブル、頂点重複テーブルを更新する。具体的には、稜線ＤＥの中点をＰ（Ｘｐ，Ｙｐ，Ｚｐ）とし、頂点座標テーブルにおいて、頂点Ｄに対応していた頂点ＩＤ“０”の頂点の座標値（Ｘｄ，Ｙｄ，Ｚｄ）を中点Ｐの座標値（Ｘｐ，Ｙｐ，Ｚｐ）に変更する。また、頂点座標テーブルにおいて、頂点Ｅに対応していた頂点ＩＤ“６”の頂点の座標値（Ｘｅ，Ｙｅ，Ｚｅ）を削除する。
【００８１】
頂点ＩＤ“６”の頂点は、頂点ＩＤ“０”の頂点の座標（Ｘｐ，Ｙｐ，Ｚｐ）と同一であり、重複することになる。このため、頂点重複テーブルにおいて、頂点座標テーブルから削除された頂点ＩＤ“６”の重複先を頂点ＩＤ“０”に設定する。また、頂点更新テーブルにおいて、頂点重複テーブルが更新された頂点ＩＤ“０”、“６”の更新フラグを“１”に設定する。この結果、各テーブルは、図２３に示した状態から図２４に示したような状態に変更される。
【００８２】
全ての削除対象ポリゴンについて処理を終えたら、頂点構成テーブル更新手段８０が、更新された頂点重複テーブルに基づいて頂点構成テーブルを更新し、頂点が重複するポリゴンを削除する（ステップＳ２７）。まず、頂点重複テーブルに基づき、頂点構成テーブルを更新する。具体的には、頂点重複テーブルにおいて重複先が設定されている頂点ＩＤについて、頂点構成テーブルの頂点ＩＤを重複先の頂点ＩＤに置き換える。例えば、図２５（ｃ）に示す頂点重複テーブルでは、頂点ＩＤ“６”の重複先が“０”、頂点ＩＤ“１０”の重複先が“９”であるので、頂点構成テーブルに記録されている全ての頂点ＩＤ“６”を頂点ＩＤ“０”に置き換え、全ての頂点ＩＤ“１０”を頂点ＩＤ“９”に置き換える。この結果、頂点構成テーブルは、図２５（ａ）に示すように、頂点ＩＤ“６”と頂点ＩＤ“１０”が存在しない状態に更新される。
【００８３】
次に、頂点が重複するポリゴンを削除する。具体的には、頂点構成テーブルにおいて、１つのポリゴンに対応する３つの頂点のうち、２つの頂点の頂点ＩＤが同一である場合に、そのポリゴンの情報を頂点構成テーブルから削除する。図２５（ａ）の例では、ポリゴンＰ５では頂点ＩＤ“０”、ポリゴンＰ６では頂点ＩＤ“０”、ポリゴンＰ１０では頂点ＩＤ“９”、ポリゴンＰ１６では頂点ＩＤ“９”がそれぞれ重複している。このため、ポリゴンＰ５、Ｐ６、Ｐ１０、Ｐ１６の情報を頂点構成テーブルから削除する。また、頂点座標テーブルにおいて、頂点座標が記録されていない頂点に関する情報を削除する。この結果、頂点構成テーブルは、図２６（ａ）に示すような状態となる。
【００８４】
頂点座標配列データが更新されたら、制御手段１０が、削減後の全ポリゴン数が、設定された目標値以下となっているか否かを判定する（ステップＳ２８）。目標値としては、事前に任意の値を設定することができる。判定の結果、全ポリゴン数が目標値以下となっていない場合は、ステップＳ２２に戻って、頂点重複テーブルを再び初期化し、ステップＳ２３において、頂点更新テーブルを再び初期化する。この結果、頂点重複テーブル、頂点更新テーブルは、それぞれ図２６（ｃ）、図２６（ｄ）に示すような状態となる。頂点重複テーブルと頂点更新テーブルを再初期化した後、ステップＳ２４において、削減後に残っている全てのポリゴンに対して、既に算出されている各ポリゴンの面積および周長を基に、面積と周長の積の最小値を再度算出し、新たに算出された最小値に前記係数αを乗算した所定のしきい値を用いて、削除対象ポリゴン設定手段５０が、残っているポリゴンの中から、他のポリゴンを削除対象ポリゴンとして選定する。
【００８５】
ステップＳ２８における判定の結果、全ポリゴン数が目標値以下である場合は、ポリゴン削減処理を終了する。図１２のフローチャートに従った処理を実行することにより、各グループについて設定された目標値以下になるように、ポリゴン数が削減されることになる。ポリゴンモデルの構造化処理が行われていなければ、削除対象ポリゴンの削除対象辺の２つの補正対象頂点と同一座標となる頂点の探索に大きな処理負荷を要することになる。ポリゴンモデルの構造化処理を行うことにより、処理負荷を大幅に削減することが可能となる。
【００８６】
図１２に示した処理を終えたら、最後に、ポリゴンモデル更新手段９０が、頂点座標配列データを更新する。具体的には、頂点座標テーブルを参照して、頂点構成テーブルの各頂点ＩＤを頂点座標に置き換える。図２６（ｄ）に示した頂点構成テーブル、図２６（ｂ）に示した頂点座標テーブルからは図２７に示したような頂点座標配列データが得られることになる。これを図１１に示した元の頂点座標配列データと置き換える。この結果、図１０、図２２に示したポリゴン群は、図２８に示すようなポリゴン群に変更されることになる。図１０、図２２に示したポリゴン群において設定された削除対象ポリゴンＰ６，Ｐ１０だけでなく、ポリゴンＰ５，Ｐ１６が削除されていることがわかる。これにより、ステップＳ２０のポリゴン削減処理が終了する。
【００８７】
＜２．５．小体積グループ削除処理＞
次に、ステップＳ４０におけるグループ合成手段１００による隔離された小体積グループの削除処理について説明する。ここで、隔離された小体積グループの削除処理の説明に先立ち、隔離された小体積グループの削除処理を行う理由について説明しておく。図２９は、ポリゴン削減処理により生じる問題を説明するための図である。図２９（ａ）は、ステップＳ１０のグループ分類処理を終えた直後の２グループのポリゴン群を示している。図２９（ａ）の例では、頂点Ｄ、Ｅ、Ｋを含む左側のグループと、頂点Ｄ’、Ｅ’、Ｋ’を含む右側のグループが存在する。頂点ＤとＤ’、頂点ＥとＥ’、頂点ＫとＫ’がそれぞれ３Ｄプリンタの解像度以下である場合、この２グループのポリゴン群を３Ｄプリンタで出力すると、頂点ＤとＤ’、頂点ＥとＥ’、頂点ＫとＫ’は樹脂等により連結され、２グループは連結した物体として出力される。すなわち、２グループは物理的には接着されている。
【００８８】
図２９（ａ）に示した右側のグループのポリゴン群に対して、ステップＳ２０におけるポリゴン削減処理を実行すると、図２９（ｂ）に示すように、ポリゴンＥ’ＦＫ’が削除対象ポリゴンとして選定される。さらに、図２９（ｃ）に示すように、削除対象ポリゴンＥ’ＦＫ’の最短稜線Ｅ’Ｆが削除対象辺として選定され、辺Ｅ’Ｆの中点Ｏが算出される。そして、図２９（ｄ）に示すように、削除対象ポリゴンＥ’ＦＫ’と、削除頂点共有ポリゴンＤ’Ｅ’Ｆが削除される。
【００８９】
図２９（ａ）と図２９（ｄ）を比較するとわかるように、左右の２つのグループは、ポリゴン削減処理前に比べて、ポリゴン削減処理後には、より離れてしまっている。図２９（ｄ）に示した２グループの距離が、３Ｄプリンタの解像度を超えると、２グループのポリゴン群を３Ｄプリンタで出力した場合、２グループは別の物体として出力される。２つのグループがともに大きな物体となれば、問題はない。しかし、一つが非常に小さな物体となる場合、例えば、一方のグループが本体であって、他方のグループが本体に付属するヒゲのようなものである場合、本体と別に出力されたヒゲは、３Ｄプリンタ内に紛れ込み、機械詰まり等の問題を引き起こす恐れがある。このような問題を除去するため、上述のヒゲのような隔離された小体積グループの削除処理を行うのである。
【００９０】
図３０は、グループ合成手段１００による隔離された小体積グループの削除処理（ステップＳ４０）の詳細を示すフローチャートである。グループ合成手段１００は、まず、グループごとに全ポリゴンの全頂点を対象として、Ｘ、Ｙ、Ｚぞれぞれの方向の最大値および最小値を算出する（ステップＳ４１）。この結果、グループｇに対して、最小値Ｘｍｉｎ（ｇ）、最大値Ｘｍａｘ（ｇ）、最小値Ｙｍｉｎ（ｇ）、最大値Ｙｍａｘ（ｇ）、最小値Ｚｍｉｎ（ｇ）、最大値Ｚｍａｘ（ｇ）が得られる。次に、グループ合成手段１００は、グループごとに外接直方体の体積を算出する（ステップＳ４２）。具体的には、以下の〔数式３〕に従った処理を実行することにより、グループｇの外接直方体の体積Ｖ（ｇ）を算出する。
【００９１】
〔数式３〕
Ｖ（ｇ）＝（Ｘｍａｘ（ｇ）−Ｘｍｉｎ（ｇ））×（Ｙｍａｘ（ｇ）−Ｙｍｉｎ（ｇ））×（Ｚｍａｘ（ｇ）−Ｚｍｉｎ（ｇ））
【００９２】
グループ合成手段１００は、全グループに対してステップＳ４１、Ｓ４２の処理を実行する。これにより、全グループについて、外接直方体の体積が得られる。次に、グループ合成手段１００は、全グループの中で外接直方体の体積が最小となるグループを特定する（ステップＳ４３）。そして、グループ合成手段１００は、特定したグループと他の全てのグループとの最小間隔を算出する（ステップＳ４４）。ステップＳ４４における最小間隔の算出処理の詳細については後述する。
【００９３】
続いて、グループ合成手段１００は、ステップＳ４４において算出された最小間隔が所定のしきい値以上か否かを判定する（ステップＳ４５）。最小間隔と比較するための所定のしきい値としては、３Ｄプリンタ出力時に０．２ｍｍに相当する数値とすることが好ましい。判定の結果、最小間隔が所定のしきい値以上である場合は、体積最小のグループが、最も近いグループとも所定の距離以上離れていることになるので、グループ合成手段１００は、体積最小のグループを削除する処理を行う（ステップＳ４６）。ステップＳ４５における判定の結果、最小間隔が所定のしきい値未満である場合は、体積最小のグループが、最も近いグループと３Ｄプリンタの解像度以下に近接していることになるので、グループ合成手段１００は、体積最小のグループを削除する処理を行わない。
【００９４】
ステップＳ４６におけるグループの削除を行った場合、削除を行わなかった場合のいずれにおいても、そのグループについては処理済みとする。そして、ステップＳ４３に戻って、未処理のグループの中で外接直方体の体積が最小となるグループを特定する。ステップＳ４３〜ステップＳ４６の処理を繰り返し、全てのグループに対する処理を終えたら、隔離された小体積グループの削除処理を終了する。
【００９５】
ステップＳ４４における最小間隔の算出処理の詳細について説明する。図３１は、図３０のステップＳ４４における最小間隔の算出処理の詳細を示すフローチャートである。グループ合成手段１００は、初期値として、最小間隔を示す変数ｄｍｉｎに所定以上の値を設定する。ｄｍｉｎに与える初期値としては、グループ間の最小間隔としてあり得ない大きな値を設定しておけばよい。そして、グループ合成手段１００は、全グループの中から、グループｇより体積が大きい１つのグループｈを抽出する（ステップＳ４７）。この結果、グループｈについて、最小値Ｘｍｉｎ（ｈ）、最大値Ｘｍａｘ（ｈ）、最小値Ｙｍｉｎ（ｈ）、最大値Ｙｍａｘ（ｈ）、最小値Ｚｍｉｎ（ｈ）、最大値Ｚｍａｘ（ｈ）が得られる。次に、グループ合成手段１００は、Ｘ、Ｙ、Ｚの各方向別に２グループ間の差分の最小値を算出する（ステップＳ４８）。具体的には、以下の〔数式４〕に従った処理を実行することにより、Ｘ方向の差分の最小値ｄｘ、Ｙ方向の差分の最小値ｄｙ、Ｚ方向の差分の最小値ｄｚを算出する。
【００９６】
〔数式４〕
ｄｘ＝Ｍｉｎ｛｜Ｘｍａｘ（ｇ）−Ｘｍｉｎ（ｈ）｜，｜Ｘｍｉｎ（ｇ）−Ｘｍａｘ（ｈ）｜｝
ｄｙ＝Ｍｉｎ｛｜Ｙｍａｘ（ｇ）−Ｙｍｉｎ（ｈ）｜，｜Ｙｍｉｎ（ｇ）−Ｙｍａｘ（ｈ）｜｝
ｄｚ＝Ｍｉｎ｛｜Ｚｍａｘ（ｇ）−Ｚｍｉｎ（ｈ）｜，｜Ｚｍｉｎ（ｇ）−Ｚｍａｘ（ｈ）｜｝
【００９７】
続いて、グループ合成手段１００は、算出された方向別最小値のうち最大のものをグループ間最小値とし、グループ間最小値を、既に処理済みの他のグループ間における全グループ間最小値と比較し、小さい方を全グループ間最小値とする処理を行う（ステップＳ４９）。具体的には、ステップＳ４８において算出された３つの方向別最小値ｄｘ、ｄｙ、ｄｚのうち、最大の値をとるグループ間最小値Ｍａｘ（ｄｘ、ｄｙ、ｄｚ）とｄｍｉｎを比較し、Ｍａｘ（ｄｘ、ｄｙ、ｄｚ）＜ｄｍｉｎの場合、ｄｍｉｎ＝Ｍａｘ（ｄｘ、ｄｙ、ｄｚ）とする。これにより、ｄｍｉｎには、グループ間最小値の中で最小の値となる全グループ間最小値が記録されることになる。
【００９８】
ステップＳ４７、Ｓ４８、Ｓ４９の処理を繰り返して実行し、グループｇより体積が大きい全てのグループに対する処理を終えたら、図３１に示した最小間隔の算出処理を終了し、図３０のステップＳ４５に進んで、最小間隔ｄｍｉｎとしきい値との比較を行うことになる。
【００９９】
＜２．６．ポリゴンモデルの合成＞
図９の処理概要に示したように、ステップＳ５０においては、グループ合成手段１００が、全グループのポリゴンモデルの合成を行って更新された領域別ポリゴンモデルを得る。また、ステップＳ７０においては、ポリゴンモデル合成手段１１０が、全領域の領域別ポリゴンモデルの合成を行って合成ポリゴンモデルを得る。いずれの合成処理も、ヘッダ部を除く領域別ポリゴンモデルを単純に結合し、ヘッダ部に合算されたポリゴン数を記録することにより行われる（ただし、これはバイナリ形式のＳＴＬフォーマットの場合で、ＡＳＣＩＩ形式のＳＴＬフォーマットではヘッダ部を除く領域別ポリゴンモデルを単純に結合するだけである）。合成ポリゴンモデル記憶手段１３０に記憶された合成ポリゴンモデルを３Ｄプリンタ７に出力することにより、３Ｄプリンタ７は効率的に立体物を造形することができる。造形された立体物については、ほとんど劣化は見られなかった。
【０１００】
＜３．変形例等＞
以上、本発明の好適な実施形態について説明したが、本発明は上記実施形態に限定されず、種々の変形が可能である。例えば、上記実施形態では、処理対象のポリゴンを三角形としたが、四角形以上の多角形であってもよい。
【０１０１】
また、上記実施形態では、ポリゴンに対して面積と周長を算出し、面積と周長の積が所定のしきい値未満となるポリゴンを削除対象ポリゴンとして選定するようにしたが、面積が所定のしきい値未満となるポリゴンを削除対象ポリゴンとして選定するようにしてもよい。この場合、所定のしきい値としては、面積の最小値に所定の係数を乗算した値を用いることができ、所定の係数としては、１以上の値をとる係数αを用いることができる。係数αは、１．５≦α≦２．５の範囲であることが好ましい。
【０１０２】
また、上記実施形態では、複数のグループの更新されたポリゴンモデルに対して、外接直方体の体積が所定のしきい値未満であるものを削除するようにしたが、削除しないようにすることも可能である。
【符号の説明】
【０１０３】
１・・・ＣＰＵ（Central Processing Unit）
２・・・ＲＡＭ（Random Access Memory）
３・・・記憶装置
４・・・キー入力Ｉ／Ｆ
５・・・データ入出力Ｉ／Ｆ
６・・・表示部
７・・・３Ｄプリンタ（立体物造形装置）
７ａ・・・データ処理部
７ｂ・・・出力部
１０・・・制御手段
２０・・・グループ分類手段
３０・・・ポリゴンモデル構造化手段
４０・・・テーブル定義手段
５０・・・削除対象ポリゴン設定手段
６０・・・頂点座標テーブル更新手段
７０・・・頂点重複テーブル更新手段
７５・・・頂点更新テーブル更新手段
８０・・・頂点構成テーブル更新手段
９０・・・ポリゴンモデル更新手段
１００・・・グループ合成手段
１１０・・・ポリゴンモデル合成手段
１２０・・・領域別ポリゴンモデル記憶手段
１３０・・・合成ポリゴンモデル記憶手段
２００・・・立体物造形用データ削減装置
【図１】
0006565494.tif 000002
【図２】
0006565494.tif 000003
【図３】
0006565494.tif 000004
【図４】
0006565494.tif 000005
【図５】
0006565494.tif 000006
【図６】
0006565494.tif 000007
【図７】
0006565494.tif 000008
【図８】
0006565494.tif 000009
【図９】
0006565494.tif 000010
【図１０】
0006565494.tif 000011
【図１１】
0006565494.tif 000012
【図１２】
0006565494.tif 000013
【図１３】
0006565494.tif 000014
【図１４】
0006565494.tif 000015
【図１５】
0006565494.tif 000016
【図１６】
0006565494.tif 000017
【図１７】
0006565494.tif 000018
【図１８】
0006565494.tif 000019
【図１９】
0006565494.tif 000020
【図２０】
0006565494.tif 000021
【図２１】
0006565494.tif 000022
【図２２】
0006565494.tif 000023
【図２３】
0006565494.tif 000024
【図２４】
0006565494.tif 000025
【図２５】
0006565494.tif 000026
【図２６】
0006565494.tif 000027
【図２７】
0006565494.tif 000028
【図２８】
0006565494.tif 000029
【図２９】
0006565494.tif 000030
【図３０】
0006565494.tif 000031
【図３１】
0006565494.tif 000032