目次第 1 章緒言研究背景研究目的脊柱の構造... 3 第 2 章腰痛診断の現状腰痛症の種類従来の姿勢変化の計測方法... 4 第 3 章提案する推定法の概要... 6 第 4 章精度検証... 8

Size: px

Start display at page:

Download "目次第 1 章緒言研究背景研究目的脊柱の構造... 3 第 2 章腰痛診断の現状腰痛症の種類従来の姿勢変化の計測方法... 4 第 3 章提案する推定法の概要... 6 第 4 章精度検証... 8"

ちとらいせき
5 years ago
Views:

1 平成 24 年度修士論文椎間板負荷の非侵襲的な推定指導教員井上喜雄教授芝田京子准教授高知工科大学大学院工学研究科基盤工学専攻知能機械システム工学コース学籍番号藤井涼

2 目次第 1 章緒言研究背景研究目的脊柱の構造... 3 第 2 章腰痛診断の現状腰痛症の種類従来の姿勢変化の計測方法... 4 第 3 章提案する推定法の概要... 6 第 4 章精度検証精度検証の実験内容精度検証の実験結果第 5 章椎間板負荷の推定腰椎部の曲率半径の推定椎骨間の隙間の推定椎間板負荷の推定第 6 章椎間板負荷の推定実験椎間板負荷の推定実験内容椎間板負荷の推定結果第 7 章結言謝辞参考文献

3 第 1 章緒言 1.1 研究背景日本では腰痛症に悩む患者が増加傾向にある. その背景として, 高齢者の総人口に占める割合が増加したこと. 生活環境や労働環境の変化によりデスクワークなど長時間同じ姿勢で仕事をする人が増えたことが挙げられ, 若年層でも腰痛症で悩まされる人が増えていることから年齢を問わない問題となっている. 腰痛症を引き起こす原因には姿勢の崩れからくるもの, 前屈みになるなど腰椎部にかかる負荷が大きい姿勢を連続して行うなど様々ある. 腰痛症の中でも代表的なものに腰椎椎間板ヘルニアがある. 発症する原因の一つとして, 脊柱を構成する椎骨間の緩衝材としての役割を果たす椎間板に過度な負荷が加わり, 椎間板が突出状態となる. この状態の椎間板が神経を圧迫することで痛みや痺れといった症状として現われる. 現在までに報告されている椎間板にかかる負荷を測定する有用な手法はスウェーデンの整形外科医である Nachemson (1) らが行った椎間板にセンサを直接挿入し様々な姿勢時での椎間板にかかる圧力を測定した方法がある. この方法ならば, 直接椎間板に加わる圧力を測定するため高い精度が得られ, 姿勢が変化していても測定可能である. しかし, 侵襲的な測定は外科手術を要するため限定された研究機関でしか行えない. 椎間板付近には人体に重要な脊髄神経が通っており, これを傷つける危険性があるため安全性に欠けた測定法であることから広く用いることは困難である. また, 非侵襲的な方法としてレントゲンや CT,MRI など様々あるが, 静止した姿勢での計測しか行うことができず, 定量的なデータの取得は難しいなどといった問題点が挙げられる. 1.2 研究目的様々な姿勢や動作中での腰椎部の椎間板に加わる負荷を知ることができれば腰痛の診断や予防法, 治療法の進歩に非常に有用であり, 整形外科分野の発展に大いに役立つものと考えられる. 本研究は,3 次元動作解析装置を用いて身体に負担をかけない非侵襲的な測定法を提案する. 腰椎部の椎骨端部が背中の体表面に近い部分に存在していることに着目し, 椎骨端部の動きを 3 次元動作解析装置を使用することで外部から測定, その動きから椎骨間の隙間を推定し, 椎間板負荷を評価する方法を考える. これにより安全かつ容易に動作中や運動中の椎間板負荷を測定することが可能になる. 本稿では提案する椎間板負荷の推定法を述べる. 次いで提案する手法を用いて推定した椎間板負荷と Nachemson の実験結果を比較検討することで, 本手法の有効性を検討する. 2

