Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7754432B2 - 嚥下障害の判定装置および判定方法 - Google Patents
[go: Go Back, main page]

JP7754432B2 - 嚥下障害の判定装置および判定方法 - Google Patents

嚥下障害の判定装置および判定方法

Info

Publication number
JP7754432B2
JP7754432B2 JP2023551872A JP2023551872A JP7754432B2 JP 7754432 B2 JP7754432 B2 JP 7754432B2 JP 2023551872 A JP2023551872 A JP 2023551872A JP 2023551872 A JP2023551872 A JP 2023551872A JP 7754432 B2 JP7754432 B2 JP 7754432B2
Authority
JP
Japan
Prior art keywords
dysphagia
analysis
phoneme
subject
degree
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2023551872A
Other languages
English (en)
Other versions
JPWO2023054632A5 (ja
JPWO2023054632A1 (ja
Inventor
康宏 大宮
真秀 岩原
萌子 谷内
恵美子 増田
茂樹 河畑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
PST CORPORATION,INC.
SOMPO CARE INC.
Original Assignee
PST CORPORATION,INC.
SOMPO CARE INC.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by PST CORPORATION,INC., SOMPO CARE INC. filed Critical PST CORPORATION,INC.
Publication of JPWO2023054632A1 publication Critical patent/JPWO2023054632A1/ja
Publication of JPWO2023054632A5 publication Critical patent/JPWO2023054632A5/ja
Application granted granted Critical
Publication of JP7754432B2 publication Critical patent/JP7754432B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/09Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being zero crossing rates
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/15Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being formant information
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/18Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being spectral information of each sub-band
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/24Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being the cepstrum
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/48Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
    • G10L25/51Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for comparison or discrimination
    • G10L25/66Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for comparison or discrimination for extracting parameters related to health condition

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Computational Linguistics (AREA)
  • Signal Processing (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Surgery (AREA)
  • Dentistry (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physiology (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Epidemiology (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Description

本開示は嚥下障害の判定装置および判定方法(医療行為を除く)に関する。
飲み込む力が低下した状態である嚥下障害により、食事がうまく取れず、のどに詰まる窒息や、飲み込んだ水分や食べ物が気管に入ってしまう誤嚥のリスクが高まることが知られている。このように嚥下障害を放置しておくと、低栄養状態が続きフレイルや要介護状態となる原因となったり、窒息や誤嚥性肺炎など生命の危機を招いたりすることもある。このため、嚥下障害は、早めに判定し、適切な対策を施すのが望ましい。
従来、被験者の喉の動きを検出することにより、嚥下障害を判定する装置が提案されている。特許文献1は、被験者の頸部に二軸加速度計を位置決めして嚥下障害の有無を判定する装置を提供する。
特許第5977255号公報
しかしながら、特許文献1は加速度計を被験者の頸部に両面テープで貼付する等、被験者に身体的負担を強いるものであった。また、センサ装置の取り付け位置がずれないようにするなど、被験者に対する検査を実施する検査担当者にも負担が多い。このため、被験者や検査担当者の負担を軽減するなど、より簡易に嚥下障害の判定を可能にする装置を研究・開発する余地が残っていた。
本開示は、このような状況に鑑みてなされたものであり、より簡易に嚥下障害の判定が可能な判定装置を提供することを目的とする。
本発明者らは鋭意研究を行った結果、被験者の発話による音声データと、その被験者の嚥下障害の進行の度合いとの間には相関があることを発見した。また、嚥下障害の進行の度合いに基づき分類したグループ間において、所定の音響特徴量に有意な差があることを発見した。これらの発見に基づき、発明者らは音声解析に基づく嚥下障害の判定装置を作成することに成功した。すなわち、本開示の判定装置の詳細は以下の通りである。
[1]音声解析により嚥下障害を判定する判定装置であって、前記判定装置は:
被験者が発話した音声データを入力する入力手段;
前記入力手段で入力した音声データを解析する解析手段;
前記解析手段による解析結果に基づいて被験者の嚥下障害を判定する判定手段、
を備える判定装置。
[2]前記解析手段は、前記音声データを以下の式F(a)に示す音響特徴量を使用して音声解析することを特徴とする、[1]に記載の判定装置:

(式中、gは嚥下障害の有無や進行の度合いを判定する線形または非線形モデルであり、xは前記音声データとして入力するフレーズに固有の係数であり、f(n)は音響パラメータであって、前記音響パラメータは、フォルマント周波数、メル周波数ケプストラム、周波数スペクトラム、音量のエンベロープ、波形の変動情報、ゼロ交差率、ハースト指数、および閉鎖開放から声帯振動が始まるまでの時間、から成る群から1つまたは複数が選択される)。
[3]前記解析手段は、前記音声データをフォルマント周波数またはメル周波数ケプストラムに基づき作成された音響特徴量を使用して音声解析することを特徴とする、[1]または[2]に記載の判定装置。
[4]前記判定装置は、前記解析結果を入力した場合に嚥下障害の進行の度合いを出力するように機械学習処理が施されている、[1]~[3]のいずれか1つに記載の判定装置。
[5]前記判定装置は、被験者の前記音声データの入力を同一のフレーズを使用して少なくとも2回行い、前記同一のフレーズにより解析された前記解析結果の差分または平均値に基づいて、嚥下障害の進行の度合いを判定する、[1]~[4]のいずれか1つに記載の判定装置。
[6]前記入力手段は、被験者が前記嚥下障害の進行の度合いに関連するフレーズを発話した音声データを入力する、[1]~[5]のいずれか1つに記載の判定装置。
[7]前記入力手段は、被験者が

の少なくとも1つの音を含む発話をした音声データを入力する、[1]~[6]のいずれ
か1つに記載の判定装置。
[8]音声解析により嚥下障害を判定する方法であって、前記方法は:
被験者が発話した音声データを入力する入力工程;
前記入力工程で入力した音声データを解析する解析工程;
前記解析工程による解析結果に基づいて被験者の嚥下障害を判定する判定工程、
を備える、方法。
[9][1]~[7]のいずれか1つに記載の判定装置の各手段としてコンピュータを機能させるためのプログラム。
本開示によれば、音声解析により嚥下障害を判定する判定装置を提供することができる。また、音声解析により嚥下障害を判定する判定方法を提供することができる。さらに、判定装置は、被験者の嚥下障害の有無を判定するだけではなく、嚥下障害の進行の度合いを判定することができる。
図1は判定装置の構成を説明する図である。 図2は判定装置により実行される処理の一例を示すフローチャートである。 図3は入出力部12に入力する音声データの一例を示す図である。 図4は本願発明の一実施例の結果を示すROC曲線である。 図5は本願発明の一実施例の結果を示すROC曲線である。 図6は本願発明の一実施例の結果を示す箱ヒゲ図である。 図7は、同じ発話内容で複数回の発話をした場合のフォルマント解析の値のばらつきを示す図である。 図8はフォルマント解析の値と嚥下障害の進行の度合いとの関係を示す図である。 図9Aは13次メル周波数ケプストラム解析の値と嚥下障害の進行の度合いとの関係を示す図である。 図9Bは13次メル周波数ケプストラム解析の値と嚥下障害の進行の度合いとの関係を示す図である。 図10は音韻クラスと活性化される音素との関係を表す図である。 図11は嚥下障害がない健常者が発した「ぱ」と中等度以上の嚥下障害がある人が発した「ぱ」とに対して音素解析を実行した場合の結果である。 図12は音響特徴量による嚥下障害の重症度による分類性能を検証した箱ヒゲ図である。 図13は音響特徴量による嚥下障害の重症度による分類性能を検証した箱ヒゲ図である。 図14は音響特徴量による嚥下障害の重症度による分類性能を検証した箱ヒゲ図である。 図15は音響特徴量による嚥下障害の重症度による分類性能を検証した箱ヒゲ図である。 図16は音響特徴量による嚥下障害の重症度による分類性能を検証した箱ヒゲ図である。 図17は音響特徴量による嚥下障害の重症度による分類性能を検証した箱ヒゲ図である。 図18は音声データを示す図である。 図19は音声の強度解析結果と嚥下障害の評価結果との相関関係を表す実験結果である。 図20は各音響特徴量の箱ヒゲ図の結果である。 図21は各音響特徴量の箱ヒゲ図の結果である。 図22は音素解析および強度解析の手法を含めて作成された音響特徴量セットを用いた機械学習モデルにより、軽度以上の嚥下障害を判定した際のROC曲線の結果である。 図23は音素解析および強度解析の手法を含めて作成された音響特徴量セットを用いた機械学習モデルにより、中等度以上の嚥下障害を判定した際のROC曲線の結果である。
以下、本開示を実施するための形態について、図面を参照して詳細に説明するが、以下に記載する構成要件の説明は、本開示の一実施態様としての一例であり、これらの内容に限定されるものではない。
<第1実施形態>
本開示の嚥下障害の判定装置は、主な構成要素として解析手段と判定手段を備える。解析手段は嚥下障害の進行の度合いを解析可能な音響特徴量(以下、「F(a)」と称することがある。)を使用して音響解析を行う。また、判定手段は、解析手段により取得した解析結果を入力した場合に嚥下障害の進行の度合いを出力するような機械学習処理が施されていてもよい。
判定装置(以下、「情報処理装置」と称することがある。)の構成について図1を用いて説明する。情報処理装置10は、全体の動作を制御する制御部11、各種の入出力を行う入出力部12、各種データやプログラム等を記憶する記憶部13、外部との通信を行う通信部14、および各ブロック同士が相互通信可能なように接続する内部バス15、を備える。
情報処理装置10は、例えばコンピュータであり、スマートフォン、PDA、タブレット、又はノートパソコンなどの、被験者が携帯可能な装置であってもよいし、被験者が携帯せずに設置位置に固定されたコンピュータでもよい。PDAは、Personal Digital Assistantの略称である。
制御部11は、例えばCPU、MCU又はMPUと呼ばれる装置であり、例えば記憶部13に記憶されたプログラムが実行される。CPUはCentral Processing Unitの略称である。MCUはMicro Controller Unitの略称である。MPUはMicro Processor Unitの略称である。
入出力部12は、情報処理装置10を操作する被験者に対する入出力を行う装置である。入出力部12は、ディスプレイ、キーボード、マウス、ボタン、タッチパネル、プリンタ、マイク、及びスピーカなどによる、情報や信号の入出力を行う。本実施例では、入出力部12は、少なくともマイクの役割を果たし、このマイクによって音声データを入力する。また、本実施例では、入出力部12は、少なくともディスプレイの役割を果たし、このディスプレイに、後述する嚥下障害の判定結果を表示する。
記憶部13は、例えば、ROM、RAM、HDD、又はフラッシュメモリといった装置であり、制御部11で実行するプログラムや各種データを記憶する。ROMは、Read Only Memoryの略称である。RAMは、Random Access Memoryの略称である。HDDは、Hard Disk Driveの略称である。
通信部16は、外部との通信を行う。通信部16による通信は、有線通信でもよいし、無線通信でもよい。通信部16による通信は、如何なる通信方式でもよい。制御部11は、通信部16によって、音声データなどの各種データの送受信を行うことが出来る。制御部11は、後述する嚥下障害の判定結果を、通信部16によって外部機器に送信するようにしてもよい。
次に図2を用いて判定装置(情報処理装置)10により実行される処理の一例を示す。まず、ステップS201において、制御部11は入出力部12によって被験者の音声データを入力する。次に、ステップS202において、算出部(または解析部)で音声データから音響特徴量を算出する。次に、ステップS203において、推定部(または判定部)で嚥下障害の有無や進行の度合いを推定(または判定)する。次に、推定結果(または判定結果)を入出力部12に出力して、フローを終了する。
なお、ステップS201の入出力部12はマイクを使用してもよい。被験者はマイクに向けて発話を行い、音声データを入力する。音声データは予め録音しておいた音声データを使用してもよい。
<音声入力に際し選択されるフレーズ>
音声入力に際し選択されるフレーズは、被験者の嚥下障害の進行の度合いを音声解析するのに適切なフレーズである。嚥下障害の進行が進むと、舌の動く程度、舌の前後位置、顎の開き具合、歯の噛み合わせ状態、歯の数、唾液の分泌量、筋肉の衰え等に起因し、咽喉頭や口腔内での音の共鳴の具合なども影響をうける。嚥下障害の進行の度合いを解析するのに適切なフレーズとは、上記の嚥下障害の進行に伴い変化する音の共鳴の具合などを発見し易いフレーズである。
例えば、音声入力に際し選択されるフレーズとして

の音を含むフレーズが好ましい。より詳細な説明は図3の説明と併せて後述する。
図3は、入出力部12に入力する音声データの一例を示す。図3は

の発音を含むフレーズをまとめた一例であるため、本開示の音声入力をこれらのフレーズに限定するものではない。被験者は、図3に示すPhrase01(Ph01)~Phrase10(Ph10)のうち、少なくとも1つを音声入力すべきフレーズとして選択できる。もちろん、Phrase01(Ph01)~Phrase10(Ph10)のうち幾つかを組合わせて音声入力してもよい。
図3において、Ph01は「ぱ」と発話する音声データである。Ph02は「ま」と発話する音声データである。Ph03は「た」と発話する音声データである。Ph04は「ら」と発話する音声データである。Ph05は「か」と発話する音声データである。Ph06は「ご」と発話する音声データである。Ph07は「パンダのたからもの」と発話する音声データである。Ph08は「たまご」と発話する音声データである。Ph09は「バナナバナナバナナバナナバナナ」と出来るだけ速く5回以上繰り返して発話する音声データである。Ph10は「きものきものきものきものきもの」と出来るだけ速く5回以上繰り返して発話する音声データである。フレーズPh01~Ph10と、嚥下障害との関連性についてさらに説明する。
<<Ph01、Ph02>>
Ph01、Ph02は、「ぱ(pa)」、「ま(ma)」の発音であり、唇を閉じる動きが必要で、嚥下機能としては咀嚼の時に口腔内の食べ物をこぼさないようにし、口腔内圧を高めることによって飲み込む際の食事の移送に関係する。
その他の発音の例としては、

等がある。
<<Ph03、Ph04>>
Ph03、Ph04は、「た(ta)」、「ら(ra)」の発音である。「た(ta)」は舌の先を使う動きであり、嚥下機能としては、咀嚼の機能、送り込み動作(口の中の水分や食物を喉の奥に移動する動き)に関係する。「ら(ra)」は舌の先を比較的にスムーズに動かす必要があり、舌の動きのスムーズさをみられる音である。「た」と同じように舌の先を使い、咀嚼の機能、送り込み動作(口の中の水分や食物を喉の奥に移動する動き)に関係する。
その他の発音の例としては、

等がある。
<<Ph05、Ph06>>
Ph05、Ph06は、「か(ka)」、「ご(go)」の発音である。両方とも舌の奥を使う動きであり、嚥下機能としては、送り込み動作や咽頭内圧を高める動作を行い、食物の移送に関係する。
その他の音の例としては、

等がある。
<<Ph07、Ph08>>
Ph07(パンダのたからもの)、Ph08(たまご)は唇および舌による音の組合せによる評価を行っている。
<<Ph09、Ph10>>
Ph09(バナナバナナバナナバナナバナナ)、Ph10(きものきものきものきものきもの)は、鼻腔共鳴の評価、唇および舌による組合せの評価、ディアドコキネシスの評価、リズムの評価を行っている。
以上のようなフレーズを音声入力して音声データを取得し、ステップS202の音声解析を行う。音声解析の際に算出部(解析部)において音響特徴量の算出を行う。音響特徴量について以下に詳しく述べる。
<音響特徴量>
音響特徴量とは、音声解析すべき音声データの特徴を、定量的に表すための数値パラメータである。本開示では、ステップS203において、判定部で嚥下障害の有無や進行の度合いを判定するが、音響特徴量に基づく音声解析の結果を入力として判定処理を出力するように機械学習処理が施されるため、精度向上のためにも音響特徴量の選択は重要である。
本開示において、音響特徴量F(a)は、以下の式で示すことができる。
式中、gは嚥下障害の有無や進行の度合いを判定する線形または非線形モデルであり、xは音声データとして入力するフレーズに固有の係数であり、f(n)は音響パラメータであって、フォルマント周波数、メル周波数ケプストラム、周波数スペクトラム、音声のエンベロープ、波形の変動情報、ゼロ交差率、ハースト指数、および閉鎖開放から声帯振動が始まるまでの時間、から成る群から1つまたは複数が選択される。また、同じフレーズの音声データが2つ以上ある場合には、その平均値または差分、3回以上の場合にはバラツキ(分散や標準偏差)や中央値などを含むことができる。また、音響特徴量は、数値に大きな開きがあるため、それぞれを正規化しても良い。さらに、3群以上の嚥下障害の有無や進行の度合いを判定する場合は、特徴量を2つ以上に分割しても良い。
音響パラメータの種類には、以下のようなものがある。
(1) 音声のエンベロープ(アタックタイム、ディケイタイム、サステインレベル、リリースタイム)
(2) 波形の変動情報(Shimmer、Jitter、倍音対雑音比 HNR:Harmonics to Noise Ratio、信号雑音比 SNR:Signals to Noise Ratio)
(3) ゼロ点交差率
(4) ハースト指数
(5) 閉鎖開放から声帯振動が始まるまでの時間(VOT:Voice Onset Time)
(6) メル周波数ケプストラムのある係数に関する発話内分布の統計値(第1四分位点、中央値、第3四分位点、95パーセント点、98パーセント点、算術平均、幾何平均、第3四分位点と中央値の差など)
(7) 周波数スペクトラムの変化の速さにおける発話内分布の統計値(第1四分位点、中央値、第3四分位点、95パーセント点、98パーセント点、算術平均、幾何平均、第3四分位点と中央値の差など)
(8) メル周波数ケプストラムのある係数の時間変化に関する発話内分布の統計値(第1四分位点、中央値、第3四分位点、95パーセント点、98パーセント点、算術平均、幾何平均、第3四分位点と中央値の差など)
(9) メル周波数ケプストラムのある係数の時間変化の時間変化に関する発話内分布の統計値(第1四分位点、中央値、第3四分位点、95パーセント点、98パーセント点、算術平均、幾何平均、第3四分位点と中央値の差など)
(10)周波数スペクトラム90パーセントロールオフの発話内時間変化における2次回帰近似に対する二乗誤差
(11)周波数スペクトラム重心の発話内時間変化における2次回帰近似に対する算術誤差その他、ピッチレート、有声音である確率、任意範囲の周波数パワー、音階、話速(一定時間におけるモーラ数)、ポーズ・間、音量など
(12)任意のフォルマント周波数(第1フォルマント、第2フォルマント、第3フォルマント、第4フォルマント、第5フォルマント、第6フォルマント、・・・)に関する発話内分布の統計値(第1四分位点、中央値、第3四分位点、95パーセント点、98パーセント点、算術平均、幾何平均、第3四分位点と中央値の差など)
(13)任意のフォルマント周波数(第1フォルマント、第2フォルマント、第3フォルマント、第4フォルマント、第5フォルマント、第6フォルマント、・・・)の時間変化に関する発話内分布の統計値(第1四分位点、中央値、第3四分位点、95パーセント点、98パーセント点、算術平均、幾何平均、第3四分位点と中央値の差など)
(14)任意のフォルマント周波数(第1フォルマント、第2フォルマント、第3フォルマント、第4フォルマント、第5フォルマント、第6フォルマント、・・・)の時間変化の時間変化に関する発話内分布の統計値(第1四分位点、中央値、第3四分位点、95パーセント点、98パーセント点、算術平均、幾何平均、第3四分位点と中央値の差など)
上述してきた音響特徴量の中から、嚥下障害の進行の度合いと相関性を有する音響特徴量に基づいて、あるいは嚥下障害の程度で分類を行った際の各群間における音響特徴量の差に基づいて、ステップS203では判定の処理を実行する。上記の音響特徴量を用いて判定を行った例を図4~6を用いて説明する。
被験者は65歳以上の高齢者であり、被験者の嚥下機能は、嚥下機能テストの結果や誤嚥の度合いから言語聴覚士が確認をし、被験者をそれぞれ嚥下障害がない健常者と、軽度の嚥下障害がある人と、中等度以上の嚥下障害がある人に分類した。
(実施例1)
図4は、被験者が図3に示す10種類のフレーズを2回ずつ読み上げた音声データを入力とし、上述した音響パラメータ(1)~(14)のうち7つを用いて特定のプログラムを作成した場合、嚥下障害の有無の分類性能を検証したROC曲線である。横軸は「1-特異度」を示し、縦軸は感度を示す。AUCは0.941であり、十分な分類性能を有していることが確認できた。
(実施例2)
図5は、被験者が図3に示す10種類のフレーズを2回ずつ読み上げた音声データを入力とし、上述した音響パラメータ(1)~(14)のうち1つの音響パラメータを用いて音響特徴量を算出し、その値の平均値を用いて作成した特定のプログラムの嚥下障害の進行の度合い(中等度以上の嚥下障害であるか否か)に関する分類性能を検証したROC曲線である。横軸は「1-特異度」を示し、縦軸は感度を示す。AUCは0.981であり、十分な分類性能を有していることが確認できた。また、図4、5の結果から判断して、プログラムを組合わせて処理することにより、まず図4に係るプログラムを用いて嚥下障害の有無を判定し、次いで嚥下障害が有りと判定された被験者群をさらに図5に係るプログラムを用いて嚥下障害の進行の度合い(軽度または中程度以上)を判定することもできる。
(実施例3)
図6は、被験者が図3に示す10種類のフレーズを2回ずつ読み上げた音声データを入力とし、上述した音響パラメータ(1)~(14)のうち1つの音響パラメータを用いて音響特徴量を算出し、その値の平均値を用いて作成した特定のプログラムの嚥下障害の進行の度合いに関する分類性能を検証した箱ヒゲ図である。図6では、健常者、軽度嚥下障害、および中等度以上の嚥下障害を区別する。横軸は健常、軽度嚥下障害、中等度以上の嚥下障害を示し、縦軸は3群判定のスコアの分布を示す。図6からも、上記プログラムが十分な分類性能を有していることが確認できた。
上記図4~6およびそれら図面に関する説明から、特定の音響特徴量を用いて判定を行う判定装置は十分に実用に耐えうることが確認できた。次に、音声データと音響パラメータとの対応関係に関して、幾つか具体例を用いて説明する。なお繰り返しとなるが以降の具体例は、本願発明を実施する上での一例であり、本願発明がこれらの実施例に限定されるものではない。
(実施例4)
図7は、音声データをフォルマント解析した結果について、同じ発話内容で複数回の発話をした場合のばらつきと、被験者の嚥下障害の進行の度合いとの関係を示す図である。図7では、発話内容としてPh07「パンダのたからもの」の音声データを用いている。横軸は、健常、軽度嚥下障害、中等度以上の嚥下障害を有する被験者の発話に関する時間軸であり、縦軸は第1フォルマントf1に係る値である。それぞれ1回目の発話と2回目の発話に関し、複数発話ごとにまとめ、その順番を時間軸に沿ってプロットしている。被験者はそれぞれが2回ずつ発話しており、最初の8人は健常者(グラフ左)であり、次の17人は軽度嚥下障害を有し(グラフ中央)、次の18人は中等度以上の嚥下障害を有している(グラフ右)。
図7から、健常者の場合、1回目の発話と2回目の発話とでf1の値の差分が小さいことが分かる。これに対して、軽度の嚥下障害がある人及び中等度以上の嚥下障害がある人の場合、1回目の発話と2回目の発話とでf1の値の差分が大きい。このことから、ステップS201で入力した音声データがPh07の場合、ステップS203では、f1により、今回の被験者が健常者であるか、それとも軽度の嚥下障害がある人及び中等度以上の嚥下障害がある人であるかを判定することが出来る。
たとえば、ステップS203では、健常者のf3又はf5をあらかじめ記憶しておき、被験者のf3又はf5が、健常者のf3又はf5に対して閾値以上の乖離がある場合に、その被験者には嚥下障害があると判定することが出来る。図7では、Ph07のf1についてのみ示したが、他の音声データ及び他のフォルマント解析結果(例えばf1以外のf2~f5のいずれか)並びにその他の周波数解析についても有効性を確認することが出来る。また、図7では、2回の発話でのf1の値の差分を対象にしたが、他の発話回数を採用することも出来る。例えば3回以上の発話でのフォルマント解析結果の最大値と最小値との差分により、嚥下障害の程度(嚥下障害の有無、及び重症度)を判定するようにしてもよい。
(実施例5)
図8は、図3に示した各音声データと、フォルマント周波数に基づく音響特徴量の値とを比較し、かつ嚥下障害の進行の度合いを照らし合わせた表である。横方向の項目は、「健常vs軽度」、「健常vs中等度以上」、「軽度vs中等度以上」の各項目を示す。縦方向の項目は、発話内容を示す。表中、「***」はP値<0.01、「**」はP値<0.03333、「*」はP値<0.05、「ns」は有意差がなかったことをそれぞれ示している。なお、本願では、P値<0.1であれば有意差ありとみなしている。また、f1は第1フォルマント、f2は第2フォルマント、f3は第3フォルマント、f4は第4フォルマント、f5は第5フォルマントである。有意差の評価には、2群比較での有意差の検定ではt検定(対応なし、片側)を、3群以上の多重比較検定ではボンフェローニの多重比較検定を用い、有意水準には10%を設定した。なお、群比較による有意差の評価には、今回用いたt検定を含むパラメトリック検定、ノンパラメトリック検定、比率による検定、分散比による検定など、いずれを用いても良い。
図8の結果から、Ph01は、「健常vs軽度」の項目においてf3、f5で有意差を示した。また「健常vs中等度以上」の項目においてf1、f5で有意差を示した。また、「軽度vs中等度以上」の項目において、f1、f2で有意差を示した。
また、Ph02は、「健常vs軽度」の項目においてf4、f5で有意差を示した。また「健常vs中等度以上」の項目においてf3~f5で有意差を示した。また、「軽度vs中等度以上」の項目において、f3、f4で有意差を示した。
また、Ph03は、「健常vs軽度」の項目においてf3で有意差を示した。また「健常vs中等度以上」の項目においてf1、f2で有意差を示した。また、「軽度vs中等度以上」の項目において、f1~f4で有意差を示した。
また、Ph04は、「健常vs中等度以上」の項目においてf2~f4で有意差を示した。また、「軽度vs中等度以上」の項目において、f1~f4で有意差を示した。Ph04は「健常vs軽度」の項目においていずれのフォルマントも有意差を示さなかったが、他の2群を比較することにより「健常vs軽度」の判定をすることは可能である。
また、Ph05は、「健常vs軽度」の項目においてf2、f5で有意差を示した。また「健常vs中等度以上」の項目においてf1で有意差を示した。また、「軽度vs中等度以上」の項目において、f1、f2で有意差を示した。
また、Ph06は、「健常vs軽度」の項目においてf3、f5で有意差を示した。また「健常vs中等度以上」の項目においてf5で有意差を示した。Ph06は「軽度vs中等度以上」の項目においていずれのフォルマントも有意差を示さなかったが、他の2群を比較することにより「軽度vs中等度以上」の判定をすることは可能である。
また、Ph07は、「健常vs軽度」の項目においてf3、f5で有意差を示した。また「健常vs中等度以上」の項目においてf1で有意差を示した。また、「軽度vs中等度以上」の項目において、f1、f3で有意差を示した。
また、Ph08は、「健常vs軽度」の項目においてf1、f3、f5で有意差を示した。また「健常vs中等度以上」の項目においてf1で有意差を示した。また、「軽度vs中等度以上」の項目において、f1、f3、f4で有意差を示した。
また、Ph09は、「健常vs軽度」の項目においてf2~f4で有意差を示した。また「健常vs中等度以上」の項目においてf1、f4、f5で有意差を示した。また、「軽度vs中等度以上」の項目において、f3、f4で有意差を示した。
また、Ph10は、「健常vs軽度」の項目においてf3、f5で有意差を示した。また「健常vs中等度以上」の項目においてf1、f3、f5で有意差を示した。また、「軽度vs中等度以上」の項目において、f1、f3で有意差を示した。
以上の結果から、Ph01~Ph10は、全てのフレーズにおいて、健常と軽度嚥下障害、健常と中等度以上の嚥下障害、軽度嚥下障害と中等度以上の嚥下障害とを区別できることが確認できた。従って、被験者が

の少なくとも1つの音を含む発話をした音声データを入力すれば、健常と軽度嚥下障害、健常と中等度以上の嚥下障害、軽度嚥下障害と中等度以上の嚥下障害とを区別できることが確認できた。
フォルマント解析は従来、口唇から舌の範囲を第1、第2フォルマントで解析する傾向にあるため、第1、第2フォルマントまでを重視し、第3フォルマント以降を軽視する傾
向にあった。しかし、本発明者らは、特に嚥下障害の進行の度合いを判定する上で口唇よりも奥の咽頭腔の形状や、その複雑な形状変化の連動は、第3フォルマント以降を確認しなければ分からないことを音声解析の研究によって明らかにした。特に高齢者の嚥下障害においては、咽頭腔の形状変化もあり、f3以降の音声解析を行うことは嚥下障害の進行の度合いを判定する上で有効である。
(実施例6)
さらに、図9Aおよび図9Bは、図3に示した各音声データを2回発話したデータの差分をとり13次メル周波数ケプストラム解析した結果の、嚥下障害の程度(嚥下障害の有無、及び重症度)の判定における有効な特徴量の1つを示した表である。
図9Aの表は図3に示した各音声データを2回発話した音声データの2発話の平均値を音響パラメータとして用いたものである。図9Bの表は図3に示した各音声データを2回発話した音声データの2発話間の差分を音響パラメータとして用いたものである。それぞれ、横方向の項目には、13次メル周波数ケプストラム係数または動的特徴量の平均値、最大値、最小値、範囲値、平均最小値、傾きを示す。縦方向の項目は発話内容を示す。表中、「〇」(マル)は嚥下障害の有無の判定に有効であることを示す。二重のマルは、嚥下障害の進行の度合いの判定に有効であることを示す。
図9Aおよび図9Bから、音声解析として13次メル周波数ケプストラム解析を行った場合も、健常と軽度嚥下障害、健常と中等度以上の嚥下障害、軽度嚥下障害と中等度以上の嚥下障害とを区別して判定できることが確認できた。
<第2実施形態>
次に、第2実施形態について説明する。なお、第2実施形態の判定装置の構成は、第1実施形態と同様の構成となるため、同一符号を付して説明を省略する。
第2実施形態の判定装置は、被験者の嚥下障害を判定する際に、音素解析を利用する点が第1実施形態と異なる。
第2実施形態の判定装置の解析手段は、被験者が発話した音声データに対して音素解析を行うことにより音響特徴量を作成する。そして、第2実施形態の判定装置の判定手段は、解析手段による解析結果である音響特徴量を使用して音声解析を実行することにより、嚥下障害の進行の度合いを判定する。
(実施例7)
図10は、音韻クラスと活性化される音素との関係を表す図である。なお、同一の音素であっても、複数の音韻クラスに属する場合がある。例えば、図10に示されているように、音素/a/は、「vocalic」「back」「open」「voiced」に属し、音素/p/は「consonantal」「stop」「labial」に属する。また、図11は、嚥下障害がない健常者は発した「ぱ」と、中等度以上の嚥下障害がある人が発した「ぱ」とに対して音素解析を実行した場合の結果である。
図11の(A),(C),(E),(G),及び(I)は、健常者が発した「ぱ」に対する音素解析結果であり、図11の(B),(D),(F),(H),及び(J)は、中等度以上の嚥下障害がある人が発した「ぱ」に対する音素解析結果である。なお、図11の各グラフは横軸が時刻(s)を表し、縦軸が各音韻クラスの音素事後確率(Phonological posteriors)を表す。図11に示されている中央部の波形が、被験者から発せられた音声の音声データに相当する。
図11の(A)と(B)とは、音韻クラスが「vocalic」「back」「consonantal」「anterior」である場合の音素解析結果である。また、図11の(C)と(D)とは、音韻クラスが「open」「nasal」「close」「stop」である場合の音素解析結果である。また、図11の(E)と(F)とは、音韻クラスが「continuant」「flap」「lateral」「trill」である場合の音素解析結果である。また、図11の(G)と(H)とは、音韻クラスが「voice」「labial」「strident」「dental」である場合の音素解析結果である。また、図11の(I)と(J)とは、音韻クラスが「velar」「pause」である場合の音素解析結果である。
図11に示されているように、健常者と中等度以上の嚥下障害がある人とでは、音素解析結果である音韻クラスの波形において大きな差異がある。このため、音素解析結果は、嚥下障害の有無を判定する際の音響特徴量として有用であると考えられる。
図12~図17に、音韻クラスの音素解析結果である、各音響特徴量による嚥下障害の有無の分類性能を検証した箱ヒゲ図を示す。なお、図中の音韻クラス(例えば、「consonantal」「close」「dental」「velar」「stop」「anterior」「back」「continuant」「open」「labial」)の後に記載されている「mean」は平均を表し、「median」は中央値を表し、「std」は標準偏差を表す。また、図12~図17内の「healthy」は健常者を表し、「mild」は嚥下障害が軽度である人を表し、「severe」は嚥下障害が中等度以上である人を表す。また、図12~図17内の縦軸は、平均、中央値、及び標準偏差等の各統計量の値を表す。具体的には、図12~図17内の縦軸は、被験者があるフレーズを発した際の、各時刻の音韻クラスの音響特徴量の値の平均、中央値、及び標準偏差等である。
図12は、被験者が「ぱ」を発した際の音響特徴量の統計量である。図12の上段に示されているように、健常者と嚥下障害が軽度以上である人とを分別する際には、その一例として、音韻クラス「close」の標準偏差である「close_std」と、音韻クラス「consonantal」の標準偏差である「consonantal_std」が有用な音響特徴量であることが分かる。また、図12の下段に示されているように、健常者及び嚥下障害が軽度である人と嚥下障害が中等度以上である人とを分別する際には、その一例として、音韻クラス「close」の平均である「close_mean」と、音韻クラス「close」の中央値である「close_median」が有用な音響特徴量であることが分かる。
図13も、被験者が「ぱ」を発した際の音響特徴量の統計量である。図13に示されているように、健常者及び嚥下障害が軽度である人と嚥下障害が中等度以上である人を分別する際には、その一例として、音韻クラス「dental」の平均値である「dental_mean」と、音韻クラス「dental」の中央値である「dental_median」と、音韻クラス「velar」の中央値である「velar_median」とが有用な音響特徴量であることが分かる。
また、図14及び図15は、被験者が「ら」を発した際の音響特徴量の統計量である。図14に示されているように、健常者と、嚥下障害が軽度である人と、嚥下障害が中等度以上である人を分別する際には、その一例として、音韻クラス「close」の平均値である「close_mean」と中央値である「close_median」と、音韻クラス「dental」の標準偏差である「dental_std」と、音韻クラス「stop」の平均値である「stop_mean」が有用な音響特徴量であることが分かる。
また、図15に示されているように、健常者と、嚥下障害が軽度である人と、嚥下障害が中等度以上である人を分別する際には、その一例として、音韻クラス「velar」の平均値である「velar_mean」と中央値である「velar_median」が有用な音響特徴量であることが分かる。
また、図16及び図17は、被験者が「たまご」を発した際の音響特徴量の統計量である。図16及び図17については詳細な説明は省略するものの、図12~図15と同様に、各音韻クラスの音響特徴量は、嚥下障害の度合いを判定するのに有用である。
以上図12~図17に例示されているように、各音韻クラスの統計量は、嚥下障害の度合いを判定するのに有用である。このため、健常と軽度嚥下障害、健常と中等度以上の嚥下障害、軽度嚥下障害と中等度以上の嚥下障害とを区別して判定できることが確認できる。
このため、第2実施形態の判定装置の解析手段は、被験者が発した音声の音声データに対して音素解析を行うことにより、音声データに含まれる音素を特定する。そして、第2実施形態の判定装置の解析手段は、音声データに含まれる音素が属する音韻クラスを特定する。なお、同一の音素であっても、複数の音韻クラスに属する場合がある。例えば、図10に示されているように、音素/a/は、「vocalic」「back」「open」「voiced」に属し、音素/p/は「consonantal」「stop」「labial」に属する。このため、第2実施形態の判定装置の解析手段は、時系列データである音声データの各時刻の音韻クラスの音素事後確率(Phonological posteriors)を計算する。そして、第2実施形態の判定装置の判定手段は、解析手段によって得られた各時刻の音韻クラスの音素事後確率の統計量(例えば、各時刻の音素事後確率の平均、中央値、及び標準偏差等)に基づいて、被験者の嚥下障害の度合いを判定する。このため、第2実施形態の判定装置によれば、被験者が発した音声の音声データを音素解析することにより、被験者の嚥下障害の度合いを精度良く判定することができる。
<第3実施形態>
次に、第3実施形態について説明する。なお、第3実施形態の判定装置の構成は、第1実施形態と同様の構成となるため、同一符号を付して説明を省略する。
第3実施形態の判定装置は、被験者の嚥下障害を判定する際に、音声の強度解析を利用する点が第1実施形態及び第2実施形態と異なる。
第3実施形態の判定装置の解析手段は、被験者が発話した音声データに対して音声の強度解析を行うことにより音響特徴量を作成する。そして、第3実施形態の判定装置の判定手段は、解析手段による解析結果である音響特徴量を使用して音声解析を実行することにより、嚥下障害の進行の度合いを判定する。
(実施例8)
図18は、音声データを示す図である。図18の横軸は時刻を表し、縦軸は音声の強度を表す。なお、図18中の黒丸は、音声データのピーク点を表す。
図19の上段の4行は、被験者が発したフレーズPh09「ばななばなな・・・ばなな」及びPh10「きものきものきもの・・・きもの」についての音声の強度解析結果と嚥下障害の評価結果との相関関係を表す実験結果である。図19の「peak_ave」は音声データの強度のピーク点の平均を表し、「peak_sd」は音声データの強度のピーク点の標準偏差を表し、「peak_span_ave」は音声データの強度のピーク点間の間隔の平均を表し、「peak_span_sd」は音声データの強度のピーク点間の間隔の標準偏差を表している。
図19は、ボンフェローニ法による検定結果のうちの「peak_ave」「peak_sd」「peak_span_ave」「peak_span_sd」を用いた、健常者と嚥下障害が軽度以上である人の分別、および健常者及び嚥下障害が軽度である人と嚥下障害が中等度以上である人の分別の結果である。図19に示されているように、Ph09の「peak_sd」を用いた健常者及び嚥下障害が軽度である人と嚥下障害が中等度以上である人の分別が有効であること、Ph09の「peak_span_sd」を用いた健常者及び嚥下障害が軽度である人と嚥下障害が中等度以上である人の分別が有効であること、Ph10の「peak_ave」を用いた健常者及び嚥下障害が軽度である人と嚥下障害が中等度以上である人の分別が有効であること、Ph10の「peak_sd」を用いた健常者及び嚥下障害が軽度である人と嚥下障害が中等度以上である人の分別が有効であることがわかる。
図20及び図21は、各音響特徴量の箱ヒゲ図である。
図20の左側は、フレーズPh09についての音響特徴量「peak_sd」の箱ヒゲ図であり、嚥下障害の症状との関連が見られ、特に嚥下障害が中等度以上である人の分別に有効であることが確認できた。
図20の右側は、フレーズPh09についての音響特徴量「peak_span_sd」の箱ヒゲ図であり、嚥下障害の症状との関連が見られ、特に嚥下障害が中等度以上である人の分別に有効であることが確認できた。
図21は、フレーズPh10についての音響特徴量「peak_sd」の箱ヒゲ図であり、嚥下障害の症状との関連が見られ、特に嚥下障害が中等度以上である人の分別に有効であることが確認できた。
以上、図19~図21に示されているように、音声データの強度解析結果は、嚥下障害の度合いを判定するのに有用である。このため、健常と軽度嚥下障害、健常と中等度以上の嚥下障害、軽度嚥下障害と中等度以上の嚥下障害とを区別して判定できることが確認できる。
このため、第3実施形態の判定装置の解析手段は、被験者が発した音声の音声データに対して強度解析を行うことにより、強度解析に関する音響特徴量を生成する。例えば、第3実施形態の判定装置の解析手段は、強度解析に関する音響特徴量として、上述したようなデータを生成する。そして、第3実施形態の判定装置の判定手段は、解析手段によって得られた音響特徴量に基づいて、被験者の嚥下障害の度合いを判定する。このため、第3実施形態の判定装置によれば、被験者が発した音声の音声データを強度解析することにより、被験者の嚥下障害の度合いを精度良く判定することができる。なお、上記実施例では、被験者が発するフレーズとしてフレーズPh09「ばななばなな・・・ばなな」及びPh10「きものきものきもの・・・きもの」を例示したが、これに限らず、所定フレーズの繰り返しを含むようなものであれば、どのようなものであってもよい。例えば「ぱたかぱたかぱたか」「ぱぱぱぱ」「たたたた」といったフレーズであってもよい。このようなフレーズは、鼻腔共鳴(鼻にかかる)を発生させるフレーズであり、嚥下障害の度合いを判定するのに有用であると考えられる。また、繰り返しを含むフレーズを採用することにより、繰り返しによって発生するリズムの一定性を判定することができるため、これも、嚥下障害の度合いを判定するのに有用であると考えられる。
なお、第2実施形態の音素解析と第3実施形態の強度解析とを組み合わせて、被験者の嚥下障害の度合いを判定するようしにてもよい。この場合には、判定装置の解析手段は、音声データに対して音素解析を行うことにより作成された音響特徴量と、音声データに対して強度解析を行うことにより作成された音響特徴量とを用いて、被験者の嚥下障害の度合いを判定する。これにより、音素解析及び強度解析の何れか一方のみを用いる場合よりも、被験者の嚥下障害の度合いを精度良く判定することができる。なお、音響特徴量の組み合わせは、音素解析に関する音響特徴量と強度解析に関する音響特徴量との組み合わせに限定されず、フォルマント周波数、メル周波数ケプストラム、周波数スペクトラム、音声のエンベロープ、波形の変動情報、ゼロ交差率、ハースト指数、及び閉鎖開放から声帯振動が始まるまでの時間等を含めたいかなる組み合わせであってもよい。
(実施例9)
図22及び図23は、音声データに対して音素解析を行うことにより作成された音響特徴量と、音声データに対して強度解析を行うことにより作成された音響特徴量とを用いて、被験者の嚥下障害の度合いを判定した場合のROC曲線図である。図22及び図23は、音響特徴量セットを用いた機械学習モデルにより、軽度以上の嚥下障害(図22)又は中等度以上の嚥下障害(図23)を判定した際のROC曲線の結果である。より詳細には、図22及び図23は、フォルマント周波数に基づき作成された音響特徴量、メル周波数ケプストラムに基づき作成された音響特徴量、音素解析を行うことにより作成された音響特徴量、及び強度解析を行うことにより作成された音響特徴量の4つの音響特徴量を用いて、被験者の嚥下障害の度合いを判定した場合のROC曲線図である。図22は、軽度以上の嚥下障害を判定した場合のROC曲線図である。また、図23は、中等度以上の嚥下障害を判定した場合のROC曲線図である。図22に示される軽度以上の嚥下障害を判定した場合のAUCは0.9954(P値<0.001)であり、十分な分類性能を有していることが確認できた。また、図23に示される中等度度以上の嚥下障害を判定した場合のAUCは1.000(P値<0.001)であり、十分な分類性能を有していることが確認できた。このため、音声データに対して音素解析を行うことにより作成された音響特徴量と、音声データに対して強度解析を行うことにより作成された音響特徴量とを用いることにより、被験者の嚥下障害の度合いを精度良く判定することができることが確認された。
以上、本開示の好ましい実施例について説明したが、本開示は、上述の実施例に限定されない。本開示の目的は、上述の実施例の機能を実現するプログラムコード(コンピュータプログラム)を格納した記憶媒体をシステムあるいは装置に供給し、供給されたシステムあるいは装置のコンピュータが記憶媒体に格納されたプログラムコードを読み出し実行することによっても達成される。この場合、記憶媒体から読み出されたプログラムコード自体が上述した実施例の機能を実現することになり、そのプログラムコードを記憶した記憶媒体は本開示を構成することになる。また、上述した実施形態では、コンピュータがプログラムを実行することにより、各処理部として機能するものとしたが、処理の一部または全部を専用の電子回路(ハードウェア)で構成するようにしても構わない。本開示は、説明された特定の実施例に限定されるものではなく、各実施例の各構成同士の置換を含み、特許請求の範囲に記載された本開示の趣旨の範囲内において、種々の変形・変更が可能である。
2021年9月29日に出願された日本国特許出願2021-159606号の開示は、その全体が参照により本明細書に取り込まれる。本明細書に記載された全ての文献、特許出願、および技術規格は、個々の文献、特許出願、および技術規格が参照により取り込まれることが具体的かつ個々に記された場合と同程度に、本明細書中に参照により取り込まれる。

Claims (11)

  1. 音声解析により嚥下状態を判定する判定装置であって、前記判定装置は:
    被験者が発した音声の時系列データである音声データの入力を受け付ける入力手段;
    前記入力手段で受け付けた音声データを解析する解析手段;
    前記解析手段による解析結果に基づいて被験者の嚥下状態を判定する判定手段、
    を備え、
    前記解析手段が、前記音声データに対して音素解析を行うことにより、前記音声データに含まれる音素を特定し、前記音声データに含まれる音素が属する音韻クラスを特定し、前記音声データの各時刻の音韻クラスの音素事後確率(Phonological posteriors)を計算し、
    前記判定手段が、前記解析手段によって得られた前記各時刻の音韻クラスの音素事後確率に基づいて、前記被験者の嚥下障害の度合いを判定する、
    判定装置。
  2. 前記解析手段は、更に、前記音声データを以下の式F(a)に示す音響特徴量を計算し、
    前記判定手段は、前記解析手段によって得られた、前記各時刻の音韻クラスの音素事後確率と、前記音響特徴量とに基づいて、前記被験者の嚥下障害の度合いを判定する、
    請求項1に記載の判定装置:

    (式中、gは嚥下状態の有無や進行の度合いを判定する線形または非線形モデルであり、xは前記音声データとして入力するフレーズに固有の係数であり、f(n)は音響パラメータであって、前記音響パラメータは、フォルマント周波数、メル周波数ケプストラム、周波数スペクトラム、音声のエンベロープ、波形の変動情報、ゼロ交差率、ハースト指数、および閉鎖開放から声帯振動が始まるまでの時間、から成る群から1つまたは複数が選択される)。
  3. 前記解析手段は、更に、前記音声データフォルマント周波数またはメル周波数ケプストラムに基づき作成された音響特徴量を取得し、
    前記判定手段は、前記解析手段によって得られた、前記各時刻の音韻クラスの音素事後確率と、前記音響特徴量とに基づいて、前記被験者の嚥下障害の度合いを判定する、
    請求項に記載の判定装置。
  4. 前記判定手段は、前記音韻クラスの音素事後確率を入力した場合に前記被験者の嚥下障害の度合いを出力するように機械学習処理が施されている機械学習モデルに対して、前記解析手段によって得られた前記音韻クラスの音素事後確率を入力することにより、前記機械学習モデルから出力される前記被験者の嚥下障害の度合いを取得し、前記被験者の嚥下障害の度合いを判定する、
    請求項に記載の判定装置。
  5. 前記入力手段は、更に、前記被験者が同一のフレーズを2回発した際の前記音声データの各々を受け付け、
    前記判定手段は、前記音響特徴量に基づいて前記被験者の嚥下障害の度合いを判定する際に、前記被験者によって発せられた2回の前記同一のフレーズの各々に対応する、2つの前記音声データから得られた前記音響特徴量の差分または平均値に基づいて、前記被験者の嚥下障害の度合いを判定する、
    請求項または請求項に記載の判定装置。
  6. 前記入力手段は、前記被験者が前記嚥下障害の度合いに関連するフレーズを発話した際の前記音声データの入力を受け付ける
    請求項1または請求項2に記載の判定装置。
  7. 前記入力手段は、前記被験者が

    の少なくとも1つの音を含む発話をした際の前記音声データの入力を受け付ける
    請求項1または請求項2に記載の判定装置。
  8. 前記解析手段は、更に、前記音声データに対して音声の強度解析を行うことにより作成された音響特徴量を取得し、
    前記判定手段は、前記解析手段によって得られた、前記各時刻の音韻クラスの音素事後確率と、前記音響特徴量とに基づいて、前記被験者の嚥下障害の度合いを判定する、
    請求項に記載の判定装置。
  9. 前記音声データは、前記被験者が所定フレーズの繰り返しを発話した際の音声データである、
    請求項1または請求項2に記載の判定装置。
  10. コンピュータが音声解析により嚥下状態を判定する方法であって、前記方法は:
    被験者が発した音声の時系列データである音声データの入力を受け付ける入力工程;
    前記入力工程で受け付けた音声データを解析する解析工程;
    前記解析工程による解析結果に基づいて被験者の嚥下状態を判定する判定工程、
    を備え、
    前記解析工程が、前記音声データに対して音素解析を行うことにより、前記音声データに含まれる音素を特定し、前記音声データに含まれる音素が属する音韻クラスを特定し、前記音声データの各時刻の音韻クラスの音素事後確率(Phonological posteriors)を計算し、
    前記判定工程が、前記解析工程によって得られた前記各時刻の音韻クラスの音素事後確率に基づいて、前記被験者の嚥下障害の度合いを判定する、
    処理をコンピュータが実行する方法。
  11. 請求項1または請求項2に記載の判定装置の各手段としてコンピュータを機能させるためのプログラム。
JP2023551872A 2021-09-29 2022-09-29 嚥下障害の判定装置および判定方法 Active JP7754432B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2021159606 2021-09-29
JP2021159606 2021-09-29
PCT/JP2022/036558 WO2023054632A1 (ja) 2021-09-29 2022-09-29 嚥下障害の判定装置および判定方法

Publications (3)

Publication Number Publication Date
JPWO2023054632A1 JPWO2023054632A1 (ja) 2023-04-06
JPWO2023054632A5 JPWO2023054632A5 (ja) 2024-06-19
JP7754432B2 true JP7754432B2 (ja) 2025-10-15

Family

ID=85782906

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2023551872A Active JP7754432B2 (ja) 2021-09-29 2022-09-29 嚥下障害の判定装置および判定方法

Country Status (2)

Country Link
JP (1) JP7754432B2 (ja)
WO (1) WO2023054632A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20250322842A1 (en) * 2022-04-18 2025-10-16 Panasonic Intellectual Property Management Co., Ltd. Oral function evaluation device, oral function evaluation system, and oral function evaluation method
CN120256885B (zh) * 2025-06-05 2025-09-16 南昌大学第一附属医院 一种基于吞咽训练仪的吞咽效果预测方法及系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019225242A1 (ja) 2018-05-23 2019-11-28 パナソニックIpマネジメント株式会社 摂食嚥下機能評価方法、プログラム、摂食嚥下機能評価装置および摂食嚥下機能評価システム

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5182773A (en) * 1991-03-22 1993-01-26 International Business Machines Corporation Speaker-independent label coding apparatus
JP3450411B2 (ja) * 1994-03-22 2003-09-22 キヤノン株式会社 音声情報処理方法及び装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019225242A1 (ja) 2018-05-23 2019-11-28 パナソニックIpマネジメント株式会社 摂食嚥下機能評価方法、プログラム、摂食嚥下機能評価装置および摂食嚥下機能評価システム

Also Published As

Publication number Publication date
JPWO2023054632A1 (ja) 2023-04-06
WO2023054632A1 (ja) 2023-04-06

Similar Documents

Publication Publication Date Title
Montaña et al. A Diadochokinesis-based expert system considering articulatory features of plosive consonants for early detection of Parkinson’s disease
Byrd Articulatory vowel lengthening and coordination at phrasal junctures
JP5120826B2 (ja) 発音診断装置、発音診断方法、記録媒体、及び、発音診断プログラム
Lu et al. Speech tasks and interrater reliability in perceptual voice evaluation
JP7754432B2 (ja) 嚥下障害の判定装置および判定方法
Mittapalle et al. Glottal flow characteristics in vowels produced by speakers with heart failure
Drugman et al. Tracheoesophageal speech: A dedicated objective acoustic assessment
CN106157212A (zh) 一种基于ema的发音障碍中文评估方法
JP7316596B2 (ja) 口腔機能可視化システム、口腔機能可視化方法およびプログラム
CN110870765A (zh) 声门闭合实时测量和视听反馈技术的嗓音治疗仪器及方法
Nissen et al. First and second language tongue movements in Spanish and Korean bilingual speakers
Yang et al. Acoustic development of vowel production in native Mandarin-speaking children
Liker et al. Tongue-Palate Contact Timing during/s/and/z/in English
Searl Magnitude and variability of oral pressure in tracheoesophageal speech
Osfar Articulation of whispered alveolar consonants
Brogan Demystifying Salvadoran [s θ]
WO2022254973A1 (ja) 口腔機能評価方法、プログラム、口腔機能評価装置および口腔機能評価システム
Peña The production and perception of stød in Danish
JP7745214B2 (ja) 音声特徴量算出方法、音声特徴量算出装置、及び、口腔機能評価装置
Penney et al. Electroglottographic analysis of coda voicelessness in Australian English
Olson Acoustic properties of bilabial trills in Mangbetu
JP7592948B2 (ja) 口腔機能評価方法、口腔機能評価装置及びプログラム
Jesus et al. Attributes Associated with Consonantal Place and Voicing in Whispered Speech
Luo The Articulatory and Acoustic Variability in Putonghua Onset/r
Haque et al. Effect of Inter-Speaker Hypernasality Variation in Cleft Palate on Acoustic Vowel Space

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240610

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240610

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20240626

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20240626

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240819

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20250218

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20250417

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20250507

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20250707

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20250827

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20250909

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20250922

R150 Certificate of patent or registration of utility model

Ref document number: 7754432

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150