JP5435543B2 - 人体の部位別の質量および重心の測定方法 - Google Patents

Description

本発明は、人体の部位別の質量および重心の測定方法に関する。さらに詳しくは、統計データに頼ることなく個人ごとに測定できる人体の部位別の質量および重心の測定方法に関する。

医療分野やスポーツ運動学分野などにおいて、運動時に人体の関節に働くモーメント（関節モーメント）を推定することは非常に重要となっている。例えば、推定された関節モーメントから、筋骨格モデルを用いて筋張力を推定することなどが行われている。
関節モーメントを推定するためには、人体の部位別の質量、重心および慣性モーメント
を測定することが必要である。測定した質量などと、運動時に人体の各部位に働く加速度および関節まわりの角加速度とから関節モーメントを推定することができるからである。
従来、人体の部位別の質量、重心および慣性モーメントを測定するには、個人ごとに測定する方法がないため、解剖などから得られた統計的なデータを用いて推定していた。
しかるに、人体の部位別の質量などは個体差があるため、本来は個人ごとに測定すべきであり、統計的なデータを用いると、個人ごとの正確な関節モーメントを推定することができないという問題がある。

これに対して、人体の部位別の質量を個人ごとに測定する従来技術として特許文献１がある。特許文献１の測定方法は、頭部、体幹部、左右の各腕部、左右の各脚部の６つの部位の各々の質量を測定する部位別質量測定方法であって、各部位の位置に対応するように板状の質量測定装置を床に配置し、その上に被測定者が仰臥位または伏臥位で横たわり、各部位の質量を測定するものである。

しかるに、特許文献１の測定方法では、各部位が関節で繋がっている他の部位の影響を受け、関節に力が入っている時などに、各部位の質量を正確に測定することができないという問題がある。
また、部位ごとに測定装置を用意する必要があるため、小さい部位の質量を測定することが困難であるという問題がある。
さらに、各部位の重心や慣性モーメントを測定できず、これのみでは関節モーメントを推定することができないという問題がある。

特開２００７−２１２４１１号公報

本発明は上記事情に鑑み、個人ごとの部位別の質量および重心を測定する方法を提供することを目的とする。

