WO2022220209A1 - 脂質濃度計測装置、プログラム、及び、方法 - Google Patents

Abstract

被検体へ光を照射する照射部と、前記照射部からの光が被検体内の脂質粒子により散乱し、前記被検体から放出された光を受光して、受光強度を検出する受光部と、前記被検体の心拍による脈動に応じた受光強度の時間変化による波形から、1以上の0.5Hz以上2.0Hz以下の周期の1周期波形の領域を判定し、前記判定された1周期波形の領域の波形面積を算出し、前記波形面積から脂質濃度を算出する制御部と、を有する。

Description

脂質濃度計測装置、プログラム、及び、方法

　この明細書に記載の実施形態は、脂質濃度計測装置、プログラム、及び、方法に関する。

　従来の血流データに基づく脂質濃度計測方法は、レーザー光など位相の揃った光源を用いることが必須であった。（例えば、特許文献１参照）。

国際公開第2018/151022号

　脂質濃度の絶対値計測は、測定結果の信頼性を得るためには重要な情報である。しかしながら、脂質濃度の絶対値を非侵襲的に計測する場合において、肌の色や血管の深さや皮膚の厚みなどが脂質濃度の絶対値計測における阻害要因となっていた。

　また、これらの阻害要因の情報を計測するために、様々な計測を組み合わせることも想定されるが、そのような場合、装置の大型化や高額化を避けることは難しくなり、簡易計測を実現することは難しい。

　さらに、特許文献１は、静脈の計測技術である。静脈には心拍による脈動（拍動）は観察されないことが知られている。したがって、特許文献１の装置は、心拍の周期性（例えば、0.5～2.0Hz）に基づく計測を行う装置ではない。

　本発明は、このような従来の課題を解決するためになされた発明であって、小型で簡易な計測装置で脂質濃度の絶対値計測を可能とする装置、プログラム、及び、装置の作動方法を提供するものである。

　本発明の脂質濃度計測装置は、被検体へ光を照射する照射部と、照射部からの光が被検体内の脂質粒子により散乱し、被検体から放出された光を受光して、受光強度を検出する受光部と、被検体の心拍による脈動に応じた受光強度の時間変化による波形から、波形の内の複数の１周期波形の領域のうち、0.5Hz以上2.0Hz以下の周期の１周期波形の領域を判定し、判定された１周期波形の領域の波形面積を算出し、波形面積から脂質濃度を算出する、制御部と、を有する。

　本発明の脂質濃度計測プログラムは、コンピュータに、被検体へ光を照射する処理と、照射部からの光が被検体内の脂質粒子により散乱し、被検体から放出された光を受光して、受光強度を検出する処理と、被検体の心拍による脈動に応じた受光強度の時間変化による波形から、波形の内の複数の１周期波形の領域のうち、0.5Hz以上2.0Hz以下の周期の１周期波形の領域を判定し、判定された１周期波形の領域の波形面積を算出し、波形面積から脂質濃度を算出する処理とを実行させる。

　本発明の脂質濃度計測方法は、被検体へ光を照射する処理と、照射部からの光が被検体内の脂質粒子により散乱し、被検体から放出された光を受光して、受光強度を検出する処理と、被検体の心拍による脈動に応じた受光強度の時間変化による波形から、波形の内の複数の１周期波形の領域のうち、0.5Hz以上2.0Hz以下の周期の１周期波形の領域を判定し、判定された１周期波形の領域の波形面積を算出し、波形面積から脂質濃度を算出する処理と、を含む。

脂質濃度計測装置の概略図 脂質濃度計測装置のブロック図 受光強度の時間変化の測定結結果および積分範囲の模式図 脂肪負荷前後の波形比較の図 脂質濃度計測装置による測定結果と脂質濃度変化量との相関図 脂質濃度が低いときの模式図 脂質濃度が高いときの模式図 実施形態の脂質濃度計測ジョブのフローチャート

　以下に実施形態を図面を用いて説明する。

　図１は、実施形態の脂質濃度計測装置１の構成例を概略的に示す図である。図１に示すように脂質濃度計測装置１は、照射部２、受光部３、及び、制御部４を有する。また、照射部２による生体（被検体）上の照射位置を照射位置２１とし、受光部３による生体上の受光位置を受光位置３１とする。

　照射部２は照射光を生体に照射する。照射部２は、制御部４により、照射する光の波長や照射強度が制御されてもよい。実施形態では、照射部２はLED（Light Emitting Diode）（810nm）である。

　実施形態の照射部２は、光の連続的な照射や光のパルス状の照射等の光を照射する時間の長さを任意に調整することができる。