4 1.3 脊柱の構造図 1.1 に示すように人間の脊柱は椎骨と椎間板が連続的に積み重なる構造となっている.7 個の頸椎,12 個の胸椎,5 個の腰椎で構成されており,S 字カーブを描くように湾曲しているため, 頭部の重量を分散し支えることができている. それらの椎骨は硬く剛体と見なすことができ, 前方では緩衝剤の役割を担う椎間板, 後方では左右の椎間関節により連結されている. また, 屈曲や伸展により, 椎骨の位置が変化することで椎骨間の隙間が変化し, それに伴い椎骨間に存在している椎間板も変形する. 前湾が増加すると椎間関節への負荷が増加し, 逆に前湾が減尐すると椎間板への負荷が増加する. 椎骨間の隙間が狭くなると, 椎間板に大きな負荷が加わることから, 椎骨間の隙間変化がわかれば, 椎間板の変形が推測でき, 椎間板負荷が評価できると考えられる. 図 1.1 脊柱 (2) 3

5 第 2 章腰痛診断の現状 2.1 腰痛症の種類症状の多い腰痛症の種類を簡単に挙げる. 腰椎椎間板ヘルニア腰部脊柱管狭窄症腰椎分離症腰椎すべり症姿勢性腰痛症急性腰痛症 ( ぎっくり腰 ) 変形性腰椎症腰椎捻挫腰椎椎間板ヘルニアや脊柱管狭窄症は日常ごく普通に遭遇し, 頻度の高い疾患 (3) であり, その発症概念も明瞭である. 姿勢性腰痛症や変形性腰椎症などは現在, 疾患概念, 診断基準などに多尐不明確な部分があり, 診断には注意を要する. また, 運動麻痺を伴う腰痛は緊急度の高い腰痛と分類 (4) でき, 腰椎椎間板ヘルニアがこれにあたる. 腰痛の原因には脊柱構造の異常から内臓疾患が起因と様々だが, ほとんどの場合脊柱構造の異常である. また, 手術を要し腰痛を治療できるのは僅かであり, ほとんどの症状は保存治療や姿勢の改善により痛みを取り除くことができる. このことから姿勢と腰痛には大きな関係があると考えられる. 2.2 従来の姿勢変化の計測方法現在行われている画像による診断方法, また姿勢変化による椎間板への負荷を直接計測する方法を示す. X 線撮影 ( レントゲン撮影 ) (2) X 線撮影では,X 線照射装置とフィルムの間に体を置き, 焼き付けて画像化する.X 線撮影には種々な特殊撮影があるが, 単純 X 線撮影が画像診断の基本であり, 最も重要である. 骨組織を扱う整形外科では最も基本的な検査であり, 骨の形態や椎体骨梁の微細な変化や骨密度などの観察に適している. CT( コンピュータ断層撮影 ) (2) 人体組織の X 線吸収に差があることを利用して細かい X 線をビーム状に体軸に水平に照射し, 高感度の検出器で検出して X 線吸収係数をコンピュータで画像に再構築したものである. 特に横断面での広がりを見るのに適する. 4

6 MRI( 核磁気共鳴画像法 ) (2) 脊椎, 脊髄疾患の診断に関しては, 単純 X 線撮影に次ぐ検査として有用である. 原理としては, 核磁気共鳴現象を利用して生体内の内部情報を画像化する方法である. 断層画像という点では X 線,CT と一見よく似た画像が得られるが,CT とは全く異なる物質の物理的性質に着目した撮影法であるゆえに,CT で得られない情報が多く得られる.CT では主に骨に関する情報を得られるのに対し,MRI では, 神経, 椎間板, 骨髄の情報を低侵襲で得ることができる, 近年 3D-MRI(MRI ミエログラフィー ) によって神経根の連続した描出も可能になっている. 腰椎椎間板ヘルニア, 脊柱管狭窄症, 脊椎, 脊髄腫瘍や炎症性疾患のスクリーニングとして優れている. 現在, これらの方法が多くの腰痛診断で利用されており, 体の内部情報を得るには非常に適している. しかし,X 線撮影では放射線を浴びるため, 人体への影響から実験には制約があり, 姿勢変化を計測するには適していない. また,CT や MRI は基本静止した状態での計測となるため, 限られた姿勢での計測しか行えないという欠点が挙げられる. 椎間板内圧計測法世界的に有名なスウェーデンの整形外科医である Nachemson らが行った研究で第 3 腰椎と第 4 腰椎の間の椎間板に直接電極針を挿入して内圧を計測した方法がある. 椎間板内圧の計測結果を図 2.1 に示す. この方法であれば, 直接椎間板に加わる圧力を測定するため高い精度が得られるため信頼性が高く, 姿勢が変化していても測定可能である. しかし, 外科手術を要する侵襲的な計測であるため, 限られた研究機関でしか行えない. 身体を傷つける危険性があり, 動作中の計測は危険を伴うため一般的ではない. 図 2.1 椎間板内圧の計測結果 (2) 5