第１発明の人体の部位別の質量および重心の測定方法は、被測定者の部位に、該部位の質量が集中していると仮定する測定点を設定する第１ステップと、前記測定点の基準姿勢からの変位x_i（iは測定点の番号を示す。）と、任意の点周りの被測定者の体全体により加えられるモーメントMとを、前記測定点を設定した部位のみ姿勢を変化させながら断続的に測定する第２ステップと、
M/g=Σm_ix_i＋m’x’
（gは重力加速度、m_iは測定点の質量、m’は測定点以外の被測定者の質量、x’はm’の重心の前記任意の点からの距離、Σは測定点の総和を示す。）
の式を用いて、前記第２ステップにより得られたデータから、前記部位の測定点における質量m_iを、M/gを従属変数、x_iを独立変数、m_iを偏回帰変数とした重回帰分析により推定する第３ステップと、前記測定点の質量m_iの総和から前記部位の質量を求める第４ステップと、前記測定点の質量m_iと位置とから前記部位の重心を求める第５ステップとを順に行うことを特徴とする。
第２発明の人体の部位別の質量および重心の測定方法は、第１発明において、前記第１ステップにおいて、１つの部位に対して２つの測定点を設定することを特徴とする。
第３発明の人体の部位別の質量および重心の測定方法は、第１発明において、前記第１ステップにおいて、１つの部位に対して３つの測定点を設定することを特徴とする。
第４発明の人体の部位別の質量および重心の測定方法は、第１発明において、前記第１ステップにおいて、１つの部位に対して４つ以上の測定点を同一平面上にならないように設定することを特徴とする。
第５発明の人体の部位別の質量および重心の測定方法は、第１発明において、前記第１ステップにおいて、１つの部位にのみ測定点を設定することを特徴とする。
第６発明の人体の部位別の質量および重心の測定方法は、第１発明において、前記第１ステップにおいて、複数の部位に測定点を設定することを特徴とする。
第７発明の人体の部位別の質量および重心の測定方法は、第１発明において、前記第２ステップにおいて、鉛直床反力を測定し、前記第３ステップにおいて、前記鉛直床反力が基準値よりも大きいもしくは小さい場合は、重回帰分析に用いるデータから取り除くことを特徴とする。
第８発明の人体の部位別の質量および重心の測定方法は、被測定者の部位に、該部位の質量が集中していると仮定する測定点を設定する第１ステップと、前記測定点の基準姿勢からの変位x _i （iは測定点の番号を示す。）と、前記基準姿勢の被測定者の床反力の重心位置を基準点とし、該基準点周りの被測定者の体全体により加えられるモーメントMとを、前記測定点を設定した部位のみ姿勢を変化させながら断続的に測定する第２ステップと、
M/g=Σm_ix_i
（gは重力加速度、m _i は測定点の質量、Σは測定点の総和を示す。）
の式を用いて、前記第２ステップにより得られたデータから、前記部位の測定点における質量m_iを、M/gを従属変数、x _i を独立変数、m _i を偏回帰変数とした重回帰分析により推定する第３ステップと、前記測定点の質量m _i の総和から前記部位の質量を求める第４ステップと、前記測定点の質量m _i と位置とから前記部位の重心を求める第５ステップとを順に行うことを特徴とする。
第９発明の人体の部位別の質量および重心の測定方法は、被測定者の部位に、該部位の質量が集中していると仮定する測定点を設定する第１ステップと、前記測定点の基準姿勢からの変位x _i （iは測定点の番号を示す。）と、任意の点周りの被測定者の体全体により加えられるモーメントMとを、前記測定点を設定した部位のみ姿勢を変化させながら断続的に測定する第２ステップと、
(dM/dt)/g=Σm_i(dx_i/dt)
（dM/dtはモーメントMの時間微分、dx_i/dtは変位x_iの時間微分、gは重力加速度、m _i は測定点の質量、Σは測定点の総和を示す。）
の式を用いて、前記第２ステップにより得られたデータから、前記部位の測定点における質量m_iを、(dM/dt)/gを従属変数、(dx _i /dt)を独立変数、m _i を偏回帰変数とした重回帰分析により推定する第３ステップと、前記測定点の質量m _i の総和から前記部位の質量を求める第４ステップと、前記測定点の質量m _i と位置とから前記部位の重心を求める第５ステップとを順に行うことを特徴とする。
第１０発明の人体の部位別の質量および重心の測定方法は、被測定者の全ての部位に、該部位の質量が集中していると仮定する測定点を設定する第１ステップと、前記測定点の基準姿勢からの変位x _i （iは測定点の番号を示す。）と、任意の点周りの被測定者の体全体により加えられるモーメントMとを、前記測定点を設定した部位の姿勢を変化させながら断続的に測定する第２ステップと、
M/g=Σm_ix_i+(F_z-Σm_i)x_N
（gは重力加速度、m _i は測定点の質量、F_zは鉛直床反力、Nは測定点の数、ΣはN番目の測定点以外の測定点の総和を示す。）
の式を用いて、前記第２ステップにより得られたデータから、前記部位の測定点における質量m_iを、M/gを従属変数、x _i を独立変数、m _i を偏回帰変数とした重回帰分析により推定する第３ステップと、前記測定点の質量m _i の総和から前記部位の質量を求める第４ステップと、前記測定点の質量m _i と位置とから前記部位の重心を求める第５ステップとを順に行うことを特徴とする。

第１発明によれば、被測定者の部位の質量および重心を直接測定することができるので、統計的データに頼ることなく正確に、個人ごとの部位別の質量および重心を測定することができる。
第２発明によれば、１つの部位に対して２つの測定点が設定されているので、測定点の質量の測定結果から、２つの測定点を結ぶ線上のどの位置に重心があるかを求めることができる。
第３発明によれば、１つの部位に対して３つの測定点が設定されているので、測定点の質量の測定結果から、３つの測定点からなる平面上のどの位置に重心があるかを求めることができる。
第４発明によれば、１つの部位に対して４つ以上の測定点が設定されているので、測定点の質量の測定結果から、３次元空間上のどの位置に重心があるかを求めることができる。
第５発明によれば、１つの部位にのみ測定点が設定されているので、測定装置が単純となり、測定データから質量および重心を求めることも容易となる。
第６発明によれば、複数の部位に測定点が設定されているので、短時間で多数の部位の質量および重心を求めることができる。
第７発明によれば、被測定者が十分静的でない場合の測定データを重回帰分析から取り除くことができるので、より正確に測定点の質量を推定することができる。
第８発明によれば、Mを基準点周りの被測定者の体全体により加えられるモーメントとすることで、偏回帰変数を１つ消去することができるので、重回帰分析の精度を向上させることができる。
第９発明によれば、微分することにより測定点以外の被測定者の質量m’の項を消去することができ、偏回帰変数を１つ消去することができるので、重回帰分析の精度を向上させることができる。
第１０発明によれば、m_Nの項を消去することにより、偏回帰変数を１つ消去することができるので、重回帰分析の精度を向上させることができる。

本発明に係る測定方法に用いる代表的な測定装置の説明図である。 前腕にマーカーを取り付ける際の説明図である。 手掌にマーカーを取り付ける際の説明図である。 頭部にマーカーを取り付ける際の説明図である。 本発明に係る測定方法の説明図である。 同測定方法の説明図である。測定部位以外の人体を省略している。

（測定装置）
つぎに、本発明に係る測定方法に用いる代表的な測定装置について図面に基づき説明する。
図１に示すように、本発明に係る測定方法に用いる代表的な測定装置は、主に床反力計と３次元動作解析装置とで構成される。

床反力計は、ひずみゲージや圧電素子などで構成されており、床反力計の上に被測定者が乗ると、被測定者に加わる床反力の大きさや向き、あるいはその床反力の重心（加重中心）を測定することができる荷重センサである。床反力計には、鉛直、前後、左右方向の床反力を測定することができるものもあるが、本発明に用いる床反力計は、少なくとも鉛直方向の床反力と加重中心を測定することができるものであればよい。
ここで、床反力とは人体が床から受ける力を意味し、静止時の力の大きさはその人体が体重により床を押す力の大きさと等しくなる。また、床反力は床と人体との接触面に分布を有するが、その床反力の分布の重心が加重中心である。

３次元動作解析装置は、主に３次元動作解析カメラと画像解析を行う電子計算機とからなり、被測定者の体に取り付けたマーカーの３次元位置を測定する３次元位置測定装置である。一般的には、複数台の３次元動作解析カメラでマーカーを取り付けた被測定者の動作を撮影し、その動画を電子計算機に組み込まれたソフトウェアで解析することで各マーカーの位置、速度、加速度などを求めるものである。

なお、３次元位置の座標は任意に設定すればよいが、説明を簡単にするため、以下では床反力計の表面を基準として、その表面に水平にx軸、y軸を、垂直にz軸を設定する。原点は３次元位置測定装置にあわせて任意に設定すればよく、例えば床反力計の角などに設定すればよい。

本発明に係る測定方法では、床反力計と３次元動作解析装置とから同じタイミングで断続的にデータを取得する必要がある。これは後述するように、被測定者の１つの姿勢に対して床反力と３次元位置とを測定する必要があり、また、複数の姿勢で測定する必要があるからである。
そのため、床反力計および３次元動作解析装置を制御し、データを収集する制御装置（図示せず）と、収集したデータを後述する方法で解析する電子計算機（図示せず）などが、必要に応じて用いられる。

（１つの部位の測定方法）
つぎに、本発明の実施形態を図面に基づき説明する。
１つの部位の質量および重心を測定する方法が基本となるので、はじめにその実施形態について説明する。

