(57)【特許請求の範囲】
【請求項１】
コンピュータグラフィックまたはリバースエンジニアリングにおいて自由曲面を表す曲面式を求める方法であって、該自由曲面を測定することによって得られた該自由曲面を示す測定データを、コンピュータを用いて点群データに分割し、該点群データにおける５点ごとの点列を４階差分商を用いて平滑化し、該平滑化された５点ごとの点列からなる点群データの各中央点の３次パラメトリック式を求め、該３次パラメトリック式に基づいて求めた５点ごとの点列からなる局所的点列データの中央点での曲面の１階、２階微分ベクトルを求め、これら微分ベクトルを４隅の微分ベクトルとする５次ベジエ曲面を生成することにより前記曲面式を得ることを特徴とする方法。
【請求項２】
前記点群データは、４辺格子上の格子点データであり、前記平滑化は、５点ごとの点列に対する局所的平滑基準値を、前記４辺格子の一方をｕ方向、他の一方をｖ方向とする空間において、前記４辺格子を構成するｕ方向のパラメータ線とｖ方向のパラメータ線に対して求めた前記４階差分商の合計とし、前記点群データ全体における全体平滑基準値を全格子点について求めた前記局所的平滑基準値の和とし、前記ｕ方向およびｖ方向のパラメータ線を同時に平滑化するものである請求項１に記載の方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、自由曲面を表す曲面式を求める方法に関するものであり、特に、複雑な自由曲面から構成された意匠形状を計測した測定データに曲面式を当てはめる方法に関するものである。
【０００２】
【従来の技術】
測定技術や情報技術の進展に伴って、自動車や家電製品を中心に意匠形状を有する製品開発のために、作成した立体モデルから計算機の内部モデルを作成するリバースエンジニアリング技術が利用されている（例えば、青山英紀：形状モデリングのためのリバースエンジニアリング、設計工学，３２，１０（１９９７）３８１．、およびＴ．  Ｖａｒａｄｙ、Ｒ．  Ｍａｒｔｉｎ  ａｎｄ  Ｊ．  Ｃｐｘ：Ｒｅｖｅｒｓｅ  Ｅｎｇｉｎｅｅｒｉｎｇ  ｏｆ  Ｇｅｏｍｅｔｒｉｃ  Ｍｏｄｅｌｓ    Ａｎ  Ｉｎｔｒｏｄｕｃｔｉｏｎ、Ｃｏｍｐｕｔｅｒ  Ａｉｄｅｄ  Ｄｅｓｉｇｎ、２９，４（１９９７）  ２５５．等参照）。。これらの形状は複雑な自由曲面より成り、また、測定データは大量でかつ誤差を含んでいるが（例えば、後藤孝行、三好隆志、高谷裕浩、枝光毅彦：曲率を考慮したＢ−スプライン曲線の測定データへのあてはめ、精密工学会誌，６０，７（１９９４）９６４、および倉賀野哲造、斉藤勝、黒田満、古川進：測定点データに基づ自由曲線・曲面のモデル化法、精密工学会誌、６５，３（１９９９）３８１．等参照）、現在市販されているシステムでは、そのまま金型を製作できるような曲率変化の滑らかな高品質のモデルを作成することは難しかった。
【０００３】
一方、コンピュータグラフィックス（ＣＧ）の分野では、これほどの精度を必要とせず、多面体での表示を目的としていることが多いため、三角形格子を中心とした技術が研究されている（例えば、鈴木宏正：ＣＧにおける形状再構成技術、設計工学，３２，１０（１９９７）３９２．等参照）。この中でも、測定誤差を含んだデータを処理するため、多量のデータを処理するため、多量のデータを取り扱う技術と共に、誤差の平準化の処理も研究されているが（例えば、Ｍ．  Ｅｃｋ  ａｎｄ  Ｈ．  Ｈｏｐｐｅｒ：  Ａｕｔｏｍａｔｉｃ  Ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ  ｏｆ  Ｂ−Ｓｐｌｉｎｅ  Ｓｕｒｆａｃｅ  ｏｆ  Ａｒｂｉｔｒａｒｙ  Ｔｏｐｏｌｏｇｉｃａｌ  Ｔｙｐｅ、Ｐｒｏｃ．  ＡＣＭ  ＳＩＧＧＲＡＰＨ’９６、ＡＣＭ  Ｐｒｅｓｓ（１９９６）３２５．、神島泰章、鈴木宏正、金井崇、木村文彦：測定点群からなる三角形メッシュ生成（第２報）、精密工学会誌、６４，１０（１９９８）１４６１．および鈴木宏正、神島泰章、金井崇、木村文彦：測定点群からなる三角形メッシュ生成（第１報）、精密工学会誌、６４，９（１９９８）１３１４．等参照）高品質な曲面生成を狙っていない。
【０００４】
これに対し、意匠形状に対する自由曲面を取り扱う計算機支援形状設計の分野では、古くは穂坂氏の論文（Ｍ．  Ｈｏｓａｋａ：  Ｔｈｅｏｒｙ  ｏｆ  Ｃｕｒｖｅｓ  ａｎｄ  Ｓｕｒｆａｃｅ  Ｓｙｎｔｈｅｓｉｓ  ａｎｄ  Ｔｈｅｉｒ  Ｓｍｏｏｔｈ  ｆｉｔｔｉｎｇ、Ｉｎｆｏｒｍａｔｉｏｎ  Ｐｒｏｃｅｓｓｉｎｇ  ｉｎ  Ｊａｐａｎ，９（１９６９）６０．）にある薄板とバネを用いたモデルのように、弾性エネルギを最小化することにより滑らかな曲面を生成する研究が盛んに行われている。その後、弾性エネルギなどの各種評価関数を用いて、全体としての最適化を図る研究が多数行われている（例えば、Ｊ．  Ｈｏｓｃｈｅｋ  ａｎｄ  Ｄ．  Ｌａｓｓｅｒ：  Ｆｕｎｄａｍｅｎｔａｌｓ  ｏｆ  Ｃｏｍｐｕｔｅｒ  Ａｉｄｅｄ  Ｇｅｏｍｅｔｒｉｃ  Ｄｅｓｉｇｎ，  Ａ  Ｋ  Ｐｅｔｅｒｓ（１９９３等参照）。
【０００５】
【発明が解決しようとする課題】
しかしながら、これらの従来の方法では誤差の中に異常データを含むとその影響が結果に表れ、測定データを対象にすると必ずしも良い結果が得られない。
【０００６】
そこで、自動車の意匠形状などの高品質曲面を対象として曲率変化に関する評価基準を導入し、全体の中からこれが最大となる点を逐次取りだしては局所的に平準化を繰り返し計算する方法が研究されている。この方法を最初に提案したのは、ファリン（Ｆａｒｉｎ）氏等である（例えば、Ｇ．  Ｆａｒｉｎ，Ｎ．  Ｓａｐｉｄｉｓ  ａｎｄ  Ｊ．  Ｗｏｒｓｅｙ：  Ｆｉｒｉｎｇ  Ｃｕｂｉｃ  Ｂ−Ｓｐｌｉｎｅ  Ｃｕｒｖｅｓ、Ｃｏｍｐｕｔｅｒ  Ａｉｄｅｄ  Ｇｅｏｍｅｔｒｉｃ  Ｄｅｓｉｇｎ，４（１９８７）９１．およびＮ．  Ｓａｐｉｄｉｓ  ａｎｄ  Ｇ．  Ｆａｒｉｎ：  Ａｕｔｏｍａｔｉｃ  Ｆａｉｒｒｉｎｇ  Ａｌｇｏｒｉｔｈｍ  ｆｏｒ  Ｂ−Ｓｐｌｉｎｅ  Ｃｕｒｖｅｓ、Ｃｏｍｕｔｅｒ  Ａｉｄｅｄ  Ｄｅｓｉｇｎ，２２，  ２（１９９０）１２１．等参照）。ファリン氏等は、入力点列をＢスプライン（Ｂ−ｓｐｌｉｎｅ）曲線として当てはめた後に、接続点での３階微分値、又は、曲率の微分値の差が最大の点を抽出し、その点を除去し再挿入することを繰り返すことにより、曲率プロット（線長を横軸に曲率を縦軸に示したグラフ）が滑らかとなる曲線を得ている。その後、Ｂスプライン（Ｂ−ｓｐｌｉｎｅ）曲線を用いず点群を平準化する方法として、離散的な曲率値を計算して、この差を局所的な評価基準とする方法が、エック（Ｅｃｋ）氏等によって示された（例えば、Ｍ．  Ｅｃｋ  ａｎｄ  Ｒ．  Ｊａｐｅｒｔ：  Ａｕｔｏｍａｔｉｃ  Ｆａｉｒｉｎｇ  ｏｆ  Ｐｏｉｎｔ  Ｓｅｔｓ、ｉｎ  Ｄｅｓｉｇｎｉｎｇ  Ｆａｉｒ  Ｃｕｒｖｅｓ  ａｎｄ  Ｓｕｒｆａｃｅｓ、ＳＩＡＭ（１９９４）４５．等参照）。また、最近では、ワグナー（Ｗａｇｎｅｒ）氏により、ロボットの経路を平滑化する手法として、曲率の代わりに点列の４階差分値を用いて、これを最小化する方法が提案されている（例えば、Ｍ．Ｇ．Ｗａｇｎｅｒ：  Ａｆｆｉｎｅ  Ｉｎｖａｒｉａｎｔ  Ｆａｉｒｉｎｇ  ｏｆ  Ｐｏｉｎｔ  Ｓｅｔ，  Ｐｒｏｃ．  ＣＩＩＳＴ’９８，  ＣＳＲＥＡ  Ｐｒｅｓｓ，（１９９８）３７０．等参照）。以上の局所的な平滑化を繰り返す処理法は、誤差に対して影響されず曲率変化レベルでの高品質の曲面データに対しては、その適用法が示されていない。
【０００７】
したがって、従来においては、例えば、リバースエンジニアリング等において、複雑な自由曲面から構成された意匠形状を再現したいような場合、その意匠形状を計測して得た点群データから先ず曲線を作成し、それら作成した曲線から面を作成するといった方法でモデリングを行っているのが先端技術であった。しかし、このような方法では、計測した意匠形状に近似した滑らかな曲面形状を再現して高品質なモデリングを行うには不十分である上、形状データの作成に時間と工数を要するという点でも問題が残っていた。
【０００８】
本発明の目的は、前述したような従来技術の問題点を解消しうるような自由曲面を表す曲面式を求める方法を提供することである。
【０００９】
【課題を解決するための手段】
本発明によれば、コンピュータグラフィックまたはリバースエンジニアリングにおいて自由曲面を表す曲面式を求める方法であって、該自由曲面を測定することによって得られた該自由曲面を示す測定データを、コンピュータを用いて点群データに分割し、該点群データにおける５点ごとの点列を４階差分商を用いて平滑化し、該平滑化された５点ごとの点列からなる点群データの各中央点の３次パラメトリック式を求め、該３次パラメトリック式に基づいて求めた５点ごとの点列からなる局所的点列データの中央点での曲面の１階、２階微分ベクトルを求め、これら微分ベクトルを４隅の微分ベクトルとする５次ベジエ曲面を生成することにより前記曲面式を得ることを特徴とする方法が提供される。ここでいう点群データは、測定点群、画像データ、色空間値、等の測定可能な値を持ったデータを総称するものである。
【００１０】
本発明の一つの実施の形態によれば、前記点群データは、４辺格子上の格子点データであり、前記平滑化は、５点ごとの点列に対する局所的平滑基準値を、前記４辺格子の一方をｕ方向、他の一方をｖ方向とする空間において、前記４辺格子を構成するｕ方向のパラメータ線とｖ方向のパラメータ線に対して求めた前記４階差分商の合計とし、前記点群データ全体における全体平滑基準値を全格子点について求めた前記局所的平滑基準値の和とし、前記ｕ方向およびｖ方向のパラメータ線を同時に平滑化するものである。
【００１３】
【発明の実施の形態】
次に、添付図面に基づいて、本発明の実施の形態および実施例について本発明をより詳細に説明する。
本発明者等は、本発明の方法の開発過程において、４階差分商を用いた曲線データの平滑化法の理論を明確化し、曲率レベルの滑らかさを得ることが出来る事を示し、添付図面の図１のフロー図に略示するように、この方法を、格子データを対象とする曲面の平滑化法に拡張し（ブロック１）、さらに、平滑化された格子データに対して曲面式を当てはめて（ブロック２）、曲面の曲率（Ｃ２）接続を行う（ブロック３）という手順をとると有効であることを確認した。以下、その実行例について説明することにより、その有効性を示す。
【００１４】
（１）曲線データの平滑化
平滑化の手法
先ず、４階差分商を用いて局所的に点列データを平滑化する方法について説明する。前述のファリン氏の方法では、曲率プロットの滑らかな曲線を得るため曲率の変化に着目し、Ｂスプライン（Ｂ−ｓｐｌｉｎｅ）曲線の接続点における両側セグメントの、曲率微分値の差を計算して、局所的な評価基準値（Ｌｏｃａｌ  Ｆａｉｒｎｅｓｓ）とし、これらを足し会わせた全体の評価値（Ｇｌｏｂａｌ  Ｆａｉｒｎｅｓｓ）が減少しなくなるまで、局所基準値の最大の点を繰り返し修正して曲線形状を平滑化している。局所的な点位置の修正法としては、その点を除去して挿入することを行っている。これに対してワーグナ氏は、曲率の微分値に相当する３階微分値の差を点列の４階差分値で表現し，これを局所平滑基準値として、この値が最大となる点を４階差分値が０となる位置に修正する方法を提案した。両者に共通するのは、３階微分ベクトルの差を滑らかさ（滑らかでないこと）の評価基準とすることと、局所的に形状を修正しようとする点である。この平滑化結果は、実用的にほぼ満足するものが得られており、本手法でも基本的にこの考え方を採用し、さらに、局所的な平滑化として、両者とも３次パラメトリック曲線を用いていることに着目して、点の修正法を４階差分商から理論的に明らかにする。
【００１５】
点列が接続点で３階微分値が連続であることは、局所的にその部分を３次式で近似できることと同値である。したがって、点列を考えた場合に５点ごとの区間が３次式で表せるとすると、ニュートン（Ｎｅｗｔｏｎ）氏の補間公式（Ｊ．  Ｈｏｓｃｈｅｋ  ａｎｄ  Ｄ．  Ｌａｓｓｅｒ：  Ｆｕｎｄａｍｅｎｔａｌｓ  ｏｆ  Ｃｏｍｐｕｔｅｒ  Ａｉｄｅｄ  Ｇｅｏｍｅｔｒｉｃ  Ｄｅｓｉｇｎ  ，Ａ  Ｋ  Ｐｅｔｅｒｓ（１９９３）参照）より、中央点での４階差分商は０となる。この差分商の絶対値を局所基準値として考え、この値が大きな点を、一番滑らかでない、部位とする。入力
0003610295.tif 000002
0003610295.tif 000003
0003610295.tif 000004
0003610295.tif 000005
ｋ区間の差分であり、
0003610295.tif 000006
0003610295.tif 000007
0003610295.tif 000008
である。式（１）の係数
0003610295.tif 000009
は、４階差分商の値を点列の間隔に会わせるスケールファクタであり、区間の大きさの３乗を掛けている。全体平滑化基準Ｆは、この局所平滑基準を全点列に対して加えて、
0003610295.tif 000010
となる。
【００１６】
曲線の平滑化法としては、従来の局所的平滑化の繰り返し処理と同様に行う。
0003610295.tif 000011
が０になるように、点位置を修正する。これを全体平滑基準が減少しなくなるまで繰り返す。
【００１７】
修正する点の位置は、４階差分商の値が０になることより、次の式より計算される。
0003610295.tif 000012
0003610295.tif 000013
0003610295.tif 000014
0003610295.tif 000015
0003610295.tif 000016
0003610295.tif 000017
ラメータが等間隔の場合には、式（８）は簡単になり、
0003610295.tif 000018
となる。これは、幾何学的にも図２に示すように簡潔に表される。
【００１８】
端点処理と終了条件
前述の平滑化手法により、ほぼ滑らかな曲線が得られるが、二つの問題が残っている。一つは曲線の、両側端点部の処理であり、他の一つはどの様に繰り返し計算を終了するかの終了条件の問題である。
第１の問題は、点列の修正は中央点を対象とするため、端の２点に誤差が含まれていると、これが修正されずに必ず端から３点目が移動するため、図３に示すように端２点で定まる接線方向が固定され、曲線全体の形状が影響を受けることである。この問題に対しては、端が３点目の局所平滑基準値が最大の場合には、端側の２点と中央の点の３点を夫々４階差分商を０とすることで修正した場合の全体平滑基準値を算出し、これが一番小さい点を選択し修正することにする。図４の例の場合には、２点目の全体平滑基準値Ｆが最小となるので、この２番目の点が選ばれて修正される。
【００１９】
第２の問題は、繰り返し計算を何処で止めるかである。単純に全体平滑基準値が減少しなくなる時点で計算を止めると、平滑化が十分行われない場合があり、これをやりすぎると図５に示すようにコーナ部のアールがやせてくる。そこで、特異点を除去するために、まず全体平滑基準値が減少しなくなるまで計算し、そこからは指定された誤差の最大値（トレーランス）以上移動した時点で、その前の全体平滑基準値が減少しない点にさかのぼって計算を終了することにする。
【００２０】
実行例
前述のアルゴリズムを用いて計算した結果を図６の（ａ）に示す。これは、スパイラル曲線に対して、ランダムな誤差を付加して作成した入力点列を平滑化した結果である。図６の（ｂ）および（ｃ）に示す曲率プロット、曲率プロファイル（曲線に沿って曲率中心の１／１０をプロットしたもの）ともに、入力点に比べて大幅に改善されていることが分かる。ここで、入力点の曲率値は、３点を用いて離散的に計算されている。
【００２１】
図７の（ａ）は凹凸形状を含む複数のベジエ（Ｂｅｚｉｅｒ）  曲線に対して、同様に誤差を付加したデータを平滑化した例である。この場合には、図７の（ｂ）および（ｃ）に示されるように、誤差による多数の曲率の反転が除去され、滑らかな曲率プロット、曲率プロファイルが得られている。なお、平滑化曲線の曲率値は両者とも、後述する曲線式への当てはめ結果を使用して計算している。
【００２２】
（２）格子データの平滑化
平滑化手法
前述した曲線での平滑化手法を曲面へ拡張する。曲面のデータは、４辺形格子上に整列された格子点を扱う。入力点Ｐｉｊとして、格子を構成する夫々の曲線のパラメータ方向に対して、差分商を定義する。
0003610295.tif 000019
0003610295.tif 000020
0003610295.tif 000021
0003610295.tif 000022
0003610295.tif 000023
【００２３】
以上を用いて曲線の場合と同様に、局所平滑基準と全体平滑基準を定める。局所平滑基準はｕ方向のパラメータ線とｖ方向のパラメータ線に対するものの合計とし、全体平滑基準は全格子点に対する和とする。
0003610295.tif 000024
0003610295.tif 000025
0003610295.tif 000026
【００２４】
格子の平滑化では、曲線の場合と同様に局所平滑基準の最大の点について、これを０とするように点位置を修正し、全体平滑基準が減少しなくなるまで繰り返す。点位置の修正では単純にｕ方向のパラメータ線とｖ方向のパラメータ線を交互に修正したり、両者の中間値をとるのでは、面として修正できず滑らかな結果を得ることができない。そこで，５×５の格子点について、両方向のパラメータ線を同時に平滑化することを考える。図８示すように、２５点のうち１６点を使って中央の９点を一度に修正する。中央の点以外はそれぞれのパラメータ線に沿って、曲線の４階差分商が０になる位置に修正し、中央点については、求められたそれぞれの４点を用いて、曲線の４階差分商が０となる位置に修正する。両者の位置は、同一の１６点からどちらの方向を先に計算するかであるから、各点についての係数は両者とも同一であり、修正点は同一点となる。
【００２５】
本手法の評価
単純に曲線としての平滑化を組み合わせた方法と本手法との比較を、例で示し、本手法の有効性を示す。図９は、双３次ベジエ曲線より１９×１２＝２２８点の格子点を計算したものであり、これに格子間隔の１／１０の誤差をランダムに与え、さらに意図的に異常に大きな誤差を数点与えたものが図１０である。このデータに対して本手法を適用した結果が図１１である。両方向のパラメータ線ともに滑らかとなっていることが分かる。
【００２６】
これに対し、図１２はｕ方向のパラメータ線とｖ方向のパラメータ線を交互に修正したものであり、図１３は両方向で修正位置の中間位置を修正位置としたものである。前者では、後に平滑化されたｖ方向については滑らかとなっているが、ｕ方向のパラメータ線は滑らかになっていない。一方後者の場合には、両方向の修正位置が逆向きとなると繰り返し計算を行っても点の位置は修正されず、滑らかな結果が得られない。これから、格子データを対象として一度に修正をする本手法の有効性が実証される。
【００２７】
（３）曲面式の当てはめ
平滑化された格子データに対して、曲面式を当てはめる。格子データはパラメータ線に沿って平滑化されているので、まず、これを曲面式に当てはめる。平滑化では５点ごとの点列に対し、４階差分商が０となるように中央点を移動しているので、これからはほぼ３次式で近似できる。したがって、５点のデータを使って、中央のセグメントを３次パラメトリック式で表す。
【００２８】
５点の点列（Ｐｉ，…，Ｐ  ｉ＋４）に対するニュートン（Ｎｅｗｔｏｎ）の補間公式は、式（２）の差分ベクトルＰｉｊを用いると、
0003610295.tif 000027
である。ここで、
0003610295.tif 000028
である。ただし、ｎ０（ｔ）＝１  である。
【００２９】
４階差分商の大きさはほぼ０であるので、３階差分商までを使って中央セグメン
0003610295.tif 000029
線ベクトルとして、
0003610295.tif 000030
を得る。後ろの４点を使うと、
0003610295.tif 000031
0003610295.tif 000032
同様に終点での接線ベクトルを求めて、ハルミッテ（Ｈｅｒｍｉｔｅ）  補間によりパラメトリック３次曲線を求める（図１４参照）。
【００３０】
以上の局所的のセグメントについての当てはめをすべてのｕ、ｖパラメータ線について求め、これらを４境界線としたパッチ式を、ク−ンズ（Ｃｏｏｎｓ）  の双１次ブレンドにより計算する。続いて、求められた２ｘ２のパッチについて、エイケン（Ａｉｔｋｅｎ）  の式でブレンドして中央点での曲面の１階、２階微分ベクトルを求める。最後に、これらの微分ベクトルを４隅の微分ベクトルとする双５次ベジエ（Ｂｅｚｉｅｒ）  曲面を生成する。この曲面は、パッチ間で２階微分ベクトルを共有しているので曲率（Ｃ２）連続となっている（図１５参照）。
【００３１】
（４）実行例
前述したようなアルゴリズムで生成した曲面の例を示す。図１６は、凹凸形状を有する曲面について、ランダムに発生した誤差と幾つかの大きな誤差のある４３×１９＝８１９点格子データを平滑化し、曲面式を当てはめたものである。繰り返し回数は５４０回である。図１６の（ａ）が入力データであり、図１６の（ｂ）が平滑化結果である。図１６の（ｂ）では、断面線を計算して示してある。この図から、滑らかな曲面が得られていることが分かる。
【００３２】
前述したように、本発明の方法は、リバースエンジニアリングにおいて高品質なモデリングを行うことができるように、誤差を持つ格子上の測定データを平滑化するため、４階差分商をもとに局所平滑基準を計算し、この総和である全体平滑基準が最小となるように、繰り返し点列を平滑化していく手法を、例えば、５×５の格子点ついて、ｕ、ｖパラメータ線を同時に考える事により、格子データの平滑化法に拡張し、その結果を曲面式として当てはめることにより、曲率変化がほぼ滑らかな曲率（Ｃ２）連続な曲面形状を得ることができるようにしたものである。
【００３３】
なお、前述したような本発明の方法は、コンピュータによる演算処理を通して実行するのにも適したものである。
【００３４】
【発明の効果】
（１）点群データから当てはめた曲面式で、例えば、リバースエンジニアリングにて高品質なモデリングを行える。
（２）製品形状を高品質でデジタル化できるため、点群データから曲線を作成し、曲線から曲面を作成する従来の方法と比べ、製品形状データの作成の時間と工数を短縮することができる。
（３）製品と設計データの誤差評価の自動化を計れ、精度向上に役立つ。
（４）試作モデルを設計データに容易に変換できる。
（５）製品の点群データから当てはめた曲面式から、製品と誤差の少ない高品質な金型を作製できる。すなわち、点群データから、そのまま金型を製作できる曲率変化の曲率（Ｃ２）連続な滑らかな高品質のモデルを作製できる。
【図面の簡単な説明】
【図１】本発明の方法の手順を概略的に示す図である。
【図２】曲線の平滑化における中心点の修正の態様を示す図である。
【図３】曲線の平滑化における終点の接線の影響を示す図である。
【図４】曲線の平滑化における終点の修正について説明するための図である。
【図５】曲線の平滑化における繰り返し計算の終了時について説明するための図である。
【図６】曲線の平滑化結果の例を示す図である。
【図７】凹凸形状を含む複数のベジエ曲線に対して誤差を付加したデータを平滑化した結果を示す図である。
【図８】格子データの中央９点の修正の態様を説明するための図である。
【図９】双３次ベジエ曲面のオリジナル例を示す図である。
【図１０】図９の双３次ベジエ曲面のオリジナル例に意図的に誤差を与えた例を示す図である。
【図１１】図１０のデータに対して本発明の方法を適用した結果を示す図である。
【図１２】図１０のデータに対してｕ方向のパラメータ線とｖ方向のパラメータ線を交互に修正した場合の結果例を示す図である。
【図１３】図１０のデータに対してｕ方向およびｖ方向の両方向で修正位置の中間位置を修正位置とした結果例を示す図である。
【図１４】曲線のローカルフィッティングを説明するための図である。
【図１５】クーンズのパッチの双１次ブレンドを説明するための図である。
【図１６】本発明の方法の実行例における結果を説明するための図である。
【符号の説明】
１  曲面の平滑化
２  曲面の当てはめ
３  曲面の接続
【図１】
0003610295.tif 000033
【図２】
0003610295.tif 000034
【図３】
0003610295.tif 000035
【図４】
0003610295.tif 000036
【図５】
0003610295.tif 000037
【図６】
0003610295.tif 000038
【図７】
0003610295.tif 000039
【図８】
0003610295.tif 000040
【図９】
0003610295.tif 000041
【図１０】
0003610295.tif 000042
【図１１】
0003610295.tif 000043
【図１２】
0003610295.tif 000044
【図１３】
0003610295.tif 000045
【図１４】
0003610295.tif 000046
【図１５】
0003610295.tif 000047
【図１６】
0003610295.tif 000048