第 3 章提案する推定法の概要本研究ではより安全で動作中の計測も可能である 3 次元動作解析装置を用いる. 図 3.1 に 3 次元動作解析装置を示す. この装置は動作を計測したい対象にマーカを貼り付け, 複数台のカメラで撮影することにより,3 次元座標を取得することができる装置である.

7 第 3 章提案する推定法の概要本研究ではより安全で動作中の計測も可能である 3 次元動作解析装置を用いる. 図 3.1 に 3 次元動作解析装置を示す. この装置は動作を計測したい対象にマーカを貼り付け, 複数台のカメラで撮影することにより,3 次元座標を取得することができる装置である. 本研究では緊急性の高い腰痛であり, 疾患頻度が高く, 発症メカニズムも明瞭である腰椎椎間板ヘルニアを研究対象とする. 腰椎部の椎骨端部が背中の体表面に近い部分に存在していることに着目する. 図 3.2 に示すように腰椎部の椎骨端部に対応する位置にマーカを貼り付け, 複数台のカメラによりマーカを撮影する. 撮影された画像データを信号処理し, それぞれのマーカの 3 次元座標を計算する. それを信号処理することにより腰椎部の姿勢変化を計測する. 3 次元動作解析装置を用いる場合には被験者の動作による皮膚の微妙な動きがマーカに影響を及ぼすため, 腰椎部の全体的な姿勢変化から椎骨間の平均的な隙間変化を推定することを考える.3 次元動作解析装置により得られた各マーカの座標から腰椎部全体としての姿勢変化を最小 2 乗法を用いて円にフィットさせ, その曲率半径から腰椎部全体の平均的な隙間変化を推定し, 椎間板負荷を評価する. 図次元動作解析装置 6

8 図次元動作解析装置を用いた計測法の概要図 7

9 第 4 章 3 次元動作解析装置による計測の精度検証 4.1 精度検証の実験内容本研究において椎間板にかかる負荷を精度よく推定するには,3 次元動作解析装置により腰椎部の形状を精度良く推定できることが重要となる. そこで, 曲率半径が既知の計測モデルを図 4.1 に示すように配置し,3 次元動作解析装置を用いて精度検証を行う. 図 4.2 に示す計測モデルは曲率半径が 200mm,400mm,600mm,800mm,1000mmにカットした木材であり, これを用いた検証となる. 曲率半径の算出方法としては円の最小 2 乗法 ( 第 5 章 ) を用いる. 図 4.1 精度検証における実験概要図曲率半径 200mm 曲率半径 400mm 8

10 曲率半径 600mm 曲率半径 800mm 曲率半径 10000mm 図 4.2 計測モデル 9

腰椎部の湾曲の円弧と同等の長さである 170mmになるよう等間隔にマーカを 5 個貼り付ける. その様子を図 4.3 に示す.

11 直径の異なるマーカでの精度検証図 4.2 に示した計測モデルを用いて,3 次元動作解析装置で測定する. 実際の木材の曲率半径と最小 2 乗法により算出した曲率半径を比較する. 腰椎部の湾曲の円弧と同等の長さである 170mmになるよう等間隔にマーカを 5 個貼り付ける. その様子を図 4.3 に示す. また, 直径の異なるマーカを使用することにより算出した曲率半径に違いがあるか検証する. 図 4.4 に検証に使用したマーカを示す. 左から直径 4mm,6mm,10mmのマーカである. 図 4.3 円弧長さが 170mm のマーカ貼り付け位置図 4.4 検証に使用したマーカ 10

6 に示すように対象物に対して高さが異なるカメラを配置した 3 次元的に撮影した場合 ( 条件 2 とする ) を比較する.