　照射部２は、波長が固定された光源を用いてもよい。照射部２は、波長が異なる複数の光源あるいは複数の波長の光を混合したものであってもよい。照射部２は、例えば、蛍光灯、LED、レーザー、白熱灯、HID、ハロゲンランプ等である。照射部２の照度は、制御部４により制御されてもよい。

　実施形態の受光部３は、受光位置３１において生体内から生体外に放出される光を受光して、光強度を検出する。実施形態の受光部３は、フォトダイオードである。受光部３は、フォトダイオードに限られず、CCDやCMOSでもよい。受光部３は、波長を可入射波長に設定し、その波長を受光できるものでもよい。実施形態では、受光部３にフォトダイオードを用い、サンプリングレートは２ｍｓに設定した。

　次に、脂質濃度計測装置１の制御系の構成について説明する。図２は実施形態の脂質濃度計測装置１のブロック図である。システムバス１４２を介して、ＣＰＵ（Central Processing Unit）１４１、ＲＯＭ（Read Only Memory）１４３、ＲＡＭ（Random Access Memory）１４４、記憶部１４５、外部Ｉ／Ｆ（Interface）１４６、照射部２、及び、受光部３が接続される。ＣＰＵ１４１とＲＯＭ１４３とＲＡＭ１４４とで制御部４を構成する。

　ＲＯＭ１４３は、ＣＰＵ１４１により実行されるプログラムや閾値を予め記憶する。

　ＲＡＭ１４４は、ＣＰＵ１４１が実行するプログラムを展開するエリアと、プログラムによるデータ処理の作業領域となるワークエリアなどの様々なメモリエリア等を有する。ＲＡＭ１４４は、図５に示すような、事前に生体を用いて計測した波形領域の波形面積データと脂質濃度との関係を示す校正データを記憶する記憶エリアを有する。

　図３は、健常人の空腹時の測定結果である。

　図3は、脂質濃度計測器を用いて、被験者の指先で脈動を計測した結果である。図３の縦軸は、計測器で計測された光の吸収強度であり、横軸は時間である。図３に示すように、約1秒に1回（すなわち、約1Hz）の周期性を持つ波形が得られた。これは、心拍による脈動を反映している。また、図中に示した部分（以下、1周期波形領域とする）の面積は、積分分析の積分範囲を示すものである。1周期波形領域の波形面積は、以下の数式１により算出できる。

　なお、実施形態では、面積として、1周期波形領域（１回の周期（１つのピーク））の波形面積を求めているが、複数の周期の波形領域の平均の面積を波形面積としてもよい。また、波形領域の数をn=3以上とり、中央値をとったり、最大値と最小値を除いて平均する、SDを求めSD範囲から外れたものを除外してから平均するなどして波形面積を求めても良い。

　図３に示すように、波形面積を求める1周期波形領域は、ベースライン（図３のａ）を差し引いたのち、山なりの1周期の波形が示すラインとベースラインで包囲した領域（図３のｂ）である。脂質濃度自体でもベースラインが上がるが、測定部位などによる影響でもベースラインの高さが変わってしまうため、規格化としてベースラインを設定するのがよい。

　また、ベースラインａは1周期波形領域のピーク波形の山の起点と終点を直線で結び、Ｘ軸と平行とするのがよい。ピークの定義は、心拍による変化であり、約１Hzの周期となる。生体計測で想定される測定レンジとして徐脈（心拍数はおおよそ40）から激しい運動時（心拍数はおおよそ120）までを考慮すると、心拍の脈動の変化に伴う、0.5Hz以上2.0Hz以下の周期の波が本計測における1周期波形の領域であると判定できる。このベースラインａの起点の受光強度と終点の受光強度は同じ値になることが望ましいが、測定条件によって、必ずしも起点の受光強度と終点の受光強度とが同じ値の形状の1周期波形領域のピークが得られるとは限らない。したがって、複数の１周期波形領域を計測し、起点の受光強度と終点の受光強度の比が最も1に近い１周期波形領域を判定し、当該判定された１周期波形領域について波形面積を算出することでもよい。さらには、起点の受光強度と終点の受光強度の比が0.8～1.2となる複数の１周期波形領域において面積を求め、複数の面積の平均を求めて波形面積とし、この波形面積から脂質濃度を求めることでもよい。

　受光強度の時間変化を測定した結果、空腹時、脂肪負荷後とも、約1秒に1回の受光強度の脈動が確認できる。このときの受光強度は、受光により得られた電圧値である。なお、ここで言う電圧値は厳格なmv単位ではないので、a.u.単位で表記した。

　図４は、脂肪負荷前（空腹時）と脂肪を負荷した際の脈動の波形である。脂肪負荷前（空腹時）については、図３中の領域ｂに示す1周期波形領域（初回1回の脈動）を拡大したものである。図4に示したように、脂肪負荷前（空腹時）と脂肪負荷後（脂肪負荷後）で、1周期波形領域の波形は異なる。

　図５は、図４で得られた1周期波形領域の数式１による積分値（波形面積）と脂質濃度の絶対値の相関を示した校正データ（線形近似データ）である。図に示したように、脂質濃度の上昇に伴い、1周期波形領域の波形面積の拡大が確認できる。

　図６と図７は、血中の脂質粒子が多い場合と少ない場合の模式図である。図６が脂質粒子が少なく、図７が脂質粒子が多い場合である。脂質粒子が増加すると、血液内で光の反拡散回数が増えることで、光路長が増すことになる。そのため、血液による吸収が増大し、図5のような正相関になったと考えられる。

　周期波形は、脈動に伴う血液の増減を反映しており、心臓の拍動により血液が拍出された際の、測定部位における血液の一時的な増加によるものである。そして光計測では、血液による光の吸収を計測しており、血液量が増加すると、吸収が増加する。そして、脈動に応じ周期波形（１周期波形など）として、観測される。このとき、血液内部に散乱体である脂質粒子が増加すると、血液内部で反射する回数が増えるため、結果的に光路長が長くなり、血液による吸収の影響を受けやすくなる。そのため、周期波形領域の面積は、脂質粒子の増加に伴い大きくなる。

　記憶部１４５は、処理に必要なデータを記憶する。記憶部１４５は、例えば、HDD(Hard Disk Drive)などである。記憶部１４５は、図５に示すような、事前に生体を用いて計測した波形領域の波形面積データと脂質濃度との関係を示す校正データを記憶してもよい。

　外部Ｉ／Ｆ１４６は、例えばクライアント端末（ＰＣ）などの外部装置と通信するためのインターフェースである。外部Ｉ／Ｆ１４６は、外部装置とデータ通信を行うインターフェースであれば良く、たとえば、外部装置にローカルに接続する機器（ＵＳＢメモリ等）であっても良いし、ネットワークを介して通信するためのネットワークインターフェースであっても良い。

　以上のような構成を備える脂質濃度計測装置１において、予め設定されているプログラムに基づいて、脂質濃度計測装置１は脂質濃度計測ジョブを実行する。図８は、実施形態の脂質濃度計測ジョブのフローチャートである。

　実施形態の脂質濃度計測装置は、上記構成を有する照射部２により、被検体へ光を照射する（ＳＴＥＰ１０１）。受光部３が、被検体から放出される受光強度を検出する（ＳＴＥＰ１０２）。制御部４が、被検体の心拍による脈動に応じた受光強度の時間変化による波形から、波形の内の0.5Hz以上2.0Hz以下の周期の１以上の１周期波形の領域を判定し、判定された１周期波形の所定の波形領域を上記数式１により積分し、波形領域の波形面積を算出する（ＳＴＥＰ１０３）。

　制御部４は、複数の１周期波形の領域の面積を平均して前記波形面積を算出してもよい。制御部４は、複数の１周期波形の領域のうち、当該１周期波形の起点の受光強度と終点の受光強度の比が最も1に近い１周期波形の領域を判定し、判定された１周期波形の領域の波形面積を算出してもよい。制御部４は、複数の前記１周期波形の領域のうち、当該１周期波形の起点の受光強度と終点の受光強度の比が0.8以上1.2以下となる、１以上の１周期波形の領域を判定し、判定された１以上の１周期波形の領域の面積を平均して波形面積を算出してもよい。

　制御部４は、0.5Hz以上2.0Hz以下の周期の1周期波形領域の波形面積から脂質濃度を算出する（ＳＴＥＰ１０４）。制御部４は、ＳＴＥＰ１０３で算出された1周期波形領域の波形面積と、事前に計測した1周期波形領域の波形面積とＴＧ濃度の絶対値との関係を示す校正データと、から、ＴＧ濃度（脂質濃度）の絶対値を算出する。

　次に、実施形態の脂質濃度計測プログラムについて説明する。なお、当該プログラムは記憶媒体に格納されてもよい。

　実施形態の脂質濃度計測プログラムは、装置のコンピュータに、被検体へ光を照射する処理と、被検体から放出される受光強度を検出する処理と、被検体の心拍による脈動に応じた受光強度の時間変化による波形から、波形の内の0.5Hz以上2.0Hz以下の周期の１以上の１周期波形の領域を判定し、判定された１周期波形の所定の波形領域を上記数式１により積分し、波形領域の波形面積を算出する処理と、1周期波形領域の波形面積から脂質濃度の絶対値を算出する処理とを実行させる。

