JP4109075B2

JP4109075B2 - 球体の回転特性と飛行特性の測定方法及び球体の回転特性と飛行特性の測定装置

Info

Publication number: JP4109075B2
Application number: JP2002297580A
Authority: JP
Inventors: 勝彦植田; 正秀大貫; 光範三木
Original assignee: Ｓｒｉスポーツ株式会社; 光範三木
Priority date: 2002-10-10
Filing date: 2002-10-10
Publication date: 2008-06-25
Anticipated expiration: 2022-10-10
Also published as: JP2004132831A; US20040076311A1; US7346194B2

Description

【０００１】
【発明の属する技術分野】
本発明は、球体の回転特性と飛行特性の測定方法及び球体の回転特性と飛行特性の測定装置に関し、詳しくは、ゴルフボール等の球体の三次元姿勢及び位置を特定し、球体の回転数、回転軸方向、飛行経路や飛行速度等を精度良く測定するものである。
【０００２】
【従来の技術】
従来から、ゴルフボール等の各種球体の回転量等の測定に対して様々な測定方法および装置が提案されている。
例えば、表面に反射テープを貼り付けた球体や、表面に光を反射しない黒色に塗装した領域を有する球体に投光し、球体の回転により得られる反射光の光量変化から球体の回転量を測定する方法がある。しかし、この方法では光量のみを測定し、球体の輪郭や姿勢変位等を測定しないため、球体の回転軸の方向等を特定することはできない。よって、表面にマークを付けた回転飛行する球体を、所定の時間間隔をあけて複数枚撮影し、撮影された各画像におけるマークの変位状況から球体の回転量や回転軸角度を求めることが多い。
【０００３】
このようなマークを付した撮影画像より球体の回転量等を求める装置や方法の具体例として、特許登録第２８１０３２０号の球体の回転量測定装置、特開平１０−１８６４７４号のゴルフボールの運動測定方法、特許登録第２９５０４５０号の移動するスポーツ用物体の飛行特性を測定する装置等が挙げられる。
【０００４】
本出願人が提案している特許登録第２８１０３２０号の測定装置は、図１２（Ａ）（Ｂ）に示すように、マークＰ、Ｑを付した中心点Ｃを有する球体Ｔを二回撮影して二枚の二次元画像Ｇ１、Ｇ２を得て、各画像における球体半径を単位半径として規定すると共に、各画像毎に球体のマークＰ、Ｑ、Ｐ’、Ｑ’や中心点Ｃ、Ｃ’の三次元座標を画像上の二次元座標から算出している。これら算出された三次元座標を三次元ベクトルとして、二つの画像Ｇ１、Ｇ２間でベクトル移動量を求め、回転量および回転軸方向を算出測定している。
【０００５】
また、特開平１０−１８６４７４号の測定方法は、図１３（Ａ）に示すように、打撃時のクラブの動きを検知するセンサー２で撮影時期を判断し、第一及び第二カメラ１Ａ、１Ｂで時間をあけてボールＢ１を撮影することにより図１３（Ｂ）に示す二個のボールＢ１、Ｂ１’が写ったボール画像Ｇ３を得ている。この二次元のボール画像Ｇ３を上記と同様に処理し、ボールの回転量等を算出測定している。
【０００６】
さらに、図１４（Ａ）の特許登録第２９５０４５０号の測定装置４は、２台の同期したカメラ５Ａ、５Ｂにより、一枚の画像にマークＢａを付したボールＢ２、Ｂ２’が二個写った画像を夫々撮影し、カメラ５Ａ、５Ｂ間の視野の関係と関連付けて三角測量に近似した原理で三次元座標を導いている。これにより図１４（Ｂ）に示すボールＢ２等の三次元領域図を得て、ボールの各種特性を測定している。なお、このようにして三次元座標を導く手法は、ＤＬＴ（Direct Liner Transformation）法として知られている。
【０００７】
さらには、特許第３１８５８５０号では、スポーツ用物体の飛行特性を測定表示するモニター装置が提案されている。具体的には、シャッター手段等を用い、１枚の画像を少なくとも１台のカメラを用いて撮影しており、少なくとも１台のカメラは予め対象となる空間に置かれた既知の座標を持った所定数以上の点を用いて校正される必要がある。既知の三次元空間座標とカメラのフィルム面に投影した二次元座標との関係を利用して、６変数（球体の質量重心座標３成分と基準座標からの回転量３つ）をテーラーの定理の線形化を８度の反復計算して解くことができ、ボール画像の半径等を基準とせずに、実空間での長さを基準としてボール上固定された点（マーク）の三次元座標を得ることができる。従って、ボール画像の輪郭を正確に撮影する必要がなくなり、高速で移動するボールの静止画を得るための高速シャッター手段による明るさ不足等の問題が緩和される等、撮影設備の自由度を向上させたものである。
【０００８】
また、本出願人は、特許第２６２６９６４号で、複数の平行な光線を投光する投光装置と、その装置からの光線の入光・遮光を検知する複数の受光装置を有し、インパクト後から既知である投光・受光装置が存在する一平面までの距離・遮光された光線の位置、光線間の寸法、及びボールの径を用いて、光線の遮光時間を計測することによって、ボール中心の高さ位置や振れ位置を決定することができ、打出角や振れ角を測定することができる球状物体の位置測定及び打出角測定方法を提案している。
【０００９】
【特許文献１】
特許第２８１０３２０号公報
【００１０】
【特許文献２】
特開平１０−１８６４７４号公報
【００１１】
【特許文献３】
特許第２９５０４５０号公報
【００１２】
【特許文献４】
特許第３１８５８５０号公報
【００１３】
【特許文献５】
特許第２６２６９６４号公報
【００１４】
【発明が解決しようとする課題】
しかしながら、図１２（Ａ）（Ｂ）の測定装置および図１３（Ａ）（Ｂ）の測定方法では、測定にかかる演算にボール画像の半径を用いるので、画像から得られる半径寸法の精度が、算出される三次元ベクトルの精度を左右してしまい、測定の基となる画像を高精度で撮影した上に、撮影画像からも高精度で球体半径を求める必要がある。しかし、高速で飛行するボール画像の静止画像を得るには、例えば、高速度シャッターを備える高速度カメラを使用することになるが、高速度シャッターはシャッターが開いている時間が僅かなので充分な光量を得にくい。
【００１５】
よって、撮影された画像は、カメラの対向位置となるボール中心付近は、比較的鮮明な画像が得られるものの、球体であるボールの輪郭を鮮明に捉えることは困難であり、このことはライティング等を工夫しても改善されにくい。その結果撮影された画像のボール輪郭は不鮮明となり、ボール画像より読みとられる半径の寸法精度等も悪化し、ボールの回転量等に対する測定精度も自然と低下せざるを得ない問題がある。
【００１６】
一方、図１４（Ａ）（Ｂ）に示す測定は、ボール画像の半径を用いることなく、実空間の長さを基準にしてボール表面のマークに対する三次元座標を得るため、ボール輪郭を明確に撮影する必要もなく光度不足の問題も比較的生じず、測定設備の負担も軽減できる利点がある。しかし、ボール表面に付されたマークの三次元座標を高精度で求めるには、各マークを正確に読みとれるように画像においてボールを比較的大きく撮影する必要がある。このようにボールを大きく写そうとすると、二つのボール画像の撮影する時間間隔を短くする必要があり、自然と二つのボール画像間の回転量は小さくなる。
【００１７】
一方、ボールの回転量の測定精度を高めるには、マークの移動距離を大きくしてマーク位置の変位を大きくすること、即ち、撮影された２つのボール画像間でボールが回転した量を大きくする必要があり、上記ボール画像を大きくすることに相反する条件が必要となる。
【００１８】
よって、ボール画像を大きくすれば、マークの三次元座標は高精度で測定できても、各画像間の位置変化が小さいためボールの回転量の測定に対しては精度が落ちてしまう問題がある。一方、ボールの回転量を大きくなるように撮影すると、時間間隔をあけて２つのボール画像を撮影しなければならず、ボールの回転量に対する測定精度は上昇するが、ボール画像は小さくなるので、マークの三次元座標の測定精度が落ちてしまい、いずれにしても、マークの三次元座標およびボールの回転量の双方を両立して高精度で測定できない問題がある。
【００１９】
なお、上記問題に対しては、測定に要する装置を二セット用意して、一セット目の装置で一方のボール画像を撮影すると共に、時間間隔をあけて二セット目の装置で他方のボール画像を撮影し、マークの三次元座標とボール回転量の両方を高精度で測定することも想定し得る。しかし、測定においては両セット間の計四台の高速カメラを関連させて校正（キャリブレーション）する必要がある上に、装置構成も非常に複雑になりコストも大幅に上昇してしまい、装置を二セット使用するの現実的には困難である。
【００２０】
さらに、撮影されたボール表面のマークの移動量からボールの回転量等を算出測定する際は、一つ目のボール画像における特定のマークが、二つ目のボール画像において何れのマークに対応するかを見きわめる必要がある。しかし、ボールの回転軸方向が予想でき、かつ、画像間のボールの回転変化が少ない場合は、上記見きわめは比較的容易であるが、回転軸の方向が測定毎に大きく変化するなどの理由で、回転軸方向が予想できない場合やボールの回転量が大きい場合等は、上記見きわめが非常に困難となり、コンピュータ等の自動認識等ではボールの回転量等の測定も不可能になるおそれがある。また、人間が上記見きわめを行う場合であっても時間を要したり、見きわめを誤るおそれもある。
【００２１】
その上、一つ目のボール画像では写っていたマークが、二つ目のボール画像では、ボールの回転により裏面側に回ってしまい画像上のボール表面に表れない場合には、測定自体が不可能となることより、カメラの測定方向や測定対象のボールの回転方向等にも制限が生じ、自由に最適な状況で測定できにくい問題もある。
【００２２】
さらには、特許第３１８５８５０号では、１つ目のボール画像上に見えていたマークが、ボールの回転によって、２つ目のボール画像上に写らなくなる。即ち、ボールの裏側に１つ目のボール画像上のマークが回ってしまう場合には、計測不能となるという問題がある。よって、カメラでの撮影方向や測定するボールの回転方向にも制約が生じ可搬性に欠けるという問題がある。
【００２３】
また、特許第２６２６９６４号は、飛行中の球状物体の通過位置を非接触で広範囲かつ高精度に測定できるものの、球状物体のスピン量等を算出するには未だ改善の余地がある。
【００２４】
本発明は、上記した種々の問題に鑑みてなされたものであり、コンピュータにより飛行する回転球体の回転数や回転軸方向等の回転特性、及び、飛行経路や飛行速度等の飛行特性を精度良く測定することを課題としている。
【００２５】
【課題を解決するための手段】
上記課題を解決するため、本発明は、表面に複数のマークを付した球体の回転飛行時に所定の時間間隔をあけて複数回撮影して、該球体の複数の二次元映像を得る一方、該球体と同様に表面にマークを付した仮想球体をコンピュータの三次元空間座標上に作成し、該仮想球体の任意姿勢及び任意位置を仮想球体の基準姿勢及び基準位置として設定すると共に、
少なくとも１台の撮影手段を用いて、三次元座標と二次元座標の間の関係を予め導出しておき、
上記関係を用いて、上記仮想球体の三次元空間座標上のマーク位置を、二次元的な画像上に変換することにより二次元の仮想マーク座標値を求めると共に、上記映像上に存在する球体の二次元の映像マーク座標値を求め、
上記仮想マーク座標値と、上記映像マーク座標値とが一致するように、上記仮想球体を、上記基準姿勢及び基準位置に対して姿勢変位操作し、該姿勢変位操作の量により、球体の各二次元映像毎に、球体の三次元姿勢及び位置を特定し、
上記二次元映像毎に特定された、ある時間の三次元姿勢及び位置と、他の時間の三次元姿勢及び位置とにより、球体の回転特性と飛行特性を演算していることを特徴とする球体の回転特性と飛行特性の測定方法を提供している。
【００２６】
このように、本発明は、上記のような方法で球体の三次元姿勢及び位置を特定しているため、ある時間の三次元姿勢及び位置と、他の時間の三次元姿勢及び位置との関係から球体の回転数と回転軸方向等の回転特性、飛行経路や飛行速度等の飛行特性を容易かつ高精度で測定することができる。よって、外力等を受けずに回転している球体であるゴルフボールのスピン量や、打球の打出角や振れ角の測定に最適であり、特に打球時の左右の振れ角の測定精度を向上することができる。
【００２７】
具体的には、少なくとも１台の撮影手段を用いて、三次元座標と二次元座標の間の関係を予め導出して把握している。このため、コンピュータ上に作成した仮想球体上のマークの位置関係から、２次元的な画像上のマークの位置関係を求めることができる。よって、上記関係を用いて、上記仮想球体の三次元空間座標上のマーク位置を、二次元的な画像上に変換し、二次元の仮想マーク座標値を求めることができる。よって、従来にはできなかった二次元の仮想マークと三次元の仮想マークの対応が可能となり、二次元の仮想球体でも三次元の情報を持たせることができ、奥行き方向の情報を把握することができる。従来の方法では、奥行き方向に関する考慮がなされておらず単に投影しているだけであり、全て二次元画像は平面にあった。また、従来は投影する方向は一方向であったが、本発明は、仮想球体の位置によって投影方向が異なっている。
よって、球体のエッジがはっきり見えていない場合や、球体映像が画面の端部に撮影されたり、小さく撮影されたような場合でも、精度良く測定することができる。
【００２８】
また、上記のように仮想マーク座標値と、映像マーク座標値とにより、仮想球体と実際の球体との姿勢及び位置を関連付けることで、実際の球体の姿勢や位置を、仮想三次元空間における基準姿勢及び基準位置に対する相対座標値で表現することが可能となる。これに伴い、球体の回転特性や飛行特性の解析にも上記仮想空間の基準姿勢及び基準位置に対する相対座標値を用いることが可能となり、球体の回転特性や飛行特性の解析をコンピュータ上で容易に自動演算することができる。なお、基準位置とは、三次元空間と二次元画像との関係を求めた時に定義した三次元座標系の原点を意味する。
【００２９】
また、本発明では、球体の輪郭データを用いず、撮影された各二次元映像の球体の複数のマークから球体の姿勢及び位置を特定している。このため、各映像において球体に付されたマークさえ明確に確認できれば、球体の輪郭が不明瞭であっても球体の姿勢及び位置等に関する測定を高精度で行うことができる。なお、撮影された二次元映像より球体のマークをコンピュータに自動認識させるためには、撮影映像に対して白と黒のみで表示する二値化処理を行うことが考えられる。
【００３０】
また、ある時間の三次元姿勢及び位置を他の時間の三次元姿勢及び位置に一致させる際の回転操作に関する回転行列を算出することにより、球体の回転特性を求めている。上記球体の回転量等は、時間間隔が特定可能な、ある時間の球体が有する回転行列に、他の時間の球体が有する回転行列の逆行列を掛けて算出している。このように算出して得られたある時間の球体の姿勢から、他の時間の球体の姿勢への回転行列の数値を判断すれば、回転量も判断でき、また、回転軸のベクトル即ち回転軸方向も容易に求めることができる。
【００３１】
さらに、ある時間の三次元姿勢及び位置を他の時間の三次元姿勢及び位置に一致させる際の球体の重心座標の移動量と移動方向を算出することにより、球体の飛行特性を求めている。
【００３２】
上記三次元座標と二次元座標の間の関係を予め導出する際に、少なくとも６つ以上の三次元座標を用いていることが好ましい。
これにより、三次元空間と二次元空間の関係を求めることができる。即ち、実空間とカメラ上の画像との関係を明確にすることができる。なお、６つより少ないと変数に対して既知数の量が少ないために三次元と二次元関係を示す連立方程式を解くことができないためである。
【００３３】
仮想球体の三次元空間座標上のマーク位置を、二次元的な画像上に変換することにより二次元の仮想マーク座標値を求める方法としては、コンピュータ上に構築した基準姿勢及び基準位置の球体のマークを、既知の実空間三次元座標と二次元的な画像上の座標の関係を利用して、二次元の画像上に投影すること等が挙げられる。
【００３４】
上記姿勢変位操作における設計変数として、上記球体の三次元重心位置及び仮想三次元空間の直交座標軸である横軸、縦軸、垂直軸の各軸回りの球体の回転角度の計６個を用いていることが好ましい。
これにより、球体の奥行き情報を算出することができるので、回転特性と飛行特性を同時に算出することができる。
【００３５】
上記姿勢変位操作は、上記仮想球体の移動操作、回転操作であり、上記基準姿勢及び基準位置に対する姿勢変位操作の量を遺伝的アルゴリズムという最適化手法に応じた演算で求めていることが好ましい。
【００３６】
コンピュータ上に形成した仮想球体は基準姿勢及び基準位置を任意に設定するため、この基準姿勢及び基準位置は、撮影された二次元映像の球体と、大きさ、位置、および、姿勢が相異している。よって、球体の姿勢及び位置を仮想空間における仮想球体の基準姿勢及び三次元と二次元の関係式を求めた三次元空間での基準位置に対する相対座標値で表示可能にするため、仮想球体が二次元映像の球体の姿勢及び位置に一致するように姿勢を変位する操作が必要となる。
【００３７】
この姿勢変位は、コンピュータ上で仮想球体を平行移動等の移動、回転操作により行うが、これら平行移動等における変位操作の量を適正に決定することにより、仮想球体を実際の球体の姿勢及び位置へと正確に一致させることができ、それに伴い、球体の姿勢及び位置を三次元と二次元の関係式を求めた三次元空間での基準姿勢及び基準位置に対する相対座標値により高精度で特定できる。本発明では、上記変位操作量の決定を最適化の問題として捉え、最適化問題の解決手法の一つである遺伝的アルゴリズムという計算機用アルゴリズムを用いて演算を行うのが好ましい。
【００３８】
上記遺伝的アルゴリズムによる演算解法は、シミュレーティド・アニーリング法や傾斜法、他の線形計画法といった他の最適化手法と異なる特徴を有しており、例えば、演算の対象となる変数を二進数、文字列、ベクトル等にコーディング（コード化）することや、演算結果を評価するための適合度を、目的関数を基に設定すること等の特徴を有している。
【００３９】
遺伝的アルゴリズムは、多数の個体から最適な解を演算により同時探索する確率的な解法であり、演算された解を最適化されたものにすることができる。よって、コンピュータ上の仮想球体を、球体上の複数のマークを参照しながら二次元映像の球体の姿勢及び位置と一致させるのに有用である。
【００４０】
上記遺伝的アルゴリズムに応じた演算では、上記移動操作する量を上記仮想球体の三次元的な移動操作に関する三つの独立した変数として、上記回転操作する量を上記仮想球体の三次元的な回転操作に関する三つの独立した変数として、それぞれ設定し、これら六変数の各数値に基づき上記各操作を行っている。このように定められた六変数が基準姿勢及び三次元と二次元の関係式を求めた三次元空間での基準位置に対する相対座標値に相当し、よって、これら六変数の数値を求めれば球体の姿勢及び位置を特定することができる。
【００４１】
なお、上記移動操作に関する三つの変数としては球体の重心位置（（ｘ，ｙ，ｚ）の３つ）が挙げられる。さらに、上記三次元的な回転操作に関する三つの独立変数としては、仮想三次元空間の直交座標軸である横軸、縦軸、垂直軸の各軸回りの回転角度が挙げられる。
【００４２】
上記球体表面の複数のマークの位置関係は、回転操作に対して４回以下の回転対称性を有していることが好ましい。
即ち、上記球体表面の複数のマークは、３６０度を４以下の自然数で割った角度のみで球体中心を通る軸を中心に球体を回転させた場合、回転の前後で一定の方向から見た球体表面の全てのマークが一致する位置に付けている。これは、ゴルフボールやテニスボール等の球技用のボールの回転速度等の測定に特に有用である。
【００４３】
このようにして球体表面に複数のマークを付加すると、９０度以下の角度で球体を回転させた際に、回転の前後で球体表面上に見える全てのマークが一致することがないため、球体のとり得る姿勢の数を少数に限定できる。なお、上記所要角度の回転による前後で一定方向からの球体表面の全マークが一致することを、以下、（回転）対称と呼ぶ。
【００４４】
具体的には、３６０度、１８０度、１２０度、９０度の各角度毎の回転になる。例えば、３６０度回転で対称とは、上記球体表面の複数のマークがランダムに付けられている場合であり、一回転で一度しか対称にならず球体の姿勢を一義的に特定でき、演算の精度も最も高くできる。１８０度の場合は、一回転で二度対称となる位置が存在するため、球体の回転姿勢を特定する場合、二通りの姿勢が想定され一義的に姿勢を特定できない反面、同じパターンを２度繰り返すことができ、マークを付ける際の効率が高い。以降、１２０度の場合は一回転で三度対称となる位置が存在し、９０度の場合は一回転で４度対称となる位置が存在することになり、球体にマークを付ける負担も段階毎に軽減される。
【００４５】
また、上述したようにマークを付しても、球体表面の一部分にマークが偏っている場合には、マークが測定視方向の裏側に回り、球体の姿勢が特定できない場合も想定されるため、上記球体表面には、表面全体にわたり１０個以上１００個以下のマークを付けている。上記範囲としているのは、、１０個より少なければ、観察方向から読みとれるマークが少数になりすぎて、姿勢を特定できなくなるおそれが生じるからであり、一方、１００個より多いと演算にかかる時間が長くなり、効率が悪化するためである。
【００４６】
なお、測定対象の球体としては、ゴルフボール、テニスボール等の各種球技用のボールが挙げられる。これらの球技用のボールの使用時の回転特性や飛行特性等の解析に有用である。
【００４７】
また、ボール映像の輪郭を明確に撮影することは根本的に困難であることから、ボール映像の輪郭のデータを用いずに回転特性を測定することが好ましい。特に、撮影映像の端部でボールが撮影された場合、ボールの右あるいは左部がストロボの影響で陰になってしまうために、できるだけ輪郭部分を用いずに回転特性を測定することが好ましい。
【００４８】
また、本発明は、表面に複数のマークを付した球体が回転しながら飛行する時に、該球体を撮影できる撮影手段と、
上記撮影により得られた球体の二次元映像を記録する記録手段と、
上記球体と同様のマークを付した仮想球体を三次元座標空間上に作成し、該仮想球体と上記球体の二次元映像に基づいて、球体の三次元姿勢及び位置を特定し、球体の回転特性と飛行特性を求める演算手段とを備え、
上記演算手段は、少なくとも１台の撮影手段を用いて、三次元座標と二次元座標の間の関係を予め導出可能な座標変換プログラムを有し、
かつ、上記二次元映像上の球体の二次元の映像マーク座標値が、上記仮想球体の三次元座標上のマーク位置を変換することにより求められる二次元の仮想マーク座標値と一致するように、仮想球体を姿勢変位操作し、上記仮想球体の基準姿勢及び基準位置に対する姿勢変位操作の量により、上記球体の三次元姿勢及び位置を特定する姿勢認識プログラムを有してることを特徴とする球体の回転特性と飛行特性の測定装置を提供している。
【００４９】
本発明の球体の回転特性と飛行特性の測定装置は、上記撮影手段と、上記記録手段と、上記座標変換プログラム及び姿勢認識プログラムを有した演算手段を備えているため、上述した本発明の球体の回転特性と飛行特性の測定方法を実行するのに最適であり、球体の姿勢、回転特性、飛行特性を精度良く求めることができる。
【００５０】
上記演算手段は、上記基準姿勢及び基準位置に対する姿勢変位操作の量を遺伝的アルゴリズムに応じて演算する最適化プログラムを有していることが好ましい。このように、演算手段において、遺伝的アルゴリズムに応じて演算する最適化プログラムを有していることで、上述した遺伝的アルゴリズムを用いた上記測定方法による演算を容易に行うことができ、演算精度を高めることができる。
【００５１】
上記撮影手段は、回転や移動している球体を所定の時間間隔をあけて複数回撮影可能な構成としていることが好ましい。
具体的には、実物の球体の二次元映像を撮影するためには、スチールカメラ、ＣＣＤカメラ等の映像を取り込むための撮影手段とフィルムや画像メモリ等の映像を記録するための記録手段が必要である。その他、短時間の発光の輝度差を利用して瞬間の静止状態のボール画像を得るためにマイクロフラッシュも用いられる。なお、１０ｍ／ｓ以上の高速で移動する球体が測定対象であっても、高速シャッターやマイクロフラッシュを用いることで、ブレのない静止画像を得ることができる。
【００５２】
上記撮影手段としては２台のカメラを用いることが好ましい。その間隔をキャリブレーション無しに容易に可変とすることが好ましい。２台のカメラの間隔や撮影時間間隔を比較的長くすることで角度変化を大きくし、回転数の測定精度を高めることができる。また、１つ目のボール映像上の特定のマークが、２つ目のボール映像上のどのマーク位置に対応するのかを認識しやすくし、自動化を容易にすることができる。なお、二次元映像の撮影手段は、一台でも複数台でも適応可能である。
【００５３】
また、複数の映像を得るには、少なくとも二回球体の映像を取り込む必要があるが、一台のカメラで複数回マイクロフラッシュを発光させるなどして、一枚の映像に複数の球体を写してもよく、また、一台のカメラで複数枚の映像にそれぞれ一つずつ球体を写すようにしてもよい。さらに、複数のカメラを用いて所要の時間間隔をあけて各カメラ毎に二次元映像を撮影してもよい。ただし、複数のカメラで撮影した場合は、各映像に対する共通の座標軸を校正等により特定することが必要である。
【００５４】
また、撮影された映像より球体に付されたマークの位置を得る手段として、記録手段にフィルムを用いた場合では、現像した写真に直接定規等を当てて測定する原始的な方法や、フィルムや写真の画像をコンピュータに取り込めるスキャナー等の画像入力装置を用いてコンピュータ画面上で測定用カーソルを合わせて測定する方法もある。さらに、記録手段に画像メモリを用いた場合は、画像メモリ内のデータをコンピュータ画面上に読み出し、上記同様測定用カーソルで測定する方法等がある。即ち、球体映像における各マークの二次元位置さえ得られれば、特に測定手段は限定されず上記のいずれの方法を用いてもよい。
【００５５】
【発明の実施の形態】
以下、本発明の実施形態を図面を参照して説明する。
図１（Ａ）は、本発明の球体の回転特性と飛行特性の測定方法を実行する球体の回転特性と飛行特性の測定装置１０の構成を示す。
上記測定装置１０は、球体を各方向から撮影できる撮影手段と、撮影により得られた球体の二次元映像を記録する記録手段と、球体と同様の仮想球体を三次元座標空間上に作成し、該仮想球体と上記球体の二次元画像に基づいて、球体の三次元姿勢及び位置を特定し、球体の回転特性と飛行特性を求める演算手段とを備えている。
【００５６】
また、演算手段は、少なくとも１台の撮影手段を用いて、三次元座標と二次元座標の間の関係を予め導出可能な座標変換プログラムを有し、かつ、二次元映像上の球体の二次元の映像マーク座標値が、仮想球体の三次元座標上のマーク位置を変換することにより求められる二次元の仮想マーク座標値と一致するように、仮想球体を姿勢変位操作し、仮想球体の基準姿勢及び基準位置に対する姿勢変位操作の量により、球体の三次元姿勢及び位置を特定する姿勢認識プログラムを有している。
【００５７】
図１（Ｂ）に示すように、測定装置１０は、撮影手段としてＣＣＤカメラ１１を使用すると共に二台のマイクロフラッシュ１２−１、１２−２も併用している。測定対象となる球体をＣＣＤカメラ１１により撮影し、順次得られる映像を記録手段である画像メモリ１５ｂに蓄積する構成としている。画像メモリ１５ｂは、演算手段であるコンピュータ１５中に設置されている。
【００５８】
映像撮影時に、所要の時間間隔をあけてマイクロフラッシュ１２−１、１２−２を一台づつ発光させてから画像メモリ１５ｂへの蓄積を終了させることで、１枚の二次元映像に球体が二個写った画像を画像メモリ１５ｂに得ている。なお、マイクロフラッシュ１２−１、１２−２を発光させるタイミングには、光電管スイッチ１６−１、１６−２とリターダー１７によるトリガー信号を利用している。
【００５９】
また、測定装置１０に含まれるコンピュータ１５は、上述した画像メモリ１５ｂ以外にも中央演算処理装置（ＣＰＵ）、メモリ、記憶装置であるハードディスク１５ａ等を備え、ハードディスク１５ａは、上記座標変換プログラム、上記姿勢認識プログラムや、基準姿勢及び基準位置に対する姿勢変位操作の量を遺伝的アルゴリズムに応じて演算する最適化プログラムを有している。
【００６０】
上記姿勢認識プログラム、座標変換プログラム、最適化プログラムは、基本的には図２のフローチャートに基づいて組まれている。
即ち、測定対象の球体を基に仮想三次元空間内に仮想球体を形成し、仮想球体に実物の測定対称と同様にマークを付け、この形成された仮想球体の任意の姿勢及び位置を基準姿勢及び基準位置として設定する。
少なくとも１台の撮影手段を用いて予め導出した三次元座標と二次元座標の間の関係を用いて、仮想球体の三次元空間座標上のマーク位置を、二次元的な画像上に変換することにより二次元の仮想マーク座標値を求める。また、上記映像上に存在する球体の二次元の映像マーク座標値を求める。
仮想マーク座標値と、映像マーク座標値とが一致するように、仮想球体を、基準姿勢及び基準位置に対して姿勢変位操作（移動、回転等）し、姿勢変位操作の量により、球体の各二次元映像毎に、球体の三次元姿勢及び位置を特定する。
二次元映像毎に特定された、ある時間の三次元姿勢及び位置と、他の時間の三次元姿勢及び位置とにより、球体の回転特性と飛行特性を演算する。
なお、上述した処理以外には撮影された二次元画像を白と黒の表示に変換する二値化処理も含まれている。
【００６１】
上記のような仮想球体の姿勢変位操作に伴う変位量は、遺伝的アルゴリズムによる演算を用いて特定するように、遺伝的アルゴリズムに応じて演算する最適化プログラムにプログラミングされている。三次元姿勢の特定は、球体の回転特性と飛行特性を求めるに当たり撮影された球体映像毎に行っている。
【００６２】
以下、球体の回転特性と飛行特性の測定装置を用いた球体の回転特性と飛行特性の測定方法について更に詳しく説明する。
【００６３】
先ず、測定対象の球体としてゴルフボールを使用し、図３（Ａ）（Ｂ）（Ｃ）に示すようにゴルフボール２０の表面に黒色のマーク２１を多数付けている。
【００６４】
この場合、３６０度を１で割った角度で、ゴルフボール２０を回転軸Ｊを中心に回転させた時、回転の前後で一度だけ対称になるようにマーク２１を付けている。即ち、回転操作に対して１回の回転対称性を有しており、ゴルフボール２０の全表面でランダムにマーク２１を付すことになる。
【００６５】
また、上記のようにマーク２１を付したゴルフボール２０に対し三次元測定器等を利用して、各マークの三次元における座標値を読みとる。この際、読み取りに用いる座標は、図４に示すように極座標（動径方向座標値ｒ、子午線方向角度θ、方位角方向角度φ）を採用している。
【００６６】
次に、図５に示すように、コンピュータ１５内の姿勢認識プログラムにより、コンピュータ上の仮想三次元座標空間Ｓに、実際のゴルフボール２０と同様の仮想ゴルフボール３０を形成し、この仮想ゴルフボール３０の表面にもゴルフボール２０と同様となるように、上記読みとった極座標（ｒ、θ、φ）の数値の位置にマーク３１を付している。なお、上述した姿勢認識プログラムでは直交座標に変換して演算を行っている。また、このように形成された仮想ゴルフボール３０の任意の姿勢を基準姿勢として設定する。
【００６７】
次に、少なくとも１台の撮影手段を用いて、三次元座標と二次元座標の間の関係を予め導出して把握し、仮想球体のマーク位置からカメラのフィルム面（二次元）のマーク位置関係を求める方法について説明する。即ち、空間座標からカメラの座標に変換する。
【００６８】
図６に三次元実空間上の座標と二次元的な画像上の座標の位置関係を示す。以下に変数の詳細な説明を示す。
実空間上に任意の座標軸を設定し、点Ｐ（Ｘ，Ｙ，Ｚ）、レンズ中心の座標Ｏ（Ｘ_０，Ｙ_０，Ｚ_０）とする。
点Ｐを含み光軸に垂直な面とレンズの中心との距離をＬ、ファイル面とレンズ中心との距離をＦとする。
フィルム面上にも任意の座標軸を設定し、点Ｐの像であるＱのフィルム面上での座標が（Ｕ，Ｖ）であるとし、光軸とフィルム面との交点の座標を（Ｕ_０，Ｖ_０）とする。
レンズ中心Ｏを原点とし、Ｚ軸が光軸と平行で、かつＸ・Ｙ軸がＵ・Ｖ軸と平行である座標系（Ｘ’Ｙ’Ｚ’）を考えた場合、実空間座標（ＸＹＺ座標）からＸ’Ｙ’Ｚ’への変換行列をＭ（ｍ_ｉｊ：ｉ＝１〜３、ｊ＝１〜３）とする。ただし、この変換行列はオイラーの定理に従うとするので、Ｍは下記の数式１のようになる。
【００６９】
【数１】