12 カメラ配置による精度検証測定に用いるカメラの配置の違いにより算出する曲率半径に影響があるかを検証する. 図 4.5 に示すような対象物に対してカメラの高さを 4 台とも一定にした 2 次元的に撮影する場合 ( 条件 1 とする ). 図 4.6 に示すように対象物に対して高さが異なるカメラを配置した 3 次元的に撮影した場合 ( 条件 2 とする ) を比較する. 円弧の長さは 170mm になるようにマーカを 5 個貼り付ける. 図次元的に撮影する場合のカメラ配置 ( 条件 1) 図次元的に撮影する場合のカメラ配置 ( 条件 2) 11

13 一致率 [%] 4.2 精度検証の実験結果表 4.1 に最小 2 乗法により算出した曲率半径をマーカの直径別に示し, 図 4.7 に一致率を示す. なお, 一致率は正確値に対する算出値の比をパーセンテージで表すこととする. 表 4.1 直径の異なるマーカを使用した場合の曲率半径算出結果木材の曲率半径曲率半径 ( マーカ直径 4mm) 曲率半径 ( マーカ直径 6mm) 曲率半径 ( マーカ直径 10mm) [mm] マーカ直径 4mm 6mm 10mm 木材の曲率半径 [mm] 図 4.7 直径の異なるマーカを使用した場合の一致率 12

14 表 4.2 に計測に用いるカメラの配置を変更した場合での算出した各曲率半径を示し, 図 4.8 に一致率を示す. 表 4.2 カメラ配置の条件を変えた場合の曲率半径算出結果木材の曲率半径曲率半径 ( 条件 1) 曲率半径 ( 条件 2) [mm] 図 4.8 カメラ配置の条件を変えた場合の一致率 13

15 曲率半径が大きくなるにつれて実際の曲率半径と算出した曲率半径との誤差が大きくなる. これは図 4.9 に示すように曲率半径が大きくなるということは, 円周の長さに対する測定した円弧長さの比が減尐するため, マーカの位置座標を読み取る精度が大きく影響したと考える. マーカの直径の違いによる一致率を比較すると, 直径が小さくなるにつれて一致率は下がる傾向がある. これはキャプチャー領域に対してマーカが小さくなると, 位置座標を読み取りにくくなることが原因であると考えられる. しかし, 撮影時におけるマーカの干渉, 椎骨端部への貼り付けなどを考えると 4mm または 6mm が理想的である. 精度検証の両マーカの結果より, 本研究では木材の曲率半径との一致率が高い 6mm のマーカを使用することとする. カメラ配置の条件を変更したことによる精度の大きな向上は見ることができなかった. マーカを立体的に撮影するという観点から高さが異なるカメラを配置した条件 2 の場合が若干ながら精度が良い. この結果から条件 2 のカメラ配置でこれ以降の実験を行う. 椎間板負荷の推定手順については第 5 章で詳しく述べるが, 本手法を用いたことによる曲率半径と椎間板負荷の関係を図 4.10 に示す. 椎間板負荷は曲率半径が大きくなるにつれて, 一定に近づく. この図より曲率半径が大きい場合には実際の腰椎部の曲率半径と誤差が生じても, 椎間板負荷の推定には影響が尐ないと考えられる. 各精度実験の結果から決定した条件で,3 次元動作解析装置を用いて腰椎部の曲率半径を推定し, 第 6 章で述べる椎間板負荷の推定実験を行うこととする. 図 4.9 曲率半径別の円周に対する円弧長さ 14

16 図 4.10 曲率半径と椎間板負荷の関係 15

17 第 5 章提案する椎間板負荷の推定法 5.1 腰椎部の曲率半径の推定姿勢変化を計測するために腰椎部の曲率半径を推定する.3 次元動作解析装置により得られた 5 個のマーカの位置座標データから最小 2 乗法を用いて円にフィットさせる. 円の最小 2 乗法 (5) の手順を以下に示す. 円の一般式は中心座標を a, b, 半径を r とすると x a y b r (5.1) 式 (5.1) を = 0の式に変形し, 点の座標 x を代入し, その式の 2 乗の総和を求める. 2 x a y b 2 2, i y i 2 i r 0 (5.2) i このまま未知数 a,b,r に関して式 (5.2) を偏微分しても a,b,r は 4 次関数となるので, 未知数の偏微分 = 0 の時が最小とは限らない. そこで, 式 (5.2) を展開し, 次式のように置く. ただし, 2 2 i 2 xi y Ax By C 0 (5.3) i i A 2a (5.4) B 2b (5.5) C (5.6) a b r 式 (5.3) に関しては A,B,C について偏微分すると A A A B A C xi B xi yi C xi xi xi yi xi yi B yi C yi xi yi yi 2 2 xi B yi C1 xi yi (5.7) (5.8) (5.9) 式 (5.7)(5.8)(5.9) を解くことによりA,B,Cが求まり, 式 (5.4) (5.5) (5.6) より円の中心座標 a, b, 円の半径 r を求めることができる. 16

