CN101836104A - 活体测量装置 - Google Patents

活体测量装置 Download PDF

Info

Publication number
CN101836104A
CN101836104A CN200880112934A CN200880112934A CN101836104A CN 101836104 A CN101836104 A CN 101836104A CN 200880112934 A CN200880112934 A CN 200880112934A CN 200880112934 A CN200880112934 A CN 200880112934A CN 101836104 A CN101836104 A CN 101836104A
Authority
CN
China
Prior art keywords
light
absorption coefficient
measured position
wavelength
medium
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
CN200880112934A
Other languages
English (en)
Other versions
CN101836104B (zh
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.)
Hamamatsu Photonics KK
Original Assignee
Hamamatsu Photonics KK
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 Hamamatsu Photonics KK filed Critical Hamamatsu Photonics KK
Publication of CN101836104A publication Critical patent/CN101836104A/zh
Application granted granted Critical
Publication of CN101836104B publication Critical patent/CN101836104B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47Scattering, i.e. diffuse reflection
    • G01N21/4795Scattering, i.e. diffuse reflection spatially resolved investigating of object in scattering medium
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0082Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes
    • A61B5/0091Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes for mammography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/43Detecting, measuring or recording for evaluating the reproductive systems
    • A61B5/4306Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
    • A61B5/4312Breast evaluation or disorder diagnosis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00Features of devices classified in G01N21/00
    • G01N2201/08Optical fibres; light guides
    • G01N2201/0846Fibre interface with sample, e.g. for spatial resolution

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Engineering & Computer Science (AREA)
  • Gynecology & Obstetrics (AREA)
  • Reproductive Health (AREA)
  • Optics & Photonics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明涉及活体测量装置,通过向受检者被测量部位照射光并检测漫射光而取得被测量部位内部信息,其具备保持光透过性介质的容器、具有固定于容器的多个光出射端并向浸于介质的被测量部位照射波长不同的第1和第2光的光照射部、具有固定于容器的多个光检测端并检测来自被测量部位的漫射光的光检测部、基于来自光检测部的输出信号计算内部信息的运算部;第1光的波长为被测量部位吸收系数与介质吸收系数的平均值实质相等的波长,第2光的波长为被测量部位吸收系数大于介质吸收系数平均值的波长;运算部基于与第1光的漫射光相关的输出信号计算内部信息,并基于与第2光的漫射光相关的输出信号计算被测量部位与介质之间的边界信息。

Description