[第１ステップ マーカーの設定]
まず、被測定者の測定する１つの部位に３次元動作解析装置用のマーカーを取り付ける。後述するように、このマーカーの位置にその部位の質量が集中していると仮定して、各マーカー位置の質量を測定し、各マーカー位置の質量の総和からその部位の質量を、各マーカー位置の質量の分布からその部位の重心を求めるのである。
本実施形態において、マーカー位置は特許請求の範囲に記載の「測定点」に相当する。

部位の重心を求めるためには、１つの部位に対して少なくとも２つのマーカーを取り付けることが必要である。
図２に示すように、例えば前腕を測定する場合、前腕は梁状の部材であり、その重心は梁の軸上にあると予想されるので、肘付近と手首付近の２か所にマーカーを取り付ければよい。マーカー位置の質量の測定結果から、２つのマーカーを結ぶ線上のどの位置に重心があるかを求めることができる。

また、図３に示すように、手掌のような平面状の部位を測定する場合、その重心は平面上にあると予想されるので、マーカーを３か所取り付ければよい。マーカー位置の質量の測定結果から、３つのマーカーからなる平面上のどの位置に重心があるかを求めることができる。

あるいは、図４に示すように、頭部のような立体的な形状の部位を測定する場合、マーカーを少なくとも４か所、同一平面上にならないように取り付ければよい。マーカーを４か所取り付けることにより、マーカー位置の質量の測定結果から、３次元空間上のどの位置に重心があるかを求めることができる。

もちろん、梁状あるいは平面状の部位に対してマーカーを４か所取り付けてもよい。また、１つの部位に対して４か所以上マーカーを取り付けてもよい。

なお、マーカーを３か所取り付け、測定回数を２回にすることでも３次元空間上の重心位置を求めることができる。この場合、まず測定する部位にマーカーを３か所取り付け、その３つのマーカーからなる平面上（例えば鉛直面）の重心位置を求める。つぎに、マーカーの位置を変更して、再度、３つのマーカーからなる新たな平面上（例えば水平面）の重心位置を求める。１回目に求めた重心位置と２回目に求めた重心位置から３次元空間上の重心を求めることが可能である。
つまり、平面状の部位に限らずとも、立体的な形状の部位に対して３つのマーカーを取り付けて、３次元空間上の重心位置を求めることも可能である。

測定する部位によって、適したマーカーの数や位置は異なるので、質量および重心の測定精度を考慮し、各部位に対して最適なマーカーの取り付け方を選択すればよい。

[第２ステップ マーカー位置およびモーメントの測定]
つぎに、被測定者が床反力計の上に乗り、加速度が重力加速度に比べて十分小さくなるようにマーカーを取り付けた部位のみをゆっくり動かし、同時にマーカーの位置を３次元動作解析装置で、床反力を床反力計により測定する。マーカーの位置と床反力の測定データは、部位の姿勢を変えながら断続的に収集する。すなわち、測定したい部位の姿勢を様々に変化させ、各姿勢に対する床反力とマーカーの位置を測定するのである。

なお、図１に示すように被測定者は測定中、床反力計の上で椅子に座った状態でもよいし、立った状態でもよい。

[第３ステップ マーカー位置の質量の推定]
つぎに、収集した測定データから各マーカー位置の質量を推定する。
図５および図６に示すように、測定する部位の測定開始時の姿勢を基準姿勢とし、その基準姿勢での被測定者の床反力の重心位置（加重中心）を基準点とする。なお、基準姿勢は測定開始時の姿勢とする必要はなく、測定の都合に合わせて任意のタイミングに測定した任意の姿勢に設定してもよい。後述するとおり、基準姿勢からのマーカー位置の変位、基準点からの加重中心の変位が測定できればよいからである。特許請求の範囲に記載の「基準姿勢」は、任意の基準となる姿勢を意味する。

姿勢が変化したときの加重中心の、基準点からのx軸方向の変位をx_G、y軸方向の変位をy_Gとし、そのときの鉛直床反力をF_zとすると、基準点まわりの被測定者の体全体により加えられるモーメントのx成分M_x、y成分M_yは次式で与えられる。

また、各マーカー位置の質量をm_i、各マーカーの基準姿勢からのx軸方向の変位をx_i、y軸方向の変位をy_iとすると次式も成り立つ。