18 5.2 椎骨間の隙間の推定腰椎部の曲率半径を推定後, 椎骨間の隙間を推定する. 図 5.1 に腰椎部の概要を示す. 椎骨は力学的にはピンジョイントと等価な形で結合されていると考える.2 次元的な姿勢変化ではピンジョイント部分の回転と椎骨関節での弾性変形により決まると考える. 本研究では椎間関節での弾性変形は微小な変形であるため, 無視する. ここでは, 図のように円弧が外側に凹の場合での椎間板負荷の推定法について説明する. 図 5.1 椎骨間の隙間の推定図 5.1 でマーカの中心からジョイント部の位置までの距離 d, ジョイント部の位置から椎間板の中心位置までの距離を e, 円の最小 2 乗法で推定した腰椎部の曲率半径を r とする.d,e はどの椎骨においても一定とする. ジョイント部の位置での椎骨の幅を v とする. 本研究では, 椎間関節の弾性変形を無視しているので一定とする. 椎間板の位置での椎骨の幅を w とする. このような仮定をおき, ジョイント部までの長さ r v を r v r d (5.10) と定義すれば, r v および v を用いて 1 個の椎骨に対応する角度 v が 17

19 v 2sin 1 v 2rv (5.11) より得られる. この式は曲率半径が小さくなれば, 角度 v が大きくなり, 椎間板の位置での隙間が広がること, すなわち圧縮力が小さくなることを意味している. 逆に曲率半径が大きくなれば, 角度 v が小さくなり, 椎間板の位置での隙間が小さくなること, すなわち圧縮力が大きくなることを意味している. 次に椎間板位置までの長さ r w を r w r d e (5.12) と定義すれば, r w および w を用いて椎間板位置での椎骨の幅に対応する角度 w が w 2sin 2 w r w (5.13) より得られる.2 つの角度の差より椎間板位置に対応する角度を v w (5.14) により計算し, 椎間板位置の半径である r w を乗じて r w (5.15) を計算すれば, 椎骨間の隙間の円弧長さを表現することができる. また, 円弧が外側に凸の場合は円弧の中心が腰椎よりも前方になることから d,e の符号が反対になる. したがって, 円弧が外側に凹の場合に展開した推定手順に関する式において, d を d,e を e に置き換え, 外側に凹の推定手順を用いれば椎間板負荷を評価できる. 18

20 5.3 椎間板負荷の推定椎骨間の隙間を推定後, 椎間板負荷の推定を行う. 椎間板に負荷がかかっていない状態での曲率半径を r 0 とすればそれに対応する 0 が得られ, 無負荷状態から計測した状態への椎間板位置での隙間変化は 0 (5.16) となる. 最も簡単な場合として, 椎間板の負荷が隙間変化に線形のばね定数を乗じたもので表現できるのであれば, に対応する等価なばね定数を k と置けば, 負荷 p は p k (5.17) で表現できる. ここで, 立位姿勢の場合に添え字 s をつけて表現すれば, p s k s (5.18) となる. p を p s で正規化すれば, 無次元の負荷が k k s s (5.19) で表現される. したがって, かなりの仮定があるが正規化した負荷であれば, ばね定数の大きさに影響を受けないので, 隙間に関する情報のみで椎間板負荷を推定することが可能となる. 19

静止した姿勢, 動作中の姿勢において測定を行い, 椎間板負荷を推定する. 図 6.1 に実験風景を, 図 6.2 に実験時のマーカの貼り付け位置を示す.

