(57)【特許請求の範囲】
【請求項１】
二次元平面上の作業領域内に多数の母点を定義する母点定義部と、
前記作業領域内に定義された個々の母点について、隣接する母点との間にそれぞれ連結線を定義する連結線定義部と、
前記作業領域に対応する参照領域内にスカラー場を定義するスカラー場定義部と、
前記個々の連結線上に、前記スカラー場上の当該連結線の位置に対応する位置に定義されたスカラー値に応じた物理的特性をもった仮想ばねを定義する仮想ばね定義部と、
個々の母点を質点とみなし、各連結線の両端に位置する一対の母点が当該連結線上に定義された仮想ばねで連結されているとみなしたときに、各仮想ばねによって作用する力に基く各質点の運動シミュレーションを実行するシミュレーション実行部と、
前記運動シミュレーションによって運動中の母点の所定時点における位置を変位母点の位置と決定する変位母点決定部と、
各変位母点についての連結線上にその中点を求め、当該連結線に対して直交し、当該連結線と同じ長さをもつ皺構成線分を、当該皺構成線分の中点が当該連結線の中点に重なるように配置する皺構成線分配置部と、
前記皺構成線分の集合体に基づいて、皺画像を生成する皺画像生成部と、
を備えることを特徴とする皺画像生成装置。
【請求項２】
請求項１に記載の皺画像生成装置において、
仮想ばね定義部が、ばね定数ｋおよび自然長Ｌを物理的特性として有する仮想ばねを定義し、前記ばね定数ｋおよび前記自然長Ｌのいずれか一方もしくは双方を、当該仮想ばねが配置される連結線の位置に対応する位置に定義されたスカラー値に基づいて決定することを特徴とする皺画像生成装置。
【請求項３】
請求項２に記載の皺画像生成装置において、
スカラー場定義部が、参照領域内に配列された画素の集合からなる参照画像の画像データを格納しており、個々の画素の画素値を当該画素位置についてのスカラー値とすることにより、スカラー場の定義を行い、
仮想ばね定義部が、参照領域内の対応位置に連結線を描画したときに、当該連結線上もしくはその近傍に位置する画素の画素値の和もしくは平均に基づいて仮想ばねの物理的特性値を決定することを特徴とする皺画像生成装置。
【請求項４】
二次元平面上の作業領域内に多数の母点を定義する母点定義部と、
前記作業領域内に定義された個々の母点について、隣接する母点との間にそれぞれ連結線を定義する連結線定義部と、
前記作業領域に対応する参照領域内にスカラー場を定義するスカラー場定義部と、
前記個々の連結線上に、所定の物理的特性をもった仮想ばねを定義する仮想ばね定義部と、
個々の母点を前記スカラー場上の当該母点の位置に対応する位置に定義されたスカラー値に応じた質量をもった質点とみなし、各連結線の両端に位置する一対の母点が当該連結線上に定義された仮想ばねで連結されているとみなしたときに、各仮想ばねによって作用する力に基く各質点の運動シミュレーションを実行するシミュレーション実行部と、
前記運動シミュレーションによって運動中の母点の所定時点における位置を変位母点の位置と決定する変位母点決定部と、
各変位母点についての連結線上にその中点を求め、当該連結線に対して直交し、当該連結線と同じ長さをもつ皺構成線分を、当該皺構成線分の中点が当該連結線の中点に重なるように配置する皺構成線分配置部と、
前記皺構成線分の集合体に基づいて、皺画像を生成する皺画像生成部と、
を備えることを特徴とする皺画像生成装置。
【請求項５】
請求項４に記載の皺画像生成装置において、
スカラー場定義部が、参照領域内に配列された画素の集合からなる参照画像の画像データを格納しており、個々の画素の画素値を当該画素位置についてのスカラー値とすることにより、スカラー場の定義を行い、
シミュレーション実行部が、参照領域内の対応位置に母点を配置したときに、当該母点上の画素の画素値、または当該母点の近傍に位置する複数の画素の画素値の和もしくは平均に基づいて、当該母点の質量を決定することを特徴とする皺画像生成装置。
【請求項６】
請求項１〜５のいずれかに記載の皺画像生成装置において、
母点定義部が、縦横格子状に母点を規則的に初期配置して定義することを特徴とする皺画像生成装置。
【請求項７】
請求項６に記載の皺画像生成装置において、
母点定義部が、規則的に初期配置した母点の各位置を、乱数を用いて変動させることにより、揺らぎの成分をもった配置がなされた母点の定義を行うことを特徴とする皺画像生成装置。
【請求項８】
請求項６または７に記載の皺画像生成装置において、
連結線定義部が、初期配置時の格子配列において、縦横４方向もしくは斜めも含めた８方向に隣接する母点との間に、それぞれ連結線の定義を行うことを特徴とする皺画像生成装置。
【請求項９】
請求項１〜８のいずれかに記載の皺画像生成装置において、
スカラー場定義部が、線状模様をもった参照画像によってスカラー場を定義することを特徴とする皺画像生成装置。
【請求項１０】
請求項９に記載の皺画像生成装置において、
スカラー場定義部が、
空間周波数を示すパラメータＦを設定するとともに、角度の範囲「θ１〜θ２」、幅の範囲「Ｕ１〜Ｕ２」、長さの範囲「Ｈ１〜Ｈ２」をパラメータとして設定するパラメータ設定部と、
参照領域内にＸＹ二次元座標系を定義し、−０．５〜＋０．５の範囲内のランダム値Rndと整数ｉおよびｊを用いて、Ｘ＝（Ｆ×ｉ）＋（Ｆ×Rnd）およびＹ＝（Ｆ×ｊ）＋（Ｆ×Rnd）なる演算によって求まる座標値（Ｘ，Ｙ）で示される位置に核点Ｐ（Ｘ，Ｙ）をプロットする処理を、ｉの値を１から所定の最大値Ｉまで１ずつ変え、ｊの値を１から所定の最大値Ｊまで１ずつ変えた各組み合わせについて実行し、合計（Ｉ×Ｊ）個の核点Ｐを設定する核点設定部と、
個々の核点Ｐについて、θ１〜θ２内のランダム値θ、Ｕ１〜Ｕ２内のランダム値Ｕ、Ｈ１〜Ｈ２内のランダム値Ｈを定め、核点Ｐを中心点として、基準となる座標軸に対して角度θをなす方向を向き、幅Ｕ，長さＨをもった参照帯状領域を定める処理を行い、合計（Ｉ×Ｊ）個の参照帯状領域を設定する参照帯状領域設定部と、
前記各参照帯状領域の外部の画素には背景部分を示す画素値を与え、前記各参照帯状領域の内部の画素には線状模様部分を示す画素値を与えることにより、線状模様をもった参照画像を生成する参照画像生成部と、
を有することを特徴とする皺画像生成装置。
【請求項１１】
請求項１０に記載の皺画像生成装置において、
参照画像生成部が、内側から外側に向かって画素値がなだらかに変化するように、各参照帯状領域の内部の画素に画素値を与えることを特徴とする皺画像生成装置。
【請求項１２】
請求項１〜８のいずれかに記載の皺画像生成装置において、
スカラー場定義部が、実在の物体表面の画像を取り込む画像取込装置を有し、取り込んだ画像もしくはこれを加工して得られた画像によってスカラー場を定義することを特徴とする皺画像生成装置。
【請求項１３】
請求項１〜１２のいずれかに記載の皺画像生成装置において、
シミュレーション実行部が、両母点間の距離Ｄと当該両母点間に連結される仮想ばねの自然長Ｌとの差を（Ｄ−Ｌ）とし、当該両母点間に当該仮想ばねに起因したＦ＝−ｋ（Ｄ−Ｌ）なる力が作用するものとし、各母点の質量をｍとしてＦ＝ｍαなる式に基づいて各母点の運動加速度αを求め、各母点の運動シミュレーションを実行することを特徴とする皺画像生成装置。
【請求項１４】
請求項１３に記載の皺画像生成装置において、
シミュレーション実行部が、第ｋ番目の時点ｔｋにおいて、各母点のそれぞれについて、その時点の仮想ばねに起因して作用する力の合力を求め、第ｋ番目の時点ｔｋから第（ｋ＋１）番目の時点ｔ（ｋ＋１）に至るまでの時間δの間における各母点の前記合力に基づく移動経路を求め、時点ｔ（ｋ＋１）における各母点の位置を決定する処理を、ｋ＝０からｋを１ずつ増加させながら繰り返し実行し、
変位母点決定部が、ｋの値が所定値Ｋに達した時点における各母点の位置を変位母点の位置と決定することを特徴とする皺画像生成装置。
【請求項１５】
請求項１〜１４のいずれかに記載の皺画像生成装置において、
シミュレーション実行部が、作業領域の輪郭近傍に位置する母点を不動の固定点として、運動シミュレーションを実行することを特徴とする皺画像生成装置。
【請求項１６】
請求項１〜１５のいずれかに記載の皺画像生成装置において、
シミュレーション実行部が、特定の両母点間に作用する力Ｆが所定のしきい値Ｆthを超えた場合に、当該特定の両母点間を連結する仮想ばねが破断したものと判断し、以後、当該特定の両母点間には仮想ばねが存在しないものとして、運動シミュレーションを実行することを特徴とする皺画像生成装置。
【請求項１７】
請求項１〜１６のいずれかに記載の皺画像生成装置において、
皺画像生成部が、各皺構成線分に幅をもたせる処理を行い、皺構成帯状領域を形成し、各皺構成帯状領域の外部の画素には背景部分を示す画素値を与え、前記各皺構成帯状領域の内部の画素には皺部分を示す画素値を与えることにより皺画像を生成することを特徴とする皺画像生成装置。
【請求項１８】
請求項１７に記載の皺画像生成装置において、
皺画像生成部が、幅の範囲「Ｗ１〜Ｗ２」をパラメータとして設定し、個々の皺構成線分について、それぞれＷ１〜Ｗ２内のランダム値Ｗを定め、幅Ｗをもった皺構成帯状領域を形成することを特徴とする皺画像生成装置。
【請求項１９】
請求項１７または１８に記載の皺画像生成装置において、
皺画像生成部が、内側から外側に向かって画素値がなだらかに変化するように、各皺構成帯状領域の内部の画素に画素値を与えることを特徴とする皺画像生成装置。
【請求項２０】
請求項１〜１９のいずれかに記載の皺画像生成装置としてコンピュータを機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、皺画像生成装置に関し、特に、皮革の表面に特有の革シボパターンを構成する多数の皺からなる画像を人為的に作成するのに適した装置に関する。
【背景技術】
【０００２】
皮革は、財布、鞄、靴などの身のまわりの物から、ソファーなどの家具、そして自動車の内装材に至るまで、様々な製品に利用されている。この皮革は、多数の皺からなる独特の模様を有しており、この皮革特有の凹凸模様は、一般に「革シボ」と呼ばれている。最近は、天然の皮革を用いる代わりに、様々な合成材料を用いて作成した人工皮革も盛んに利用されるようになってきているが、人工皮革では、表面に革シボパターンを人工的に形成する必要がある。
【０００３】
このような革シボパターンは、自然界の生物によって生み出された自然な風合いをもった独特のパターンであるため、従来は、天然皮革から抽出したパターンを転用して作成するのが一般的であった。具体的には、天然皮革の表面を写真撮影し、当該撮影画像の中から意匠的に優れた部分を切り出し、これを繰り返し配置して、繋ぎ目を目立たなくするシームレス処理を施すことにより、必要な大きさの革シボパターンを作成することができる。
【０００４】
一方、最近では、ＣＧ技術の発達により、コンピュータを用いた演算処理のみで、自然な風合いをもった革シボパターンの皺画像を作成する技術も提案されている。たとえば、下記の非特許文献１には、ランダムに配置した母点に基づくボロノイ分割を利用して革シボパターンを疑似的に生成する手法が開示されており、下記の特許文献１には、同様の手法で爬虫類の皮膚や植物の組織などの表現に適した不定形セルからなる皺模様を生成する手法が開示されている。また、下記の非特許文献２には、生物の皮膚細胞をパーティクルに見立て、パーティクルの配向を制御するポテンシャル場を定義し、Runge-Kutta法を用いて皮革細胞の成長シミュレーションを行い、細胞形状をメタボールで表現することにより革シボパターンを疑似的に生成する手法が開示されている。
【特許文献１】特開平９−２６５５４２号公報
【非特許文献１】石井、安田、横井、鳥脇、”表面の微細形状に注目した皮膚の質感表現の一手法”、情処論誌、Vol. 32、No. 5、pp.645-654 (1991)
【非特許文献２】宮田、坂口、今尾、須崎、”パーティクルとメタボールを用いた皮革テクスチャの生成法”、情処研報、Vol. 2006、No. 119、pp.13-18 (2006)
【発明の開示】
【発明が解決しようとする課題】
【０００５】
上述した天然皮革の撮影画像を利用して革シボパターンを作成する従来の手法は、撮影時の照明環境と天然皮革の凹凸模様の位置関係により、撮影画像上に光沢や陰影が付加されてしまうため、天然皮革の凹凸模様を忠実に再現することが困難である。また、革シボパターンの特徴は、用いた天然皮革の特徴に依存したものになり、用途に応じた嗜好に合致した革シボパターン等の皺画像を作成することもできない。
【０００６】
これに対して、前掲の各文献に記載されたＣＧ技術を用いた手法では、パラメータの設定値を変えることにより、ある程度、好みに応じた革シボパターン等の自然な風合いをもった皺画像を作成することが可能になる。しかしながら、設定可能なパラメータの数は限られているため、自由度の高い皺画像（特に異方性をもった画像）の作成は困難である。
【０００７】
そこで本発明は、ＣＧ技術を用いた手法により、自然な風合いをもちつつ、非常に自由度の高い革シボパターン等の皺画像を生成できる皺画像生成装置を提供することを目的とする。
【課題を解決するための手段】
【０００８】
(1)   本発明の第１の態様は、皺画像生成装置において、
二次元平面上の作業領域内に多数の母点を定義する母点定義部と、
作業領域内に定義された個々の母点について、隣接する母点との間にそれぞれ連結線を定義する連結線定義部と、
作業領域に対応する参照領域内にスカラー場を定義するスカラー場定義部と、
個々の連結線上に、スカラー場上の当該連結線の位置に対応する位置に定義されたスカラー値に応じた物理的特性をもった仮想ばねを定義する仮想ばね定義部と、
個々の母点を質点とみなし、各連結線の両端に位置する一対の母点が当該連結線上に定義された仮想ばねで連結されているとみなしたときに、各仮想ばねによって作用する力に基く各質点の運動シミュレーションを実行するシミュレーション実行部と、
運動シミュレーションによって運動中の母点の所定時点における位置を変位母点の位置と決定する変位母点決定部と、
各変位母点についての連結線上にその中点を求め、当該連結線に対して直交し、当該連結線と同じ長さをもつ皺構成線分を、当該皺構成線分の中点が当該連結線の中点に重なるように配置する皺構成線分配置部と、
皺構成線分の集合体に基づいて、皺画像を生成する皺画像生成部と、
を設けるようにしたものである。
【０００９】
(2)   本発明の第２の態様は、上述の第１の態様に係る皺画像生成装置において、
仮想ばね定義部が、ばね定数ｋおよび自然長Ｌを物理的特性として有する仮想ばねを定義し、ばね定数ｋおよび自然長Ｌのいずれか一方もしくは双方を、当該仮想ばねが配置される連結線の位置に対応する位置に定義されたスカラー値に基づいて決定するようにしたものである。
【００１０】
(3)   本発明の第３の態様は、上述の第２の態様に係る皺画像生成装置において、
スカラー場定義部が、参照領域内に配列された画素の集合からなる参照画像の画像データを格納しており、個々の画素の画素値を当該画素位置についてのスカラー値とすることにより、スカラー場の定義を行い、
仮想ばね定義部が、参照領域内の対応位置に連結線を描画したときに、当該連結線上もしくはその近傍に位置する画素の画素値の和もしくは平均に基づいて仮想ばねの物理的特性値を決定するようにしたものである。
【００１１】
(4)   本発明の第４の態様は、皺画像生成装置において、
二次元平面上の作業領域内に多数の母点を定義する母点定義部と、
作業領域内に定義された個々の母点について、隣接する母点との間にそれぞれ連結線を定義する連結線定義部と、
作業領域に対応する参照領域内にスカラー場を定義するスカラー場定義部と、
個々の連結線上に、所定の物理的特性をもった仮想ばねを定義する仮想ばね定義部と、
個々の母点をスカラー場上の当該母点の位置に対応する位置に定義されたスカラー値に応じた質量をもった質点とみなし、各連結線の両端に位置する一対の母点が当該連結線上に定義された仮想ばねで連結されているとみなしたときに、各仮想ばねによって作用する力に基く各質点の運動シミュレーションを実行するシミュレーション実行部と、
運動シミュレーションによって運動中の母点の所定時点における位置を変位母点の位置と決定する変位母点決定部と、
各変位母点についての連結線上にその中点を求め、当該連結線に対して直交し、当該連結線と同じ長さをもつ皺構成線分を、当該皺構成線分の中点が当該連結線の中点に重なるように配置する皺構成線分配置部と、
皺構成線分の集合体に基づいて、皺画像を生成する皺画像生成部と、
を設けるようにしたものである。
【００１２】
(5)   本発明の第５の態様は、上述の第４の態様に係る皺画像生成装置において、
スカラー場定義部が、参照領域内に配列された画素の集合からなる参照画像の画像データを格納しており、個々の画素の画素値を当該画素位置についてのスカラー値とすることにより、スカラー場の定義を行い、
シミュレーション実行部が、参照領域内の対応位置に母点を配置したときに、当該母点上の画素の画素値、または当該母点の近傍に位置する複数の画素の画素値の和もしくは平均に基づいて、当該母点の質量を決定するようにしたものである。
【００１３】
(6)   本発明の第６の態様は、上述の第１〜第５の態様に係る皺画像生成装置において、
母点定義部が、縦横格子状に母点を規則的に初期配置して定義するようにしたものである。
【００１４】
(7)   本発明の第７の態様は、上述の第６の態様に係る皺画像生成装置において、
母点定義部が、規則的に初期配置した母点の各位置を、乱数を用いて変動させることにより、揺らぎの成分をもった配置がなされた母点の定義を行うようにしたものである。
【００１５】
(8)   本発明の第８の態様は、上述の第６または第７の態様に係る皺画像生成装置において、
連結線定義部が、初期配置時の格子配列において、縦横４方向もしくは斜めも含めた８方向に隣接する母点との間に、それぞれ連結線の定義を行うようにしたものである。
【００１６】
(9)   本発明の第９の態様は、上述の第１〜第８の態様に係る皺画像生成装置において、
スカラー場定義部が、線状模様をもった参照画像によってスカラー場を定義するようにしたものである。
【００１７】
(10)  本発明の第１０の態様は、上述の第９の態様に係る皺画像生成装置において、
スカラー場定義部が、
空間周波数を示すパラメータＦを設定するとともに、角度の範囲「θ１〜θ２」、幅の範囲「Ｕ１〜Ｕ２」、長さの範囲「Ｈ１〜Ｈ２」をパラメータとして設定するパラメータ設定部と、
参照領域内にＸＹ二次元座標系を定義し、−０．５〜＋０．５の範囲内のランダム値Rndと整数ｉおよびｊを用いて、Ｘ＝（Ｆ×ｉ）＋（Ｆ×Rnd）およびＹ＝（Ｆ×ｊ）＋（Ｆ×Rnd）なる演算によって求まる座標値（Ｘ，Ｙ）で示される位置に核点Ｐ（Ｘ，Ｙ）をプロットする処理を、ｉの値を１から所定の最大値Ｉまで１ずつ変え、ｊの値を１から所定の最大値Ｊまで１ずつ変えた各組み合わせについて実行し、合計（Ｉ×Ｊ）個の核点Ｐを設定する核点設定部と、
個々の核点Ｐについて、θ１〜θ２内のランダム値θ、Ｕ１〜Ｕ２内のランダム値Ｕ、Ｈ１〜Ｈ２内のランダム値Ｈを定め、核点Ｐを中心点として、基準となる座標軸に対して角度θをなす方向を向き、幅Ｕ，長さＨをもった参照帯状領域を定める処理を行い、合計（Ｉ×Ｊ）個の参照帯状領域を設定する参照帯状領域設定部と、
各参照帯状領域の外部の画素には背景部分を示す画素値を与え、各参照帯状領域の内部の画素には線状模様部分を示す画素値を与えることにより、線状模様をもった参照画像を生成する参照画像生成部と、
を有するようにしたものである。
【００１８】
(11)  本発明の第１１の態様は、上述の第１０の態様に係る皺画像生成装置において、
参照画像生成部が、内側から外側に向かって画素値がなだらかに変化するように、各参照帯状領域の内部の画素に画素値を与えるようにしたものである。
【００１９】
(12)  本発明の第１２の態様は、上述の第１〜第８の態様に係る皺画像生成装置において、
スカラー場定義部が、実在の物体表面の画像を取り込む画像取込装置を有し、取り込んだ画像もしくはこれを加工して得られた画像によってスカラー場を定義するようにしたものである。
【００２０】
(13)  本発明の第１３の態様は、上述の第１〜第１２の態様に係る皺画像生成装置において、
シミュレーション実行部が、両母点間の距離Ｄと当該両母点間に連結される仮想ばねの自然長Ｌとの差を（Ｄ−Ｌ）とし、当該両母点間に当該仮想ばねに起因したＦ＝−ｋ（Ｄ−Ｌ）なる力が作用するものとし、各母点の質量をｍとしてＦ＝ｍαなる式に基づいて各母点の運動加速度αを求め、各母点の運動シミュレーションを実行するようにしたものである。
【００２１】
(14)  本発明の第１４の態様は、上述の第１３の態様に係る皺画像生成装置において、
シミュレーション実行部が、第ｋ番目の時点ｔｋにおいて、各母点のそれぞれについて、その時点の仮想ばねに起因して作用する力の合力を求め、第ｋ番目の時点ｔｋから第（ｋ＋１）番目の時点ｔ（ｋ＋１）に至るまでの時間δの間における各母点の合力に基づく移動経路を求め、時点ｔ（ｋ＋１）における各母点の位置を決定する処理を、ｋ＝０からｋを１ずつ増加させながら繰り返し実行し、
変位母点決定部が、ｋの値が所定値Ｋに達した時点における各母点の位置を変位母点の位置と決定するようにしたものである。
【００２２】
(15)  本発明の第１５の態様は、上述の第１〜第１４の態様に係る皺画像生成装置において、
シミュレーション実行部が、作業領域の輪郭近傍に位置する母点を不動の固定点として、運動シミュレーションを実行するようにしたものである。
【００２３】
(16)  本発明の第１６の態様は、上述の第１〜第１５の態様に係る皺画像生成装置において、
シミュレーション実行部が、特定の両母点間に作用する力Ｆが所定のしきい値Ｆthを超えた場合に、当該特定の両母点間を連結する仮想ばねが破断したものと判断し、以後、当該特定の両母点間には仮想ばねが存在しないものとして、運動シミュレーションを実行するようにしたものである。
【００２４】
(17)  本発明の第１７の態様は、上述の第１〜第１６の態様に係る皺画像生成装置において、
皺画像生成部が、各皺構成線分に幅をもたせる処理を行い、皺構成帯状領域を形成し、各皺構成帯状領域の外部の画素には背景部分を示す画素値を与え、各皺構成帯状領域の内部の画素には皺部分を示す画素値を与えることにより皺画像を生成するようにしたものである。
【００２５】
(18)  本発明の第１８の態様は、上述の第１７の態様に係る皺画像生成装置において、
皺画像生成部が、幅の範囲「Ｗ１〜Ｗ２」をパラメータとして設定し、個々の皺構成線分について、それぞれＷ１〜Ｗ２内のランダム値Ｗを定め、幅Ｗをもった皺構成帯状領域を形成するようにしたものである。
【００２６】
(19)  本発明の第１９の態様は、上述の第１７または第１８の態様に係る皺画像生成装置において、
皺画像生成部が、内側から外側に向かって画素値がなだらかに変化するように、各皺構成帯状領域の内部の画素に画素値を与えるようにしたものである。
【００２７】
(20)  本発明の第２０の態様は、上述の第１〜第１９の態様に係る皺画像生成装置としてコンピュータを機能させるためのプログラムを用意するようにしたものである。
【発明の効果】
【００２８】
本発明の皺画像生成装置では、多数の母点間に仮想ばねを連結し、コンピュータ上で力学的シミュレーションを行うことにより、革シボパターン等の皺画像を生成することができる。各母点間に仮想ばねに基づく力を作用させるシミュレーションを行うため、生物の皮膚の成長過程に似た現象をシミュレートすることができるので、ＣＧ技術を用いた手法を採りながら、天然皮革などの自然な風合いをもつ皺画像の生成が可能になる。また、力学的シミュレーションのパラメータとなるばねの物理的特性や母点の質量を、二次元画像として与えられるスカラー場によって設定することができるため、非常に自由度の高い皺画像の生成が可能になる。
【発明を実施するための最良の形態】
【００２９】
以下、本発明を図示する実施形態に基づいて説明する。
【００３０】
＜＜＜  §１．本発明に係る皺画像生成手法の基本概念  ＞＞＞
本発明に係る皺画像生成装置は、主として天然皮革に特有な革シボパターンを、コンピュータを用いた演算処理によって生成するために開発されたものである。したがって、この演算処理のアルゴリズムの根底には、革シボパターンに固有の特徴を表現する上で好都合な要素が盛り込まれている。そこで、まず、本発明の着想を得るきっかけとなった一般的な革シボパターンに固有の特徴を簡単に述べておく。
【００３１】
既に述べたとおり、「革シボパターン」とは、皮革特有の凹凸模様、すなわち、多数の皺からなる独特の模様であるが、その実体は、動物の種類によって様々であり、また、同じ動物でも皮膚の部位によって様々である。ただ、一般的には、細かな不定形のセル状要素の集合体として把握できるパターンであり、個々のセルの境界部分に皺が生じることになる。
【００３２】
図１は、革シボパターンの基本構成例を示す平面図である。この平面図は、説明の便宜上、出願人が作成した人為的な図であり、個々のセルがほぼ完全な閉領域を構成している。すなわち、図の白い部分が個々のセルを構成し、黒い部分がその境界部分となっている。もちろん、実際の天然皮革の表面に現れる革シボパターンは、このように個々のセルが完全な閉領域を構成するものにはならず、白い部分と黒い部分とが、かなり乱雑に入り乱れたものになるが、ここでは、図に白い部分として示される皮丘部と、図に黒い部分として示される皮溝部とを対比する上で、個々のセル（皮丘部）とその境界部分（皮溝部）という明確な区分けがなされた例を示すことにする。
【００３３】
また、実際の革シボパターンは、平面的な模様ではなく、凹凸構造をもった立体模様である。すなわち、図示の黒い部分は溝を形成し、図示の白い部分は丘を形成する。白い部分（セルの部分）を皮丘部と呼び、黒い部分（セルの境界部分）を皮溝部と呼ぶのは、このような立体構造において、両者が溝と丘の関係をなすためである。また、図１にも例示したとおり、このような皮溝部がある程度の長さにわたって連なって見える箇所は、一般に血筋と呼ばれている。
【００３４】
なお、図１では、便宜上、白黒の二値画像を示したが、実際の天然皮革の凹凸構造は、凹と凸の２通りの高さをもつわけではなく、なだらかな高低をもつ構造であり、たとえば、皮溝部の部分であっても、浅い部分から深い部分に至るまで様々な深さを有している。このような連続的な高低情報を示すには、図１のような二値画像ではなく、濃淡の情報をもった階調画像が用いられる。
【００３５】
このように、革シボパターンに、皮丘部や皮溝部が生じるのは、皮革表面の生体構造に由来する。図２は、天然皮革表面の一般的な構造を示す拡大斜視図である。皮革は生物の皮膚の表面部分であり、図示のとおり、皮丘と皮溝とによって構成されている。そして、通常、皮溝の交差部分に汗腺や毛根が位置し、そこから体毛が育成されることになる。図１に示す皮丘部および皮溝部は、図２に示す皮丘および皮溝に対応する部分である。
【００３６】
天然皮革の革シボパターンに自然の風合いが感じられるのは、個々の皮丘部や皮溝部の形状・大きさ・向きなどに、自然の揺らぎ要素が盛り込まれているためと考えられる。生物の皮膚の各部位によって、このような揺らぎが生じる理由は、生物学的な探究のテーマであるが、いずれにしても、このような揺らぎが生じるのは、生物の成長過程において、各部分の成長速度がそれぞれ異なるためと考えられる。本願発明者が知る限り、皮膚の各部位によって成長速度に差が生じる正確なメカニズムは、現段階では解明されていないようである。ただ、一般論的には、何らかの要因で皮膚細胞内の保湿力に差が生じ、保湿力の高い部位は芳しい成長を遂げるが、保湿力の低い部位は緩慢な成長を遂げる、とする説が唱えられている。
【００３７】
天然皮革の表面に生じる革シボパターンを、人為的に生成するためには、本来、このような生物学的な成長過程を詳細に把握し、当該生物学的な成長過程をコンピュータ上でシミュレートするべきである。しかしながら、現時点では、このような生物学的な成長過程の正確なメカニズムは不明であり、そのようなシミュレーションを行うことは非常に困難である。そこで、本願発明者は、このような生物学的な皮膚の成長過程を、より単純な物理学的現象に置き換えることにより、革シボパターンに類似したパターンを人為的に発生できるのではないかと考えた。
【００３８】
本願発明者が着目した物理学的現象は、相互にばねで連結された多数の質点の挙動である。たとえば、規則的に配列された多数の質点を用意し、個々の質点について、その四方もしくは八方に隣接する別な質点との間にばねを連結したものとする。このとき、個々のばねを強いばねにしたり、弱いばねにしたりすれば、強いばねで連結された２つの質点は相互に接近しようとするが、弱いばねで連結された２つの質点は、結果的に相互の距離が離されてしまうことになろう。つまり、ある領域については強いばねを用い、別なある領域については弱いばねを用いる、というように、領域ごとにばねの条件を変えてやれば、質点を密集させたり、離散させたり制御することができる。
【００３９】
本願発明者は、このように、相互にばねで連結された多数の質点の挙動に関する物理学的なシミュレーションは、前述した生体の皮膚の成長過程のシミュレーションと似た結果を生む可能性があるのではないかと考えた。すなわち、生体の皮膚の成長過程において、保湿力の高い部位は著しい成長を遂げ、保湿力の低い部位は緩慢な成長を遂げる、とすれば、これを質点とばねのモデルに置き換え、保湿力をばね定数として取り扱えば、類似した結果が得られるのではないかと考えたのである。
【００４０】
このような基本概念に基づいて質点とばねのモデルについてのシミュレーションを行った結果、自然な風合いを有する疑似的な革シボパターンの皺画像を生成することに成功した。しかも、この手法では、生物学的なモデルにおける「保湿力分布」に相当するパラメータとなる「ばね定数分布」を二次元スカラー場（二次元画像）として用意することができるので、この二次元スカラー場を変えることにより、非常に自由度の高い皺画像の生成が可能になる。
【００４１】
＜＜＜  §２．本発明に係る皺画像生成装置の基本構成  ＞＞＞
図３は、本発明に係る皺画像生成装置の基本構成を示すブロック図である。図示のとおり、この装置は、母点定義部１０、連結線定義部２０、仮想ばね定義部３０、スカラー場定義部４０、シミュレーション実行部５０、変位母点決定部６０、皺構成線分配置部７０、皺画像生成部８０によって構成されている。もっとも、実際には、これらの各構成要素は、コンピュータに専用のプログラムを組み込むことによって実現されるものである。なお、図３に示す実線の矢印は、各構成要素間の情報の流れを示す。また、破線の矢印は、§６で述べる変形例における情報の流れを示している。
【００４２】
母点定義部１０は、二次元平面上の作業領域内に多数の母点を定義する機能を有している。ここで定義した個々の母点は、所定の質量をもった質点として取り扱われることになる。連結線定義部２０は、この作業領域内に定義された個々の母点について、隣接する母点との間にそれぞれ連結線を定義する機能を有している。
【００４３】
仮想ばね定義部３０は、この個々の連結線上に仮想ばねを定義する処理を行うが、このとき、スカラー場定義部４０によって定義されたスカラー場を参照して、個々の仮想ばねの物理的特性を決定する。スカラー場定義部４０は、作業領域に対応する参照領域内に二次元のスカラー場を定義する機能を有している。ここで定義される二次元のスカラー場は、二次元平面上の個々の位置に何らかのスカラー値を特定することができるような場であれば、どのような場であってもかまわないが、実用上は、後述するように、二次元画像を用いてスカラー場の定義がなされる。ここで定義されるスカラー場は、§１で述べた生物学的なモデルにおける「保湿力分布」を示す場である。
【００４４】
仮想ばね定義部３０は、連結線定義部２０によって定義された個々の連結線上に、スカラー場定義部４０によって定義されたスカラー場上の当該連結線の位置に対応する位置に定義されたスカラー値に応じた物理的特性をもった仮想ばねの定義を行うことになる。具体的には、後述するように、スカラー値に応じて、仮想ばねの「ばね定数ｋ」もしくは「自然長Ｌ」が設定されることになる。§１で述べた生物学的なモデルとの関係においては、個々の場所における「保湿力」に応じて、仮想ばねの物理的特性が設定されることになる。
【００４５】
シミュレーション実行部５０は、母点定義部１０が定義した個々の母点を質点とみなし、各連結線の両端に位置する一対の母点が当該連結線上に定義された仮想ばねで連結されているとみなしたときに、各仮想ばねによって作用する力に基く各質点の運動シミュレーションを実行する。そして、変位母点決定部６０は、この運動シミュレーションによって運動中の母点の所定時点における位置を変位母点の位置と決定する。
【００４６】
皺構成線分配置部７０は、個々の変位母点の位置に基づいて、作業領域内に皺構成線分を配置する処理を実行する。具体的には、各変位母点についての連結線上にその中点を求め、当該連結線に対して直交し、当該連結線と同じ長さをもつ皺構成線分を、当該皺構成線分の中点が当該連結線の中点に重なるように配置する。そして、皺画像生成部８０は、この皺構成線分の集合体に基づいて、皺画像を生成する。生成された皺画像は、皺画像データとして出力される。
【００４７】
以上、図３のブロック図を参照しながら、本発明に係る皺画像生成装置の基本構成を説明したが、図にブロックとして示される個々の構成要素の詳細な動作については、以下に述べる§３以降で実例を挙げながら説明する。
【００４８】
＜＜＜  §３．母点および連結線の定義  ＞＞＞
ここでは、図３に示す母点定義部１０によって定義される母点および連結線定義部２０で定義される連結線について、実例を挙げながら、より詳細な説明を行う。
【００４９】
図４は、図３に示す装置の母点定義部１０で定義された母点Ｍの一例を示す平面図である。上述したとおり、母点定義部１０は、二次元平面上の作業領域内に多数の母点を定義する機能を有している。図４に示す例では、ハッチングを施して示す作業領域Ａ内に、５行５列の規則的な格子点配置をもった母点Ｍ（白丸で示す）が定義された例が示されている。すなわち、図示の例では、縦横の寸法がａである正方形からなる作業領域Ａ内に、ピッチｄの間隔で合計２５個の母点Ｍが配置されている。
【００５０】
ここでは、図示の便宜上、５行５列という小規模な母点配置を行った例を示すが、実際には、たとえば、１００行１００列といったより大きな配列上に母点の定義が行われる。ピッチｄは、図１に示す皮シボパターンにおける皮丘部のピッチを左右するパラメータとなるので、最終的に生成させたい皺画像に応じて、適当な寸法に設定すればよい。また、作業領域Ａは、皺画像の形成領域に対応する領域であるので、その寸法ａは、最終的に生成させたい皺画像のサイズに応じて、適当な寸法に設定すればよい。ピッチｄおよび寸法ａを定めれば、母点定義部１０は、作業領域Ａ内に多数の母点Ｍを自動的に定義することができる。
【００５１】
なお、図示の例では、作業領域Ａを正方形状の領域にしているが、作業領域Ａは長方形であってもよいし、矩形以外の任意の形状であってもかまわない。ピッチｄも縦横で同一にする必要はなく、縦ピッチと横ピッチとを別個に設定してもよい。あるいは、六角形を平面状に隙間なく敷き詰め、各六角形の頂点位置に母点を配置するような形態をとってもかまわない。また、図４に示す例では、縦横格子状に母点Ｍを規則的に配置しているが、母点Ｍは必ずしも規則的に配置する必要はなく、乱数を用いてランダムな位置に母点を配置するようにしてもかまわない。
【００５２】
要するに、母点定義部１０は、所定の二次元作業領域内に、全体的にほぼ一様に分布するように、多数の母点Ｍを定義することができれば、どのようなアルゴリズムに基づいて母点の定義を行ってもかまわない。
【００５３】
もちろん、図４に示す例のように、縦横格子状に母点Ｍを規則的に初期配置して定義した後、この規則的に初期配置した母点の各位置を、乱数を用いて変動させることにより、揺らぎの成分をもった配置がなされた母点の定義を行うことも可能である。たとえば、−０．５〜＋０．５の間のランダムな値を示す乱数Rndを生成し、１つの母点Ｍについて、横方向の変動成分ΔＸを、ΔＸ＝(ｄ×Rnd）と定義し、縦方向の変動成分ΔＹを、ΔＹ＝(ｄ×Rnd）と定義し、規則的に初期配置された母点Ｍの位置を、右方向にΔＸだけ移動させ（ΔＸが負の場合は左方向へ移動させる）、上方向にΔＹだけ移動させると（ΔＹが負の場合は下方向へ移動させる）、規則的な初期配置がランダムに乱されることになる。その結果、作業領域Ａを全体的に観察した場合、母点Ｍはほぼ一様に分布しているが、個々の母点は規則的ではなく、ランダムな揺らぎの成分をもった配置をとることになる。
【００５４】
後述する§５の力学的な運動シミュレーションを行うことにより、個々の母点Ｍは変位することになるので、母点定義部１０によって規則的な母点の定義を行ったとしても、結局、変位後の母点の分布の規則性は失われることになる。しかしながら、母点定義部１０によって定義された各母点の配置は、運動シミュレーションにおける最初の母点位置を与えることになるので、より揺らぎ成分の強い皺画像を生成したい場合には、母点定義部１０で母点を定義する段階から、不規則配置をとる母点を定義しておくのが好ましい。
【００５５】
なお、図４では、説明の便宜上、平面図として各母点Ｍの配置を示したが、母点定義部１０による母点定義の実際の作業は、個々の母点についての座標値を決定する作業ということになる。すなわち、１つの母点Ｍの座標値を、Ｍ（ｘ，ｙ）のようなＸ座標値とＹ座標値との組み合わせで表現することにすれば、図４に示す母点定義は、全２５個の母点のそれぞれについて、Ｍ（ｘ，ｙ）のような座標値を算出する処理として実行され、個々の座標値は、各母点の位置を示すデータとして、母点定義部１０内に格納されることになる。具体的には、全２５個の母点をＭ１〜Ｍ２５とすれば、母点定義部１０内には、Ｍ１（ｘ，ｙ），Ｍ２（ｘ，ｙ），... ，Ｍ２５（ｘ，ｙ）といったデータが用意されることになる。
【００５６】
続いて、連結線定義部２０の機能を説明する。この連結線定義部２０は、作業領域Ａ内に定義された個々の母点Ｍについて、隣接する母点との間にそれぞれ連結線を定義する処理を実行する。図５は、図３に示す装置の連結線定義部２０で定義された連結線Ｃの一例を示す平面図である。図に実線で示す線分が連結線Ｃを示しており、この例では、縦横の他、斜め方向に隣接する母点との間にも連結線Ｃが定義されている。
【００５７】
図６は、特定の着目母点Ｍ０（ここでは、便宜上、二重丸で示す）について定義された８本の連結線Ｃ１〜Ｃ８を示す平面図である。図示のとおり、この例では、着目母点Ｍ０の縦横４方向に隣接する４個の母点Ｍ２，Ｍ４，Ｍ５，Ｍ７と、斜め４方向に隣接する４個の母点Ｍ１，Ｍ３，Ｍ６，Ｍ８という合計８個の隣接母点との間に、それぞれ連結線Ｃ１〜Ｃ８が定義されている。
【００５８】
もちろん、連結線の定義は、縦横４方向に隣接する４個の母点Ｍ２，Ｍ４，Ｍ５，Ｍ７との間の４本だけを定義するようにしてもよいし、図示の例のように、斜めも含めた８方向に隣接する母点との間の８本を定義するようにしてもよい。あるいは、母点定義部１０による母点定義の形態が、平面状に隙間なく敷き詰めた六角形の頂点位置に母点を配置するような形態を採った場合は、各六角形の辺をそのまま連結線とする定義を行うこともできる。
【００５９】
また、母点定義部１０による母点定義の方法として、まず、母点を規則的に格子状に並べて初期配置とした後、乱数を用いて各母点の位置を変動させ、揺らぎの成分をもった配置をなす母点を得る方法を採用した場合は、初期配置時の格子配列において、縦横４方向もしくは斜めも含めた８方向に隣接する母点との間に、それぞれ連結線の定義を行うようにすればよい。また、母点定義部１０による母点定義の方法として、最初から全くランダムな座標値をもつ母点を発生させる方法を採った場合には、各母点間に規則的な格子を前提とした縦横斜めの関係を定義することはできないが、そのような場合には、個々の母点をそれぞれ着目母点Ｍ０とし、当該着目母点Ｍ０から所定の半径をもった円内に位置する他の母点を隣接母点と認識し、当該着目母点Ｍ０と当該隣接母点との間に連結線を定義するようにすればよい。
【００６０】
なお、図５では、説明の便宜上、各母点Ｍ間を結ぶ線分として各連結線Ｃを示したが、連結線定義部２０による実際の連結線定義作業は、このような連結線を構成する図形としての線分を発生させる作業ではなく、「ある特定の母点と別なある特定の母点との間に、概念的な連結線が存在する」ことを示すデータを作成する作業になる。たとえば、図６に示す例の場合、「連結線Ｃ１を定義する」作業は、「連結線Ｃ１」という図形としての線分を発生させたり、描画したりする作業ではなく、「母点Ｍ０と母点Ｍ１との間に、連結線Ｃ１が存在する」ことを示すデータ、たとえば、「Ｃ１＝Ｍ０，Ｍ１」のようなデータを作成して格納する作業ということになる。
【００６１】
結局、図６に示す連結線Ｃ１〜Ｃ８を定義する作業は、「Ｃ１＝Ｍ０，Ｍ１」，「Ｃ２＝Ｍ０，Ｍ２」，「Ｃ３＝Ｍ０，Ｍ３」，... ，「Ｃ８＝Ｍ０，Ｍ８」といったデータを作成し、これを格納する作業ということになる。このように、各連結線は、２つの母点の連結関係を示す概念的な線であるから、後述する運動シミュレーションにより、個々の母点の位置が変動したとしても、各連結線それ自体は何ら変わることはない。別言すれば、図５に実線で示す各連結線Ｃは、個々の母点Ｍが移動すると、伸びたり縮んだりする仮想の線ということになる。
【００６２】
＜＜＜  §４．スカラー場と仮想ばねの定義  ＞＞＞
続いて、図３に示す装置における仮想ばね定義部３０およびスカラー場定義部４０の機能について詳述する。
【００６３】
スカラー場定義部４０は、作業領域Ａに対応する参照領域内にスカラー場を定義する機能を果たす。このスカラー場の実体は、実質的には、二次元画像というべきものである。ただ、一般に「画像」とは、本来、ディスプレイ画面や紙面といった媒体上に提示し、人間が視覚的に把握することを前提とした対象物を指す言葉であるのに対し、スカラー場定義部４０によって定義されるべきスカラー場は、後述する力学的な運動シミュレーションにおけるパラメータを設定するために利用されるスカラー値の二次元平面上の分布を示すものであるため、ここでは、「画像」という用語の代わりに「スカラー場」という用語を用いている。
【００６４】
図７は、図３に示す装置のスカラー場定義部４０で定義されたスカラー場を示すための参照画像Ｒの一例を示す平面図である。この参照画像Ｒは、一辺の寸法がａの正方形の参照領域内の画像であり、この参照領域の大きさは、図４に示す作業領域Ａの大きさに一致する。したがって、図７の参照画像Ｒを、図４の作業領域Ａに重ねると、両者はぴったりと一致し、作業領域Ａ内の任意の１点は、参照画像Ｒ内の特定の１点に対応することになる。
【００６５】
図７では、参照画像Ｒが二値画像から構成される単純な例が示されている。すなわち、この参照画像Ｒは、ハッチングを施した楕円状の特定領域Ｅと白地の背景領域とによって構成されており、特定領域Ｅ内の画素には画素値”１”が与えられ、背景領域内の画素には画素値”０”が与えられている。別言すれば、この参照画像Ｒは、二次元平面上の各位置に、”０”か”１”かのいずれかのスカラー値を与える二次元スカラー場を定義していることになる。このスカラー値は、§１で述べた生物学的なモデルにおける「保湿力」を示す値になる。たとえば、スカラー値”０”が定義された領域は、高い保湿力が与えられた部分であり、細胞が十分に成長する部分となり、スカラー値”１”が定義された領域は、低い保湿力しか与えられなかった部分であり、細胞の成長が不十分な部分となる（もちろん、これと逆の定義を行ってもよい）。
【００６６】
一方、仮想ばね定義部３０は、連結線定義部２０で定義された個々の連結線Ｃ上に、それぞれ所定の物理的特性をもった仮想ばねの定義を行う。たとえば、図６には、着目母点Ｍ０に関して定義された８本の連結線Ｃ１〜Ｃ８が示されているが、これらの各連結線Ｃ１〜Ｃ８上には、図８に示すように、それぞれ仮想ばねＳ１〜Ｓ８が定義されることになる。結局、図５に実線で示されている個々の連結線Ｃ上に、それぞれ独立した仮想ばねＳが定義されることになる。
【００６７】
一般に「ばね」は、「自然長Ｌ」と「ばね定数ｋ」という物理的特性を有している。そこで、ここでは、個々の仮想ばねＳごとに、それぞれ別個独立した物理的特性の定義を行うことにする。すなわち、仮想ばね定義部３０は、特定の連結線Ｃ上に仮想ばねを定義する際に、スカラー場定義部４０によって定義されたスカラー場上の当該連結線Ｃの位置に対応する位置に定義されたスカラー値に応じた物理的特性をもった仮想ばねの定義を行うのである。
【００６８】
図９は、各仮想ばねＳの物理的特性を決定するために、図５に示す連結線Ｃが定義された作業領域Ａ上に、図７に示す参照画像Ｒを重ね合わせた状態を示す平面図である。このように重ね合わせれば、各連結線Ｃ上に位置するスカラー場のスカラー値を認識することができる。たとえば、図９に太線で示す連結線Ｃ９は、ハッチングを施して示す特定領域Ｅ内に位置するので、この連結線Ｃ９上には、スカラー値”１”に応じた物理的特性をもった仮想ばねＳ９が定義される。これに対して、同じく図９に太線で示す連結線Ｃ１０は、特定領域Ｅ外に位置するので、この連結線Ｃ１０上には、スカラー値”０”に応じた物理的特性をもった仮想ばねＳ１０が定義される。
【００６９】
たとえば、仮想ばねの物理的特性として、「ばね定数ｋ」の値をスカラー値に応じて定義することにし、スカラー値”０”に応じた「ばね定数ｋ」として、ｋ＝４を与え、スカラー値”１”に応じた「ばね定数ｋ」として、ｋ＝８を与えることにすれば、連結線Ｃ９上の仮想ばねＳ９のばね定数ｋは８になり、連結線Ｃ１０上の仮想ばねＳ１０のばね定数ｋは４になる。
【００７０】
なお、図９に示す例を見ればわかるとおり、特定領域Ｅの内側部分（スカラー値”１”）から外側部分（スカラー値”０”）に跨がる連結線も多数存在するので、そのような連結線上の仮想ばねに対しては、ｋ＝４またはｋ＝８のいずれか一方のばね定数を与えることはできない。このような場合は、スカラー値”１”の領域に係る部分とスカラー値”０”に係る部分との按分比に応じて、中間的なばね定数を与えるようにすればよい。
【００７１】
実際には、スカラー場定義部４０には、参照領域内に配列された画素の集合からなる参照画像Ｒの画像データが格納されており、個々の画素の画素値を当該画素位置についてのスカラー値とすることにより、スカラー場の定義が行われている。したがって、たとえば、図７に示す参照画像Ｒも、それぞれ所定の画素値（”１”または”０”）をもつ画素の集合からなるデータである。したがって、仮想ばね定義部３０は、参照領域内の対応位置に連結線を描画したときに、当該連結線上もしくはその近傍に位置する画素の画素値の平均に基づいて、当該連結線上に配置する仮想ばねの物理的特性値（ここで述べる例の場合は、ばね定数ｋの値）を決定すればよい、
【００７２】
図１０は、作業領域Ａに重ね合わせた参照画像Ｒ内の画素の画素値に基づいて、各仮想ばねの物理的特性を決定する具体的な方法の一例を示す平面図である。図に示す細かな升目は、参照画像Ｒを構成する個々の画素を示している。ここで、たとえば、母点Ｍ１１とＭ１２との間の連結線Ｃ１１上に配置する仮想ばねＳ１１の物理的特性値は、連結線Ｃ１１上の６個の画素（連結線Ｃ１１に接触している斜線ハッチングを施した６個の画素）の画素値の平均に基づいて決定すればよい。同様に、母点Ｍ１３とＭ１５との間の連結線Ｃ１２上に配置する仮想ばねＳ１２の物理的特性値は、連結線Ｃ１２上の６個の画素（連結線Ｃ１２に接触しているドットハッチングを施した６個の画素）の画素値の平均に基づいて決定すればよい。
【００７３】
たとえば、連結線Ｃ１２上の６個の画素のうち、３個の画素が画素値”０”、残りの３個の画素が画素値”１”であった場合、その平均値である”０．５”を連結線Ｃ１２の位置に定義されたスカラー値であるものとして取り扱い、スカラー値”０．５”に基づいて、仮想ばねＣ１２のばね定数ｋを決定すればよい。上述した例のように、スカラー値”０”に対してはｋ＝４なるばね定数を与え、スカラー値”１”に対してはｋ＝８なるばね定数を与えるようにした場合は、スカラー値「０〜１」を、ｋの値「４〜８」に線形対応させ、仮想ばねＣ１２のばね定数ｋとしては、スカラー値”０．５”に対応する「６」を与えるようにすればよい。
【００７４】
なお、連結線上に位置する画素の画素値の平均に基づいて仮想ばねのばね定数ｋを決定する代わりに、連結線上に位置する画素の画素値の和に基づいて仮想ばねのばね定数ｋを決定するようにしてもよい。ただ、このように和を用いると、連結線上に位置する画素の数も、スカラー値を決定する要因の１つになるので、一般的に、長い連結線については、より大きなスカラー値が与えられることになる。このように、画素の数に左右されずにスカラー値を決定したい場合は、上述の例のように、和ではなく、平均を用いるようにすればよい。もっとも、図１０に示す例の場合、連結線Ｃ１１の全長に比べて、連結線Ｃ１２の全長の方が長くなっているが、連結線Ｃ１２が画素配列に対して斜め方向の線であるため、いずれも連結線上に位置する画素の数は６個となっている。
【００７５】
また、図１０に示す例では、連結線上に位置する画素の画素値のみを参照するようにしているが、連結線上に位置する画素のみならず、連結線の近傍に位置する画素（たとえば、画素の中心位置と連結線との距離が所定の値以下となるような画素）の画素値の和もしくは平均に基づいて仮想ばねの物理的特性値を決定するようにすることも可能である。
【００７６】
図１１は、図３に示す装置の仮想ばね定義部３０において定義される仮想ばねＳの概念図であり、図１１(a) は、仮想ばねＳが単独で置かれている状態、図１１(b) は、この仮想ばねＳを２つの母点Ｍ１０，Ｍ２０間に接続した状態を示す。図１１(a) に示すように、仮想ばねＳが単独で置かれている状態では、ばねに対する外力は加わっていないので、仮想ばねＳは自然長Ｌの状態を維持している。本発明における仮想ばねＳは、この自然長Ｌとばね定数ｋとを物理的特性値としてもつ仮想のばねである。
【００７７】
このような仮想ばねＳを、２つの母点Ｍ１０，Ｍ２０間に接続すると、図１１(b) に示すように、母点Ｍ１０，Ｍ２０には、それぞれ力Ｆが作用することになる。この力Ｆの大きさは、母点Ｍ１０，Ｍ２０間の距離をＤとすると、フックの法則により、Ｆ＝−ｋ（Ｄ−Ｌ）で示される。すなわち、作用する力の大きさは、ばねの伸び量（Ｄ−Ｌ）に比例する。ここで、ｋはばね定数、Ｌは仮想ばねＳの自然長であり、ｋの前のマイナス符号は、変位が生じた方向と逆向きの力が生じることを示している。図１１に示す例の場合、Ｄ＞Ｌの例であるため、母点Ｍ１０，Ｍ２０間に作用する力の向きは、白矢印で示すとおり互いに引っ張り合う方向になるが、Ｄ＜Ｌの場合、力の向きは逆向きになる。
【００７８】
前述したとおり、ここで述べる実施形態の場合、個々の仮想ばねＳのばね定数ｋは、スカラー場上の対応する位置に定義されたスカラー値に応じて決定される。一方、個々の仮想ばねＳの自然長は、次のように定めている。
【００７９】
まず、図５に示すような母点の初期配置状態（格子状に規則的配置がなされた状態）において、母点間を水平に連結する連結線および垂直に連結する連結線（図６に示す例における連結線Ｃ２，Ｃ４，Ｃ５，Ｃ７）の上に配置される仮想ばね（以下、「縦横連結ばね」という）については、その長さＬを０に設定している。もちろん、自然長Ｌ＝０になるようなばねは実在しないが、ここで定義する仮想ばねは、後述する運動シミュレーションを実施するためにコンピュータ上で定義されるばねであるので、自然長Ｌ＝０なる設定を行っても支障はない。したがって、もし個々の母点Ｍが、図５に示す例のように、一定ピッチｄで規則的に配置されている状態であったとすると、すべての「縦横連結ばね」の伸び量は「ｄ（母点間隔）−０（ばねの自然長）＝ｄ」ということになり、各母点Ｍが、１本の「縦横連結ばね」から受ける力は、Ｆ＝−ｋｄになる。
【００８０】
一方、図５に示すような母点の初期配置状態において、母点間を斜めに連結する連結線（図６に示す例における連結線Ｃ１，Ｃ３，Ｃ６，Ｃ８）の上に配置される仮想ばね（以下、「斜め連結ばね」という）については、その長さＬを「（ｄ×√２）−ｄ」に設定している。したがって、もし個々の母点Ｍが、図５に示す例のように、一定ピッチｄで規則的に配置されている状態であったとすると、すべての「斜め連結ばね」の伸び量は「（ｄ×√２）（母点間隔）−（（ｄ×√２）−ｄ）（ばねの自然長）＝ｄ」ということになり、各母点Ｍが、１本の「斜め連結ばね」から受ける力は、Ｆ＝−ｋｄになる。
【００８１】
結局、上述した例のように、「縦横連結ばね」の自然長Ｌ＝０、「斜め連結ばね」の自然長Ｌ＝「（ｄ×√２）−ｄ」という設定を行えば、もし個々の母点Ｍが、図５に示す例のように、一定ピッチｄで規則的に配置されている状態であったとすると、各母点Ｍが、１本の「仮想ばね」から受ける力は、すべてＦ＝−ｋｄになり、ばね定数ｋが等しければ、どのばねからも同じ力を受けることになる。したがって、このような設定を行えば、一定ピッチｄで規則的に配置されている母点Ｍに対して作用する力の大きさは、専ら、ばね定数ｋに依存して定まることになる。そして、ここで述べる実施形態の場合、各仮想ばねＳのばね定数ｋは、スカラー場定義部４０で定義されたスカラー場（すなわち、図７に示すような参照画像Ｒの画素値分布）に基づいて決められるので、スカラー場定義部４０においてどのようなスカラー場を定義するかによって、各母点Ｍに対して作用する力を制御することが可能になる。
【００８２】
既に述べたとおり、仮想ばね定義部３０の機能は、個々の連結線Ｃ上に、スカラー場上の当該連結線Ｃの位置に対応する位置に定義されたスカラー値に応じた物理的特性をもった仮想ばねを定義することにある。ここで仮想ばねの物理的特性とは、上述したとおり、自然長Ｌとばね定数ｋであるから、スカラー値に応じて定める物理的特性として、ばね定数ｋを用いる代わりに、自然長Ｌを用いるようにしてもよい。
【００８３】
たとえば、すべての仮想ばねＳのばね定数ｋを一定値に固定し、「縦横連結ばね」の基準自然長Ｌstd＝０、「斜め連結ばね」の基準自然長Ｌstd＝「（ｄ×√２）−ｄ」と設定した上で、対応位置に定義されたスカラー値に応じて、これら基準自然長Ｌstdを増減するようにして自然長Ｌを定め（負の値になってもかまわない）、個々の仮想ばねＳの定義を行うようにしてもよい。このような実施形態を採った場合、一定ピッチｄで規則的に配置されている母点Ｍに対して作用する力の大きさは、専ら、各仮想ばねの自然長Ｌに依存して定まることになる。したがって、この自然長Ｌを、スカラー場定義部４０で定義されたスカラー場に基づいて決めるようにすれば、スカラー場定義部４０においてどのようなスカラー場を定義するかによって、各母点Ｍに対して作用する力を制御することが可能になる。
【００８４】
もちろん、スカラー値に応じて定める物理的特性として、ばね定数ｋと自然長Ｌの双方を用いるようにしてもかまわない。この場合、ばね定数ｋと自然長Ｌの双方が、スカラー場定義部４０で定義されたスカラー場に基づいて決定されることになり、やはりスカラー場定義部４０においてどのようなスカラー場を定義するかによって、各母点Ｍに対して作用する力を制御することが可能になる。要するに、仮想ばね定義部３０は、ばね定数ｋおよび自然長Ｌを物理的特性として有する仮想ばねを定義し、ばね定数ｋおよび自然長Ｌのいずれか一方もしくは双方を、当該仮想ばねが配置される連結線の位置に対応する位置に定義されたスカラー値に基づいて決定すればよい。
【００８５】
＜＜＜  §５．運動シミュレーション  ＞＞＞
ここでは、図３に示す装置におけるシミュレーション実行部５０および変位母点決定部６０の機能について詳述する。シミュレーション実行部５０の基本機能は、母点定義部１０で定義された個々の母点Ｍを質点とみなし、連結線定義部２０で定義された各連結線Ｃの両端に位置する一対の母点が、仮想ばね定義部３０によって当該連結線Ｃ上に定義された仮想ばねＳで連結されているとみなしたときに、各仮想ばねＳによって作用する力に基く各質点の運動シミュレーションを実行することであり、変位母点決定部６０の基本機能は、この運動シミュレーションによって運動中の母点Ｍの所定時点における位置を変位母点の位置と決定することである。
【００８６】
そこで、まず、ある１つの母点に着目し、当該着目母点に作用する力を求める手法を説明する。図１２は、図８に示す特定の着目母点Ｍ０に対して、その八方に定義された仮想ばねＳ１〜Ｓ８によって加わる力を示す平面図である。図示のとおり、着目母点Ｍ０に対しては、仮想ばねＳ１によって、隣接母点Ｍ１の方向に向かう力Ｆ１が作用する。もちろん、この場合、隣接母点Ｍ１に着目すれば、仮想ばねＳ１によって、母点Ｍ０の方向に向かう力（図のＦ１と逆方向の力）が作用していることになる。
【００８７】
図１１(b) に示すとおり、１本の仮想ばねＳに起因して、当該仮想ばねＳの両端に位置する２つの母点Ｍ１０，Ｍ２０間に働く力Ｆの大きさは、両母点間の距離をＤとし、当該両母点間に連結される仮想ばねＳの自然長をＬとすれば、その差（Ｄ−Ｌ）を用いて、Ｆ＝−ｋ（Ｄ−Ｌ）なる式で求められる。ここで行う運動シミュレーションは、図１１(b) に示すような系において、両母点Ｍ１０，Ｍ２０間に仮想ばねＳに起因したＦ＝−ｋ（Ｄ−Ｌ）なる力が作用するものとし、更に、各母点の質量をｍとして、ニュートンの運動方程式「Ｆ＝ｍα」に基づいて各母点の運動加速度αを求め、各母点の運動をシミュレートするものである。
【００８８】
実際には、図１２に示すように、ある１つの着目母点Ｍ０に対しては、図示のとおり、その八方向に向かって、力Ｆ１〜Ｆ８が作用することになるので、着目母点Ｍ０に関する力学的な運動シミュレーションを行う上では、これら八方向への力Ｆ１〜Ｆ８の合力を求めればよい。図１３は、図３に示す装置のシミュレーション実行部５０において、着目母点Ｍ０に対して行われる動作シミュレーションの様子を示す平面図である。
【００８９】
シミュレーションの開始時点ｔ０において、各母点はすべて静止状態であるものとすると、図１３(a) に示すように、静止状態の着目母点Ｍ０（質量ｍをもった質点）に対して、ΣＦなる力が作用し、その結果、ΣＦ＝ｍαで示される加速度αが作用することになる。ここで、ΣＦは、図１２に示す八方向への力Ｆ１〜Ｆ８の合力であり、力Ｆ１〜Ｆ８をベクトルとして表した場合、これら８個のベクトルのベクトル和として求まる値である。
【００９０】
ここで、たとえば、全母点Ｍの質量ｍに一定値（たとえば、ｍ＝１）を与えるようにすれば、結局、加速度αの大きさは合力ΣＦの大きさにのみ依存して定まる量になり、加速度αの向きは合力ΣＦの向きに一致する。そこで、あるステップ時間δを設定し、シミュレーションの開始時点ｔ０から時間δだけ経過した時点ｔ１までの着目母点Ｍ０の運動をシミュレートする。すなわち、時点ｔ０における着目母点Ｍ０の位置、速度（＝０）と、その時点に加わった合力ΣＦに基づいて生じる加速度αを用いて、時間δだけ経過した時点ｔ１における着目母点Ｍ０の新たな位置と速度Ｖとを求める。同様の処理を、すべての母点をそれぞれ着目母点として実行する。
【００９１】
かくして、時点ｔ０から時点ｔ１に至る期間の運動シミュレーションが完了し、すべての母点の時点ｔ１における新たな位置が決定する。なお、シミュレーションの開始時点ｔ０では、すべての母点は静止状態であり、その速度成分はいずれも０であったが、時点ｔ１では、すべての母点はそれぞれ固有の速度成分を有していることになる。したがって、シミュレーション実行部５０は、時点ｔ１において、各母点の位置（たとえば、ＸＹ座標値）と速度（たとえば、速度ベクトルＶ）とをデータとして格納していることになる。
【００９２】
もちろん、時点ｔ０〜ｔ１に至る期間における各母点の運動を厳密に捉えると、隣接する母点の位置も連続的に変化し、その結果、着目母点Ｍ０に作用する合力ΣＦの大きさや向きも連続的に変化してゆくので、着目母点Ｍ０に加わる加速度αも連続的に変化してゆくことになる。しかしながら、コンピュータ上のシミュレーションとしては、個々の母点の加速度αを、時点ｔ０において算出し、その後、時点ｔ１に至るまでは、当該加速度αが一定であるものとして、上述のような運動シミュレーションを行っても問題はない。
【００９３】
続いて、時点ｔ１から時間δだけ経過した時点ｔ２までの各母点の運動を、同様の方法でシミュレートする。すなわち、図１３(b) に示すように、時点ｔ１における着目母点Ｍ０の位置、速度Ｖと、当該時点ｔ１に加わった合力ΣＦ（これは、時点ｔ１において、各母点Ｍが新たな位置まで移動し、その結果、各仮想ばねＳに伸縮が生じ、仮想ばねに起因する力に変動が生じたものとして、新たに算出された合力である）に基づいて生じる加速度α（時点ｔ１における加速度）を用いて、時間δだけ経過した時点ｔ２における着目母点Ｍ０の新たな位置と速度Ｖとを求めるのである。同様の処理を、すべての母点をそれぞれ着目母点として実行する。
【００９４】
かくして、時点ｔ１から時点ｔ２に至る期間の運動シミュレーションが完了し、すべての母点の時点ｔ２における新たな位置が決定する。この時点ｔ２でも、すべての母点はそれぞれ固有の速度成分を有していることになる。したがって、シミュレーション実行部５０は、今度は、時点ｔ２における、各母点の位置と速度とをデータとして格納することになる。
【００９５】
このようにして、ステップ時間δごとに、各母点の新たな位置および速度を求めてゆく処理を、所定回数だけ繰り返し実行してゆけば、各母点の位置は徐々に変化してゆくことになる。こうして、予め設定された所定回数だけ上記処理を繰り返したら、その時点における最終的な母点の位置を、変位母点の位置と決定する。
【００９６】
なお、上記シミュレーションに相当する現実の物理的現象を考えると、各ばねは所定のダンパー係数をもった実在の物体であり、運動には様々な抵抗が加わるため、ばねで連結された多数の母点からなる系全体は、やがて安定した状態で収束することになり、最終的には各母点は静止状態となる。ところが、上記シミュレーションでは、各仮想ばねＳは、各時点での実際の長さと自然長との差に基づいて力を発生させる概念的な要素であり、個々の母点Ｍは、体積をもたない質点であるため、ステップ時間δごとの運動シミュレーションを繰り返したとしても、運動が収束することはない。したがって、予め所定の実行回数を定めておき、当該実行回数分だけ繰り返し処理を完了したら、その時点における最終的な母点の位置を、変位母点の位置とするようにしている。
【００９７】
もちろん、各仮想ばねにそれぞれダンパー係数を定義し、各母点を所定の形状および体積をもった物体として取り扱い、様々な抵抗要素を考慮したシミュレーションを実行し、全体の系がある程度収束した状態になるまで、シミュレーションを続行するようにしてもかまわない。ただ、実用上は、上述した例のように、単純な系でシミュレーションを実行し、所定の実行回数分だけ繰り返した時点で、その時点における最終的な母点の位置を、変位母点の位置とする方が、演算負担を軽くする上で好ましい。
【００９８】
要するに、上述した実施形態では、シミュレーション実行部５０は、第ｋ番目の時点ｔｋにおいて、各母点のそれぞれについて、その時点の仮想ばねに起因して作用する力の合力を求め、第ｋ番目の時点ｔｋから第（ｋ＋１）番目の時点ｔ（ｋ＋１）に至るまでの時間δの間における各母点の合力に基づく移動経路を求め、時点ｔ（ｋ＋１）における各母点の位置を決定する処理を、ｋ＝０からｋを１ずつ増加させながら繰り返し実行することになる。そして、変位母点決定部６０は、ｋの値が所定値Ｋに達した時点における各母点の位置を変位母点の位置と決定することになる。
【００９９】
図１４は、こうして変位母点決定部６０によって決定された変位母点の位置を示す平面図である。この図１４に示す例は、図５に示す規則的配置がなされた母点Ｍについて、上述した運動シミュレーションを行うことにより得られた最終結果を示している。なお、この例では、作業領域Ａの輪郭近傍に位置する母点を不動の固定点として、運動シミュレーションを実行しているため、作業領域Ａの輪郭上の１６個の母点の位置は、初期設定の状態のままとなっており、内部の９個の母点の位置のみが変動している。
【０１００】
もちろん、全母点を変動の対象としたシミュレーションを行ってもかまわないが、最終的に利用する皺画像は、実用上、その輪郭線が矩形であることが好ましい。そこで、ここで述べる実施形態の場合、作業領域Ａの輪郭近傍に位置する母点を不動の固定点として取り扱い、最終的に得られる皺画像の輪郭線が矩形となるよう配慮している。このように、作業領域Ａの輪郭近傍に位置する母点を固定点とする取り扱いを行うと、シミュレーションの途中で、母点が作業領域Ａを超えて発散してしまう現象を防止するメリットも得られる。
【０１０１】
なお、前述したとおり、各母点Ｍ間を連結する連結線Ｃは、２つの母点の連結関係を示す概念的な線であるから、上述した運動シミュレーションにより、個々の母点の位置が変動したとしても、各連結線Ｃそれ自体は何ら変わることはない。すなわち、図１４に示す各連結線Ｃは、個々の母点Ｍの移動に伴い伸縮したことになる。
【０１０２】
＜＜＜  §６．運動シミュレーションの変形例  ＞＞＞
ここでは、§５で述べた運動シミュレーションの変形例をいくつか述べておく。
【０１０３】
（１）  母点の質量ｍをスカラー値に応じて決定する変形例
これまで述べた実施形態では、仮想ばね定義部３０によって、個々の連結線Ｃ上に、スカラー場定義部４０で定義されたスカラー場を利用して、対応位置のスカラー値に応じた物理的特性をもった仮想ばねを定義していた。この場合、スカラー値に応じて定める物理的特性としては、ばね定数ｋを用いてもよいし、自然長Ｌを用いてもよいが、いずれにせよ、スカラー値に応じて特性が定められる対象は、「仮想ばね」であった。
【０１０４】
すなわち、§５で述べた運動シミュレーションでは、すべての母点に同一の質量ｍを与え、ある母点に作用する合力ΣＦに基づいて、ΣＦ＝ｍαなる式を用いて、当該母点に作用する加速度αを求め、時間δ経過後の当該母点の位置を求めていた。この場合、仮想ばねの物理的特性（ばね定数ｋもしくは自然長Ｌ）が、スカラー値に応じて定められるため、各母点に作用する力Ｆがスカラー場の影響を受け、最終的な変位母点の位置が、スカラー場の影響を受ける結果となる。
【０１０５】
これに対して、ここで述べる変形例では、スカラー値に応じて定められる対象を、「仮想ばねの物理的特性」ではなく、「母点の質量ｍ」にするのである。この場合、仮想ばね定義部３０が定義する仮想ばねの物理的特性は、スカラー場とは全く無関係に定めてかまわない。したがって、図３において、仮想ばね定義部３０は、スカラー場定義部４０内に定義されたスカラー場を参照する必要はない（図３のブロック４０から３０へ向かう矢印は不要になる）。
【０１０６】
ただ、１つの母点に作用するすべての力が釣り合ってしまうと（たとえば、図１２において、力Ｆ１〜Ｆ８の大きさがすべて等しくなってしまうと）、合力ΣＦが０となり、母点が移動しなくなってしまう。したがって、仮想ばね定義部３０は、各母点に作用する力が釣り合わないように、各仮想ばねの物理的特性（ばね定数ｋや自然長Ｌ）が互いに適宜異なるように設定するようにする（たとえば、各ばね定数ｋを、乱数を用いて決定すればよい）。
【０１０７】
一方、シミュレーション実行部５０は、運動シミュレーションを行う上で、個々の母点をスカラー場定義部４０で定義されたスカラー場上の、当該母点のシミュレーション開始時の位置に対応する位置に定義されたスカラー値に応じた質量をもった質点とみなし、各連結線の両端に位置する一対の母点が当該連結線上に定義された仮想ばねで連結されているとみなしたときに、各仮想ばねによって作用する力に基く各質点の運動シミュレーションを実行する。図３に示すブロック４０からブロック５０に向かう破線の矢印は、ここで述べる変形例を実施する場合に必要な情報の流れを示しており、スカラー場定義部４０で定義されたスカラー場の情報が、シミュレーション実行部５０におけるシミュレーション実行時に参照されることを示している。
【０１０８】
具体的には、図９に示すように、作業領域Ａ上に参照画像Ｒを重ね、ある母点の質量ｍを、当該母点上に位置する画素の画素値、または当該母点の近傍に位置する複数の画素の画素値の和もしくは平均に基づいて決定するようにすればよい。なお、各母点の質量ｍは、シミュレーション開始時の当該母点の位置に基づいて一度決定した後は不変とする取り扱いをしてもよいし、その代わりに、ステップ時間δごとの新たな母点位置に基づいてその都度更新するようにしてもかまわない。現実の物理現象としては、時々刻々と質量が変化してゆく物体は考えにくいが、本発明でのシミュレーションでは、そのような取り扱いをしても問題はない。
【０１０９】
このように、各母点の質量ｍを、スカラー場定義部４０で定義されたスカラー場（すなわち、図７に示すような参照画像Ｒの画素値分布）に基づいて決めるようにすれば、スカラー場定義部４０においてどのようなスカラー場を定義するかによって、各母点Ｍの運動を制御することが可能になり、§５で述べた実施形態と同様に、スカラー場の影響を受けた変位母点の分布を得ることができるようになる。
【０１１０】
もちろん、必要があれば、「仮想ばねの物理的特性（ばね定数ｋや自然長Ｌ）」と「母点の質量ｍ」との双方を、スカラー場定義部４０で定義されたスカラー場に基づいて決めるようにしてもよい。
【０１１１】
（２）  ばねの破断を考慮する変形例
これまで述べた実施形態では、仮想ばねＳは、どのような大きな力が加わろうと、常に、加えられた力に応じた伸びを生じる仮想のばねであった。これに対して、ここで述べる変形例では、各仮想ばねＳを、所定のしきい値Ｆth（各ばねごとに異なる値を設定してもよい）を超える力が作用した場合には破断するばねとして取り扱うのである。すなわち、シミュレーション実行部５０は、特定の両母点間に作用する力Ｆが所定のしきい値Ｆthを超えた場合に、当該特定の両母点間を連結する仮想ばねＳが破断したものと判断し、以後、当該特定の両母点間には仮想ばねが存在しないものとして、運動シミュレーションを実行するようにすればよい。なお、前述したとおり、両母点間を連結する連結線Ｃは、２つの母点の連結関係を示す概念的な線であるから、仮想ばねＳが破断したとしても、連結線Ｃ自体はそのまま残ることになる。
【０１１２】
このように、ある特定の母点間を連結する仮想ばねが破断し、当該両母点間には仮想ばねが存在しないものとすると、当該両母点間の距離は、運動シミュレーションを実行してゆく過程で増大してゆく可能性が高い。したがって、運動シミュレーションが完了し、変位母点決定部６０によって、各母点の最終的な位置が決定した段階において、仮想ばねが破断した両母点間の距離、すなわち、当該両母点を連結する連結線Ｃの長さは、非常に長くなっている可能性が高い。後の§７のプロセスで述べるように、このように非常に長い連結線の存在は、非常に長い皮溝部を形成する要因になり、図１に示すような血筋を形成する要因になる。
【０１１３】
＜＜＜  §７．皺画像の生成  ＞＞＞
ここでは、図３に示す装置における皺構成線分配置部７０および皺画像生成部８０の機能について詳述する。まず、皺構成線分配置部７０の基本機能は、図１に示す皮シボパターンにおける皮溝部の芯となるべき皺構成線分を作業領域Ａ上に配置する処理を行うことにある。この皺構成線分は、運動シミュレーションによって変位した後の個々の母点間を連結する連結線Ｃ（たとえば、図１４に示す各連結線）を利用して作成される。
【０１１４】
皺構成線分配置部７０による皺構成線分の配置アルゴリズムは、各変位母点Ｍについての連結線Ｃ上にその中点Ｇを求め、当該連結線Ｃに対して直交し、当該連結線Ｃと同じ長さをもつ皺構成線分Ｔを、当該皺構成線分Ｔの中点が当該連結線Ｃの中点Ｇに重なるように配置するというものである。図１５は、このような皺構成線分Ｔの配置原理を示す平面図である。
【０１１５】
いま、図１５に示すように、２つの母点Ｍ１０，Ｍ２０（いずれも、§５，§６で述べた運動シミュレーションによって変位した後のもの）の間に連結線Ｃが存在するものとし、この連結線Ｃに基づいて、皺構成線分Ｔを作成して配置することを考える。この場合、まず、連結線Ｃ上にその中点Ｇを求める。そして、この連結線Ｃと同じ長さをもつ皺構成線分Ｔを作成し、これを連結線Ｃに対して直交するように配置する。しかも、この皺構成線分Ｔの中点が連結線Ｃの中点Ｇに重なるように配置するのである。いわば、皺構成線分Ｔは、２点Ｍ１０，Ｍ２０を結ぶ連結線Ｃの垂直二等分線に相当する。
【０１１６】
皺構成線分配置部７０は、変位母点決定部６０によって位置が決定された変位母点に関するすべての連結線Ｃについて、それぞれ上述のアルゴリズムに基づいて、皺構成線分Ｔを作成して配置する処理を行う。図１６は、図１４に示す変位母点に基づいて配置された皺構成線分Ｔを示す平面図である。すなわち、図１６に示されている各皺構成線分Ｔは、図１４に示されている連結線Ｃのいずれか１つに対応している。たとえば、皺構成線分Ｔ１は連結線Ｃ１に対応し、皺構成線分Ｔ２は連結線Ｃ２に対応し、皺構成線分Ｔ３は連結線Ｃ３に対応している。そして、ある皺構成線分Ｔと、これに対応する連結線Ｃとは、互いに同一の長さを有し、互いにその中点において直交する関係にある。図１７は、図１６から変位母点Ｍを削除した平面図であり、作業領域Ａ内に、皺構成線分Ｔのみが配置された状態を示している。なお、皺構成線分Ｔによっては、その一部分が作業領域Ａの輪郭線を食み出して外側へと伸びるものもあるが、ここに示す実施形態は、皺構成線分Ｔのうち、作業領域Ａの輪郭線を食み出た部分はカットしている。
【０１１７】
前述した§６（２）では、「ばねの破断を考慮する変形例」を説明した。この変形例の場合、ある特定の母点間を連結する仮想ばねが破断すると、当該両母点を連結する連結線Ｃの長さは、非常に長くなっている可能性が高い。そのため、当該連結線Ｃに対応する皺構成線分Ｔも非常に長いものになる。そして、このような長い皺構成線分Ｔは、図１に示すような血筋を形成する要因になるのである。したがって、部分的に血筋を含む皮シボパターンを生成するためには、この「ばねの破断を考慮する変形例」を実施し、部分的にばねの破断が生じるような条件設定を行い、運動シミュレーションを行うようにするとよい。
【０１１８】
さて、こうして配置された多数の皺構成線分Ｔは、皮シボパターンにおける皮溝部の芯となるべき線分であり、この芯の部分に肉付けを行うことにより、皮溝部（皺部分）の形成が行われる。図３に示す皺画像生成部８０は、図１７に示すような皺構成線分Ｔの集合体に基づいて、皺画像を生成する機能を果たす。
【０１１９】
各皺構成線分Ｔに対して肉付けして、皮溝部を形成するには、各皺構成線分Ｔに幅をもたせる処理を行い、皺構成帯状領域Ｖを形成し、各皺構成帯状領域Ｖの外部の画素には背景部分を示す画素値を与え、各皺構成帯状領域Ｖの内部の画素には皺部分を示す画素値を与える処理を行えばよい。このような画素値をもった画素の集合体が、本発明に係る装置を用いて生成する目的となる皺画像（皮シボパターン）である。
【０１２０】
図１８(a) は、１本の皺構成線分Ｔに幅Ｗをもたせる処理を行い、皺構成帯状領域Ｖを形成する最も簡単な方法の一例を示す平面図である。この方法では、皺構成線分Ｔを中心線として、幅Ｗをもった矩形図形を生成し、当該矩形図形の内部（図にハッチングを施す部分）の画素には皺部分を示す画素値を与え、外部の画素には背景部分を示す画素値を与えればよい。幅Ｗは任意に設定することができ、たとえば、１画素分の幅に設定すれば、非常に細い皺構成帯状領域Ｖが形成される。
【０１２１】
なお、幅Ｗは、個々の皺構成線分Ｔごとにバラバラな値に設定するのが好ましい。そのためには、幅の範囲「Ｗ１〜Ｗ２」をパラメータとして設定し、個々の皺構成線分Ｔについて、それぞれＷ１〜Ｗ２内のランダム値Ｗを定め、それぞれランダムな幅Ｗをもった皺構成帯状領域Ｖを形成するようにすればよい。
【０１２２】
図１８(a) に示す例では、皺構成帯状領域Ｖの内部か外部かを識別するための２通りの画素値しか用いられていないため、生成される皺画像は二値画像になる。そこで、階調をもった皺画像を生成するためには、内側から外側に向かって画素値がなだらかに変化するように、各皺構成帯状領域Ｖの内部の画素に画素値を与えるようにすればよい。図１８(b) は、皺構成線分Ｔの位置を中心として、上下方向になだらかに変化するような画素値を定義した例である。いわば、上下方向にグラデーションがかかった皺構成帯状領域Ｖが形成されている。
【０１２３】
このように階調をもった二次元の皺画像は、階調値を高さ（深さ）を示す情報として取り扱えば、実質的に三次元の形状を表現した三次元画像としての機能を果たすことができる。一般に、皮シボパターンは、凹凸をもったエンボス加工などにも利用されるが、階調をもった二次元皺画像を生成すれば、凹凸をもったエンボス加工にもそのまま利用することが可能である。たとえば、図１８(b) の上下方向に関する画素値の分布が、ガウス分布となるように設定しておけば、凹凸をもったエンボス加工により、深さの分布がガウス分布となるような皮溝部を形成することができるようになり、非常に自然な深さ分布を得ることができる。
【０１２４】
もちろん、図１８(a) ，(b) に示す例は、皺構成線分Ｔに対して肉付けして、皺構成帯状領域Ｖ（皮溝部などの皺の部分となる領域）を形成する単純な方法の一例を示すものであり、実際には、この他にも様々な方法を利用して、皺構成帯状領域Ｖを形成することができる。また、１本の線分を加工して、帯状の閉領域を形成するアルゴリズムも種々のものが知られており、このようなアルゴリズムを利用した画像処理ソフトウエアも種々の製品が市販されている。
【０１２５】
図１９は、このような市販の画像処理ソフトウエアを利用して、図１７に示す各皺構成線分Ｔに肉付け処理を施し、皺構成帯状領域Ｖを形成した一例を示す平面図である。この画像処理ソフトウエアは、１本の線からなる線画に対して毛筆調の肉付けを施すためのアルゴリズムを用いた処理を行う機能を有しており、各皺構成帯状領域Ｖは、それぞれ毛筆で描いた線のような形態をなしている。もちろん、皺構成線分Ｔに対して肉付けして、皺構成帯状領域Ｖを形成する手法としては、この他にも種々のアルゴリズムを利用することが可能である。要するに、皺画像生成部は、各皺構成線分Ｔに対して肉付けして、帯状の領域Ｖを形成する処理を行うことができれば、どのようなアルゴリズムで肉付けを行ってもかまわない。
【０１２６】
この図１９に示す例は、もともと図４に示すような２５個の母点配置から出発して得られた皺画像であるが、わずか２５個（そのうちの１６個は固定点）の母点を用いた単純な実施例であるにもかかわらず、ある程度自然な風合いをもった皺画像となっており、皮シボパターンとして十分に利用可能である。
【０１２７】
図２０は、天然皮革の表面に現れる皮シボパターンを写真撮影し、これに画像処理（二値画像に変換する処理）を施して得られた画像を示す拡大平面図である。別言すれば、天然の皮革に表現された皺画像ということになる。これに対して、図２１は、図３に示す本発明に係る皺画像生成装置を用いて人為的に作成した皮シボパターンに画像処理（二値画像に変換する処理）を施して得られた画像を示す拡大平面図である。この図２１に示す皺画像の生成には、１００行１００列に配置された１００００個の母点が用いられている。図２０と図２１とを比較すれば、本発明に係る装置によって生成された皺画像は、天然皮革の革シボパターンに比べても遜色のない意匠性を有していることがわかる。しかも、用いるスカラー場（参照画像Ｒ）を変更することにより、異なる革シボパターンを容易に生成することができ、非常に自由度の高い皺画像の生成が可能である。
【０１２８】
＜＜＜  §８．好ましいスカラー場の例  ＞＞＞
最後に、図３に示すスカラー場定義部４０におけるスカラー場の定義方法について、いくつかの実例を示しておく。既に述べたとおり、スカラー場定義部４０で定義されるスカラー場とは、実質的に、作業領域Ａに対応する参照領域と同じサイズの参照画像Ｒと言うべきものであり、スカラー場定義部４０は、このような参照画像Ｒを定義する機能を有していればよい。
【０１２９】
もっとも、参照画像Ｒは、仮想ばねの物理的特性値もしくは母点に与えるべき質量を決定するために用いられるものであり、所定の位置に所定の値（スカラー値）が定義されていれば、本発明に係る装置に利用可能である。したがって、原理的には、世の中に存在するあらゆる画像が、本発明において参照画像Ｒとして利用可能である。ただ、自然な風合いをもった皺画像を得るためには、より適切な画像を用いるのが好ましい。
【０１３０】
本発明に係る装置に参照画像Ｒとして用いるのに適した第１の候補画像は、天然の皮革などの表面画像である。そこで、スカラー場定義部４０に、実在の物体表面の画像を取り込む画像取込装置を設けておき、取り込んだ画像もしくはこれを加工して得られた画像を参照画像Ｒとして利用するようにするとよい。たとえば、牛や豚などの天然皮革の画像を画像取込装置を用いて取り込み、トリミング処理やフィルタリング処理などを施したものを参照画像Ｒとして利用すれば、自然な風合いをもった皺画像を得るために好ましい。取り込んだ画像それ自身は、皺画像として利用されるものではなく、あくまでもスカラー場を与える参照画像として用いられるものであるから、照明むらや影が混在していても問題はない。
【０１３１】
もちろん、スカラー場定義部４０に、人為的にスカラー場を定義する機能（すなわち、
人為的に何らかの参照画像Ｒを生成する機能）をもたせておくこともできる。上述したとおり、原理的には、どのような画像でも、本発明において参照画像Ｒとして利用することが可能である。しかしながら、本願発明者が、種々の画像を参照画像Ｒとして用いて試したところ、自然な風合いをもった革シボパターンを生成する上では、線状模様をもった参照画像Ｒを用いるのが好ましいことが判明した。特に、ある共通した方向を向いた複数の帯状領域から構成される画像を参照画像Ｒとして用いると、革シボパターンとして好ましい皺画像が得られることがわかった。
【０１３２】
そこで、以下、このような線状模様をもった参照画像Ｒを作成する機能をスカラー場定義部４０にもたせた実施形態を述べる。図２２は、このような実施形態に係るスカラー場定義部４０の基本構成を示すブロック図である。図示のとおり、このスカラー場定義部４０は、パラメータ設定部４１、核点設定部４２、参照帯状領域設定部４３、参照画像生成部４４によって構成されている。もっとも、これらの各構成要素も、実際には、コンピュータに所定のプログラムを組み込むことによって構成される。
【０１３３】
パラメータ設定部４１は、空間周波数を示すパラメータＦを設定するとともに、角度の範囲「θ１〜θ２」、幅の範囲「Ｕ１〜Ｕ２」、長さの範囲「Ｈ１〜Ｈ２」をパラメータとして設定する機能を有する。これらのパラメータは予め既存値として記憶させておいてもよいし、オペレータの入力操作に基づいて、適宜設定できるようにしておいてもよい。ここで、パラメータＦは、生成される線状模様の空間周波数（配置ピッチに相当）を決めるためのものである。また、角度の範囲「θ１〜θ２」、幅の範囲「Ｕ１〜Ｕ２」、長さの範囲「Ｈ１〜Ｈ２」は、生成される線状模様の向き、幅、長さの範囲を決めるためのものである。上述したとおり、自然な風合いをもった皺画像を得るための参照画像Ｒを生成する上では、ある共通した方向を向いた複数の帯状領域を配置するのが効果的であるので、角度の範囲「θ１〜θ２」は、あまり大きくしないようにするのが好ましい（具体的には、θ１−θ２＜３０°程度がよい）。
【０１３４】
核点設定部４２は、参照領域内にＸＹ二次元座標系を定義し、多数の核点Ｐ（Ｘ，Ｙ）をプロットする処理を行う。ここで、核点Ｐとは、線状模様を形成するための個々の帯状領域の核となる点である。
【０１３５】
図２３は、この核点設定部４２の具体的な処理手順を示す流れ図である。まず、ステップＳ１において、変数ｉおよびｊを初期値０に設定する。変数ｉはＸ軸方向に関する配置位置を示すパラメータであり、変数ｊはＹ軸方向に関する配置位置を示すパラメータである。続くステップＳ２では、変数ｊを１だけ増加させ、更にステップＳ３では、変数ｉを１だけ増加させる。次のステップＳ４では、核点Ｐを配置するためのＹ座標値が、Ｙ＝（Ｆ×ｊ）＋（Ｆ×Rnd）なる演算式で求められ、ステップＳ５では、核点Ｐを配置するためのＸ座標値が、Ｘ＝（Ｆ×ｉ）＋（Ｆ×Rnd）なる演算式で求められる。ここで、Ｆは、パラメータ設定部４１で設定された空間周波数を示すパラメータであり、Rndは、−０．５〜＋０．５の範囲内の乱数値である。たとえば、Ｆ＝１０に設定されていた場合、ｉ＝１，ｊ＝１では、Ｘは５〜１５の範囲内の値をとり、Ｙも５〜１５の範囲内の値をとる。
【０１３６】
こうして、座標値Ｘ，Ｙが決まると、ステップＳ６，Ｓ７を経て、ステップＳ８へと進み、核点Ｐ（Ｘ，Ｙ）のプロットが行われる。そして、ステップＳ３へと戻り、変数ｉの値が１だけ更新され、ステップＳ４，Ｓ５へと進む。今度は、ｉ＝２，ｊ＝１であるから、Ｘは１５〜２５の範囲内の値をとり、Ｙは５〜１５の範囲内の値をとるので、ステップＳ８では、この範囲内の座標位置に次の核点Ｐ（Ｘ，Ｙ）のプロットが行われる。再びステップＳ３へと戻り、変数ｉの値が１だけ更新され、ステップＳ４，Ｓ５へと進む。今度は、ｉ＝３，ｊ＝１であるから、Ｘは２５〜３５の範囲内の値をとり、Ｙは５〜１５の範囲内の値をとるので、ステップＳ８では、この範囲内の座標位置に次の核点Ｐ（Ｘ，Ｙ）のプロットが行われる。
【０１３７】
やがて、Ｘの値が最大値Ｘmaxを超えると（この時点でｉ＝最大値Ｉとなっている）、ステップＳ７からステップＳ９へと分岐し、変数ｉ＝０にリセットされ、ステップＳ２からの処理が繰り返される。今回は、ｊ＝２となるので、Ｙは１５〜２５の範囲内の値をとることになる。こうして、最後に、Ｙの値が最大値Ｙmaxを超えると（この時点でｊ＝最大値Ｊとなっている）、ステップＳ６を経て処理は終了する。図２４は、このような処理によって、ＸＹ二次元座標系上に複数の核点Ｐがプロットされた状態を示す平面図である。上述の手順により核点の設定を行えば、ほぼ空間周波数Ｆで示されるピッチで、ほぼ一様に分布するように、ランダムに核点Ｐを配置することが可能になる。
【０１３８】
このように、図２２に示す核点設定部４２の基本機能は、−０．５〜＋０．５の範囲内のランダム値Rndと整数ｉおよびｊを用いて、Ｘ＝（Ｆ×ｉ）＋（Ｆ×Rnd）およびＹ＝（Ｆ×ｊ）＋（Ｆ×Rnd）なる演算によって求まる座標値（Ｘ，Ｙ）で示される位置に核点Ｐ（Ｘ，Ｙ）をプロットする処理を、ｉの値を１から所定の最大値Ｉまで１ずつ変え、ｊの値を１から所定の最大値Ｊまで１ずつ変えた各組み合わせについて実行し、合計（Ｉ×Ｊ）個の核点Ｐを設定することである。
【０１３９】
参照帯状領域設定部４３は、こうして設定された個々の核点Ｐについて、それぞれ参照帯状領域Ｂを設定する処理を行う。具体的には、個々の核点Ｐについて、θ１〜θ２内のランダム値θ、Ｕ１〜Ｕ２内のランダム値Ｕ、Ｈ１〜Ｈ２内のランダム値Ｈを定め、核点Ｐを中心点として、基準となる座標軸（たとえば、Ｘ軸）に対して角度θをなす方向を向き、幅Ｕ，長さＨをもった参照帯状領域Ｂを定める処理を行い、合計（Ｉ×Ｊ）個の参照帯状領域を設定する処理が行われる。
【０１４０】
図２５は、図２４に示す個々の核点Ｐについて、それぞれ参照帯状領域Ｂを設定した状態を示す平面図である。図に描かれた個々の矩形が、各核点Ｐを中心として設定された参照帯状領域Ｂである。この参照帯状領域Ｂは、Ｘ軸に対して角度θをなす方向を向き、幅Ｕ、長さＨをもった矩形であるが、角度θはθ１〜θ２内のランダム値として定められ、幅ＵはＵ１〜Ｕ２内のランダム値として定められ、長さＨはＨ１〜Ｈ２内のランダム値として定められるため、それぞれバラバラな形状を有し、バラバラな向きを向いている。ただ、θ１−θ２＜３０°程度に設定しておけば、向きに関するバラつきは、角度３０°程度に抑えられ、ある程度共通した方向を向いた複数の帯状領域が配置されることになる。
【０１４１】
参照画像生成部４４は、各参照帯状領域Ｂの外部の画素には背景部分を示す画素値（たとえば、”０”）を与え、各参照帯状領域Ｂの内部の画素には線状模様部分を示す画素値（たとえば、”１”）を与えることにより、線状模様をもった参照画像を生成する処理を行う。結局、図２５に示すような参照画像Ｒが二値画像として生成されることになる。なお、二値画像ではなく、階調画像を生成する場合には、参照画像生成部４４が、内側から外側に向かって画素値がなだらかに変化するように、各参照帯状領域Ｂの内部の画素に画素値を与えるようにすればよい。
【０１４２】
このような線状模様をもった参照画像Ｒを、スカラー場定義部４０内で自動生成し、これをスカラー場として用いるようにすると、革シボパターンとして好ましい皺画像を得ることができる。
【図面の簡単な説明】
【０１４３】
【図１】革シボパターンの基本構成例を示す平面図である。
【図２】天然皮革表面の一般的な構造を示す拡大斜視図である。
【図３】本発明に係る皺画像生成装置の基本構成を示すブロック図である。
【図４】図３に示す装置の母点定義部１０で定義された母点Ｍの一例を示す平面図である。
【図５】図３に示す装置の連結線定義部２０で定義された連結線Ｃの一例を示す平面図である。
【図６】特定の着目母点Ｍ０について定義された８本の連結線Ｃ１〜Ｃ８を示す平面図である。
【図７】図３に示す装置のスカラー場定義部４０で定義されたスカラー場を示すための参照画像Ｒの一例を示す平面図である。
【図８】図３に示す装置の仮想ばね定義部３０において、図６に示す各連結線Ｃ１〜Ｃ８上に定義された仮想ばねＳ１〜Ｓ８を示す平面図である。
【図９】各仮想ばねの物理的特性を決定するために、図５に示す連結線Ｃが定義された作業領域Ａ上に、図７に示す参照画像Ｒを重ね合わせた状態を示す平面図である。
【図１０】作業領域Ａに重ね合わせた参照画像Ｒ内の画素の画素値に基づいて、各仮想ばねの物理的特性を決定する具体的な方法の一例を示す平面図である。
【図１１】図３に示す装置の仮想ばね定義部３０において定義される仮想ばねＳの概念図である。
【図１２】図８に示す特定の着目母点Ｍ０に対して、その八方に定義された仮想ばねＳ１〜Ｓ８によって加わる力を示す平面図である。
【図１３】図３に示す装置のシミュレーション実行部５０において、着目母点Ｍ０に対して行われる動作シミュレーションの様子を示す平面図である。
【図１４】図３に示す装置の変位母点決定部６０によって決定された変位母点の位置を示す平面図である。
【図１５】図３に示す装置の皺構成線分配置部７０による皺構成線分Ｔの配置原理を示す平面図である。
【図１６】図３に示す装置の皺構成線分配置部７０によって、図１４に示す変位母点に基づいて配置された皺構成線分Ｔを示す平面図である。
【図１７】図１６の皺構成線分Ｔのみを示す平面図である。
【図１８】図３に示す装置の皺画像生成部８０による皺画像生成処理の原理を示す平面図である。
【図１９】図３に示す装置の皺画像生成部８０によって、図１７に示す皺構成線分Ｔに基づいて生成された皺画像を示す平面図である。
【図２０】天然皮革の表面に現れる皮シボパターンを写真撮影し、これに画像処理を施して得られた二値画像を示す拡大平面図である。
【図２１】本発明に係る皺画像生成装置を用いて人為的に作成した皮シボパターンに、画像処理を施して得られた二値画像を示す拡大平面図である。
【図２２】図３に示す装置のスカラー場定義部４０の具体的な構成例を示すブロック図である。
【図２３】図２２に示す核点設定部４２の処理手順を示す流れ図である。
【図２４】図２２に示す核点設定部４２によって設定された核点Ｐの一例を示す平面図である。
【図２５】図２２に示す参照帯状領域設定部４３によって設定された参照帯状領域Ｂの一例を示す平面図である。
【符号の説明】
【０１４４】
１０：母点定義部
２０：連結線定義部
３０：仮想ばね定義部
４０：スカラー場定義部
４１：パラメータ設定部
４２：核点設定部
４３：参照帯状領域設定部
４４：参照画像生成部
５０：シミュレーション実行部
６０：変位母点決定部
７０：皺構成線分配置部
８０：皺画像生成部
Ａ：作業領域
ａ：作業領域Ａ／参照領域Ｒの寸法
Ｂ：参照帯状領域
Ｃ，Ｃ１〜Ｃ１２：連結線
Ｄ：母点間距離
ｄ：母点の初期配置の格子間隔
Ｅ：特定領域
Ｆ，Ｆ１〜Ｆ８：仮想ばねによって作用する力
Ｇ：連結線Ｃの中点
Ｈ：参照帯状領域の長さ
ｋ：ばね定数
Ｌ：仮想ばねの自然長
Ｍ，Ｍ１〜Ｍ８，Ｍ１０〜Ｍ２０：母点
Ｍ０：着目母点
ｍ：母点の質量
Ｐ：核点
Ｒ：参照領域（参照画像）
Ｓ，Ｓ１〜Ｓ８：仮想ばね
Ｓ１〜Ｓ９：流れ図の各ステップ
Ｔ，Ｔ１〜Ｔ３：皺構成線分
Ｕ：参照帯状領域の幅
Ｖ：皺構成帯状領域／速度
Ｗ：皺構成帯状領域の幅
α：加速度
θ：参照帯状領域の傾き
【図１】
0004905282.tif 000002
【図２】
0004905282.tif 000003
【図３】
0004905282.tif 000004
【図４】
0004905282.tif 000005
【図５】
0004905282.tif 000006
【図６】
0004905282.tif 000007
【図７】
0004905282.tif 000008
【図８】
0004905282.tif 000009
【図９】
0004905282.tif 000010
【図１０】
0004905282.tif 000011
【図１１】
0004905282.tif 000012
【図１２】
0004905282.tif 000013
【図１３】
0004905282.tif 000014
【図１４】
0004905282.tif 000015
【図１５】
0004905282.tif 000016
【図１６】
0004905282.tif 000017
【図１７】
0004905282.tif 000018
【図１８】
0004905282.tif 000019
【図１９】
0004905282.tif 000020
【図２０】
0004905282.tif 000021
【図２１】
0004905282.tif 000022
【図２２】
0004905282.tif 000023
【図２３】
0004905282.tif 000024
【図２４】
0004905282.tif 000025
【図２５】
0004905282.tif 000026