ここでiは部位におけるマーカーの番号であり、その部位にNか所のマーカーを取り付けた場合、i=1,2,・・・,Nである。また、gは重力加速度である。

式３、式４においてM_x、M_yは床反力計の測定データ、すなわちx_G、y_G、F_zより得られ、x_i、y_iは３次元動作解析装置の測定データより得られる。したがって、M_x/g、M_y/gを従属変数、x_i、y_iを独立変数、m_iを偏回帰係数とすれば、重回帰分析により各マーカー位置の質量m_iを推定することができる。

ここで重回帰分析とは、独立変数X_i（i=1,2,・・・,N N:独立変数の数）と従属変数Yとの間に次式の関係が成り立つと予想される場合に、

Yの実測値y_k（k=1,2,・・・,n n:測定データ数）と、X_iの実測値から式５で与えられる予測値Y_k（k=1,2,・・・,n n:測定データ数）との差（残差）e_kの二乗和Qが最小になるような偏回帰係数b_iを求める方法である。

また、重回帰分析で求められた偏回帰係数の確からしさは、重相関係数などの値で表わされる。

式３、式４からは重回帰分析によりそれぞれ独立に各マーカー位置の質量m_iを推定するため、１つのマーカー位置の質量に対して、式３から推定された値と、式４から推定された値の２つが存在することになる。理想的には、式３の推定値と式４の推定値は一致するはずであるが、実際には、一致することはあまりない。そのため、式３の重回帰分析と式４の重回帰分析で得られる残差の二乗和Qを比較し、Qが小さい方の式で得られた推定値m_iをマーカー位置の質量として採用すればよい。また、重相関係数などの指標から、どちらの式で得られた推定値m_iを採用するかを判断してもよい。

あるいは、床反力計の加重位置の測定がx軸方向のみ可能であり、３次元動作解析装置がx軸方向のみ測定可能な場合であれば式４を用い、逆に、床反力計の加重位置の測定がy軸方向のみ可能であり、３次元動作解析装置がy軸方向のみ測定可能な場合であれば式３を用いるというように、測定装置の制限により用いる式を選択してもよい。

また、式３、式４を独立に重回帰分析を行うのではなく、式３、式４合わせて重回帰分析を行ってもよい。
すなわち式３の残差を

とし、式４の残差を

としたときに、e_xkの二乗和とe_ykの二乗和Qが最小になるような偏回帰係数m_iを求めてもよい。

この場合は、一意にマーカー位置の質量m_iを推定することができる。

ところで、鉛直床反力F_zは被測定者が加速度運動をしない限り以下の式を満たす。

ここでm’は測定する部位以外の被測定者の質量である。被測定者が十分静的であればF_zは、被測定者の体重により床を押す力の大きさと等しくなり一定のはずである。F_zは被測定者の体重により得られるので、例えば事前に測定した被測定者の体重の±1%を基準値として、F_zの値がその基準値よりも大きいもしくは小さい測定データに関しては、被測定者が十分静的でないとして、前記の重回帰分析に使用するデータから取り除くという操作をしてもよい。この操作により、より正確にマーカー位置の質量m_iを推定することができる。
特許請求の範囲に記載の「基準値」は、被測定者が十分静的でない場合の測定データを取り除くための、鉛直床反力F_zに対する任意の値を意味する。

なお、本実施形態では、M_x、M_yを基準点まわりの被測定者の体全体により加えられるモーメントとしたが、この代わりに、M_x、M_yを任意の点周りの被測定者の体全体により加えられるモーメントとしてもよい。この任意の点には、例えば床反力計の角に設定した３次元座標の原点などを用いることができる。
この場合、式３、式４には、測定する部位以外の被測定者の質量m’の項が現れ、

となる。ここで、x’、y’は測定する部位以外の被測定者の質量m’の前記任意の点からの距離である。被測定者は測定する部位のみ動かし、測定する部位以外の部分は動かさない。つまり、x’、y’は一定のはずである。したがってm’x’、m’y’を式５におけるb₀として式１２および式１３を用いて重回帰分析をすれば、マーカー位置の質量m_iを推定することができる。