【００７０】
実空間上の点Ｐとその像Ｑとの関係を直接表すためには、両者の関係を同じ座標系で表さなければならない。そこで、上記で定義した実空間から光軸をＺ軸とし、画像軸Ｕ及びＶ軸に平行であるＸ及びＹ軸をもつ座標系への変換行列Ｍを用いて、数式２が得られる。
【００７１】
【数２】

【００７２】
ＯＰとＯＱは互いに同一直線上にあり、その長さの比がＬ：Ｆであることから、数式３が得られる。
【００７３】
【数３】

【００７４】
数式２を数式３へ代入することによって、３式のうちの１式からＬを求め、そのＬを他の２式に代入することにより、以下の数式４が得られる。
【００７５】
【数４】

【００７６】
この数式４における変数をまとめ、Ｕ，Ｖ，Ｘ，Ｙ，Ｚについて整理すれば、数式５が得られる。
【００７７】
【数５】

【００７８】
１１個の係数を求めるためには、カメラの台数によるが、カメラが２台の場合、６組の（Ｘ，Ｙ，Ｚ）と（Ｕ，Ｖ）が得られれば数式５に代入することによって算出できるが、測定誤差の影響を抑えるために、６組以上の既知の実空間座標を測定することが望ましい。また、カメラ１台では１２ヶ以上の（Ｘ，Ｙ，Ｚ）と（Ｕ，Ｖ）が必要で、３台では４ヶ以上となる。
【００７９】
以上の方法によって、既知の実空間座標と二次元的な画像上の座標の関係を把握することができる。この関係を利用して、図７に示すように、コンピュータ上の仮想的な球体のマークの座標をカメラのフィルム面上の二次元的な座標を求めることができる。即ち、仮想ゴルフボール３０の三次元空間座標上のマーク３０の位置を、二次元的な画像上に変換することができる。
【００８０】
次に、仮想的な三次元球体上のマークを二次元座標へ投影した場合、フィルム画像上にマークが存在するかしないか（マークが表に存在するか裏に存在するか）を判定する方法について説明する。
【００８１】
図２に示すフローチャートで示したマークが表あるいは裏のどちらに存在するのかを判定する方法を以下に示す。上記数式５で定義した１１個の係数を求めることができ、それらの係数を利用してカメラレンズの三次元座標を計算する。
【００８２】
図８に示すように、レンズ中心Ｏ（Ｘ_０，Ｙ_０，Ｚ_０）とし、仮想球体重心Ｐ（Ｇ_Ｘ，Ｇ_Ｙ，Ｇ_Ｚ）とすると、ベクトルＯＰは、数式６のようになる。
【００８３】
【数６】

