ホーム > 特許ランキング > エバ・ジャパン 株式会社
知財求人 - 知財ポータルサイト「IP Force」
▶ エバ・ジャパン 株式会社の特許一覧 ▶ 国立大学法人東京農工大学の特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2022-10-04
(45)【発行日】2022-10-13
(54)【発明の名称】生体情報算出装置及び生体情報算出方法
(51)【国際特許分類】
A61B 5/1455 20060101AFI20221005BHJP
A61B 5/00 20060101ALI20221005BHJP
【ＦＩ】
A61B5/1455
A61B5/00 G
【請求項の数】
5
(21)【出願番号】P 2021029206
(22)【出願日】2021-02-25
(65)【公開番号】P2022130179
(43)【公開日】2022-09-06
【審査請求日】2021-03-05
【国等の委託研究の成果に係る記載事項】（出願人による申告）令和２年度、国立研究開発法人日本医療研究開発機構、長寿・障害総合研究事業 長寿科学研究開発事業「ＷＥＢカメラを用いた脈拍・呼吸・酸素飽和度の非接触型遠隔リアルタイムモニタリングシステムの開発」委託研究開発、産業技術力強化法第１７条の適用を受ける特許出願
(73)【特許権者】
【識別番号】510241915
【氏名又は名称】エバ・ジャパン 株式会社
(73)【特許権者】
【識別番号】504132881
【氏名又は名称】国立大学法人東京農工大学
(74)【代理人】
【識別番号】100120868
【弁理士】
【氏名又は名称】安彦 元
(72)【発明者】
【氏名】野呂 直樹
(72)【発明者】
【氏名】高良 洋平
(72)【発明者】
【氏名】西舘 泉
【審査官】▲高▼原 悠佑
(56)【参考文献】
【文献】特表２０２１－５０３３２３（ＪＰ，Ａ）
【文献】特表２００２－５２４１７７（ＪＰ，Ａ）
【文献】特開２０１６－０３０２１４（ＪＰ，Ａ）
【文献】特開２０００－３５０７０２（ＪＰ，Ａ）
【文献】米国特許出願公開第２０１６／０１６６１５０（ＵＳ，Ａ１）
【文献】中国特許出願公開第１０８０９０４１０（ＣＮ，Ａ）
【文献】特開２０１７－１８９５３６（ＪＰ，Ａ）
【文献】特開２０１２－０９５９４０（ＪＰ，Ａ）
【文献】川瀬達也，他7名，ＲＧＢ画像に基づく静脈閉塞に対する皮膚血流応答の評価，２０１１年春季 ＜第５８回＞応用物理学関係連合講演会［講演予稿集］，26p-BH-6，公益社団法人応用物理学会，2011年
【文献】KHATUN, F. et al.，Transcutaneous monitoring of hemoglobin derivatives during methemoglobinemia in rats using spectral diffuse reflectance，Journal of Biomedical Optics，2021年02月13日，Vol.26, No.3，033708，<DOI:https://doi.org/10.1117/1.JBO.26.3.033708>
【文献】NISHIDATE, I. et al.，Estimation of Melanin and Hemoglobin Using Spectral Reflectance Images Reconstructed from a Digital RGB Image by the Wiener Estimation Method，sensors，2013年06月19日，Vol.13, No.6，pp.7902-7915，<DOI:https://doi.org/10.3390/s130607902>
【文献】NISHIDATE, I. et al.，Noninvasive imaging of human skin hemodynamics using a digital red-green-blue camera，Journal of Biomedical Optics，2011年08月01日，Vol.16, No.8，086012，<DOI:https://doi.org/10.1117/1.3613929>
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ５／１４５５－５／１４６４
(57)【特許請求の範囲】
【請求項1】
生体である被測定対象の画像情報から前記被測定対象の１又は複数の経皮的動脈血酸素飽和度を含む生体情報を算出する生体情報算出装置であって、前記画像情報をＲＧＢ形式で取得する取得部と、
取得した前記画像情報の表色をＲＧＢ形式からＸＹＺ形式又はＬａｂ形式へと変換式を用いて変換する表色系変換部と、
前記ＸＹＺ形式又は前記Ｌａｂ形式の前記表色から酸素化血液色素濃度及び脱酸素化血液色素濃度を、推定式を用いて推定する色素濃度推定部と、
前記酸素化血液色素濃度及び前記脱酸素化血液色素濃度の少なくとも１つのパワースペクトル密度を算出する高速フーリエ変換部と、
前記パワースペクトル密度のピークから前記１又は複数の前記生体情報を抽出する生体情報抽出部と、
を備え、
前記経皮的動脈血酸素飽和度は、前記酸素化血液色素濃度の前記パワースペクトル密度における心拍由来の振動成分のピークと、前記脱酸素化血液色素濃度の前記パワースペクトル密度における心拍由来の振動成分のピークと、の比を所定の関数に入力することで出力として得られる、
生体情報算出装置。
【請求項2】
前記生体情報は、さらに脈拍数を含み、前記脈拍数は、前記酸素化血液色素濃度又は前記脱酸素化血液色素濃度の前記パワースペクトル密度における心拍由来の振動成分がピークとなった際の周波数から算出される、請求項１に記載の生体情報算出装置。
【請求項3】
前記生体情報は、さらに組織酸素飽和度を含み、前記組織酸素飽和度は、前記酸素化血液色素濃度を前記酸素化血液色素濃度及び前記脱酸素化血液色素濃度の和で除することで算出される、請求項２に記載の生体情報算出装置。
【請求項4】
前記生体情報は、さらに呼吸数を含み、前記呼吸数は、前記酸素化血液色素濃度又は前記脱酸素化血液色素濃度の前記パワースペクトル密度における呼吸由来の振動成分がピークとなった際の周波数から算出される、請求項３に記載の生体情報算出装置。
【請求項5】
生体である被測定対象の画像情報から前記被測定対象の１又は複数の経皮的動脈血酸素飽和度を含む生体情報を算出する生体情報算出装置が行う生体情報算出方法であって、前記画像情報をＲＧＢ形式で取得する第１のステップと、
取得した前記画像情報の表色をＲＧＢ形式からＸＹＺ形式又はＬａｂ形式へと変換式を用いて変換する第２のステップと、
前記ＸＹＺ形式又は前記Ｌａｂ形式の前記表色から酸素化血液色素濃度及び脱酸素化血液色素濃度を、推定式を用いて推定する第３のステップと、
前記酸素化血液色素濃度及び前記脱酸素化血液色素濃度の少なくとも１つのパワースペクトル密度を算出する第４のステップと、
前記パワースペクトル密度のピークから前記１又は複数の前記生体情報を抽出する第５のステップと、
を備え、
前記経皮的動脈血酸素飽和度は、前記酸素化血液色素濃度の前記パワースペクトル密度における心拍由来の振動成分のピークと、前記脱酸素化血液色素濃度の前記パワースペクトル密度における心拍由来の振動成分のピークと、の比を所定の関数に入力することで出力として得られる、
生体情報算出方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、生体情報算出装置及び生体情報算出方法に関する。
【背景技術】
【0002】
例えば特許文献１に開示された装置では、カメラと、２つの列状ＬＥＤ照明装置（以下、これを照明システムと呼んでもよい）と、制御部と、処理部と、を備え、哺乳類の経皮的動脈血酸素飽和度ＳｐＯ２を対象となる哺乳類に装置を接触させずに測定する。
【先行技術文献】
【特許文献】
【0003】
【文献】ＵＳ２０１７／０１２７９８８Ａ１号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
しかしながら特許文献１に記載されたような従来の手法においては、生体情報の測定には撮影環境を整えるために専用の照明システムが必要であるという課題がある。
【0005】
本発明の実施の形態の一態様は、スマートフォンなどの一般的なカメラ機能に用いられる機器を流用し汎用的な構成で生体情報を算出する生体情報算出装置を提供することを目的とする。
【課題を解決するための手段】
【0006】
生体である被測定対象の画像情報から被測定対象の１又は複数の生体情報を算出する生体情報算出装置であって、画像情報をＲＧＢ形式で取得する取得部と、取得した画像情報の表色をＲＧＢ形式からＸＹＺ形式又はＬａｂ形式へと変換する表色系変換部と、ＸＹＺ形式又はＬａｂ形式の表色から酸素化血液色素濃度及び脱酸素化血液色素濃度を推定する色素濃度推定部と、酸素化血液色素濃度及び脱酸素化血液色素濃度の少なくとも１つのパワースペクトル密度を算出する高速フーリエ変換部と、パワースペクトル密度のピークから１又は複数の生体情報を抽出する生体情報抽出部と、を備える生体情報算出装置を提供することを目的とする。
【0007】
生体である被測定対象の画像情報から前記被測定対象の１又は複数の生体情報を算出する生体情報算出装置が行う生体情報算出方法であって、画像情報をＲＧＢ形式で取得する第１のステップと、取得した画像情報の表色をＲＧＢ形式からＸＹＺ形式又はＬａｂ形式へと変換する第２のステップと、ＸＹＺ形式又はＬａｂ形式の表色から酸素化血液色素濃度及び脱酸素化血液色素濃度を推定する第３のステップと、酸素化血液色素濃度及び脱酸素化血液色素濃度の少なくとも１つのパワースペクトル密度を算出する第４のステップと、パワースペクトル密度のピークから１又は複数の生体情報を抽出する第５のステップと、を備える生体情報算出方法を提供することを目的とする。
【発明の効果】
【0008】
本発明の実施の形態の一態様によれば、スマートフォンなどの一般的なカメラ機能に用いられる機器を流用し汎用的な構成で生体情報を算出する生体情報算出装置を実現できる。
【図面の簡単な説明】
【0009】
【発明を実施するための形態】
【0010】
以下図面を用いて、本発明の実施の形態の一態様を詳述する。
【0011】
（本実施の形態）
まず図１を用いて利用段階における生体情報算出装置１について説明する。図１は、本実施の形態による生体情報算出装置１の構成を示すブロック図である。生体情報算出装置１は、生体である被測定対象の画像情報から測定対象の１又は複数の生体情報を算出する。
まず図１を用いて利用段階における生体情報算出装置１について説明する。図１は、本実施の形態による生体情報算出装置１の構成を示すブロック図である。生体情報算出装置１は、生体である被測定対象の画像情報から測定対象の１又は複数の生体情報を算出する。
【0012】
生体情報算出装置１は、中央演算装置２、主記憶装置３及び補助記憶装置４を備える。例えばスマートフォンは、生体情報算出装置１及び撮像装置１０を備える。中央演算装置２は、例えばＣＰＵ（Central Processing Unit）であって、主記憶装置３に記憶されたプログラムを呼び出すことで処理を実行する。主記憶装置３は、例えばＲＡＭ（Random Access Memory）であって、後述の取得部５、表色系変換部６、色素濃度推定部７、高速フーリエ変換部８及び生体情報抽出部９といったプログラムを記憶する。
【0013】
補助記憶装置４は、例えばＳＳＤ（Solid State Drive）やＨＤＤ（Hard Disk Drive）であって、変換係数テーブルＴＢ１や推定式係数テーブルＴＢ２や関数係数テーブルＴＢ３といったテーブルを記憶する。
【0014】
取得部５は、カメラなどの撮像装置１０から例えば１０秒間人の顔が撮像された複数の時系列的に連続する画像情報をＲＧＢ形式で取得する。表色系変換部６は、画像情報の表色をＲＧＢ形式からＸＹＺ形式又はＬａｂ形式へと変換式を用いて変換する。
【0015】
色素濃度推定部７は、ＸＹＺ形式又はＬａｂ形式の表色から酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲを、推定式を用いて推定する。高速フーリエ変換部８は、酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲの少なくとも１つのパワースペクトル密度を高速フーリエ変換により算出する。
【0016】
生体情報抽出部９は、パワースペクトル密度のピークＰＳＤＨｂＯ，ＰＳＤＨｂＲから１又は複数の生体情報をする、又は／及び、酸素化血液色素濃度ＣＨｂＯと酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲの和との比から生体情報を抽出する。ここでピークとは、パワースペクトル密度が極大かつ最大となることを指す。
【0017】
例えば、生体情報は、経皮的動脈血酸素飽和度ＳｐＯ２を含み、経皮的動脈血酸素飽和度ＳｐＯ２は、酸素化血液色素濃度ＣＨｂＯのパワースペクトル密度における心拍由来の振動成分のピークＰＳＤＨｂＯと、脱酸素化血液色素濃度ＣＨｂＲのパワースペクトル密度における心拍由来の振動成分のピークＰＳＤＨｂＲと、の比を所定の関数に入力することで出力として得られる。
【0018】
また例えば、生体情報は、脈拍数ＨＲを含み、脈拍数ＨＲは、酸素化血液色素濃度ＣＨｂＯ又は脱酸素化血液色素濃度ＣＨｂＲのパワースペクトル密度における心拍由来の振動成分がピークＰＳＤＨｂＯ，ＰＳＤＨｂＲとなった際の周波数ｆｐＨｂＯ，ｆｐＨｂＲから算出される。
【0019】
また例えば、生体情報は、組織酸素飽和度ＳｔＯ２を含み、組織酸素飽和度ＳｔＯ２は、酸素化血液色素濃度ＣＨｂＯを酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲの和で除することで算出される。
【0020】
また例えば、生体情報は、呼吸数ＢＲを含み、呼吸数ＢＲは、酸素化血液色素濃度ＣＨｂＯ又は脱酸素化血液色素濃度ＣＨｂＲのパワースペクトル密度における呼吸由来の振動成分がピークとなった際の周波数ｆｂＨｂＯ，ｆｂＨｂＲから算出される。
【0021】
次に図２を用いて表色系変換部６が変換の際に使用する変換式、色素濃度推定部７が推定の際に使用する推定式及び生体情報抽出部９が経皮的動脈血酸素飽和度ＳｐＯ２を算出する際に使用する所定の関数の詳細について説明する。
【0022】
図２は、本実施の形態による補助記憶装置４に記憶されている情報を示す。補助記憶装置４には、変換式、推定式及び所定の関数の係数を保持する変換係数テーブルＴＢ１、推定式係数テーブルＴＢ２及び関数係数テーブルＴＢ３が情報として記憶されている。
【0023】
表色系変換部６は、下記（１）式に示す変換式を用いてＲＧＢ形式からＸＹＺ形式へと表色を変換する。なおＲＧＢ形式からＬａｂ形式への変換は国際照明委員会が定義するＸＹＺ形式からＬａｂ形式への変換式を用いて行う。ＲＧＢ形式からＬａｂ形式への変換に必要なパラメータは変換係数テーブルＴＢ１に記憶されていてもよいものとする。
………（１）
パラメータＴ１１，Ｔ１２，Ｔ１３，Ｔ２１，Ｔ２２，Ｔ２３，Ｔ３１，Ｔ３２，Ｔ３３は、変換係数テーブルＴＢ１に記憶されており、例えば国際照明委員会が定義する数値とする。
【数1】
パラメータＴ１１，Ｔ１２，Ｔ１３，Ｔ２１，Ｔ２２，Ｔ２３，Ｔ３１，Ｔ３２，Ｔ３３は、変換係数テーブルＴＢ１に記憶されており、例えば国際照明委員会が定義する数値とする。
【0024】
色素濃度推定部７は、下記（２）式に示す推定式を用いて酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲをＸＹＺ形式の表色から推定する。
………（２）
【数2】
【0025】
なお色素濃度推定部７は、下記（３）式に示す推定式を用いて酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲをＬａｂ形式の表色から推定してもよい。
………（３）
なおメラニン色素濃度Ｃｍは、酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲをＸＹＺ形式又はＬａｂ形式の表色から推定する際に算出される。
【数3】
なおメラニン色素濃度Ｃｍは、酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲをＸＹＺ形式又はＬａｂ形式の表色から推定する際に算出される。
【0026】
パラメータＬ１１，Ｌ１２，Ｌ１３，Ｌ１０，Ｌ２１，Ｌ２２，Ｌ２３，Ｌ２０，Ｌ３１，Ｌ３２，Ｌ３３，Ｌ３０，Ｌ４１，Ｌ４２，Ｌ４０，Ｌ５１，Ｌ５２，Ｌ５０，Ｌ６１，Ｌ６２，Ｌ６０は、光伝搬モンテカルロシミュレーションを利用して求められる。
【0027】
具体的には例えば、説明変数（Ｘ，Ｙ，Ｚ）及び目的変数（Ｃｍ，ＣＨｂＯ，ＣＨｂＲ）の組み合わせを３００組用意し、式（２）で示す重回帰式を用いて重回帰分析を行うことでパラメータＬ１１，Ｌ１２，Ｌ１３，Ｌ１０，Ｌ２１，Ｌ２２，Ｌ２３，Ｌ２０，Ｌ３１，Ｌ３２，Ｌ３３，Ｌ３０は、算出される。
【0028】
また具体的には例えば、説明変数（ａ，ｂ）及び目的変数（Ｃｍ，ＣＨｂＯ，ＣＨｂＲ）の組み合わせを３００組用意し、式（３）で示す重回帰式を用いて重回帰分析を行うことでパラメータＬ４１，Ｌ４２，Ｌ４０，Ｌ５１，Ｌ５２，Ｌ５０，Ｌ６１，Ｌ６２，Ｌ６０は、算出される。
【0029】
生体情報抽出部９は、下記（４）式に示す所定の関数を用いて、経皮的動脈血酸素飽和度ＳｐＯ２を算出する。なお％はｖｏｌ．％と同意義とする。
………（４）
パラメータＡ，Ｂ，Ｃは、予め測定したＳｐＯ２とＰＳＤＨｂＯ／ＰＳＤＨｂＲとの組み合わせから決定されるものとする。なお（４）式で示したように、所定の関数は、適切な曲線で２つの値を関係づけるようなフィッテイング関数であって、例えば指数関数とする。
【数4】
パラメータＡ，Ｂ，Ｃは、予め測定したＳｐＯ２とＰＳＤＨｂＯ／ＰＳＤＨｂＲとの組み合わせから決定されるものとする。なお（４）式で示したように、所定の関数は、適切な曲線で２つの値を関係づけるようなフィッテイング関数であって、例えば指数関数とする。
【0030】
生体情報抽出部９は、下記（５）式を用いて組織酸素飽和度ＳｔＯ２を算出する。
………（５）
【数5】
【0031】
次に図４～図６を用いて、高速フーリエ変換部８が行う高速フーリエ変換、及び生体情報抽出部９が行う生体情報抽出の詳細について説明する。図３は、本実施の形態による一定時間における酸素化血液色素濃度を示すグラフ１５である。図４は、図３のグラフ１５に対して高速フーリエ変換を行うことで得られるグラフ２０である。
【0032】
【0033】
図３及び図５に示すように、酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲは、時系列的に連続する画像における一定範囲内のＸＹＺ形式又はＬａｂ形式の表色の例えば平均値から、（２）式又は（３）式を用いて推定される。
【0034】
画像における一定範囲とは、例えば人の顔における額部分とする。額部分を対象としているため、本実施の形態においては、色素濃度推定部７は、脳における酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲを推定する。
【0035】
また図４及び図６に示すように、パワースペクトル密度の振動成分は、呼吸由来成分と心拍由来成分とに分けることができる。パワースペクトル密度において周波数が５Ｈｚ未満の帯域を呼吸由来帯域とし、呼吸由来帯域における振動成分を呼吸由来成分と呼ぶ。
【0036】
図４や図６に示すようなパワースペクトル密度において周波数が０．７～４Ｈｚの帯域を心拍由来帯域とし、心拍由来帯域における振動成分を心拍由来成分と呼ぶ。なお呼吸由来帯域は、雑音の影響を除いた帯域である１Ｈｚ以下の帯域が特に好ましい。
【0037】
【0038】
【0039】
次に図７を用いて、本実施の形態における生体情報算出方法について説明する。図７は、本実施の形態による生体情報算出処理の処理手順を示すフローチャートである。まず取得部５は、画像情報をＲＧＢ形式で取得する（Ｓ１）。
【0040】
次に表色系変換部６は、取得した画像情報の表色をＲＧＢ形式からＸＹＺ形式又はＬａｂ形式へと変換式を用いて変換する（Ｓ２）。次に色素濃度推定部７は、ＸＹＺ形式又はＬａｂ形式の表色から酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲを、推定式を用いて推定する（Ｓ３）。
【0041】
次に高速フーリエ変換部８は、酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲの少なくとも１つのパワースペクトル密度を算出する（Ｓ４）。次に生体情報抽出部９は、パワースペクトル密度のピークＰＳＤＨｂＯ，ＰＳＤＨｂＲから１又は複数の生体情報を抽出する、又は／及び、酸素化血液色素濃度ＣＨｂＯと酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲの和との比から生体情報を抽出する（Ｓ５）。
【0042】
以上のように本実施の形態によれば、ＲＧＢ形式からＸＹＺ形式又はＬａｂ形式の表色系へと表色を変換し、ＸＹＺ形式又はＬａｂ形式の表色から色素濃度を推定し、色相濃度に対して高速フーリエ変換を行うことで、スマートフォンなどの一般的なカメラ機能に用いられる機器を流用し汎用的な構成で生体情報を算出することが可能となる。
【0043】
また本実施の形態では、例えば経皮的動脈血酸素飽和度ＳｐＯ２や組織酸素飽和度ＳｔＯ２や呼吸数ＢＲといった複数の生体情報を算出することができる。複数の生体情報を参照することで、生体情報算出装置１の使用者は、病態を鑑別することが可能となる。
【0044】
例えば具体的な病態として、脳などの臓器に低酸素状態がある場合には、呼吸数ＢＲが増え、経皮的動脈血酸素飽和度ＳｐＯ２や組織酸素飽和度ＳｔＯ２が低下する。またコロナウィルスに感染した場合、感染者が若年者の場合は経皮的動脈血酸素飽和度ＳｐＯ２が低下する前に呼吸数ＢＲが増え、感染者が高齢者の場合は呼吸数が増える前に経皮的動脈血酸素飽和度ＳｐＯ２が上昇する。
【0045】
また低血圧性ショックの場合や交感神経活動が急激に上昇するような場合といった、末梢灌流が大きく低下するような場合には、経皮的動脈血酸素飽和度ＳｐＯ２は測定不能又は不正確となるが、組織酸素飽和度ＳｔＯ２は正確である。このため、生体情報のいずれかの値が測定不能などになった場合においても生体情報算出装置１の使用者は、病態を鑑別することが可能となる。
【0046】
また本実施の形態においては、取得部５、表色系変換部６、色素濃度推定部７、高速フーリエ変換部８及び生体情報抽出部９は、プログラムとしたがこれに限らず論理回路でもよい。
【0047】
また取得部５、表色系変換部６、色素濃度推定部７、高速フーリエ変換部８、生体情報抽出部９、変換係数テーブルＴＢ１、推定式係数テーブルＴＢ２及び関数係数テーブルＴＢ３は１つの装置に実装されておらず、ネットワークで接続された複数の装置に分散して実装されていてもよい。
【0048】
本実施の形態においては、人の顔が画像として取得されるものとしたが、これに限らず、露出している肌を画像として取得してもよい。人の顔のほかに人の手のひらの画像が取得されてもよいし、同一人物における顔部及び手のひらの２つの画像が取得されてもよい。なお露出している肌の顔や手のひら以外の例として、手の甲や足なども挙げられる。
【0049】
また対象は人に限らなくてもよいものとする。なお手のひらの画像が取得される場合は、生体情報算出装置１は、末梢組織の経皮的動脈血酸素飽和度ＳｐＯ２及び組織酸素飽和度ＳｔＯ２を算出する。
【0050】
同一人物において複数個所の生体情報が算出される場合、例えば手のひらの画像情報を基に算出された経皮的動脈血酸素飽和度ＳｐＯ２は通常時と同じ正常値であって、額の画像情報を基に算出された組織酸素飽和度ＳｔＯ２が通常時と違う異常値であるような場合を鑑別することが可能となる。
【0051】
同一人物においてある部分の経皮的動脈血酸素飽和度ＳｐＯ２が正常値で、別の部分の組織酸素飽和度ＳｔＯ２が異常値であるような場合、別の部分への血液循環の低下を鑑別することが可能となる。血液循環の低下の原因は、血管狭窄や血管閉塞や出血や圧迫障害などが考えられる。
【0052】
本実施の形態においては、（２）式又は（３）式のように推定式の次数は１としたがこれに限らず、（１０）式又は（１１）式のように推定式は２次以上の多項式としてもよい。なお推定式の次数は３でもよいものとする。
……（１０）
……（１１）
【数10】
【数11】
【0053】
２次以上の多項式とした場合においてもパラメータＮ１０，Ｎ１１，Ｎ１２，Ｎ１３，Ｎ１４，Ｎ１５，Ｎ１６，Ｎ１７，Ｎ１８，Ｎ１９，Ｎ２０，Ｎ２１，Ｎ２２，Ｎ２３，Ｎ２４，Ｎ２５，Ｎ２６，Ｎ２７，Ｎ２８，Ｎ２９，Ｎ３０，Ｎ３１，Ｎ３２，Ｎ３３，Ｎ３４，Ｎ３５，Ｎ３６，Ｎ３７，Ｎ３８，Ｎ３９，Ｎ４０，Ｎ４１，Ｎ４２，Ｎ４３，Ｎ４４，Ｎ４５，Ｎ５０，Ｎ５１，Ｎ５２，Ｎ５３，Ｎ５４，Ｎ５５，Ｎ６０，Ｎ６１，Ｎ６２，Ｎ６３，Ｎ６４，Ｎ６５は、光伝搬モンテカルロシミュレーションを利用して求められる。推定式の次数を高くすることでＸＹＺ形式又はＬａｂ形式の表色から酸素化血液色素濃度ＣＨｂＯ及び脱酸素化血液色素濃度ＣＨｂＲをより正確に推定することが可能となる。
【0054】
本実施の形態においては、撮像装置１０は、スマートフォンに含まれる場合を例に挙げたが、これに限らず、ウェブカメラなどとしてもよい。なおウェブカメラにはリングＬＥＤなどが装着されていてもよいものとする。
【0055】
本実施の形態においては、経皮的動脈血酸素飽和度ＳｐＯ２を算出するフィッテイング関数は、例として（４）式のような指数関数を挙げたが、これに限らず、下記（１２）式のような２次関数であってもよい。
………（１２）
パラメータＤ，Ｅ，Ｆ，Ｇは、予め測定したＳｐＯ２とＰＳＤＨｂＯ／ＰＳＤＨｂＲとの組み合わせから決定されるものとする。
【数12】
パラメータＤ，Ｅ，Ｆ，Ｇは、予め測定したＳｐＯ２とＰＳＤＨｂＯ／ＰＳＤＨｂＲとの組み合わせから決定されるものとする。
【符号の説明】
【0056】
１……生体情報算出装置、２……中央演算装置、３……主記憶装置、４……補助記憶装置、５……取得部、６……表色系変換部、７……色素濃度推定部、８……高速フーリエ変換部、９……生体情報抽出部、１０……撮像装置。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】