(57)【要約】  在宅等で非医療者が努力呼吸に特徴的な呼吸補助筋の動員を、簡便に検出かつ客観的に表現する装置、方法、プログラムを提供する。人の頸部を含む動画像から呼吸補助筋が存在する領域を抽出し、その領域内の形状的特徴を抽出する。抽出された形状特徴が、呼吸に由来するか否かを判断し、呼吸に由来するものである場合には、その形状的特徴の変化から、呼吸状態を表す指標を算出する。
2020226182.tif 000002
【特許請求の範囲】
【請求項１】
人の頸部を含む動画像から呼吸補助筋が存在する呼吸補助筋存在領域を抽出する呼吸補助筋領域検出装置と、
前記呼吸補助筋存在領域中の形状的特徴を抽出する動画像処理装置と、
前記形状的特徴の時系列的変動から呼吸補助筋の動員を判断する呼吸動作判定装置と、
を備えた呼吸状態検出装置。
【請求項２】
前記呼吸補助筋領域検出装置は、胸鎖乳突筋及び頸窩あるいは鎖骨上窩に基づいて前記呼吸補助筋存在領域を決定する、請求項１に記載の呼吸状態検出装置。
【請求項３】
前記呼吸補助筋領域検出装置は、前記動画像から人体の輪郭線を抽出するとともに、前記輪郭線上の人体の左右の首付根位置及び肩峰位置を抽出し、抽出された前記左右の首付根位置及び肩峰位置に基づいて呼吸補助筋存在領域を抽出する、請求項１または２に記載の呼吸状態検出装置。
【請求項４】
前記動画像処理装置は、前記呼吸補助筋存在領域の形状的特徴として、鎖骨上窩および／または頸窩部に対する胸鎖乳突筋の隆起を抽出するとともに、
前記呼吸動作判定装置は、前記胸鎖乳突筋の隆起の時系列的変動の振幅が所定値以上であり、かつ前記胸鎖乳突筋の隆起の時系列的変動の周期が所定の範囲内であるとき、呼吸補助筋が動員されたと判断する、請求項１乃至３のいずれか一項に記載の呼吸状態検出装置。
【請求項５】
前記動画像処理装置は、前記呼吸補助筋存在領域の形状的特徴が抽出されない場合、前記動画像中の肩部を含む動画像から、当該肩部の動きの時系列的変動を抽出するとともに、
前記呼吸動作判定装置は、前記肩部の動きの時系列的変動の振幅と前記肩部の動きの時系列的変動の周期に基づいて、前記肩部の動きの時系列的変動により、呼吸補助筋動員を判断する、請求項１乃至４のいずれか一項に記載の呼吸状態検出装置。
【請求項６】
前記動画像処理装置および前記呼吸動作判定装置の出力情報を基に呼吸状態を評価する呼吸評価演算装置をさらに備える、請求項１乃至５のいずれか一項に記載の呼吸状態検出装置。
【請求項７】
前記呼吸評価演算装置が評価した前記呼吸状態を継続的に記録する記録装置をさらに備える、請求項６に記載の呼吸状態検出装置。
【請求項８】
人の頸部を含む動画像から呼吸補助筋が存在する呼吸補助筋存在領域を抽出する呼吸補助筋領域検出ステップと、
前記呼吸補助筋存在領域中の形状的特徴を抽出する動画像処理ステップと、
前記形状的特徴の時系列的変動から呼吸補助筋の動員を判断する呼吸動作判定ステップと、
を含む呼吸状態検出方法。
【請求項９】
コンピュータを動作させるプログラムであって、
人の頸部を含む動画像から呼吸補助筋が存在する呼吸補助筋存在領域を抽出する呼吸補助筋領域検出ステップと、
前記呼吸補助筋存在領域中の形状的特徴を抽出する動画像処理ステップと、
前記形状的特徴の時系列的変動から呼吸補助筋の動員を判断する呼吸動作判定ステップと、
を含む呼吸状態検出方法をコンピュータに実行させる呼吸状態検出プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、非接触で被験者の呼吸状態を検出する呼吸状態検出装置、呼吸状態検出方法および呼吸状態検出プログラムに関する。
【背景技術】
【０００２】
正常安静時の呼吸では、吸気は主に横隔膜の収縮によって行われ、呼気は筋肉を用いず、伸展された肺の受動的反跳によって行われる。一方、心不全などで肺機能の低下が起き、酸素摂取能が低下すると、その代償反応として頸部にある胸鎖乳突筋、斜角筋などの呼吸補助筋が吸気時に動員されて胸腔容積を増加させる努力呼吸が行われる。また、ＣＯＰＤなどの呼吸器患者では横隔膜の筋力が低下するため、常に胸鎖乳突筋を動員した努力呼吸様式となり、胸鎖乳突筋の肥大がみられる。このように、努力呼吸時には呼吸補助筋が動員されるような特徴的な動きがあることは医療者の間ではよく知られている。この呼吸補助筋の活動を計測する手段としては、これまで、筋電計測による方法が提案されている（非特許文献１、非特許文献２）。【０００３】
一方、非接触の呼吸計測技術としては、特許文献１がある。ここでは、両肩の中央部を領域の上限、右肩部を領域の右限、左肩部を領域の左限、及び腰中央部を領域の下限として定められる矩形領域を呼吸主要領域として確定し、その領域の距離情報から人の呼吸波形を生成する方法が提案されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１７−２１７２９８号公報
【非特許文献】
【０００５】
【非特許文献１】「表面筋電図を使用した補助呼吸筋の活動分析による努力呼吸の評価」，人工呼吸，  ３０−１，  （２０１３）
【非特許文献２】「吸気流量の増加を伴う頻呼吸時における吸気筋活動の評価法の検討」，人工呼吸，  ３４−２，  （２０１７）
【発明の概要】
【０００６】
これまで努力呼吸の有無の判断は、通常、専門医の視診によっており、非医療者がその状況を客観的に表現する手段が存在しなかった。上述のように、呼吸補助筋活動を計測する方法としては、従来、筋電計測を用いるものが開発されている。しかし、筋電によって筋活動を計測するには、電極の貼り付けを常に正しい位置に行う必要があること、また、その評価には異なる活動時の筋電振幅に対する比率を求める必要があることから、専門的知識が不可欠である。従って、医療者の補助なく呼吸補助筋の筋電計測を実施するのは困難であり、在宅で、いつでも、どこでも簡便に計測できる方法ではなかった。
【０００７】
また、従来開発されている非接触の呼吸計測装置は、人が静止している場合の呼吸数計測には有効な方法ではあるが、努力呼吸時に局所に現れる特徴的な動きを計測する方法ではなかった。
【０００８】
すなわち、在宅等で心不全等の慢性疾患の増悪の予兆を客観的な手法で捉えることは重要と認識されつつも、非医療者が努力呼吸に特徴的な動きを簡便に計測する装置、方法、プログラムはこれまで開発されていなかった。
【０００９】
上記課題を解決するために、本発明に係る呼吸状態検出装置は、人の頸部を含む動画像から呼吸補助筋が存在する呼吸補助筋存在領域を抽出する呼吸補助筋領域検出装置と、呼吸補助筋存在領域中の形状的特徴を抽出する動画像処理装置と、形状的特徴の時系列的変動から呼吸補助筋の動員を判断する呼吸動作判定装置とを備える。
【００１０】
呼吸補助筋領域検出装置は、胸鎖乳突筋及び頸窩あるいは鎖骨上窩に基づいて呼吸補助筋存在領域を決定することが好ましい。
【００１１】
呼吸補助筋領域検出装置は、動画像から人体の輪郭線を抽出するとともに、輪郭線上の人体の左右の首付根位置及び肩峰位置を抽出し、抽出された左右の首付根位置及び肩峰位置に基づいて呼吸補助筋存在領域を抽出することが好ましい。
【００１２】
また、動画像処理装置は、呼吸補助筋存在領域の形状的特徴として、鎖骨上窩および／または頸窩部に対する胸鎖乳突筋の隆起の時系列的変動を抽出するとともに、呼吸動作判定装置は、胸鎖乳突筋の隆起の時系列的変動の振幅が所定値以上であり、かつ胸鎖乳突筋の隆起の時系列的変動の周期が所定の範囲内であるとき、形状的特徴の時系列的変動が呼吸に由来すると判断することが好ましい。
【００１３】
ここで、呼吸補助筋存在領域の形状的特徴が抽出されない場合、動画像処理装置は、さらに、動画像中の肩部を含む動画像から、当該肩部の動きの時系列的変動幅を利用して呼吸状態を評価する機能を有してもよい。具体的には、肩部の動きの時系列的変動の振幅と肩部の動きの時系列的変動の周期に基づいて、肩部の動きの時系列的変動により、呼吸補助筋動員を判断することが好ましい。
【００１４】
また、本発明に係る呼吸状態検出装置は、動画像処理装置および呼吸動作判定装置の出力情報を基に呼吸状態を評価する呼吸評価演算装置、さらに、評価した呼吸状態を記録する結果記録装置を備え、継続的に得られた呼吸状態の評価値から平常の状態とそうでない状態を判別する機能を有してもよい。
【００１５】
また、本発明に係る呼吸状態検出装置は、さらに、呼吸評価演算装置の評価結果を数値、または／および、グラフとして表示すること、および／または、計測結果記録装置に記録されたデータを呼吸評価演算装置の評価結果と比較して表示する機能を有しても良い。
【００１６】
本発明に係る呼吸状態検出方法は、人の頸部を含む動画像から呼吸補助筋が存在する領域を抽出する呼吸補助筋領域検出ステップと、呼吸補助筋存在領域中の形状的特徴を抽出する動画像処理ステップと、形状的特徴の時系列的変動から呼吸補助筋の動員を判断する呼吸動作判定ステップと、を含む。
【００１７】
本発明に係る呼吸状態検出プログラムは、コンピュータを動作させるプログラムであって、人の頸部を含む動画像から呼吸補助筋が存在する領域を抽出する呼吸補助筋領域検出ステップと、呼吸補助筋存在領域中の形状的特徴を抽出する動画像処理ステップと、形状的特徴の時系列的変動から呼吸補助筋の動員を判断する呼吸動作判定ステップと、を含む呼吸動作判定方法をコンピュータに実行させるプログラムである。
【００１８】
本発明の呼吸状態検出装置、方法、プログラムを用いれば、在宅等で非医療者が努力呼吸に特徴的な呼吸補助筋の動員を簡便に検出できる。
【図面の簡単な説明】
【００１９】
【図１】本発明の装置の設置例である。
【図２】人の頸部の形状図である。
【図３】第一の実施形態に係る呼吸状態検出装置のシステムの構成図である。
【図４】第一の実施形態に係る呼吸補助筋領域の決定フロー図である。
【図５】第一の実施形態に係る呼吸補助筋領域検出の説明図である。
【図６】第一の実施形態に係る胸鎖乳突筋指標算出の説明図である。
【図７】第一の実施形態に係る胸鎖乳突筋指数（ＳＣＭ＿ｉｎｄｅｘ）の時系列データ算出フロー図である。
【図８】第一の実施形態に係る呼吸補助筋動員判定フロー図である。
【図９】種々の呼吸状態における鎖骨上窩の奥行の時間的変動を表すタイミングチャートである。
【図１０】第二の実施形態に係る呼吸状態検出装置の呼吸補助筋動員判定フロー図である。
【図１１】第三の実施形態に係る呼吸状態検出装置のシステム構成図である。
【発明を実施するための形態】
【００２０】
以下、図面を参照しながら本発明の実施の形態を詳細に説明する。全図面に渡り、対応する構成要素には共通の参照符号を付す。
【００２１】
図１は本発明の実施形態に係る呼吸状態検出装置２００の設置例である。慢性心不全、慢性呼吸不全患者Ｐは、状態が悪化すると軽い労作でも息切れや呼吸困難が起こり、努力呼吸が出現しやすくなる。前述の通り、努力呼吸の吸気時には呼吸補助筋が動員される特徴的な動きがあるため、これを日常生活の中で、手間なく、かつ、患者の負担なく捉えられれば、早期に増悪予兆を把握できると考えられる。そのために、労作が起こりやすい場所、例えば、トイレ３００、脱衣所等で呼吸補助筋付近の体の動きを捉えて指標化し、これを日々、患者Ｐに負荷をかけることなく継続的に記録することで、通常と異なるデータとなった場合に患者Ｐを観察する医療関係者や患者Ｐに対して医療行為を行う病院４００等に対してアラートＡを出すようなシステムが考えられる。本実施形態に係る呼吸状態検出装置２００の目的の一つは、このような呼吸専門医が努力呼吸を視診で判断する際の呼吸補助筋の動員を画像センサーと情報処理技術で指標化し、誰もが簡単に判断できるようにするものである。
【００２２】
一方、呼吸専門医は肩の動きにも注目する。すなわち、肩を上下させて息をする肩呼吸では、安静呼吸時に比べ、大きな肩の動きが呼吸周期で見られる。この場合も呼吸補助筋が動員されたと判定できる。
【００２３】
［第一の実施形態に係る呼吸状態検出の原理］
図２は人の頸部の形状図である。努力呼吸の吸気時に活動する呼吸補助筋は、頸部から鎖骨上部にかけて存在する。呼吸補助筋には胸鎖乳突筋３０１、斜角筋など色々あるが、視診で分かりやすいのは胸鎖乳突筋３０１である。この筋肉は両側頭部から胸骨中央上部の胸骨柄に伸び、頭部の動きの他、胸骨と鎖骨の挙上にかかわっている。呼吸に関連する動きは後者であり、胸骨の中心線を基準に左右同時に緊張する。胸鎖乳突筋３０１は、頸部の耳の下から胸骨切痕３０２へ向かう斜めの領域にあり、重度の努力呼吸の吸気時には、左右の胸鎖乳突筋３０１が隆起し、胸骨上部の頸窩部３０３や、胸鎖乳突筋３０１と鎖骨３０４との間の鎖骨上窩３０５で陥没が見られる。そこで、例えば、胸鎖乳突筋３０１と頸窩部３０３あるいは鎖骨上窩３０５との深度方向の差が呼吸周期に合わせ変動するか否かが確認できれば、呼吸補助筋動員の有無が判定できる。
【００２４】
［第一の実施形態に係る呼吸状態検出装置］
図３は第一の実施形態に係る呼吸状態検出装置のシステム構成図である。呼吸状態検出装置２００は、メモリ１と、制御部２と、通信インターフェース（Ｉ／Ｆ）３と、を備え、これらは内部バス４により接続されている。通信Ｉ／Ｆ３は、病院４００（図１参照）等に設置された端末やサーバ（図示せず）に情報を送信することができる。制御部２は、プロセッサ（ＣＰＵ）を含む。制御部２は、メモリ１に格納された制御プログラムを実行することにより、呼吸状態検出装置２００の制御を行なう。制御部２は、呼吸補助筋領域検出装置２０１と、動画像処理装置２０２と、呼吸動作判定装置２０３と、呼吸評価演算装置２０４と、を備えている。呼吸補助筋領域検出装置２０１、動画像処理装置２０２、呼吸動作判定装置２０３及び呼吸評価演算装置２０４は、プログラムモジュールとして、一つのプロセッサで実現されてもよいし、それぞれ一つまたは複数のプロセッサと、周辺回路とを有する構成としてもよい。
【００２５】
呼吸補助筋領域検出装置２０１は呼吸補助筋が存在する呼吸補助筋存在領域を入力画像から検出する。呼吸補助筋領域検出装置２０１は、胸鎖乳突筋３０１及び頸窩部３０３あるいは鎖骨上窩３０５に基づいて呼吸補助筋存在領域を決定することができる。呼吸補助筋領域検出装置２０１には、人体の頸部を含む動画像が入力される。呼吸補助筋の動員を判定するには、この動画像から頸部にある胸鎖乳突筋３０１の隆起、すなわち、奥行方向の形状変化を捉える必要がある。胸鎖乳突筋の隆起は、胸鎖乳突筋３０１に対する頸窩部３０３の奥行方向の形状変化、あるいは胸鎖乳突筋３０１に対する鎖骨上窩３０５の奥行方向の形状変化を観測すればよい。呼吸補助筋領域検出装置２０１に入力する動画像としては、Ｌｉｇｈｔ  Ｃｏｄｉｎｇ方式やＴｏＦ（Ｔｉｍｅ  ｏｆ  Ｆｌｉｇｈｔ）画像センサーを搭載したデバイス（例えば、３Ｄセンサー１００）の出力を用いることができる。これらは、直接的に奥行方向の形状データを動画像として出力できるため、そのままの形式でデータを入力すればよい。なお、他に、ステレオカメラ像や単眼の影付き動画像を用いても良いが、その場合は、奥行情報を算出するための周知の変換が必要な場合がある。いずれにしても、以下では、奥行方向のデータを含む動画像が入力されたものとして説明する。
【００２６】
呼吸補助筋領域検出装置２０１は、図４に示す手順を用いて、入力された人体の頸部を含む動画像から呼吸補助筋が存在する領域を検出する。
まず呼吸補助筋領域検出装置２０１は、入力された人体の頸部を含む動画像から、周知の方法により図５に示す人の上半身の輪郭４０１（点線）を求める（ステップＳ５０１）。
次に、呼吸補助筋領域検出装置２０１は、ステップＳ５０１で求めた輪郭線上の首付根位置４０２と肩峰の位置４０３を取得する（ステップＳ５０２）。首付根位置４０２や肩峰４０３は、人の上半身の輪郭で曲率の変化が大きい箇所であり、その画像的特徴から容易に検出できる。例えば、輪郭を多角形近似した頂点から選ぶこともできる。
【００２７】
そして、呼吸補助筋領域検出装置２０１は、ステップＳ５０２で取得した左右の首付根位置４０２及び肩峰４０３の４点から、左右の胸鎖乳突筋、鎖骨上窩、頸窩がそれぞれ含まれる頸部の計測領域を決定する（ステップＳ５０３）。産業技術総合研究所人体寸法データベース１９９１−９２によると、頚椎点と頸窩点の間隔が最大４ｃｍ、肩峰点と頸窩点の間隔が最大３ｃｍ、胸骨丙の幅が大きく見積もっても６ｃｍ程度であることから、例えば縦方向（または、体軸に沿った方向）は、首付根位置の高さから３ｃｍ下の位置４０４を上辺、肩峰位置高さ４０５を下辺とし、横方向（または、体軸に垂直な方向）は左右辺を左右両肩峰点の中点から左右に３ｃｍとする矩形領域４０８を設定することができる。なお上述の各数値はあくまで一例であり、これに限られるものではなく、少なくとも、頸窩、または、胸骨上窩が含まれる領域が含まれていればよい。この矩形領域４０８を左右両肩峰点の中点で分離することで、２つの計測領域４０６および４０７が得られる。これにより、呼吸補助筋が存在する２つの計測領域４０６および４０７（呼吸補助筋存在領域）それぞれの領域内で左右夫々の胸鎖乳突筋の状態を確認することができる。
【００２８】
なお、上記は輪郭４０１から求める方法だが、これ以外にも頸窩部付近の鎖骨等の形状特徴を周知の方法で捉えることにより計測領域を決定してもよい。
【００２９】
動画像処理装置２０２は動画像処理装置であって、呼吸補助筋領域検出装置２０１で決定した呼吸補助筋存在領域における形状的特徴として胸鎖乳突筋突出度を算出する。前述の通り、努力呼吸の吸気時には左右両方の胸鎖乳突筋が隆起し、相対的に鎖骨上窩、頸窩部が窪む。医療者が視診する場合は、この胸鎖乳突筋の突出度合い、すなわち、より目立つ突出であるか否かを直感的に判断すると考えられる。そこで動画像処理装置２０２は、図６に示すある動画像の１フレーム内の計測領域の横方向のスキャンライン６０５に沿った横断面６０６及び図７に示す手順を用い、頸窩部左側の計測領域６０４の胸鎖乳突筋６０１突出度算出例について説明する。
【００３０】
まず動画像処理装置２０２は、スキャンライン６０５に沿って左側計測領域６０４をＸ方向（横方向）に走査し、その線上の頸窩部６０３の深度：Ｆｚｙ、および、そのフレーム内で最も隆起した箇所の深度（深度の最小箇所）Ｓｚｙとの差（Ｆｚｙ−Ｓｚｙ）を求め（ステップＳ７０１）、かつ、その値を頸窩部の水平位置（Ｘ座標）：Ｆｘｙと最も隆起した箇所の水平位置：Ｓｘｙとの距離（｜Ｆｘｙ−Ｓｘｙ｜）で除す（ステップＳ７０２）。この値は、上下方向の位置がｙにおける胸鎖乳突筋６０１の突出度であり、値が大きければ、それだけ目立つ突出があることになる。さらに、これを計測領域の下（ｙｍｉｎ）から上（ｙｍａｘ）まで積算する（ステップＳ７０３）ことで、突出の縦方向の広がりを示すことができる。ここで求めた積算値を、そのフレームにおける胸鎖乳突筋指数（ＳＣＭ＿ｉｎｄｅｘ）と呼ぶことにする。数式で表現すると式１のようになる。
【００３１】
【数１】
2020226182.tif 000003
但し、
ｙｍａｘ：  計測領域上端のＹ座標
ｙｍｉｎ：  計測領域下端のＹ座標
Ｆｘｙ：  Ｙ座標がｙの時の頸窩部のＸ座標値
Ｆｚｙ：  Ｙ座標がｙの時の頸窩部の深度（Ｚ座標値）
Ｓｘｙ：  Ｙ座標がｙの時の最大突出部のＸ座標値
Ｓｚｙ：  Ｙ座標がｙの時の最大突出部の深度（Ｚ座標値）
【００３２】
これを左右両方の計測領域４０６および４０７で毎フレームごと行うことで、鎖骨上窩および／または頸窩部に対する胸鎖乳突筋の隆起の時系列的変動をＳＣＭ＿ｉｎｄｅｘの時系列データとして表すことができる（ステップＳ７０４）。なお、動画像処理装置２０２は、計測領域内の画像的特徴点をＡＫＡＺＥ特徴量等を用いて求め、その画像的特徴点の毎フレームでの移動から計測領域を継続的に追跡する構成とすることも可能である。
【００３３】
呼吸動作判定装置２０３は、動画像処理装置２０２で求めたＳＣＭ＿ｉｎｄｅｘ（胸鎖乳突筋指数）の時間的変動が、呼吸補助筋動員によるものであるか否かを判定する。前述の通り、努力呼吸時に胸鎖乳突筋の隆起が見られるのは左右の胸鎖乳突筋が一定以上の大きさで同時に働き、かつ、その周期が呼吸周期に一致する場合である。また、隆起が現れるのは特に吸気時であるため、呼吸の位相の前半部分に大きな隆起が現れる特徴を持つ。
【００３４】
図８は、呼吸動作判定装置２０３で実行される処理の流れを示したフローチャートである。まず、呼吸動作判定装置２０３は、動画像処理装置２０２で求めた左右の計測領域それぞれのＳＣＭ＿ｉｎｄｅｘ時系列データの振幅と周期を周波数解析、ピーク検出等周知の方法により求める（ステップＳ８０１）。
【００３５】
次に呼吸動作判定装置２０３は、左右の計測領域それぞれのＳＣＭ＿ｉｎｄｅｘの振幅が所定値以上であり、かつ、この周期が呼吸周期であるか否かを判定する（ステップＳ８０２）。この所定値については、病状が安定している場合の安静呼吸時に計測したＳＣＭ＿ｉｎｄｅｘの振幅と周期を用いれば良い。ＳＣＭ＿ｉｎｄｅｘの振幅が所定値以下の場合は呼吸以外の体動や計測ノイズによるものと考えられるので、この場合はＳＣＭ＿ｉｎｄｅｘの変動は呼吸に由来しないと判断することができる。また、例えば、一定区間のＳＣＭ＿ｉｎｄｅｘの時系列データに自己相関解析を行ってデータの周期性を解析し、求めた周期が成人健常者の一般的な呼吸周期（例えば２〜６秒）に入るかどうかを判定することで、ＳＣＭ＿ｉｎｄｅｘの変動が呼吸に由来するか否かを判定することができる。また、振幅は、呼吸周期の波形のピーク、すなわち、山と谷を求め、その差から求められる。
【００３６】
ステップＳ８０２の判定結果がＮｏの場合は、ステップＳ８０５に進み、呼吸補助筋の動員はないものと判断する。
【００３７】
ステップＳ８０２の判定結果がＹｅｓの場合は、ステップＳ８０３に進み、左右のＳＣＭ＿ｉｎｄｅｘ時系列データの相関係数ｒを求め、計算されたｒが、予め決められた一定値以上か否かを判定する（ステップＳ８０４）。例えば、ｒが０．７以上であれば、左右の胸鎖乳突筋が同時に働いており、呼吸補助筋が動員されたと判定することができる。ステップＳ８０４の判定結果がＮｏの場合は、ステップＳ８０６へ進み、呼吸補助筋の動員はないものと判断する。ステップＳ８０４の判定結果がＹｅｓの場合は呼吸補助筋の動員がされているものと判断する（ステップＳ８０７）。この場合は、好ましくは、呼吸評価演算装置２０４において呼吸補助筋の動員度合いを評価する。
【００３８】
呼吸評価演算装置２０４は、呼吸補助筋の動員の度合いから、ある時間の努力呼吸を指標化する。ここでは、例えば、呼吸動作判定装置２０３によって呼吸補助筋の動員が判定された場合、動画像処理装置２０２で求めたＳＣＭ＿ｉｎｄｅｘの計測時間：Ｔにわたる平均値（努力呼吸指標）を求める。数式で表現すると式２のようになる。
【００３９】
【数２】
2020226182.tif 000004
但し、
ＳＣＭ＿ｉｎｄｅｘ＿ｌｅｆｔ（ｔ）：時刻ｔにおける左側ＳＣＭ＿ｉｎｄｅｘ
ＳＣＭ＿ｉｎｄｅｘ＿ｒｉｇｈｔ（ｔ）：時刻ｔにおける右側ＳＣＭ＿ｉｎｄｅｘ
【００４０】
例えば、日々、この努力呼吸指標が予め定めた値よりも大きい場合や増加傾向にある場合は、努力呼吸が増す、すなわち、普段よりも息苦しいと感じている可能性がある。
【００４１】
上記の説明では胸鎖乳突筋指数（ＳＣＭ＿ｉｎｄｅｘ）の時間的変動から努力呼吸の有無を判定する例を示した。次に、図３に示した呼吸状態検出装置２００を用いて、鎖骨上窩の奥行方向の変動量から努力呼吸の有無を判定する方法について説明する。被験者が努力呼吸を行う場合には、鎖骨上窩が凹むことが分かっている。従って、鎖骨上窩の奥行方向の変動量を観測することにより、被験者が努力呼吸を行っているか否かを判定することができる。
【００４２】
まず、被験者の胸鎖乳突筋及び鎖骨上窩を含む頸部を観察できるように、被験者の正面に図３に示す３Ｄセンサー１００を設置し、動画像を撮像する。
【００４３】
次に、呼吸補助筋領域検出装置２０１が、３Ｄセンサー１００が撮像した、被験者の頸部を含む動画像から呼吸補助筋が存在する呼吸補助筋存在領域を抽出する。具体的には、図５に示すように、撮像した画像を人体形状のデータベースと照合し、呼吸補助筋存在領域（４０６、４０７）を抽出する。
【００４４】
次に、動画像処理装置２０２が、図４に示したフローに基づいて、呼吸補助筋存在領域（４０６、４０７）中の形状的特徴を抽出する。具体的には、図６に示すように、胸鎖乳突筋６０１及び鎖骨上窩６０２の測定点を決定する。次に、被験者が、通常呼吸、腹式呼吸、胸式呼吸、無呼吸、及び努力呼吸を行った場合の胸鎖乳突筋６０１及び鎖骨上窩６０２の位置を、３Ｄセンサー１００を用いて測定する。
【００４５】
ここで、「通常呼吸」とは、被験者が安静時に行う呼吸をいう。「腹式呼吸」とは、一般的には胸郭（肋骨等からなる籠状の骨格）をなるべく動かさずに行う呼吸をいう。「胸式呼吸」とは、胸の周辺で行う呼吸をいう。「無呼吸」とは、気道の空気の流れが止まった状態をいう。「努力呼吸」とは、安静時呼吸では使用されない呼吸筋を動員して行う呼吸であり、具体的には、呼吸困難のため、吸気時に胸鎖乳突筋等の呼吸補助筋を動かしたり、呼気時に内肋間筋や腹筋等を動かしたりして、努力的に行なう呼吸をいう。
【００４６】
次に、呼吸動作判定装置２０３が、動画像処理装置２０２が抽出した形状的特徴の時系列的変動から呼吸補助筋の動員の有無を判断する。具体的には、図７に示したフローと同様の手順により、３Ｄセンサー１００を設置した位置から、各呼吸状態における被験者の胸鎖乳突筋６０１までの距離と、鎖骨上窩６０２までの距離との差を、３Ｄセンサー１００を用いて測定し、胸鎖乳突筋６０１の位置に対する鎖骨上窩６０２の相対的な位置を算出する。観察時間は、各呼吸状態それぞれについて約１０秒間とした。
【００４７】
図９（ａ）〜（ｅ）に、図３に示した呼吸状態検出装置２００を用いて測定した、種々の呼吸状態における鎖骨上窩の奥行方向の変動量の時間的変化を表すタイミングチャートを示す。鎖骨上窩の奥行方向の変動量とは、胸鎖乳突筋の位置に対する鎖骨上窩の相対的な位置の変化をいう。従って、被験者の上体の揺動等により胸鎖乳突筋の位置と鎖骨上窩の位置が同じだけ変動する場合は、鎖骨上窩の奥行方向の変動量はないものとする。図９（ａ）〜（ｅ）において、（ａ）は通常呼吸、（ｂ）は腹式呼吸、（ｃ）は胸式呼吸、（ｄ）は無呼吸、（ｅ）は努力呼吸にそれぞれ対応する。尚、（ｅ）は、努力呼吸時に生じる呼吸補助筋の動員を健常人で示すため、被検者が鼻すすりを行った状態での検知結果を示す。
【００４８】
各タイミングチャートの縦軸は、鎖骨上窩の奥行方向の変動量［ｍｍ］を表し、横軸は時間［ｓｅｃ］を表す。図９（ａ）〜（ｄ）に示すように、努力呼吸以外の場合は、鎖骨上窩の奥行方向の変動量ｄ1〜ｄ4は最大で３〜５［ｍｍ］程度である。これに対して、図９（ｅ）に示すように、努力呼吸の場合は、鎖骨上窩の奥行方向の変動量ｄ5は１０［ｍｍ］程度であり、努力呼吸以外の場合と比べて大きな変動が観測されることが分かる。ここで、努力呼吸以外の呼吸で観測される鎖骨上窩の奥行方向の変動量の最大値と、努力呼吸で観測される鎖骨上窩の奥行方向の変動量の大きさとの間の所定の値（例えば、７［ｍｍ］）を閾値として設定することにより、鎖骨上窩の奥行方向の変動量の大きさが、閾値未満であれば努力呼吸ではないと判定し、閾値以上であれば努力呼吸であると判定することができる。
【００４９】
このように、鎖骨上窩の奥行方向の変動量を観測することにより、被験者が努力呼吸を行っているか否かを判定することができる。なお、図９（ｅ）に示した努力呼吸の場合、図８に示した胸鎖乳突筋指数（ＳＣＭ＿ｉｎｄｅｘ）を用いた判定方法においても呼吸補助筋の動員有りと判断され、努力呼吸であると判定することができる。
【００５０】
図９（ａ）〜（ｅ）は、胸鎖乳突筋６０１（図６参照）に対する鎖骨上窩６０２の奥行方向の変動量を観測した例を示したが、図６に示すように胸鎖乳突筋６０１に対する頸窩部６０３の奥行方向の変動量を観測することにより、努力呼吸の有無を判定するにしてもよい。
【００５１】
本実施形態によれば、呼吸専門医が努力呼吸を視診で判断する際の呼吸補助筋の動員を画像センサーと情報処理技術で指標化し、誰もが簡単に判断できるようになるため、在宅等で非医療者が努力呼吸に特徴的な呼吸補助筋の動員を簡便に計測かつ客観的に表現でき、在宅慢性疾患患者の病態をより正確かつ効率的に医療者に伝えることが可能になる。従って、医療者は、この情報を他のモニタリング情報やカルテ情報に照らし合わせて、患者の急性増悪の予兆を早期に捉えることが可能となる。
【００５２】
［第二の実施形態に係る呼吸状態検出原理］
上述の第一の実施形態では、呼吸補助筋の動員を人体の頸部を含む動画像から判定した。一方、頸部の肉付きや姿勢により、頸部まわりの動画像では呼吸補助筋の動員が判定しにくい場合がある。本実施形態は、このような場合に対応するために考えられたものであり、ここでは、努力呼吸の有無を肩の動きから判定する装置をさらに加える。前述の通り、努力呼吸時には、肩が呼吸と同期して大きく上下することが知られている。但し、肩は随意的にも上下に動かすことができる。従って、肩の動きから努力呼吸の有無を判定するには、まず、肩の動きが呼吸由来のものであるか否かを判定し、かつその動きの大きさが通常の呼吸時よりも大きいか否かを判断する必要がある。安静にもかかわらず、通常よりも大きな肩の動きがあれば、努力呼吸が発生している可能性が疑われる。すなわち、第一の実施形態では呼吸動作判定装置は頸部の胸鎖乳突筋の動員を判定したが、第二の実施形態では、胸鎖乳突筋の動員が判定できなかった場合に、上記のような肩―頸相対距離による呼吸動作判定を行う。
【００５３】
［第二の実施形態に係る呼吸状態検出装置］
第二の実施形態に係る呼吸状態検出装置２００´は第一の実施形態に係る呼吸状態検出装置２００とハードウェア構成としては同一であるため、呼吸状態検出装置２００´のシステム構成図は図３で代用する。第二の実施形態に係る呼吸状態検出装置２００´は、呼吸補助筋領域検出装置２０１´と、動画像処理装置２０２´と、呼吸動作判定装置２０３´と、呼吸評価演算装置２０４´とを備えている。第二の実施形態に係る呼吸状態検出装置２００´の呼吸補助筋領域検出装置２０１´、動画像処理装置２０２´、呼吸動作判定装置２０３´及び呼吸評価演算装置２０４´は、それぞれ第一の実施形態に係る呼吸状態検出装置２００の呼吸補助筋領域検出装置２０１、動画像処理装置２０２、呼吸動作判定装置２０３、呼吸評価演算装置２０４と同様の処理を実行するとともに、胸鎖乳突筋の動員が判定できなかった場合に、肩の動きから呼吸努力の有無を判定する機能を付加したものである。以下では、上記付加機能について説明する。
【００５４】
呼吸補助筋領域検出装置２０１´は呼吸補助筋が存在する領域を入力画像から検出する。第一の実施形態とは異なり、呼吸補助筋領域検出装置２０１´に入力される入力動画像は頸部だけでなく、肩部も含んだものとなる。
【００５５】
動画像処理装置２０２´は、肩の上下動を、例えば、上述の方法で求めた肩峰４０３と首付根位置４０２との上下方向の相対距離４０９の変動から求める。なお、体格による差をなくすため、この相対距離は、実際の寸法を用いるのではなく、例えば、首付根幅（左右の首付根位置４０２の間隔）や肩峰幅（左右の肩峰４０３の間隔）等の寸法で除し、正規化してもよい。
【００５６】
図１０は、呼吸動作判定装置２０３´で実行される処理を説明するためのフローチャートである。呼吸動作判定装置２０３´は、第一の実施形態と同様に動画像処理装置２０２´が求めた左右のそれぞれのＳＣＭ＿ｉｎｄｅｘ時系列データの振幅と周期を求め（ステップＳ８０１）、ＳＣＭ＿ｉｎｄｅｘの振幅が一定以上であり、かつ、この周期が呼吸周期であるか否かを判定する（ステップＳ８０２）。例えば、通常呼吸周期は４秒程度であり、頻呼吸の場合でも、その周期は２秒程度であるので、一定区間のＳＣＭ＿ｉｎｄｅｘの時系列データの周期性の解析に自己相関解析を行い、求めた周期が１〜６秒 に入るかどうかで呼吸周期であるか否かの判定できる。ステップＳ８０２の判定結果がＹｅｓの場合、左右のＳＣＭ＿ｉｎｄｅｘ時系列データの相関係数ｒを求める（ステップＳ８０３）。これが、予め決められた一定値以上か否かを判定する（ステップＳ８０４）。例えば、ｒが０．７以上であれば、呼吸補助筋が動員されたと判定できる。そしてステップＳ８０４の判定結果がＹｅｓの場合は呼吸補助筋の動員がされているものと判断する（ステップＳ８０７）。
【００５７】
第一の実施形態の処理と異なるのは、ステップＳ８０２またはステップＳ８０４の判定結果がＮｏの場合、すなわちＳＣＭ＿ｉｎｄｅｘの振幅が一定以下、または、ＳＣＭ＿ｉｎｄｅｘの周期が呼吸周期と判定されなかった場合、あるいは、左右のＳＣＭ＿ｉｎｄｅｘの相関係数が一定値以下だった場合である。この場合は、人体の頸部を含む動画像から呼吸補助筋の動員を判定することができないため、肩の動きから呼吸努力の有無を判定する。
【００５８】
動画像処理装置２０２´で求めた、ある時刻の肩峰４０３と首付根位置４０２との上下方向の相対距離４０９を、Ｄｓ（ｔ）とすると、肩の動きが呼吸補助筋の動員によるものか否かは、ＳＣＭ＿ｉｎｄｅｘ同様、この相対距離Ｄｓ（ｔ）が周期的に変動し、かつ、その周期が呼吸周期と同等であるか否かを調べること、および、努力呼吸時は肩の動きが一定以上の大きさであることから、その振幅が一定値以上かどうかを調べることで、判定できる。
【００５９】
まず、呼吸動作判定装置２０３´は、Ｄｓ（ｔ）の時系列データの計測時間での平均振幅（ΔＤ  ｍｅａｎ）と周期を算出する（ステップＳ９０１）。具体的には、第一実施例と同様、この周期と振幅は、波形の周波数解析、ピーク検出等周知の方法により求められる。
【００６０】
次に、呼吸動作判定装置２０３´は、この振幅が予め決められた一定値（Ｄ０）以上の大きさで、かつ、周期が呼吸周期なのかどうかを判定する（ステップＳ９０２）。呼吸周期の判定方法はステップＳ８０２と同様、自己相関解析から周期を求めて呼吸周期と比較する方法等がある。
【００６１】
ステップＳ９０２の判定結果がＮｏの場合は、ステップＳ９０３に進み、呼吸補助筋の動員はないものと判断する。ステップＳ９０２の判定結果がＹｅｓの場合は呼吸補助筋の動員がされているものと判断する（ステップＳ８０７）。この場合は、好ましくは、呼吸評価演算装置２０４´において呼吸補助筋の動員度合いを評価する。第二の実施例では、以下の式３のような処理により、努力呼吸指標を求める。
【００６２】
【数３】
2020226182.tif 000005
【００６３】
ここで、ａはＳＣＭ＿ｉｎｄｅｘ、すなわち、呼吸補助筋から求める努力呼吸指標とレベルを合わせるため、指標を正規化するための係数である。予め集めた人のデータを統計的に解析して求めた値を用いることができる。本実施形態によれば、呼吸専門医が努力呼吸を視診で判断する際の呼吸補助筋の動員を画像センサーと情報処理技術で指標化し、誰もが簡単に判断できるようになるため、在宅等で非医療者が努力呼吸に特徴的な呼吸補助筋の動員を簡便に計測かつ客観的に表現でき、在宅慢性疾患患者の病態をより正確かつ効率的に医療者に伝えることが可能になる。
【００６４】
［第三の実施形態に係る呼吸状態検出装置］
上記、第一および第二の実施形態では、算出した努力呼吸指標の利用についての機能は含まれていない。しかし、実際の利用場面を考えると、この努力呼吸指標を患者ごとに記録し、そのトレンドの変化から良い状態が続いているのか、または、増悪を起こしそうなのか等の目安が示されると良い。本実施形態は、そのような目的のために、算出した努力呼吸指標を経時的に記録する記録装置である結果記録装置（１００１）をさらに加えたものである（図１１）。このような、日々の努力呼吸指標を、例えば、数値、または／および、グラフとしてディスプレイ５００等に表示すること、および／または、過去、結果記録装置に記録されたデータと呼吸評価演算装置の出力した評価結果、すなわち、現在の努力呼吸指標とを比較して表示することができれば、患者の状態をより客観的にかつ直感的にわかりやすく判断することが可能となる。
【００６５】
以上の説明において、被験者の努力呼吸の有無の判定を図３及び図１１に示す呼吸状態検出装置２００（２００´）で行う例を示したが、このような例には限られない。即ち、呼吸状態検出装置２００（２００´）が観測した被験者のデータを通信Ｉ／Ｆ３によりサーバに送信し、サーバ側が被験者の努力呼吸の有無の判定するようにしてもよい。
【００６６】
以上説明したように、本発明によれば、在宅等で非医療者が努力呼吸に特徴的な呼吸補助筋の動員を簡便に計測かつ客観的に表現できるようになるため、在宅における慢性呼吸疾患患者の病態をより正確かつ効率的に医療者に伝えることが可能になる。
【図１】
2020226182.tif 000006
【図２】
2020226182.tif 000007
【図３】
2020226182.tif 000008
【図４】
2020226182.tif 000009
【図５】
2020226182.tif 000010
【図６】
2020226182.tif 000011
【図７】
2020226182.tif 000012
【図８】
2020226182.tif 000013
【図９】
2020226182.tif 000014
【図１０】
2020226182.tif 000015
【図１１】
2020226182.tif 000016
【国際調査報告】
2020226182.tif 000017
2020226182.tif 000018
2020226182.tif 000019
2020226182.tif 000020