【００８４】
ベクトルＯＰと球の径Ｒを用いて、レンズ中心から球体への接線長さｌは、数式７で表される。
【００８５】
【数７】

【００８６】
仮想的な球体の斜線部分を二次元に投影した場合、実際に撮影できるマーク（表側）と定義し、それ以外のマークを裏側と定義する。レンズ中心と仮想的な球体上のマークの長さをＬとすると、Ｌは仮想的な球体の重心位置と既知である各マークの位置関係により決まる。
【００８７】
Ｌ＜ｌ：球体上のマークは実際に測定できる画像の方向（表側）であり、ＧＡ計算に利用するマーク。
Ｌ≧ｌ：球体上のマークは実際に測定できない画像の方向（裏側）であり、ＧＡ計算に利用しないマーク。
【００８８】
以上の方法により、仮想的な球体上の全てのマークでＧＡの計算に用いることができるマークであるか否かを判定することができる。しかし、レンズの中心位置Ｏ（Ｘ_０，Ｙ_０，Ｚ_０）を求める必要がある。そこで、上記数式４と数式５を用いてレンズの中心位置を決定する。これらの式を比較すると、以下のような関係式（数式８）が成り立つ。
【００８９】
【数８】

【００９０】
数式８は変数が１２個存在するが、１１式しかないので、このままでは、Ｘ_０，Ｙ_０，Ｚ_０を計算できない。そこで、変換行列Ｍの行列式を計算し、数式９の関係式を追加することによって、１４式から１４個の変数を求めることができ、Ｘ_０，Ｙ_０，Ｚ_０を求めることができる。
【００９１】
【数９】