活体测量装置
技术领域
本发明涉及活体测量装置。
背景技术
作为对头部或乳房等活体的内部信息进行非侵入的测量的装置,提出了利用活体的光吸收特性而得到内部信息的装置(例如,参照专利文献1)。在这样的测量装置中,从规定的照射位置对作为测量对象的活体的部位照射光,并在规定的检测位置对在该部位的内部散射并传播的光进行检测,从而由其强度及时间波形等测定结果,能够得到该部位的内部信息,即与存在于该部位的内部的肿瘤等光吸收体相关的信息。
此外,专利文献1中记载了:使介质(以下,称为光学干涉(interface)材料)介于光的照射位置以及检测位置与被测量部位之间,通过防止在被测量部位的表面的光的反射及散射等从而可以提高测量精度,其中该介质在照射光的波长(测量波长)下具有与被测量部位大致相同的光学特性。
专利文献1:日本专利特许第3771364号公报
发明内容
在使如上所述的光学干涉材料介于光的照射位置与被测量部位之间的情况下,由于在测量波长下被测量部位与介质具有大致相同的光学特性,因而能够高精度地测量光的照射位置或检测位置与存在于被测量部位的内部的光吸收体之间的相对位置。但是,由于难以区别光学干涉材料和被测量部位,因而难以测量该光吸收体相对于被测量部位整体的存在位置。因此,难以进行与由X射线检查装置或MRI、超声波诊断装置等其它医用图像诊断装置得到的测量结果的比较研究讨论,或者难以进行用于取得肿瘤的组织标本的活组织检查(biopsy)(穿刺针)。
本发明正是鉴于以上的问题,其目的在于,在通过向受检者的被测量部位照射光并检测漫射光从而取得被测量部位的内部信息的活体测量装置中,除了被测量部位的内部信息之外,还取得肿瘤等光吸收体相对于被测量部位整体的存在位置。
解决问题的方法
为了解决上述问题,本发明的活体测量装置的特征在于,是通过向受检者的被测量部位照射光并检测漫射光从而取得被测量部位的内部信息的活体测量装置,具备:保持光透过性介质的容器、具有被固定于容器的多个光出射端并向浸于介质的被测量部位照射波长彼此不同的第1和第2光的光照射部、具有被固定于容器的多个光检测端并检测来自被测量部位的漫射光的光检测部、以及基于来自光检测部的输出信号计算出内部信息的运算部;第1光的波长为被测量部位的吸收系数与介质的吸收系数的平均值实际上相等的波长,第2光的波长为被测量部位的吸收系数大于介质的吸收系数的平均值的波长;运算部基于与第1光的漫射光相关的输出信号而计算出内部信息,并基于与第2光的漫射光相关的输出信号而计算出被测量部位与介质之间的边界信息。
在上述的活体测量装置中,能够进行以下的两种测量,即内部信息测量以及轮廓测量。所谓内部信息测量,是指光照射部朝向浸于介质的被测量部位照射第1光,光检测单元检测其漫射光,运算部基于该检测结果计算出内部信息的测量。所谓轮廓测量,是指光照射部朝向浸于介质的被测量部位照射第2光,光检测部检测其漫射光,运算部基于该检测结果计算出被测量部位与介质之间的边界信息(即被测量部位的轮廓信息)的测量。
在第1光的波长(设为λ1)下,被测量部位的吸收系数(设为μab)与介质的吸收系数(设为μai)实际上相等(即μab(λ1)=μai(λ1))。因此,在使用第1光的内部信息测量中,可以防止被测量部位的表面上的光的反射及散射等,从而能够高精度地测量以光的照射位置以及检测位置为基准的肿瘤等光吸收体的位置以及大小。另一方面,在第2光的波长(设为λ2)下,被测量部位的吸收系数(μab)大于介质的吸收系数(μai)(即μab(λ2)>μai(λ2))。因此,在使用第2光的轮廓测量中,能够高精度地测量以光的照射位置以及检测位置为基准的被测量部位的轮廓信息。并且,通过对内部信息测量的结果以及轮廓测量的结果进行综合,从而得到光吸收体相对于被测量部位整体的存在位置。
这样,根据上述的活体测量装置,在通过向受检者的被测量部位照射光并检测漫射光从而取得被测量部位的内部信息的装置中,不仅能够高精度地测量肿瘤等光吸收体的有无以及大小,而且能够高精度地测量光吸收体相对于被测量部位整体的存在位置。
此外,活体测量装置也可以为:介质对第1光的吸收系数大于介质对第2光的吸收系数。在被测量部位的吸收系数μab(λ1)以及μab(λ2)彼此相等的情况下,通过这样使介质的吸收系数μai(λ1)>μai(λ2),从而能够得到μab(λ1)=μai(λ1)且μab(λ2)>μai(λ2)的关系。或者是,即使在被测量部位的吸收系数μab(λ1)大于μab(λ2)的情况下,通过以介质的吸收系数的差(μai(λ1)-μai(λ2))大于被测量部位的吸收系数的差(μab(λ1)-μab(λ2))的方式调整介质的光吸收特性或设定波长λ1、λ2,也能得到μab(λ1)=μai(λ1)且μab(λ2)>μai(λ2)的关系。因此,能够适宜地进行内部信息测量以及轮廓测量。
此外,活体测量装置也可以为:介质对第1光的吸收系数与介质对第2光的吸收系数实际上相等。在被测量部位的吸收系数μab(λ1)小于μab(λ2)的情况下,也可以以这样的使介质的吸收系数为μai(λ1)=μai(λ2)的方式调整介质的光吸收特性或设定波长λ1、λ2。由此,能够得到μab(λ1)=μai(λ1)且μab(λ2)>μai(λ2)的关系。因此,能够适宜地进行内部信息测量以及轮廓测量。
发明的效果
根据本发明,在通过向受检者的被测量部位照射光并检测漫射光从而取得被测量部位的内部信息的活体测量装置中,除了被测量部位的内部信息之外,还能够取得肿瘤等光吸收体相对于被测量部位整体的存在位置。
附图说明
图1为示意活体测量装置的一个实施方式的系统构成图。
图2为示意图1所示的活体测量装置的使用状态的图。
图3为图1所示的活体测量装置具备的容器周边的构成图。
图4为示意光源的构成例的图。
图5为示意光学干涉材料以及被测量部位的吸收系数与光的波长之间的关系一例的图。
图6为示意(a)(b)容器的内部中的吸收系数的空间分布的概念图,并示意内部信息测量结果。
图7为示意(a)(b)容器的内部中的吸收系数的空间分布的概念图,并示意轮廓测量结果。
图8为示意(a)(b)容器的内部中的吸收系数的空间分布的概念图,并示意对内部信息测量结果以及轮廓测量结果进行综合后的结果。
图9为示意光透过吸收系数均匀的介质的形态的图。
图10为示意光透过吸收系数不均匀的介质的形态的图。
图11为示意作为变形例的、光学干涉材料以及被测量部位的吸收系数与光的波长之间的关系的另一个例子的曲线图。
符号的说明
10…活体测量装置、12…容器、14…运算/控制部、16…光出射/检测端、18…遮光板、20…光学干涉材料、22…光源、24…光开关、26…光源用光纤、28…出射用光纤、30…光检测器、32…快门(shutter)、34…检测用光纤、36…信号处理电路、38…显示部、40a、40b…光源、42…光开关、B…被测量部位。
具体实施方式
以下,参照附图,对活体测量装置的实施方式进行详细说明。在此,在图的说明中,对同一要素标记同一符号,省略重复的说明。
图1为示意活体测量装置的一个实施方式的系统构成图。图2为示意图1所示的活体测量装置的使用状态的图。图3为图1所示的活体测量装置具备的容器周边的构成图。本实施方式的活体测量装置10是通过向作为测量对象的受检者的被测量部位B照射光并检测漫射光(返回光)从而取得被测量部位B的内部信息(肿瘤的有无等)的装置。在此,在本实施方式中,如图2所示,作为被测量部位B假定女性的乳房。
活体测量装置10具备:保持作为光透过性介质的光学干涉材料20并导入被测量部位B的容器12、向容器12的内部照射波长彼此不同的第1和第2光的光照射部、检测通过来自光照射部的光的照射而从被测量部位B产生的漫射光的光检测部、以及基于来自光检测部的输出信号计算被测量部位B的吸收系数的空间分布并计算出被测量部位B的内部信息的运算/控制部14(运算部)。
容器12具有能够充分地容纳被测量部位B的大小,且为在上表面具有开口部的圆筒状或者半球状。在容器12的内表面上,在相互不同的位置上三维配置有n(n为2以上的整数)个光出射/检测端16,从而构成测量部(构架(gantry))。从n个光出射/检测端16各自所具有的光出射端朝向被测量部位B依次射出光。该光在被测量部位B中被吸收和漫射,从被测量部位B射出的漫射光入射到n个光出射/检测端16各自的光检测端。在此,在本实施方式中,配置了n组由一组光出射端以及光入射端构成的光出射/检测端16,然而也可以在相互不同的位置独立地设置光出射端以及光检测端。
容器12由遮光性的材料构成,防止光从光出射/检测端16以外的地方入射到容器12的内部。此外,容器12的开口部可以由可拆卸的遮光板18覆盖。遮光板18被安装于容器12的开口部,则防止了光从开口部入射到容器12的内部。优选在容器12的内表面上进行用于降低漫射光的反射的处理。例如,容器12的内表面可以由用黑色染料着色的被阳极氧化处理(黑色耐酸铝处理)的铝材形成。或者,容器12的内表面也可以由黑色的树脂材料形成。
容器12的内部在活体测量时被光学干涉材料20充满。光学干涉材料20是通过充填被测量部位B与容器12之间的缝隙而起到降低被测量部位B的表面上的光学特性的不连续性的作用的液状介质。光学干涉材料20被构成为:在某个波长λ1下,使光学干涉材料20的吸收系数μai(λ1)与被测量部位B的吸收系数的平均值μab(λ1)实际上相等。此外,进一步优选光学干涉材料20被构成为:使其散射系数、折射率、旋光度、偏光度等光学特性中一个以上的特性与被测量部位B的散射系数的平均值、折射率的平均值、旋光度的平均值、偏光度的平均值等实际上相等。作为被测量部位B为活体的情况下的光学干涉材料20的一个例子,使用如下的溶液,即在折射率与活体大致相等的水中,针对散射系数而添加二氧化硅、脂肪乳剂(Intralipid)等;针对吸收系数而添加在规定波长下具有特有的吸收系数的墨水等;针于旋光度、偏光度而添加葡萄糖或果糖等,从而使该溶液的光学特性接近于被测量部位B。在此,所谓“实际上相等”,是指由测量精度等观点出发相同或者被认为相同。
本实施方式的光照射部由上述的n个光出射/检测端16各自所具有的光出射端、光源22、以及光开关24构成。作为光源22,例如能够使用激光二极管。本实施方式的光源22被构成为可以射出第1光(波长λ1)、以及波长与第1光不同的第2光(波长λ2)。作为这样的构成的例子,能够使用波长可变激光或使用波长选择器而切换使用波长,或者是如图4所示,能够使用具有分别产生波长λ1、λ2(λ1≠λ2)的光的光源40a、40b、以及通过选择性地切换光源40a、40b而选择波长的光开关42的光源。此外,作为第1和第2光的各自的波长λ1、λ2,由活体的透过率与应定量的吸收体部分的吸收系数之间的关系等,优选700nm~900nm左右的近红外线区域的波长。此外,为了缩短测量时间并减轻受检者的负担,优选第1和第2光以微小的时间差而被输出。
在此,图5为示意光学干涉材料20以及被测量部位B的吸收系数μai、μab与光的波长λ1、λ2之间的关系一例的图。在图5中,用实线表示的曲线G1表示了光学干涉材料20的吸收系数μai与光的波长之间的关系的一个例子(光学干涉材料20的吸收系数μai随着波长变大而变大的情况)。此外,用虚线表示的曲线G2~G4表示被测量部位B的吸收系数的平均值μab与光的波长之间的关系的一个例子(被测量部位B的吸收系数μab随着波长变大而变大的情况、几乎没有变化的情况、以及随着波长变大而变小的情况的三个例子)。
如图5所示,第1光的波长λ1被设定为使被测量部位B的吸收系数的平均值μab与光学干涉材料20的吸收系数μai大致相等(即μab(λ1)=μai(λ1))的波长。此外,第2光的波长λ2被设定为使被测量部位B的吸收系数的平均值μab大于光学干涉材料20的吸收系数μai(即μab(λ2)>μai(λ2))的波长。在图5所示的例子中,曲线G1的斜率大于曲线G2~G4的斜率,因而波长λ2被设定为比波长λ1短。
再次参照图1~图3。上述的第1和第2光例如作为连续光而从光源22被射出。从光源22被射出这些光从光出射/检测端16向被测量部位B照射。光开关24为1输入n输出的光开关,从光源22经由光源用光纤26而输入光,并将该光分别依次供给到上述n个光出射/检测端16。即光开关24逐一地依次选择连接于各光出射/检测端16的n条出射用光纤28,并光学地连接该出射用光纤28与光源22。
本实施方式的光检测部由上述的n个光出射/检测端16各自所具有的光检测端、分别对应于n个光出射/检测端16的n个光检测器30、以及被配置于各光检测器的输入部之前的n个快门(shutter)32构成。在n个光检测器30的各自上,经由检测用光纤34而输入入射于各光出射/检测端16的来自被测量部位B的漫射光。光检测器30根据到达所对应的光出射/检测端16的漫射光的光强度而生成模拟信号。作为光检测器30,除了光电倍增管(PMT:Photomultiplier Tube)之外,还能够使用光电二极管、雪崩光电二极管(avalanche photodiode)、PIN光电二极管等各种元件。优选光检测器30具有能够充分地检测波长λ1、λ2的波长成分的光谱灵敏度特性。此外,在来自被测量部位B的漫射光较为微弱时,优选使用高灵敏度或者高增益的光检测器。信号处理电路36连接于光检测器30的信号输出端,信号处理电路36对从光检测器30输出的模拟信号进行A/D变换并基于漫射光的光强度生成数字信号,并向运算/控制部14提供该数字信号。
运算/控制部14是基于从信号处理电路36提供的数字信号,进行与被测量部位B的内部信息、以及被测量部位B和光学干涉材料20之间的边界信息相关的解析运算的运算部。运算/控制部14利用例如具有CPU(Central Processing Unit)等的运算部以及存储器等的存储部的电脑而加以实现。此外,运算/控制部14还可以具有控制光源22的发光、光开关24的动作以及快门32的开闭的功能。此外,在运算/控制部14上连接有记录/显示部38,并可以使运算/控制部14中的运算结果、即被测量部位B的内部信息、以及被测量部位B和光学干涉材料20之间的边界信息可视化。
被测量部位B的内部信息的计算、即内部信息测量,以例如以下的方式进行。如图1所示,在容器12的内部充满了光学干涉材料20的状态下,将遮光板18安装于容器12。然后,从n个光出射/检测端16分别向容器12的内部依次照射第1光(波长λ1),并利用n个光检测器30,经由n个光出射/检测端16,检测通过了光学干涉材料20而漫射的光。此外,与此分开进行如下操作,即如图2所示,在容器12的内部充满了光学干涉材料20的状态下将被测量部位B浸于光学干涉材料20。然后,从n个光出射/检测端16分别向容器12的内部依次照射第1光(波长λ1),并利用n个光检测器30,经由n个光出射/检测端16,检测通过了光学干涉材料20而漫射的光。通过在运算/控制部14中比较这样得到的两个检测结果,从而对容器12的内部中的吸收系数的空间分布进行运算,得到与肿瘤等吸收体的位置以及形状相关的信息(内部信息)。
图6为示意由上述运算计算出的、容器12的内部中的吸收系数的空间分布的概念图。在图6(a)中,圆形的框表示容器12。此外,在容器12的内侧,吸收系数越大的区域,其着色越浓,存在于中央附近的吸收系数比较大的区域A1表示存在于被测量部位B的内部的肿瘤等吸收体。此外,在图6(b)中,纵轴表示吸收系数,横轴表示图6(a)的假想线C上的位置。在上述运算中,得到了区域A1(吸收体)相对于容器12整体的位置以及形状等。
此外,被测量部位B和光学干涉材料20之间的边界信息的计算、即轮廓测量,例如可以通过以下方式进行。如图2所示,在容器12的内部充满了光学干涉材料20的状态下,将被测量部位B浸于光学干涉材料20。然后,从n个光出射/检测端16分别向容器12的内部依次照射第2光(波长λ2),并利用n个光检测器30,经由n个光出射/检测端16,检测通过了光学干涉材料20而漫射的光。通过在运算/控制部14中比较这样得到的检测结果以及在前面所说明的在将遮光板18安装于容器12的状态下得到的检测结果,从而对容器12的内部中的吸收系数的空间分布进行运算,得到与被测量部位B的位置以及形状相关的信息(轮廓信息)。
图7为示意由上述运算计算出的、容器12的内部中吸收系数的空间分布的概念图。在图7(a)中,存在于中央附近的吸收系数比较大的区域A2表示被测量部位B所占据的区域。此外,在图7(b)中,纵轴表示吸收系数,横轴表示图7(a)的假想线C上的位置。在上述运算中,得到了区域A2(被测量部位)相对于容器12整体的位置以及形状等。
此后,运算/控制部14通过对图6以及图7所示的运算结果进行综合,从而如图8(a)、(b)所示,提供表示肿瘤等光吸收体相对于被测量部位B整体的存在位置的图像信息等。
接着,针对由运算/控制部14得到的吸收系数分布的运算方法,说明其基本原理。图9为示意光透过吸收系数均匀的介质的形态的图。图10为示意光透过吸收系数不均匀的介质的形态的图。为了便于说明,将作为散射吸收体的介质设想为二维放大的正方形,并将该介质分割为N(=25)个相等大小的正方形的体积单元(由于为二维,因而实际为面积单元)。设想在各个体积单元内吸收系数一定,且设想斜线部等所表示的体积单元与其它的体积单元相比,其吸收系数不同。
在从图9所示的吸收系数均匀(吸收系数为μa)的介质的一点,向介质内部投射光,并在一点检测输出光的情况下,使用入射光量I、各体积单元的影响度Wj(j=1~N)、表示由于散射以及反射等而使入射光射出到介质的外部的比率的衰减常数Dsr,将检测光量S表示为式
(1):
S=Dsr·I·exp{-μa(W1+W2……+WN)}…(1)
在此,所谓各体积单元的影响度,是指在从某一点投射光并在某一点检测光的情况下,伴随各体积单元的吸收系数的变化的检测光量的变化比率,在后面会叙述其具体的计算方法。
其次,如图10所示,使用基准的吸收系数μa、以及各体积单元的吸收系数相对于μa的变化量Δμaj(j=1~N),将在每个体积单元中具有不同的吸收系数μaj(j=1~N)的介质的各体积单元的吸收系数表示为式(2):
μaj=μa+Δμaj(j=1,2,…N)…(2)
如果衰减常数Dsr在吸收系数均匀的情况下没有变化,那么该情况下的检测光量O表示为式(3):
O=Dsr·I·exp{-[W1a+Δμa1)+W2a+Δμa2)+…+WNa+ΔμaN)]}
=S·exp{-[W1Δμa1+W2Δμa2+…+WNΔμaN]}…(3)
因此,通过在式(3)的两边取对数,导出以下的式(4):
ln S - ln O = ( W 1 Δμ a 1 + W 2 Δμ a 2 + . . . + W N Δμ aN )
= Σ j = 1 N W j Δμ aj . . . ( 4 )
在此,式(4)为从吸收系数均匀的介质的一点投射并在一点输出的光的检测光量S(以下称为基准光量S)、从吸收系数不均匀的介质的一点投射并在一点输出的光的检测光量O(以下称为测定光量O)、各体积单元内的影响度Wj(j=1~N)、以及相对于各体积单元的吸收系数相对于μa的变化量Δμaj(j=1~N)的函数。上述变量中,通过测量而得到基准光量S以及测定光量O,通过计算而得到各体积单元的影响度Wj(j=1~N)(在后面会详细叙述),因此未知数仅为N个各体积单元的吸收系数相对于μa的变化量Δμaj(j=1~N)。因此,通过针对不同的光投射点·光检测点组联立N个式(4)所示的方程式,可以求得N个Δμaj,并能够计算出介质的吸收系数的空间分布。
具体而言,将第i组(i=1~N)的光投射点·光检测点组中的基准光量设为Si、检测光量设为Oi、各体积单元的影响度设为Wij(j=1~N),则式(4)被表示为式(5):
ln S i - ln O i = Σ j = 1 N W ij Δμ aj . . . ( 5 )
在此,列出关于全部的i的式(5),以行列式表示,则为式(6):
Figure GPA00001111891900104
因此,可以如式(7)所示求得N个Δμaj、即介质的吸收系数的空间分布。
Figure GPA00001111891900111
在此,对各体积单元的影响度Wij(j=1~N)的求取方法进行说明。入射到各体积单元的连续光(光束)的稳定光漫射方程式被表示为式(8):
ΔΦ-μaD-1Φ=0 …(8)
在此,Φ:光束(每单位体积的光密度)
μa:各体积单元的光吸收系数
μs:各体积单元的光的各向同性散射系数
D:各体积单元的漫射系数
Figure GPA00001111891900112
此外,设介质的内部与外部的边界条件为式(9):
ΦBL=0          …(9)
其中,设下标BL表示介质的内部与外部的边界。此外,式(9)与光在该边界被完全吸收的条件、例如介质的周围被涂为全黑的状态等价。
使用式(8)以及式(9),针对各光投射点·光检测点组、即第i组(i=1~N)的光投射点·光检测点组进行的光的透过模拟(以下称为第1模拟),从而计算出检测光强度。其中,在第1模拟中,假定介质具有一定的吸收系数μa,此外假定在上述式(8)中为完全漫射,并使容器12的大小大于1/μs。将由第1模拟得到的第i组(i=1~N)光投射点·光检测点组中的检测光量设为di0
接着,使用式(8)以及式(9),进行第2模拟。在第2模拟中,假定介质的一个体积单元具有与吸收系数μa不同的吸收系数μa+Δμa,针对各光投射点·光检测点组进行光的透过模拟。例如通过使Δμa=0.01[mm-1],从而计算出检测光强度。对于第i组(i=1~N)光投射点·光检测点组,设在使第j组(j=1~N)体积单元的吸收系数发生了变化的情况下的检测光量为dij
使用通过上述第1模拟以及第2模拟而计算出的检测光量,从而将各体积单元的影响度Wij表示为式(10):
Wij=μa -1ln(dio/dij)…(10)
因此,由式(10)求得Wij,其结果,由式(7)计算出吸收系数的空间分布。
此外,针对式(4)所示的方程式的个数与体积单元的个数相等的情况进行了说明,但是,在方程式的个数少于体积单元的个数的情况、或者方程式的个数多于体积单元的个数的情况下,通过使用奇异值分解法等,能够将奇异问题变换为非奇异问题,因而也可以求得吸收系数的空间分布。
对本实施方式的活体测量装置10所达到的效果进行说明。如上所述,在活体测量装置10中,能够进行内部信息测量以及轮廓测量。在内部信息测量中,利用使得被测量部位B的吸收系数μab与光学干涉材料20的吸收系数μai实际上相等的波长λ1的光(第1光),测量容器12的内部的吸收系数的空间分布。因此,在该内部信息测量中,防止了被测量部位B的表面上的光的反射及散射等,从而能够高精度地测量以光的照射位置以及检测位置为基准的肿瘤等光吸收体的位置以及大小。另一方面,在轮廓测量中,利用使得被测量部位B的吸收系数μab大于光学干涉材料20的吸收系数μai的波长λ2的光(第2光),测量容器12的内部的吸收系数的空间分布。因此,在该轮廓测量中,能够高精度地测量以光的照射位置以及检测位置为基准的被测量部位B的轮廓。并且,通过对内部信息测量的结果以及轮廓测量的结果进行综合,可以测量光吸收体相对于被测量部位B整体的存在位置。
这样,根据本实施方式的活体测量装置10,在通过向受检者的被测量部位B照射光并检测漫射光从而取得被测量部位B的内部信息的装置中,不仅能够高精度地测量肿瘤等光吸收体的有无以及大小,而且能够高精度地测量光吸收体相对于被测量部位B整体的存在位置。
此外,如图5所示,光学干涉材料20对波长λ1的光(第1光)的吸收系数μai(λ1)可以大于光学干涉材料20对波长λ2的光(第2光)的吸收系数μai(λ2)。在波长λ1下的被测量部位B的吸收系数μab(λ1)以及波长λ2下的被测量部位B的吸收系数μab(λ2)彼此相等的情况(图5所示的曲线G3的情况)下,或者在被测量部位B的吸收系数μab(λ1)小于μab(λ2)的情况(图5所示的曲线G4的情况)下,通过这样设定光学干涉材料20的吸收系数μai,可以适宜地得到如下关系:在波长λ1下被测量部位B和光学干涉材料20的吸收系数μab(λ1)、μai(λ1)实际上相等且在波长λ2下被测量部位B的吸收系数μab(λ2)大于光学干涉材料20的吸收系数μai(λ2)。或者是,即使在被测量部位B的吸收系数μab(λ1)大于μab(λ2)的情况(图5所示的曲线G2的情况)下,通过以光学干涉材料20的吸收系数差(μai(λ1)-μai(λ2))大于被测量部位B的吸收系数差(μab(λ1)-μab(λ2))的方式进行光学干涉材料20的光吸收特性的调整以及波长λ1、λ2的设定,也能够适宜地实现上述关系。因此,能够适宜地进行内部信息测量以及轮廓测量。
(变形例)
图11为用于说明上述实施方式的变形例的曲线图,示意了光学干涉材料20以及被测量部位B的吸收系数μai、μab与光的波长λ1、λ2之间的关系的另一个例子的图。在图11中,用实线表示的曲线G5示意了光学干涉材料20的吸收系数μai与光的波长之间的关系的一个例子(光学干涉材料20的吸收系数μai无论波长的长短几乎没有变化的情况)。此外,用虚线表示的曲线G6示意了被测量部位B的吸收系数的平均值μab与光的波长之间的关系的一个例子(被测量部位B的吸收系数μab随着波长变大而变大的情况。
在光学干涉材料20以及被测量部位B具有如图11所示的光吸收特性的情况下,第1光的波长λ1被设定为使被测量部位B的吸收系数的平均值μab与光学干涉材料20的吸收系数μai大致相等的波长(即μab(λ1)=μai(λ1))。此外,第2光的波长λ2被设定为使被测量部位B的吸收系数的平均值μab大于光学干涉材料20的吸收系数μai的波长(即μab(λ2)>μai(λ2))。即波长λ2被设定为比波长λ1长。
这样,在被测量部位B的吸收系数μab(λ1)小于μab(λ2)的情况(即μab(λ1)<μab(λ2))下,也可以设定使光学干涉材料20的吸收系数μai(λ1)和μai(λ2)实际上相等。由此,能够适宜地得到如下关系:在波长λ1下被测量部位B和光学干涉材料20的吸收系数μab(λ1)、μai(λ1)实际上相等且在波长λ2下被测量部位B的吸收系数μab(λ2)大于光学干涉材料20的吸收系数μai(λ2)。因此,能够适宜地进行内部信息测量以及轮廓测量。
本发明的活体测量装置并不局限于上述的实施方式以及变形例,可以有其它的各种变形。例如,在上述实施方式中,作为光源例示了激光二极管,然而作为光源也能够使用固体激光、色素激光、或者气体激光。或者是,也可以使用通过波长选择器对来自LED或白色光源的光进行了波长选择而得到的光。
此外,在上述实施方式中,论述了使用连续光作为照射光时的解析方法,然而本发明也可以使用采用了时间分辨光谱法(TRS法:TimeResolved Spectroscopy)的活体测量装置、或者采用了调相谱检测技术(PMS法:Phase Modulation Spectroscopy)的活体测量装置,其中,时间分辨光谱法使用脉冲光,并由漫射光的时间分辨波形取得被检测部位的内部信息,调相谱检测技术使用位相调制光,从而由漫射光的光强度以及相位信息取得被检测部位的内部信息。
此外,在上述实施方式中,作为被测量部位例示了乳房,但是,本发明的活体测量装置也可以适用于测量头、手、脚、躯干等其它的活体部位。
此外,作为充满容器内的介质(光学干涉材料),不局限于液状的介质,也可以使用例如随着时间的流逝而固化的液体、凝胶状的材质的介质。通过使用由随着时间的流逝而固化的材质构成的介质,可以固定被测量部位,并提高测量精度,同时由于能够在舒适的姿势下进行测量,因而能够减轻受检者的负担。

Claims (3)

1.一种活体测量装置,其特征在于,
是通过向受检者的被测量部位照射光并检测漫射光从而取得所述被测量部位的内部信息的活体测量装置,
所述活体测量装置具备:
保持光透过性介质的容器;
具有被固定于所述容器的多个光出射端,并向浸于所述介质的所述被测量部位照射波长彼此不同的第1和第2光的光照射部;
具有被固定于所述容器的多个光检测端,并检测来自所述被测量部位的漫射光的光检测部;以及
基于来自所述光检测部的输出信号计算出所述内部信息的运算部,
所述第1光的波长为所述被测量部位的吸收系数与所述介质的吸收系数的平均值实际上相等的波长,
所述第2光的波长为所述被测量部位的吸收系数大于所述介质的吸收系数的平均值的波长,
所述运算部基于与所述第1光的所述漫射光相关的所述输出信号而计算出所述内部信息,并基于与所述第2光的所述漫射光相关的所述输出信号而计算出所述被测量部位与所述介质之间的边界信息。
2.如权利要求1所述的活体测量装置,其特征在于,
所述介质对所述第1光的吸收系数大于所述介质对所述第2光的吸收系数。
3.如权利要求1所述的活体测量装置,其特征在于,
所述介质对所述第1光的吸收系数与所述介质对所述第2光的吸收系数实际上相等。
CN2008801129347A 2007-10-24 2008-10-10 活体测量装置 Active CN101836104B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2007-276770 2007-10-24
JP2007276770A JP5420163B2 (ja) 2007-10-24 2007-10-24 生体計測装置
PCT/JP2008/068485 WO2009054278A1 (ja) 2007-10-24 2008-10-10 生体計測装置

Publications (2)

Publication Number Publication Date
CN101836104A true CN101836104A (zh) 2010-09-15
CN101836104B CN101836104B (zh) 2012-01-25

Family

ID=40579383

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008801129347A Active CN101836104B (zh) 2007-10-24 2008-10-10 活体测量装置

Country Status (5)

Country Link
US (1) US8279443B2 (zh)
EP (1) EP2204645B1 (zh)
JP (1) JP5420163B2 (zh)
CN (1) CN101836104B (zh)
WO (1) WO2009054278A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107613850A (zh) * 2015-04-09 2018-01-19 通用医疗公司 用于时间分辨漫射相关谱的系统和方法
CN108209867A (zh) * 2016-12-15 2018-06-29 松下知识产权经营株式会社 摄像装置

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11444248B2 (en) * 2016-11-30 2022-09-13 Kyushu University, National University Corporation Organic electro-luminescent element and bioinstrumentation device

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0760091B1 (de) * 1994-05-19 1999-11-03 Roche Diagnostics GmbH Verfahren und vorrichtung zur bestimmung eines analyten in einer biologischen probe
DE4417639A1 (de) * 1994-05-19 1995-11-23 Boehringer Mannheim Gmbh Verfahren zur Bestimmung eines Analyten in einer biologischen Probe
EP0808124B1 (en) * 1995-01-03 2003-04-16 Non-Invasive Technology, Inc. Optical coupler for in vivo examination of biological tissue
US5813988A (en) * 1995-02-03 1998-09-29 Research Foundation Time-resolved diffusion tomographic imaging in highly scattering turbid media
JP3844815B2 (ja) * 1996-08-30 2006-11-15 浜松ホトニクス株式会社 散乱体の吸収情報計測方法及び装置
US6687532B2 (en) * 1997-12-12 2004-02-03 Hamamatsu Photonics K.K. Optical CT apparatus and image reconstructing method
JP3771364B2 (ja) * 1997-12-12 2006-04-26 浜松ホトニクス株式会社 光ct装置及び画像再構成方法
JP3950243B2 (ja) * 1998-11-05 2007-07-25 浜松ホトニクス株式会社 散乱吸収体の内部情報の計測方法及び装置
DE10050232A1 (de) * 2000-10-11 2002-05-02 Karlsruhe Forschzent Hochauflösender Ultraschalltomograph
FR2900043B1 (fr) * 2006-04-24 2008-07-04 Commissariat Energie Atomique Procede d'imagerie optique par fluorescence de tissus biologiques, notamment pour delimiter des regions d'interet des tissus a analyser par tomographie
WO2008011112A2 (en) * 2006-07-19 2008-01-24 University Of Connecticut Method and apparatus for medical imaging using combined near-infrared optical tomography, fluorescent tomography and ultrasound

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107613850A (zh) * 2015-04-09 2018-01-19 通用医疗公司 用于时间分辨漫射相关谱的系统和方法
CN108209867A (zh) * 2016-12-15 2018-06-29 松下知识产权经营株式会社 摄像装置

Also Published As

Publication number Publication date
WO2009054278A1 (ja) 2009-04-30
EP2204645A4 (en) 2012-12-12
CN101836104B (zh) 2012-01-25
US8279443B2 (en) 2012-10-02
JP2009103619A (ja) 2009-05-14
EP2204645B1 (en) 2013-12-04
US20110063624A1 (en) 2011-03-17
EP2204645A1 (en) 2010-07-07
JP5420163B2 (ja) 2014-02-19

Similar Documents

Publication Publication Date Title
JP3844815B2 (ja) 散乱体の吸収情報計測方法及び装置
JP3662376B2 (ja) 内部特性分布の計測方法および装置
US6549284B1 (en) Calibration methods and systems for diffuse optical tomography and spectroscopy
US7428434B2 (en) Quantitative broadband absorption and scattering spectroscopy in turbid media by combined frequency-domain and steady state methodologies
US9277866B2 (en) Method and apparatus for analysis of turbid media via single-element detection using structured illumination
US7440659B2 (en) Depth-resolved reflectance instrument and method for its use
Gulsen et al. Design and implementation of a multifrequency<? xpp qa?> near-infrared diffuse optical tomography system
Kacprzak et al. Time-resolved optical imager for assessment of cerebral oxygenation
US20110105865A1 (en) Diffuse reflectance spectroscopy device for quantifying tissue absorption and scattering
JP3660761B2 (ja) 散乱体の吸収情報の計測方法及び装置
US7818154B2 (en) Monte Carlo based model of fluorescence in turbid media and methods and systems for using same to determine intrinsic fluorescence of turbid media
JPH0749307A (ja) 散乱吸収体計測方法及び装置
JPH11337476A (ja) 散乱吸収体の内部特性分布の計測方法及び装置
Patwardhan et al. Quantitative diffuse optical tomography for small animals using an ultrafast gated image intensifier
US6704110B2 (en) Method and apparatus for measuring internal information of scattering medium
CN101836104B (zh) 活体测量装置
Song et al. Statistical analysis of nonlinearly reconstructed near-infrared tomographic images. II. Experimental interpretation
JPH11311569A (ja) 内部特性分布の計測方法及び装置
JP4077477B2 (ja) 散乱体の吸収情報計測方法及び装置
JP4077476B2 (ja) 散乱体の吸収情報計測方法及び装置
JP3878943B2 (ja) 散乱体の吸収情報の計測方法及び装置
JP4077475B2 (ja) 散乱体の吸収情報計測方法及び装置
JP2012005556A (ja) 生体機能情報の光計測装置
Qin et al. Determination of the optical properties of turbid materials by hyperspectral diffuse reflectance
RU2637102C1 (ru) Устройство для спектрофотометрической оценки уровня кровенаполнения поверхностных слоев тканей и органов человека in vivo

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant