"ゴルフクラブのヘッドの挙動を計測する方法であって、クラブヘッドのフェース面に少なくとも3個のマークを付するステップと、複数台のカメラにより、移動する上記クラブヘッドの少なくとも2時点において、上記マークの複数の2次元データを得るステップと、上記複数の2次元データから、マークの3次元座標上の位置を特定するステップと、3個以上のマークの3次元座標上の位置データからフェース面を特定するステップと、フェース面とゴルフボールとの位置関係を時系列に更新することにより、フェース面とゴルフボールとの接触時点を特定するステップとを含むクラブヘッドの挙動の計測方法。上記少なくとも3個のマークの位置が等速直線運動をしていると仮定し、この3つの位置によって構成される面の法線ベクトルであって、ゴルフボールの中心を通る法線ベクトルを用いて、この面がゴルフボールと接触する時点を求める請求項1に記載のクラブヘッドの挙動の計測方法。上記フェース面とゴルフボールとの位置関係を、上記3次元座標の原点をゴルフボールの中心位置に設定することにより特定する請求項1又は2に記載のクラブヘッドの挙動の計測方法。上記フェース面の3個以上のマークの3次元座標上の位置データから、フェース面の2次元座標を設定する請求項1から3のうちのいずれかに記載のクラブヘッドの挙動の計測方法。フェース面とゴルフボールとの接触時点における、接触点の3次元座標上のデータを得るステップと、この接触点の3次元座標上のデータをフェース面の2次元座標に転換するステップとを含む請求項4に記載のクラブヘッドの挙動の計測方法。上記マークの3次元座標上の位置ベクトルと、この位置ベクトルから得られる基準軸ベクトルとを用いてクラブヘッドの動ロフト角を算出する請求項1に記載のクラブヘッドの挙動の計測方法。上記マークの3次元座標上の位置ベクトルを用いてクラブヘッドのフェース角を算出する請求項1又は6に記載のクラブヘッドの挙動の計測方法。上記少なくとも2時点における、マークの3次元座標上の位置ベクトルを用いて、クラブヘッドのブロー角を算出する請求項1、6、7のうちのいずれかに記載のクラブヘッドの挙動の計測方法。上記少なくとも2時点における、マークの3次元座標上の位置ベクトルを用いて、クラブヘッドの進入角を算出する請求項1、6、7、8のうちのいずれかに記載のクラブヘッドの挙動の計測方法。上記クラブヘッドのクラウン部のトップ部に、帯状のマークを付するステップと、上記2時点における上記帯状マークの各位置データから、上記2時点間の帯状マークの回転角度を算出し、回転行列を得るステップと、上記2時点のうちの一方の時点における3個以上のマークの位置データと、他方の時点における、上記一方の時点におけるより少ない個数のマークの位置データとから、上記回転行列を用いて、他方の時点における未取得のマークの位置データを推定するステップとをさらに含んでいる請求項1から9のうちのいずれかに記載のクラブヘッドの挙動の計測方法。上記クラブヘッドのクラウン部のトップ部に、互いに離間した少なくとも2個のマークを付するステップと、上記2時点におけるトップ部の上記2個のマークを結ぶ仮想直線の各位置データから、上記2時点間の仮想直線の回転角度を算出し、回転行列を得るステップと、上記2時点のうちの一方の時点における、3個以上のマークの位置データと、他方の時点における、上記一方の時点におけるより少ない個数のマークの位置データとから、上記回転行列を用いて、他方の時点における未取得のマークの位置データを推定する推定するステップとをさらに含んでいる請求項1から9のうちのいずれかに記載のクラブヘッドの挙動の計測方法。ゴルフクラブのヘッドの挙動を計測するシステムであって、ゴルフクラブのヘッドのフェース面に付された少なくとも3個のマークと、このクラブヘッドの挙動を連続撮影する右カメラ及び左カメラと、情報処理装置とを備えており、上記右カメラが飛球方向の前方右側に配置されており、左カメラが飛球方向の前方左側に配置されており、上記カメラが同期させられて連続撮影可能にされており、上記情報処理装置が、連続撮影されたクラブ画像データから、上記マークの3次元座標上の位置を特定し、この位置データからフェース面を特定し、フェース面とゴルフボールとの位置関係を時系列に更新することにより、フェース面とゴルフボールとの接触時点を特定するクラブヘッドの挙動の計測システム。上記右カメラ及び左カメラの各光軸が、XYZの3次元直交座標系のXZ平面において、水平線に対して30°以上60°以下の角度をなしており、上記左カメラの光軸が、XY平面において、X軸に対して0°以上35°以下の角度をなしており、上記右カメラの光軸が、XY平面において、X軸に対してマイナス35°以上0°以下の角度をなしており、上記左右のカメラの光軸同士が、XY平面において、20°以上90°以下の角度をなしている請求項12に記載のクラブヘッドの挙動の計測システム。ゴルフクラブのヘッドの挙動を連続撮影する上カメラを備えており、この上カメラが、打撃されるゴルフボールの上方に配置されており、上記上カメラの光軸が、ゴルフボールの中心を通る鉛直線に対して80°以上100°以下の角度をなしている請求項12又は13に記載のクラブヘッドの挙動の計測システム。上記フェース面上のマークのうちの第一のマークがフェース面上の縦仮想直線よりトウ側に配置され、第二のマークが上記縦仮想直線よりヒール側に配置され、上記縦仮想直線はフェース面上をフェース面の中点を通って上下に延びる直線であり、上記第一及び第二のマークの一方を基準マークとして、第三のマークが、この基準マークの上方又は下方に配置されており、上記基準マークと第三のマークとを結ぶ仮想直線と、第一のマークと第二のマークとを結ぶ仮想直線とが、85°以上95°以下の角度をなしている請求項12から14のうちのいずれかに記載のクラブヘッドの挙動の計測システム。上記特許第2826697号公報に開示された装置では、ボールを打撃するインパクト以前のクラブヘッドの挙動は計測されるが、クラブヘッドとボールとの位置関係は計測も算出もされない。従って、クラブヘッドにおけるボールの打点を推定することはできない。上記特許第2950450号公報は、移動物体としてのゴルフボールの複数の2次元データから3次元データを得るDLT法を解説している。この公報に開示されたシステムでは、ゴルフボールの飛跡、速度、飛行の初期部分における回転は得られるが、クラブヘッドとボールとの位置関係は得られない。従って、クラブヘッドにおけるボールの打点を推定することはできない。上記特許第4109076号公報及び特許第4307511号公報に開示された方法及び装置では、コンピュータ上で作成されたクラブヘッドのCADデータを用いて、対象となるクラブヘッドの挙動が求められる。上記CADデータとクラブヘッドの実物の形状とは完全には一致しないため、測定誤差が発生しうる。また、クラブヘッドのCADデータを必要とするため、測定は容易ではない。この方法及び装置では、クラブヘッドとボールとの位置関係は計測も算出もされない。従って、クラブヘッドにおけるボールの打点を推定することはできない。上記特開2007-167549公報には、インパクト時のフェース角、ブロー角、ライ角、打点位置トウ-ヒール方向、上下方向と記載されてはいるが、その方法は具体的に開示されていない。もちろん、打点を定量的に得ることはできない。本発明は以上の現況に鑑みてなされたものであり、本発明の目的は、ゴルフクラブのインパクト直前又はインパクト時の挙動を計測する方法の提供にある。本発明に係るゴルフクラブのヘッドの挙動を計測する方法は、クラブヘッドのフェース面に少なくとも3個のマークを付するステップと、複数台のカメラにより、移動する上記クラブヘッドの少なくとも2時点において、上記マークの複数の2次元データを得るステップと、上記複数の2次元データから、マークの3次元座標上の位置を特定するステップと、3個以上のマークの3次元座標上の位置データからフェース面を特定するステップと、フェース面とゴルフボールとの位置関係を時系列に更新することにより、フェース面とゴルフボールとの接触時点を特定するステップとを含んでいる。好ましくは、上記少なくとも3個のマークの位置が等速直線運動をしていると仮定し、この3つの位置によって構成される面の法線ベクトルであって、ゴルフボールの中心を通る法線ベクトルを用いて、この面がゴルフボールと接触する時点を求める。好ましくは、上記フェース面とゴルフボールとの位置関係を、上記3次元座標の原点をゴルフボールの中心位置に設定することにより特定する。原点は自在に設定可能であるが、算定式を簡略化することによって計算時間を短縮するためには、原点をゴルフボールの中心位置に設定するのが好ましい。好ましくは、上記フェース面の3個以上のマークの3次元座標上の位置データから、フェース面の2次元座標を設定する。好ましくは、フェース面とゴルフボールとの接触時点における、接触点の3次元座標上のデータを得るステップと、この接触点の3次元座標上のデータをフェース面の2次元座標に転換するステップとをさらに含んでいる。好ましくは、上記マークの3次元座標上の位置ベクトルと、この位置ベクトルから得られる基準軸ベクトルとを用いてクラブヘッドの動ロフト角を算出する。好ましくは、上記マークの3次元座標上の位置ベクトルを用いてクラブヘッドのフェース角を算出する。好ましくは、上記少なくとも2時点における、マークの3次元座標上の位置ベクトルを用いて、クラブヘッドのブロー角を算出する。好ましくは、上記少なくとも2時点における、マークの3次元座標上の位置ベクトルを用いて、クラブヘッドの進入角を算出する。好ましくは、上記クラブヘッドのクラウン部のトップ部に、帯状のマークを付するステップと、上記2時点における上記帯状マークの各位置データから、上記2時点間の帯状マークの回転角度を算出し、回転行列を得るステップと、上記2時点のうちの一方の時点における、3個以上のマークの位置データと、他方の時点における、上記一方の時点におけるより少ない個数のマークの位置データとから、上記回転行列を用いて、他方の時点における未取得のマークの位置データを推定するステップとをさらに含んでいる。又は、上記クラブヘッドのクラウン部のトップ部に、互いに離間した少なくとも2個のマークを付するステップと、上記2時点におけるトップ部の上記2個のマークを結ぶ仮想直線の各位置データから、上記2時点間の仮想直線の回転角度を算出し、回転行列を得るステップと、上記2時点のうちの一方の時点における3個以上のマークの位置データと、他方の時点における、上記一方の時点におけるより少ない個数のマークの位置データとから、上記回転行列を用いて、他方の時点における未取得のマークの位置データを推定するステップとをさらに含むのが好ましい。本発明に係るゴルフクラブのヘッドの挙動を計測するシステムは、ゴルフクラブのヘッドのフェース面に付された少なくとも3個のマークと、このクラブヘッドの挙動を連続撮影する右カメラ及び左カメラと、情報処理装置とを備えており、上記右カメラが飛球方向の前方右側に設置されており、左カメラが飛球方向の前方左側に設置されており、上記カメラが同期させられて連続撮影可能にされており、上記情報処理装置が、連続撮影されたクラブ画像データから、上記マークの3次元座標上の位置を特定し、この位置データからフェース面を特定し、フェース面とゴルフボールとの位置関係を時系列に更新することにより、フェース面とゴルフボールとの接触時点を特定するように構成されている。好ましくは、上記右カメラ及び左カメラの各光軸が、Z軸が鉛直であるXYZの3次元直交座標系のXZ平面において、水平線に対して30°以上60°以下の角度をなしており、上記左カメラの光軸が、XY平面において、X軸に対して0°以上35°以下の角度をなしており、上記右カメラの光軸が、XY平面において、X軸に対してマイナス35°以上0°以下の角度をなしており、上記左右のカメラの光軸同士が、XY平面において、20°以上90°以下の角度をなしている。好ましくは、ゴルフクラブのヘッドの挙動を連続撮影する上カメラを備えており、この上カメラが、打撃されるゴルフボールの上方に設置されており、上記上カメラの光軸が、ゴルフボールの中心を通る鉛直線に対して80°以上100°以下の角度をなしている。好ましくは、上記フェース面上のマークのうちの第一のマークが、フェース面上をフェース面の中点を通って上下に延びる縦仮想直線よりトウ側に配置され、第二のマークが、上記縦仮想直線よりヒール側に配置されており、上記第一及び第二のマークの一方を基準マークとして、第三のマークが、この基準マークの上方又は下方に配置されており、上記基準マークと第三のマークとを結ぶ仮想直線と、第一のマークと第二のマークとを結ぶ仮想直線とが、85°以上95°以下の角度をなしている。本発明に係るクラブヘッドの挙動の計測方法によれば、インパクト時及びフェース面とゴルフボールとの接触点を推定することができる。以下、適宜図面が参照されつつ、好ましい実施形態に基づいて本発明が詳細に説明される。図1及び図2は、本発明に係るクラブヘッドの挙動の計測方法を実施するための計測システムの一実施形態を概略的に示す図である。この計測システム1により、ゴルファーPにスイングされているゴルフクラブの挙動が計測される。この計測では、スイング中のゴルフクラブの3次元の直交座標が時系列的に計測される。ゴルフクラブのヘッド2上の三箇所以上のポイントが3次元座標上で計測される。3次元座標の計測には、DLTDirectLinearTransformation法と呼ばれる方法が用いられるのが好ましい。このDLT法は、例えば、特許第2950450号公報、特開2004-344418号公報等に記載されている。DLT法は、バイオメカニクスの分野においても一般的に用いられている。DLT法は、異なる方向から見た複数の画像を用いて3次元座標を得る方法である。すなわち、対象物の異なる方向からの複数の2次元座標の位置データの総当たり的な組み合わせから、一つの3次元座標の位置データを得る方法である。上記複数の2次元座標の位置データは、もちろん、同一時点で得られたもの同時に撮影されたものである必要がある。このDLT法は、カメラの配置に関する制約が少なく、汎用性が高い。また、DLT法は、実空間におけるカメラの位置、光軸の方向、レンズの焦点距離等のカメラ定数に関する情報が不要とされうる利点を有する。DLT法では、3次元座標中の既知の点、すなわち、実空間での座標値が知られている点コントロールポイントの画像に基づいて3次元座標が再構築される。コントロールポイントの設定は、通常は、位置関係が定量的に決定された6個以上のマークが形成されたリファレンスフレームキャリブレーションフレームが用いられる。挙動の計測においては、クラブヘッド2の所定箇所にマークが取り付けられる。マークの取付部位については後述される。このマークの挙動が計測される。クラブヘッド2の挙動解析においては、マークの挙動の計測結果に基づいて解析がなされる。このマークの3次元座標が前述のDLT法によって得られる。図1及び図2に示される実施形態は、DLT法を用いた3次元挙動の計測システム1を例示している。この計測システム1は、複数のカメラ3、情報処理装置4、ゴルフクラブ5、及び、ゴルフクラブ5に付された複数のマークMを備えている。カメラ3は高速で連続撮影が可能な高速度カメラである。この計測システム1において、ゴルファーPの前後方向がY軸方向とされ、Y軸方向に垂直で且つ鉛直の方向がZ軸方向とされ、YZ両軸に垂直な、ゴルフボール打撃方向ボール飛行方向がX軸方向とされている。この3次元直交座標の原点は、定位置に載置されたゴルフボール以下、単にボールという6の中心に設定されている。X軸は、ボール6の飛んでいく向きがプラスである。Y軸は、ゴルファーPの立ち位置側がプラスである。Z軸は、上方がプラスである。この3次元座標において、クラブヘッド2の位置及び姿勢が特定される。複数のカメラ3は、情報処理装置4に接続されている。情報処理装置4は、図示されない制御部を有する。この制御部は、複数のカメラ3を同期して撮影可能に制御する。情報処理装置4は、DLT法による解析プログラムが記憶された記憶部と、演算部とを備える。典型的な記憶部は、ハードディスクである。典型的な演算部は、CPUである。複数のカメラ3は、互いに異なる位置に設置されている。本実施形態では、3台のカメラ3が用いられている。3台のカメラ3は同期してシャッターが切られる。1台のカメラ第一カメラと呼ぶ3Aは、インパクト直前のクラブヘッドの挙動及びクラブヘッド2のトップ部を撮影することができるように、ボール載置位置のほぼ真上に設置されている。他の2台のカメラ第二カメラ及び第三カメラと呼ぶ3B、3Cは、フェース面を撮影することができるように、ボールの飛球方向に設置されている。具体的には、以下のとおりである。なお、カメラ3の数は3台には限定されない。4台以上のカメラ3が用いられてもよい。第一カメラ3Aは、X軸上の原点ボール6の中心からやや後方の位置の上方に設置されている。本実施形態では、原点から80mm後方であり、床面から4000mm上方の位置X=-80mm、Y=0、Z=4000mmに設置されている。第一カメラ3Aは上カメラ3Aと呼べる。第二カメラ3Bは、X軸上の原点から4000mm前方で、Y軸プラス方向に1250mmの位置の、4000mm上方の位置X=4000mm、Y=1250mm、Z=4000mmに設置されている。第二カメラ3Bは左カメラ3Bと呼べる。第三カメラ3Cは、X軸上の原点から4000mm前方で、Y軸マイナス方向に1250mmの各位置の、4000mm上方の位置X=4000mm、Y=-1250mm、Z=4000mmに設置されている。第三カメラ3Cは右カメラ3Cと呼べる。第二及び第三のカメラ3B、3Cはそれぞれ、X軸方向に1000mmの範囲で移動可能にされている。すなわち、2台のカメラ3B、3Cはそれぞれ、原点から3000~4000mmの範囲で位置調節可能にされている。上記各寸法は例示であり、それらには限定されない。上記第二及び第三のカメラ3B、3Cは、その光軸同士が互いに直交する配置にされると、DLT法の数値解法で得られる解の精度が高くなる。しかし、光軸同士の交差角が大きくなるに伴い、撮影対象であるクラブヘッドのフェース面と光軸とのなす角度が小さくなる。その結果、複数のマークの映像が近接し、識別されにくくなる。逆に、光軸同士の交差角が小さくなるに伴い、DLT法の数値解法で得られる解の精度が低くなる。以上に説明された数値解析の精度の観点から、カメラ3の設置位置は、以下のとおり決定されるのが好ましい。第一カメラ3Aは、その光軸LAが、ボール6の中心を通る鉛直線VLに対して80°以上100°以下の角度をなしている。第二及び第三のカメラ3B、3Cはそれぞれ、その光軸LB、LCの方向が、ボール6の載置位置から地面床面に対して30°以上60°以下の角度をなして設置される。また、第二カメラ3Bは、その光軸LBの方向が、平面視でX軸に対して0°以上35°Y軸のプラス方向以下の角度をなして設置される。第三カメラ3Cは、その光軸LCの方向が、平面視でX軸に対してマイナス35°Y軸のマイナス方向以上0°未満の角度をなして設置される。ただし、平面視で第二カメラ3Bの光軸LBと第三カメラ3Cの光軸LCとのなす角度が、20°以上90°以下に設定される。X軸方向に第一カメラ3Aを挟んで、前後それぞれに第一ストロボ7Aが取り付けられている。第一ストロボ7Aは、Y=0、Z=4000mmにおける、X=-30mm及び-130mmの位置のそれぞれに取り付けられている。また、X=2000mm、Y=0、Z=4000mmの位置に2台の第二ストロボ7Bが取り付けられている。第一ストロボ7Aはほぼ鉛直方向下向きに照射する。第二ストロボ7Bはボール6近傍を照射する。これらのストロボ7は、カメラ3の作動と同期して照射する。第二ストロボ7Bは、効果的な照明の観点から、ボール6の載置位置から地面床面に対して40°以上80°以下の角度をなして設置される。ボール6の後方には、3台のカメラ3A、3B、3C及びストロボ7A、7Bの作動のタイミングを決定するトリガー装置が設置されている。このトリガー装置は、2台の光センサ8A、8Bから構成されている。光センサ8A、8Bは、クラブヘッド2の移動経路であるX軸を挟んで対向配置された検出光照射器9A、9Bと、検出光を受光しうる受光器10A、10Bとを有している。第一光センサ8AはX=-130mmの位置にあり、第二光センサ8BはX=-30mmの位置にある。スイングされたゴルフクラブ5のヘッド2が検出光を遮蔽したとき、受光器10A、10Bがこれを検知する。各受光器10A、10Bによる遮蔽検知により、それぞれの検知時点が基準時点となる。この基準時点に基づいてストロボ7及びカメラ3が作動する。図3及び図4には、クラブヘッド2におけるマークMの形成位置が示されている。図3にはウッド系のクラブヘッド2Wが示されている。図3aはクラブヘッド2Wの平面図であり、図3bはクラブヘッド2Wの正面図である。フェース面12Wに4つのマークMA、MB、MC、MDが形成されている。これらのマークとして、全てプリズムタイプの反射マークが採用されている。フェース面12Wには、フェース面の中点であるフェースセンタFoが設定されている。このフェースセンタFoを通ってフェース面12W上を上下方向に延びる縦仮想直線をYw軸とする。Yw軸に対して垂直にフェースセンタFoを通ってフェース面12W上をトウ-ヒール方向に延びる横仮想直線をXw軸とする。4つのマークMA、MB、MC、MDは、フェース面12W上の長方形の四頂点A、B、C、Dに位置している。2辺AC、BDはYw軸に平行であり、他の2辺AB、CDはXw軸に平行である。仮想直線ADと仮想直線BCとの交点がフェースセンタFoとなる。上記フェースセンタFoは、後述する位置データの取得のために妥当と思われる任意の位置に設定することが可能である。例えば以下のとおりである。第一には、フェース面の幅と高さとからフェースセンタが特定されうる。第二には、フェース面の面積重心がフェースセンタとして設定されうる。第三には、クラブヘッド全体の重心位置とフェース面との位置関係からフェースセンタが特定されうる。一例として、フェース面の幅と高さとからフェースセンタを設定する方法が以下に説明される。まず、クラブヘッドが、所定のリアルロフト角フェース角が0及びライ角となるように水平面に設置される。すなわち、ゴルフクラブが上記姿勢に固定される。この姿勢におけるクラブヘッドのフェース面の最トウ側点フェース最トウ側点が特定され、最ヒール側点フェース最ヒール側点が特定される。フェース最トウ側点とフェース最ヒール側点との水平方向離間距離がフェース幅として特定される。フェース幅の中央点を通り、水平面に垂直で且つその中央位置でフェース面に垂直な面垂直面が特定される。この垂直面とフェース面とのフェース面上の交線が特定される縦ライン。フェース面上のこの縦ラインの最上側点フェース最上側点が特定される。次いで、フェース面上の縦ラインの最下側点フェース最下側点が特定される。フェース最上側点とフェース最下側点とを結ぶ仮想直線の中央点がフェースセンタと設定される。このフェースセンタ設定方法は一例である。ウッド系クラブヘッド2Wには、そのクラウン部の頂部トップ部13Wに帯状のマークMEが形成されている。図3aに示されるように、トップ部13Wは、クラウン部の上面におけるフェース面12Wに近接した部位である。帯状のマークMEは、トップ部13W上をトウ-ヒール方向に延在している。この帯状マークMEは、後述されるように、フェース面12Wのマークのうち2つしか撮影されなかった場合、撮影されなかったマークの位置を推定するために用いられる。この帯状マークMEは、前述した第一カメラ3Aによって容易に撮影されうる。図4にはサンドウェッジと呼ばれるロフト角の大きいアイアン系クラブのヘッド2Aが示されている。図4aはヘッド2Aの平面図であり、図4bはヘッド2Aの正面図である。フェース面12Aに3つのマークMA、MB、MCが形成されている。これらのマークには、ガラスビーズタイプの反射テープが反射マークとして採用されている。ロフト角が小さいクラブの場合は、反射マークによって反射された光が撮影されにくいので、プリズムタイプの反射マークを用いるのが好ましい。一方、ロフト角が大きいクラブの場合は、反射マークによって反射された光が撮影されやすいので、とくにプリズムタイプの反射マークを用いる必要がない。従って、ロフト角が大きいクラブの場合は、より低コストのガラスビーズタイプの反射テープ等を用いるのが好ましい。具体的には、プリズムタイプの反射マークは、ロフト角が25°以下のクラブにとって好ましい。ロフト角が20°以下のクラブにとってはさらに好ましい。ロフト角が15°以下のクラブにとっては特に好ましい。一方、コスト削減の観点から、ガラスビーズタイプの反射テープ等は、ロフト角が16°以上のクラブにとって好ましい。ロフト角が21°以上のクラブにとってはさらに好ましい。ロフト角が26°以上のクラブにとっては特に好ましい。上記マークMA、MB、MCの全てが、フェース面12Aのブラスト処理された中央側部分に形成されている。トウ-ヒール方向両端側では、フェース面A自体が反射しやすいため、マークの識別が難しくなる可能性があるからである。フェース面12AにはフェースセンタFoが設定されている。上記ウッド系クラブヘッド2Wにおけると同様にして、横仮想直線であるXa軸及び縦仮想直線であるYa軸が設定される。3つのマークMA、MB、MCは、フェース面12A上の直角三角形の三頂点A、B、Cに位置している。底辺ABはXa軸に一致している。フェースセンタFoは底辺ABの中点である。辺ACはYa軸に平行である。アイアン系クラブのヘッド2AのフェースセンタFoも、前述した各種の設定方法によって設定されうる。また、以下の方法によっても設定されうる。まず、クラブごとに予め決まっているライ角となるように、対象のクラブの姿勢が固定される。フェース面のフェース材例えば、図4bにおけるA及びCを通る直線AC線と、Bを通り且つ上記AC線に平行な線とに囲まれた範囲の幅方向中央点が特定され、この中央点を通る鉛直線中央縦ラインが特定される。フェース材上の中央縦ラインの最上側点フェース最上側点が特定される。次いで、フェース材上の中央縦ラインの最下側点フェース最下側点が特定される。フェース最上側点とフェース最下側点とを結ぶ仮想直線の中央点がフェースセンタと設定される。このフェースセンタ設定方法は一例である。アイアン系クラブヘッド2Aには、そのクラウン部のトップ部13Aに2つのマークMF、MGが形成されている。図4aに示されるように、トップ部13Aは、クラウン部の上面におけるフェース面12Aに近接した部位である。上記2つのマークMF、MGは、トップ部13A上をトウ-ヒール方向に互いに離間して形成されている。この2つのマークMF、MGは、後述されるように、フェース面12Aのマークのうち2つしか撮影されなかった場合、撮影されなかったマークの位置を推定するために用いられる。この目的のために、上記2つのマークMF、MGを結ぶ仮想直線が利用される。この2つのマークMF、MGは、前述した第一カメラ3Aによって容易に撮影されうる。本実施形態においては、ウッド系クラブヘッド2Wのフェース面12には4つのマークが形成されているが、4つには限定されない。3つでもよく、5つ以上でもよい。アイアン系クラブヘッド2Aのフェース面12には3つのマークが形成されているが、3つには限定されない。4つ以上でもよい。フェース面上のマークは面を構成しうる少なくとも3つの位置に必要である。すなわち、一直線上の配置ではない少なくとも3つのマークが必要である。これは、クラブヘッド2のフェース面を特定する必要があるからである。各マークMA、MB、MC、MDの位置も、上記実施形態におけるものには限定されない。フェース面12上のマークのうちのいずれか一のマーク第一のマークが、上記Ya軸及びYw軸以下、Ya軸で代表させるよりトウ側に配置され、他のいずれか一のマーク第二のマークが、Ya軸よりヒール側に配置され、さらに他のいずれか一のマーク第三のマークが、上記第一及び第二のマークのうちの一方のマーク基準マークというの上方又は下方に配置され、この基準マークと第三のマークとを結ぶ仮想直線と、第一のマークと第二のマークとを結ぶ仮想直線とは、85°以上95°以下の角度をなしているのが好ましい。最も好ましい角度は90°である。基準マークとしては、本実施形態では、図3b及び図4bにおけるマークMAが該当する。このような第一、第二、第三のマークが、下記アからオのうちの少なくとも一つの要件を満たしているのがさらに好ましい。ア上記基準マークは、Ya軸よりトウ側に配置されるのが好ましい。一般的にフェース面はヒール側よりトウ側の方が広いため、複数のマークを広く分散配置することができるからである。その結果、隣接するマークが共にボールに隠れてしまう事態が回避されうる。イ上記第一のマーク及び第二のマークは、それらを結ぶ仮想直線がYa軸によって2等分されるように配置されるのが好ましい。フェースセンタFoが容易に計算されうるからである。さらに、ボールがフェース面の中央付近に衝突したとき、マークがボールに隠れてしまう事態が回避されうるからである。ウ上記第一のマーク及び第二のマークは、それらを結ぶ仮想直線がYa軸と直交するように配置されているのが好ましい。3次元座標が2次元座標に変換された際に、打点がフェースセンタFoからの上下左右の各方向の距離によって定義されうるからである。その結果、フェース面のどの部位にボールが衝突したかがイメージされやすい。エ上記第一のマークと第二のマークとの離間距離両マークを結ぶ仮想直線の実長さ。以下同じは、40mm以上が好ましく、45mm以上がさらに好ましい。隣接するマークが共にボールに隠れてしまう事態が回避されうるからである。この離間距離は、フェース面の大きさの範囲内で設定されるが、通常、80mm以下、さらには70mm以下に設定されうる。本実施形態では、上記離間距離は、ウッド系クラブヘッド2W及びアイアン系クラブヘッド2Aともに60mmにされている。オ上記基準マークと第三のマークとの離間距離は、22mm以上が好ましく、24mm以上がさらに好ましい。隣接するマークが共にボールに隠れてしまう事態が回避されうるからである。この離間距離は、フェース面の大きさの範囲内で設定されるが、通常、50mm以下、さらには40mm以下に設定されうる。本実施形態では、上記離間距離は、ウッド系クラブヘッド2Wでは30mm、アイアン系クラブヘッド2Aでは25mmにされている。以上に説明された計測システム1により、DLT法が用いられて、クラブヘッド2の位置、姿勢等の3次元データが時系列で得られる。得られた上記3次元データに基づいて、フェース面12とボール6との3次元座標上の位置関係が時系列に更新され、フェース面とボールとの接触時点インパクト時が特定される。最終的に、インパクト時のクラブヘッド2の挙動が推定される。上記挙動とは、例えば、クラブヘッド2の動ロフト角、フェース角開き角、ブロー角、進入角、フェース面12におけるボール6の打点等である。図5、図6から図13を参照しつつ、インパクト時におけるクラブヘッド2の挙動の推定の一例を説明する。図5において、前述した第一から第三の3台のカメラ3A、3B、3Cによってインパクト直前のクラブヘッドが同時に撮影されるSTEP1。この撮影は、例えば2つの時点のそれぞれにおいて行われる。次いで、撮影された画像から、フェース面12の複数のマークの位置が、画像処理によって抽出されるSTEP2。上記複数のマークとは、ウッド系クラブヘッドでは4個のマーク、アイアン系クラブヘッドでは3個のマークを意味している。しかし、ウッド系クラブヘッドについても、3個のマークが採用されてもよい。3台のカメラのいずれの撮影画像においても、各マークの位置が抽出された場合STEP3、各位置は、3台計測用の補正式DLT法において規定される近似式によって3次元座標に変換されるSTEP4。3台のうちの2台のカメラの撮影画像においてのみ、各マークの位置が抽出された場合STEP5、各位置は、2台計測用の補正式DLT法において規定される近似式によって3次元座標に変換されるSTEP6。又は、例えば第二及び第三の2台のカメラのみが使用された場合でもSTEP5、各位置は、2台計測用の補正式DLT法において規定されるによって3次元座標に変換されるSTEP6。このDLT法において規定される補正式は、情報処理装置4に記憶されている。3台計測用の補正式は、第一カメラ3A、第二カメラ3B及び第3カメラ3Cで撮影されたコントロールポイントの画像に基づいて予め求められている。第一カメラ3A及び第二カメラ3Bの2台計測用の補正式は、第一カメラ3A及び第二カメラ3Bで撮影されたコントロールポイントの画像に基づいて予め求められている。第一カメラ3A及び第三カメラ3Cの2台計測用の補正式は、第一カメラ3A及び第三カメラ3Cで撮影されたコントロールポイントの画像に基づいて予め求められている。第二カメラ3B及び第三カメラ3Cの2台計測用の補正式は、第二カメラ3B及び第三カメラ3Cで撮影されたコントロールポイントの画像に基づいて予め求められている。このSTEP1からSTEP6の2次元座標から3次元座標への変換は公知のDLT法によってなされるため、詳細な説明は省略される。2時点のそれぞれにおけるフェース面上の3つ以上のマークの3次元位置情報が取得できたか否かが判定されるSTEP7。上記3つ以上の3次元位置情報が取得できた場合には、後述されるように2時点のフェース面の位置から、インパクト時のフェース面12の位置が推定され、フェース面12とボール6との接触点の3次元座標上の位置が計算されるSTEP8。STEP8においては、フェース面12とボール6との3次元座標上の位置関係が時系列に更新されることにより、フェース面12とボール6との接触時点を特定され、接触点が算出される。2時点のいずれか一方において3つ以上のマークの3次元位置情報が取得できているが、他方においては上記一方において取得されたマーク位置情報より少ない個数の位置情報しか得られていない場合にはSTEP9、他方において取得されなかったマークの3次元座標位置が推定されるSTEP10。具体的には、ウッド系クラブヘッド2Wのフェース面12Wに形成された4つのマークMについて、最初の時点においては4つのマークMの3次元位置情報が取得できたが、後の時点においては2つのマークがボール6の陰になって撮影されず、その3次元位置情報が取得できなかった場合であるSTEP9。この場合には、トップ部13Wの帯状マークMEの2時点間の回転角に基づいて、後の時点において取得できなかった2つのマークの3次元位置情報が推定されるSTEP10。この推定方法については後述される。2時点のフェース面12の位置から、インパクト時のフェース面12の位置が推定され、フェース面12とボール6との接触点の3次元座標上の位置が計算されるSTEP8。算出されたフェース面12とボール6との接触点の3次元座標が、2次元座標であるフェース面座標に変換され、インパクト時の打点を含むクラブヘッドの挙動が推定されるSTEP11。[フェース面とボールとの接触点の推定]図6及び図7を参照しつつ、以下に、上記STEP8におけるフェース面12とボール6との接触点の3次元座標上の位置の計算について説明される。図6に示されように、例示されたクラブヘッドのフェース面上の3箇所A、B、Cには、それぞれマークが形成されている。上記接触点の3次元座標位置の計算の基本的な考え方は、上記3つの点A、B、Cが等速直線運動をしていると仮定し、3つの点A、B、Cによって構成される面がボールと接触する時点を求め、その時点での上記位置A、B、C及び上記接触点Qの各座標を求めることである。前述した2時点のうちの最初の時点をt=0とし、後の時点をt=1とする。図示のごとく、t時点における点A、B、Cの座標位置がAt、Bt、Ctとされる。t=0時点における点A、B、Cの座標位置がA0、B0、C0とされる。t=1時点における点A、B、Cの座標位置がA1、B1、C1とされる。上記3次元座標の原点0、0、0は、図6に示されるボール6の中心点である。t=0とt=1との時間間隔はTとされる。この時間Tは遮光時間間隔の計測値である。遮光時間間隔とは、クラブヘッド2が、図2に示される第一光センサ8Aの検出光を遮蔽してから、第二光センサ8Bの検出光を遮蔽するまでの時間間隔である。以上から、上記座標位置At、Bt、Ctは下記の式1で表される。At=A1A1-A0/T×δtBt=B1B1-B0/T×δt1Ct=C1C1-C0/T×δt上記式1において、δtは計算時間間隔である。δtの初期値は例えば100μsマイクロ秒に設定され、終了値が例えば300μsに設定され、増加ステップが例えば1μsに設定される。ボールの半径が例えばrに設定される。図7のフローチャートにおいて、上記式1によって座標位置At、Bt、Ctが計算されるSTEP21。時刻tにおけるABベクトル及びACベクトルが計算される。そして、上記両ベクトルの外積NNx、Ny、Nzが計算されるSTEP22。このベクトルの外積Nは法線ベクトルとなる。この法線ベクトルは単位ベクトルである。この法線ベクトルはフェース面に垂直なベクトルフェース面ベクトルである。上記At、Bt、Ctと、法線ベクトルNとを用いて、下記の式2により、ボールの中心0、0、0とフェース面との最短距離rrが計算されるSTEP23。rr=-Nx×Atx-Ny×Aty-Nz×Atz2上記最短距離rrが、ボールの半径r以下になったか否かが判定されるSTEP24。半径r以下になったとき、その時点のδtがインパクト時のδtであると決定され、インパクト時の座標位置At、Bt、Ctが決定されるSTEP25。次いで、インパクト時の3次元座標上の打点Qx、Qy、Qzが下記の式3によって決定されるSTEP26。Qx=-Nx×rQy=-Ny×r3Qz=-Nz×r上記STEP24において、上記最短距離rrがボールの半径r以下にならない間、及び、δtの累積が300μsにならない間は、δtがδt1μsに置き換えられてSTEP28、上記STEP21からSTEP24が繰り返される。以上の演算によって求められた打点Qx、Qy、Qzは3次元座標上の位置である。従って、この打点は、図5中のSTEP11に従って、フェース面の2次元座標に変換される必要がある。上記3次元座標上の打点の2次元座標への変換は、図8及び図9を参照しつつ、アイアン系クラブヘッド2Aの場合と、ウッド系クラブヘッド2Wの場合とに分けて、以下に説明される。[アイアン系クラブヘッド2Aの打点の2次元変換]インパクト時の3つのマークの位置A、B、Cは、それぞれXa、Ya、Za、Xb、Yb、Zb、Xc、Yc、Zcとされ、打点はXq、Yq、Zqとされる。この打点Xq、Yq、Zqは、前述した打点Qx、Qy、Qzと同一位置である。前述のとおり、アイアン系のフェース面12Aのマーク位置は、AB間の中点がフェースセンタFoとなるように定められている図4、図8。図8は、FX軸及びFY軸からなるフェース面12A上の2次元直交座標を示している。フェースセンタFo及びB点の3次元座標位置が既知であり、打点Qの3次元座標位置も既知である。これらにより、FoQベクトル及びFoBベクトルそれぞれの大きさ|FoQベクトル|、|FoBベクトル|、並びに、FoQベクトルとFoBベクトルとのなす角度α図8は、下記の式4、式5、式6によって求められる。なお、sqrtは、内の平方根を示す。また、数式中に示されるベクトルには、その上部に矢印→を付してベクトルである旨を表すこととする。打点QのFX軸フェース面上のY軸上の位置は、|FoQベクトル|cosαによって求まる。FY軸フェース面上のY軸とFoQベクトルとのなす角度βは、90°-αであるから、打点QのFY軸上の位置は、|FoQベクトル|cosβによって求まる。このように、図5中のSTEP11における打点Qの2次元座標への変換は容易である。[ウッド系クラブヘッド2Wの打点の2次元変換]図9は、FX軸及びFY軸からなるフェース面12W上の2次元直交座標を示している。インパクト時の4つのマークの位置A、B、C、Dは、それぞれXa、Ya、Za、Xb、Yb、Zb、Xc、Yc、Zc、Xd、Yd、Zdと置かれ、打点はXq、Yq、Zqと置かれる。この打点Xq、Yq、Zqは、前述した打点Qx、Qy、Qzと同一位置である。前述のとおり、ウッド系のフェース面12Wのマーク位置は、長方形ABCDの中点がフェースセンタFoとなるように定められている図3、図9。A、B、C、Dの各点が既知であるため、フェースセンタFoの座標位置Fox、Foy、Fozは、下記の式7によって求められる。Fox=XaXb/2XcXd/2/2Foy=YaYb/2YcYd/2/27Foz=ZaZb/2ZcZd/2/2以上の各値により、FoQベクトル及びFoBベクトルそれぞれの大きさ|FoQベクトル|、|FoBベクトル|並びに、FoQベクトルとFoBベクトルとのなす角度α図9は、下記の式8、式9、式10によって求められる。打点QのFX軸上の位置は、|FoQベクトル|cosαによって求まる。図9中のFY軸とFoQベクトルとのなす角度βは、90°-αであるから、打点QのFY軸上の位置は、|FoQベクトル|cosβによって求まる。このように、図5中のSTEP11における打点Qの2次元座標への変換は容易である。[動ロフト角の推定]図10を参照しつつ、インパクト直前及びインパクト時の動ロフト角θの推定が以下に説明される。図10は、フェース面に設定された2次元座標と、3次元座標とを示している。動ロフト角とは、スイング等によるクラブヘッドの移動に応じて変化するロフト角である。ここで、図10に示されるように、動ロフト角θとは、X軸とABベクトルとの外積AJと、ACベクトルとのなす角度θである。従って、AJベクトルが動ロフト角θの基準軸となる。ここで、ABベクトル=Xb-Xa、Yb-Ya、Zb-Zaであり、ACベクトル=Xc-Xa、Yc-Ya、Zc-Zaであり、X軸ベクトル=1、0、0である。ABベクトルとX軸ベクトルとの外積を求めることによりAJベクトルが算出される。そして、動ロフト角θは、下記の式11によって求められる。Xc<Xaの場合、動ロフト角θはプラスとなる。[フェース角の推定]図11を参照しつつ、インパクト直前及びインパクト時のフェース角開き角γの推定が以下に説明される。図11は、フェース面に設定された2次元座標と、3次元座標とを合わせて示している。ここで、フェース角γとは、BAベクトルのYZ平面に投影されたベクトルBGと、BAベクトルとのなす角度γである。ここで、BAベクトル=Xa-Xb、Ya-Yb、Za-Zbである。BAベクトルのYZ平面へ投影されたベクトルBG=0.0、Ya-Yb、Za-Zbである。そして、フェース角γは、下記の式12によって求められる。Xa<Xbの場合、フェース角γはプラスであり、いわゆるオープンである。それ以外であればフェース角γはマイナスであり、いわゆるクローズである。上記フェース角の推定方法と同様な手法により、動的ライ角を推定することも可能である。動的ライ角を推定する際には、BAベクトルが投影されるのはYZ平面ではなくXY平面であり、BAベクトルのXY平面に投影されたベクトルと、BAベクトルとのなす角度と、クラブの所定のライ角とにより、動的ライ角が求められる。[ブロー角の推定]次に、インパクト直前及びインパクト時のブロー角ηの推定が以下に説明される。ここで、ブロー角ηとは、前述したt=0及びt=1におけるフェースセンタのベクトルFot0-Fot1のXY平面に投影されたベクトルFot0’-Fot1’と、フェースセンタベクトルFot0-Fot1とのなす角度ηである。まず、t=0及びt=1におけるフェースセンタの移動距離が求められる。フェースセンタFoの座標はXaXb/2、YaYb/2、ZaZb/2である。t=0におけるFoの座標とt=1におけるFoの座標とから、上記フェースセンタの移動距離が求まる。フェースセンタベクトルFot0-Fot1は、下記の式13で表される。フェースセンタのベクトルのXY平面に投影されたベクトルFot0’-Fot1’は、下記の式14で表される。従って、下記の式15により、ブロー角ηが求められる。Za1Zb1/2>Za0Zb0/2の場合、ブロー角はプラスであり、いわゆるアッパーブローである。それ以外ではブロー角はマイナスであり、いわゆるアンダーブローである。[進入角の推定]次に、インパクト直前及びインパクト時の進入角ξの推定が以下に説明される。ここで、進入角ξとは、t=0及びt=1におけるフェースセンタ同士を結ぶ線分の座標のXY平面に投影された座標と、X軸とのなす角度ξである。t=0において、XY平面に投影されたフェースセンタFoxy0は、下記の式16で表される。t=1において、XY平面に投影されたフェースセンタFoxy1は、下記の式17で表される。従って、下記の式18により、進入角ξが求められる。Foxy0=Xa0Xb0/2,Ya0Yb0/2,016Foxy1=Xa1Xb1/2,Ya1Yb1/2,017tanξ=Ya1Yb1/2-Ya0Yb0/2/Xa1Xb1/2-Xa0Xb0/218Ya1Yb1/2>Ya0Yb0/2の場合、進入角はマイナスであり、いわゆるインサイドアウトである。それ以外では進入角はプラスであり、いわゆるアウトサイドインである。以上のとおり、上記計測システム1を用いたクラブヘッドの挙動計測方法により、インパクト直前及びインパクト時のクラブヘッドの挙動が計測され、座標上に特定され、推定される。すなわち、クラブヘッドの位置関係が定量化される。[未取得位置の推定]次に、図5中のSTEP9及びSTEP10における、2時点のうちの1時点で取得されなかったマークの位置の推定について、以下に説明される。図12には、t=0及びt=1の2時点におけるウッド系クラブヘッドのフェース面12Wが示されている。t=0においては、4つのマークの位置A、B、C、Dが全て抽出された。t=1においては、トウ側の2つのマークがボールの陰に隠れたために撮影されず、この2つの位置A’、C’が抽出されなかった。ヒール側の2つの位置B’、D’は抽出された。図12には、抽出されなかった位置A’、C’は白地で示されている。トップ部13Wに形成された帯状マークMEは、第一カメラ3A図1及び図2参照によってほぼ確実に撮影されうる。帯状マークMEの長手方向の端部位置はそれぞれE点及びF点とされる。上記2位置A’、C’の推定の基本的な考え方は、t=0におけるDCベクトルDCx0、DCy0、DCz0及びBAベクトルBAx0、BAy0、BAz0を用いて、t=1におけるD’C’ベクトルDCx1、DCy1、DCz1及びB’A’ベクトルBAx1、BAy1、BAz1を求めることである。その際、帯状マークMEの端部位置のFEベクトルの、t=0からt=1への回転角度が用いられる。第一に、t=0におけるFEベクトルFEx0、FEy0、FEz0及びt=1におけるF’E’ベクトルFEx1、FEy1、FEz1が求められる。ここで、である。第二に、t=0におけるFEベクトルと、t=1におけるF’E’ベクトルとの外積ベクトルVVx、Vy、Vzが求められる。ここで、Vx、Vy、Vzはそれぞれ、下記の式19、式20、式21によって求められる。第三に、t=0におけるFEベクトルと、t=1におけるF’E’ベクトルとのなす角度φが、下記の式22によって求められる。なお、acosは、アークコサインを意味する。第四に、t=0からt=1への帯状マークMEの回転行列線分FEの回転行列Mtが、下記の式23によって求められる。第五に、t=0におけるDCベクトルDCx0、DCy0、DCz0及びBAベクトルBAx0、BAy0、BAz0が求められる。ここで、である。上記DCベクトル及びBAベクトルのそれぞれに上記回転行列Mtが乗ぜられる。その結果、t=1における推定されたD’C’ベクトル及びB’A’ベクトルが求められる。下記の式24及び式25のとおりである。取得することができた位置B’に、上記B’A’ベクトルを加えることにより、未取得であった位置A’t=1におけるマークMAの位置が推定される。また、取得することができた位置D’に、上記D’C’ベクトルを加えることにより、未取得であった位置C’t=1におけるマークMCの位置が推定される。このようにして得られた推定値を用いることにより、図5中のSTEP8以降の演算が実行される。その結果、t=0及びt=1の両時点において3つ以上の位置が抽出し得た場合と同様に、インパクト直前及びインパクト時におけるクラブヘッドの挙動が推定されうる。図13に示されるウッド系クラブヘッドでは、t=0においては、4つのマークの位置A、B、C、Dが全て抽出された。t=1においては、図12とは異なり、ヒール側の2つのマークがボールの陰に隠れたために撮影されず、この2つの位置B’、D’が抽出されなかった。この場合でも、未取得の位置B’、D’を推定するために、図12を参照しつつ説明された手法と同様の手法が用いられる。従って、詳細な説明は省略される。図13に示される例については、ベクトルの向きが、図12におけるのとは逆になる。すなわち、t=0におけるCDベクトルCDx0、CDy0、CDz0及びABベクトルABx0、ABy0、ABz0を用いて、t=1におけるC’D’ベクトルCDx1、CDy1、CDz1及びA’B’ベクトルABx1、ABy1、ABz1が求められる。その際、帯状マークMEの端部位置のベクトルEFの、t=0からt=1への回転角度が用いられる。以上説明された未取得位置の推定では、各マークの3次元座標上の位置データが用いられたが、かかる手法には限定されない。各マークの2次元座標上の位置データが用いられてもよい。具体的には、例えば、クラブヘッドの真上に位置する第一カメラ3Aによって撮影されたデータから得られた2次の回転行列が用いられてもよい。計算時間の短縮が重要視される場合には2次元位置データが用いられ、精度が重要視される場合には3次元位置データが用いられるのがよい。"