(57)【要約】  一実施形態に係る図形情報処理装置は、有限の境界線を有する曲面のｕｖ座標系のパラメータと、該曲面の制御点に関するデータとを含む入力データを受け付ける受付部と、曲面上の一点での第１の曲率及び第２の曲率を用いて定義される拡張パラメータを用いた目的関数を記憶部から読み出し、ｕｖ座標系のパラメータを該目的関数に適用することで、該目的関数が最小値を取る制御点の位置を算出する算出部と、算出部により算出された制御点の位置に基づいて曲面を整形する整形部と、整形部により整形された曲面を示す図形情報を出力する出力部とを備える。
2013183396.tif 000002
【特許請求の範囲】
【請求項１】
有限の境界線を有する曲面のｕｖ座標系のパラメータと、該曲面の制御点に関するデータとを含む入力データを受け付ける受付部と、
前記曲面上の一点での第１の曲率及び第２の曲率を用いて定義される拡張パラメータを用いた目的関数を記憶部から読み出し、前記ｕｖ座標系のパラメータを該目的関数に適用することで、該目的関数が最小値を取る前記制御点の位置を算出する算出部と、
前記算出部により算出された制御点の位置に基づいて前記曲面を整形する整形部と、
前記整形部により整形された曲面を示す図形情報を出力する出力部と
を備える図形情報処理装置。
【請求項２】
前記目的関数が前記拡張パラメータと前記曲面の第１基本量とを用いて定義される、
請求項１に記載の図形情報処理装置。
【請求項３】
前記拡張パラメータをΩとした場合に前記目的関数が下記式で表され、
【数１】
2013183396.tif 000045
ここで、Ｉは、曲面の第１基本量で表される行列であり、Ωｕ２は、前記ｕｖ座標系におけるｕ方向の前記拡張パラメータの微分値を２乗した値であり、Ωｖ２は、前記ｕｖ座標系におけるｖ方向の前記拡張パラメータの微分値を２乗した値であり、Ｗは重みである、
請求項２に記載の図形情報処理装置。
【請求項４】
前記第１の曲率が最大曲率であり、前記第２の曲率が最小曲率である、
請求項３に記載の図形情報処理装置。
【請求項５】
前記拡張パラメータが、前記最大曲率をｋｍａｘとし、前記最小曲率をｋｍｉｎとし、α及びβを任意の数として、Ω＝ｋｍａｘ−αｋｍｉｎ−βで示される、
請求項４に記載の図形情報処理装置。
【請求項６】
前記拡張パラメータが、前記最大曲率をｋｍａｘとし、前記最小曲率をｋｍｉｎとし、α及びβを任意の数として、Ω＝（ｋｍａｘ−α＋ｋｍｉｎ−β）／２で示される、
請求項４に記載の図形情報処理装置。
【請求項７】
前記第１の曲率が前記ｕｖ座標系におけるｕ方向の曲率であり、前記第２の曲率が該ｕｖ座標系におけるｖ方向の曲率である、
請求項３に記載の図形情報処理装置。
【請求項８】
前記拡張パラメータが、前記ｕ方向の曲率の微分値をｋｕとし、前記ｖ方向の曲率の微分値をｋｖとし、α及びβを任意の数として、Ω＝ｋｕ−αｋｖ−βで示される、
請求項７に記載の図形情報処理装置。
【請求項９】
前記拡張パラメータが、前記ｕ方向の曲率の微分値をｋｕとし、前記ｖ方向の曲率の微分値をｋｖとし、α及びβを任意の数として、Ω＝（ｋｕ−α＋ｋｖ−β）／２で示される、
請求項７に記載の図形情報処理装置。
【請求項１０】
前記拡張パラメータが、前記ｕ方向の曲率の微分値をｋｕとし、前記ｖ方向の曲率の微分値をｋｖとし、α及びβを任意の数として、Ω＝ｋｕ−α＋ｋｖ−βで示される、
請求項７に記載の図形情報処理装置。
【請求項１１】
前記制御点に関するデータが曲面式であり、
前記算出部が前記曲面式により決まる前記制御点の位置を算出し、
前記整形部が、前記算出部により算出された制御点の位置を前記曲面式に適用することで前記曲面を整形する、
請求項１〜１０のいずれか一項に記載の図形情報処理装置。
【請求項１２】
前記受付部が、前記制御点の位置の変更に対する拘束条件を更に受け付け、
前記算出部が、前記拘束条件を満たす範囲内で、前記目的関数が最小値を取る前記制御点の位置を算出する、
請求項１〜１１のいずれか一項に記載の図形情報処理装置。
【請求項１３】
図形情報処理装置により実行される図形情報処理方法であって、
有限の境界線を有する曲面のｕｖ座標系のパラメータと、該曲面の制御点に関するデータとを含む入力データを受け付ける受付ステップと、
前記曲面上の一点での第１の曲率及び第２の曲率を用いて定義される拡張パラメータを用いた目的関数を記憶部から読み出し、前記ｕｖ座標系のパラメータを該目的関数に適用することで、該目的関数が最小値を取る前記制御点の位置を算出する算出ステップと、
前記算出ステップにおいて算出された制御点の位置に基づいて前記曲面を整形する整形ステップと、
前記整形ステップにおいて整形された曲面を示す図形情報を出力する出力ステップと
を含む図形情報処理方法。
【請求項１４】
有限の境界線を有する曲面のｕｖ座標系のパラメータと、該曲面の制御点に関するデータとを含む入力データを受け付ける受付部と、
前記曲面上の一点での第１の曲率及び第２の曲率を用いて定義される拡張パラメータを用いた目的関数を記憶部から読み出し、前記ｕｖ座標系のパラメータを該目的関数に適用することで、該目的関数が最小値を取る前記制御点の位置を算出する算出部と、
前記算出部により算出された制御点の位置に基づいて前記曲面を整形する整形部と、
前記整形部により整形された曲面を示す図形情報を出力する出力部と
をコンピュータに実行させる図形情報処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の一形態は、曲面を生成若しくは修正する図形情報処理装置、図形情報処理方法、及び図形情報処理プログラムに関する。
【背景技術】
【０００２】
工業製品の設計、作画、グラフィック・デザイン、フォント作成、あるいは手書き文字の入力などにおいて、自由曲線および自由曲面は基本的なデザイン要素である。そのうち曲線については、デザインに好適な曲線（人が美しいと感じる曲線、あるいは対数型美的曲線とも言う）を生成するための図形情報処理装置が下記特許文献１に記載されている。この図形情報処理装置は、まず、線図形を複数のセグメントに分割し、続いて、各セグメントの曲率半径または曲率のべき乗が該セグメントの曲線長の一次関数で表されるように各セグメントにおける曲線を整形する。そして、この装置は、セグメントの一端点での曲率と該一端点に接続している他のセグメントの一端点での曲率との差が所定の閾値以下になるように、複数のセグメントを整形する。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】国際公開第２００９／１４８１５７号パンフレット
【発明の概要】
【発明が解決しようとする課題】
【０００４】
一方、曲面については、上記特許文献に記載の手法を単純に拡張することで、デザインに好適な曲面を生成することが考えられる。しかし、この単純な拡張では、理想的なパラメータを設定したり、二つの主曲率に対して適切な目的関数を設定したりする必要があるので、複雑な曲面のデザインにその手法を適用することが非常に困難か又は不可能である。そこで、様々な曲面をデザインに好適な態様に整形できるような図形処理の手法が要請されている。
【課題を解決するための手段】
【０００５】
本発明の一形態に係る図形情報処理装置は、有限の境界線を有する曲面のｕｖ座標系のパラメータと、該曲面の制御点に関するデータとを含む入力データを受け付ける受付部と、曲面上の一点での第１の曲率及び第２の曲率を用いて定義される拡張パラメータを用いた目的関数を記憶部から読み出し、ｕｖ座標系のパラメータを該目的関数に適用することで、該目的関数が最小値を取る制御点の位置を算出する算出部と、算出部により算出された制御点の位置に基づいて曲面を整形する整形部と、整形部により整形された曲面を示す図形情報を出力する出力部とを備える。【０００６】
本発明の一形態に係る図形情報処理方法は、図形情報処理装置により実行される図形情報処理方法であって、有限の境界線を有する曲面のｕｖ座標系のパラメータと、該曲面の制御点に関するデータとを含む入力データを受け付ける受付ステップと、曲面上の一点での第１の曲率及び第２の曲率を用いて定義される拡張パラメータを用いた目的関数を記憶部から読み出し、ｕｖ座標系のパラメータを該目的関数に適用することで、該目的関数が最小値を取る制御点の位置を算出する算出ステップと、算出ステップにおいて算出された制御点の位置に基づいて曲面を整形する整形ステップと、整形ステップにおいて整形された曲面を示す図形情報を出力する出力ステップとを含む。【０００７】
本発明の一形態に係る図形情報処理プログラムは、有限の境界線を有する曲面のｕｖ座標系のパラメータと、該曲面の制御点に関するデータとを含む入力データを受け付ける受付部と、曲面上の一点での第１の曲率及び第２の曲率を用いて定義される拡張パラメータを用いた目的関数を記憶部から読み出し、ｕｖ座標系のパラメータを該目的関数に適用することで、該目的関数が最小値を取る制御点の位置を算出する算出部と、算出部により算出された制御点の位置に基づいて曲面を整形する整形部と、整形部により整形された曲面を示す図形情報を出力する出力部とをコンピュータに実行させる。【０００８】
このような形態によれば、曲面上の一点での二種類の曲率を用いて定義される拡張パラメータに着目して、当該拡張パラメータを用いた目的関数が用いられる。そして、この目的関数が最小値を取るような曲面の制御点の位置が算出され、その制御点の位置に基づいて曲面が整形される。このように、拡張パラメータに基づいて曲面の制御点を調整することで、様々な曲面をデザインに好適な態様に整形することが可能になる。【０００９】
別の形態に係る図形情報処理装置では、目的関数が拡張パラメータと曲面の第１基本量とを用いて定義されてもよい。【００１０】
さらに別の形態に係る図形情報処理装置では、拡張パラメータをΩとした場合に目的関数が下記式（４１）で表され、ここで、当該式（４１）において、Ｉは、曲面の第１基本量で表される行列であり、Ωｕ２は、ｕｖ座標系におけるｕ方向の拡張パラメータの微分値を２乗した値であり、Ωｖ２は、ｕｖ座標系におけるｖ方向の拡張パラメータの微分値を２乗した値であり、Ｗは重みであってもよい。【００１１】
さらに別の形態に係る図形情報処理装置では、第１の曲率が最大曲率であり、第２の曲率が最小曲率であってもよい。【００１２】
さらに別の形態に係る図形情報処理装置では、拡張パラメータが、最大曲率をｋｍａｘとし、最小曲率をｋｍｉｎとし、α及びβを任意の数として、Ω＝ｋｍａｘ−αｋｍｉｎ−βで示されてもよい。【００１３】
さらに別の形態に係る図形情報処理装置では、拡張パラメータが、最大曲率をｋｍａｘとし、最小曲率をｋｍｉｎとし、α及びβを任意の数として、Ω＝（ｋｍａｘ−α＋ｋｍｉｎ−β）／２で示されてもよい。【００１４】
さらに別の形態に係る図形情報処理装置では、第１の曲率がｕｖ座標系におけるｕ方向の曲率であり、第２の曲率が該ｕｖ座標系におけるｖ方向の曲率であってもよい。【００１５】
さらに別の形態に係る図形情報処理装置では、拡張パラメータが、ｕ方向の曲率の微分値をｋｕとし、ｖ方向の曲率の微分値をｋｖとし、α及びβを任意の数として、Ω＝ｋｕ−αｋｖ−βで示されてもよい。【００１６】
さらに別の形態に係る図形情報処理装置では、拡張パラメータが、ｕ方向の曲率の微分値をｋｕとし、ｖ方向の曲率の微分値をｋｖとし、α及びβを任意の数として、Ω＝（ｋｕ−α＋ｋｖ−β）／２で示されてもよい。【００１７】
さらに別の形態に係る図形情報処理装置では、拡張パラメータが、ｕ方向の曲率の微分値をｋｕとし、ｖ方向の曲率の微分値をｋｖとし、α及びβを任意の数として、Ω＝ｋｕ−α＋ｋｖ−βで示されてもよい。【００１８】
さらに別の形態に係る図形情報処理装置では、制御点に関するデータが曲面式であり、算出部が曲面式により決まる制御点の位置を算出し、整形部が、算出部により算出された制御点の位置を曲面式に適用することで曲面を整形してもよい。【００１９】
さらに別の形態に係る図形情報処理装置では、受付部が、制御点の位置の変更に対する拘束条件を更に受け付け、算出部が、拘束条件を満たす範囲内で、目的関数が最小値を取る制御点の位置を算出してもよい。
【発明の効果】
【００２０】
本発明の一側面によれば、様々な曲面をデザインに好適な態様に整形することができる。
【図面の簡単な説明】
【００２１】
【図１】変分原理を説明するための図である。
【図２】実施形態に係るＣＡＤシステムのハードウェア構成を示す図である。
【図３】実施形態に係るＣＡＤシステムの機能構成を示すブロック図である。
【図４】ｕｖ座標系の定義の一例を示す図である。
【図５】制御点の設定の一例を示す図である。
【図６】図３に示すＣＡＤシステムの動作を示すフローチャートである。
【図７】実施形態に係る図形情報処理プログラムの構成を示す図である。
【発明を実施するための形態】
【００２２】
以下、添付図面を参照しながら本発明の実施形態を詳細に説明する。なお、図面の説明において同一又は同等の要素には同一の符号を付し、重複する説明を省略する。
【００２３】
本実施形態では、本発明に係る図形情報処理装置をＣＡＤ（Computer Aided Design）システム１０に適用する。ＣＡＤシステム１０は、平面上に表現された三次元モデルの図形を規定するデータを受け付け、そのデータで表現される図形を整形して出力するコンピュータシステムである。より具体的には、ＣＡＤシステム１０はモデル図形に含まれる曲面を整形することで、デザインに好適な曲面（対数型美的曲面）を生成する。
【００２４】
ＣＡＤシステム１０の構成を具体的に説明する前に、そのＣＡＤシステム１０により生成される美的曲面の理論を説明する。この理論は、変分原理に基づく対数型美的曲面の定式化である。
【００２５】
理解を容易にするために、まず曲線について説明する。平面曲線である対数型美的曲線は、曲率対数分布グラフの傾きαが０でない場合に次式（１）を満足する。
【数１】
2013183396.tif 000003
ここで、ρは曲率半径、ｓは曲線長、ｃ及びｄは定数である。α＝０であれば次式（２）が得られる。
【数２】
2013183396.tif 000004
ここで、ｃ０及びｃ１は定数である。また、その曲率対数分布グラフは直線となる。したがって、次式（３）で示される関係が得られる。
【数３】
2013183396.tif 000005
【００２６】
この式（３）は、座標系（ｌｏｇρ，ｌｏｇ（ρｄｓ／ｄρ））の２点を結ぶ最短経路として対数型美的曲線が定義されることを意味している。変分原理によれば、２点間の最短経路は直線であることが示される。
【００２７】
変分計算では、ある積分量を最小又は最大にする問題（停留積分）を扱う。例えば次式（４）を設定する。
【数４】
2013183396.tif 000006
ここで、ｙはｘの関数であり、ｙｘはｙのｘに関する導関数である。ｙは未知である。積分量Ｊが停留値を取るための条件式は下記偏微分方程式（５）となる。
【数５】
2013183396.tif 000007
この式（５）はオイラー方程式と呼ばれ、ｄｙ／ｄｘ≠０であれば次式（６）と等価である。
【数６】
2013183396.tif 000008
【００２８】
もし、ｆ＝ｆ（ｙ，ｙｘ）のように変数ｘが明示的に現われない問題では、上記式（６）は次式（７）で示され、したがって、次式（８）の関係が得られる。
【数７】
2013183396.tif 000009
【数８】
2013183396.tif 000010
【００２９】
オイラー方程式の最も簡単な応用は、ｘｙ平面上の２点間の最短距離を決める問題である。距離の微小要素は次式（９）で定義されるので、距離Ｊは次式（１０）で示される。
【数９】
2013183396.tif 000011
【数１０】
2013183396.tif 000012
【００３０】
この式（１０）と式（４）とを比較すると次式（１１）が得られる。
【数１１】
2013183396.tif 000013
【００３１】
この式（１１）を式（８）に代入すると次式（１２）が得られる。ここで、Ｃは定数である。
【数１２】
2013183396.tif 000014
【００３２】
したがって、定数ａが存在して次式（１３）が成立し、さらに次式（１４）が得られる。
【数１３】
2013183396.tif 000015
【数１４】
2013183396.tif 000016
ここで、ｂも定数である。これら二つの定数ａ，ｂは、２点（ｘ１，ｙ１）及び（ｘ２，ｙ２）を通るように決定される。
【００３３】
式（１）においてρα＝σと定義すると次式（１５）が得られる。
【数１５】
2013183396.tif 000017
この式（１５）は、図１に示すように、横軸をｓとし縦軸をσ＝ραとする２次元空間で２点（ｓ１，σ１）及び（ｓ２，σ２）を直線で（すなわち最短距離で）結ぶと対数型美的曲線が得られることを示している。この場合には、次式（１６）が最小化される。
【数１６】
2013183396.tif 000018
ここで、σｓ２＝（ｄσ／ｄｓ）２である。
【００３４】
α＝１ならば、その曲線は対数（等角）らせんである。この場合には次式（１７）が成立する。この式（１７）は、横軸をｓとし縦軸をρとする２次元空間で２点（ｓ１，ρ１）及び（ｓ２，ρ２）を最短で結ぶ曲線が対数らせんであることを意味する。
【数１７】
2013183396.tif 000019
【００３５】
一方、α＝−１ならば、その曲線はクロソイド曲線である。この場合には次式（１８）が成立する。ここで、κ＝１／ρは曲率である。この式（１８）は、横軸をｓとし縦軸をκとする２次元空間で２点（ｓ１，κ１）及び（ｓ２，κ２）を最短で結ぶ曲線がクロソイド曲線であることを意味する。
【数１８】
2013183396.tif 000020
【００３６】
曲線が、曲線長ｓではなく一般的なパラメータｔで与えられるとすると、次式（１９）が成り立つ。
【数１９】
2013183396.tif 000021
【００３７】
したがって、式（１６）は次式（２０）となる。
【数２０】
2013183396.tif 000022
【００３８】
したがって、次式（２１）が成立する。
【数２１】
2013183396.tif 000023
【００３９】
式（１６）に対する式（５）を求めると、次式（２２），（２３）より、当該式（５）は次式（２４）となる。
【数２２】
2013183396.tif 000024
【数２３】
2013183396.tif 000025
【数２４】
2013183396.tif 000026
【００４０】
したがって、次式（２５）が得られる。
【数２５】
2013183396.tif 000027
【００４１】
対数型美的曲線では、ραはｓの１次式で与えられるので次式（２６）が成立し、したがって、次式（２７）が得られる。
【数２６】
2013183396.tif 000028
【数２７】
2013183396.tif 000029
この式（２７）は式（２５）と一致する。
【００４２】
ρ＝１／κ、ρｓ＝−κｓ／κ２、ρｓｓ＝−（κｓｓκ＋２κｓ）／κ３であるから、式（２５）は曲率κを用いて次式（２８）で表すこともできる。
【数２８】
2013183396.tif 000030
【００４３】
式（２５）の解を求める。この式はｓを独立変数とする非線形２階常微分方程式であるが、ｓを陽に含まないので、この式を、ρを独立変数としｓを従属変数とする式に置き換える。ρｓ＝１／ｓρ，ρｓｓ＝−ｓρρ／ｓρ３であるから、式（２５）は次式（２９）となる。
【数２９】
2013183396.tif 000031
ここでｔ＝ｓρとすると、次式（３０）が得られる。
【数３０】
2013183396.tif 000032
【００４４】
この式（３０）はアーベルの常微分方程式（y'(x)=f(x)+g(x)y(x)+h(x)y(x)2+k(x)y(x)3）であることに留意する必要がある。したがって、次式（３１）が得られる。
【数３１】
2013183396.tif 000033
【００４５】
この式（３１）においてα≠０であれば次式（３２）が得られ、α＝０であれば次式（３３）が得られる。
【数３２】
2013183396.tif 000034
【数３３】
2013183396.tif 000035
これら二つの式（３２），（３３）は共に、対数型美的曲線に一致する。
【００４６】
以上説明した変分原理の考え方を曲面に拡張する。曲線長ｓに曲面の面積Ａを対応させ、曲率κにガウス曲率Κを対応させる。曲面上の一点での最大曲率及び最小曲率をそれぞれｋｍａｘ，ｋｍｉｎとすると、ρ＝κ−αであり、Κ＝ｋｍａｘｋｍｉｎであるから、拡張パラメータΩをΩ＝ｋｍａｘ−αｋｍｉｎ−βと定義する。
【００４７】
次式（３４）をｓ＝ｓ（ｔ）で再パラメータ化すると、次式（３５）が得られる。
【数３４】
2013183396.tif 000036
【数３５】
2013183396.tif 000037
【００４８】
したがって、次式（３６）が得られる。
【数３６】
2013183396.tif 000038
【００４９】
ここでλｃ（ｔ）を
【数３７】
2013183396.tif 000039
と定義すると、次式（３７）が得られる。
【数３８】
2013183396.tif 000040
【００５０】
この式（３７）を曲面に拡張して次式（３８）を定義する。
【数３９】
2013183396.tif 000041
ここで、Ｉは曲面の第１基本量Ｅ，Ｆ，Ｇで表される行列であり、次式（３９）で定義され、局所的な曲面の面積の２乗に対応する。
【数４０】
2013183396.tif 000042
第１基本量Ｅ，Ｆ，Ｇは、曲面式をＳ＝Ｓ（ｕ，ｖ）とすると以下のように定義される。なお、ＳｕはＳのｕに関する導関数であり、ＳｖはＳのｖに関する導関数であり、“・”は内積を表す。
Ｅ＝Ｓｕ・Ｓｕ
Ｆ＝Ｓｕ・Ｓｖ
Ｇ＝Ｓｖ・Ｓｖ
【００５１】
したがって、曲率を線形化する対数型美的曲線に対応する曲面の目的関数は次式（４０）で示される。すなわち、この式（４０）を最小化する曲面が対数型美的曲面である。ここで、Ωｕ２は、ｕｖ座標系におけるｕ方向の拡張パラメータΩの微分値を２乗した値であり、Ωｖ２は、ｕｖ座標系におけるｖ方向の拡張パラメータΩの微分値を２乗した値である。
【数４１】
2013183396.tif 000043
【００５２】
この目的関数は、スカラー定数である重みＷを用いて次式（４１）で定義されてもよい。この重みＷはデザインパラメータとして設計者により任意に設定されてよい。Ｗ＝１であれば式（４１）は式（４０）と一致する。
【数４２】
2013183396.tif 000044
【００５３】
次に、ＣＡＤシステム１０の構成を説明する。ＣＡＤシステム１０のハードウェア構成を図２に示す。この図に示すように、ＣＡＤシステム１０は、オペレーティングシステムやアプリケーション・プログラムなどを実行するＣＰＵ１０１と、ＲＯＭ及びＲＡＭで構成される主記憶部１０２と、ハードディスクなどで構成される補助記憶部１０３と、ネットワークカードなどで構成される通信制御部１０４と、キーボードやマウスなどの入力装置１０５と、ディスプレイなどの出力装置１０６とを備えている。
【００５４】
後述するＣＡＤシステム１０の各機能的構成要素は、ＣＰＵ１０１又は主記憶部１０２の上に所定のソフトウェアを読み込ませ、ＣＰＵ１０１の制御の下で通信制御部１０４や入力装置１０５、出力装置１０６などを動作させ、主記憶部１０２又は補助記憶部１０３におけるデータの読み出し及び書き込みを行うことで実現される。処理に必要なデータやデータベースは主記憶部１０２又は補助記憶部１０３内に格納される。なお、図２ではＣＡＤシステム１０が１台のコンピュータで構成されているように示しているが、ＣＡＤシステム１０の機能を複数台のコンピュータに分散させてもよい。
【００５５】
図３に示すように、ＣＡＤシステム１０は、機能的構成要素として受付部１１、算出部１２、整形部１３、及び出力部１４を備えている。
【００５６】
受付部１１は、対数型美的曲面の生成に必要な入力データを受け付ける機能要素である。入力データは、曲面の有限の境界線を示す点列データと、曲面のｕｖ座標系における一組のパラメータｕ，ｖと、内挿に用いる曲面式と、曲面の形状を調整するための上記パラメータα，βと、式（４１）における重みＷと、制御点の位置の変更に対する拘束条件とを含む。
【００５７】
曲面の境界線の形状は限定されない。例えば、受付部１１は四つの曲線から成る境界線を受け付けてもよく、この場合には各曲線が上記特許文献１で示される美的曲線であってもよい。あるいは、受付部１１は四つの直線から成る境界線、曲線及び直線の双方を含む境界線、円形の境界線、若しくは自由曲線から成る境界線を受け付けてもよい。
【００５８】
ｕｖ座標系は、境界線が４曲線から成り立っているのであれば、例えば図４に示すように定義される。パラメータｕ，ｖはそれぞれ正規化されて０から１の間の値を取る。
【００５９】
曲面式は、生成しようとする曲面を規定する制御点の個数（自由度）を決めるために入力される関数であり、制御点に関するデータに相当する。この曲面式も限定されない。例えば、ユーザは、比較的単純な曲面を生成したいのであれば双３次ベジエ曲面を示す曲面式を指定すればよいし、より複雑な曲面を生成したいのであればｍ×ｎ次ＮＵＲＢＳ（Non-Uniform Rational B-Spline（非一様有利Ｂスプライン））曲面を指定すればよい。受付部１１は指定された曲面式を示すデータを受け付ける。
【００６０】
パラメータα及びβは曲面の形状を調整するためのパラメータであり、それぞれ−∞から＋∞の範囲で任意に設定可能である。
【００６１】
拘束条件の種類は限定されない。例えば、受付部１１は、整形後の曲面が境界線を通過するという条件、境界線における接平面の指定、隣接する曲面との間で曲率が連続するという条件（曲率の連続性）などを拘束条件として受け付ける。
【００６２】
受付部１１はこれらのデータを受け付けて算出部１２に出力する。なお、受付部１１はユーザがスタイラスペン、キーボード、又はマウスにより直接入力したデータを受けてもよいし、ユーザ指示に基づいて所定の記憶部（例えば、補助記憶部１０３や外部の記憶装置など）から読み出したデータを受け付けてもよい。
【００６３】
算出部１２は、式（４１）で示される目的関数が最小値を取る制御点の位置を算出する機能要素である。算出部１２は受付部１１から取得した入力データを用いてこの計算を実行する。この算出部１２は、上述した理論が組み込まれたプログラムがＣＰＵ１０１上で実行されることで実現される。以下では、上記理論にそのまま従って目的関数を用いる第１の手法と、処理時間を削減するために拡張パラメータΩの変形を用いる第２〜第５の手法とについて説明する。
【００６４】
［第１の手法］
まず、算出部１２は曲面式から制御点の個数を特定する。図５の例では、算出部１２は１６個の制御点（Ｐ００，Ｐ０１，…，Ｐ３３）を特定している。続いて、算出部１２は曲面式で定義される曲面上の一点での最大曲率ｋｍａｘ及び最小曲率ｋｍｉｎを算出し、これらの曲率とパラメータα，βとから、拡張パラメータΩ＝ｋｍａｘ−αｋｍｉｎ−βを算出する。したがって、この拡張パラメータは曲面上の一点での最大曲率及び最小曲率で定義される。さらに、算出部１２は曲面の第１基本量Ｅ，Ｆ，Ｇを求めることで、式（３９）で示される行列Ｉを得る。
【００６５】
続いて、算出部１２は式（４１）で示される目的関数を所定の記憶部（例えば、補助記憶部１０３や外部の記憶装置など）から読み出す。なお、目的関数はプログラム内に記述されていてもよい。そして、算出部１２は、パラメータｕ，ｖをその目的関数に適用（代入）することで、拘束条件下において当該目的関数が最小値を取る１以上の制御点の位置を算出する。例えば、隣接する曲面との間の曲率の連続性を示す拘束条件が入力された場合には、算出部１２は、当該隣接曲面を示す図形情報を上記記憶部から読み出す。そして、算出部１２は隣接曲面との境界線における２曲面の曲率の差が一致するか又は所定の閾値以下である範囲内で、目的関数を最小化する１以上の制御点の位置を求める。
【００６６】
算出部１２はこのような手順で求めた各制御点の位置を算出結果として整形部１３に出力する。また、算出部１２は入力データも整形部１３に出力する。
【００６７】
［第２の手法］
ガウス曲率Κの代わりに平均曲率Ｈ＝（ｋｍａｘ＋ｋｍｉｎ）／２を用いることを考えて、式（４１）における拡張パラメータをΩ＝（ｋｍａｘ−α＋ｋｍｉｎ−β）／２で置き換えてもよい。この場合には、算出部１２はこの変形例に係る目的関数（式（４１））が最小値を取る１以上の制御点の位置を算出する。
【００６８】
［第３の手法］
式（４１）における拡張パラメータをΩ＝ｋｕ−αｋｖ−βで置き換えてもよい。ここで、ｋｕはｕｖ座標系におけるｕ方向の曲率の微分値であり、ｋｖはｕｖ座標系におけるｖ方向の曲率の微分値である。算出部１２はこの変形例に係る目的関数（式（４１））が最小値を取る１以上の制御点の位置を算出する。
【００６９】
［第４の手法］
上記第２及び第３の手法を応用して、式（４１）における拡張パラメータをΩ＝（ｋｕ−α＋ｋｖ−β）／２で置き換えてもよい。この場合には、算出部１２はこの変形例に係る目的関数（式（４１））が最小値を取る１以上の制御点の位置を算出する。
【００７０】
［第５の手法］
二次元モデルの場合と同じ目的関数を得るために、式（４１）における拡張パラメータをΩ＝ｋｕ−α＋ｋｖ−βで置き換えてもよい。これは、上記第４の手法において分母を１にしたことを意味する。例えば、平面曲線をそれに直交する方向に掃引した場合には、ｋｍｉｎ＝（ｋｍｉｎ）ｕ＝（ｋｍｉｎ）ｖ＝０、且つ（ｋｍａｘ）ｖ＝０であるので、分母を１とする必要がある。
【００７１】
このように、上記第１〜第５の手法で互いに異なる点は、拡張パラメータの定義のみである。
【００７２】
整形部１３は、算出部１２により算出された制御点の位置に基づいて曲面を整形する機能要素である。整形部１３は、算出部１２により算出された制御点の位置を、入力データの一部である曲面式に適用（代入）することで対数型美的曲面を生成する。これにより、入力データで指定された曲面が対数型美的曲面に整形される。整形部１３は、整形後の曲面を示す点列データを図形情報として出力部１４に出力する。
【００７３】
出力部１４は、整形部１３から入力された図形情報をモニタやプリンタなどに出力する機能要素である。これにより、ＣＡＤシステム１０のユーザはデザインに好適な曲面（対数型美的曲面）を視認できる。
【００７４】
次に、図６を用いて、ＣＡＤシステム１０の処理を説明するとともに本実施形態に係る図形情報処理方法について説明する。
【００７５】
まず、受付部１１が必要な入力データを受け付ける（ステップＳ１１、受付ステップ）。具体的には、受付部１１は曲面の境界線のデータと、パラメータｕ，ｖと、曲面式と、パラメータα，βと、制御点の拘束条件とを受け付ける。続いて、算出部１２が、式（４１）で示される目的関数が最小値を取る制御点の位置を算出する（ステップＳ１２、算出ステップ）。具体的には、算出部１２はパラメータｕ，ｖを目的関数に適用（代入）することで、拘束条件下において当該目的関数が最小値を取る制御点の位置を算出する。続いて、整形部１３が算出された制御点の位置に基づいて曲面を整形する（ステップＳ１３、整形ステップ）。そして、出力部１４が整形された曲面（対数型美的曲面）を出力する（ステップＳ１４、出力ステップ）。
【００７６】
次に、図７を用いて、コンピュータをＣＡＤシステム１０として機能させるための図形情報処理プログラムＰ１を説明する。
【００７７】
図形情報処理プログラムＰ１は、メインモジュールＰ１０、受付モジュールＰ１１、算出モジュールＰ１２、整形モジュールＰ１３、及び出力モジュールＰ１４を備えている。
【００７８】
メインモジュールＰ１０は、上述した一連の図形情報処理を統括的に制御する部分である。受付モジュールＰ１１、算出モジュールＰ１２、整形モジュールＰ１３、及び出力モジュールＰ１４を実行することにより実現される機能はそれぞれ、上記の受付部１１、算出部１２、整形部１３、及び出力部１４の機能と同様である。
【００７９】
図形情報処理プログラムＰ１は、例えば、ＣＤ−ＲＯＭやＤＶＤ−ＲＯＭ、半導体メモリ等の有形の記録媒体に固定的に記録された上で提供される。また、図形情報処理プログラムＰ１は、搬送波に重畳されたデータ信号として通信ネットワークを介して提供されてもよい。
【００８０】
以上説明したように、本実施形態によれば、曲面上の一点での二種類の曲率を用いて定義される拡張パラメータΩに着目して、当該拡張パラメータを用いた目的関数（式（４１））が用いられる。そして、この目的関数が最小値を取るような曲面の制御点の位置が算出され、その制御点の位置に基づいて曲面が整形される。このように、拡張パラメータに基づいて曲面の制御点を調整することで、様々な曲面をデザインに好適な態様に整形することが可能になる。ＣＡＤシステム１は様々な物品のデザインやグラフィック・デザインにおいて利用でき、設計者はこのシステム１を用いることで、見る者が美しいと感じる曲面を簡単に得ることができる。
【００８１】
以上、本発明をその実施形態に基づいて詳細に説明した。しかし、本発明は上記実施形態に限定されるものではない。本発明は、その要旨を逸脱しない範囲で様々な変形が可能である。
【００８２】
上記実施形態では、受付部１１は入力データの一部として拘束条件を受け付けたが、この拘束条件は省略可能である。
【符号の説明】
【００８３】
１０…ＣＡＤシステム（図形情報処理装置）、１１…受付部、１２…算出部、１３…整形部、１４…出力部、Ｐ１…図形情報処理プログラム、Ｐ１０…メインモジュール、Ｐ１１…受付モジュール、Ｐ１２…算出モジュール、Ｐ１３…整形モジュール、Ｐ１４…出力モジュール。
【図１】
2013183396.tif 000046
【図２】
2013183396.tif 000047
【図３】
2013183396.tif 000048
【図４】
2013183396.tif 000049
【図５】
2013183396.tif 000050
【図６】
2013183396.tif 000051
【図７】
2013183396.tif 000052
【国際調査報告】
2013183396.tif 000053
2013183396.tif 000054