21 第 6 章椎間板負荷の推定実験 6.1 椎間板負荷の推定実験内容今回は 2 次元的な姿勢変化で重量物を持たない場合を対象とした測定を行う. 被験者の腰椎部の椎骨端部に対応する位置にマーカを貼り付け,4 台のカメラを用いて測定する. また, この実験におけるサンプリング周波数は 100Hz である. 静止した姿勢, 動作中の姿勢において測定を行い, 椎間板負荷を推定する. 図 6.1 に実験風景を, 図 6.2 に実験時のマーカの貼り付け位置を示す. 実験には 3 人の被験者を用いる.3 人の被験者には実験内容を説明し, 同意を得た上で協力して頂いた. 表 6.1 に各被験者の身体情報を示す. 図 6.1 実験風景図 6.2 マーカの貼り付け位置 20

表 6.1 被験者の身体データ性別身長 (cm) 体重 (kg) BMI 体脂肪率 (%) 備考被験者 1 男 168 58.6 20.7 20.4 腰痛あり被験者 2 男 167 73.3 26.3 27 腰痛なし被験者 3 男 181 73.2 22.4 15.7 腰痛なし実験時に使用する装置本研究での椎間板負荷の推定実験において, 上体の前屈角度により測定する姿勢を決定する.

22 表 6.1 被験者の身体データ性別身長 (cm) 体重 (kg) BMI 体脂肪率 (%) 備考被験者 1 男腰痛あり被験者 2 男腰痛なし被験者 3 男腰痛なし実験時に使用する装置本研究での椎間板負荷の推定実験において, 上体の前屈角度により測定する姿勢を決定する. そのため, 上体を前傾した時の前屈角度を測る必要がある. これまでは医療用の角度計とデジタル角度計を併用することで上体の前屈角度を測ってきた. しかし, これでは動作中の測定時に前屈角度を測ることができない. そこで今回の椎間板負荷の推定実験では図 6.3 に示す Xsens 社製の姿勢センサ MTx を用いてリアルタイムに前屈角度を測定する. 姿勢センサでの測定と 3 次元動作解析装置の撮影を同期させ, 以降に述べる姿勢時の実験を行う. 図 6.3 姿勢センサ 21

4 静止時の測定の測定姿勢動作中の測定動作中の測定では, 低速動作と高速動作の測定を行う.

23 静止時の測定重量物を持たない静止時の姿勢の測定を行う.Nachemson が実験で行った測定姿勢の内, 図 6.4 に示す 4 つ姿勢の測定を行う. 各姿勢において静止状態 (1 秒間 ) の姿勢を測定する. 図 6.4 静止時の測定の測定姿勢動作中の測定動作中の測定では, 低速動作と高速動作の測定を行う. 低速動作は楽に立った姿勢から前屈を行う動作を測定する. 計測開始から 2 秒間は静止する.2 秒経過後から前屈動作を行う. 撮影時間は 7 秒間とする 0 秒 2 秒 7 秒図 6.5 立位姿勢からの前屈動作 22

24 高速動作では咳を行った動作の測定する. 立位姿勢をとり, 測定開始 4 秒後に 2 回,7 秒後に 1 回咳を行う. 撮影時間は 10 秒間とする. 0 秒 4 秒 7 秒 10 秒図 6.6 立位姿勢からの咳動作 23

25 6.2 椎間板負荷の推定結果静止時の測定での椎間板負荷の推定結果椎間板負荷の推定に用いる腰椎部の曲率半径は 1 秒間の平均値を用いる. 指標とする Nachemson と同様に楽に立った姿勢で正規化を行い, 椎間板負荷を評価する. 各被験者の椎間板の物性値を調査するのは困難であるため, 比較対象である Nachemson の実験結果のスケールに合わせるように決定した. そのため, 今回は無負荷状態の隙間部分の円弧長さ 0 は各被験者に対して仮定の値を用いる. そして, 椎間板負荷の推定結果と Nachemson の実験結果と比較する. 表 6.2 に各被験者の腰椎部の曲率半径を推定した結果を示す. 表 6.2 各被験者の腰椎部の曲率半径測定姿勢被験者 1 被験者 2 被験者 [mm] 表 6.3 に各被験者に使用した無負荷状態での隙間部分の円弧長さ 0 の仮定値を示す. 表 6.3 各被験者の無負荷状態での仮定値無負荷状態での隙間部分の円弧長さ 0 被験者被験者被験者 [mm] 図 6.6 から図 6.8 に各被験者の椎間板負荷の推定結果を示す.Nachemson の実験結果と同様に楽に立った姿勢を基準姿勢 (100) とし, 各姿勢の椎間板負荷を示す. 24