しかし、この場合偏回帰係数の数が式３、式４に比べて１つずつ増えるので、重回帰分析の精度が悪くなる可能性がある。この点、M_x、M_yを基準点まわりのモーメントとすることで、式５におけるb₀の項を消去した式３、式４を用いるほうが好ましい。

あるいは、m’x’、m’y’が一定であることから、式１２および式１３を時間で微分し、

として、m’の影響を排除した式を用いて重回帰分析を行ってもよい。式１４、式１５では式３、式４と同様にb₀に相当する項が消去されるので偏回帰係数の数を少なくすることができ、重回帰分析に好適である。

なお、式１４、式１５は、荷重センサや３次元位置測定装置で直接測定される値が、時間で微分したものである場合にも好適である。それらの測定値を積分する必要なく重回帰分析を行うことができるためである。

以上の式３、式４、式１２、式１３、式１４、式１５は測定に用いる測定装置などの種々の測定条件によって適した式を選択すればよい。

[第４ステップ 部位の質量の推定]
つぎに、推定されたマーカー位置の質量m_iから、マーカーを取り付けた部位の質量mを求める。部位の質量mはマーカー位置の質量m_iの総和で求められる。

[第５ステップ 部位の重心の推定]
つぎに、マーカー位置の質量m_iとマーカー位置とから、その部位の重心x_g、y_g,、z_gを求める。重心は以下の式で与えられる。

x_i、y_i、z_iは３次元動作解析装置で測定した値を用いてもよいし、被測定者の部位を基準としてマーカーを取り付けた位置を別途測定した値を用いてもよい。

以上の手順で、１つの部位の質量と重心を測定することができる。 被測定者の部位の質量および重心を直接測定することができるので、統計的データに頼ることなく正確に、個人ごとの部位別の質量および重心を測定することができる。
また、１つの部位のみの測定では、測定装置が単純となり、測定データから質量および重心を求めることも容易となる。
他の部位の質量を測定するためには、マーカーを付け替えて以上の方法を繰り返すことで測定することもできる。

（複数の部位の測定方法）
つぎに、複数の部位の質量および重心を同時に測定する方法について説明する。

[第１ステップ マーカーの設定]
まず、被測定者の複数の部位に３次元動作解析装置用のマーカーを取り付ける。被測定者の体のすべての部位にマーカーを取り付け、すべての部位を同時に測定することも可能であるが、マーカーの数が多いと、マーカー位置の質量の測定精度が落ちるため、上腕と前腕、大腿と下腿というように、部位を限定する方が好ましい。各部位のマーカーの取り付け方は前述の方法と同様である。

[第２ステップ マーカー位置およびモーメントの測定]
前述の方法と同様に、被測定者が床反力計の上に乗り、加速度が重力加速度に比べて十分小さくなるようにマーカーを取り付けた部位のみをゆっくり動かし、同時にマーカーの位置を３次元動作解析装置で、床反力を床反力計により測定する。マーカーの位置と床反力の測定データは、部位の姿勢を変えながら断続的に収集する。

[第３ステップ マーカー位置の質量の推定]
つぎに、収集した測定データから各マーカー位置の質量を推定する。
複数の部位にマーカーが取り付けられているので、前記の式３、式４の代わりに次式を用いる。

ここでjはマーカーを取り付けた部位の番号であり、Mか所の部位にマーカーを取り付けた場合、j=1,2,・・・,Mである。すなわち、m_jiはj番目部位のi番目のマーカー位置の質量を意味し、x_ji、y_jiはj番目部位のi番目マーカーの基準姿勢からのx方向の変位およびy方向の変位を意味する。式２０、式２１を用いて前述と同様の方法で重回帰分析を行うことにより、各マーカー位置の質量m_jiを推定することができる。

もちろん、式１２、式１３、式１４、式１５についても同様に変換し、その式を用いて重回帰分析を行うことにより、各マーカー位置の質量m_jiを推定することもできる。