　実施形態の脂質濃度計測プログラムは、装置のコンピュータに、複数の１周期波形の領域の面積を平均して前記波形面積を算出する処理を実行させてもよい。実施形態の脂質濃度計測プログラムは、装置のコンピュータに、複数の１周期波形の領域のうち、当該１周期波形の起点の受光強度と終点の受光強度の比が最も1に近い１周期波形の領域を判定し、判定された１周期波形の領域の波形面積を算出する処理を実行させてもよい。実施形態の脂質濃度計測プログラムは、装置のコンピュータに、複数の前記１周期波形の領域のうち、当該１周期波形の起点の受光強度と終点の受光強度の比が0.8以上1.2以下となる、１以上の１周期波形の領域を判定し、判定された１以上の１周期波形の領域の面積を平均して波形面積を算出する処理を実行させてもよい。

　なお、実施形態では、照射部と受光部と制御部を一体の装置として構成したが、これに限られず、例えば、照射部として、携帯端末（スマートフォン、タブレット、モバイルＰＣ）などのユーザ装置に備わる光源を使用し、受光部として携帯端末（スマートフォン、タブレット、モバイルＰＣ）などのユーザ装置に備わるセンサ（ＣＭＯＳ等）を使用し、制御部をユーザ装置とネットワーク接続したサーバー装置等に設置してもよい。

　実施形態の脂質濃度計測装置は、被検体へ、光を照射する照射部と、被検体から放出される受光強度を検出する受光部を有するユーザ装置に、通信可能に接続する。脂質濃度計測装置は、ユーザ装置から送信された受光強度の時間変化による波形を上記数式１により積分して1周期波形領域の波形面積を算出し、1周期波形領域の波形面積から脂質濃度の絶対値を算出する制御部を有する。

　以上、実施形態を説明したが、この実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。この新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。この実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。

１：脂質濃度計測装置
２：照射部
３：受光部
４：制御部

Claims (8)

1. 　被検体へ光を照射する照射部と、
   　前記照射部からの光が前記被検体内の脂質粒子により散乱し、前記被検体から放出された光を受光して、受光強度を検出する受光部と、
   　前記被検体の心拍による脈動に応じた前記受光強度の時間変化による波形から、1以上の0.5Hz以上2.0Hz以下の周期の１周期波形の領域を判定し、
   　前記判定された１周期波形の領域の波形面積を算出し、
   　前記波形面積から脂質濃度を算出する、制御部と、
   を有する脂質濃度計測装置。
2. 　前記１周期波形の領域は、前記１周期波形と、前記１周期波形の起点の受光強度と終点の受光強度を直線で接続したベースラインと、により包囲された領域である、請求項１に記載の脂質濃度計測装置。
3. 　前記制御部は、
   　複数の前記１周期波形の領域の面積を平均して前記波形面積を算出する、
   請求項１に記載の脂質濃度計測装置。
4. 　前記制御部は、
   　複数の前記１周期波形の領域のうち、前記１周期波形の起点の受光強度と終点の受光強度の比が最も1に近い１周期波形の領域を判定し、前記判定された１周期波形の領域の波形面積を算出する、
   請求項１に記載の脂質濃度計測装置。
5. 　前記制御部は、
   　複数の前記１周期波形の領域のうち、前記１周期波形の起点の受光強度と終点の受光強度の比が0.8以上1.2以下となる１周期波形の領域を判定し、前記判定された１以上の１周期波形の領域の面積を平均して波形面積を算出する、請求項１に記載の脂質濃度計測装置。
6. 　前記制御部は、
   　事前に計測した前記波形面積と前記脂質濃度との関係を示す校正データを記憶し、前記校正データから、前記脂質濃度を算出する、請求項１から５のいずれかに記載の脂質濃度計測装置。
7. 　コンピュータに、
   　被検体へ光を照射する処理と、
   　前記照射部からの光が前記被検体内の脂質粒子により散乱し、前記被検体から放出された光を受光して、受光強度を検出する処理と、
   　前記被検体の心拍による脈動に応じた前記受光強度の時間変化による波形から、1以上の0.5Hz以上2.0Hz以下の周期の１周期波形の領域を判定し、
   　前記判定された１周期波形の領域の波形面積を算出し、
   　前記波形面積から脂質濃度を算出する処理と、
   を実行させる脂質濃度計測プログラム。
8. 　被検体へ光を照射する処理と、
   　前記照射部からの光が前記被検体内の脂質粒子により散乱し、前記被検体から放出された光を受光して、受光強度を検出する処理と、
   　前記被検体の心拍による脈動に応じた前記受光強度の時間変化による波形から、1以上の0.5Hz以上2.0Hz以下の周期の１周期波形の領域を判定し、
   　前記判定された１周期波形の領域の波形面積を算出し、
   　前記波形面積から脂質濃度を算出する処理と、
   を含む脂質濃度計測方法。

Images