CN105814435A - 光学活性物质的浓度计算系统、光学活性物质的浓度计算系统的制作方法、计算机可读介质及程序 - Google Patents
光学活性物质的浓度计算系统、光学活性物质的浓度计算系统的制作方法、计算机可读介质及程序 Download PDFInfo
- Publication number
- CN105814435A CN105814435A CN201580003067.3A CN201580003067A CN105814435A CN 105814435 A CN105814435 A CN 105814435A CN 201580003067 A CN201580003067 A CN 201580003067A CN 105814435 A CN105814435 A CN 105814435A
- Authority
- CN
- China
- Prior art keywords
- optically active
- active substance
- concentration
- value
- optical rotation
- 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
Links
Classifications
-
- 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/21—Polarisation-affecting properties
-
- 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/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/483—Physical analysis of biological material
- G01N33/487—Physical analysis of biological material of liquid biological material
- G01N33/49—Blood
- G01N33/492—Determining multiple analytes
-
- 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/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N2021/3129—Determining multicomponents by multiwavelength light
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Hematology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Ecology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Urology & Nephrology (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
数据处理单元(30)获得与多个波长λ1,λ2,λ3,…对应的旋光度αM1,αM2,αM3,…,并用非线性最小二乘法根据获得的成对的旋光度和波长(λ1:αM1,λ2:αM2:λ3:αM3,…)以及公式计算第一光学活性物质的浓度C,该公式用第一非线性函数和第二非线性函数的总和表示,第一非线性函数表示第一光学活性物质的旋光色散,第二非线性函数包括用于限定第二光学活性物质的旋光色散特性的固有值作为未知值。
Description
技术领域
本发明涉及一种光学活性物质的浓度计算系统、光学活性物质的浓度计算系统的制作方法、计算机可读介质及程序。
背景技术
PTL1(专利文献1)公开了一种尿分析方法。在该尿分析方法中,测量了尿的旋光角度(旋光角度的范围是已知的),并且确定光学活性物质的浓度C[kg/dl]在(A-Ah)/(α×L)≤C≤(A-Al)/(α×L)的范围内。尿的旋光角度用除了浓度未知的光学活性物质之外的受干扰光学活性物质(interferenceoptically-activesubstance)表示。“A”表示测量出的尿的旋光角度[deg]。“Ah”表示用受干扰光学活性物质表示的旋光角度的最大值[deg]。“Al”表示用受干扰光学活性物质表示的旋光角度的最小值[deg]。“α”表示光学活性物质的旋光率[deg/cm·dl/kg]。“L”表示测量出的光程长度[cm]。
引用列表
专利文献
[PTL1]JP-A-9-138231
发明内容
技术问题
本发明的目的在于提供一种如下的光学活性物质的浓度计算系统等:在使用光学活性物质的旋光度来计算光学活性物质的浓度的情况下,即使当测量目标物质中所包含的所有种类光学活性物质均未知时,该系统也可以计算出想要得到的光学活性物质的浓度。
解决问题的技术方案
[1]根据本发明的一个方面,提供了一种光学活性物质的浓度计算系统,其是根据表示旋光度的公式计算光学活性物质的浓度的系统。所述公式包括:第一函数,其表示第一光学活性物质的旋光度的波长相关性;以及第二函数,其表示第二光学活性物质的旋光度的波长相关性。在所述第一函数中,所述第一光学活性物质的浓度具有未知值,并且用于限定所述第一光学活性物质的旋光色散特性的固有值是已知值或是一定受限范围内的未知值。在所述第二函数中,用于限定所述第二光学活性物质的旋光色散特性的固有值是未知值。使用最小二乘法根据所述公式以及测量目标的分别与多个波长对应的旋光度来计算所述第一光学活性物质的所述浓度。
[2]在根据第[1]项所述的光学活性物质的浓度计算系统中,所述第二光学活性物质包含多种光学活性物质,所述第二函数是如下的函数:以单个项表示所述多种光学活性物质的旋光度的波长相关性的集合,并且所述第二函数是如下的函数:所述单个项中用于限定所述多种光学活性物质的旋光色散特性的固有值用未知值来表示。
[3]在根据第[1]或[2]项所述的光学活性物质的浓度计算系统中,在所述第一函数中用于限定所述第一光学活性物质的旋光色散特性的至少一些固有值是一定受限范围内的未知值。
[4]在根据第[1]至[3]项中任一项所述的光学活性物质的浓度计算系统中,在所述第二函数中用于所述第二光学活性物质的至少一些固有值是一定受限范围内的未知值。
[5]在根据第[1]至[4]项中任一项所述的光学活性物质的浓度计算系统中,在改变所述第一函数中所包含的未知值和所述第二函数中所包含的未知值中的至少一个未知值的初始值的同时,根据最小二乘法多次计算所述第一光学活性物质的所述浓度。当多次计算所述第一光学活性物质的所述浓度时,计算出的所述第一光学活性物质的与使目标函数的值最大的初始值对应的浓度不作为所述第一光学活性物质的所述浓度。所述目标函数用所述最小二乘法中的旋光度的理论值与获得的旋光度之差的平方和表示。
[6]在根据第[1]至[4]项中任一项所述的光学活性物质的浓度计算系统中,在对为所述第一函数中所包含的未知值限制一定范围的限制条件以及为所述第二函数中所包含的未知值限制一定范围的限制条件中的至少一个限制条件进行改变的同时,根据最小二乘法多次计算所述第一光学活性物质的所述浓度。当多次计算所述第一光学活性物质的所述浓度时,计算出的所述第一光学活性物质的与使目标函数的值最大的限制条件对应的浓度不作为所述第一光学活性物质的所述浓度。所述目标函数用所述最小二乘法中的旋光度的理论值与获得的旋光度之差的平方和表示。
[7]在根据第[1]至[6]项中任一项所述的光学活性物质的浓度计算系统中,对传输通过人的眼球中的房水的光而言,获得所述测量目标的分别与所述多个波长对应的旋光度。
[8]根据本发明的另一个方面,提供一种光学活性物质的浓度计算系统的制作方法。所述制作方法包括:获得分别与多个光线的波长对应的旋光度的组合;对旋光度的每个组合而言,使用最小二乘法根据公式和获得的旋光度的组合来计算第一光学活性物质的浓度,所述公式包括:第一函数,其表示所述第一光学活性物质的旋光度的波长相关性;以及第二函数,其表示第二光学活性物质的旋光度的波长相关性;以及将所述第一光学活性物质的彼此相关的旋光度的组合和浓度存储在存储单元中,所述第一光学活性物质的所述浓度与旋光度的组合对应。在所述第一函数中,所述第一光学活性物质的所述浓度具有未知值,并且用于限定所述第一光学活性物质的旋光色散特性的固有值是已知值或是一定受限范围内的未知值。在所述第二函数中,用于限定所述第二光学活性物质的旋光色散特性的固有值是未知值。
[9]根据本发明的另一个方面,提供了一种计算机可读介质,其是非临时性的并存储使计算机执行根据表示旋光度的公式计算光学活性物质的浓度的功能的程序。所述公式包括:第一函数,其表示第一光学活性物质的旋光度的波长相关性;以及第二函数,其表示第二光学活性物质的旋光度的波长相关性。在所述第一函数中,所述第一光学活性物质的浓度具有未知值,并且用于限定所述第一光学活性物质的旋光色散特性的固有值是已知值或是一定受限范围内的未知值。在所述第二函数中,用于限定所述第二光学活性物质的旋光色散特性的固有值是未知值。使用最小二乘法根据所述公式以及测量目标的分别与多个波长对应的旋光度来计算所述第一光学活性物质的所述浓度。
[10]根据本发明的另一个方面,提供了一种使计算机执行根据表示旋光度的公式计算光学活性物质的浓度的功能的程序。所述公式包括:第一函数,其表示第一光学活性物质的旋光度的波长相关性;以及第二函数,其表示第二光学活性物质的旋光度的波长相关性。在所述第一函数中,所述第一光学活性物质的浓度具有未知值,并且用于限定所述第一光学活性物质的旋光色散特性的固有值是已知值或是一定受限范围内的未知值。在所述第二函数中,用于限定所述第二光学活性物质的旋光色散特性的固有值是未知值。使用最小二乘法根据所述公式以及测量目标的分别与多个波长对应的旋光度来计算所述第一光学活性物质的所述浓度。
本发明的有益效果
根据第[1]项的构造,即使当测量目标物质中所包含的所有种类光学活性物质均未知时,也可以计算出想要得到的光学活性物质的浓度。
根据第[2]项的构造,与多个光学活性物质的旋光度的波长相关性以多个项表示的情况相比,可以减少公式中所包含的未知值的数量。
根据第[3]和[4]项的构造中的每一个构造,与未知值的范围不受限的情况相比,更加提高了光学活性物质的浓度的计算精度。
根据第[5]和[6]项的构造中的每一个构造,与没有执行多次计算的情况相比,更加提高了光学活性物质的浓度的计算精度。
根据第[7]项的构造,非侵入性地得到眼球的房水中所包含的特定光学活性物质的浓度。
根据第[8]项的构造,与没有存储波长和旋光度的组合的对应关系以及与该组合对应的光学活性物质的浓度的情况相比,可以更快地计算光学活性物质的浓度。
根据第[9]项的构造,即使当测量目标物质中所包含的所有种类光学活性物质均未知时,也可以计算出想要得到的光学活性物质的浓度。
根据第[10]项的构造,即使当测量目标物质中所包含的所有种类光学活性物质均未知时,也可以计算出想要得到的光学活性物质的浓度。
附图说明
图1是示出了实施例所应用的光学活性物质的浓度计算系统的构造实例的视图。
图2是示出了旋光率[α]与波长λ的相关性的视图。
图3是示出了在测量目标物质包含作为光学活性物质的葡萄糖和白蛋白的情况下观测旋光度αM与葡萄糖的旋光度αg和白蛋白的旋光度αa的总和对应的视图。
图4是示出了以使用非线性最小二乘法来应用函数的方式得到基于关于多个波长λ的观测旋光度αM的函数g(λ)的方法的视图。
图5是示出了从观测旋光度αM计算出葡萄糖浓度Cg的计算实例1的视图;这里,图5中的(a)示出了波长λ与观测旋光度αM之间的关系;图5中的(b)示出了葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系。
图6是示出了使用公式(10)从观测旋光度αM计算出葡萄糖浓度Cg的计算实例2的视图;这里,图6中的(a)示出了波长λ与观测旋光度αM之间的关系;图6中的(b)示出了葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系。
图7是示出了使用公式(10)并使用与图6中的410nm至470nm的短波长区域不同的波长区域中的波长λ来从观测旋光度αM得到葡萄糖浓度Cg的情况的视图;这里,图7中的(a)示出了660nm至850nm的长波长区域中的波长λ与观测旋光度αM之间的关系;图7中的(b)示出了图7中的(a)的葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系;图7中的(c)示出了短波长区域和长波长区域中的波长λ与观测旋光度αM之间的关系;图7中的(d)示出了图7中的(c)的葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系。
图8是示出了使用公式(11)从观测旋光度αM得到葡萄糖浓度Cg的计算实例3的视图;这里,图8中的(a)示出了波长λ与观测旋光度αM之间的关系;图8中的(b)示出了葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系。
图9是用于将在使用公式(11)来计算葡萄糖浓度(计算值)Cgc时的不对未知值设置限制条件的情况与对未知值设置限制条件的情况进行比较的视图;这里,图9中的(a)示出了不对未知值设置限制条件的情况;图9中的(b)和图9中的(c)示出了对未知值设置限制条件的情况。
图10是用于将在使用公式(11)来计算葡萄糖浓度(计算值)Cgc时所设定的葡萄糖浓度(计算值)Cgc的初始值进行比较的视图;这里,图10中的(a)示出了初始值是100ml/dl的情况;图10中的(b)示出了初始值是300ml/dl的情况;图10中的(c)示出了初始值是500ml/dl的情况。
图11是示出了从在限制条件A和B下计算出的光学活性物质的浓度C得到优选的计算值的实例的视图;这里,图11中的(a)示出了在限制条件A的情况下葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系;图11中的(b)示出了在限制条件B的情况下的葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系;图11中的(c)示出了针对具有某个葡萄糖浓度Cgr的测量目标物质计算出的目标函数的值与葡萄糖浓度Cgc之间的关系。
图12是示出了数据处理单元的概况的视图。
图13是示出了实例1中的数据处理单元的功能单元的视图。
图14是示出了实例2中的数据处理单元的功能单元的视图。
图15是示出了存储在公式存储单元中的公式以及存储在固有值存储单元中的固有值的实例的视图;这里,图15中的(a)示出了由公式存储单元存储的公式(算法)的实例;图15中的(b)示出了由固有值存储单元存储的固有值的实例。
图16是示出了实例2中的数据处理单元的硬件构造的实例的视图。
图17是示出了实例3中的数据处理单元的功能单元的视图。
图18是示出了存储在LUT存储单元中的LUT的实例的视图。
图19是示出了在对葡萄糖浓度测量装置应用光学活性物质的浓度计算系统的情况下在UI单元中所包含的显示器显示光学活性物质的浓度的显示实例的视图;这里,图19中的(a)示出了以数值进行显示的显示实例;图19中的(b)示出了以条形图(水平)进行显示的显示实例;图19中的(c)示出了以标记“OK/NG”进行显示的显示实例;图19中的(d)示出了使用数值除了显示葡萄糖浓度之外还显示其他光学活性物质的浓度的显示实例。
具体实施方式
接下来,将参考附图对根据本发明的实施例进行描述。
(光学活性物质的浓度计算系统1)
图1是示出了实施例所应用的光学活性物质的浓度计算系统1的实例的视图。
光学活性物质具有使得进行照射的线偏振光的偏振面旋转的旋光性。这里,在本实施例中,偏振面指的是对于线偏振光而言,电场振动的平面。
图1所示的光学活性物质的浓度计算系统1利用线偏振光照射包含光学活性物质的测量目标物质13,并通过测量在线偏振光传输通过测量目标物质13之后所观察到的偏振面的旋转角度(观测旋光度αM)来计算测量目标物质13中所包含的光学活性物质的浓度。
这里,在测量目标物质13包含一种光学活性物质的情况下所观测到的旋光度用旋光度α表示,而在测量目标物质13包含多种光学活性物质的情况下所观测到的旋光度表示为旋光度αM。如上所述,观测旋光度αM被设定成测量目标物质13中所包含的光学活性物质的旋光度α的总和。测量目标物质13所包含的光学活性物质的每个旋光度α反应了测量目标物质13所包含的相应光学活性物质的浓度。
即使当测量目标物质13包含多种光学活性物质时,也可以知道想要得到的光学活性物质的浓度。
想要得到的光学活性物质的浓度指的是在测量目标物质13包含多种光学活性物质情况下使用者想要知道的光学活性物质的浓度。想要得到的光学活性物质的浓度还指的是作为UI单元40(见图2)等的显示目标的光学活性物质的浓度。这种光学活性物质可以表示成第一光学活性物质,而其余光学活性物质可以表示成第二光学活性物质。可以提供多种第一光学活性物质。第二光学活性物质是除了第一光学活性物质之外的光学活性物质。然而,在提供除了第一光学活性物质之外的多种光学活性物质的情况下,可以不需要对多种光学活性物质进行区分,并且可以将除了第一光学活性物质之外的多种光学活性物质的一部分或全部的集合设定成第二光学活性物质。
光学活性物质的浓度计算系统1包括测量单元10、控制单元20、数据处理单元30和UI单元40。测量单元10连接至控制单元20。控制单元20连接至数据处理单元30。UI单元40与控制单元20和数据处理单元30连接。连接可以是有线或无线的。
测量单元10利用线偏振光照射包含光学活性物质的测量目标物质13,并测量在线偏振光传输通过测量目标物质13之后所观察到的偏振面的旋转角度(观测旋光度αM)。
控制单元20控制测量单元10测量观测旋光度αM并将测量数据发送至数据处理单元30。通过将测量中所使用的波长λ和观测旋光度αM组合起来而得到测量数据。
数据处理单元30由从控制单元20接收的测量数据计算测量目标物质13中所包含的光学活性物质的浓度,上述测量数据是波长λ和观测旋光度αM的组合。数据处理单元30借助于数值计算处理执行计算并将计算出的浓度发送至UI单元40。
UI单元40包括输入装置和输出装置。输入装置接收数据输入或来自使用者的指令,并包括键盘等。输出装置向使用者显示处理结果等,并包括显示器等。
使用者借助于UI单元40的例如键盘等输入装置向控制单元20给出对测量单元10的操作指示。使用者向数据处理单元30输入公式或固有值(将在下文中进行描述)。
使用者借助于UI单元40的例如显示器等输出装置从数据处理单元30得到测量单元10的操作状态或想要得到的光学活性物质的浓度。
测量单元10包括光源11、偏光器12、补偿器14、分析器15和光接收元件16。光源11发射具有预定波长的光。偏光器12从光源11发射的光中提取出预定偏振面的线偏振光。补偿器14使传输通过测量目标物质13的线偏振光的偏振面旋转。分析器15使预定偏振面的线偏振光传输通过分析器15。光接收元件16接收传输通过分析器15的光。这些部件构成一个光学系统。
在图1所示的测量单元10中,光源11与偏光器12之间、偏光器12与测量目标物质13之间、测量目标物质13与补偿器14之间、补偿器14与分析器15之间以及分析器15与光接收元件16之间分别设置有位于圆圈中的箭头。每个圆圈中的箭头均表示从光的行进方向观察到的偏振形式。
这里,假设测量目标物质13包含多种光学活性物质。
光源11可以是具有宽带宽的例如发光二极管(LED)或灯等光源。光源11可以是具有窄带宽的例如激光等光源。能够以具有至少两个以上的波长的光进行照射的光源被用作光源11。在例如测量目标物质13中所包含的光学活性物质的旋光度α可以近似于鲁德(Drude)单项式(将在下文中进行描述)的区域中的波长范围被用作光的波长范围。波长范围的实例是400nm至900nm。
这里,假设光源11发射的光包括具有如图1所示的随机偏振面的光。光源11发射的光可以是线偏振光。在该情况下,可以不需要偏光器12(接下来将要描述)。
偏光器12是例如尼科耳(Nicol)棱镜等。具有随机偏振面的光入射到偏光器12,并且偏光器12使具有预定偏振面的线偏振光透过。在图1中,作为实例,假设透过了与纸面平行的偏振面的线偏振光。
就传输通过偏光器12的线偏振光而言,测量目标物质13中所包含的光学活性物质使线偏振光的偏振面旋转。在图1中,假设偏振面足足旋转了角度αM(观测旋光度αM)。
补偿器14是例如磁光元件,诸如使用石榴石等的法拉第元件。补偿器14利用磁场使线偏振光的偏振面旋转。
与偏光器12类似,分析器15使预定偏振面的线偏振光透过。
光接收元件16是例如硅二极管等的光接收元件。光接收元件16输出与光的强度对应的输出信号。
上述测量单元10仅是实例,并可以包括例如反射镜、透镜、波长板和棱镜等其他光学元件。
接下来,将描述由测量单元10执行的观测旋光度αM的测量方法的实例。
首先,在没有放置(移除)测量目标物质13的状态下,光学系统中的补偿器14和分析器15设定为使光接收元件16的输出信号最小。光学系统由光源11、偏光器12、补偿器14、分析器15和光接收元件16构成。在没有放置测量目标物质13的状态下,透过偏光器12的线偏振光的偏振面与透过光检测器15的光的偏振面垂直。
在图1中,偏光器12的偏振面以及光透过分析器15之前的偏振面与纸面更严格地平行。然而,在补偿器14使偏振面旋转的情况下,在光透过分析器15之前的偏振面可以相对于与纸面平行的表面倾斜。也就是说,在没有将测量目标物质13置于测量单元10中的情况下,补偿器14和分析器15可以设定为使光接收元件16的输出信号最小。
然后,将测量目标物质13置于测量单元10中。测量目标物质13中所包含的光学活性物质使偏振面旋转。因此,从光接收元件16输出的输出信号偏离最小值。相应地,设定对补偿器14施加的磁场使得从光接收元件16输出的输出信号最小。也就是说,补偿器14使偏振面旋转,并且使偏振面与透过分析器15的光的偏振面垂直。
因补偿器14而旋转的偏振面的角度与测量目标物质13中所包含的光学活性物质所引起的观测旋光度αM对应。已预先知道了对补偿器14施加的磁场的大小与旋转的偏振面的角度之间的关系。因此,从对补偿器14施加的磁场的大小中发现观测旋光度αM。
描述了使用补偿器14的实例作为得到观测旋光度αM的方法。然而,可以通过使用除了补偿器14之外的部件来得到观测旋光度αM。图1示出了正交偏光器方法(使用补偿器),该方法为用于测量旋转角度(旋光度)的最基本的测量方法。然而,可以应用例如旋转分析器方法或法拉第调制方法和光学延迟调制方法等其他测量方法。
利用来自光源11的具有多个波长λ1,λ2,λ3,…的光照射测量目标物质13,并且针对多个波长λ1,λ2,λ3,…中的每一个波长,得到每个观测旋光度αM1,αM2,αM3,…。以该方式,在提供了多个波长λ和观测旋光度αM的情况下,执行如下的表示:波长λ1,λ2,λ3,…和观测旋光度αM1,αM2,αM3,…。类似地,对其他情况应用这种表示。
控制单元20控制光源11发射的光的波长λ1,λ2,λ3,…的设定(切换),并控制开启/关断(On/Off)。控制单元20设定对补偿器14施加的磁场以使光接收元件16的输出信号最小。控制单元20向数据处理单元30发送测量数据(λ1:αM1,λ2:αM2,λ3:αM3,…)。测量数据是波长λ与通过利用对补偿器14施加的磁场且利用波长λ而测量出的观测旋光度αM的组合。
数据处理单元30根据测量数据(λ1:αM1,λ2:αM2,λ3:αM3,…)(为波长λ与观测旋光度αM的组合)来计算测量目标物质13中所包含的光学活性物质的浓度。
活体中所包含的大部分物质是具有旋光性的光学活性物质。因此,可以将利用了光学活性物质的旋光性的浓度测量方法应用于对活体中所包含的光学活性物质的浓度的测量。
例如,血液中的葡萄糖浓度称为血糖水平并被广泛地用作糖尿病等的指标。具有与血清的成分大致相同的成分的房水还包含如下的许多种光学活性物质:这些光学活性物质具有旋光性并包含葡萄糖、例如白蛋白和球蛋白等蛋白质以及抗坏血酸。血液中的葡萄糖浓度与房水中的葡萄糖浓度之间存在的相关性是已知的。
房水具有高透明性。因此,如果允许利用旋光性测量房水的葡萄糖浓度,则可以使用非侵入性的葡萄糖浓度测量方法。
然而,旋光度α根据光学活性物质的种类而变化是已知的。如果混合有多种光学活性物质的测量目标物质13被线偏振光照射,则不会观察到每种光学活性物质的旋光度α,而是观察到受到测量目标物质中所包含的所有光学活性物质的旋光度α影响的观测旋光度αM。
例如,在利用旋光性得到葡萄糖浓度的情况下,需要使与设定作为目标的葡萄糖对应的信号(浓度、旋光度α等)与观测旋光度αM分开。
在以下描述中,就测量目标物质13中所包含的多种光学活性物质而言,将描述从测量出的观测旋光度αM计算设定为要得到的光学活性物质的浓度的方法。
(旋光度α)
首先,将描述旋光度α。
关于包含某个单种光学活性物质的测量目标物质13,对波长λ而言的旋光度α用公式(1)表示。也就是说,旋光度α(deg)用旋光率[α](deg/(dm·g/ml))、光程长度L(dm)和浓度C(g/ml)的乘积表示。旋光率[α]不包括光学活性物质的浓度C或光程长度L。旋光率[α]是恒定温度下的光学活性物质固有的常量。
α=[α]·L·C公式(1)
图2是示出了旋光率[α]与波长λ的相关性的视图。旋光率[α]的波长相关性称为旋光色散(opticalrotatorydispersion)。
在长度比处于最大点与最小点之间的长度长的波长区域中旋光率[α]可以用单调递减或单调递增的鲁德单项式表示。鲁德单项式是表示光学活性物质的旋光色散的函数的实例。鲁德单项式是用公式(2)表示的非线性函数。在图2中,在超过最大点的波长区域中,鲁德单项式示出为单调递减的函数。这里,“λ”是变量,而“A”和“λ0”是光学活性物质固有的常量(用于限定光学活性物质中的旋光色散特性的固有值)。
在测量目标物质13包含多种光学活性物质的情况下,观测旋光度αM用公式(3)表示。也就是说,对包含多种光学活性物质的测量目标物质13进行观察而得到的观测旋光度αM描述成:将光学活性物质的旋光度αi(旋光率[α]用公式(2)的鲁德单项式表示)加起来。也就是说,观测旋光度αM用表示各光学活性物质的旋光度的波长相关性的函数的总和表示。作为实例,函数用鲁德单项式和浓度的乘积表示,或用鲁德单项式、浓度和光程长度等的乘积表示。就光学活性物质的固有值Ai和λi而言,提供了已知的情况和未知的情况。
这里,在测量目标物质13包含多种光学活性物质的情况下,就每种光学活性物质而言,执行如下表示:旋光度αi、浓度Ci、旋光率[αi]、固有值Ai和λi。可以使用与光学活性物质的名称相关的符号替代下标“i”。例如,在葡萄糖的情况下,可以使用下标“g”等。因为光程长度L由测量目标物质13确定,所以光学活性物质的光程长度L彼此相同。
在测量目标物质13仅包含葡萄糖和白蛋白作为光学活性物质的情况下,观测旋光度αM用公式(4)表示。也就是说,观测旋光度αM用指示葡萄糖的旋光度的波长相关性的函数以及指示白蛋白的旋光度的波长相关性的函数的总和来表示。这里,Ag和λg是葡萄糖的固有值,而Aa和λa是白蛋白的固有值。葡萄糖与浓度Cg(葡萄糖浓度Cg)对应,而白蛋白与浓度Ca(白蛋白浓度Ca)对应。
图3是示出了在测量目标物质13包含作为光学活性物质的葡萄糖和白蛋白的情况下观测旋光度αM是葡萄糖的旋光度αg与白蛋白αa的旋光度的总和的视图。
因为如图3所示葡萄糖表现出顺时针的旋光性(右旋)而白蛋白表现出逆时针的旋光性(左旋),所以葡萄糖的旋光度αg是正的并且相对于波长λ为单调递减。白蛋白的旋光度αa是负的并且相对于波长λ为单调递增。
如图3所示,观测旋光度αM由葡萄糖的旋光度αg与白蛋白的旋光度αa的总和表示。
在图3的实例中,葡萄糖对观测旋光度αM的影响比白蛋白对观测旋光度αM的影响大。如公式(1)所示,旋光度α由光程长度L、旋光率[α]和浓度C的乘积确定。因此,在图3中,与白蛋白相比,葡萄糖具有较大的旋光率[α]和/或较高的浓度C。
如图3所示,一般来说,在房水中,葡萄糖对观测旋光度αM的影响比白蛋白对观测旋光度αM的影响大,并且与其他光学活性物质相比,葡萄糖和白蛋白具有较大的旋光率[α]和较高的浓度C。
(光学活性物质的浓度Ci的计算方法)
接下来,就包含多种光学活性物质的测量目标物质13而言,将描述从观测旋光度αM计算想要得到的光学活性物质的浓度Ci的方法。在实施例中,如公式(3)所示,至少想要得到的光学活性物质的浓度Ci设定为具有使用表示测量目标物质13的旋光度的波长相关性的公式(理论公式)而得到的未知值,并且利用非线性最小二乘法计算浓度Ci。
例如,在假设房水是测量目标物质13的情况下,可以认为的是,至少15种以上的右旋和左旋的光学活性物质彼此混合在房水中。房水中的光学活性物质的实例包括抗坏血酸、乳酸、葡萄糖、丙氨酸、精氨酸、半胱氨酸、谷氨酸、组氨酸、亮氨酸、异亮氨酸、赖氨酸、丝氨酸、缬氨酸、白蛋白和球蛋白。
因此,可以认为的是,根据由所有约15种光学活性物质的旋光度αi的总和表示的公式(如公式(3)所示),应用计算例如葡萄糖等特定光学活性物质的浓度Ci的方法。
然而,因为房水位于人的眼球中并且在眼泪、尿等的情况下不能容易地收集排出至身体外部的物质,所以对所包含的光学活性物质种类的精确检验是困难的且是不明确的。因为可以认为房水中所包含的光学活性物质的种类或浓度比率也会因个人的生活习惯或药的剂量而变化,所以可以认为:即使当根据由所有约15种光学活性物质的旋光度αi的总和表示的公式来单单计算例如葡萄糖等特定光学活性物质的浓度Ci时,也难以以期望的精度进行计算。
因此,在实施例中,如下所述进行数值计算。得到特定光学活性物质的旋光度αi以及其余光学活性物质的集合的旋光度αx的总和。使包含多种光学活性物质的测量目标物质13的观测旋光度αM近似于得到的总和,并且将其余光学活性物质的固有值设定成未知值。
也就是说,如公式(5)所示,观测旋光度αM近似于特定光学活性物质的旋光度αi和其余光学活性物质的集合的旋光度αx的总和。也就是说,观测旋光度αM设定成以下两项的总和:包括鲁德单项式且表示特定光学活性物质的旋光色散的函数之和;以及X(λ,A),其表示其余光学活性物质的旋光度的波长相关性的集合的旋光色散。
这里,“n”表示测量目标物质13中所包含的光学活性物质的全部种类。“k”表示在测量目标物质13中所包含的所有光学活性物质之中的特定光学活性物质数。因此,满足1≤k<n。例如,在测量目标物质13包含15种光学活性物质的情况下,满足n=15。就特定光学活性物质而言,以高精度计算出对测量目标物质13的观测旋光度αM的影响程度大的物质(旋光率的绝对值大的物质)的浓度。因此,如果设定成包含影响程度最大的物质,则提高了计算出的物质浓度的精度。在多种物质被选定作为特定光学活性物质的情况下,可以按影响程度大的物质的次序选择多种物质。
特定光学活性物质中的至少一种可以是想要得到浓度的光学活性物质。想要得到浓度的光学活性物质的浓度设定为具有未知值。如下文将要描述的,想要得到浓度的光学活性物质的固有值设定为具有已知值或一定受限范围内的未知值。就X(λ,A)而言,其余光学活性物质的至少一个固有值被设定成未知值。
如下文将要描述的,在已知与多种其余光学活性物质的旋光度αx成正比的X(λ,A)对观测旋光度αM具有小影响的情况下,X(λ,A)可以被设定成“0”。
例如,在测量目标物质13设定为包含15种光学活性物质并且葡萄糖和白蛋白设定作为15种光学活性物质中的想要得到浓度的光学活性物质的情况下,将公式(6)用作公式(5)。公式(6)中的第一项与葡萄糖的旋光度αg对应。公式(6)中的第二项与白蛋白的旋光度αa对应。公式(6)中的第三项与除了葡萄糖和白蛋白之外的其余光学活性物质的集合的旋光度αx对应。第三项中的Ax和λx与用于对其余光学活性物质的固有值的集合进行限定的常量对应。Cx与其余光学活性物质的浓度对应。
例如,在得到了葡萄糖的浓度Cg的情况下,第一项与第一函数(非线性函数)对应,而第二项和第三项对应于第二函数(非线性函数)。在得到了葡萄糖的浓度Cg和白蛋白的浓度Ca的情况下,第一项和第二项对应于第一函数,而第三项与第二函数对应。也就是说,与第一光学活性物质的浓度相关的函数是第一函数的实例。与第二光学活性物质的浓度相关的函数是第二函数的实例。
图4是示出了以使用非线性最小二乘法来应用函数的方式基于关于多个波长λ的观测旋光度αM得到函数g(λ)的方法。函数g(λ)是表示测量目标物质13的旋光度的公式(理论公式),并与公式(5)对应。然而,这里,对函数g(λ)执行一般设定。
如图2和3所示,观测旋光度αM对波长λ的波长相关性是非线性的。因此,使用非线性最小二乘法来执行函数g(λ)的函数应用。也就是说,在对想要得到的浓度Ci等的未知值应用数值时,得到了使平方和最小的函数g(λ)。
如图4所示,假设提供了测量值(分别与波长λ(λ1,λ2,λ3和λ4)对应的观测旋光度αMj(αM1,αM2,αM3和αM4))。如果提供了相应的测量值,则函数g(λj)与测量值(观测旋光度αMj)之差rj用公式(7)表示。差rj的平方和用公式(8)表示。平方和的一般形式用公式(9)表示。公式(9)是表示真值(理论值)与计算值(将在下文中进行描述)之差的目标函数。
r1=g(λ1)-αM1
r2=g(λ2)-αM2
r3=g(λ3)-αM3
r4=g(λ4)-αM4公式(7)
r1 2+r2 2+r3 2+r4 2公式(8)
在非线性最小二乘法中,设定未知值并应用函数,使得函数g(λj)与作为测量值的观测旋光度αMj之差的平方和最小。例如将列文伯格-马夸尔特(Levenberg-Marquardt)法、拟牛顿(quasi-Newton)法和共轭梯度法用作该方法。
这里,使用列文伯格-马夸尔特法来计算想要得到的光学活性物质的浓度C。为了改进高斯-牛顿法的解的收敛不稳定性而开发出列文伯格-马夸尔特法,并且列文伯格-马夸尔特法被广泛地用作非线性最小二乘法的算法。因为列文伯格-马夸尔特法是已知的方法,所以将省略对列文伯格-马夸尔特法的描述。
多个观测旋光度αM(即,对两个以上的波长λ而言的观测旋光度αM)用作测量数据。随着波长λ的数量增加,提高了得到的未知值的精度。然而,如果波长λ的数量大,则需花费一定时间执行数值计算。因此,可以根据情况来选择波长λ的数量。
这里,使用了列文伯格-马夸尔特法,但可以应用其他算法。
替代这种最小二乘法,可以通过解联立方程来计算想要得到的光学活性物质的浓度。在联立方程中,需要使测量中所使用的波长λ的数量与未知值的数量匹配。在采用非线性最小二乘法的列文伯格-马夸尔特法的函数应用中,即使当测量中所使用的波长λ的数量不与未知值的数量匹配时,也可以计算出想要得到的光学活性物质的浓度C。也就是说,测量中所使用的波长λ的数量可以小于或大于未知值的数量,只要使用最小二乘法即可。
在以下描述中,将描述通过仿真来计算葡萄糖浓度Cg的计算实例。
这里,计算混合有包括葡萄糖在内的19种光学活性物质的测量目标物质13中的葡萄糖的浓度(计算值)Cgc。以使用列文伯格-马夸尔特法的非线性最小二乘法来从对多个波长λ而言的观测旋光度αM进行计算。将计算出的浓度Cgc与测量目标物质13中所包含(混合)的葡萄糖的浓度(真值)Cgr进行比较。
光学活性物质的种类和浓度Ci被设定在可识别的范围内,以使测量目标物质13近似于人的房水。因为身体中存在变化,所以设定了具有可能性的上限、下限及中间这三个条件。在图5至图8(将在下文中进行描述)中,三个条件用符号○、□和Δ表示。
(计算实例1)
图5是示出了由观测旋光度αM计算葡萄糖浓度Cg的计算实例1的视图。图5中的(a)是示出了波长λ与观测旋光度αM之间的关系的视图。图5中的(b)是示出了葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系的视图。
在计算实例1中,使用公式(6)来计算混合有包括葡萄糖在内的19种光学活性物质的测量目标物质13中的葡萄糖浓度Cg。
这里,在公式(6)中,将葡萄糖浓度Cg、白蛋白浓度Ca、白蛋白的固有值Aa以及固有值λx设定成未知值,并且将X(λ,A)中的Ax和Cx的乘积设定成一个未知值,其中固有值λx为通过集合除了葡萄糖和白蛋白之外的其余光学活性物质而得到的X(λ,A)中的常量。
葡萄糖的固有值Ag和λg以及白蛋白的固有值λa被设定成已知值,并分别被设定成1.72×107、150和250。
在测量观测旋光度αM中使用图5中的(a)所示的四个波长λ。对波长λ1(410nm)而言得到观测旋光度αM1(2.89mdeg)。对波长λ2(420nm)而言得到观测旋光度αM2(2.79mdeg)。对波长λ3(440nm)而言得到观测旋光度αM3(2.60mdeg)。对波长λ4(470nm)而言得到观测旋光度αM4(2.33mdeg)。观测旋光度αM的分辨能力可以测量到两个小数位。
由上述结果可以了解的是,在从波长λ1(410nm)至波长λ4(470nm)的波长范围中,观测旋光度αM单调递减。
通过将葡萄糖浓度Cg的初始值设定成500mg/dl来开始计算。
如图5中的(b)所示,在以使用列文伯格-马夸尔特法的非线性最小二乘法计算出的葡萄糖浓度(计算值)Cgc与在葡萄糖浓度(真值)Cgr是0mg/dl至400mg/dl的范围内的葡萄糖浓度(真值)Cgr之间得到正比关系。三个条件(○、□和Δ)中的差(变化量)比葡萄糖浓度Cg的值小。也就是说,可以以高精度从观测旋光度αM计算混合有包括葡萄糖在内的19种光学活性物质的测量目标物质13的葡萄糖浓度Cg。
如上所述,在例如房水等测量目标物质13中所包含的其余光学活性物质的种类或浓度比率不明确的情况下,很显然可以认为浓度计算是困难的。然而,在计算实例中,X(λ,A)中的固有值被设定成未知值,因此,可以以期望的精度计算出浓度。
在上述计算中,虽然白蛋白固有值Aa被设定成未知值,但是白蛋白固有值Aa也可以被设定成已知值。例如,在白蛋白浓度Ca本身就是想要得到的浓度的情况下,使用已知值作为白蛋白固有值Aa。在白蛋白浓度Ca本身不是想要得到的浓度的情况下,可以将白蛋白浓度Ca与白蛋白固有值Aa的乘积设定成一个未知值。
也就是说,在上述计算中,未知值的数量被设定成5。然而,未知值的数量可以小于或大于5。
波长λ的数量被设定成4。然而,可以设置多个波长λ。波长λ的数量可以设定成2或3。波长λ的数量可以大于4。
此外,在数值计算中,需要设定未知值的初始值。在实例1中,葡萄糖浓度Cg被设定成500mg/dl(作为初始值)。然而,葡萄糖浓度Cg可以大于或小于500mg/dl。
(计算实例2)
在计算实例2中,在公式(6)中将除了葡萄糖和白蛋白之外的其余光学活性物质的旋光度αx(X(λ,A))设定成“0”。也就是说,使用公式(10)替代公式(6)。公式(10)是公式(4)。
图6是示出了使用公式(10)从观测旋光度αM计算葡萄糖浓度Cg的计算实例2的视图。图6中的(a)是示出了波长λ与观测旋光度αM之间的关系的视图。图6中的(b)是示出了葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系的视图。
葡萄糖浓度Cg、白蛋白浓度Ca和白蛋白的固有值Aa被设定成具有未知值,并且在葡萄糖的固有值Ag和λg以及白蛋白的固有值λa中使用已知值。因为Aa和Ca是简单的乘积,所以可以将组合设定成一个未知值。
在测量观测旋光度αM中使用图6中的(a)所示的四个波长。这四个波长与图5中的(a)所示的四个波长一样。然而,对波长λ而言的观测旋光度αM的分辨能力被设定为具有一个小数位。这是因为假设观测旋光度αM为亚毫米级别,并且因为澄清了如下情况的影响:即,难以得到观测旋光度αM的分辨能力的情况。
如图6中的(b)所示,在以使用列文伯格-马夸尔特法的非线性最小二乘法计算出的葡萄糖浓度(计算值)Cgc与在葡萄糖浓度(真值)Cgr是0mg/dl至600mg/dl的范围内的葡萄糖浓度(真值)Cgr之间得到正比关系。三个条件(○、□和Δ)中的差(变化量)比葡萄糖浓度Cg的值小。也就是说,在公式(6)中,即使当将除了葡萄糖和白蛋白之外的其余光学活性物质的集合的旋光度αx(X(λ,A))设定成“0”并且将观测旋光度αM的分辨能力设定成具有一个小数位时,也可以以高精度从混合有包括葡萄糖在内的19种光学活性物质的测量目标物质13的观测旋光度αM计算葡萄糖浓度Cg。
虽然白蛋白的固有值Aa被设定成未知值,但白蛋白的固有值Aa也可以被设定成已知值。在白蛋白的固有值Aa和λa被设定成已知值的情况下,并且在房水中所包含的光学活性物质的种类或浓度比率因个人的生活习惯或药的剂量而变化的情况下,应考虑不能以期望的精度计算浓度的可能性。在这种情况下,可以使用计算实例1中的公式(6)或计算实例3中的公式(11)(将在下文中进行描述)。
在对除了房水之外的其他测量目标应用公式(10)的情况下,可以按对测量目标物质13的观测旋光度的影响程度大的物质(旋光率的绝对值大的物质)的次序选择多种物质,并且可以将其他物质的旋光度设定成“0”。
图7是示出了使用公式(10)并使用与图6中的410nm至470nm的短波长区域不同的波长区域的波长λ来从观测旋光度αM得到葡萄糖浓度Cg的情况的视图。图7中的(a)示出了660nm至850nm的长波长区域中的波长λ与观测旋光度αM之间的关系。图7中的(b)示出了图7中的(a)的葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系。图7中的(c)示出了短波长区域和长波长区域中的波长λ与观测旋光度αM之间的关系。图7中的(d)示出了图7中的(c)的葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系。
如图7中的(a)所示,对波长λ5(600nm)而言得到观测旋光度αM5(1.2mdeg)。对波长λ6(690nm)而言得到观测旋光度αM6(1.1mdeg)。对波长λ7(780nm)而言得到观测旋光度αM7(0.9mdeg)。对波长λ8(850nm)而言得到观测旋光度αM8(0.7mdeg)。旋光度αM被设定成具有一个小数位。
如图7中的(a)所示,在从波长λ5(660nm)至波长λ8(850nm)的波长范围中观测旋光度αM单调递减。然而,观测旋光度αM的变化比在使用具有图6中的(a)所示的短波长区域的光的情况下观测旋光度αM的变化小。
如图7中的(b)所示,在葡萄糖浓度(计算值)Cgc与在葡萄糖浓度(真值)Cgr是0mg/dl至600mg/dl的范围内的葡萄糖浓度(真值)Cgr之间得到正比关系。然而,三个条件(○、□和Δ)中的差(变化量)比图6中的(b)所示的情况大。可以认为的是,这是因为对波长λ而言的观测旋光度αM的变化是小的。也就是说,在图2中,这是因为长波长区域对应于旋光率[α]的变化小的波长λ区域。
与之相反,在图7的(c)中,使用图6中的(a)所示的对波长λ3(440nm)而言的观测旋光度αM3(2.6mdeg)以及对波长λ4(470nm)而言的观测旋光度αM4(2.3mdeg),并且使用图7中的(a)所示的对波长λ7(780nm)而言的观测旋光度αM7(0.9mdeg)以及对波长λ8(850nm)而言的观测旋光度αM8(0.7mdeg)。在波长λ3(440nm)至波长λ8(850nm)的波长范围中观测旋光度αM单调递减。观测旋光度αM的变化比在使用图6中的(a)所示的短波长区域的情况下以及在使用图7中的(a)所示的长波长区域的情况下的观测旋光度αM的变化大。
如图7中的(d)所示,与图6中的(b)和图7中的(b)所示的情况类似,在葡萄糖浓度(计算值)Cgc与在葡萄糖浓度(真值)Cgr是0mg/dl至600mg/dl的范围内的葡萄糖浓度(真值)Cgr之间得到正比关系。三个条件(○、□和Δ)中的差(变化量)比图7中的(b)所示的情况以及图6中的(b)所示的情况小。可以认为的是,这是因为对波长λ而言的观测旋光度αM的变化是大的。也就是说,如果波长λ选择为被包括在从短波长区域至长波长区域的范围内使得观测旋光度αM大幅变化,则提高了想要得到的光学活性物质的浓度Ci的计算精度。
为了提高浓度Ci的计算精度,光源11可以构造成由不同的半导体材料形成的半导体激光器,从而选择例如多个波长λ以包括可见光波长区域中的从短波长区域至长波长区域的范围。有源层由氮化镓的混晶(GaN、AlGaInN、GaInN等)形成的半导体激光器可以用作具有蓝光区域(波长λ为400nm至500nm)的光源11。有源层由磷化镓的混晶(AlGaInP、GaInP等)形成的半导体激光器可以用作具有红光区域(波长λ为600nm至700nm)的光源11。有源层由砷化镓的混晶(AlGaAs、GaAs等)形成的半导体激光器可以用作具有红光区域或红外线区域(波长λ在800nm附近)的光源11。
这里所述的半导体激光器是示例性的,并且可以选择和使用与波长λ对应的半导体激光器。
半导体激光器是紧凑的并具有窄带宽,并因此有助于减小光学活性物质的浓度计算系统1的尺寸。
(计算实例3)
在计算实例3中,使用公式(11)。通过用除了葡萄糖之外的其余光学活性物质的集合所表述的旋光度αx替代公式(10)中的白蛋白的旋光度αa来得到公式(11)。除了葡萄糖之外的其余光学活性物质的集合被设定成浓度Cx,而除了葡萄糖之外的其余光学活性物质的集合被设定成常量Ax和λx,因此,使用单个项。
图8是示出了使用公式(11)从观测旋光度αM得到葡萄糖浓度Cg的计算实例3的视图。图8中的(a)示出了波长λ与观测旋光度αM之间的关系。图8中的(b)示出了葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系。
葡萄糖浓度Cg、除了葡萄糖之外的其余光学活性物质的集合的浓度Cx以及除了葡萄糖之外的其余光学活性物质的集合的常量Ax和λx被设定成未知值。1.72×107和150分别用作葡萄糖的固有值Ag和λg。因为Ax和Cx是简单的乘积,所以将Ax和Cx的组合作为一个未知值进行计算。以这种方式将固有值和浓度的乘积作为未知值的情况是固有值为未知值的实例。
在测量观测旋光度αM中使用图8中的(a)所示的四个波长。这四个波长与图7中的(c)所示的四个波长一样。采用了使得观测旋光度αM的差变大的波长区域。
如图8中的(b)所示,在葡萄糖浓度(计算值)Cgc与在葡萄糖浓度(真值)Cgr是0mg/dl至600mg/dl的范围内的葡萄糖浓度(真值)Cgr之间得到正比关系。三个条件(○、□和Δ)中的差(变化量)是小的。也就是说,在公式(10)中,即使当白蛋白的旋光度αa被设定成除了葡萄糖之外的其余光学活性物质的集合的旋光度αx时,也可以以高精度从混合有包括葡萄糖在内的19种光学活性物质的测量目标物质13的观测旋光度αM计算葡萄糖浓度Cg。在对除了房水之外的其他测量目标物质13应用公式(11)的情况下,可以选择对测量目标物质13的观测旋光度αM具有最大程度影响的物质(旋光率的绝对值为最大的物质),并且可以将其余光学活性物质的集合设定成旋光度αx。
在上述计算中,葡萄糖的固有值Ag和λg被设定成1.72×107和150。然而,可以指定固有值Ag和λg的范围。例如,在固有值因例如温度等环境波动而变化的情况下,可以使用在预先考虑波动的范围内的未知值(即,一定受限范围内的未知值)。例如,固有值λg可以在145至155的范围内。这与对未知值设置限制条件的情况(将在下文中进行描述)类似。
在想要得到白蛋白浓度Ca而不是葡萄糖浓度Cg的情况下,可以使用公式(12)替代公式(11)。
如上所述,可以对假设被包含在房水(其被设定作为测量目标物质13的实例)中的葡萄糖应用公式(6)、公式(10)和公式(11)。
在以上描述中,描述了使用最小二乘法来计算浓度的实例。然而,可以通过解联立方程来计算浓度。在该情况下,可以将波长λ的数量设定为与联立方程中的未知值的数量相同。
因为公式(11)简单,所以用于公式(11)的计算时间比用于公式(6)的计算时间短。因此,在想要减少计算时间的情况下,可以不使用公式(6)而使用公式(11)。就其他测量目标物质13而言,可以根据想要得到的光学活性物质的种类或测量目标物质13中所包含的光学活性物质的种类来改变和使用公式(6)、公式(10)和公式(11)中的光学活性物质。
在实施例中,如公式(6)或公式(11)所示,用于限定除了想要得到浓度的光学活性物质之外的光学活性物质的旋光色散特性的固有值被设定成未知值。因此,就房水而言,针对除了想要得到浓度的光学活性物质之外的光学活性物质的种类或浓度比率不明确的测量目标物质13,计算想要得到浓度的光学活性物质的浓度C。
如公式(5)所示,观测旋光度αM用数量比测量目标物质13中所包含的光学活性物质的数量小的非线性函数的总和表示。也就是说,就除了想要得到浓度的光学活性物质之外的多种光学活性物质而言,使用如下的函数:其中,旋光度的波长相关性的集合用单个项表示。因此,不需要利用分别与测量目标物质13中所包含的所有光学活性物质对应的非线性函数的总和来表示公式。因此,减少了公式中所包含的未知值的数量。因此,例如,可以减少计算时间。
即使当使用非线性最小二乘法来执行函数应用并且因此未知值的数量比测量中所使用的波长λ的数量大时,也可以计算出想要得到的光学活性物质的浓度C。因此,如果将测量中所使用的波长λ的数量设定成大于未知值的数量,则与数量小的情况相比提高了浓度计算的精度。
在以上描述中,旋光色散近似于公式(2)所示的鲁德单项式。然而,旋光色散可以近似于单调递减或单调递增的其他非线性函数。此外,如果不使用从长波长区域至短波长区域的大范围内的波长并且例如仅使用长波长区域中的波长或仅使用短波长区域中的波长,则旋光率的倾斜比使用大范围内的波长的情况更接近直线。因此,在这种情况下,可以通过设定为接近线性函数并通过使用线性最小二乘法来计算浓度。
(改进计算值的相关性的方法)
需要改进计算值与真值的相关性(正比关系)。因此,将对改进计算值与真值的相关性的方法进行描述。
利用例如列文伯格-马夸尔特法和拟牛顿法等梯度法(其用作不搜索函数的最小值而搜索某个区域中的最小值的数值计算算法)来进行局部优化计算。因此,如果不能进行适当区域中的计算,则得到的解(计算值)与真值具有大的差异。即使当使用给定的初始值(初始估计的值)时,得到的解(计算值)也可能与真值具有大的差异。
即使在这种情况下,也为光学活性物质的浓度C的未知值、固有值等设置用于限制某个范围的限制条件。因此,计算值和真值彼此匹配或彼此接近。
也通过替换未知值的初始值(初始估计的值)使计算值和真值彼此匹配或彼此接近。
也就是说,提高了计算值的精度。
首先,将描述为光学活性物质的浓度C的未知值、固有值等设置限制条件的情况。
图9是使用公式(11)来计算葡萄糖浓度(计算值)Cgc时将不为未知值设置限制条件的情况与为未知值设置限制条件的情况进行比较的视图。图9中的(a)示出了不为未知值设置限制条件的情况。图9中的(b)和图9中的(c)示出了为未知值设置限制条件的情况。图9中的(a)至图9中的(c)示出了葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系。
在图9中的(a)至图9中的(c)中,符号◇、□、Δ、■和O表示在混合有包括葡萄糖在内的19种光学活性物质的测量目标物质13中光学活性物质的混合比不同的情况。
葡萄糖浓度Cgc的初始值(将在下文中进行描述)被设定成300mg/dl。
在图9中的(a)中不设置限制条件的情况下,在一定程度上获得真值与计算值之间的相关性。然而,计算值(绘制点)稍微宽泛地分散。
在图9中的(b)中设置了设定Cgc>0且Ax<0限制条件的情况下,与图9中的(a)的情况类似,计算值(绘制点)稍微宽泛地分散。
然而,在图9中的(c)中设置了设定Cgc>0、Ax<0且250≤λx≤300限制条件的情况下,与图9中的(a)至图9中的(b)不同的是,计算值(绘制点)几乎不分散,并且葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的相关性(正比关系)是高的。
如上所示,如果为光学活性物质的浓度C的未知值、固有值等设置限制条件,则计算值与真值之间的相关性增大。如图9中的(b)至图9中的(c)所示,在提供了许多未知光学活性物质且混合比率(浓度比率)是未知的情况下,难以给出合适的限制条件。
接下来,将描述替换光学活性物质的浓度C的未知值的初始值、固有值等的情况。这里,将描述替换光学活性物质的浓度C的初始值的情况。
图10是将使用公式(11)来计算葡萄糖浓度(计算值)Cgc时所设定的葡萄糖浓度(计算值)Cgc的初始值进行比较的视图。图10中的(a)示出了初始值是100ml/dl的情况。图10中的(b)示出了初始值是300ml/dl的情况。图10中的(c)示出了初始值是500ml/dl的情况。
图10中的(a)至图10中的(c)中的符号◇、□、Δ、■和O与图9中的(a)至图9中的(c)中情况下的符号相同。
关于未知值的限制条件被设定为Cgc>0、Ax<0且250≤λx≤300(即,与图9中的(c)的限制条件一样)。
在图10中的(a)所示的葡萄糖浓度Cgc的初始值被设定成100mg/dl的情况下,如区域I所示,许多计算值(绘制点)大大偏离葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的正比关系(1:1)。也就是说,计算值大大背离真值。
与之相反,在图10中的(b)所示的葡萄糖浓度Cgc的初始值被设定成300mg/dl的情况下,计算值(绘制点)分布成接近葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的正比关系。
此外,在图10中的(c)所示的葡萄糖浓度Cgc的初始值被设定成500mg/dl的情况下,计算值(绘制点)分布成更接近葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的正比关系(1:1)。
如图10中的(a)所示,以设定光学活性物质的浓度C的初始值的方式使计算值大大背离真值。因此,在提供了许多未知光学活性物质且混合比率(浓度比率)未知的情况下,难以给出合适的初始值。
在实施例中,对于未知值的限制条件不限于某个值或某个范围,并且光学活性物质的浓度C设定为在多个限制条件(通过改变限制条件来得到)下来计算。因此,在分别在多个限制条件下得到的计算值之中,优选的计算值被设定作为光学活性物质的浓度C。
类似地,未知值的初始值也不限于某个值,并且光学活性物质的浓度C设定为使用多个初始值(通过改变初始值而得到)来计算。因此,在分别使用多个初始值而得到的计算值之中,优选的计算值被设定作为光学活性物质的浓度C。
以多个的方式提供用于未知值的限制条件和初始值中的一者或两者。也就是说,可以提供限制条件和初始值的多个组合,并且可以计算出光学活性物质的浓度C。
同时,需要设定使真值与计算值之差小的计算值作为优选的计算值。然而,因为真值不是已知的,所以需要用于选择优选的计算值的指标。因此,这里,与真值和计算值之差对应的公式(9)被设定作为目标函数,并且优选的计算值被设定为利用目标函数的值而选择得到。因为由公式(9)表示的目标函数的值变小,所以真值与计算值之差是小的。
因此,公式(9)可以设定为目标函数,可以将对于未知值而言的多个限制条件的目标函数的值、多个初始值等彼此比较,从而可以选出优选的计算值。
图11是示出了从在限制条件A和B下计算出的光学活性物质的浓度C得到优选的计算值的实例的视图。这里,图11中的(a)示出了在限制条件A的情况下的葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系。图11中的(b)示出了在限制条件B的情况下的葡萄糖浓度(真值)Cgr与葡萄糖浓度(计算值)Cgc之间的关系。图11中的(c)是示出了针对具有某个葡萄糖浓度Cgr的测量目标物质13计算得到的目标函数的值与葡萄糖浓度Cgc之间的关系的表格。
图11中的(a)中的限制条件A与Cgc>0、Ax<0且250≤λx≤300对应。葡萄糖浓度Cgc的初始值被设定成300mg/dl。在限制条件A中,在一定程度上获得真值与计算值之间的相关性,但计算值(绘制点)稍微宽泛地分散。
图11中的(b)中的限制条件B与Cgc>0、Ax<0且250≤λx≤300对应。葡萄糖浓度Cgc的初始值被设定成300mg/dl。限制条件B与图9中的(c)一样。在限制条件B中,计算值(绘制点)不像在限制条件A中那样分散,并且葡萄糖浓度(真值)Cgr和葡萄糖浓度(计算值)Cgc具有接近正比关系的关系。
也就是说,如果使用限制条件B计算葡萄糖浓度Cgc,则与使用限制条件A的情况相比,提高了浓度计算的精度。
如图11中的(c)所示,就具有某个葡萄糖浓度Cgr的测量目标物质13而言,在限制条件A中计算出的目标函数的值是3.5×10-3,并且葡萄糖浓度是230mg/dl。与之相比,在限制条件B中计算出的目标函数的值是1.56×10-5,并且葡萄糖浓度是198mg/dl。
就限制条件A和限制条件B而言,限制条件B中的目标函数的值比限制条件A中的目标函数的值小。也就是说,真值与计算值之差是小的。因此,限制条件B下所计算出的198mg/dl可以被选择作为葡萄糖浓度Cgc。
图11中的(a)和图11中的(b)示出了:可以确认的是,如果选择限制条件B,则真值与计算值之差是小的。图11中的(c)示出了:根据目标函数的值选择更接近真值的计算值。
在以上描述中,描述了提供两个限制条件作为限制条件A和B的情况。然而,可以提供三个以上的限制条件。在这种情况下,可以选择在目标函数具有最小值的其他限制条件中所计算出的光学活性物质的浓度C。可以排除在目标函数具有最大值的限制条件中所计算出的光学活性物质的浓度C,而可以选择在其他限制条件中所计算出的光学活性物质的其他浓度C。此外,通过执行基于其他限制条件的计算而得到的结果可以用作光学活性物质的浓度C,例如,可以对在限制条件中所计算出的光学活性物质的其他浓度C求平均,并且求平均的浓度可以用作光学活性物质的浓度C。可以设定限制条件之中目标函数具有最大值的多个排除限制条件,目标函数具有最大值的限制条件之后的限制条件等。
在以上描述中,描述了对未知值应用多个限制条件并执行计算的情况。然而,以上描述可以类似地应用于对使用多个初始值执行一个未知值的计算的情况或使用限制条件和初始值的多个组合来执行未知值的计算的情况。例如,对于光学活性物质的浓度C为使用多个初始值来计算未知值的实例的情况,可以使用图10中的(a)至图10中的(c)所示的多个初始值来计算光学活性物质的浓度C,并且可以选择优选的计算值。
使用多个限制条件或多个初始值来执行参量的计算,或者使用多个限制条件和多个初始值的组合来执行计算。因此,得到与真值具有高的相关性的计算值。
这里,根据公式(11)计算光学活性物质的浓度C(葡萄糖浓度Cgc)。然而,可以使用其他公式(6)、(10)和(12)来执行计算。
(数据处理单元30)
在以下描述中,将描述数据处理单元30的实例。
图12是示出了数据处理单元30的概况的视图。数据处理单元30输出设定为从测量数据(λ1:αM1,λ2:αM2:λ3:αM3,…)和光程长度L得到的光学活性物质的浓度C。通过组合波长λ和使用波长λ测得的观测旋光度αM来得到测量数据。
<实例1>
图13是示出了实例1中的数据处理单元30的功能单元的视图。
数据处理单元30与控制单元20和UI单元40连接。数据处理单元30包括浓度计算单元101。
数据处理单元30的浓度计算单元101根据已从控制单元20发送来的测量数据(λ1:αM1,λ2:αM2,λ3:αM3,…)和光程长度L来计算光学活性物质的浓度C。通过组合波长λ与观测旋光度αM来得到测量数据。浓度计算单元101将计算出的光学活性物质的浓度C发送至UI单元40。
浓度计算单元101是例如用作专用于计算葡萄糖浓度Cgc的机器,并由例如专用集成电路(ASIC)等硬件构成。也就是说,在浓度计算单元101中,用于根据公式(6)、公式(10)、公式(11)等来计算葡萄糖浓度Cgc的算法由硬件(硬件模块)构成。
如果从控制单元20输入测量数据(λ1:αM1,λ2:αM2,λ3:αM3,…)和光程长度L,则浓度计算单元101立即计算葡萄糖浓度Cgc并输出计算出的葡萄糖浓度Cgc。
光程长度L设定为与测量数据(λ1:αM1,λ2:αM2,λ3:αM3,…)一起从控制单元20被输入至数据处理单元30。然而,可以不从控制单元20输入光程长度L。
例如,在光程设定为横跨眼睛前房的情况下的光程长度L可以设定为标准人的眼睛前房中的预定光程长度。在光程长度L具有固定值的情况下,可以预先包括浓度计算单元101。
此外,光程长度L可以设定为在医疗机构等中测量的眼睛前房中的光程长度。在上述情况中,可以从UI单元40输入光程长度L。浓度计算单元101可以从UI单元40在预先设置的多个光程长度中选择接近期望光程长度的光程长度。
以该方式,在从系统的外部输入光程长度L的情况下,数据处理单元30可以不经由控制单元20获得光程长度L。
<实例2>
图14是示出了实例2中的数据处理单元30的功能单元的视图。
数据处理单元30包括浓度计算单元111、公式存储单元112和固有值存储单元113。
公式存储单元112存储例如用于计算测量目标物质13中所包含且想要得到浓度的光学活性物质的浓度C的公式(6)、公式(10)和公式(11)等公式(算法)。公式存储单元112可以存储多个公式。
固有值存储单元113存储光学活性物质的固有值,这些光学活性物质的固有值被代入到计算光学活性物质的浓度C的公式中。固有值存储单元113可以存储分别与多个公式对应的多个固有值。
浓度计算单元111将已从控制单元20发送的测量数据(λ1:αM1,λ2:αM2,λ3:αM3,…)和光程长度L以及从固有值存储单元113读取出的固有值代入到从公式存储单元112读取的公式中。通过组合波长λ与观测旋光度αM来得到测量数据。浓度计算单元111通过使用代入结果来计算预定的光学活性物质的浓度C。
在UI单元40的显示器等中显示计算出的光学活性物质的浓度C。
可以从控制单元20输入光程长度L。然而,与实例1类似,可以包括数据处理单元30。此外,数据处理单元30可以不经由控制单元20获得光程长度L。
图15是示出了存储在公式存储单元112中的公式以及存储在固有值存储单元113中的固有值的实例的视图。图15中的(a)示出了由公式存储单元112存储的公式(算法)的实例。图15中的(b)示出了由固有值存储单元113存储的固有值的实例。这里,通过使用房水作为测量目标物质13的实例来进行描述。
首先,将描述图15中的(a)所示的公式存储单元112。
作为实例,公式存储单元112存储多个公式#1至#5。公式存储单元112设定为存储公式,但可以是执行公式的程序。即使在该情况下,也描述公式存储单元112存储公式的情况。
如上所述,公式根据想要显示在UI单元40中的光学活性物质(第一光学活性物质)以及房水中所包含的除了第一光学活性物质之外的光学活性物质(第二光学活性物质)的组合来设定浓度。因此,利用第一光学活性物质和第二光学活性物质的组合来选择公式。
在葡萄糖被设定成第一光学活性物质而除了葡萄糖之外的其余光学活性物质的集合被设定成第二光学活性物质的#1情况下,选择公式(11)。
在白蛋白被设定作为第一光学活性物质而除了白蛋白之外的其余光学活性物质的集合被设定作为第二光学活性物质的#2情况下,选择公式(12)。
在葡萄糖和白蛋白被设定作为第一光学活性物质而除了葡萄糖和白蛋白之外的其余光学活性物质的集合被设定作为第二光学活性物质的#3情况下,选择公式(6)。
在葡萄糖被设定作为第一光学活性物质而包含白蛋白以及除了葡萄糖和白蛋白之外的其余光学活性物质在内的光学活性物质被设定成第二光学活性物质的#4情况下,选择公式(6)。
在葡萄糖被设定作为第一光学活性物质而白蛋白被设定作为第二光学活性物质的#5情况下,选择公式(10)。
公式存储单元112可以存储其他组合。
接下来,将描述图15中的(b)所示的固有值存储单元113。
例如,固有值存储单元113存储葡萄糖的固有值Ag和λg以及白蛋白的固有值λa。因此,根据选自公式存储单元112的公式来读取需要的固有值。除了存储λa之外,还可以存储公式(6)中的Aa等作为白蛋白的固有值。
图16是示出了实例2中的数据处理单元30的硬件构造的实例的视图。
数据处理单元30包括中央处理单元(下文中将描述成CPU)31、RAM32、ROM33、硬盘驱动器(将在下文中描述成HDD)34、输入和输出接口(将在下文中描述成I/O)35。
CPU31、RAM32、ROM33、HDD34和I/O35经由单根总线36彼此连接。
在光学活性物质的浓度计算系统中,I/O35与控制单元20和UI单元40连接。
可以使用ASIC37等构造CPU31、RAM32、ROM33和I/O35。
浓度计算单元111与CPU31对应。公式存储单元112和固有值存储单元113对应于ROM33或HDD34。
CPU31包括执行逻辑运算和算术运算等的算术逻辑单元(ALU)。
随机存取存储器(RAM)32保持CPU31执行运算时所使用的程序或数据。一般来说,RAM32是易失性存储介质。在提供电力期间,RAM32可以读写数据并保持已写数据,而如果不提供电力,则失去数据。
ROM33是非易失性存储介质(非易失性存储器)。在提供电力期间,ROM33可以读写数据并保持已写数据,并且当不提供电力时,ROM33也可以保持已写数据。这里,ROM33可以是不能够重写数据的只读存储器(ROM),也可以是能够重写数据的闪速存储器。
HDD34是可以重写大量数据的非易失性存储器。HDD34保持程序或数据。ROM33可以保持作为固件的程序或数据。在该情况下,数据处理单元30可以不包括HDD34。
如果电力被输入到光学活性物质的浓度计算系统1,则数据处理单元30的CPU31将程序(公式)和数据(固有值)从HDD34或ROM33加载至RAM32。
CPU31将测量数据写入到RAM32中,将已加载的固有值和测量数据(λ1:αM1,λ2:αM2,λ3:αM3,…)代入到已加载的程序(公式),并计算光学活性物质的浓度C。通过组合从控制单元20发送的波长λ和使用波长λ测得的观测旋光度αM来得到测量数据。
然后,在UI单元40中所包含的显示器等中显示得到的光学活性物质的浓度C。
<实例3>
图17是示出了实例3中的数据处理单元30的功能单元的视图。
数据处理单元30包括浓度输出单元121和查找表(LUT)存储单元122。LUT存储单元122是存储单元的实例。
LUT存储单元122存储参照表(LUT)(将在下文中进行描述)。
基于通过组合从控制单元20发送来的波长λ和观测旋光度αM2而得到的测量数据(λ1:αM1,λ2:αM2,λ3:αM3,…),数据处理单元30中的浓度输出单元121参考由LUT存储单元122存储的LUT将与波长λ和观测旋光度αM的组合匹配的光学活性物质的浓度C输出至UI单元40。
在LUT中没有与波长λ和观测旋光度αM的组合匹配的组合的情况下,例如可以选择最接近的组合,或可以使用多个组合来执行补充计算。
浓度输出单元121可以具有如下构造:提取出与由LUT存储单元122存储的波长λ和观测旋光度αM的组合匹配的光学活性物质的浓度C。浓度输出单元121不需要进行计算光学活性物质的浓度C的操作的性能。
LUT存储单元122可以是能够存储LUT的存储单元,可以是例如ROM和HDD等非易失性存储器,并可以是能够重写的非易失性存储器。
图18是示出了由LUT存储单元122存储的LUT的实例的视图。这里,假设得到了房水中所包含的葡萄糖浓度Cg。在测量中使用波长λ1、λ2和λ3这三个波长。
如图18所示,葡萄糖浓度Cgc相关联地被存储在由LUT存储单元122存储的LUT中。在假设测量目标物质13被光照射的情况下,预先根据波长λ和观测旋光度αM的组合来计算葡萄糖浓度Cgc。因此,如果输入通过组合波长λ和观测旋光度αM而得到的测量数据(λ1:αM1,λ2:αM2和λ3:αM3),则数据处理单元30中的浓度输出单元121参考LUT存储单元122的LUT获得与测量数据(λ1:αM1,λ2:αM2和λ3:αM3)的组合匹配或与接近组合对应的葡萄糖浓度Cgc。浓度输出单元121将获得的葡萄糖浓度Cgc输出至UI单元40。也就是说,如果输入通过组合波长λ和观测旋光度αM而得到的测量数据,则浓度输出单元121从LUT存储单元122获得与该输入对应的葡萄糖浓度Cgc,并将获得的葡萄糖浓度Cgc输出至UI单元40。例如,在浓度计算系统的设计步骤或制造步骤中预先计算LUT存储单元122中所存储的数据并进行存储。
在使用LUT存储单元122的实例3中,因为可以读取预先计算出且被存储在LUT存储单元122中的葡萄糖浓度Cgc,所以与未使用LUT存储单元122的实例1或实例2相比,不需要根据输入的测量数据(λ1:αM1,λ2:αM2和λ3:αM3)来计算和算出葡萄糖浓度Cg。因此,与根据实例1或实例2所述且被输入的测量数据(λ1:αM1,λ2:αM2和λ3:αM3)来计算并输出葡萄糖浓度Cg的情况相比,输入测量数据(λ1:αM1,λ2:αM2和λ3:αM3)之后输出葡萄糖浓度Cg为止的时间是短的。
可以如下所述那样制作这种LUT。也就是说,提供如下处理:获得分别与多个光线的波长对应的旋光度的组合(λ1:αM1,λ2:αM2,λ3:αM3,…)的处理;对旋光度的每个组合而言,使用最小二乘法根据公式以及获得的旋光度的组合来计算想要得到的光学活性物质(第一光学活性物质)的浓度的处理,上述公式包括表示第一光学活性物质中的旋光度的波长相关性的第一函数以及表示其他光学活性物质(第二光学活性物质)中的旋光度的波长相关性的第二函数;以及将第一光学活性物质的彼此相关的旋光度的组合和浓度存储在存储单元中的处理,第一光学活性物质的浓度与旋光度的组合对应。在第一函数中,想要得到的光学活性物质的浓度具有未知值,并且用于限定第一光学活性物质的旋光色散特性的固有值被设定成已知值或一定受限范围内的未知值。在第二函数中,用于限定其他光学活性物质的旋光色散特性的固有值被设定成未知值。
在制作中,不仅可以使用最小二乘法,而且还可以使用对由数量与公式中所包含的未知值的数量相同的等式构成的联立方程求解的方法。
在波长的数量设定成m并且对于波长设定了n个观测旋光度αM的情况下,由LUT存储单元122存储的LUT存储通过组合(求和)上述常量而得到的m×n×n个光学活性物质的浓度C。
LUT存储单元122可以存储多个LUT。
可以为固定的光程长度L(例如,标准人的眼睛前房中的预定的光程长度)设置从由LUT存储单元122存储的LUT计算出的葡萄糖浓度Cgc。在该情况下,获得了考虑光程长度L的葡萄糖浓度Cgc。
可以为单位光程长度设置从由LUT存储单元122存储的LUT计算出的葡萄糖浓度Cgc。在该情况下,从UI单元40等输入光程长度L,并对光程长度L执行修正。因为用于光程长度L的修正是简单的乘法,所以浓度输出单元121可以执行处理。与根据输入的测量数据(λ1:αM1,λ2:αM2和λ3:αM3)计算葡萄糖浓度Cgc相比,短时间地进行处理。因此,在输入测量数据(λ1:αM1,λ2:αM2和λ3:αM3)之后输出葡萄糖浓度Cgc为止的期间所产生的影响是小的。
(光学活性物质的浓度C的显示实例)
在光学活性物质的浓度计算系统1中,将描述在UI单元40中所包含的显示器显示光学活性物质的浓度C的显示实例。这里,将描述显示计算出的葡萄糖浓度Cgc的实例。
图19是示出了在对葡萄糖浓度测量装置应用光学活性物质的浓度计算系统1的情况下在UI单元40中所包含的显示器显示光学活性物质的浓度C的显示实例的视图。图19中的(a)示出了以数值进行显示的显示实例。图19中的(b)示出了以条形图(水平)进行显示的显示实例。图19中的(c)示出了以标记“OK/NG”进行显示的显示实例。图19中的(d)示出了使用数值除了显示葡萄糖浓度Cg之外还显示其他光学活性物质的浓度的显示实例。
因为图19中的(a)和图19中的(d)用数值示出,所以可以容易地追踪葡萄糖浓度Cgc的变化。因为图19中的(b)用条形图(水平)示出,所以一瞥之下可以识别出葡萄糖浓度Cgc的水平。在图19中的(c)中,在数值等于或小于预定阈值的情况下设定“OK”,在数值大于预定阈值的情况下设定“NG”,并且在“NG”的情况下警告使用者。这时,“OK”可以用绿色显示而“NG”可以用红色显示。
显示实例仅是实例,并可以被组合。此外,可以应用其他显示方法。UI单元40的输出不限于例如显示器等显示视听信息的装置(即,例如提供音频信息的装置等装置)。
在为了测量血糖水平目的的情况下,可以根据血液中的葡萄糖浓度与房水中的葡萄糖浓度之间的相关性将在房水中测量出(计算出)的葡萄糖浓度Cgc转换成血液中的葡萄糖浓度并显示。
在以上描述中,描述了房水中的光学活性物质的浓度C,尤其是葡萄糖浓度Cg。
以上实施例除了应用于人的眼球房水之外还可以应用于使用光学池等的测量装置。可以提供如下的光学系统:包含光学活性物质的测量目标物质13被放入到光学池(容器)中,并且光横跨且穿过该光学池。在该情况下,光学池决定光程长度L。
由数据处理单元30执行并用于计算光学活性物质的浓度C的算法可以用作例如由通用计算机(PC)执行的程序。
例如,个人计算机(PC)的使用者在PC中所包含的显示器(显示屏幕)中指定用于计算光学活性物质的浓度C的公式、固有值和光程长度L,以及用于计算光学活性物质的浓度C的波长λ和旋光度α的至少一个组合等。如果进行了指定,则PC执行程序,并因此使用最小二乘法计算光学活性物质的浓度C,并在显示器(显示屏幕)上显示计算出的光学活性物质的浓度C。执行程序时所使用的固有值、光程长度L以及波长λ和旋光度α的组合等不一定使用真实的测量结果。例如,可以通过利用光实际照射测量目标物质来获得旋光度α。此外,也就是说,在没有用光实际照射测量目标物质而假设用光进行照射的情况下,旋光度α可以是假设的旋光度。以这种方式,在由使用者假设设定的条件下计算浓度,因此,可以在各种仿真等中使用旋光度α。可以使用对由数量与公式中所包含的未知值的数量相同的等式构成的联立方程求解的方法替代最小二乘法。
工业实用性
根据本发明的浓度计算系统等在计算测量目标物质中所包含的例如葡萄糖或白蛋白等光学活性物质的浓度时是有用的。
参考具体实施例详细地描述了本发明。然而,对本领域的技术人员显而易见的是,可以在不脱离本发明的主旨和范围的情况下进行各种修改和变型。
本发明基于且要求2014年3月20日提交的日本专利申请公开No.2014-059221和2015年3月9日提交的日本专利申请公开No.2015-046453的优先权;上述日本专利申请公开的全部内容通过引用并入本文中。
附图标记列表
1光学活性物质的浓度计算系统
10测量单元
11光源
12偏光器
13测量目标
14补偿器
15分析器
16光接收元件
20控制单元
30数据处理单元
40UI单元
101,111浓度计算单元
112公式存储单元
113固有值存储单元
121浓度输出单元
122LUT存储单元
α旋光度
[α]旋光率
αM观测旋光度
Claims (10)
1.一种光学活性物质的浓度计算系统,其根据表示旋光度的公式计算光学活性物质的浓度,其中
所述公式包括:第一函数,其表示第一光学活性物质的旋光度的波长相关性;以及第二函数,其表示第二光学活性物质的旋光度的波长相关性,
在所述第一函数中,所述第一光学活性物质的浓度具有未知值,并且用于限定所述第一光学活性物质的旋光色散特性的固有值是已知值或是一定受限范围内的未知值,
在所述第二函数中,用于限定所述第二光学活性物质的旋光色散特性的固有值是未知值,并且
使用最小二乘法根据所述公式以及测量目标的分别与多个波长对应的旋光度来计算所述第一光学活性物质的所述浓度。
2.根据权利要求1所述的光学活性物质的浓度计算系统,其中
所述第二光学活性物质包含多种光学活性物质,
所述第二函数是如下的函数:以单个项表示所述多种光学活性物质的旋光度的波长相关性的集合,并且
所述第二函数是如下的函数:所述单个项中用于限定所述多种光学活性物质的旋光色散特性的固有值用未知值来表示。
3.根据权利要求1或2所述的光学活性物质的浓度计算系统,其中
在所述第一函数中用于限定所述第一光学活性物质的旋光色散特性的至少一些固有值是一定受限范围内的未知值。
4.根据权利要求1至3中任一项所述的光学活性物质的浓度计算系统,其中
在所述第二函数中用于所述第二光学活性物质的至少一些固有值是一定受限范围内的未知值。
5.根据权利要求1至4中任一项所述的光学活性物质的浓度计算系统,其中
在改变所述第一函数中所包含的未知值和所述第二函数中所包含的未知值中的至少一个未知值的初始值的同时,根据最小二乘法多次计算所述第一光学活性物质的所述浓度,并且
当多次计算所述第一光学活性物质的所述浓度时,计算出的所述第一光学活性物质的与使目标函数的值最大的初始值对应的浓度不作为所述第一光学活性物质的所述浓度,所述目标函数用所述最小二乘法中的旋光度的理论值与获得的旋光度之差的平方和表示。
6.根据权利要求1至4中任一项所述的光学活性物质的浓度计算系统,其中
在对为所述第一函数中所包含的未知值限制一定范围的限制条件以及为所述第二函数中所包含的未知值限制一定范围的限制条件中的至少一个限制条件进行改变的同时,根据最小二乘法多次计算所述第一光学活性物质的所述浓度,并且
当多次计算所述第一光学活性物质的所述浓度时,计算出的所述第一光学活性物质的与使目标函数的值最大的限制条件对应的浓度不作为所述第一光学活性物质的所述浓度,所述目标函数用所述最小二乘法中的旋光度的理论值与获得的旋光度之差的平方和表示。
7.根据权利要求1至6中任一项所述的光学活性物质的浓度计算系统,其中
对传输通过人的眼球中的房水的光而言,获得所述测量目标的分别与所述多个波长对应的旋光度。
8.一种光学活性物质的浓度计算系统的制作方法,所述方法包括:
获得分别与多个光线的波长对应的旋光度的组合;
对旋光度的每个组合而言,使用最小二乘法根据公式和获得的旋光度的组合来计算第一光学活性物质的浓度,所述公式包括:第一函数,其表示所述第一光学活性物质的旋光度的波长相关性;以及第二函数,其表示第二光学活性物质的旋光度的波长相关性,以及
将所述第一光学活性物质的彼此相关的旋光度的组合和浓度存储在存储单元中,所述第一光学活性物质的所述浓度与旋光度的组合对应,其中
在所述第一函数中,所述第一光学活性物质的所述浓度具有未知值,并且用于限定所述第一光学活性物质的旋光色散特性的固有值是已知值或是一定受限范围内的未知值,并且
在所述第二函数中,用于限定所述第二光学活性物质的旋光色散特性的固有值是未知值。
9.一种计算机可读介质,其是非临时性的并存储使计算机执行根据表示旋光度的公式计算光学活性物质的浓度的功能的程序,其中
所述公式包括:第一函数,其表示第一光学活性物质的旋光度的波长相关性;以及第二函数,其表示第二光学活性物质的旋光度的波长相关性,
在所述第一函数中,所述第一光学活性物质的浓度具有未知值,并且用于限定所述第一光学活性物质的旋光色散特性的固有值是已知值或是一定受限范围内的未知值,
在所述第二函数中,用于限定所述第二光学活性物质的旋光色散特性的固有值是未知值,
使用最小二乘法根据所述公式以及测量目标的分别与多个波长对应的旋光度来计算所述第一光学活性物质的所述浓度。
10.一种使计算机执行根据表示旋光度的公式计算光学活性物质的浓度的功能的程序,其中
所述公式包括:第一函数,其表示第一光学活性物质的旋光度的波长相关性;以及第二函数,其表示第二光学活性物质的旋光度的波长相关性,
在所述第一函数中,所述第一光学活性物质的浓度具有未知值,并且用于限定所述第一光学活性物质的旋光色散特性的固有值是已知值或是一定受限范围内的未知值,
在所述第二函数中,用于限定所述第二光学活性物质的旋光色散特性的固有值是未知值,
使用最小二乘法根据所述公式以及测量目标的分别与多个波长对应的旋光度来计算所述第一光学活性物质的所述浓度。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014-059221 | 2014-03-20 | ||
JP2014059221 | 2014-03-20 | ||
JP2015-046453 | 2015-03-09 | ||
JP2015046453A JP5783342B1 (ja) | 2014-03-20 | 2015-03-09 | 光学活性物質の濃度算出システム、光学活性物質の濃度算出システムの製造方法及びプログラム |
PCT/JP2015/058350 WO2015141804A1 (ja) | 2014-03-20 | 2015-03-19 | 光学活性物質の濃度算出システム、光学活性物質の濃度算出システムの製造方法、コンピュータ読み取り媒体及びプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105814435A true CN105814435A (zh) | 2016-07-27 |
CN105814435B CN105814435B (zh) | 2018-12-28 |
Family
ID=54144766
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580003067.3A Active CN105814435B (zh) | 2014-03-20 | 2015-03-19 | 光学活性物质的浓度计算系统、光学活性物质的浓度计算系统的制作方法、计算机可读介质及程序 |
Country Status (4)
Country | Link |
---|---|
US (1) | US9851293B2 (zh) |
JP (1) | JP5783342B1 (zh) |
CN (1) | CN105814435B (zh) |
WO (1) | WO2015141804A1 (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10702197B2 (en) * | 2016-07-26 | 2020-07-07 | The Texas A&M University System | Dual amplitude modulation and polarization frequency modulation as well as compensation for noninvasive glucose monitoring |
JP6841146B2 (ja) | 2017-04-21 | 2021-03-10 | 富士ゼロックス株式会社 | 眼球の光計測装置及び眼球の光計測方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4128458A1 (de) * | 1991-08-28 | 1993-03-04 | Siemens Ag | Verfahren und vorrichtung zur quantitativen bestimmung optisch aktiver substanzen |
JP2002082047A (ja) * | 2000-09-08 | 2002-03-22 | Matsushita Electric Ind Co Ltd | 尿検査装置 |
WO2007044486A1 (en) * | 2005-10-13 | 2007-04-19 | Baylor University | Spectroscopic determination of sucrose |
CN201147318Y (zh) * | 2007-10-12 | 2008-11-12 | 中国人民解放军第三军医大学野战外科研究所 | 正交双偏振光无创血糖测量装置 |
JP2012112907A (ja) * | 2010-11-26 | 2012-06-14 | Global Fiber Optics Co Ltd | 旋光成分分析装置および旋光成分分析方法ならびに旋光度の温度特性または波長特性測定装置 |
CN102755167A (zh) * | 2011-04-29 | 2012-10-31 | 财团法人工业技术研究院 | 非侵入式血糖监测装置与方法以及生化分子的分析方法 |
JP2014032146A (ja) * | 2012-08-06 | 2014-02-20 | Seiko Epson Corp | 濃度測定装置及び濃度測定装置の制御方法 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3481671A (en) * | 1966-07-27 | 1969-12-02 | Melpar Inc | Apparatus and method for obtaining optical rotatory dispersion measurements |
JP3014633B2 (ja) | 1995-11-16 | 2000-02-28 | 松下電器産業株式会社 | 尿検査方法 |
EP1300670A3 (en) * | 1995-11-16 | 2003-10-29 | Matsushita Electric Industrial Co., Ltd. | Method of measuring angle of rotation and polarimeter |
EP0880688B1 (en) * | 1996-01-31 | 2003-08-27 | Gary W. Yanik | Optical activity detector for use with optically active compounds |
US5788632A (en) * | 1996-03-19 | 1998-08-04 | Abbott Laboratories | Apparatus and process for the non-invasive measurement of optically active compounds |
US6327037B1 (en) * | 1997-11-12 | 2001-12-04 | Chien Chou | Optical rotation angle polarimeter |
US6620622B1 (en) * | 1999-06-29 | 2003-09-16 | Matsushita Electric Industrial Co., Ltd. | Method of polarimetry and method of urinalysis using the same |
JP4507372B2 (ja) * | 2000-09-08 | 2010-07-21 | パナソニック株式会社 | 尿検査装置 |
US6885882B2 (en) * | 2002-05-28 | 2005-04-26 | Cote Gerard L. | Method and apparatus for non-invasive glucose sensing through the eye |
JP5006049B2 (ja) * | 2004-01-13 | 2012-08-22 | ザ・ユニバーシティ・オブ・トレド | 非侵襲複屈折補償感知旋光計 |
WO2008058184A2 (en) * | 2006-11-07 | 2008-05-15 | Deka Products Limited Partnership | Method and apparatus for determining concentration using polarized light |
JP2008209204A (ja) * | 2007-02-26 | 2008-09-11 | Citizen Holdings Co Ltd | 濃度測定装置 |
JP2017023328A (ja) * | 2015-07-21 | 2017-02-02 | 富士ゼロックス株式会社 | 光学活性物質の濃度算出システム及びプログラム |
-
2015
- 2015-03-09 JP JP2015046453A patent/JP5783342B1/ja not_active Expired - Fee Related
- 2015-03-19 CN CN201580003067.3A patent/CN105814435B/zh active Active
- 2015-03-19 WO PCT/JP2015/058350 patent/WO2015141804A1/ja active Application Filing
-
2016
- 2016-04-26 US US15/138,646 patent/US9851293B2/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4128458A1 (de) * | 1991-08-28 | 1993-03-04 | Siemens Ag | Verfahren und vorrichtung zur quantitativen bestimmung optisch aktiver substanzen |
JP2002082047A (ja) * | 2000-09-08 | 2002-03-22 | Matsushita Electric Ind Co Ltd | 尿検査装置 |
WO2007044486A1 (en) * | 2005-10-13 | 2007-04-19 | Baylor University | Spectroscopic determination of sucrose |
CN201147318Y (zh) * | 2007-10-12 | 2008-11-12 | 中国人民解放军第三军医大学野战外科研究所 | 正交双偏振光无创血糖测量装置 |
JP2012112907A (ja) * | 2010-11-26 | 2012-06-14 | Global Fiber Optics Co Ltd | 旋光成分分析装置および旋光成分分析方法ならびに旋光度の温度特性または波長特性測定装置 |
CN102755167A (zh) * | 2011-04-29 | 2012-10-31 | 财团法人工业技术研究院 | 非侵入式血糖监测装置与方法以及生化分子的分析方法 |
JP2014032146A (ja) * | 2012-08-06 | 2014-02-20 | Seiko Epson Corp | 濃度測定装置及び濃度測定装置の制御方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2015141804A1 (ja) | 2015-09-24 |
CN105814435B (zh) | 2018-12-28 |
US20160238518A1 (en) | 2016-08-18 |
US9851293B2 (en) | 2017-12-26 |
JP5783342B1 (ja) | 2015-09-24 |
JP2015194482A (ja) | 2015-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chatani et al. | Water molecular system dynamics associated with amyloidogenic nucleation as revealed by real time near infrared spectroscopy and aquaphotomics | |
Gil | Review on Mueller matrix algebra for the analysis of polarimetric measurements | |
EP3403067A1 (en) | System and method for multi-parameter spectroscopy | |
Arteaga et al. | Pseudopolar decomposition of the Jones and Mueller-Jones exponential polarization matrices | |
TW201643408A (zh) | 對位裝置及對位方法 | |
Van Eeckhout et al. | Depolarizing metrics for plant samples imaging | |
JP2007040805A (ja) | 撮像偏光計測方法 | |
US20180228415A1 (en) | Method and system for sensing glucose concentration | |
Pham et al. | Extraction of effective parameters of anisotropic optical materials using a decoupled analytical method | |
US20170020385A1 (en) | Concentration calculation system of optically active substance and non-transitory computer readable medium storing program | |
Luber | Solvent effects in calculated vibrational Raman optical activity spectra of α-helices | |
TWI810247B (zh) | 用於校準和使用攝影機來偵測樣品中的分析物之方法及系統 | |
Le Gratiet et al. | Review on complete Mueller matrix optical scanning microscopy imaging | |
Bakshi et al. | Fourier transform infrared spectroscopy of peptides | |
WO2015154105A1 (en) | Device for non-invasive determination of glucose concentration in human organism | |
Ma et al. | Understanding the effect of urea on the phase transition of poly (N-isopropylacrylamide) in aqueous solution by temperature-dependent near-infrared spectroscopy | |
CN105814435A (zh) | 光学活性物质的浓度计算系统、光学活性物质的浓度计算系统的制作方法、计算机可读介质及程序 | |
Wang et al. | Principle and implementation of stokes vector polarization imaging technology | |
Graziano et al. | On the intactness of hydrogen bonds around nonpolar solutes dissolved in water | |
Leng et al. | Mueller matrix imaging with a spatially modulated polarization light source | |
US11607136B2 (en) | System and method for conformal vision | |
CN114305406B (zh) | 浓度感测系统与方法 | |
Otsuki | Multiple scattering in turbid media containing chiral components: a Monte Carlo simulation | |
CN110613461A (zh) | 葡萄糖浓度的感测系统与感测方法 | |
JP2007286011A (ja) | 光学特性計測装置および方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: Tokyo Patentee after: Fuji film business innovation Co.,Ltd. Address before: Tokyo Patentee before: Fuji Xerox Co.,Ltd. |