JPWO2018143119A1 - 脂質計測装置及びその方法 - Google Patents

脂質計測装置及びその方法 Download PDF

Info

Publication number
JPWO2018143119A1
JPWO2018143119A1 JP2018565525A JP2018565525A JPWO2018143119A1 JP WO2018143119 A1 JPWO2018143119 A1 JP WO2018143119A1 JP 2018565525 A JP2018565525 A JP 2018565525A JP 2018565525 A JP2018565525 A JP 2018565525A JP WO2018143119 A1 JPWO2018143119 A1 JP WO2018143119A1
Authority
JP
Japan
Prior art keywords
light
reachable range
lipid
living body
light reachable
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.)
Granted
Application number
JP2018565525A
Other languages
English (en)
Other versions
JP6894087B2 (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.)
Medical Photonics Co Ltd
Original Assignee
Medical Photonics Co Ltd
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 Medical Photonics Co Ltd filed Critical Medical Photonics Co Ltd
Publication of JPWO2018143119A1 publication Critical patent/JPWO2018143119A1/ja
Application granted granted Critical
Publication of JP6894087B2 publication Critical patent/JP6894087B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4869—Determining body composition
    • A61B5/4872—Body fat
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B10/00—Other methods or instruments for diagnosis, e.g. instruments for taking a cell sample, for biopsy, for vaccination diagnosis; Sex determination; Ovulation-period determination; Throat striking implements
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/14546—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring analytes not otherwise provided for, e.g. ions, cytochromes
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4866—Evaluating metabolism
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47—Scattering, i.e. diffuse reflection
    • G01N21/4795—Scattering, i.e. diffuse reflection spatially resolved investigating of object in scattering medium
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • G16H40/63—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for local operation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Optics & Photonics (AREA)
  • Obesity (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Databases & Information Systems (AREA)
  • Business, Economics & Management (AREA)
  • General Business, Economics & Management (AREA)
  • Data Mining & Analysis (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

【課題】 計測者の熟練がなくとも、容易に非侵襲脂質計測が可能な装置及び方法を提供する。【解決手段】 生体外から生体内に向けて、生体の所定の部位に所定の光強度で光を照射する照射部と、生体から放出される光強度に基づいて生体内の光到達範囲を検出する光強度検出部と、光到達範囲に基づいて所定の光到達範囲パラメータを算出し、光到達範囲パラメータに基づいて生体内の脂質濃度を算出する制御部と、を有する。【選択図】 図１

Description

本発明は、脂質計測装置及びその方法に関する。
食後高脂血症は、動脈硬化のリスクファクターとして注目されている。非空腹時の中性脂肪濃度が高くなると、冠動脈疾患のイベント発症リスクが高くなることが報告されている。
食後高脂血症の診断は、食後6〜8時間の血中の脂質濃度変化を観測する必要がある。つまり、食後の高脂血状態を計測するためには、被験者を6〜8時間拘束し、複数回の採血が必要である。そのため、食後高脂血症の診断は臨床研究の域を出ず、食後高脂血症の診断を臨床現場で実施することは、現実的ではなかった。
このような課題を解決する手法が、特許文献１に開示されている。特許文献１の手法によれば、非侵襲脂質計測により、採血を無くすことができる。これにより医療機関のみならず家庭でも血中脂質を計測できるようになる。即時的なデータ取得を可能とすることで、時間的に連続した血中脂質を計測することが可能となる。
国際公開第2014/087825号公報
しかしながら、特許文献１に示すような非侵襲の脂質計測手法では、最適計測部位の決定に計測者の熟練を要し、計測誤差の要因となっていた。
また、生体を透過する光は、皮膚、筋肉、血液などの影響により、その光強度が減衰する。生体中の特定の物質濃度を検出するためには、測定対象物以外の影響を極力取り除くことが望ましい。
一方で、精度はS/N（signal/noise）比で表現されるように、測定対象物質のシグナル強度をより強く検出することで、測定精度向上も可能といえる。
特許文献１に示す計測手法は、１次元（ライン状）の検出であったが、静脈、筋、骨などにより光拡散が不均一であるため、計測時における測定機器の位置ズレや脱着などによる同一場所での計測が困難である。したがって、精度よく計測するためには、計測者の熟練が必要であった。
本発明は、このような従来の問題を解決するためになされたもので、計測者の熟練がなくとも、容易に非侵襲脂質計測を可能とする装置及び方法を提供することである。
本発明の脂質計測装置は、生体外から生体内に向けて、生体の所定の部位に所定の光強度で光を照射する照射部と、生体から放出される光強度に基づいて生体内の光到達範囲を検出する光強度検出部と、光到達範囲に基づいて所定の光到達範囲パラメータを算出し、光到達範囲パラメータに基づいて生体内の脂質濃度を算出する制御部とを有する。
また、本発明の脂質計測方法は、生体外から生体内に向けて、生体の所定の部位に所定の光強度で光を照射する照射工程と、生体から放出される光強度に基づいて生体内の光到達範囲を検出する光強度検出工程と、光到達範囲に基づいて所定の光到達範囲パラメータを算出するパラメータ算出工程と、光到達範囲パラメータに基づいて生体内の脂質濃度を算出する脂質濃度算出工程と、を有する。
また、本発明の脂質計測装置は、生体外から生体内に向けて、生体の所定の部位に所定の光強度で光を照射する照射部と、生体から放出される光強度に基づいて生体内の光到達範囲を検出する光強度検出部と、光強度検出部により検出された光到達範囲を送信する通信部とを有するユーザー装置に、通信可能に接続される脂質計測装置であって、ユーザー装置から送信された光到達範囲に基づいて所定の光到達範囲パラメータを算出し、光到達範囲パラメータに基づいて生体内の脂質濃度を算出する制御部を有する。
本発明の脂質計測装置及び方法によれば、計測者の熟練がなくとも、容易に非侵襲脂質計測が可能となる。
実施形態の脂質計測装置の構成を示す図である。 血中脂質による光の散乱を示す図である。 実施形態の脂質計測装置の制御系の構成を示す図である。 円状の光到達範囲を示す図である。 歪んだ形状の光到達範囲を示す図である。 実施形態の脂質計測装置の作動方法のフローチャートである。 実施形態の脂質計測システムの構成を示す図である。 実施形態の脂質計測装置の制御系の構成を示す図である 光到達範囲を撮像した結果を示す図である。 光到達範囲を撮像した結果を示す図である。 静脈近傍の光到達範囲を撮像した結果を示す図である。 脂質濃度の変化と光到達範囲の面積の変化とを比較した図である。 脂質濃度の変化と光到達範囲の面積の変化とを比較した図である。 光の最小到達距離と脂質濃度との間の関連性を示す図である。 光の最小到達距離と脂質濃度との間の関連性を示す図である。 光の到達体積と脂質濃度との間の関連性を示す図である。 光の到達体積と脂質濃度との間の関連性を示す図である。 図２と異なる照射部と光強度検出部の配置を示す図である。 図１５に示す照射部と光強度検出部の配置により撮影した例を示す図である。 図１５に示す照射部と光強度検出部の配置により計測した結果を示す図である。 図１５に示す照射部と光強度検出部の配置により計測した結果を示す図である。
以下、実施形態である脂質計測装置及びその作動方法について、図を参照して詳細に説明をする。
図１は、実施形態の脂質計測装置の構成を示す図である。
図１に示すように、実施形態の脂質計測装置１００は、生体の外から生体内に向けて生体の所定の部位に光を照射する照射部１０１と、生体から放出される光を受光して、光強度に基づいて生体内の光到達範囲Ｆを検出する光強度検出部１０２と、光強度検出部１０２により検出された光到達範囲Ｆに基づき光到達範囲パラメータを算出するパラメータ算出し、光到達範囲パラメータに基づき脂質濃度を算出する制御部１０３とを有する。
照射部１０１は、生体の所定の部位の生体外から生体内に向けて、所定の照射位置に光を照射するための光源を有する。実施形態の照射部１０１は、照射する光の波長を調整することができる。照射部１０１は、波長範囲を血漿の無機物によって光が吸収される波長範囲以外に調整できる。照射部１０１は、血液の細胞成分によって光が吸収される波長範囲以外に調整できる。ここで、血液の細胞成分とは、血中の赤血球、白血球及び血小板である。血漿の無機物とは、血中の水及び電解質である。
照射部１０１が照射する光の波長範囲は、血漿の無機物により光を吸収する波長範囲を考慮して約1400nm以下、及び、約1500nm〜約1860nmとするのが好ましい。さらに、照射部１０１が照射する光の波長範囲は、血液の細胞成分によって光が吸収される波長範囲を考慮して約580nm〜約1400nm、及び、約1500nm〜約1860nmとするのがより好ましい。
照射部１０１に用いられる波長範囲を上記範囲とすることにより、後述する光強度検出部１０２により検出される光において、血漿の無機物による光の吸収の影響、及び、血液の細胞成分により光の吸収の影響を抑制する。これにより、物質を特定するほどの吸収は存在せず、吸収による光エネルギー損失は無視できるほど小さくなる。そのため、血中の光は血中の脂質による散乱によって遠くまで伝搬し、体外へ放出される。
実施形態の照射部１０１は、光の連続的な照射や光のパルス状の照射等の光を照射する時間長さを任意に調整することができる。照射部１０１は、照射する光の強度または光の位相を任意に変調することができる。
照射部１０１は、波長が固定された光源を用いてもよい。照射部１０１は、波長が異なる複数の光源あるいは複数の波長の光を混合したものであってもよい。
光強度検出部１０２は、生体から生体外に放出される光を受光して、その光強度を検出し、生体内の光到達範囲Ｆを検出する。
図２は、血中脂質による光の散乱を示す図である。図２に示すように、照射部１０１から、生体Ｄ表面の照射位置（図中のＥ）に照射された光（図中のＢ）は、リポ蛋白等の脂質が存在する深さに達したのち、生体Ｄ内の血液中の脂質（図中のＡ）によって反射する。さらに、照射された光は血中の脂質による光の散乱を経て、後方散乱光（図中のＣ）が生体から放出される。光強度検出部１０２は後方散乱光Ｃの光強度を検出する。
なお、図２では、照射部１０１の先端は生体Ｄに接触しているが、図１３に示すように照射部１０１の先端は生体Ｄから所定の距離、離れていてもよい。
図２に示すように、照射部１０１の照射位置Ｅから、所定レベルの光強度が得られる範囲（以下、光到達範囲Ｆとする）の外周までの距離を光到達距離ｌとする。
計測対象であるリポ蛋白は、アポ蛋白等に覆われた球状構造をしている。リポ蛋白は血中において固体のような状態で存在する。リポ蛋白は、光を反射する性質を有する。特に、粒子径や比重の大きいカイロミクロン（CM）やVLDL等は中性脂肪（TG）を多く含み、光をより散乱させ易い特性を有する。よって、光強度検出部１０２により検出される光強度には、リポ蛋白による光の散乱の影響が含まれる。
光強度検出部１０２は、ＣＣＤやＣＭＯＳ等の受光素子であってもよい。また、光強度検出部１０２は、受光素子をアレイ状に配置したものでもよく、同心円状に配置してもよい。受光素子数を少なくする場合には、受光素子を照射位置Ｅを中心に十字状、V字状に配置してもよく、直線上に配置し移動や回転をさせて計測してもよい。
また、図２では、光強度検出部１０２が、照射部１０１の真上に配置されるが、これに限られず、光到達範囲Ｆを検出できる位置であればよい。
次に、脂質計測装置１００の制御系の構成について説明する。図３は実施形態の脂質計測装置１００のブロック図である。システムバス１０９を介して、ＣＰＵ（Central Processing Unit）１０４、ＲＯＭ（Read Only Memory）１０５、ＲＡＭ（Random Access Memory）１０６、記憶部１０７、外部Ｉ／Ｆ（Interface）１０８、照射部１０２、及び、光強度検出部１０２が接続される。ＣＰＵ１０４とＲＯＭ１０５とＲＡＭ１０６とで制御部（コントローラー）１０３を構成する。
ＲＯＭ１０５は、ＣＰＵ１０４により実行されるプログラムや閾値を予め記憶する。
ＲＡＭ１０６は、ＣＰＵ１０４が実行するプログラムを展開するエリアと、プログラムによるデータ処理の作業領域となるワークエリアなどの様々なメモリエリア等を有する。
記憶部１０７は、予め用意された、静的パラメータ及び動的パラメータの適切な数値範囲のデータを記憶する。記憶部１０７は、ＨＤＤ（Hard Disk Drive）や、フラッシュメモリや、ＳＳＤ（Solid State Drive）等の、不揮発性に記憶する内部メモリーでよい。
外部Ｉ／Ｆ１０８は、例えばクライアント端末（ＰＣ）などの外部装置と通信するためのインターフェースである。外部Ｉ／Ｆ１０８は、外部装置とデータ通信を行うインターフェースであれば良く、たとえば、外部装置にローカルに接続する機器（ＵＳＢメモリ等）であっても良いし、ネットワークを介して通信するためのネットワークインターフェイスであっても良い。
制御部１０３は、光強度検出部１０２により検出された光到達範囲Ｆに基づき、光到達範囲パラメータを算出する。
光到達範囲Ｆの検出は、２値化法を採用すればよい。光強度検出部１０２が検出した光強度を0〜255の256段階とし、光強度検出部１０２は、光強度のしきい値を254に設定し、255の場合の光到達範囲Ｆとする。
また、照射部１０１から光強度検出部１０２が遠くなるほうが、散乱の影響をよく反映するため、上記のしきい値に限定されず、値を下げてもよい。この場合、実際の計測においては外乱光の影響を受けやすくなるため、装置形状、遮光の程度、受光素子の感度により適時設定を行うことが好ましい。
また受光部にPDなどの素子を用いる場合、AD値や電圧値を用いてもよく、計測に用いる測定範囲は、照射強度、受光素子の感度、遮光の程度により適時設定を行うことが好ましい。
本実施例では、暗室にて外乱光が無視できる状態で行った場合は、強度が１０程度のノイズが入り込むことから、強度１１以上の場合についても検討した。
図４は、図２のＸ方向から見た、生体表面の光到達範囲Ｆを示す図である。図に示すように、毛細管のみの場合には、照射光は照射位置Ｅを中心とし光到達距離ｌを半径とした円状に拡散し、光到達範囲Ｆは生体表面で円状となる。
制御部１０３は、光到達範囲パラメータとして、光到達範囲Ｆにおける照射位置Ｅから光到達範囲の外周（外縁）までの距離（光到達距離ｌとする）を算出する。
また、制御部１０３は、光到達範囲パラメータとして、光到達範囲Ｆの面積（光到達面積Ｓとする）を算出する。なお、光到達面積Ｓは光到達距離ｌから算出してもよい。光到達面積Ｓは、しきい値255の画素数から算出してもよい。計測誤差を平均化するため、光到達面積Ｓを光到達最大距離と光到達最小距離から楕円面積として算出してもよい。
また、制御部１０３は、光到達範囲パラメータとして、光到達範囲Ｆの体積（光到達体積Ｖとする）を算出する。なお、光到達体積Ｖは、光到達体積V＝（4/3π×a×b×c）/2と、計算することもできる。
式中のa,b,cは、それぞれ球体の座標軸x,y,z方向のそれぞれが９０°に交わる半径である。また、光到達範囲に歪みがない場合は、a=b=cとなるため、図２では、ｌ=ｒであり、光到達体積Vは（4/3π×ｌ3）/2）となる。
したがって、光到達範囲パラメータは、光到達面積Ｓ、光到達距離ｌ、最小光到達距離ｌ2、光到達面積Ｓと最小光到達距離ｌ2、最大光到達距離ｌ1と最小光到達距離ｌ2の比又は差、光到達体積Ｖ、光到達体積Ｖと最小光到達距離ｌ2、もしくは、それらの組み合わせであってもよい。
制御部１０３は、算出された光到達範囲パラメータ（光到達距離ｌ、及び、光到達面積Ｓ等）に基づいて、血中の脂質濃度を算出する。
血中の脂質濃度の変化に伴い、照射光の拡散の面積が小さくなる。これは、血中の脂質粒子による光の散乱増加に伴い、光の拡散距離が低下するためであると判断できる。したがって、脂質濃度算出部１０４は、光到達距離ｌ若しくは光到達面積Ｓから、血中の脂質濃度を算出する。この手法は、特に毛細血管などの情報のみでも計測可能であることから、計測部位に依存しない。
図１２Ｂに示すように、脂質濃度変化量と光到達面積Ｓの間には、相関係数0.875という、良好な関係が得られていることから、少なくとも個人内変動においては、あらかじめ定めた相関係数から算出することができる。
また、制御部１０３は、光到達範囲パラメータから散乱係数を算出した後、脂質濃度を算出してもよい。臨床現場において、濃度と濁度とは同義で使われることがあり、本発明における濃度には濁度が含まれる。よって、制御部１０３は、その算出結果として、濃度のみならず単位量当たりの粒子数やホルマジン濁度あるいは散乱係数とすることができる。
図５は、図２のＸ方向から見た、生体表面の光到達範囲Ｆを示す図である。図に示すように、静脈を介した場合、照射部１０１からの光は同心円状に拡散せず、光到達範囲Ｆは、生体表面で、最大光到達距離ｌ1と最小光到達距離ｌ2を有する歪みを生じた形状になる。ここで、制御部１０３は、最小光到達距離ｌ2から血中の脂質濃度を算出する。この手法は、静脈を介した場合に計測可能な手法である。
また、制御部１０３は、光到達面積Ｓと最小光到達距離ｌ2とから脂質濃度を算出してもよい。これにより、静脈を含む計測部位でも、静脈および毛細血管の情報を総合的に取得できる。
また、制御部１０３は、最大光到達距離ｌ1と最小光到達距離ｌ2の比又は差を取ることで、より静脈情報としての精度を高めてもよい。さらに、制御部１０３は、最大光到達距離ｌ1と最小光到達距離ｌ2から、光到達範囲Ｆの楕円率を求め、または、光到達範囲Ｆの楕円面積を求めることにより、静脈情報としての精度を高めてもよい。
以上のような構成を備える脂質計測装置１００において、予め設定されているプログラムに基づいて、脂質計測装置１００は脂質計測処理を実行する。図６は、実施形態の脂質計測処理のフローチャートである。
照射工程（Ｓ１０１）では、照射部１０１が、生体の照射位置に対して連続光を照射する。
光強度検出工程（Ｓ１０２）では、光強度検出部１０２が、照射位置周辺における生体から放出される光強度を検出し、この光強度に基づいて生体内の光到達範囲Ｆを検出する。光強度検出工程で検出された光到達範囲Ｆは、パラメータ算出工程へと送られる。
パラメータ算出工程（Ｓ１０３）では、制御部１０３が、光到達範囲Ｆに基づき、所定の光到達範囲パラメータを算出する。光到達範囲パラメータは、光到達範囲Ｆの面積Ｓ、もしくは、光到達範囲Ｆの体積Ｖ、あるいは、光到達範囲Ｆにおける照射位置Ｅから光到達範囲Ｆの外周（外縁）までの距離ｌでもよい。光到達範囲パラメータは、最小光到達距離ｌ2のみ、光到達面積Ｓと最小光到達距離ｌ2、光到達体積Ｖと最小光到達距離ｌ2もしくは、最大光到達距離ｌ1と最小光到達距離ｌ2の比又は差、もしくは、それらの組み合わせであってもよい。算出した光到達範囲パラメータは、脂質濃度算出工程へと送られる。
脂質濃度算出工程（Ｓ１０４）では、制御部１０３が、光到達範囲パラメータに基づいて、血中の脂質濃度を算出する。脂質濃度算出工程では、光到達範囲パラメータから散乱係数を算出した後、脂質濃度を算出してもよい。
以上説明したように、本実施形態の脂質計測装置及び方法によれば、生体から放出される光強度の２次元情報を取得することで、静脈情報及び毛細血管情報を取得し、計測者の熟練がなくとも、容易に非侵襲脂質計測が可能となる。
次に、他の実施形態の脂質計測装置について説明をする。なお、他の実施形態の脂質計測装置の構成は、上記実施形態の脂質計測装置の構成と共通する部分もあるため、相違する部分を主に説明する。
上記実施形態では、照射部１０１と光強度検出部１０２と制御部１０３とを一体として構成した例を示したが、これに限られず、照射部１０１と光強度検出部１０２をユーザー装置として構成し、制御部１０３を、ユーザー装置に接続したサーバー装置に設けたシステムとしてもよい。
図７は、実施形態の脂質計測システムの構成を示す図である。システムは、脂質計測装置２００と、アクセスポイント３００と、ユーザー装置４００とを有する。
脂質計測装置２００は、ユーザー装置４００から送信された光強度に基づいて所定の処理を行い、脂質濃度を算出するための装置であり、具体的には、パーソナルコンピュータや、装置の台数や送受信するデータ量によってはサーバー装置が適宜用いられる。
ユーザー装置４００は、ユーザーが所持する装置であり、単独の装置である場合もあり、スマートフォン、携帯電話、腕時計等に搭載される場合もある。また、照射部４０１、光強度検出部４０２、通信部４０４として、スマートフォンや携帯電話に備わるカメラや照明、通信機能等、を使用してもよい。
ユーザー装置４００は、光を照射する照射部４０１と光強度検出部４０２と通信部４０４とを有する。通信部４０４は、光強度検出部４０２で検出された光強度を送信する。照射部４０１と光強度検出部４０２の機能・動作については上述した。
脂質計測装置２００は、通信部２０４ａと制御部２０３とを有する。通信部２０４は、通信部４０４から送信された光強度をアクセスポイント３００を介して受信し、制御部２０３へ送信する。
次に、脂質計測装置２００の制御系の構成について説明する。図８は実施形態の脂質計測装置２００のブロック図である。システムバス２０９を介して、ＣＰＵ（Central Processing Unit）２０４、ＲＯＭ（Read Only Memory）２０５、ＲＡＭ（Random Access Memory）２０６、記憶部２０７、及び、通信部（外部Ｉ／Ｆ（Interface））２０８が接続される。ＣＰＵ２０４とＲＯＭ２０５とＲＡＭ２０６とで制御部（コントローラー）２０３を構成する。
ＲＯＭ２０５は、ＣＰＵ２０４により実行されるプログラムや閾値を予め記憶する。
ＲＡＭ２０６は、ＣＰＵ２０４が実行するプログラムを展開するエリアと、プログラムによるデータ処理の作業領域となるワークエリアなどの様々なメモリエリア等を有する。
記憶部２０７は、予め用意された、静的パラメータ及び動的パラメータの適切な数値範囲のデータを記憶する。記憶部２０７は、ＨＤＤ（Hard Disk Drive）や、フラッシュメモリや、ＳＳＤ（Solid State Drive）等の、不揮発性に記憶する内部メモリーでよい。
通信部（外部Ｉ／Ｆ）２０８は、例えばクライアント端末（ＰＣ）などの外部装置と通信するためのインターフェースである。外部Ｉ／Ｆ２０８は、外部装置とデータ通信を行うインターフェースであれば良く、たとえば、外部装置にローカルに接続する機器（ＵＳＢメモリ等）であっても良いし、ネットワークを介して通信するためのネットワークインターフェイスであっても良い。制御部２０３の機能・動作については上述した。
なお、実施形態では、ユーザー装置４００から脂質計測装置２００へ、アクセスポイント３００を介して光強度を送信したが、これに限られず、ユーザー装置４００と脂質計測装置２００とが、アクセスポイントを介さずに直接接続し、有線通信や無線通信等の手段により光強度を送信してもよい。
以下に、本発明の実施例について説明するが、本発明は、下記の実施例に限定されない。
本実施例の脂質計測装置は、照射光が生体内の血中脂質に反射し、散乱して、生体から放出される光強度の２次元情報を取得することで、静脈情報及び毛細血管情報を取得し、計測者の熟練がなくとも、容易に非侵襲脂質計測を可能とする。
図９は、LED（照射部１０１）を生体の皮膚に直接当て、赤外線カメラ（光強度検出部１０２）を使用し撮影した結果を示す図である。図９に示すように、LED（照射部１０１）からの照射光は、生体内で同心円状に光は拡散することが確認される。
図１０は、生体の皮膚の同一部位で脂肪負荷後（血液濁度上昇後）に計測した結果を示す図である。
図１０では、図９と同様に、LED（照射部１０１）からの照射光は、生体内で同心円状に拡散するが、図９と比較して、光の周辺への広がりが小さくなっていることが確認できる。ここで計測したデータは、目視で静脈が見えない部分を計測している。
図１１は、前腕部の静脈近傍にて計測した結果を示す図である。静脈では血液による光の減衰と思われ現象が確認され、同心円ではなく歪みを生じた拡散をしていることが確認できる。
これらの得られる情報から、以下の手法により脂質濃度の算出が可能となる。
（１）光の拡散の光到達面積Ｓから脂質濃度を算出する手法（手法１）
（２）静脈による光の拡散の光到達範囲Ｆの歪みから脂質濃度を算出する手法（手法２）（３）光の拡散の光到達面積Ｓから脂質濃度を算出する手法（手法３）
以下に、上記各手法により脂質濃度を算出する方法について説明する。
（１）光の拡散の光到達面積Ｓから脂質濃度を算出する手法（手法１）
この手法は、静脈上以外の部分を分析することで毛細血管などの情報のみでも計測可能であることから、計測部位に依存しない。また、簡易的に光到達範囲あるいは光到達面積を光到達距離として分析してもよい。
図１２は、脂肪負荷試験による脂質濃度の変動と光到達面積Ｓを比較した図である。図１２Ａは、脂肪負荷時におけるTG変化量の時間変化と光到達面積の時間変化をプロットしたものである。図１２Ｂは、図１２ＡのTG変化量と光到達面積の相関を示したものである。図１２Ａに見られるように、脂質濃度の増加に伴い、光到達面積Ｓが減少することが確認できる。これは、脂質粒子による散乱増加に伴い、光の拡散距離が低下したためと判断できる。図１２Ｂから、TG変化量と光到達面積との間には相関度0.875の相関があることがわかる。
（２）静脈による光の拡散の光到達範囲Ｆの歪みから脂質濃度を算出する手法（手法２）
静脈を介した場合、光は同心円状に拡散せず、光到達範囲Ｆは歪みを生じた形状になる。ここで、光入射地点から光到達点の最大光到達距離l1と最小光到達距離l2を比較した。
図１３は、最小光到達距離ｌ2と脂質濃度との間の関連性を示す図である。図１３Ａは、脂肪負荷時におけるTG変化量の時間変化と最小光到達距離ｌ2の時間変化をプロットしたものである。図１３Ｂは、図１３ＡのTG変化量と最小光到達距離ｌ2の相関を示したものである。図１３Ａに見られるように、脂質濃度の増加に伴い、最小光到達距離ｌ2が減少することが確認できる。これは、脂質粒子による散乱増加に伴い、光の拡散距離が低下したためと判断できる。図１３Ｂから、TG変化量と最初光到達距離ｌ2との間には相関度0.877の相関があることがわかる。
（３）光の拡散の光到達体積Ｖから脂質濃度を算出する手法（手法３）
この手法は、特に毛細血管などの情報のみでも計測可能であることから、計測部位に依存しない。
図１４は、脂肪負荷試験による脂質濃度の変動と光到達体積Ｖを比較した図である。図１４Ａは、脂肪負荷時におけるTG変化量の時間変化と光到達面積の時間変化をプロットしたものである。図１４Ｂは、図１４ＡのTG変化量と光到達面積の相関を示したものである。図１４Ａに見られるように、脂質濃度の増加に伴い、光到達体積Ｖが減少することが確認できる。これは、脂質粒子による散乱増加に伴い、光の拡散距離が低下したためと判断できる。図１４Ｂから、TG変化量と光到達体積Ｖとの間には相関度0.851の相関があることがわかる。
また、手法１と手法２を合わせ、静脈を含む計測部位でも、光到達面積Ｓを算出し、静脈および毛細血管の情報を総合的に取得することも可能である。
また、手法２においては最大光到達距離l1と最小光到達距離l2の比又は差を取ることで、より静脈情報としての精度を高めることが可能である。さらに、手法２においては最大光到達距離l1と最小光到達距離l2から楕円率を求め、または楕円面積より静脈情報としての精度を高めることが可能である。
また、静脈情報の正確性を高める場合は、照射部１０１による照射光の入射点を増やし、複数点からの情報により静脈位置を特定することも可能である。
図１５は、図２と異なる照射部１０１と光強度検出部１０２の配置を示す図であり、図１６は、図１５に示す方法により撮影した例を示す図である。
図１６は、照射部１０１にレーザー使用し、レーザーを広範囲に照射し、レーザーのスペックルを計測することで、毛細血管（光到達深度１ｍｍ程度）の血流を計測したものである。
光の到達深度は、光源の光量を調節するなどでも良い。
図１７は、体温、脈などの影響を考慮し、計測時は同一姿勢で、安静状態で計測した結果を示す図である。
図１７Ａは、脂肪負荷時におけるTG変化量の時間変化と流量の時間変化をプロットしたものである。図１７Ｂは、図１７ＡのTG変化量と流量の相関を示したものである。図１７Ａに見られるように、脂質濃度の増加変化に伴い、流量が減少することが確認できる。図１７Ｂから、TG変化量と流量との間には相関度0.757の相関があることがわかる。この結果からも、静脈以外の血液情報から、脂質濃度を算出することが可能であることがわかった。
参考文献に記載の方法などを使用した静脈情報と比較し、代謝情報をより正確に得ることができる。また、光源を接触した場合と、非接触の場合を比較することで、静脈のみの情報を得ることも可能となる。
１００ 脂質計測装置
１０１ 照射部
１０２ 光強度検出部
１０３ 制御部
【０００８】
＝（４／３π×ａ×ｂ×ｃ）／２と、計算することもできる。
［００４５］
式中のａ，ｂ，ｃは、それぞれ球体の座標軸ｘ，ｙ，ｚ方向のそれぞれが９０°に交わる半径である。また、光到達範囲に歪みがない場合は、ａ＝ｂ＝ｃとなるため、図２では、ｌ＝ｒであり、光到達体積Ｖは（４／３π×ｌ３）／２）となる。
［００４６］
したがって、光到達範囲パラメータは、光到達面積Ｓ、光到達距離ｌ、最小光到達距離ｌ２、光到達面積Ｓと最小光到達距離ｌ２、最大光到達距離ｌ１と最小光到達距離ｌ２の比又は差、光到達体積Ｖ、光到達体積Ｖと最小光到達距離ｌ２、もしくは、それらの組み合わせであってもよい。さらには、取得データ全体ではなく最適部位のみの光到達面積Ｓ、光到達距離ｌ、最小光到達距離ｌ２、光到達面積Ｓと最小光到達距離ｌ２、最大光到達距離ｌ１と最小光到達距離ｌ２の比又は差、光到達体積Ｖ、光到達体積Ｖと最小光到達距離ｌ２、もしくは、それらの組み合わせであってもよい。
［００４７］
制御部１０３は、算出された光到達範囲パラメータ（光到達距離ｌ、及び、光到達面積Ｓ等）に基づいて、血中の脂質濃度を算出する。
［００４８］
血中の脂質濃度の変化に伴い、照射光の拡散の面積が小さくなる。これは、血中の脂質粒子による光の散乱増加に伴い、光の拡散距離が低下するためであると判断できる。したがって、脂質濃度算出部１０４は、光到達距離ｌ若しくは光到達面積Ｓから、血中の脂質濃度を算出する。この手法は、特に毛細血管などの情報のみでも計測可能であることから、計測部位に依存しない。
［００４９］
図１２Ｂに示すように、脂質濃度変化量と光到達面積Ｓの間には、相関係数０．８７５という、良好な関係が得られていることから、少なくとも個人内変動においては、あらかじめ定めた相関係数から算出することができる。
［００５０］
また、制御部１０３は、光到達範囲パラメータから散乱係数を算出した後、脂質濃度を算出してもよい。臨床現場において、濃度と濁度とは同義で使われることがあり、本発明における濃度には濁度が含まれる。よって、制御部１０３は、その算出結果として、濃度のみならず単位量当たりの粒子数やホルマジン濁度あるいは散乱係数とすることができる。
［００５１］
図５は、図２のＸ方向から見た、生体表面の光到達範囲Ｆを示す図である。図に示すように、静脈を介した場合、照射部１０１からの光は同心円状に拡散せず、光到達範囲Ｆは、生体表面で、最大光到達距離ｌ１と最小光到達距離ｌ２を有する歪みを生じた形状になる。ここで、制御部１０３は、最小光到

Claims (18)

  1. 生体外から生体内に向けて、生体の所定の部位に所定の光強度で光を照射する照射部と、
    前記生体から放出される光強度に基づいて生体内の光到達範囲を検出する光強度検出部と、
    前記光到達範囲に基づいて所定の光到達範囲パラメータを算出し、前記光到達範囲パラメータに基づいて生体内の脂質濃度を算出する制御部と、
    を有することを特徴とする脂質計測装置。
  2. 前記光到達範囲パラメータは、前記光到達範囲の面積に基づくことを特徴とする請求項１に記載の脂質計測装置。
  3. 前記光到達範囲パラメータは、前記光到達範囲における前記照射位置から前記光到達範囲の外周までの距離に基づくことを特徴とする請求項１または２に記載の脂質計測装置。
  4. 前記光到達範囲パラメータは、前記光到達範囲の体積に基づくことを特徴とする請求項１から３のいずれかに記載の脂質計測装置。
  5. 前記光到達範囲パラメータは、前記光到達範囲における前記照射位置から前記光到達範囲の外周までの最大距離と最小距離の比又は差を含むことを特徴とする請求項１から４のいずれかに記載の脂質計測装置。
  6. 前記制御部は、
    前記光到達範囲パラメータから散乱係数を算出した後、脂質濃度を算出することを特徴とする請求項１から５のいずれかに記載の脂質計測装置。
  7. 生体外から生体内に向けて、生体の所定の部位に所定の光強度で光を照射する照射工程と、
    前記生体から放出される光強度に基づいて生体内の光到達範囲を検出する光強度検出工程と、
    前記光到達範囲に基づいて所定の光到達範囲パラメータを算出するパラメータ算出工程と、
    前記光到達範囲パラメータに基づいて生体内の脂質濃度を算出する脂質濃度算出工程と、
    を有することを特徴とする脂質計測方法。
  8. 前記光到達範囲パラメータは、前記光到達範囲の面積に基づくことを特徴とする請求項７に記載の脂質計測方法。
  9. 前記光到達範囲パラメータは、前記光到達範囲における前記照射位置から前記光到達範囲の外周までの距離に基づくことを特徴とする請求項７または８に記載の脂質計測方法。
  10. 前記光到達範囲パラメータは、前記光到達範囲の体積に基づくことを特徴とする請求項７から９のいずれかに記載の脂質計測方法。
  11. 前記光到達範囲パラメータは、前記光到達範囲における前記照射位置から前記光到達範囲の外周までの最大距離と最小距離の比又は差を含むことを特徴とする請求項７から１０のいずれかに記載の脂質計測方法。
  12. 前記脂質濃度算出工程では、
    前記光到達範囲パラメータから散乱係数を算出した後、脂質濃度を算出することを特徴とする請求項７から１１のいずれかに記載の脂質計測方法。
  13. 生体外から生体内に向けて、生体の所定の部位に所定の光強度で光を照射する照射部と、
    前記生体から放出される光強度に基づいて生体内の光到達範囲を検出する光強度検出部と、前記光強度検出部により検出された光到達範囲を送信する通信部とを有するユーザー装置に、通信可能に接続される脂質計測装置であって、
    前記ユーザー装置から送信された前記光到達範囲に基づいて所定の光到達範囲パラメータを算出し、前記光到達範囲パラメータに基づいて生体内の脂質濃度を算出する制御部、を有することを特徴とする脂質計測装置。
  14. 前記光到達範囲パラメータは、前記光到達範囲の面積に基づくことを特徴とする請求項１３に記載の脂質計測装置。
  15. 前記光到達範囲パラメータは、前記光到達範囲における前記照射位置から前記光到達範囲の外周までの距離に基づくことを特徴とする請求項１３または１４に記載の脂質計測装置。
  16. 前記光到達範囲パラメータは、前記光到達範囲の体積に基づくことを特徴とする請求項１３から１５のいずれかに記載の脂質計測装置。
  17. 前記光到達範囲パラメータは、前記光到達範囲における前記照射位置から前記光到達範囲の外周までの最大距離と最小距離の比又は差を含むことを特徴とする請求項１３から１６のいずれかに記載の脂質計測装置。
  18. 前記制御部は、
    前記光到達範囲パラメータから散乱係数を算出した後、脂質濃度を算出することを特徴とする請求項１３から１７のいずれかに記載の脂質計測装置。
JP2018565525A 2017-01-31 2018-01-29 脂質計測装置及びその方法 Active JP6894087B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2017016012 2017-01-31
JP2017016012 2017-01-31
PCT/JP2018/002679 WO2018143119A1 (ja) 2017-01-31 2018-01-29 脂質計測装置及びその方法

Publications (2)

Publication Number Publication Date
JPWO2018143119A1 true JPWO2018143119A1 (ja) 2019-11-21
JP6894087B2 JP6894087B2 (ja) 2021-06-30

Family

ID=63040646

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018565525A Active JP6894087B2 (ja) 2017-01-31 2018-01-29 脂質計測装置及びその方法

Country Status (4)

Country Link
US (1) US20200245929A1 (ja)
JP (1) JP6894087B2 (ja)
TW (1) TWI750309B (ja)
WO (1) WO2018143119A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3988015A1 (en) * 2017-02-14 2022-04-27 Medical Photonics Co., Ltd. Scatterer measuring apparatus

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007117221A (ja) * 2005-10-25 2007-05-17 Hirato Koji Ｆｔｉｒを用いた非侵襲血液検査「光人間ドックシステム」
JP2010048703A (ja) * 2008-08-22 2010-03-04 Hokkaido Univ 血清脂質の測定方法及び測定装置
WO2014087825A1 (ja) * 2012-12-06 2014-06-12 国立大学法人北海道大学 非侵襲型生体脂質濃度計測器、非侵襲型生体脂質代謝機能計測器、非侵襲による生体脂質濃度計測方法および非侵襲による生体脂質代謝機能検査方法
JP6029128B1 (ja) * 2016-05-18 2016-11-24 メディカルフォトニクス株式会社 血中脂質濃度計測装置及びその作動方法
JP2017009398A (ja) * 2015-06-20 2017-01-12 関根　弘一 光学式生体情報測定装置および生体情報測定方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SG134283A1 (en) * 2006-01-24 2007-08-29 Ind Tech Res Inst Biomarkers for liver fibrotic injury
EP2841951B1 (en) * 2012-04-25 2019-12-11 Biodesy, Inc. Methods for detecting allosteric modulators of proteins

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007117221A (ja) * 2005-10-25 2007-05-17 Hirato Koji Ｆｔｉｒを用いた非侵襲血液検査「光人間ドックシステム」
JP2010048703A (ja) * 2008-08-22 2010-03-04 Hokkaido Univ 血清脂質の測定方法及び測定装置
WO2014087825A1 (ja) * 2012-12-06 2014-06-12 国立大学法人北海道大学 非侵襲型生体脂質濃度計測器、非侵襲型生体脂質代謝機能計測器、非侵襲による生体脂質濃度計測方法および非侵襲による生体脂質代謝機能検査方法
JP2017009398A (ja) * 2015-06-20 2017-01-12 関根　弘一 光学式生体情報測定装置および生体情報測定方法
JP6029128B1 (ja) * 2016-05-18 2016-11-24 メディカルフォトニクス株式会社 血中脂質濃度計測装置及びその作動方法

Also Published As

Publication number Publication date
US20200245929A1 (en) 2020-08-06
TW201837450A (zh) 2018-10-16
JP6894087B2 (ja) 2021-06-30
WO2018143119A1 (ja) 2018-08-09
TWI750309B (zh) 2021-12-21

Similar Documents

Publication Publication Date Title
US20240057883A1 (en) Deep tissue flowmetry using diffuse speckle contrast analysis
JP2022514226A (ja) 組織成分の非侵襲的検出方法、装置、システム、及びウェアラブルデバイス
US20210196135A1 (en) Blood vessel detection device and method therefor
TWI770119B (zh) 散射體測量裝置及其方法
JP6894088B2 (ja) 散乱体濃度計測装置及びその方法
JP6894087B2 (ja) 脂質計測装置及びその方法
CN106580328B (zh) 基于血红细胞浓度测量的烧伤程度与面积检测系统及方法
US20200359938A1 (en) Lipid measurement device and method therefor
JP6948069B2 (ja) 脂質計測装置及びその方法
JP2020022722A (ja) 癌発生疑い部位特定装置
US20210059575A1 (en) Lipid concentration measurement device and method therefor
JP6991634B1 (ja) 脂質濃度計測装置、プログラム、及び、方法
JP2019010346A (ja) 被検体情報取得装置および超音波探触子
JP2023004122A (ja) 高感度粒子濃度計測装置
WO2020080409A1 (ja) 粒子濃度計測装置、粒子濃度計測プログラム、及び、粒子濃度計測方法

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20190719

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200814

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210316

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210423

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: 20210511

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210524

R150 Certificate of patent or registration of utility model

Ref document number: 6894087

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150