なお、被測定者の全ての部位に３次元動作解析装置用のマーカーを取り付ける場合、式１１において、測定する部位以外の被測定者の質量m’がなくなり、

となる。式２２は以下のように変形できるので、

式２３を式２０、式２１に代入して

とすることができる。すなわち、m_MNを消去することで偏回帰係数m_jiを１つ少なくすることができ、重回帰分析の精度を向上させることができるのである。
式２４、式２５を用いて重回帰分析を行った後、m_MN以外の質量の推定値m_jiを用いて式２３からm_MNを求めることができる。

[第４、第５ステップ 部位の質量および重心の推定]
式２０、式２１もしくは式２４、式２５を用いて推定された各マーカー位置の質量m_jiから、各部位の質量および重心を求める。その方法は前述と同様であるので説明を省略する。

以上の手順で、複数の部位の質量および重心を同時に測定することができる。被測定者の部位の質量および重心を直接測定することができるので、統計的データに頼ることなく正確に、個人ごとの部位別の質量および重心を測定することができる。
また、複数の部位の質量および重心を同時に測定することができるので、短時間で多数の部位の質量および重心を求めることができる。

（他の測定装置）
以上の実施形態では、荷重センサとして床反力計を用い、３次元位置測定装置として３次元動作解析装置を用いたが、これに代わる装置を用いてもよい。
３次元位置測定装置としては、例えば、被測定者の各関節にゴニオメーターを取り付け関節の角度を測定し、その測定値と事前に測定した被測定者の各部位の長さとから各部位の３次元位置を計算して求めてもよい。この場合、床と接触する足から質量などを測定する部位までの各部位の長さと、各関節の角度を測定することが必要である。
あるいは、ジャイロ，加速度計，地磁気センサを組み合わせたセンサを被測定者が身に着け、その測定値から地球を基準とした３次元位置を求めてもよい。
いずれの３次元位置測定装置の場合においても、３次元位置の測定方法以外は、本発明に係る測定方法は同じである。

本発明に係る人体の部位別の質量および重心の測定方法は、医療分野やスポーツ運動学分野などにおいて利用される。例えば、本発明の測定方法により得られた部位別の質量および重心は、個人ごとに関節モーメントを求める際のパラメータとして利用され、ひいては筋骨格モデルを用いて個人ごとの筋張力を推定することなどに利用される。
また、人体をモデル化して運動をシミュレーションすることが様々な分野で行われているが、本発明の測定方法により得られた部位別の質量および重心は、個人ごとの運動をシミュレーションする際のパラメータとしても利用される。

Claims (10)

1. 被測定者の部位に、該部位の質量が集中していると仮定する測定点を設定する第１ステップと、
   前記測定点の基準姿勢からの変位x_i（iは測定点の番号を示す。）と、
   任意の点周りの被測定者の体全体により加えられるモーメントMとを、
   前記測定点を設定した部位のみ姿勢を変化させながら断続的に測定する第２ステップと、
   M/g=Σm_ix_i＋m’x’
   （gは重力加速度、m_iは測定点の質量、m’は測定点以外の被測定者の質量、x’はm’の重心の前記任意の点からの距離、Σは測定点の総和を示す。）
   の式を用いて、前記第２ステップにより得られたデータから、前記部位の測定点における質量m_iを、M/gを従属変数、x_iを独立変数、m_iを偏回帰変数とした重回帰分析により推定する第３ステップと、
   前記測定点の質量m_iの総和から前記部位の質量を求める第４ステップと、
   前記測定点の質量m_iと位置とから前記部位の重心を求める第５ステップとを順に行う
   ことを特徴とする人体の部位別の質量および重心の測定方法。
2. 前記第１ステップにおいて、
   １つの部位に対して２つの測定点を設定する
   ことを特徴とする請求項１記載の人体の部位別の質量および重心の測定方法。
3. 前記第１ステップにおいて、
   １つの部位に対して３つの測定点を設定する
   ことを特徴とする請求項１記載の人体の部位別の質量および重心の測定方法。