【００９２】
次に、図９（Ａ）（Ｂ）に示すように、回転装置１９により、実物のゴルフボール２０が回転する状態を撮影できる位置にＣＣＤカメラ１１を配置しており、ゴルフボール２０との距離やカメラレンズの倍率等を適宜調整している。
【００９３】
上記配置状態で、ＣＣＤカメラ１１の前をゴルフボール２０が通過する際に、マイクロフラッシュを所要の時間間隔で２回発光させ、一枚の二次元映像にゴルフボール２０が二個写った静止映像を画像メモリ１１ｂに得ている。ここで得られた二次元映像を図１０に示す。このように撮影した映像データを二値化処理用のプログラムで白と黒に二値化してから、ディスプレイ画面１８上で、各マーク２１毎に二次元の座標値を読み取ると共に、読みとった値をハードディスク１５ａに記録している。
【００９４】
姿勢認識プログラムにより、上記方法で予め導出した三次元座標と二次元座標の間の関係を用いて、仮想ゴルフボール３０の三次元空間座標上のマーク位置を、二次元的な画像上に変換することで求めた二次元の仮想マーク座標値と、上記撮影により得られた映像上に存在する実物のゴルフボール２０の二次元の映像マーク座標値とが一致するように、各々移動操作、回転操作を行っている。
これにより、仮想ゴルフボール３０の基準姿勢及び基準位置に対する姿勢及び位置を特定している。この移動操作、回転操作を行う量を、最適化手法の一種である遺伝的アルゴリズムを用いて演算している。
【００９５】
以下に、本実施形態の遺伝的アルゴリズム（ＧＡ）を用いた最適化の条件を示す。
任意に設定した三次元直行座標の各軸周りの回転角度（α、β、γ）、三次元空間での仮想ゴルフボール３０の重心座標（Ｘ_０、Ｙ_０、Ｚ_０）の計６項目を設計変数としている。この６項目の数値を六変数としている。ＧＡによる最適化では、設計変数をコード化しておく必要があるが、α、β、γ、Ｘ_０、Ｙ_０、Ｚ_０の各変数（実数値）を１０ビットの２進法に変換したものをコードとして使用している。
【００９６】
また、本実施形態のＧＡによる最適化の基本設定値を以下に示す。
個体数：５０
染色体長：６０ビット
設計変数の数：６変数（１変数あたり１０ビット）
終了条件（適合度の最大値が１０００世代にわたり変化しなければ終了）
マークの数：４２
【００９７】
コンピュータにより上記操作を行って得られた二次元投影図と、撮影した二次元投影図をパターンマッチングする。この時のパターンの類似度を下記の数式１０で定義し、これを目的関数として、その最大値を遺伝的アルゴリズムによって探索する。
【００９８】
【数１０】