26 図 6.6 静止時の椎間板負荷比 ( 被験者 1) 図 6.7 静止時の椎間板負荷比 ( 被験者 2) 25

27 図 6.8 静止時の椎間板負荷比 ( 被験者 3) 図 6.9 に Nachemson の実験結果と本手法による実験結果を比較したものを示す. 左端に Nachemson の実験結果, 右の 3 つが本手法による被験者の椎間板負荷の推定結果である. 図 6.9 Nachemson の実験結果と本手法による実験結果の比較 26

28 表 6.4 に Nachemson の実験結果と本手法による椎間板負荷の推定結果の詳細な数値を示す. 表 6.4 Nachemosn の実験結果と本手法による椎間板負荷の推定結果の数値計測姿勢 Nachemson 被験者 1 被験者 2 被験者 [%] 比較対象である Nachemosn の実験結果と本手法による椎間板負荷の推定結果を比較したところ, 全ての被験者において同様の傾向が得られた. 椎間板負荷の詳細な数値で見ると多尐の誤差があるが,Nachemson が用いた被験者で測定を行っているわけではないため, 厳密な精度に関する比較は難しい. しかし, 各姿勢の傾向が見られたことから, かなり一致した結果が得られたと考える. 測定した 4 つの姿勢での椎間板負荷の推定実験は継続して行われてきており,Nachemson の実験結果と同様の傾向は以前からも得られている. そこで,Nachemson の被験者の数値で精度の比較は難しいが, 今回協力して頂いた各被験者の椎間板負荷の推定結果を比較対象とすれば, 精度に関する議論が可能だと考える. 次に述べる動作中の測定では表 6.4 に示した各被験者の推定結果を比較対象とし, 同角度での椎間板負荷を見る. 27

29 動作中の測定での推定結果前屈動作では撮影開始から静止していた 0 秒から 2 秒時の平均的な姿勢を基準姿勢 (100) として, 姿勢変化による椎間板負荷を評価する. 図 6.10, 図 6.11 に楽に立った姿勢から前屈動作を行った被験者 1, 被験者 3 の椎間板負荷の推定結果を示す. この結果は上体の前屈角度に対する椎間板負荷の推定結果である. また, 表 6.4 に示した推定結果と比較するために推定した椎間板負荷を近似曲線で示している. 図 6.10 前屈動作における椎間板負荷の推定結果 ( 被験者 1) 図 6.11 前屈動作における椎間板負荷の推定結果 ( 被験者 3) 28

30 表 6.5 に 2 人の被験者の前屈角度が 20 時の静止時と動作時の負荷推定結果の数値を示す. 表 6.5 前屈角度 20 における静止時と動作時の負荷推定結果静止時の負荷推定結果動作時の負荷推定結果被験者被験者人の被験者から椎間板負荷の変化を見ることができた. 推定結果は両者とも前屈角度が増加していくと椎間板負荷も増加する. 静止時の実験結果からも明らかであり, その傾向が得られた. 前屈角度が 20 の場合での推定した椎間板負荷は静止時の実験, 前屈動作の実験ともにかなり一致した数値を得ることができた. 前屈動作は比較対象がなく, 精度に関する比較は行えなかったが, 姿勢センサを用いたこの実験方法ならば同被験者で精度の比較は可能と考えられる. また, 0 から 30 まで 2 人の被験者ともに椎間板負荷は線形が保たれている. このことから前屈動作の推定精度が明確になれば, 前屈角度が 20 以外の姿勢にも新しい指標ができると考えられる. 29

31 咳動作では撮影開始時の姿勢を基準姿勢 (100) として, 姿勢変化による椎間板負荷を評価する. 図 6.12 から図 6.14 に楽に立った姿勢から咳動作を行った被験者 3 人の椎間板負荷の推定結果を示す. 図 6.12 咳動作における椎間板負荷の推定結果 ( 被験者 1) 図 6.13 咳動作における椎間板負荷の推定結果 ( 被験者 2) 30