4. 前記第１ステップにおいて、
   １つの部位に対して４つ以上の測定点を同一平面上にならないように設定する
   ことを特徴とする請求項１記載の人体の部位別の質量および重心の測定方法。
5. 前記第１ステップにおいて、
   １つの部位にのみ測定点を設定する
   ことを特徴とする請求項１記載の人体の部位別の質量および重心の測定方法。
6. 前記第１ステップにおいて、
   複数の部位に測定点を設定する
   ことを特徴とする請求項１記載の人体の部位別の質量および重心の測定方法。
7. 前記第２ステップにおいて、
   鉛直床反力を測定し、
   前記第３ステップにおいて、
   前記鉛直床反力が基準値よりも大きいもしくは小さい場合は、重回帰分析に用いるデータから取り除く
   ことを特徴とする請求項１記載の人体の部位別の質量および重心の測定方法。
8. 被測定者の部位に、該部位の質量が集中していると仮定する測定点を設定する第１ステップと、
   前記測定点の基準姿勢からの変位x _i （iは測定点の番号を示す。）と、
   前記基準姿勢の被測定者の床反力の重心位置を基準点とし、該基準点周りの被測定者の体全体により加えられるモーメントMとを、
   前記測定点を設定した部位のみ姿勢を変化させながら断続的に測定する第２ステップと、
   M/g=Σm_ix_i
   （gは重力加速度、m _i は測定点の質量、Σは測定点の総和を示す。）
   の式を用いて、前記第２ステップにより得られたデータから、前記部位の測定点における質量m_iを、M/gを従属変数、x _i を独立変数、m _i を偏回帰変数とした重回帰分析により推定する第３ステップと、
   前記測定点の質量m _i の総和から前記部位の質量を求める第４ステップと、
   前記測定点の質量m _i と位置とから前記部位の重心を求める第５ステップとを順に行う
   ことを特徴とする人体の部位別の質量および重心の測定方法。
9. 被測定者の部位に、該部位の質量が集中していると仮定する測定点を設定する第１ステップと、
   前記測定点の基準姿勢からの変位x _i （iは測定点の番号を示す。）と、
   任意の点周りの被測定者の体全体により加えられるモーメントMとを、
   前記測定点を設定した部位のみ姿勢を変化させながら断続的に測定する第２ステップと、
   (dM/dt)/g=Σm_i(dx_i/dt)
   （dM/dtはモーメントMの時間微分、dx_i/dtは変位x_iの時間微分、gは重力加速度、m _i は測定点の質量、Σは測定点の総和を示す。）
   の式を用いて、前記第２ステップにより得られたデータから、前記部位の測定点における質量m_iを、(dM/dt)/gを従属変数、(dx _i /dt)を独立変数、m _i を偏回帰変数とした重回帰分析により推定する第３ステップと、
   前記測定点の質量m _i の総和から前記部位の質量を求める第４ステップと、
   前記測定点の質量m _i と位置とから前記部位の重心を求める第５ステップとを順に行う
   ことを特徴とする人体の部位別の質量および重心の測定方法。
10. 被測定者の全ての部位に、該部位の質量が集中していると仮定する測定点を設定する第１ステップと、
    前記測定点の基準姿勢からの変位x _i （iは測定点の番号を示す。）と、
    任意の点周りの被測定者の体全体により加えられるモーメントMとを、
    前記測定点を設定した部位の姿勢を変化させながら断続的に測定する第２ステップと、
    M/g=Σm_ix_i+(F_z-Σm_i)x_N
    （gは重力加速度、m _i は測定点の質量、F_zは鉛直床反力、Nは測定点の数、ΣはN番目の測定点以外の測定点の総和を示す。）
    の式を用いて、前記第２ステップにより得られたデータから、前記部位の測定点における質量m_iを、M/gを従属変数、x _i を独立変数、m _i を偏回帰変数とした重回帰分析により推定する第３ステップと、
    前記測定点の質量m _i の総和から前記部位の質量を求める第４ステップと、
    前記測定点の質量m _i と位置とから前記部位の重心を求める第５ステップとを順に行う
    ことを特徴とする人体の部位別の質量および重心の測定方法。

Images