【００９９】
上記数式１０において、Ｘｉは撮影した投影図のマークの二次元座標（原点は球体（ボール）の中心座標）、Ｙｉは遺伝的アルゴリズムで得た値によって回転させた後のボールの三次元空間での重心位置から、既知の実空間座標と二次元的な画像上の座標の関係を利用して求めた投影図のマークの二次元座標（原点は球体の中心座標）である。
また、distance(ａ,ｂ)は２つの点ａ、ｂの距離を示し、ｍｉｎ_ｊはｊ（ｊ＝１からｎ、ｎはマークの個数）に関して最小値を表し、Σはｉ（ｉ＝１からｎ、ｎはマークの個数）に関して和をとることを意味する。ただし、Ｎはマーク数である。それぞれの画像について、マーク間距離の最小値を全てのマークについて計算し、その総和を適合度とする。
【０１００】
上記目的関数の演算は揃えられた各個体全てに対して行われ、これら演算より各個体毎の適合度が求められる。
このような適合度計算を何世代にもわたりコンピュータ１５で繰り返し、適合度の最大値が１０００世代にわたり変化しなければ、その時点で演算を終了している。この適合度に最大値をもたらす個体が有している六変数の値が、仮想ゴルフボールの姿勢変位の最適解であると判断し、この個体によりゴルフボール２０の三次元姿勢及び位置を特定している。なお、遺伝的アルゴリズムの世代交代シミュレーションを終了させる方法は、上記以外の方法でも良い。
【０１０１】
三次元座標上の点の回転は３×３の回転行列を点の座標にかければ良い。この回転行列は上記α、β、γの値が決まれば決定する。よって、この回転行列を求めることにより、各画像での基準となる位置からのゴルフボールの回転量を決定することができる。最適化計算では、α、β、γを変量させることによりゴルフボールを回転させ、撮影映像とのずれが最小になるようにα、β、γを決定している。
また、上記手法で各々の映像における回転行列を求めることができるが、２つの映像間でどれだけ回転しているかと求めるためには、１つの映像の回転行列に、もう１つの映像の回転行列の逆行列をかければ良い。ここで求められた回転行列をもとに、回転軸ベクトルと回転角度を算出する。
【０１０２】
具体的には、適合度の最大値をもたらす個体が有する三次元座標空間のＸ軸、Ｙ軸、Ｚ軸回りの回転角度（α、β、γ）の数値より、三行三列の回転行列を求め、この回転行列が、ゴルフボール２０の姿勢を定める行列となっている。即ち、Ｘ軸、Ｙ軸、Ｚ軸の各軸回りに回転させる回転行列Ｒｘ、Ｒｙ、Ｒｚは、下記のように表記される。
【０１０３】
【数１１】

【０１０４】
一方、上記仮想ゴルフボールの基準姿勢及び基準位置から、姿勢認識プログラムで特定されるゴルフボール２０の姿勢及び位置を特定する行列である回転行列Ｒは、各軸回りの回転行列Ｒｘ、Ｒｙ、Ｒｚより下記のように表記される。
【０１０５】
【数１２】

【０１０６】
また、他のゴルフボール２０の映像に対しても上記と同様の演算および適合度の評価を行い、ゴルフボール２０の姿勢及び位置を特定する回転角度（α’、β’、γ’）より回転行列を求めている。
【０１０７】
以上より、上記数式１２で特定されるある時間のゴルフボール２０の姿勢を特定する回転行列をＲ１、他の時間のゴルフボール２０の姿勢を特定する回転行列をＲ２として、ある時間のゴルフボール２０の姿勢を他の時間のゴルフボール２０の姿勢へと一致させる回転行列をＲ１２とすると、回転行列Ｒ１２は下式で求められる。
【０１０８】
【数１３】

【０１０９】
上式で求められる回転行列Ｒ１２に伴う、図１１に示す座標変換される際の回転軸（座標軸）の方向（ｕ１、ｕ２、ｕ３）と回転角度Ψは、下記の数式１４、１５で求められる。
【０１１０】
【数１４】

【０１１１】
【数１５】

【０１１２】
上記求めた回転角度Ψと二つの映像間の撮影された時間間隔から、ある時間から他の時間にわたるゴルフボール２０の回転量と回転軸方法を求めている。
なお、二つ以上の映像間の回転量等を求めるには、各連続する前後の映像から上記同様にして各映像間の回転量等を求めて球体の回転量等を連続的に求めるようにしている。
【０１１３】
また、上記手法で各々の映像におけるゴルフボール重心位置の３成分の移動速度を求めることによって、打出角や振れ角等の飛行経路を含む飛行特性を算出できる。具体的には、ある時刻でのボールの重心位置Ｇｘ（ｔ）、Ｇｙ（ｔ）、Ｇｚ（ｔ）とすると、３成分の移動速度Ｖｘ、Ｖｙ、Ｖｚとボール速度Ｖは、下記の数式１６により求められる。
【０１１４】
【数１６】

【０１１５】
また、打出角及び振れ角は、下記の数式１７により求められる。ただし、ボールの進行方向をＸ軸、左右方向をＹ軸、高さ方向をＺ軸とする。
【０１１６】
【数１７】

【０１１７】
（実験１）
球体の回転特性の測定を行った。
【０１１８】
以下、本発明の実施例、比較例について詳述する。
予め既知の位置にマークを付けたゴルフボールを図９に示すような回転ステージに取り付け、ストロボ撮影によってボールの二次元映像を撮影した。ボールを既知の角度で回転ステージを回転させ、２枚写真（あるいは４枚写真）を撮影し、その２枚（４枚）の写真から求めたボールの回転量と回転ステージで読み取った回転量を比較し、精度の検証を行った。回転ステージは位置の固定をしておき、カメラの高さや位置を調整している。
【０１１９】
（ボールの撮影条件）
カメラ１の映像はカメラフィルム面座標でＸ＝０、Ｙ＝０、奥方向へＺ＝０で撮影し、カメラ２の映像を以下の４条件で撮影した。なお、ボール撮影位置の移動距離単位はｍｍとした。
（Ａ）ボールが大きく写っており、かつボールのエッジが良く見えるような好条件で撮影した。（カメラフィルム面座標でＸ＝０、Ｙ＝０、奥方向へＺ＝０（これを基準とする。））
（Ｂ）ボールが大きく写っており、ボールのエッジが片一方のみ見えにくい条件で撮影した。（カメラフィルム面座標でＸ＝１００、Ｙ＝１００、奥方向へＺ＝０）
（Ｃ）ボールが小さく写っており、かつボールのエッジが良く見えるような好条件で撮影した。（カメラフィルム面座標でＸ＝０、Ｙ＝０、奥方向へＺ＝−１００）
（Ｄ）ボールが小さく写っており、ボールのエッジが片一方のみ見えにくい条件で撮影した。（カメラフィルム面座標でＸ＝１００、Ｙ＝１００、奥方向へＺ＝−１００）
【０１２０】
（実施例１）
上記実施形態と同様の方法により、本発明の測定方法でスピン量を算出した。
【０１２１】
（比較例１）
撮影した二次元映像を基にボール映像の半径を単位半径としてマーク等のゴルフボール上に固定された点の三次元ベクトルを求め、複数の映像間での該ベクトルの回転角度からスピン量を求めた。ＧＡによる姿勢変位操作を行わなかった。
【０１２２】
（比較例２）
少なくとも１台の撮影手段を用いて、三次元座標と二次元座標の間の関係を予め導出することを行わなかった。即ち、仮想球体の三次元空間座標上のマーク位置を、二次元的な画像上に変換することを行わず、単に二次元方向へ投影することによって（視点を固定する）、変位姿勢操作を行った。その他は実施例１と同様とした。
【０１２３】
上記実施例及び比較例について、バックスピン成分、サイドスピン成分の計２通り測定した。測定された回転数と読み取り差を下記の表１、２に示す。
【０１２４】
【表１】

【０１２５】
【表２】

【０１２６】
表１、２に示すように、（Ａ）〜（Ｄ）のように撮影条件を変更しても、全般的に実施例１は読み取り値との差がかなり小さくなり、精度良くスピン量が求められていることが確認できた。実施例１では、ボールのエッジがはっきり見えていない場合や、ボール画像が画面の端部に撮影されたり、小さく撮影されても、測定精度が悪くならない事が確認できた。
即ち、ボールが大きく見える場合全ての例で精度は良く、ボールが小さく見える場合比較例１は精度が悪く、ボールが中心に見える場合全ての例で精度が良く、ボールが端部に見える場合比較例１と比較例２は精度が悪かった。
【０１２７】
（実験２）
飛行特性の測定を行った。
【０１２８】
（実施例２）
上記実施形態と同様の方法により、本発明の測定方法で振れ角、打出角、ボール速度を算出した。
【０１２９】
実打テストにおける振れ角、打出角、ボール速度をセンサを用いて正確に測定した。実施例２及びセンサによる測定結果を表３に示す。
【０１３０】
【表３】

【０１３１】
実施例２と、センサによる測定の結果がほぼ一致していた。よって、本発明は、センサによる測定とほぼ同程度で、振れ角、打出角、ボール速度を測定できることが確認できた。特に、従来の測定方法に比べ、振れ角の測定精度が極めて向上していることが確認できた。
【０１３２】
【発明の効果】
以上の説明より明らかなように、本発明の測定方法及び測定装置によれば、上記のような方法で球体の三次元姿勢及び位置を特定しているため、ある時間の三次元姿勢及び位置と、他の時間の三次元姿勢及び位置との関係から球体の回転数と回転軸方向等の回転特性、飛行経路や飛行速度等の飛行特性を容易かつ高精度で測定することができる。よって、外力を受けず回転している球体であるゴルフボールのスピン量や、打球の打出角や振れ角の測定に最適であり、特に振れ角の測定の精度を向上することができる。
【０１３３】
また、対象となる球体を適宜撮影すれば、球体の映像を基にコンピュータ内のプログラムに従って演算され、自動的に球体の回転特性や飛行特性を測定できるので、測定にかかる手間を大幅に削減することができる。
【０１３４】
さらに、本発明の測定方法では球体の輪郭のデータを用いないため、従来の測定方法と比べると大幅に測定精度を向上させることができる。よって、測定結果が撮影状況により左右されることを防止でき、フラッシュの光の当て方を詳細に調整するなどの撮影に要する手間にそれほど労力をかける必要もなく、高精度の測定が可能となり、測定に要するコストも低減できる。
【０１３５】
さらには、一枚の二次元映像から球体の姿勢及び位置を特定できるので姿勢及び位置の特定にかかる手間を削減できる。高速シャッター付きのカメラやマイクロフラッシュ等を利用することで、高速で移動する球体の回転特性や飛行特性をも求めることができる。
【図面の簡単な説明】
【図１】（Ａ）は、本発明の球体の回転特性と飛行特性の測定装置の構成を示す図、（Ｂ）は概略斜視図である。
【図２】本発明の球体の回転特性と飛行特性の測定方法のフローチャートである。
【図３】球体にマークを付した位置を示しており（Ａ）は斜視図、（Ｂ）平面図、（Ｃ）は底面図である。
【図４】球体の極座標を説明する図である。
【図５】マークを付した仮想球体の概略図である。
【図６】三次元空間上の座標と二次元的な画像上の座標の位置関係を示す図である。
【図７】仮想球体の三次元空間座標上のマーク位置を、二次元的な画像上に変換する方法を示す図である。
【図８】レンズ中心と球体との関係を示す図である。
【図９】ゴルフボール回転変位装置であり、（Ａ）は斜視図、（Ｂ）は側面図である。
【図１０】球体の二次元映像を示す図である。
【図１１】回転軸方向および回転角度を説明する図である。
【図１２】（Ａ）（Ｂ）は従来の測定方法による二次元ボール画像である。
【図１３】（Ａ）は別の従来の測定装置の概略図、（Ｂ）は（Ａ）の装置によるボール画像である。
【図１４】（Ａ）は他の別の従来の測定装置の斜視図、（Ｂ）は測定領域を通過したゴルフボールを示す三次元領域の斜視図である。
【符号の説明】
１０測定装置
１１ＣＣＤカメラ
１５コンピュータ
２０ゴルフボール
２１マーク
３０仮想ゴルフボール

Claims

表面に複数のマークを付した球体の回転飛行時に所定の時間間隔をあけて複数回撮影して、該球体の複数の二次元映像を得る一方、該球体と同様に表面にマークを付した仮想球体をコンピュータの三次元空間座標上に作成し、該仮想球体の任意姿勢及び任意位置を仮想球体の基準姿勢及び基準位置として設定すると共に、
少なくとも１台の撮影手段を用いて、三次元座標と二次元座標の間の関係を予め導出しておき、
上記関係を用いて、上記仮想球体の三次元空間座標上のマーク位置を、二次元的な画像上に変換することにより二次元の仮想マーク座標値を求めると共に、上記映像上に存在する球体の二次元の映像マーク座標値を求め、
上記仮想マーク座標値と、上記映像マーク座標値とが一致するように、上記仮想球体を、上記基準姿勢及び基準位置に対して姿勢変位操作し、該姿勢変位操作の量により、球体の各二次元映像毎に、球体の三次元姿勢及び位置を特定し、
上記二次元映像毎に特定された、ある時間の三次元姿勢及び位置と、他の時間の三次元姿勢及び位置とにより、球体の回転特性と飛行特性を演算していることを特徴とする球体の回転特性と飛行特性の測定方法。
上記三次元座標と二次元座標の間の関係を予め導出する際に、少なくとも６つ以上の三次元座標を用いている請求項１に記載の球体の回転特性と飛行特性の測定方法。
上記姿勢変位操作における設計変数として、上記球体の三次元重心位置及び仮想三次元空間の直交座標軸である横軸、縦軸、垂直軸の各軸回りの球体の回転角度の計６個を用いている請求項１または請求項２に記載の球体の回転特性と飛行特性の測定方法。
上記姿勢変位操作は、上記仮想球体の移動操作、回転操作であり、
上記基準姿勢及び基準位置に対する姿勢変位操作の量を遺伝的アルゴリズムという最適化手法に応じた演算で求めている請求項１乃至請求項３のいずれか１項に記載の球体の回転特性と飛行特性の測定方法。
上記球体表面の複数のマークの位置関係は、回転操作に対して４回以下の回転対称性を有している請求項１乃至請求項４のいずれか１項に記載の球体の回転特性と飛行特性の測定方法。
表面に複数のマークを付した球体が回転しながら飛行する時に、該球体を撮影できる撮影手段と、
上記撮影により得られた球体の二次元映像を記録する記録手段と、
上記球体と同様のマークを付した仮想球体を三次元座標空間上に作成し、該仮想球体と上記球体の二次元映像に基づいて、球体の三次元姿勢及び位置を特定し、球体の回転特性と飛行特性を求める演算手段とを備え、
上記演算手段は、少なくとも１台の撮影手段を用いて、三次元座標と二次元座標の間の関係を予め導出可能な座標変換プログラムを有し、
かつ、上記二次元映像上の球体の二次元の映像マーク座標値が、上記仮想球体の三次元座標上のマーク位置を変換することにより求められる二次元の仮想マーク座標値と一致するように、仮想球体を姿勢変位操作し、上記仮想球体の基準姿勢及び基準位置に対する姿勢変位操作の量により、上記球体の三次元姿勢及び位置を特定する姿勢認識プログラムを有してることを特徴とする球体の回転特性と飛行特性の測定装置。
上記演算手段は、上記基準姿勢及び基準位置に対する姿勢変位操作の量を遺伝的アルゴリズムに応じて演算する最適化プログラムを有している請求項６に記載の球体の回転特性と飛行特性の測定装置。
上記撮影手段は、回転や移動している球体を所定の時間間隔をあけて複数回撮影可能な構成としている請求項６または請求項７に記載の球体の回転特性と飛行特性の測定装置。