(57)【要約】【課題】最適なブラッシング動作をコンピュータシミュレーションにより導出し、可視化することができるブラッシング動作のシミュレーション方法を提供すること。【解決手段】本発明のブラッシング動作のシミュレーション方法は、本発明のブラッシング動作のシミュレーション方法は、ブラッシングステップ１と、ブラッシング計測ステップ２と、評価パラメータ計測ステップ３と、正規化処理ステップ４と、主成分分析ステップ５と、主成分抽出ステップ６と、重回帰分析ステップ７と、最適化処理ステップ８と、ブラッシング動作生成ステップ９と、可視化ステップ１０とを有することを特徴とする。【選択図】        図１
2019030578.tif 000002
【特許請求の範囲】
【請求項１】
歯ブラシを用いて特定の口腔部位に対しブラッシング動作を行わせるブラッシングステップと、
前記ブラッシングステップにおける前記歯ブラシの位置と姿勢を計測するブラッシング計測ステップと、
前記ブラッシングステップにおける前記口腔部位に対する負荷パラメータを計測する負荷パラメータ計測ステップと、
前記ブラッシングステップの後における前記口腔部位に対する評価パラメータを計測する評価パラメータ計測ステップと、
前記ブラッシング計測ステップで計測した計測データを１ストローク毎に正規化する正規化処理ステップと、
前記正規化処理ステップで処理した正規化データに対して主成分分析を行う主成分分析ステップと、
前記主成分分析ステップで抽出された主成分の中から一部又は全ての前記主成分を解析対象として抽出する主成分抽出ステップと、
前記主成分抽出ステップで抽出した前記主成分の主成分得点から、前記負荷パラメータと前記評価パラメータを予測するための重回帰分析を行う重回帰分析ステップと、
前記重回帰分析ステップで得られた重回帰式を用いて目的関数を設定し、前記目的関数を最大化又は最小化する主成分得点を導出する最適化処理ステップと、
前記最適化処理ステップで導出された前記主成分得点から、最適ブラッシング動作を再構築するブラッシング動作生成ステップと、
前記ブラッシング動作生成ステップで再構築された前記最適ブラッシング動作を表示する可視化ステップと
を有する
ことを特徴とするブラッシング動作のシミュレーション方法。
【請求項２】
前記口腔部位を歯面又は歯肉とし、
前記負荷パラメータを、前記歯ブラシを介して前記歯面又は前記歯肉に作用させるブラッシング力とし、前記評価パラメータを、前記ブラッシングステップの前後における歯垢等のデンタルプラークや歯石の除去効果を評価するプラーク等除去率とした
ことを特徴とする請求項１に記載のブラッシング動作のシミュレーション方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、最適な口腔ケアスキルを導出・提示することができるブラッシング動作のシミュレーション方法に関する。
【背景技術】
【０００２】
口腔内の常在菌により形成されるデンタルプラークは、う蝕や歯周病の原因である他、菌が歯周ポケットを通じ日常的に体内へ侵入することで心内膜炎や狭心症のリスクを高め、さらにその炎症因子によって糖尿病等の生活習慣病を助長することが知られている。加えて、高齢者がデンタルプラークによる誤嚥性肺炎で死に至るなど、全身への感染リスクが見過ごせないものとなっている。デンタルプラークは強固に歯面に粘着しており、機械的なブラッシングによる除去が不可欠である。
ブラッシングの位置、姿勢、部位、時間、回数、周期、力など、ブラッシングの状態を検出し、確認できる装置が提案されている（例えば特許文献１から３）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００５−１５２２１７号公報
【特許文献２】特開２０１７−０６０６６１号公報
【特許文献３】特開２０１７−０６０６６２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかし、ブラッシングスキルを正確に評価する装置や、エビデンスに基づいた最適なブラッシング手技を獲得する方法が確立されていない。
ここで、最適なブラッシングとは、例えば、歯ブラシによる歯表面（歯面）の摩耗や歯肉の損傷を可及的に減じるために、歯面に対して歯ブラシから加わる力、すなわちブラッシング力を最小にしつつ、プラーク除去率を最大化するような歯ブラシの動かし方と定義できる。あるいはこれとは逆に、ブラッシング力を最大にしつつ、プラーク除去率を最小化するような悪条件に最も適した歯ブラシの動かし方を避けること、とも定義することもできる。歯面や歯肉へ為害作用を及ぼさない範囲内を条件に付与することもできる。
【０００５】
本発明は、このような目的に対して、最適なブラッシング動作をコンピュータシミュレーションにより導出し、可視化し、これをフィードバックすることができるブラッシング動作のシミュレーション方法を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明の第１の実施の形態によるブラッシング動作のシミュレーション方法は、歯ブラシを用いて特定の口腔部位に対しブラッシング動作を行わせるブラッシングステップと、次いで計測ステップとして、ブラッシングステップにおける歯ブラシの位置と姿勢を計測するブラッシング計測ステップ、ブラッシングステップにおける口腔部位に対する負荷パラメータを計測する負荷パラメータ計測ステップ、ブラッシングステップの後における口腔部位に対する評価パラメータを計測する評価パラメータ計測ステップがまず存在する。さらに計算ステップとして、ブラッシング計測ステップで計測した計測データを１ストローク毎に正規化する正規化処理ステップ、正規化処理ステップで処理した正規化データに対して主成分分析を行う主成分分析ステップ、主成分分析ステップで抽出された主成分の中から一部又は全ての主成分を解析対象として抽出する主成分抽出ステップと、主成分抽出ステップで抽出した主成分の主成分得点から、負荷パラメータと評価パラメータを予測するための重回帰分析を行う重回帰分析ステップ、重回帰分析ステップで得られた重回帰式を用いて目的関数を設定し、目的関数を最大化又は最小化する主成分得点を導出する最適化処理ステップ、最適化処理ステップで導出された主成分得点から、最適ブラッシング動作を再構築するブラッシング動作生成ステップがあり、最後に、ブラッシング動作生成ステップで再構築された最適ブラッシング動作を表示する可視化ステップとを有するものである。
本実施の形態によれば、ブラッシングステップでは、例えばブラッシングの位置、角度、力、速度、振幅、部位、被験者などを変えることで、様々なブラッシングを計測することができる。そして、ブラッシング計測ステップにおける計測データを正規化し、正規化した複数の条件の正規化データに対して、主成分分析ステップで運動の分解を行い、分解した運動（主成分得点）を変数として、負荷パラメータと評価パラメータを予測する重回帰式をそれぞれ求め、最適化処理ステップで目的関数を設定し、その目的を満足する最適な変数（主成分得点）を最適計算で求め、得られた最適値から、動作の再構築と可視化を行うことにより、予め計測した計測データとは異なる最適なブラッシング動作を生成することができる。
【０００７】
本発明の第２の実施の形態は、第１の実施の形態によるブラッシング動作のシミュレーション方法において、口腔部位を歯面又は歯肉とし、負荷パラメータを、歯ブラシを介して歯面又は歯肉に作用させるブラッシング力とし、評価パラメータを、前記ブラッシングステップの前後における歯垢等のデンタルプラークや歯石の除去効果を評価するプラーク等除去率としたものである。
本実施の形態によれば、歯面や歯肉への障害にもなり得るブラッシング力を最少にしつつ、プラーク等除去率を最大化にできるブラッシング動作を生成することができる。
【発明の効果】
【０００８】
本発明のブラッシング動作のシミュレーション方法によれば、最適なブラッシング動作をコンピュータシミュレーションにより導出し、可視化することができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の一実施例によるブラッシング動作のシミュレーション方法を示す概念図
【図２】本実施例に用いた上顎歯列模型、模型歯、及び歯ブラシを示す図
【図３】磨き方パターンを示す図
【図４】プラークの付着状態を示す写真
【図５】正規化処理を示すグラフ
【図６】主成分分析結果を示す図
【図７】シミュレーションでＦｚｍｅａｎとＥｓを用いて最適化したブラッシング動作（実線）と被験者のブラッシング動作（破線）との比較を可視化した図
【図８】シミュレーションでＦｚｍｅａｎとＥｓを用いて最適化した動作と、最悪化した動作を再現して可視化した図
【発明を実施するための形態】
【００１０】
以下本発明の実施例について図面とともに説明する。
図１は本発明の一実施例によるブラッシング動作のシミュレーション方法を示す概念図である。
本実施例によるブラッシング動作のシミュレーション方法は、以下の処理手順で行う。
まず、測定空間内に基準座標系（xyz座標の原点と向き）を定め、その測定空間に、歯ブラシと歯列模型を設置する。測定データは、基準座標系を原点として、歯ブラシの位置と姿勢、及び歯列模型の位置と姿勢を計測する。なお、本実施例では、歯列模型はテーブルに固定され、静止状態にあるが、実際のブラッシング指導の際には、対象者の歯列の位置と姿勢が変化するため、その位置と姿勢に応じて、導出した最適な歯ブラシの位置と姿勢を座標変換してから可視化する必要がある。
歯ブラシを用いて特定の口腔部位である歯面や歯肉に対しブラッシング動作を行わせ（ブラッシングステップ１）、ブラッシングステップ１における歯ブラシの位置と姿勢を計測するとともに（ブラッシング計測ステップ２）、口腔部位である歯面や歯肉に対する負荷パラメータ（ブラッシング力）を計測する（負荷パラメータ計測ステップ３）。ブラッシングステップ１における計測は、対象者の歯型をとり（印象採得し），作成した歯列模型を用いることが好ましい。
また、ブラッシングステップ１の後における口腔部位に対する評価パラメータ（プラーク除去率）を計測する（評価パラメータ計測ステップ４）。ブラッシングステップ１において歯列模型でなく対象者の実際の口腔部位をブラッシングの対象とする場合は，評価パラメータ計測ステップ４では、例えば、対象者の口腔部位に対して，実際に歯垢染色剤で着色し、この着色の除去率を数値化し、又はUVなどを用いてブラッシングステップ１の前後における実際の歯垢や歯石を計測する。すなわち、評価パラメータ計測ステップ４では、ブラッシングステップ１の前後における歯垢等のデンタルプラークや歯石の除去効果を評価する。
ブラッシング計測ステップ２で計測した計測データは、１ストローク毎に正規化する（正規化処理ステップ５）。
正規化処理ステップ５で処理した正規化データに対して主成分分析を行う（主成分分析ステップ６）。主成分分析ステップ６で抽出された主成分の中から累積寄与率が所定値以上を満たす主成分を解析対象として抽出する（主成分抽出ステップ７）。なお、本実施例では、主成分の中から累積寄与率が所定値以上を満たす主成分を解析対象として抽出したが、全ての主成分を解析対象としてもよく、又は例えばブラッシング力及びプラーク等除去率の少なくとも一方と相関のある主成分のみを抽出、すなわち累積寄与率ではなく一部の主成分を抽出してもよい。
主成分抽出ステップ７で抽出した主成分の主成分得点から、負荷パラメータと評価パラメータを予測するための重回帰分析を行う（重回帰分析ステップ８）。
重回帰分析ステップ８で得られた重回帰式を用いて目的関数を設定し、目的関数を最大化又は最小化する主成分得点を導出する（最適化処理ステップ９）。なお、本実施例では、重回帰分析を用いたが、重回帰分析に代えて、予測式を作成できる他の多変量解析を用いることもできる。
最適化処理ステップ９で導出された主成分得点から、最適ブラッシング動作を再構築する（ブラッシング動作生成ステップ１０）。
ブラッシング動作生成ステップ１０で再構築された最適ブラッシング動作を表示する（可視化ステップ１１）。
正規化処理ステップ５から可視化ステップ１１は、コンピュータによって行わせることができる。
【実施例】
【００１１】
以下に、負荷パラメータを、歯ブラシを介して前記歯面又は前記歯肉に作用させるブラッシング力とし、評価パラメータを、プラーク除去率として本発明の実施例を説明する。
【００１２】
図２は、本実施例に用いた上顎歯列模型、模型歯、及び歯ブラシを示している。
図２（ａ）に示すように、上顎歯列模型１の模型歯２の左側中切歯唇側を被験歯面とした。歯ブラシ３の歯ブラシヘッド部３ａには、ブラッシング力を計測する力覚センサ４を設けている。なお、歯ブラシ３の位置と姿勢（３次元運動）は、磁気式３次元位置・姿勢計測システムで計測した（図示を省略）。ブラッシング力は、歯ブラシヘッド部３ａを加工し、許容値 ＦＸ＝ＦＹ＝±２．ＯＮ、ＦＺ＝２．ＯＮの小型３軸力覚センサ（ＳＰ２２−ＦＦＣ１５、Ｔｏｕｃｈｅｎｃｅ  Ｉｎｃ．）を２個組み込むことで、ブラッシング時の合力ベクトルＦと作用点Ｌを計測した（図２（ｂ））。
３名を被験者とし、実験には図２に示す上顎歯列模型１、模型歯２、そして人工プラークを用いた。
まず、人工プラークの付着量に差が生じないようにエタノールで模型歯２を拭き、歯面の前面に均一に塗布した。
【００１３】
次に、できるだけ多くのブラッシングパターンをデータベース化するために、力や振幅、速さ、磨く部位などが異なる種類の磨き方を被験者に指示した。
図３は、磨き方パターンを示している。歯頚部については、磨き方のパターンの中に動作を繰り返させることで複数回計測した。これは、人のブラッシング動作がある程度ばらつくことが予測されること、そして、歯頚部のプラーク除去率に焦点を当てて、後に最適化計算を行うためである。
ストローク回数は各磨き方に対して１６回と設定した。ブラッシング後はプラーク除去率を求めるために歯面の撮影を行った。撮影は一定の環境で行うために、LEDライトが固定された遮光BOXの中で、デジタルカメラと被写体との位置を１０５ｍｍ、ＩＳＯ感度を８０、画素数を１４１０万画素で行った。写真には画像処理ソフト（ＩｍａｇｅＪ）を用いて画像処理を行うことでプラーク除去率を算出した。
【００１４】
図４はプラークの付着状態を示す写真である。
図４（ａ）に示すように、プラーク除去率は、歯頚部と、中央部と、それらの合計との３種類とした。図４（ｂ）はブラッシング後のプラークの付着状態を示している。
【００１５】
図５は正規化処理を示している。
変数は、歯ブラシの座標（ｘ、ｙ、ｚ）、角度（α、β、γ）、座標と角度を微分し求めた速度（ｘ’、ｙ’、ｚ’）、加速度（ｘ’’、ｙ’’、ｚ’’）、角速度（α’、β’、γ’）、角加速度（α’’、β’’、γ’’）の１８パラメータ、ブラッシング時間（ｔ）の１パラメータからなる合計１９パラメータとし、これらを１ストロークごとに切り出して平均化を行い、０から１００%まで１%刻み（１０１フレーム）で表示できるように正規化を行い、主成分分析を行った。
１９１９変数（１９パラメータ×１０１フレーム）全てについて、式（１）を用いて平均０、分散１となるように正規化した。
【００１６】
【数１】
2019030578.tif 000003
【００１７】
ここで、Ｚｔは正規化された値、Ｘｔは実測値、μｔは加算平均、σｔは標準偏差である。
【００１８】
主成分分析ステップ６では、正規化処理ステップ５で処理した正規化データから６８試行×１９１９変数（１９パラメータ×１０１フレーム）の行列（この行列が後述する[データベースの値]に相当する）を作成し主成分分析を行った。主成分分析によって、ブラッシング運動を、別の視点で表現している。なお、図３に示す２５種類×３名＝75試行のデータのうち、16ストローク以上磨いてしまい解析対象から外れたデータを除く68試行を分析対象データとして用いた。
図６は、主成分分析結果を示し、主成分分析の寄与率、累積寄与率、また、各主成分（ＰＣＶ）毎のブラッシング力と除去率のそれぞれと主成分得点の相関を示している。
主成分の個数は累積寄与率＞９５%を満たす範囲とすると、主成分分析の結果より、１８個の主成分でブラッシングを表すことができる。
主成分分析を行うことで、任意の被験者データ（座標、角度、ブラッシング力など）から固有ベクトル行列を介してその被験者の主成分得点を求めることができる。また逆に、主成分得点から固有ベクトル行列を介してある被験者データ（座標、角度、ブラッシング力など）を求めることもできる。この特性を用いて最適な主成分得点を求めることで、最適なブラッシングを導き出せる。すなわち、[主成分得点]＝[固有ベクトル行列][データベースの値]の関係であり、実際に測定したデータベースに含まれていないブラッシング動作も導き出すことができる。
【００１９】
次に重回帰分析ステップ８と最適化処理ステップ９について説明する。
主成分分析によって得た主成分得点と、ブラッシング力、プラーク除去率に重回帰分析を行い、それぞれの重回帰係数を求める。
重回帰分析を用いて、主成分得点からＺ方向のブラッシング力の平均Ｆｚｍｅａｎと歯頚部のプラーク除去率Ｅｓを求めた（式２）。なお、負荷パラメータには、合力の大きさ、歯ブラシヘッド部に垂直な力の成分、接線方向成分，あるいは1ストロークにおける各成分の最大値、平均値、積算値などが考えられるが、今回はこのうちの垂直に押し付ける力の平均値を選択した。またプラーク除去率も、歯頚部、中央部、全体（それから隣接部）などがあるが、今回は歯頚部を選択した。
【００２０】
【数２】
2019030578.tif 000004
【００２１】
ここで、Eはプラーク除去率、Fはブラッシング力、 ＰＣＳは主成分得点、Ａはブラッシング力の重回帰係数、Ｂはプラーク除去率の重回帰係数、ｋは主成分数（ｋ＝１、２、・・・Ｎ）、Ｎ＝18，Ｃ0及びＤ0は定数項である。なお、Ｎ＝18は累積寄与率＞９５%を満たす主成分数である。
ブラッシング力が大きすぎると摩耗し歯のエナメル質や象牙質が侵食され、結果として欠損を生じる要因となる。そのため、歯面や歯肉へ為害作用を及ぼさないようにブラッシング力はできるだけ小さくしたい。また、プラーク除去率は高いほうがいい。
そこで目的関数に用いる式を式３とし、これを最大にするＰＣＳを求めること、すなわち、単位ブラッシング力あたりのプラーク除去率を最大化することを目的とする。また、比較対象として式３が最小となる条件も検討する。
【００２２】
【数３】
2019030578.tif 000005
【００２３】
次に、目的関数の制約条件について述べる。ブラッシング力Ｆに関しては、実測値のＦｚｍｅａｎの最大値、最小値を探索範囲とし、プラーク除去率Ｅに関しては実測値のＥｓの最大値、最小値の範囲とした。また、求めるＰＣＳにもＰＣＶごとの最大値、最小値を制約条件として適用した。
【００２４】
次に可視化ステップ１１について説明する。
最適化処理ステップ９で導出された主成分得点から、最適ブラッシング動作を再構築することで、最適ブラッシング動作を表示できる。すなわち、[最適ブラッシング動作] ＝[固有ベクトル行列]-1[導出された主成分得点]の関係となる。
図７は、最適化したブラッシング動作（実線）と被験者のブラッシング動作（破線）との比較を可視化した図である。図７では、１ストロークを１００％として、歯ブラシ３の位置と姿勢を示している。従って、図７は、１ストロークを６段階に分けて、それぞれの段階での歯ブラシ３の位置と姿勢について、最適化したブラッシング動作と被験者のブラッシング動作とを示している。また，最適化シミュレーションでは，最適化したブラッシング動作を実施する最適な１ストロークの時間も同時に得られるため，図７の描写をストロークの段階ごとではなく，一定時間ごとに描写することで，ブラッシングのスピードの可視化とそれによる比較も可能である。
Ｆｚｍｅａｎと除去率Ｅｓで最適化をした場合、元の動作よりもブラッシング動作が＋X方向に移動していることから、歯ブラシ３の中央部で磨いていて、歯ブラシ３が歯頚部側に傾いているのが分かる。また、歯ブラシ３が歯頚部側に移動していることが確認された。これは、歯頚部の除去率で最適化を行ったために歯頚部をより磨こうと移動したと考えられる。
【００２５】
図８では、シミュレーションで上記式２によるＦｚｍｅａｎとＥｓを用いて最適化した動作と、最悪化した動作を再現して可視化した。
このように、主成分分析ステップ６と、重回帰分析ステップ８と、最適化処理ステップ９とによってブラッシングの最適化シミュレーションを行える。
このように、本実施例では、１ストロークを複数の段階に分け、各段階での歯ブラシ３の位置と姿勢とを表示して、歯ブラシ３の位置と姿勢の変化の様子を任意視点で可視化することで、最適化前後の歯ブラシ３の動かし方の違いを視覚的に確認することができる。更に最適化した動作、最悪化した動作、及び実際の動作を比較して描写することもできる。
また、時間経過に伴う歯ブラシ３の動きを任意視点における動画像として描写することで、実時間におけるブラッシングの様子やスピードを可視化することができる。更に異なる目的関数を用いて最適化した複数のブラッシング動作の比較も同様である。
なお、可視化にあたっては、理想的な歯ブラシ３の動きと力ベクトルに加えて、実際の動きから予想される力ベクトルの推定値を可視化するか、又は力センサを内蔵した歯ブラシ３を用いて実際の力ベクトルを可視化することもでき、歯ブラシ３の動きと力を視覚的にフィードバックしながら個々の患者に適したブラッシング指導を行え、また個々の患者が各自で最適化したブラッシングを行うことができる。
【００２６】
本発明によれば、従来には出来なかった個々人に最適な歯磨き方法、すなわち歯面や歯肉に為害作用を及ぼさない範囲でプラーク除去を確実に行える最適なブラッシング力を使用者に提示でき、適切な口腔ケアを可能にするシステムを提供できる。
また本発明によれば、事前に作られた多種な磨き方のデータベースと、力覚センサ４付きの歯ブラシ３にてブラッシングしたデータを基に、その人に合った磨き方を最適化計算により導出し提示できる。
また、本発明によれば、歯ブラシ毎に事前にデータベースを作成しておけば、各歯ブラシの特徴やその最適な使い方、個人に適した歯ブラシの選定、実際に使用している歯ブラシに最適なブラッシング方法などを提示することができる。
また、ブラッシング力やプラーク等除去率を有限要素解析などでコンピュータシミュレーションした結果データを実測値と置き換えることもできる。
なお、本実施例ではプラーク除去率を評価パラメータとしたが、歯石の除去率、歯周病菌の除去率、歯茎のマッサージ効果などを評価パラメータとしてもよく、また歯ブラシを動かす距離やスピード、振動数（周期）などの動きの要素を負荷パラメータに加えてもよく、最適化計算に用いる負荷パラメータ、評価パラメータ、及び目的関数は自由に設定することができる。また、本実施例では、歯ブラシヘッド部に作用するブラッシング力を用いたが、歯ブラシのネック部やグリップ部を測定対象とすることもできる。
【産業上の利用可能性】
【００２７】
本発明は、個人に対してだけではなく、歯科医療従事者や口腔ケアに携わる介助者などに対しても、ブラッシングスキルの評価や改善に適している。
【符号の説明】
【００２８】
１  上顎歯列模型
２  模型歯
３  歯ブラシ
３ａ  歯ブラシヘッド部
４  力覚センサ
ステップ１  ブラッシングステップ
ステップ２  ブラッシング計測ステップ
ステップ３  負荷パラメータ計測ステップ
ステップ４  評価パラメータ計測ステップ
ステップ５  正規化処理ステップ
ステップ６  主成分分析ステップ
ステップ７  主成分抽出ステップ
ステップ８  重回帰分析ステップ
ステップ９  最適化処理ステップ
ステップ１０  ブラッシング動作生成ステップ
ステップ１１  可視化ステップ
【図１】
2019030578.tif 000006
【図２】
2019030578.tif 000007
【図３】
2019030578.tif 000008
【図４】
2019030578.tif 000009
【図５】
2019030578.tif 000010
【図６】
2019030578.tif 000011
【図７】
2019030578.tif 000012
【図８】
2019030578.tif 000013