32 図 6.14 咳動作における椎間板負荷の推定結果 ( 被験者 3) 咳動作測定では全ての被験者から咳を行ったと同時に椎間板負荷の変化を見られた. 咳を行うことで腰痛を患う例は多数報告されており, 症状の悪化の原因でもあるといわれている. このことから咳を行うことで椎間板に加わる負荷は増加すると考えられる. 今回の推定結果においてその傾向を得ることができた. 精度に関する議論はこれからの課題であるが, 咳のような高速な連続動作でも本手法は椎間板負荷を推定するのに対応可能である見通しを得られた. 31

33 第 7 章結言本研究では, 身体に負担をかけず比較的容易に椎間板負荷を測定する方法として 3 次元動作解析装置を用いた非侵襲的な推定方法を提案した. 本手法は腰椎部の椎骨端部が背中の体表面に近い部分に着目し, 椎骨の動きを 3 次元動作解析装置を用いて外部から測定する. その動きから椎骨間の隙間変化を推定し, 椎間板負荷として評価をした. 静止時の測定においては本手法による推定結果と侵襲的な測定である Nachemson の実験結果を比較し, 傾向としてはかなり一致した結果を得られた. 動作中の測定では前屈動作, 咳動作の椎間板負荷の推定をし, 動作による負荷の変化を見ることができた. 前屈動作での推定結果は前屈角度が増加するにつれて椎間板負荷も増加することが確認できた. また, 各被験者の静止時の推定結果を新たな指標とし, 精度を評価したところかなり一致した値が得られた. 咳動作では咳を行うと椎間板負荷が増加する傾向が推定結果から得ることができ, 高速な動作でも対応可能であることが確認できた. 今回の実験結果により, 本手法が椎間板負荷を推定するにあたり有効であることがわかった. 今後は, ここで示した測定姿勢についても検討を行う必要があるが, 姿勢が複雑になればなるほど比較対象がないため, 様々な姿勢での椎間板負荷を推定しても評価が難しいのが現状である. 単純な姿勢や動作であるが, 被験者を増やしデータを蓄積することで新たな指標を生み出し, 本手法の有用性をより確かなものにしていく必要があると考える. 32

34 謝辞本研究を進めるにあたり, 井上喜雄教授ならびに芝田京子准教授に多大なご指導賜り, 心から感謝の意を表したいと思います. また, 共に研究に励んだ研究室の方々に感謝致します. 33

35 参考文献 (1)LUMBAR DISC PRESSURE AND MYOELECTRIC BACK MUSCLE ACTIVITY DURING SITTIHG: I.Studies on an Experimental Chair B.J.G.Andersson, R.Örtengren, A.Nachemson, and G.Elfström,1974,pp (2) 最新腰痛症ハンドブック腰椎椎間板ヘルニアからスポーツ, 事故の治療まで : 遠藤健司, 金岡恒治,2008,pp7-81. (3) 整形外科 : 痛みへのアプローチ腰背部の痛み : 寺山和男, 片山治, 山本博司,2006,pp (4) 腰痛知る診る治す : 清水克時,2008,pp26-58 (5) 一般式による最小二乗法 ( 円の最小二乗法 ) 画像ソリューション ( 参照日 : ) 34

画像類似度測定の初歩的な手法の検証

画像類似度測定の初歩的な手法の検証島根大学総合理工学部数理情報システム学科計算機科学講座田中研究室 S539 森瀧昌志 1 目次第 1 章序論第章画像間類似度測定の初歩的な手法について.1 A. 画素値の平均を用いる手法.. 画素値のヒストグラムを用いる手法.3 C. 相関係数を用いる手法.4 D. 解像度を合わせる手法.5 E. 振れ幅のヒストグラムを用いる手法.6 F. 周波数ごとの振れ幅を比較する手法第

目次 第 1 章緒言 研究背景 研究目的 脊柱の構造... 3 第 2 章腰痛診断の現状 腰痛症の種類 従来の姿勢変化の計測方法... 4 第 3 章提案する推定法の概要... 6 第 4 章精度検証... 8

目次第 1 章緒言研究背景研究目的脊柱の構造... 3 第 2 章腰痛診断の現状腰痛症の種類従来の姿勢変化の計測方法... 4 第 3 章提案する推定法の概要... 6 第 4 章精度検証... 8