(57)【特許請求の範囲】
【請求項１】
生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出する光源と、
複数のセンサが所定の配置で規則的に配設されており、前記光源から射出され前記生体を透過した前記測定光を当該複数のセンサで検出する検出部と、
前記検出部により検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行う解析部と、
を備え、
前記検出部は、複数のレンズが格子状に規則的に配設されたマイクロレンズアレイを利用したセンサにより、前記生体を透過した測定光を検出する、測定装置。
【請求項２】
前記マイクロレンズアレイと前記複数のセンサとの間には、検出する前記測定光のクロストークを防止する遮光体が設けられる、請求項１に記載の測定装置。
【請求項３】
前記解析部は、前記測定光の光量の時間変化に基づき、前記生体の活動に伴う脈動に関する情報として、前記測定領域に対応する前記生体内部の動脈の拍動に由来する脈波形を更に算出する、請求項２に記載の測定装置。
【請求項４】
前記解析部は、予め測定された前記生体における前記脈波形と、前記測定光の光量の時間変化に基づき算出された前記脈波形と、の類似度合いに応じて、前記測定位置を特定する、請求項３に記載の測定装置。
【請求項５】
前記解析部は、前記測定光の光量の時間変化に基づき算出された前記脈波形の最大振幅を与える位置を、前記測定位置として特定する、請求項３に記載の測定装置。
【請求項６】
前記解析部は、前記生体の活動に伴う脈動に関する情報として前記脈波形の振幅を測定するための前記測定位置を特定する、請求項４に記載の測定装置。
【請求項７】
前記解析部は、前記生体の活動に伴う脈動に関する情報として前記脈波形の形状を測定するための前記測定位置を特定する、請求項４に記載の測定装置。
【請求項８】
前記解析部は、前記光源から射出される前記測定光の波長毎に、前記解析処理を実施する、請求項１に記載の測定装置。
【請求項９】
前記解析部は、算出した前記脈波形を利用して、前記動脈における酸素飽和度を更に算出する、請求項３に記載の測定装置。
【請求項１０】
前記解析部は、前記複数のセンサにおける前記測定光の検出位置と前記光源との離隔距離に応じて、前記検出結果の時間軸における位相を補正する、請求項１に記載の測定装置。
【請求項１１】
前記解析部は、前記複数のセンサにおいて前記光源からの離隔距離が同一である検出位置で得られた前記検出結果を積算する、請求項１に記載の測定装置。
【請求項１２】
前記解析部は、互いに異なる時間に対応する複数の前記検出結果について、当該時間を補正したうえで合成する、請求項１に記載の測定装置。
【請求項１３】
生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出することと、
複数のセンサが所定の配置で規則的に配設された検出部により、光源から射出され前記生体を透過した前記測定光を検出することと、
検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行うことと、
を含み、
前記検出部は、複数のレンズが格子状に規則的に配設されたマイクロレンズアレイを利用したセンサにより、前記生体を透過した測定光を検出する、測定方法。
【請求項１４】
生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出する光源と、複数のセンサが所定の配置で規則的に配設されており、前記光源から射出され前記生体を透過した前記測定光を当該複数のセンサで検出する検出部と、を備え、前記検出部は、複数のレンズが格子状に規則的に配設されたマイクロレンズアレイを利用したセンサにより、前記生体を透過した測定光を検出する測定機器と通信可能なコンピュータに、
前記検出部により検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行う解析機能
を実現させるためのプログラム。
【請求項１５】
生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出する光源と、複数のセンサが所定の配置で規則的に配設されており、前記光源から射出され前記生体を透過した前記測定光を当該複数のセンサで検出する検出部と、を備え、前記検出部は、複数のレンズが格子状に規則的に配設されたマイクロレンズアレイを利用したセンサにより、前記生体を透過した測定光を検出する測定機器と通信可能なコンピュータに、
前記検出部により検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行う解析機能
を実現させるためのプログラムが記録された記録媒体。
【請求項１６】
生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出する光源と、
複数のレンズが格子状に規則的に配設されたマイクロレンズアレイを利用したセンサにより、前記光源から射出され前記生体を透過した前記測定光を検出する検出部と、
前記検出部により検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行うとともに、記測定光の光量の時間変化に基づき、前記生体の活動に伴う脈動に関する情報として、前記測定領域に対応する前記生体内部の動脈の拍動に由来する脈波形を算出する解析部と、
を備える、測定装置。
【請求項１７】
生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出することと、
複数のレンズが格子状に規則的に配設されたマイクロレンズアレイを利用したセンサにより、光源から射出され前記生体を透過した前記測定光を検出することと、
検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行うとともに、前記測定光の光量の時間変化に基づき、前記生体の活動に伴う脈動に関する情報として、前記測定領域に対応する前記生体内部の動脈の拍動に由来する脈波形を算出することと、
を含む、測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、測定装置、測定方法、プログラム及び記録媒体に関する。
【背景技術】
【０００２】
医療分野において、脈波形（動脈の拍動波形）を測定する技術が多用されている。脈波形は、その波形形状を解析することにより、循環器系の検査（例えば、動脈硬化度測定、ストレス測定）等に用いられるほか、パルスオキシメータ（動脈血中酸素飽和濃度測定装置）等にも利用されている。
【０００３】
なかでも、採血や穿刺を伴わない非侵襲で経皮的に脈波形を測定する技術の一つである光学式容積脈波法（Ｐｈｏｔｏｐｌｅｔｈｙｓｍｏｇｒａｐｈｙ）と呼ばれる技術（例えば、以下の特許文献１を参照。）は、人体に負担をかけずに簡便に測定が可能であるため、非常に広い範囲で利用されている。
【０００４】
光学式容積脈波法を利用したパルスオキシメータでは、酸化ヘモグロビンは赤外光（例えば、波長９４０ｎｍ前後）を吸収し、還元ヘモグロビンは赤色光（例えば、波長６６０ｎｍ前後）を吸収するという特性を利用し、測定した脈波形に基づいて酸素飽和濃度を決定する。すなわち、これら２種類の波長の光を人体に照射してそれぞれの波長における吸光度を測定し、得られた吸光度に基づいて、各ヘモグロビンの相対濃度を算出し、得られた相対濃度から酸素飽和濃度を算出する。
【０００５】
このような光学式容積脈波法を利用したパルスオキシメータには、測定部位である指を光源と光検出器とで挟み込む透過型のパルスオキシメータと、指等の測定部位の片側に光源と光検出器とを配置する反射型のパルスオキシメータという、２種類のタイプが存在する。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００３−１３５４３４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
上記のような２つのタイプの光学式容積脈波法を利用したパルスオキシメータにおいて、安定した脈波形を得るためには透過型のものが有利であるが、光学的な安定を保持するためには、強く、かつ、隙間なく測定部位を挟み込むことが求められるため、長時間の使用は使用者に負担がかかってしまう。一方、反射型のパルスオキシメータは、人体の片側にのみ接触するものであるため、使用者への低負荷な測定が可能であったが、光の効率が不利であるため、十分な精度の脈波形を得ることが困難であった。
【０００８】
そこで、本開示では、上記事情に鑑みて、使用者への負荷を低下させつつ、精度良く脈波形を測定することが可能な測定装置、測定方法、プログラム及び記録媒体を提案する。
【課題を解決するための手段】
【０００９】
本開示によれば、生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出する光源と、複数のセンサが所定の配置で規則的に配設されており、前記光源から射出され前記生体を透過した前記測定光を当該複数のセンサで検出する検出部と、前記検出部により検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行う解析部と、を備える測定装置が提供される。
【００１０】
また、本開示によれば、生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出することと、複数のセンサが所定の配置で規則的に配設された検出部により、前記光源から射出され前記生体を透過した前記測定光を検出することと、検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行うことと、を含む測定方法が提供される。
【００１１】
また、本開示によれば、生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出する光源と、複数のセンサが所定の配置で規則的に配設されており、前記光源から射出され前記生体を透過した前記測定光を当該複数のセンサで検出する検出部と、を備える測定機器と通信可能なコンピュータに、前記検出部により検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行う解析機能を実現させるためのプログラムが提供される。
【００１２】
また、本開示によれば、生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出する光源と、複数のセンサが所定の配置で規則的に配設されており、前記光源から射出され前記生体を透過した前記測定光を当該複数のセンサで検出する検出部と、を備える測定機器と通信可能なコンピュータに、前記検出部により検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行う解析機能を実現させるためのプログラムが記録された記録媒体が提供される。
【００１３】
また、本開示によれば、生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出する光源と、複数のレンズが格子状に規則的に配設されたマイクロレンズアレイを利用したセンサにより、前記光源から射出され前記生体を透過した前記測定光を検出する検出部と、前記検出部により検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行うとともに、記測定光の光量の時間変化に基づき、前記生体の活動に伴う脈動に関する情報として、前記測定領域に対応する前記生体内部の動脈の拍動に由来する脈波形を算出する解析部と、を備える測定装置が提供される。
【００１４】
また、本開示によれば、生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出することと、複数のレンズが格子状に規則的に配設されたマイクロレンズアレイを利用したセンサにより、前記光源から射出され前記生体を透過した前記測定光を検出することと、検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行うとともに、前記測定光の光量の時間変化に基づき、前記生体の活動に伴う脈動に関する情報として、前記測定領域に対応する前記生体内部の動脈の拍動に由来する脈波形を算出することと、を含む測定方法が提供される。
【００１５】
本開示によれば、生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光が射出され、複数のセンサが所定の配置で規則的に配設されたセンサにより、光源から射出され生体を透過した測定光が検出され、検出された検出結果を利用し、検出された測定光の光量の時間変化に基づいて、生体の活動に伴う脈動に関する情報を測定するための測定位置が測定領域の中から特定される。
【発明の効果】
【００１６】
以上説明したように本開示によれば、使用者への負荷を低下させつつ、精度良く脈波形を測定することが可能である。
【図面の簡単な説明】
【００１７】
【図１】人体の皮膚構造モデルの例について示した説明図である。
【図２Ａ】一般的なパルスオキシメータの構成について示した説明図である。
【図２Ｂ】一般的なパルスオキシメータの構成について示した説明図である。
【図３】拡張ランベルト・ベールの法則について説明するための説明図である。
【図４】酸素化ヘモグロビン及び還元ヘモグロビンの光吸収スペクトルの例を示したグラフ図である。
【図５】生体の皮膚組織と吸光度との関係を示した説明図である。
【図６】本開示の第１の実施形態に係る測定装置の構成の一例を示したブロック図である。
【図７】同実施形態に係る測定装置が備える測定部の一例を示した説明図である。
【図８】同実施形態に係る測定装置が備える測定部の一例を示した説明図である。
【図９】同実施形態に係る測定装置が備える測定部の一例を示した説明図である。
【図１０Ａ】同実施形態に係る測定装置について説明するための説明図である。
【図１０Ｂ】同実施形態に係る測定装置について説明するための説明図である。
【図１１】同実施形態に係る測定装置について説明するための説明図である。
【図１２】同実施形態に係る測定装置について説明するための説明図である。
【図１３】本開示の実施形態に係る測定装置のハードウェア構成の一例を示したブロック図である。
【発明を実施するための形態】
【００１８】
以下に添付図面を参照しながら、本開示の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
【００１９】
なお、説明は、以下の順序で行うものとする。
（１）人体の皮膚構造モデルについて
（２）パルスオキシメータの原理について
（３）第１の実施形態
（４）本開示の実施形態に係る測定装置のハードウェア構成について
（５）まとめ
【００２０】
（人体の皮膚構造モデルについて）
本開示の実施形態に係る測定装置、測定方法、プログラム及び記録媒体について説明するに先立ち、人体の皮膚構造をモデル化した皮膚構造モデルについて、図１を参照しながら簡単に説明する。図１は、人体の皮膚構造モデルの例を示した説明図である。
【００２１】
先述のように、人体内に存在するグルコース、アルブミン、ＡＧＥｓ（糖化最終産物）、コレステロール、酸化・還元ヘモグロビン等の血中・体液成分を非侵襲の光学測定により計測する技術が開発されている。
【００２２】
測定されたデータを解析するにあたっては、人体の皮膚構造をどのようにモデル化するかが重要となってくる。このような人体の皮膚構造モデルの例としては、図１に示したような３層モデルが存在する。
【００２３】
図１に示した３層モデルは当し、真皮層、皮膚の角質層より下部に位置する皮下組織を、表皮層、真皮層、皮下脂肪の３つの層にモデル化したものである。この３層モデルにおいて、角質層は、個人差はあるものの体表から内部方向に０．０１〜０．０２ｍｍ程度に相当し、表皮層は、体表から０．０４〜０．１５ｍｍ程度に相は、体表から１〜４ｍｍ程度に相当し、皮下脂肪は、体表から数ｍｍ〜数ｃｍ程度に相当する。
【００２４】
このような皮膚構造において、真皮層には、毛細血管が存在しており、この毛細血管の内部には、酸化ヘモグロビン、還元ヘモグロビンや、各種の血中成分が存在しており、皮下脂肪には、主に脂肪細胞が存在している。従って、上記のような成分を非侵襲の光学測定により計測する場合には、どのような皮膚構造モデルを考慮するかが重要になるのである。
【００２５】
ところが、人体の皮膚構造や、皮膚構造に含まれる各種体内成分は、被測定者本人の経年変化、肥満、日焼け、人種の差、性別、体質等により変動するものであり、個人差も大きいものである。このため、皮膚構造モデルに基づく体内成分の非侵襲計測では、図１に示したような皮膚構造モデルの変動の影響を強く受けるものであると言える。
【００２６】
従って、動脈の拍動に由来する脈波形を精度よく検出するためには、被測定者それぞれにおいて、適切な検出位置を特定し、特定した検出位置において検出を行うことが好ましい。そのため、本発明者は、以上のような事情を鑑みて、使用者への負荷を低下させつつ、精度良く脈波形を測定することが可能な技術について鋭意検討を行った。その結果、本発明者は、以下で説明するような本開示の実施形態に係る測定装置及び測定方法に想到した。
【００２７】
（パルスオキシメータの原理について）
続いて、本開示の実施形態に係る測定装置、測定方法、プログラム及び記録媒体について説明するに先立ち、光学式容積脈波法を用いた一般的なパルスオキシメータの原理について、図２Ａ、図２Ｂ及び図３を参照しながら、簡単に説明する。図２Ａ及び図２Ｂは、一般的なパルスオキシメータの構成について示した説明図である。図３は、拡張ランベルト・ベールの法則について説明するための説明図である。
【００２８】
パルスオキシメータは、経皮的に動脈血の酸素飽和度（ｐｅｒｃｕｔａｎｅｏｕｓ  ｏｘｙｇｅｎ  ｓａｔｕｒａｔｉｏｎ，ＳｐＯ２と呼ばれる。）を測定する機器である。このパルスオキシメータは、その構成要素として、一般的な脈波測定装置を含む。
【００２９】
パルスオキシメータは、測定プローブが接続された吸光度測定部で、生体情報を測定する。測定プローブは、図２Ａ及び図２Ｂに示したように、光源と光検出器とで構成されており、生体による光吸収の時間変化を測定する。測定プローブにより測定された光吸収に関する測定結果は、周波数フィルタへと出力され、周波数フィルタにおいて、動脈の拍動に由来する周波数成分（以下、単に「拍動成分」とも称する。）が分離される。その後、分離された拍動成分の振幅に基づいて、酸素飽和度が算出される。算出された酸素飽和度は、ユーザへと出力されることとなる。
【００３０】
パルスオキシメータでは、生体に向けて照射する光（入射光）に、複数の波長の光を用いることが求められる。そのため、測定プローブには複数の光源が配置され、時分割で切り替えて用いられる。また、入射光には、生体内に届きやすいという理由から、赤色光から近赤外光の帯域に属する波長の光が用いられることが多い。
【００３１】
パルスオキシメータの測定プローブでは、図２Ａ及び図２Ｂに示したように、光源から入射光を生体の皮膚面に照射し、生体内で反射したり拡散したりして生体から射出された出射光を、光検出器で検出する。ここで、図２Ａに示した透過型タイプのパルスオキシメータでは、光源と光検出器とは、生体の一部（例えば指等）を挟んで対向するように設けられており、光検出器は、生体を拡散しながら透過した出射光を検出する。また、図２Ｂに示した反射型タイプのパルスオキシメータでは、光源と光検出器とは、生体の一部の同じ側に設けられており、生体内を反射したり拡散したりして略Ｕ字形状に伝播した出射光が、光検出器により検出される。このとき、入射光は、生体内に存在する動脈や静脈やその他の体組織によって一部吸収されて、出射光として観測されることとなる。
【００３２】
パルスオキシメータでは、測定された実際のデータと、着目する生体内成分（すなわち、酸素化ヘモグロビンや、還元ヘモグロビン等）に起因する光吸収量とを関係づけるために、拡張ランベルト・ベールの法則を利用する。パルスオキシメータは、生体という光を拡散させる物体（光の散乱体）に着目し、生体内の光の伝播を考慮するものであるため、散乱・拡散の効果を考慮することができない通常のランベルト・ベールの法則を利用することはできない。そのため、パルスオキシメータは、以下の式１１に示したような拡張ランベルト・ベールの法則を利用して、得られた測定データの解析を実施する。以下、拡張ランベルト・ベールの法則について、図３を参照しながら簡単に説明する。
【００３３】
【数１】
0005948836.tif 000002
【００３４】
ここで、上記式１１において、
λ：着目する光の波長
Ａ（λ）：波長λにおける吸光度
Ｉ０（λ）：散乱体に入射した波長λの光の強度
Ｉ（λ）：散乱体を透過した波長λの光の検出強度
Ｇ（λ）：波長λの光の散乱による減衰量
εｉ（λ）：物質ｉの波長λの光に関する吸光係数であり、物質に固有の値である。
Ｃｉ：物質ｉの濃度
ｌｉ：波長λの光が物質ｉを伝播する際の平均光路長
である。
【００３５】
ここで、上記拡張ランベルト・ベールの法則を、図３に示したような層構造を有する散乱体に適用することを考える。以下では、層を特定するための添え字を改めてｉと記載することとし、層ｉに含まれる物質の数を添え字ｊで表すこととする。すると、図３に示したような層構造を有する散乱体における拡張ランベルト・ベールの法則は、下記式１２及び式１３のように表すことができる。
【００３６】
【数２】
0005948836.tif 000003
【００３７】
ここで、上記式１２及び式１３において、
λ：着目する光の波長
Ａ（λ）：波長λにおける吸光度
Ｉ０（λ）：散乱体に入射した波長λの光の強度
Ｉ（λ）：散乱体を透過した波長λの光の検出強度
Ｇ（λ）：波長λの光の散乱による減衰量
εｉ（λ）：層ｉの波長λの光に関する吸光係数
Ｃｉ：層ｉに含まれる物質の濃度
ｌｉ：波長λの光が層ｉを伝播する際の平均光路長
εｉｊ（λ）：層ｉに含まれる物質ｊの波長λの光に関する吸光係数
Ｃｉｊ：層ｉに含まれる物質ｊの濃度
である。
【００３８】
ここで、着目する生体内成分の吸光係数は、予め着目する生体内成分の吸収スペクトルを測定したり、公知のデータベースからデータを取得したりすることで、特定することができる。従って、着目する生体内成分の吸光係数は、これらのデータを利用することによって、既知量として取り扱うことができる。また、上記式１２における最左辺は、パルスオキシメータにより測定された、各波長に関する実際の測定データとなる。
【００３９】
パルスオキシメータで着目する血液中のヘモグロビンは、図４に例示したように、酸素との結合の有無によって吸光度が変化し、また、観測する波長によっても吸光度が異なる。従って、複数の波長で吸光度を測定することによって、酸素と結合していない還元ヘモグロビン（Ｈｂ）と、酸素と結合したヘモグロビン（酸素化ヘモグロビン：ＨｂＯ２）の比率を求めることができる。
【００４０】
血液中に含まれる総ヘモグロビンのうち、酸素化ヘモグロビンの割合を、血中酸素飽和度と呼ぶ。生体情報として特に有用なものは動脈血の酸素飽和度ＳａＯ２（ａｒｔｅｒｉａｌ  ｏｘｙｇｅｎ  ｓａｔｕｒａｔｉｏｎ）であるが、この酸素飽和度ＳａＯ２は、以下の式１４で算出することができる。なお、先述のＳｐＯ２は、ＳａＯ２を経皮的に測定したものである。
【００４１】
【数３】
0005948836.tif 000004
【００４２】
なお、上記式１４において、
ＳａＯ２：動脈血酸素飽和度
ＣＨｂＯ２：酸素化ヘモグロビン濃度
ＣＨｂ：還元ヘモグロビン濃度
である。
【００４３】
上述の通り、パルスオキシメータにおける測定プローブの光検出器で検出される出射光は、入射光が体内における反射・散乱の過程で、体組織や血液成分による吸収を受けたものである。出射光強度を解析することでＳｐＯ２を算出することが可能であるが、ＳｐＯ２は動脈血の酸素飽和度であるから、出射光から動脈血以外による光吸収の影響を除外することが求められる。
【００４４】
入射光に光吸収をもたらす要素は、図５に例示したように、動脈血、静脈血、その他の体組織の３種類に大別することができる。このとき、出射光は、下記の式１５に示すような光吸収を受ける。
【００４５】
【数４】
0005948836.tif 000005
【００４６】
ここで、上記式１５において、
λ：波長
ε：吸光係数
Ｃ：濃度
ｄ：光路長
である。
【００４７】
また、上記式１５において、最右辺第１項は、血液以外の成分に起因する光吸収を表しており、最右辺第２項は、静脈血に起因する光吸収を表しており、最右辺第３項は、動脈血に起因する光吸収を表しており、最右辺第４項は、生体内での拡散に起因する光吸収を表している。
【００４８】
パルスオキシメータは、上記３種類の要素のうち、動脈のみに拍動性があることを利用し、動脈血の吸光をその他の要素から分離する。すなわち、上記式１５を時間微分することで、拍動性を持たない（換言すれば、時間変化の無い）静脈及びその他の体組織による光吸収の影響を除去する。この微分操作は、信号処理においては、周波数フィルタによる直流成分の除去に相当するものであり、脈波形の抽出処理に他ならない。
【００４９】
上記式１４において、ＳａＯ２を算出するための未知数は、還元ヘモグロビン濃度（ＣＨｂ）と、酸素化ヘモグロビン濃度（ＣＨｂＯ２）の２種類であるため、２つの未知数を特定するために２つの測定結果を連立させることが求められる。従って、パルスオキシメータでは、少なくとも２つの波長を用いて測定が行われる。
【００５０】
いま、波長λ１，λ２の２種類の入射光で測定を行い、出射光強度の時間変化であるΔＯＤλ１，ΔＯＤλ２を求める場合を考える。この場合、２つの波長を用いて測定された出射光強度の時間変化は、上記式１５から、以下の式１６のように表すことができる。従って、未知数であるヘモグロビン濃度（ＣＨｂ）及び酸素化ヘモグロビン濃度（ＣＨｂＯ２）は、ヘモグロビン及び酸素化ヘモグロビンの吸光係数等と、測定結果とを用いて、以下の式１７に示したように算出することが可能となる。
【００５１】
【数５】
0005948836.tif 000006
【００５２】
従って、式１７を式１４に代入すると、以下の式１８を得ることができる。ここで、下記式１８において、パラメータα、β，Φは、以下の式１９ａ〜式１９ｃの通りである。
【００５３】
【数６】
0005948836.tif 000007
【００５４】
上記式１８の最右辺から明らかなように、ＳａＯ２の値は、パラメータΦに比例する関数として与えられることがわかる。パラメータΦは、上記式１９ｃのように、波長λ１と波長λ２で測定した脈波形の振幅の比である。また、パラメータα，βは、式１９ａ及び式１９ｂや、図４に示したように、ヘモグロビンの吸光係数から理論的に算出することが可能であるが、多くの場合には、事前に実験で求めた変換表に基づく校正を行うことで求めている。そのようにすることで、ランベルト・ベールの法則の成立する条件と、実際の生体における条件との乖離も含めた補正が可能となるからである。
【００５５】
このような方法がパルスオキシメータの動作原理であり、パルスオキシメータは、２種類の波長による測定結果を利用して、動脈血の酸素飽和度ＳｐＯ２を算出する。
【００５６】
（第１の実施形態）
＜測定装置の構成について＞
続いて、図６〜図１２を参照しながら、本開示の第１の実施形態に係る測定装置及び測定方法について、詳細に説明する。図６は、本実施形態に係る測定装置１０の構成について示したブロック図である。図７〜図９は、本実施形態に係る測定装置１０が備える測定部の一例を示した説明図である。図１０〜図１２は、本実施形態に係る測定装置１０について説明するための説明図である。
【００５７】
まず、図２を参照しながら、本実施形態に係る測定装置１０の全体構成について、詳細に説明する。
本実施形態に係る測定装置１０は、測定対象物である生体Ｂの少なくとも一部を所定の波長を有する光を用いて測定し、得られた測定結果に基づいて、生体の活動に伴う脈動に関する情報を測定するための測定位置を測定領域の中から特定する。また、本実施形態に係る測定装置１０は、得られた測定結果を利用して、生体の活動に伴う脈動に関する情報（例えば、脈波形に関する情報や、酸素飽和度に関する情報や、これらの情報から二次的に算出可能な情報等）を算出することもできる。
【００５８】
この測定装置１０は、図６に示したように、生体Ｂの測定領域を測定する測定部１０１と、制御部１０３と、解析部１０５と、記憶部１０７と、を主に備える。
【００５９】
［測定部１０１について］
以下では、まず、図７〜図９を参照しながら、本実施形態に係る測定部１０１の構成について、具体的に説明する。
本実施形態に係る測定部１０１は、図７に示したように、光源１１１と、検出部１１３と、から構成される。
【００６０】
○光源について
光源１１１は、所定の波長帯域に属する少なくとも１種類の測定光を、生体Ｂの測定領域に向かって射出する。この光源１１１は、測定光の射出面が生体Ｂと対向するように、所定のフレームＦに配設される。
【００６１】
光源１１１が射出する測定光の波長は、着目する生体内成分に応じて適宜設定することが可能であるが、例えば、酸素化ヘモグロビンを測定する場合には、波長９４０ｎｍ前後の近赤外光を使用することが可能であり、還元ヘモグロビンを測定する場合には、波長６６０ｎｍ前後の赤色光を使用することが可能である。また、光源１１１が９４０ｎｍ、９５０ｎｍのような波長の光を射出することで、皮下組織に存在する脂肪に関する知見を得ることができる。また、光源１１１が５６８ｎｍ、５８０ｎｍ、６６０ｎｍ、８９０ｎｍのような波長の光を射出することで、メラニン色素に関する知見を得ることができる。また、光源１１１が１４００ｎｍ〜２２００ｎｍの波長の光を射出することで、グルコースに関する知見を得ることができる。このような複数の波長の光は、例えば、光源１１１から時分割で射出される。
【００６２】
なお、前述の各種波長は、あくまでも一例であって、本実施形態に係る測定装置１０の光源１１１が射出する光が、上記の例に限定されるわけではない。
【００６３】
このような光源１１１としては、例えば、発光ダイオード（Ｌｉｇｈｔ  Ｅｍｉｔｔｉｎｇ  Ｄｉｏｄｅ：ＬＥＤ）や小型のレーザ等を利用可能であり、このような発光デバイスが、光源１１１として１又は複数個設けられる。
【００６４】
また、光源１１１は、後述する制御部１０３により、上記測定光の射出タイミングや射出される測定光の強度等が制御される。
【００６５】
なお、光源１１１の配設されるフレームＦの形状は特に限定されるわけではないが、光源１１１と、後述する検出部１１３との間に図７に示したような壁部を設けることで、このような壁部を、光源１１１から射出される光が検出部１１３に入射することを防止する遮光壁として用いることが可能となる。
【００６６】
○検出部について
本実施形態に係る測定装置１０が備える検出部１１３は、複数のセンサが所定の配置で規則的に配設されており、光源１１１から射出され生体Ｂを透過した測定光を、複数のセンサで検出するものである。換言すれば、本実施形態に係る検出部１１３は、いわゆるマルチタップセンサにより構成されている。図７では、このような検出部１１３の一例として、マイクロレンズアレイ（Ｍｉｃｒｏ  Ｌｅｎｓ  Ａｒｒａｙ：ＭＬＡ）を利用したセンサを示している。
【００６７】
本実施形態に係る測定装置１０が備える検出部１１３は、図７に示したように、例えば、光源１１１から射出された測定光が属する波長帯域の光を透過可能な透明基板１２１と、第１遮光体１２３と、マイクロレンズアレイ１２５と、第２遮光体１２９と、アパーチャ（絞り）１３１と、センサ１３３と、を主に備える。
【００６８】
透明基板１２１は、測定対象物である生体Ｂの一部が配設される部分である。この透明基板１２１は、測定処理で利用される波長の光を透過可能な基板を用いて形成される。光源１１１から射出され、生体Ｂの内部を透過してきた測定光は、透明基板１２１を透過すると、第１遮光体１２３によって指向性が制御される。
【００６９】
第１遮光体１２３は、透明基板１２１を透過した測定光の指向性を制御する指向性制御板として機能するものであり、後述するマイクロレンズアレイ１２５において互いに隣り合うマイクロレンズ１２７の境界部に設けられている。このような第１遮光体１２３を設けることにより、各マイクロレンズ１２７に入射する測定光の指向性を制御することが可能となり、より高精度な測定を行うことが可能となる。第１遮光体１２３を通過した測定光は、マイクロレンズアレイ１２５へと導光される。
【００７０】
マイクロレンズアレイ１２５は、図７上段に示したように、受光レンズである複数のマイクロレンズ１２７から構成されており、各マイクロレンズ１２７は、所定の基板上にｘ方向及びｙ方向に沿って格子状に配列されている。各マイクロレンズ１２７は、マイクロレンズ１２７に入射した測定光を、後述するセンサ１３３へと導光する。マイクロレンズアレイ１２５は、像面湾曲が少なく深さ方向のひずみがないレンズアレイであるため、このようなマイクロレンズアレイ１２５を用いることで、良好な測定データを得ることができる。なお、マイクロレンズアレイ１２５を構成する各マイクロレンズ１２７の被写界深度は、生体Ｂが接写距離に存在している場合であっても本実施形態に係る測定装置１０で着目する皮膚構造を包括するように（例えば、体表から数ミリ〜十数ミリの深さの範囲までがフォーカスされるように）設定される。
【００７１】
なお、本実施形態に係るマイクロレンズアレイ１２５に配設されるマイクロレンズ１２７の個数は、図７上段に示した例に限定されるわけではない。本実施形態に係るマイクロレンズアレイ１２５に配設されるマイクロレンズ１２７の個数は、撮像したい生体の大きさや、センサ１３３の大きさに応じて、自由に設定することが可能である。
【００７２】
マイクロレンズアレイ１２５に入射した測定光は、マイクロレンズ１２７により集光されて、後述するセンサ１３３へと結像されることとなる。
【００７３】
ここで、マイクロレンズアレイ１２５におけるセンサ１３３側に位置する面では、互いに隣り合うマイクロレンズ１２７の境界部に、第２遮光体１２９及びアパーチャ（絞り）１３１が設けられる。この第２遮光体１２９及びアパーチャ１３１により、マイクロレンズアレイ１２５を透過した測定光の指向性を制御することが可能となり、各マイクロレンズ１２７に入射した光を、隣接するマイクロレンズ１２７に入射した光と分離することができる。これにより、本実施形態に係る測定装置１０では、センサ１３３に集光される測定光を選択することが可能となる。
【００７４】
また、本実施形態に係る測定装置１０では、以上のような各種遮光体やアパーチャを設けることにより、各マイクロレンズ１２７に入射する光の入射角度を制限して、体内散乱によって生じる各マイクロレンズ１２７間のクロストークを防止することが可能となる。また、各マイクロレンズ１２７間のクロストークを防止することで、マイクロレンズアレイ１２５に設けられた複数のマイクロレンズ１２７のうち、一部のマイクロレンズ１２７に対応するセンサ画素から得られた信号（すなわち、測定領域のうちの局所的な位置に対応する信号）を取得することが可能となり、後述するセンサ１３３によって測定されるデータの時間分解能及び空間分解能を向上させることが可能となる。
【００７５】
センサ１３３は、図７上段に示したｘｙ平面の各位置における測定光の強度を検出する。このセンサ１３３は、光検出器（Ｐｈｏｔｏ  Ｄｅｔｅｃｔｏｒ：ＰＤ）等により受光した測定光の強度を電気信号に変換して、後述する解析部１０５へと出力する。このセンサ１３３としては、フォトダイオードや、ＣＣＤ（Ｃｈａｒｇｅ  Ｃｏｕｐｌｅｄ  Ｄｅｖｉｃｅｓ）型画像センサ、ＣＭＯＳ（Ｃｏｍｐｌｅｍｅｎｔａｒｙ  Ｍｅｔａｌ  Ｏｘｉｄｅ  Ｓｅｍｉｃｏｎｄｕｃｔｏｒ）型画像センサ、有機ＥＬを受光素子としたセンサ、ＴＦＴ（Ｔｈｉｎ  Ｆｉｌｍ  Ｔｒａｎｓｉｓｔｏｒ）型画像センサや等の２次元エリアセンサを利用することができる。また単純化したモデルとしては、例えば図８に示したように、このセンサにｘ軸方向のラインセンサ等の１次元センサを実装した測定装置１０を実現することも可能である。
【００７６】
なお、１つのマイクロレンズ１２７の下には、１又は複数の画素が配置されることとなるが、１つのマイクロレンズ１２７に対応して複数の画素が設けられる場合には、マイクロレンズ１２７と被写体との距離に応じて生じる無効な画素が存在しないように、後述する制御部１０３やソフトウェアにより制御が行われる。
【００７７】
このセンサ１３３は、後述する制御部１０３により走査タイミング等が制御され、例えば図７上段における任意の位置の検出強度を解析部１０５に出力することができる。
【００７８】
以上、図７及び図８を参照しながら、本実施形態に係る測定部１０１の構成について、詳細に説明した。
【００７９】
○測定部により測定されるデータについて
次に、図９を参照しながら、本実施形態に係る測定部１０１により測定されるデータ（測定データ）について、詳細に説明する。
【００８０】
人体の皮膚は、図１を参照しながら説明したように、皮膚表面から、表皮、真皮及び皮下脂肪の３つに区分することができる。これら３つの領域は、それぞれ光学的な特性が異なっており、互いに厚さも異なっている。また、真皮の散乱係数が２７ｍｍ−１程度であり、脂肪の散乱係数が１２．６ｍｍ−１程度であることからも明らかなように、人体は、光を極めて良く散乱させる媒質である。そのため、これらの皮膚構造の厚みの比率の個人差や経年的な変化によって、皮膚構造の総合的な光学特性は大きく異なることとなる。また、表皮に含まれるメラニン色素も、光の透過率を低下させる一因となっている。
【００８１】
光源１１１から射出され生体Ｂに入射した測定光は、生体Ｂの内部を拡散しながら、図９に示したように、略Ｕ字形状に進行して、ある位置に設けられた検出部により検出されることとなる。この際、図９に模式的に示したように、光源１１１からの位置が離れた位置に存在する検出部ほど、より深い位置まで散乱して体表へと戻ってきた測定光を検出することができる。すなわち、図９において、ｘ軸方向の位置が光源１１１から離れた位置に存在するセンサ（例えば図９において右側に位置するセンサ）ほど、より深くまで浸透した測定光を検出することができる。また、測定深度の深さは、光源１１１と着目するセンサとの離隔距離をＬとすると、Ｌ／２程度となると言われている。測定光は、この光が進んだ距離（光学的距離）の長さに応じて、その光路上に存在する各種の生体内成分により特定の波長のエネルギーが吸収され、その強度が減衰していくこととなる。
【００８２】
このように、測定光の経路が略Ｕ字形状となる反射型タイプの測定装置では、動脈での血流変動を光信号として取得する際に、探索すべき光学的な深度が大きくばらつくこととなる。従って、血管位置とこの血管位置の仮想的な深さが測定光の略Ｕ字形状の経路に含まれるか否かで、良好な信号が得られるかどうかが影響されることとなる。
【００８３】
本実施形態に係る測定装置１０では、マイクロレンズアレイを用いたセンサから得られる信号を独立して取り出すことが可能である。そのため、以下で説明するように、各マイクロレンズアレイからの出力波形を分析することで、求められる深さの信号に対応するセンサと光源との距離を正確に選択することが可能となる。
【００８４】
また、図２Ａ及び図２Ｂに示したような一般的なパルスオキシメータでは、指や耳たぶで計測を行う場合、光源と光検出器とが１０ｍｍ程度の距離にあるため、光源からの光は、指や耳たぶ内で十分に拡散される。また、光検出器は一つのフォトディテクタであるため、面積的な広がりを加算した結果のみを計測している。このようにして得られる測定データは、時間軸方向での分解能を大きく劣化させる原因となる。一方、本実施形態に係る測定装置１０のような反射型の装置の場合、測定対象との距離が数ミリまで短縮でき、また、本実施形態に係る測定装置１０では、ＭＬＡ光学系で更に微小領域の変化を検出することが可能となる。測定領域として用いられる生体部位によっては、その領域によって、良い測定データ（例えば、良い脈波形など）が得られる部位と、動脈の存在しない無効領域とが存在するが、本実施形態に係る測定装置１０では、多数のマイクロレンズから得られる信号を独立に取得して処理に利用することができるため、後述する解析部１０５と制御部１０３とが互いに連携することで、血管位置の差異を自動調整することができる。これにより、本実施形態に係る測定装置１０では、非常に時間分解能の高い脈波形を抽出することが可能となる。
【００８５】
また、例えば酸素飽和濃度など、生体の活動に伴う脈動に関する情報として着目するものの種別によっては、測定光として複数の波長を有する光を利用する場合もあるが、例えば赤外光と赤色光とでは、人の皮膚において光学的に侵入可能な深度が異なることが知られている。例えば、７００〜１０００ｎｍの波長の光では光散乱が主体となるが、可視光や波長１０００ｎｍ以上の光では吸収が主体となるため、光が生体内の深くまで侵入することが困難となる。このため、例えば波長６６０ｎｍの赤色光は、波長９４０ｎｍの赤外光（近赤外光）よりも浅い部分を介してセンサに光が戻る性質がある。その結果、本実施形態に係る測定装置１０においても、ＭＬＡにおける動脈からの最適距離に、波長依存性が生じる場合がある。本実施形態に係る測定装置１０では、後述する解析部１０５によって、この距離の差をＭＬＡからの信号の振幅に基づいて解析することで、波長により異なる深度調整のためのＭＬＡ位置の最適化を、波長毎に可能としている。
【００８６】
［制御部１０３について］
再び図６に戻って、本実施形態に係る測定装置１０が備える制御部１０３について説明する。
制御部１０３は、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）等により実現される。制御部１０３は、測定部１０１に設けられた光源１１１やセンサ１３３等の駆動制御を行うことにより、測定部１０１における生体Ｂの測定処理全般を統括する。より詳細には、制御部１０３は、所定の同期信号等に基づいて、センサ１３３の走査タイミングや、情報を取得するセンサ１３３の選択等といったセンサの駆動制御を行う。また、制御部１０３は、光源１１１に対しても、測定光の射出タイミングや強度に関する駆動制御を行う。
【００８７】
制御部１０３が以上のような駆動制御を行うことで、測定部１０１の光源１１１は、時分割で異なる波長の測定光を射出することが可能となるとともに、センサ１３３上の任意の位置の測定データを時分割で取得することが可能となる。
【００８８】
制御部１０３により駆動制御された測定部１０１によって測定された測定データは、後述する解析部１０５へと出力されて、測定データの解析処理が実施される。
【００８９】
ここで、制御部１０３は、測定部１０１の制御を行うにあたり、後述する記憶部１０７に記録されている各種のプログラムやパラメータやデータベース等を参照することが可能である。
【００９０】
［解析部１０５について］
本実施形態に係る測定装置１０が備える解析部１０５は、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ等により実現される。解析部１０５は、測定部１０１の検出部１１３により検出された検出結果を利用し、検出された測定光の光量の時間変化に基づいて、生体の活動に伴う脈動に関する情報を測定するための測定位置を、測定領域の中から特定する解析処理を行う。
【００９１】
図２Ａ及び図２Ｂに示したような一般的な測定機器では、指先や耳たぶなどの各小静脈で観測される血流による時間変化群を、ひとつのセンサで加算した結果のみを取り出して解析していた。そのため、図２Ａ及び図２Ｂに示したような一般的な測定機器では、以下のような欠点が生じることとなる。
【００９２】
・無効となる静脈血が存在する領域や動脈の分布の少ない領域の信号加算に伴う信号雑音比（Ｓ／Ｎ）の劣化
・時間的に位相の異なる波形が加算される（合成される）ことによる信号波形の劣化
・広範囲のデータの積分加算による波形の時間分解能の劣化
【００９３】
一方で、本実施形態に係る測定装置１０では、測定部１０１の光学系としてマイクロレンズアレイ光学系を採用することにより、各マイクロレンズに対応するセンサ画素から、それぞれ測定部位の異なる様々な波形の信号を選択して取得することができる。そこで、本実施形態に係る解析部１０５では、各マイクロレンズと光源との離隔距離が、皮下の反射散乱によって戻ってくる深さ（すなわち、測定深度）と大きな相関があることを利用して、以下のような処理を実施する。
【００９４】
・微小領域での血流変化によって得られる適切な脈波形データを取得
・個人によって異なる血管位置や深さに応じて、適切なデータ取得位置を調整
・波長によって異なる戻り光の深さの変化の影響を個別解析し、正確な位置からの脈動の輝度変化を取得
【００９５】
また、本実施形態に係る解析部１０５では、このような微小領域での波形の結果を合成することにより、より良い脈波形の取得が可能となる。
【００９６】
脈拍などといった周期成分を求める場合であれば、本実施形態に係る解析部１０５は、脈波形の最大振幅が得られるマイクロレンズの位置からの信号を選択する方法がある。一方、この方法では波形が適切ではない場合もあるので、解析部１０５は、最適波形を求める用途では、以下の方法を利用することができる。
【００９７】
本実施形態に係る測定装置１０では、上述のように、各マイクロレンズに対応する画素から取得したデータ（測定領域中の微小領域でのデータ）を別個に取得して利用することができる。そこで、解析部１０５は、各マイクロレンズに対応する画素（又は画素群）から取得した微小領域でのデータの時間変化に着目し、生体の活動に伴う脈動に関する情報（以下、単に脈動情報とも称する。）を測定するための測定位置を測定領域の中から特定する。
【００９８】
すなわち、本実施形態に係る解析部１０５では、微小領域でのデータの時間変化に着目することで、着目する微小領域での脈波形データを算出することができる。そこで、解析部１０５は、事前に予め測定された脈波形（例えば、他の測定機器等により測定されたデータ）を利用して、この予め測定された脈波形と、それぞれの微小領域での脈波形データとの類似度合いを算出する。ここで、解析部１０５は、２つの波形の類似度合いを表す評価値として、相互相関係数を算出してもよいし、差分絶対値和（Ｓｕｍ  ｏｆ  Ａｂｓｏｌｕｔｅ  Ｄｉｆｆｅｒｅｎｃｅ：ＳＡＤ）や差分２乗和（Ｓｕｍ  ｏｆ  Ｓｑｕａｒｅｄ  Ｄｉｆｆｅｒｅｎｃｅ：ＳＳＤ）等を算出してもよい。解析部１０５は、算出した類似度合いが所定の閾値以上であった場合（すなわち、２つの波形が類似していると判定できる場合）に、そのような類似度合いを与える微小領域を、脈動情報を測定するために適した測定位置として特定することができる。
【００９９】
これにより、解析部１０５は、適切な脈波形の振幅を与える微小領域の集合や、適切な脈波形形状を与える微小領域の集合など、着目する脈動情報の種別に応じて、測定に適した測定位置を自動的に特定することができる。適切な脈波形の振幅を与える微小領域の集合は、例えば、測定装置１０による測定結果に基づいて生体の脈拍数を特定する場合に有用な測定位置として用いることが可能である。また、適切な脈波形形状を与える微小領域の集合は、例えば、脈拍数以外の脈動情報に関する知見を得る際に有用な測定位置として用いることが可能である。
【０１００】
本実施形態に係る解析部１０５は、このようにして得られた脈動情報を得るために適した測定位置に関する情報に基づいて、測定部１０１から出力される測定データを処理に利用する際のデータの選択を行うことができる。また、本実施形態に係る解析部１０５は、このようにして得られた測定位置に関する情報を制御部１０３にフィードバックしてもよい。これにより、制御部１０３では、測定部１０１の制御を効率良く行うことが可能となる。
【０１０１】
また、先だって説明したように、光源として用いる光の波長に応じて、脈動情報を得るために適した測定位置は変化する可能性があるが、本実施形態に係る解析部１０５では、上記のような解析処理を光源として用いられる光の波長毎に実施するため、脈動情報を得るために適した測定位置を測定光の波長毎に特定することが可能となる。
【０１０２】
このように、本実施形態に係る解析部１０５による測定位置の解析処理により、本実施形態に係る測定装置１０では、脈動情報を得るための測定位置を自動的に調整することが可能となる。その結果、例えば図１０Ａに示したような、指を測定領域とする場合や、図１０Ｂに示したように、動脈が広範な測定領域の一部にしか存在していないような場所を測定箇所とした場合であっても、有効な測定領域を自動的に特定することができる。その結果、本実施形態に係る測定装置１０では、生体を測定装置１０の特定の位置に接触させなかった場合であっても、適切な測定を行うことが可能となる。本実施形態に係る測定装置１０のこのような機能を、例えば携帯電話や携帯用タブレット等といった各種携帯機器に実装することで、使用者が携帯機器を持つだけで脈動情報の測定を行うことも可能となる。
【０１０３】
本実施形態に係る解析部１０５では、以上のようにして特定した測定位置でのデータを利用して解析処理を行うことで、各種の脈動情報（例えば、脈波形そのものや、酸素飽和度や、これらの情報から二次的に算出可能な情報等）を更に算出することができる。
【０１０４】
すなわち、解析部１０５は、適切な測定位置から得られた測定データの時間変化に着目し、センサにより検知された光の光量や輝度の時間推移を追うことで、動脈の拍動に由来する脈波形を得ることができる。
【０１０５】
また、解析部１０５は、上記式１８〜式１９ｃに示した演算式を利用することで、適切な測定位置から得られた測定データに基づいて、動脈血の酸素飽和度ＳｐＯ２を算出することができる。
【０１０６】
また、解析部１０５は、解析処理によって得られた脈波形データを更に利用して、ＡＩ値（Ａｕｇｍｅｎｔａｔｉｏｎ  Ｉｎｄｅｘ：脈波増大係数）や脈波伝播速度（ＰＷＶ）等といった２次的な情報を更に算出することが可能である。これらの２次的な情報を算出する方法については特に限定されるわけではなく、公知の方法を利用可能である。
【０１０７】
更に、解析部１０５は、適切な測定位置から得られた測定データを利用して、公知の方法により多変量解析を行うことで、メラニン色素の成分量や、グルコース、アルブミン、コレステロール等といった血中成分の成分量などのように、広範な物質濃度の測定を行うことも可能である。
【０１０８】
本実施形態に係る解析部１０５は、例えば図１１に示したように、光源からほぼ等距離にあるマイクロレンズ群に対応する画素から得られた検出信号を加算して、信号強度を増加させてもよい。図１１に示したように、光源からの離隔距離がほぼ等距離となる位置に対応するセンサの水平画素は、互いに同位相となる測定データを検知している。そのため、図１１に示したように、これらの水平ラインから得られた検出信号を積算することによって、信号強度を増加させることが可能となる。
【０１０９】
複数の信号を積算する（換言すれば、複数の信号を合成する）場合には、その時間軸にズレが生じていると、互いの信号が打ち消し合ってしまう場合があり、時間分解能が劣化する恐れがある。そのため、本実施形態に係る解析部１０５では、複数の信号を合成する場合には、図１２に示したように、信号波形（例えば、脈波形等）のトップ又はボトムの位置が一致するように時間軸を補正した後、信号の積算処理を実施する。
【０１１０】
ここで、図１１に示した垂直走査方向（すなわち、垂直同期方向）を、複数の光源の配置方向に対して垂直となるようにセンサを配設することで、ＣＣＤ以外にＣＭＯＳセンサでも時間軸の一致を図ることが可能となる。なお、垂直走査方向が光源の配置方向と平行となる場合には、時間軸の異なる垂直方向の画素が加算されることとなるため、時間分解能が劣化してしまう。そのため、垂直走査方向は光源の配置方向に対して垂直となるようにすることが望ましい。
【０１１１】
また、本実施形態に係る解析部１０５では、互いに同位相にある信号だけでなく、位相の異なる信号同士（例えば、図１１において、異なる垂直走査方向に位置する画素から得られた信号同士）を合成して、各種の解析処理を実施することも可能である。解析部１０５は、この場合においても、図１２で説明したような方法で時間軸を適切に補正して、信号が互いに打ち消し合わないようにした上で合成処理を実施する。
【０１１２】
以上、図１０〜図１２を参照しながら、本実施形態に係る解析部１０５について、詳細に説明した。
【０１１３】
［記憶部１０７について］
再び図６に戻って、本実施形態に係る測定装置１０が備える記憶部１０７について説明する。
記憶部１０７は、本実施形態に係る測定装置１０に設けられたＲＡＭやストレージ装置等により実現される。記憶部１０７には、解析部１０５における解析処理に用いられる光吸収スペクトルのデータや、各種のデータベースやルックアップテーブル等が格納されている。また、記憶部１０７には、本実施形態に係る測定部１０１により測定された測定データや、本実施形態に係る制御部１０３や解析部１０５が実施する処理に用いられる各種のプログラムやパラメータやデータ等が記録されていてもよい。また、記憶部１０７には、これらのデータ以外にも、測定装置１０が、何らかの処理を行う際に保存する必要が生じた様々なパラメータや処理の途中経過等を適宜記憶することが可能である。この記憶部１０７は、測定部１０１、制御部１０３、解析部１０５等の各処理部が、自由にアクセスし、データを書き込んだり読み出したりすることができる。
【０１１４】
以上、図６〜図１２を参照しながら、本実施形態に係る測定装置１０の構成について、詳細に説明した。
【０１１５】
以上説明したように、本実施形態に係る測定装置１０によれば、指や耳たぶを用いた反射型や透過型のパルスオキシメータでは捉えることが困難であった微小信号を検出することが可能となり、これらの微小信号に基づく正確で時間分解能の高い脈波形を検出することが可能となる。
【０１１６】
また、人への負担が少ない片側だけでの計測が可能な反射散乱型の照明方式により、長時間の計測を低負荷で実現できる。また、本実施形態に係る測定装置１０では、指等の測定部位を厳密に固定するメカニズムが不要となり、自動的にイメージセンサの有効な領域を用いた自動調整が可能となる。これらの効果により、血液中の成分分析をより広範囲な物質濃度の計測へと利用できるようになり、例えば、グルコース、アルブミン、コレステロールなどの微小成分への対応が可能となる。
【０１１７】
なお、本実施形態に係る制御部１０３及び解析部１０５は、本実施形態に係る測定装置１０の一部であってもよいし、測定装置１０に接続されているコンピュータ等の外部機器に実現されていてもよい。また、測定部１０１によって生成される測定データがリムーバブル記憶媒体等に格納され、この記憶媒体が測定装置１０から取り外されて、解析部１０５を有する他の装置に接続されることで、測定データが解析されてもよい。
【０１１８】
以上、本実施形態に係る測定装置１０の機能の一例を示した。上記の各構成要素は、汎用的な部材や回路を用いて構成されていてもよいし、各構成要素の機能に特化したハードウェアにより構成されていてもよい。また、各構成要素の機能を、ＣＰＵ等が全て行ってもよい。従って、本実施形態を実施する時々の技術レベルに応じて、適宜、利用する構成を変更することが可能である。
【０１１９】
なお、上述のような本実施形態に係る測定装置の各機能を実現するためのコンピュータプログラムを作製し、パーソナルコンピュータ等に実装することが可能である。また、このようなコンピュータプログラムが格納された、コンピュータで読み取り可能な記録媒体も提供することができる。記録媒体は、例えば、磁気ディスク、光ディスク、光磁気ディスク、フラッシュメモリなどである。また、上記のコンピュータプログラムは、記録媒体を用いずに、例えばネットワークを介して配信してもよい。
【０１２０】
（ハードウェア構成について）
次に、図１３を参照しながら、本開示の実施形態に係る測定装置１０のハードウェア構成について、詳細に説明する。図１３は、本開示の実施形態に係る測定装置１０のハードウェア構成を説明するためのブロック図である。
【０１２１】
測定装置１０は、主に、ＣＰＵ９０１と、ＲＯＭ９０３と、ＲＡＭ９０５と、を備える。また、測定装置１０は、更に、ホストバス９０７、ブリッジ９０９、外部バス９１１、インターフェース９１３、センサ９１４、入力装置９１５、出力装置９１７、ストレージ装置９１９、ドライブ９２１、接続ポート９２３および通信装置９２５を備える。
【０１２２】
ＣＰＵ９０１は、演算処理装置および制御装置として機能し、ＲＯＭ９０３、ＲＡＭ９０５、ストレージ装置９１９、またはリムーバブル記録媒体９２７に記録された各種プログラムに従って、測定装置１０内の動作全般またはその一部を制御する。ＲＯＭ９０３は、ＣＰＵ９０１が使用するプログラムや演算パラメータ等を記憶する。ＲＡＭ９０５は、ＣＰＵ９０１が使用するプログラムや、プログラムの実行において適宜変化するパラメータ等を一次記憶する。これらはＣＰＵバス等の内部バスにより構成されるホストバス９０７により相互に接続されている。
【０１２３】
ホストバス９０７は、ブリッジ９０９を介して、ＰＣＩ（Ｐｅｒｉｐｈｅｒａｌ  Ｃｏｍｐｏｎｅｎｔ  Ｉｎｔｅｒｃｏｎｎｅｃｔ／Ｉｎｔｅｒｆａｃｅ）バスなどの外部バス９１１に接続されている。
【０１２４】
センサ９１４は、例えば、ユーザに固有の生体情報、または、かかる生体情報を取得するために用いられる各種情報を検出する検出手段である。このセンサ９１４として、例えば、ＣＣＤ（Ｃｈａｒｇｅ  Ｃｏｕｐｌｅｄ  Ｄｅｖｉｃｅ）やＣＭＯＳ（Ｃｏｍｐｌｅｍｅｎｔａｒｙ  Ｍｅｔａｌ  Ｏｘｉｄｅ  Ｓｅｍｉｃｏｎｄｕｃｔｏｒ）等の各種の撮像素子を挙げることができる。また、センサ９１４は、生体部位を撮像するために用いられるレンズ等の光学系や光源等を更に有していてもよい。また、センサ９１４は、音声等を取得するためのマイクロフォン等であってもよい。なお、センサ９１４は、上述のもの以外にも、温度計、照度計、湿度計、速度計、加速度計などの様々な測定機器を備えていてもよい。
【０１２５】
入力装置９１５は、例えば、マウス、キーボード、タッチパネル、ボタン、スイッチおよびレバーなどユーザが操作する操作手段である。また、入力装置９１５は、例えば、赤外線やその他の電波を利用したリモートコントロール手段（いわゆる、リモコン）であってもよいし、測定装置１０の操作に対応した携帯電話やＰＤＡ等の外部接続機器９２９であってもよい。さらに、入力装置９１５は、例えば、上記の操作手段を用いてユーザにより入力された情報に基づいて入力信号を生成し、ＣＰＵ９０１に出力する入力制御回路などから構成されている。測定装置１０のユーザは、この入力装置９１５を操作することにより、測定装置１０に対して各種のデータを入力したり処理動作を指示したりすることができる。
【０１２６】
出力装置９１７は、取得した情報をユーザに対して視覚的または聴覚的に通知することが可能な装置で構成される。このような装置として、ＣＲＴディスプレイ装置、液晶ディスプレイ装置、プラズマディスプレイ装置、ＥＬディスプレイ装置およびランプなどの表示装置や、スピーカおよびヘッドホンなどの音声出力装置や、プリンタ装置、携帯電話、ファクシミリなどがある。出力装置９１７は、例えば、測定装置１０が行った各種処理により得られた結果を出力する。具体的には、表示装置は、測定装置１０が行った各種処理により得られた結果を、テキストまたはイメージで表示する。他方、音声出力装置は、再生された音声データや音響データ等からなるオーディオ信号をアナログ信号に変換して出力する。
【０１２７】
ストレージ装置９１９は、測定装置１０の記憶部の一例として構成されたデータ格納用の装置である。ストレージ装置９１９は、例えば、ＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）等の磁気記憶部デバイス、半導体記憶デバイス、光記憶デバイス、または光磁気記憶デバイス等により構成される。このストレージ装置９１９は、ＣＰＵ９０１が実行するプログラムや各種データ、および外部から取得した各種データなどを格納する。
【０１２８】
ドライブ９２１は、記録媒体用リーダライタであり、測定装置１０に内蔵、あるいは外付けされる。ドライブ９２１は、装着されている磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリ等のリムーバブル記録媒体９２７に記録されている情報を読み出して、ＲＡＭ９０５に出力する。また、ドライブ９２１は、装着されている磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリ等のリムーバブル記録媒体９２７に記録を書き込むことも可能である。リムーバブル記録媒体９２７は、例えば、ＤＶＤメディア、ＨＤ−ＤＶＤメディア、Ｂｌｕ−ｒａｙメディア等である。また、リムーバブル記録媒体９２７は、コンパクトフラッシュ（登録商標）（ＣｏｍｐａｃｔＦｌａｓｈ：ＣＦ）、フラッシュメモリ、または、ＳＤメモリカード（Ｓｅｃｕｒｅ  Ｄｉｇｉｔａｌ  ｍｅｍｏｒｙ  ｃａｒｄ）等であってもよい。また、リムーバブル記録媒体９２７は、例えば、非接触型ＩＣチップを搭載したＩＣカード（Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ  ｃａｒｄ）または電子機器等であってもよい。
【０１２９】
接続ポート９２３は、機器を測定装置１０に直接接続するためのポートである。接続ポート９２３の一例として、ＵＳＢ（Ｕｎｉｖｅｒｓａｌ  Ｓｅｒｉａｌ  Ｂｕｓ）ポート、ＩＥＥＥ１３９４ポート、ＳＣＳＩ（Ｓｍａｌｌ  Ｃｏｍｐｕｔｅｒ  Ｓｙｓｔｅｍ  Ｉｎｔｅｒｆａｃｅ）ポート等がある。接続ポート９２３の別の例として、ＲＳ−２３２Ｃポート、光オーディオ端子、ＨＤＭＩ（Ｈｉｇｈ−Ｄｅｆｉｎｉｔｉｏｎ  Ｍｕｌｔｉｍｅｄｉａ  Ｉｎｔｅｒｆａｃｅ）ポート等がある。この接続ポート９２３に外部接続機器９２９を接続することで、測定装置１０は、外部接続機器９２９から直接各種データを取得したり、外部接続機器９２９に各種データを提供したりする。
【０１３０】
通信装置９２５は、例えば、通信網９３１に接続するための通信デバイス等で構成された通信インターフェースである。通信装置９２５は、例えば、有線または無線ＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）、Ｂｌｕｅｔｏｏｔｈ（登録商標）、またはＷＵＳＢ（Ｗｉｒｅｌｅｓｓ  ＵＳＢ）用の通信カード等である。また、通信装置９２５は、光通信用のルータ、ＡＤＳＬ（Ａｓｙｍｍｅｔｒｉｃ  Ｄｉｇｉｔａｌ  Ｓｕｂｓｃｒｉｂｅｒ  Ｌｉｎｅ）用のルータ、または、各種通信用のモデム等であってもよい。この通信装置９２５は、例えば、インターネットや他の通信機器との間で、例えばＴＣＰ／ＩＰ等の所定のプロトコルに則して信号等を送受信することができる。また、通信装置９２５に接続される通信網９３１は、有線または無線によって接続されたネットワーク等により構成され、例えば、インターネット、家庭内ＬＡＮ、赤外線通信、ラジオ波通信または衛星通信等であってもよい。
【０１３１】
以上、本開示の実施形態に係る測定装置１０の機能を実現可能なハードウェア構成の一例を示した。上記の各構成要素は、汎用的な部材を用いて構成されていてもよいし、各構成要素の機能に特化したハードウェアにより構成されていてもよい。従って、本実施形態を実施する時々の技術レベルに応じて、適宜、利用するハードウェア構成を変更することが可能である。
【０１３２】
（まとめ）
以上説明したように、本開示の実施形態によれば、生体の活動に伴う脈動に関する情報を測定するために適した測定位置を特定することが可能であり、個人差に基づく測定位置の変化に柔軟に対応することが可能となる。
【０１３３】
また、マイクロレンズアレイ光学系を用いたセンサから得られた局所的な測定データを利用することで、一般的な反射型タイプや透過型タイプの測定装置よりも時間分解能の高い脈波形を得ることが可能となる。この際、マイクロレンズアレイを用いたセンサから得られた各信号を時間補正して合成することにより、たとえ指先であっても十分に正確な脈波形を測定することができる。
【０１３４】
また、本開示の実施形態に係る測定装置では、測定に適した位置を自動的に検出することが可能となるため、位置合わせのための生体の固定部材が不要となり、使用者の利便性を高めることができる。
【０１３５】
また、本開示の実施形態に係る測定装置では、空間分解能及び時間分解能に優れた脈波形データを測定することができるため、酸素飽和度ＳｐＯ２のみならず、グルコース、アルブミン、コレステロールなどの微小な血液成分の分析にも対応できるようになる。
【０１３６】
更に、本開示の実施形態に係る測定装置を複数利用することで、脈波形データを生体の複数の測定位置で検出することも可能となり、中心血圧等の測定にも利用することができる。
【０１３７】
以上、添付図面を参照しながら本開示の好適な実施形態について詳細に説明したが、本開示の技術的範囲はかかる例に限定されない。本開示の技術分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本開示の技術的範囲に属するものと了解される。
【０１３８】
なお、以下のような構成も本開示の技術的範囲に属する。
（１）
生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出する光源と、
複数のセンサが所定の配置で規則的に配設されており、前記光源から射出され前記生体を透過した前記測定光を当該複数のセンサで検出する検出部と、
前記検出部により検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行う解析部と、
を備える、測定装置。
（２）
前記検出部は、複数のレンズが格子状に規則的に配設されたマイクロレンズアレイを利用したセンサにより、前記生体を透過した測定光を検出する、（１）に記載の測定装置。
（３）
前記マイクロレンズアレイと前記複数のセンサとの間には、検出する前記測定光のクロストークを防止する遮光体が設けられる、（２）に記載の測定装置。
（４）
前記解析部は、前記測定光の光量の時間変化に基づき、前記生体の活動に伴う脈動に関する情報として、前記測定領域に対応する前記生体内部の動脈の拍動に由来する脈波形を更に算出する、（１）〜（３）の何れか１つに記載の測定装置。
（５）
前記解析部は、予め測定された前記生体における前記脈波形と、前記測定光の光量の時間変化に基づき算出された前記脈波形と、の類似度合いに応じて、前記測定位置を特定する、（４）に記載の測定装置。
（６）
前記解析部は、前記測定光の光量の時間変化に基づき算出された前記脈波形の最大振幅を与える位置を、前記測定位置として特定する、（４）に記載の測定装置。
（７）
前記解析部は、前記生体の活動に伴う脈動に関する情報として前記脈波形の振幅を測定するための前記測定位置を特定する、（５）に記載の測定装置。
（８）
前記解析部は、前記生体の活動に伴う脈動に関する情報として前記脈波形の形状を測定するための前記測定位置を特定する、（５）に記載の測定装置。
（９）
前記解析部は、前記光源から射出される前記測定光の波長毎に、前記解析処理を実施する、（１）〜（８）の何れか１つに記載の測定装置。
（１０）
前記解析部は、算出した前記脈波形を利用して、前記動脈における酸素飽和度を更に算出する、（４）〜（９）の何れか１つに記載の測定装置。
（１１）
前記解析部は、前記複数のセンサにおける前記測定光の検出位置と前記光源との離隔距離に応じて、前記検出結果の時間軸における位相を補正する、（１）〜（１０）の何れか１つに記載の測定装置。
（１２）
前記解析部は、前記複数のセンサにおいて前記光源からの離隔距離が同一である検出位置で得られた前記検出結果を積算する、（１）〜（１１）の何れか１つに記載の測定装置。
（１３）
前記解析部は、互いに異なる時間に対応する複数の前記検出結果について、当該時間を補正したうえで合成する、（１）〜（１２）の何れか１つに記載の測定装置。
（１４）
生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出することと、
複数のセンサが所定の配置で規則的に配設された検出部により、前記光源から射出され前記生体を透過した前記測定光を検出することと、
検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行うことと、
を含む、測定方法。
（１５）
生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出する光源と、複数のセンサが所定の配置で規則的に配設されており、前記光源から射出され前記生体を透過した前記測定光を当該複数のセンサで検出する検出部と、を備える測定機器と通信可能なコンピュータに、
前記検出部により検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行う解析機能
を実現させるためのプログラム。
（１６）
生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出する光源と、複数のセンサが所定の配置で規則的に配設されており、前記光源から射出され前記生体を透過した前記測定光を当該複数のセンサで検出する検出部と、を備える測定機器と通信可能なコンピュータに、
前記検出部により検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行う解析機能
を実現させるためのプログラムが記録された記録媒体。
（１７）
生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出する光源と、
複数のレンズが格子状に規則的に配設されたマイクロレンズアレイを利用したセンサにより、前記光源から射出され前記生体を透過した前記測定光を検出する検出部と、
前記検出部により検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行うとともに、記測定光の光量の時間変化に基づき、前記生体の活動に伴う脈動に関する情報として、前記測定領域に対応する前記生体内部の動脈の拍動に由来する脈波形を算出する解析部と、
を備える、測定装置。
（１８）
生体の少なくとも一部からなる測定領域に対して、所定の波長帯域に属する少なくとも１種類の測定光を射出することと、
複数のレンズが格子状に規則的に配設されたマイクロレンズアレイを利用したセンサにより、前記光源から射出され前記生体を透過した前記測定光を検出することと、
検出された検出結果を利用し、検出された前記測定光の光量の時間変化に基づいて、前記生体の活動に伴う脈動に関する情報を測定するための測定位置を前記測定領域の中から特定する解析処理を行うとともに、前記測定光の光量の時間変化に基づき、前記生体の活動に伴う脈動に関する情報として、前記測定領域に対応する前記生体内部の動脈の拍動に由来する脈波形を算出することと、
を含む、測定方法。
【符号の説明】
【０１３９】
１０    測定装置
１０１    測定部
１０３    制御部
１０５    解析部
１０７    記憶部
１１１    光源
１２１    透明基板
１２３    第１遮光体
１２５    マイクロレンズアレイ
１２７    マイクロレンズ
１２９    第２遮光体
１３１    アパーチャ（絞り）
１３３    センサ
【図１】
0005948836.tif 000008
【図２Ａ】
0005948836.tif 000009
【図２Ｂ】
0005948836.tif 000010
【図３】
0005948836.tif 000011
【図４】
0005948836.tif 000012
【図５】
0005948836.tif 000013
【図６】
0005948836.tif 000014
【図７】
0005948836.tif 000015
【図８】
0005948836.tif 000016
【図９】
0005948836.tif 000017
【図１０Ａ】
0005948836.tif 000018
【図１０Ｂ】
0005948836.tif 000019
【図１１】
0005948836.tif 000020
【図１２】
0005948836.tif 000021
【図１３】
0005948836.tif 000022