ゴルフクラブをスイングした時のシャフトの撓りを解析するゴルフクラブシャフトの3次元変形解析方法であって、シャフト上に複数のマーカーが離間して設けられたゴルフクラブのスイングを複数のカメラで撮影する工程と、前記カメラで撮影された画像から前記シャフト上のマーカーの3次元座標を算出する工程と、前記シャフト上に設けられた複数のマーカーのうち、任意の隣り合う2つのマーカーを選択して基準ベクトルを作成する工程と、前記基準ベクトルの延長線と他のマーカー間の距離を算出する工程と、を含むことを特徴とするゴルフクラブシャフトの3次元変形解析方法。前記基準ベクトルの延長線と前記マーカー間との距離をベクトル計算によって算出することを特徴とする請求項1に記載のゴルフクラブシャフトの3次元変形解析方法。しかしながら、ゴルフクラブシャフトは、撓ったり、ねじれたりして非常に複雑な挙動をすることが知られており、特許文献1に示されている2次元のシャフトの撓りを求める方法では、シャフトの詳細な撓り変形を求めることができないという問題がある。また、特許文献2に記載のシャフト挙動自動計測システムは、金属体の3次元座標を算出するものであり、特許文献2には、シャフトの撓り変形を求める方法は開示されていない。本発明は、このような問題点に鑑みてなされたものであり、本発明の目的は、ゴルフクラブシャフトの変形を詳細に解析できる3次元変形解析方法を提供することにある。上記目的を達成するため、第1の発明は、ゴルフクラブをスイングした時のシャフトの撓りを解析するゴルフクラブシャフトの3次元変形解析方法であって、シャフト上に複数のマーカーが離間して設けられたゴルフクラブのスイングを複数のカメラで撮影する工程と、前記カメラで撮影された画像から前記シャフト上のマーカーの3次元座標を算出する工程と、前記シャフト上に設けられた複数のマーカーのうち、任意の隣り合う2つのマーカーを選択して基準ベクトルを作成する工程と、前記基準ベクトルの延長線と他のマーカー間の距離を算出する工程とを含むことを特徴とする。第2の発明は、前記基準ベクトルの延長線と前記マーカー間との距離をベクトル計算によって算出することを特徴とする。第1の発明によれば、3次元でのゴルフクラブシャフトの撓りを計測できるので、ゴルフクラブシャフトの変形を詳細に解析できる。第2の発明によれば、基準ベクトルの延長線と前記マーカー間との距離をベクトル計算によって算出するので、数学的処理により容易に撓りを解析できる。本発明の実施の形態について説明する前に、まず、本発明によるゴルフクラブシャフトの3次元変形解析方法において用いられるカメラシステムについて説明する。図1は、上述したカメラシステムを用いた撮影システムの一実施形態を説明する図である。なお、ここでは、ゴルファがゴルフクラブによりゴルフボールを打撃するときのゴルフスイングを撮影する場合を例にして撮影システムを説明するが、図1に示す撮影システムは、野球のバットスイングや、テニスのラケットスイング等を撮影する場合にも、適用できるものである。図1に示す撮影システム11は、カメラシステム15と、計測装置13コンピュータからなる。カメラシステム15は、ゴルファがゴルフクラブによりゴルフボールを打撃するときのゴルフスイングをステレオ撮影法により撮影する2台の第1カメラ12A及び第2カメラ12Bと、第1カメラ12A及び第2カメラ12Bを固定する治具24からなる。計測装置13は、第1カメラ12A及び第2カメラ12Bから取得した画像から、ゴルフスイングの軌跡を取得して、表示する。第1カメラ12A及び第2カメラ12Bは、ゴルフスイングの3次元映像を取得できるように、所定間隔離間して配置されている。また、第1カメラ12A及び第2カメラ12Bは、座標が既知の対象物を複数撮影し、座標が既知の対象物の映像を取得して、第1カメラ12A及び第2カメラ12Bの間の相対的な位置関係の算出キャリブレーションが行われており、算出された相対的な位置関係を保持できるように冶具24で固定されている。図2は、カメラシステムを示す図であり、図3は、カメラシステムのカメラの配置を示す図である。第1カメラ12Aと第2カメラ12Bのレンズの光軸間の距離を1100mmとしている。第1カメラ12Aと第2カメラ12Bのレンズの光軸間の距離は、800~1500mmが好ましい。第1カメラ12A及び第2カメラ12Bは、レンズの光軸を、カメラ同士を結ぶ線分に対して直角の方向から内側に15°回転させた状態で冶具24に固定されている。第1カメラ12A及び第2カメラ12Bのレンズの光軸を内側に回転させる角度は10°~30°であることが好ましい。上述したカメラシステムでは、治具に固定された2台のカメラは、既にキャリブレーションが行われているので、計測の度にキャリブレーションを行わなくても良い。また、2台のカメラは、治具に固定されているので、撮像範囲を調整する際に2台同時に調整可能である。また、カメラシステム15は、被験者ゴルファ14によるゴルフスイングを撮像するために、被験者14が使用するゴルフクラブが、スイング中の各段階バックスイング、ダウンスイング、フォロースルーにおいて、撮像できるように、適宜被験者14から離間して配置される。ここで、上述のように、第1カメラ12A及び12Bの間の間隔は固定されているため、被験者の身長や体格等に応じて第1カメラ12A及び12Bと、被験者14との間隔を適宜調整したとしても、同一平面上で測定したスイング軌跡を比較することができる。図1では、被験者14が打撃を行うゴルフボールの中心位置を原点とする。そして、原点を通り、水平面に対して垂直な方向の軸をz軸とする。さらに、原点をとおり、被験者14の両肩を結んだ線に対して垂直であり、かつ水平面上に延在する軸をx軸とする。さらに、原点を通り、x軸及びz軸に対して垂直な方向打ち出し方向の軸をy軸とする。図4は、第1カメラ12A及び第2カメラ12Bと、ゴルフクラブの部位が描く軌跡がのっている平面スイングプレーンとの位置関係を説明する図である。ゴルフスイング撮影の際は、第1カメラ12Aのレンズの光軸撮像素子面に対する垂線22Aと、第2カメラ12Bのレンズの光軸撮像素子面に対する垂線22Bとを含む面が、計測したいゴルフクラブ26の部位が描くスイング軌跡25が乗っている平面スイングプレーン20に対して垂直となるように、カメラシステム15を被験者14の正面に配置する。図5は、カメラの設定角度を説明する図である。図5では、スイングプレーン20が、ゴルフボール21の位置において、水平面23に対して平均で48°傾いているとき、カメラのレンズの光軸撮像素子面に対する垂線22A、22Bを含む面が、スイングプレーン20に対して垂直となるようにすると、第1カメラ12A及び第2カメラ12Bは、水平面23に対して平均で42°傾けることになることを示している。スイングプレーン20の水平面23に対する角度は、被験者ゴルファ14によって、また、ゴルフクラブによって異なる。スイングプレーン20の水平面23に対する角度が30°~80°であるとすると、第1カメラ12A及び第2カメラ12Bの水平面23に対する傾き角度は、10°~60°であることが好ましい。なお、上述した実施形態では、キャリブレーションを行ったカメラ2台で、ゴルフスイングの軌跡を撮影したが、カメラ3台以上を治具24に固定し、キャリブレーションを行って、3台以上のカメラでゴルフスイングの軌跡を撮影するようにしても良い。次に、上述したカメラシステムを用いた本発明によるゴルフクラブシャフトの3次元変形解析方法の実施の形態について説明する。図6は、本発明によるゴルフクラブシャフトの3次元変形解析方法が適用される解析システムの一実施形態を説明する図である。図6に示す解析システム31は、カメラシステム15と、解析装置33コンピュータからなる。カメラシステム15は、ゴルファがゴルフクラブによりゴルフボールを打撃するときのゴルフスイングをステレオ撮影法により撮影する2台の第1カメラ12A及び第2カメラ12Bと、第1カメラ12A及び第2カメラ12Bを固定する治具24からなる。解析装置33は、第1カメラ12A及び第2カメラ12Bから取得した画像から、ゴルフクラブの撓りを解析する。第1カメラ12A及び第2カメラ12Bは、被験者ゴルファ14によるゴルフスイングを撮像する。第1カメラ12A及び第2カメラ12Bは、ゴルフスイング軌跡を3次元計測するための映像を取得できるように、所定間隔離間して配置されている。これらの第1カメラ12A及び第2カメラ12Bは、ゴルフスイングの映像を取得した後に、第1カメラ12A及び第2カメラ12Bの間の相対的な位置関係を容易に算出可能として、ゴルフスイング軌跡の3次元計測が容易になるように冶具24で固定されている。第1カメラ12Aと第2カメラ12Bのレンズの光軸間の距離は、図2に示すように、1100mmとしている。第1カメラ12Aと第2カメラ12Bのレンズの光軸間の距離は、800~1500mmが好ましい。第1カメラ12A及び第2カメラ12Bは、レンズの光軸を、図3に示すように、カメラ同士を結ぶ線分に対して直角の方向から内側に15°回転させた状態で冶具24に固定されている。第1カメラ12A及び第2カメラ12Bのレンズの光軸を内側に回転させる角度は10°~30°であることが好ましい。また、カメラシステム15は、被験者ゴルファ14によるゴルフスイングを撮像するために、被験者14が使用するゴルフクラブが、スイング中の各段階バックスイング、ダウンスイング、フォロースルーにおいて、撮像できるように、適宜被験者14から離間して配置される。ここで、上述のように、第1カメラ12A及び12Bの間の間隔は固定されているため、被験者の身長や体格等に応じて第1カメラ12A及び12Bと、被験者14との間隔を適宜調整したとしても、同一平面上で測定したスイング軌跡を比較することができる。図6では、被験者14が打撃を行うゴルフボールの中心位置を原点とする。そして、原点を通り、水平面に対して垂直な方向の軸をz軸とする。さらに、原点をとおり、被験者14の両肩を結んだ線に対して垂直であり、かつ水平面上に延在する軸をx軸とする。さらに、原点を通り、x軸及びz軸に対して垂直な方向の軸をy軸とする。y軸方向は、ゴルフスイングによりボールを打ち出す方向に略対応するため、以下、打ち出し方向とも称する。撮影範囲は、ゴルフボールを中心にして前方0.2mから後方1.0m打ち出し方向の範囲と、ゴルフボールを中心にしてボール中心でゴルファの正面方向0.1mからゴルファの正面逆方向0.4mゴルファーの体に対して垂直方向の範囲である。また、第1カメラ12A及び12Bの撮影速度秒間撮像フレーム数は10000fpsである。図7は、図6に示す解析システムにおいて使用するゴルフクラブの一例を示す図である。図6に示す被験者14は、このゴルフクラブ以下、計測用クラブと称するを用いてスイングを行う。そして、第1カメラ12A及び第2カメラ12Bはその映像を取得する。計測用クラブは、グリップ36付近のシャフト34上に、追跡用のマーカー1が設けられ、さらに、ゴルフボールを打撃する部分であるゴルフクラブヘッド35に向かってシャフト34上に所定間隔約20cmでマーカー2~5が設けられている。追跡用のマーカーは、ゴルフクラブの軌跡を追いたい点に3点以上設けられることが好ましい。このマーカーは、例えば、白色の反射テープ、球形マーカー、反射マーカーなどを付したものであり、ゴルフスイングの撮影時には、黒背景で撮影を行うことにより、画像処理によって、マーカーの軌跡を容易に追跡することができる。図8は、解析装置の概略構成を示す機能ブロック図である。解析装置33は、第1カメラ12A及び第2カメラ12Bから、インタフェースI/F29A及び29Bを介してゴルフスイングの映像を取得し、演算を行う演算部16と、解析装置33全体を制御する制御部17と、演算部16による演算結果や、第1カメラ12A及び第2カメラ12Bから取得した映像を表示する表示部18と、第1カメラ12A及び第2カメラ12Bから取得したスイング画像データを格納する記憶部19とを備える。演算部16は、3次元座標算出部27と、マーカー変位量決定部28とを備える。3次元座標算出部27は、記憶部19から取得したスイング画像データから各マーカーの3次元座標を算出し、2つのマーカーを結ぶ直線を基準ベクトルとする。マーカー変位量決定部28は、基準ベクトルの延長線とマーカーとの距離を算出する。なお、3次元座標算出部27、マーカー変位量決定部28は、ソフトウェアプログラムを実行することによって実現される。図9は、解析装置の動作を説明するフローチャートである。また。図10は、解析装置による解析結果の一例をイメージ化した図であり、図7に示すマーカー1~5が設けられた計測用クラブを、被験者14がスイングしたときの状態を第1カメラ12A及び第2カメラ12Bで撮影したとき映像を示している。演算部16の3次元座標算出部27は、第1カメラ12A及び第2カメラ12Bがステレオ撮影した映像のスイング画像データからマーカー1~5の3次元座標を算出しS101、シャフト34の実質上撓りが生じないグリップ36近くのマーカー1とマーカー2を結ぶ線分を基準ベクトルとするS102。次に、演算部16のマーカー変位量決定部28は、マーカー3から基準ベクトルの延長線上に垂線を下ろしたときの垂線と延長線との交点を求めるS103。p3は、マーカー3から基準ベクトルの延長線上に垂線を下ろしたときの垂線と延長線との交点であり、p4は、マーカー4から基準ベクトルの延長線上に垂線を下ろしたときの垂線と延長線との交点であり、p5は、マーカー5から基準ベクトルの延長線上に垂線を下ろしたときの垂線と延長線との交点である。また、l3は、マーカー1とp3との間の距離ベクトルであり、l4は、マーカー1とp4との間の距離ベクトルであり、l5は、マーカー1とp5との間の距離ベクトルである。次に、マーカー変位量決定部22は、ベクトル計算により距離d3~d5を求めるS104。d3は、マーカー3と交点p3との間の距離ベクトルであり、d4は、マーカー4と交点p4との間の距離ベクトルであり、d5は、マーカー5と交点p5との間の距離ベクトルである。このd3~d5がゴルフクラブシャフトの変形量撓りとなる。図11は、ゴルフクラブシャフトの変形量を示す図である。縦軸は、基準ベクトル延長線上でのマーカー1からの距離lnを表しており、横軸は、基準ベクトル延長線に直交する平面上での基準ベクトル延長線からの距離絶対値dnを表している。図10には、基準ベクトル延長線上でのマーカー1~5のマーカー1からの距離lnと、マーカー1~5の、方向成分を含まない、基準ベクトル延長線までの距離dn絶対値が示されており、この図から、3次元でのゴルフクラブシャフトの変形量撓りを解析でき、剛性分布等の解析が可能となる。また、その被験者ゴルファ14に最適なゴルフクラブシャフトの選択等が可能となる。