CN101980656B - 用于重构混浊介质内部的荧光图像的方法以及用于对混浊介质内部成像的设备 - Google Patents
用于重构混浊介质内部的荧光图像的方法以及用于对混浊介质内部成像的设备 Download PDFInfo
- Publication number
- CN101980656B CN101980656B CN2009801110542A CN200980111054A CN101980656B CN 101980656 B CN101980656 B CN 101980656B CN 2009801110542 A CN2009801110542 A CN 2009801110542A CN 200980111054 A CN200980111054 A CN 200980111054A CN 101980656 B CN101980656 B CN 101980656B
- Authority
- CN
- China
- Prior art keywords
- turbid medium
- fluorescence
- light
- contrast agent
- wave
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
- A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
- A61B5/4312—Breast evaluation or disorder diagnosis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0082—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes
- A61B5/0091—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes for mammography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N21/4795—Scattering, i.e. diffuse reflection spatially resolved investigating of object in scattering medium
-
- G—PHYSICS
- 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/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
- G01N21/6428—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Veterinary Medicine (AREA)
- Biochemistry (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Analytical Chemistry (AREA)
- Biomedical Technology (AREA)
- Gynecology & Obstetrics (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Reproductive Health (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本发明提供了一种用于重构混浊介质内部的荧光图像的方法。该方法包括步骤:在测量体(4)中容纳向其施用了荧光对比剂的混浊介质(1)。荧光对比剂能够在用光照射时发射第一波长范围内的光。该方法还包括:通过随后利用来自多个不同源位置的光照射混浊介质(1)并且对于每个源位置在多个检测位置检测从混浊介质(1)发出的光,在多个不同的波长(λ1,...,λk)下执行衰减测量;根据衰减测量对于所述多个不同的波长重构作为混浊介质(1)内部的位置的函数的吸收特性(μa(r,λ));对于第一波长范围的波长计算作为混浊介质(1)内部的位置的函数的吸收特性;通过随后利用来自所述多个源位置的使得荧光对比剂发射第一波长范围内的光的光照射混浊介质(1)并且对于每个源位置在所述多个检测位置检测从荧光对比剂发出的光,执行荧光测量;以及根据荧光测量使用计算的吸收特性重构混浊介质(1)内部荧光对比剂的空间分布的荧光图像。
Description
技术领域
本发明涉及用于重构混浊介质内部的荧光图像的方法以及用于对混浊介质的内部成像的设备。
背景技术
在本申请的上下文中,术语混浊介质应当被理解为表示由具有高光散射系数的材料(例如脂肪乳溶液或生物组织)组成的物质。此外,光应当被理解为表示电磁辐射,特别是具有从180nm到1400nm的范围内的波长的电磁辐射。术语“光学特性”覆盖减小的散射系数μ′s和吸收系数μa。此外,“匹配光学特性”应当被理解为具有类似的减小的散射系数μ′s和类似的吸收系数μa。
用于对混浊介质的内部成像(例如用于乳腺癌筛查)的方法是通过使用光,特别是使用近红外(NIR)内的光成像,该方法近年来变得流行。这样的方法在乳腺X光检查设备以及用于检查人体或动物体的其他部分的设备中实现。用于借助于光对混浊介质内部成像的这种方法的突出的实例是漫射光学层析术(DOT)。例如,用于对混浊介质内部成像的这种DOT设备使用照射混浊介质的光源以及用于测量通过混浊介质传输的光的一部分(即其强度)的光电检测器。控制单元被提供以用于控制扫描过程。处理单元被提供以用于根据测量的强度重构混浊介质内部的图像。一些已知的设备特别适于检查女性胸部。为了允许检查混浊介质,设备设有包含测量体并且被设置成接纳混浊介质的接纳部分。来自光源的光耦合到接纳体中并且耦合到混浊介质中。光被选择成使得它能够传播通过混浊介质。为了对女性胸部内部成像,典型地使用NIR(近红外)内的光。作为将光耦合到接纳体中的结果从混浊介质发出的散射光耦合出接纳体。耦合出接纳体的光用来重构混浊介质内部的图像。由于要检查的混浊介质的不同尺寸的原因,接纳部分的尺寸可能不完美地与混浊介质的尺寸匹配,即接纳体的边界与混浊介质之间留有空间。研究中的混浊介质的部分被填充接纳体中的空间的散射介质(耦合介质)包围。该散射介质被选择成使得散射介质的光学参数(例如吸收和散射系数)类似于混浊介质的相应光学参数。光源随后从不同的方向照射混浊介质并且光电检测器测量透射通过混浊介质的光的部分。利用从不同的方向定向到混浊介质的光执行多个这样的测量,并且基于这些测量结果(即获得的数据集),处理单元重构被检查的混浊介质的图像。
依照该方法的一个发展,执行光的衰减扫描,其中对于源位置和检测位置的多个组合检测光的衰减。在这些测量中,使用了混浊介质的内在对比度,即由于诸如氧合血红蛋白、脱氧血红蛋白、水和脂质之类的散射体和发色团的存在,不同波长的光按不同的量衰减。根据这些衰减扫描,可以重构混浊介质的吸收图像以及诸如血红蛋白浓度之类的生理参数的图像。这种技术称为漫射光学层析术(DOT)。
依照该方法的另一发展,为了所述测量施用优先地积聚在所研究的混浊介质中的病变(例如女性胸部中的癌组织)处的荧光对比剂。利用来自优选地为激光器的光源的光照射混浊介质,并且检测混浊介质发射的荧光。根据该测量,重构胸部的荧光发射的体图像,即使用外源对比度。这样,重构混浊介质中的对比剂的空间分布。这种方法称为漫射光学荧光层析术。
已经开发了适于执行衰减测量和荧光测量两者的设备。在这样的设备中,对于特定波长范围内的多个波长执行衰减测量。衰减测量中获得的数据允许重构对于所述多个波长的混浊介质的吸收特性的空间分辨图像。因此,可以根据衰减测量重构这些波长下光如何传播通过混浊介质。为了执行荧光测量,明确地激发荧光对比剂并且仅仅测量荧光对比剂发射的光。这例如通过在测量体与检测单元之间的光路上引入适当的滤光器来实现。然而,荧光对比剂发射的光包括位于一定波长范围内的波段,其可能与衰减测量探测的波长不同。因此,对于荧光测量的相关波长而言,出现了不知道光如何与混浊介质相互作用的问题。
在申请人已知的一种方法中,使用了以下近似:光与混浊介质的相互作用对于与衰减测量相关的波长以及对于与荧光测量相关的波长范围基本上相同。然而,该近似通常不成立并且该近似的质量取决于例如不同的患者以及他们的组织的成分而不同。结果,用于重构混浊介质内部的荧光图像,即用于重构混浊介质中荧光对比度的空间分布的图像的当前的方法包含以下问题:根据重构图像不能令人满意地确定混浊介质中不同位置处的荧光对比剂的浓度。
Koenig.A等人(Fluorescence diffuse optical(fDOT)system for smallanimal studies.Proceedings of the 29th Annual International Conference ofthe IEEE EMBS CitéInternationale,Lyon,France,August 23-262007.pp.2626-2629)描述了一种用于荧光漫射光学层析术的系统,其中在荧光重构之前计算光学异质图。
发明内容
本发明的目的是提供一种用于重构混浊介质内部的荧光图像的方法以及一种用于对混浊介质内部成像的设备,利用其可以根据重构图像更精确地确定混浊介质中不同位置处的荧光对比剂的浓度。此外,这将快速而高效地实现,而无需附加的硬件。
这个目的是通过依照权利要求1的用于重构混浊介质内部的荧光图像的方法来获得的。该方法包括步骤:在测量体中容纳向其施用了荧光对比剂的混浊介质。荧光对比剂能够在用光照射时发射第一波长范围内的光。该方法还包括:通过随后利用来自多个不同源位置的光照射混浊介质并且对于每个源位置在多个检测位置检测从混浊介质发出的光,在多个不同的波长下执行衰减测量。该方法还包括:根据衰减测量,对于所述多个不同的波长重构作为混浊介质内部的位置的函数的吸收特性;对于第一波长范围的波长计算作为混浊介质内部的位置的函数的吸收特性;通过随后利用来自所述多个源位置的使得荧光对比剂发射第一波长范围内的光的光照射混浊介质并且对于每个源位置在所述多个检测位置检测从荧光对比剂发出的光,执行荧光测量;以及根据荧光测量使用计算的吸收特性重构混浊介质内部荧光对比剂的空间分布的荧光图像。结果,更正确地估计具有第一波长范围的波长的光与混浊介质的相互作用,并且可以以更高的精确度重构混浊介质内对比剂的浓度和空间分布。此外,对于离散的波长执行的衰减测量可以用来实现该更高的精确度,从而无需附加的硬件。“作为位置的函数的吸收特性”表示吸收特性依赖于混浊介质中的位置而确定。(数学意义上的)函数依赖性不一定必须确立。计算的吸收特性用作重构荧光图像的步骤中的输入信息。
优选地,在重构吸收特性的步骤中,计算吸收系数。在这种情况下,已知的吸收系数波长依赖性可以用来实现对于第一波长范围的波长的有利计算。
依照一个方面,通过将混浊介质看作多种物质的线性组合并且在所述多个不同的波长下根据衰减测量确定所述多种物质的局部浓度而计算吸收系数。在这种情况下,可以以令人满意的精确度计算关于第一波长范围的波长的吸收特性,从而可以在重构的荧光图像中更正确地可视化荧光对比剂的浓度。
如果使用所述多种物质的已知光谱行为计算对于第一波长范围的波长的吸收特性,那么可以以令人满意的精确度计算关于第一范围的波长的吸收特性。结果,重构的荧光图像提供了有关不同位置处荧光对比剂的浓度的可靠信息。因此,可以提供允许将恶性病变与其他异常区分的荧光图像。
如果所述多种物质包括荧光对比剂,那么即使荧光对比剂的浓度为高,也可以根据荧光测量重构提供有关荧光对比剂浓度的令人满意的信息的荧光图像。
优选地,将重构的荧光图像作为输入反馈到重构作为混浊介质内部的位置的函数的吸收特性的步骤,并且迭代地重复外插吸收特性并且重构荧光图像的步骤。在这种情况下,即使混浊介质中出现具有高浓度的位置,也可以以提高的精确度确定荧光对比剂的浓度。
所述目的进一步通过依照权利要求7的用于对混浊介质内部成像的设备来获得。该设备包括:测量体,其适于容纳向其施用荧光对比剂的混浊介质,所述荧光对比剂在照射时能够发射第一波长范围的光;光源单元,其适于随后利用来自多个不同源位置的光照射测量体;光源单元能够有选择地发射多个不同波长的光;以及检测单元,其适于在多个不同检测位置检测从测量体发出的光。该设备还包括控制和重构单元,该控制和重构单元适于控制设备以:通过随后利用来自多个不同源位置的光照射测量体在多个不同的波长下执行衰减测量并且对于每个源位置在多个检测位置处检测从测量体发出的光;通过随后利用来自所述多个源位置的使得荧光对比剂发射第一波长范围内的光的光照射测量体并且对于每个源位置在所述多个检测位置检测从荧光对比剂发出的光,执行荧光测量;根据衰减测量对于所述多个不同的波长重构作为测量体内的位置的函数的吸收特性;对于第一波长范围的波长计算作为测量体内的位置的函数的吸收特性;以及根据荧光测量使用计算的吸收特性重构测量体内荧光对比剂的空间分布的荧光图像。结果,更正确地估计具有第一波长范围的波长的光与混浊介质的相互作用,并且可以以更高的精确度重构混浊介质内对比剂的浓度和空间分布。此外,对于若干离散的波长执行的衰减测量可以用来实现该更高的精确度,从而无需附加的硬件。如果在设备中提供可在测量体与检测单元之间的光路上引入的滤光器以用于荧光测量,那么可以以有限数量的硬件实现重构。
优选地,所述设备是医疗图像采集设备。
附图说明
根据参照附图对于实施例的详细描述将得到本发明的另外的特征和优点。
图1示意性地示出了用于对混浊介质内部成像的设备的测量体。
图2示意性地示出了图1设备中的测量体、光源单元和检测单元的设置。
具体实施方式
第一实施例
现在将参照图1和图2描述本发明的实施例。依照该实施例的用于对混浊介质内部成像的设备是用于漫射光学层析术(DOT)的设备。特别地,该设备适于检查女性胸部。这种设备的整体构造在本领域中是已知的。该设备包括受检查的个人以卧姿躺在其上的床(未示出)。床内形成开口,测量体4在其下延伸。图1中示出了测量体4。
在图1所示的设备中,要检查的混浊介质1是女性胸部。测量体4以适于接纳混浊介质1的接纳部分2为界,如图1中示意性所示。接纳部分2具有杯状形状并且设有开口3。由图1可见,要检查的混浊介质1置于测量体4中,使得它从开口3的侧面自由地悬在测量体4中。接纳部分2用来放置和稳定被检查的混浊介质1。
接纳部分2面向混浊介质3的内表面设有多个由光学引导纤维形成的如图2中示意性所示连接到光源单元6和检测单元7的光导5的末端。光源单元6包括至少一个能够有选择地发射多个不同波长的光的光源。优选地,光源单元6包括多个激光器,每个激光器发射单色的但是不同波长的光。在这种情况下,光源单元6适于使得一次来自这些光源之一的光可以有选择地耦合到测量体中。光导5的末端分布在接纳部分2的内表面上。所述设备适于使得来自光源单元6的光可以从多个不同方向(源位置)定向到混浊介质1,并且从混浊介质1发出的光可以在围绕测量体4分布的多个不同检测位置处由检测单元7检测。在该实施例中,检测单元7由多个检测器实现,其相应光导5分布在接纳部分2的内表面上。接纳部分2内表面处光导5的末端形成多个源位置和多个检测位置。在该实施例中,源位置的总数等于检测位置的总数;然而,本发明并不限于相等的数量。例如,在依照该实施例的设备中,在接纳部分2的内表面上提供了256个不同的源位置,并且提供了256个检测位置,即光导5的各末端。来自光源的光随后从这256个源位置定向到混浊介质1,并且对于每个源位置在所述256个检测位置处检测从混浊介质1发出的光。然而,本发明并不限于这些特定的数量。
如图2中示意性所示,所述设备包括控制和处理单元8,其用于控制图像的采集并且重构混浊介质1内部的图像。控制和处理单元8基于来自检测单元7的信号重构混浊介质1内部的图像。为了进行重构,使用了扫描期间采样的信号,其中光从不同的方向定向到混浊介质1。为了简单起见,不再描述用于对混浊介质内部成像的设备的这些在本领域中已知的元件。
接纳部分2还被结构化成使得接纳部分2的内表面与混浊介质1之间留有空间。为了进行检查,该空间填充有光学匹配的介质。该光学匹配的介质被选择成在要成像的混浊介质1与分布在内表面上的源和检测位置之间提供适当的光学耦合。为此目的,向该光学匹配的介质设有与要检查的混浊介质1的光学特性相似的光学特性。
现在,将描述依照该实施例的设备的操作。该设备特别地适于其中检查向其施用对比剂的混浊介质1的方法。对比剂是能够在用适当的光照射时发射第一波长范围的(荧光)光的荧光对比剂。例如,用于光学荧光层析术的典型的荧光染料发射包括大约700-950nm的光谱的波段内的光。向其施用了对比剂的混浊介质1放置在测量体4中。
然后,执行衰减测量,其中来自光源单元6的光随后从所述多个源位置定向到混浊介质。在衰减测量中,当测量衰减测量数据,即入射光的衰减时,在所述多个检测位置对于每个源位置检测用于照射的光的衰减。因此,在衰减测量中,使用混浊介质1的内在对比度。对于来自光源单元6的光的多个不同的波长λ1、...、λk测量用于照射的光的衰减。依照该实施例,这是由光源单元6中的若干不同光源完成的,这些光源发射不同波长λ1、...、λk的光。通常,这些离散波长位于与荧光对比剂的荧光光谱的第一波长范围不同的第二波长范围内。更特别地,在许多情况下,第二波长范围包含更短的波长。优选地,这些光源由发射诸如690nm和730nm之类的不同波长的单色光的激光器形成。
根据这些衰减测量,可以使用本领域中已知的重构方法重构混浊介质1的吸收特性的三维图像。由于对于不同的波长λ1、...、λk执行衰减测量,因而测量对混浊介质1的不同成分敏感。因此,可以以空间分辨的方式重构混浊介质1的不同成分的浓度。这通过(数学上)将混浊介质1建模成多种物质的线性组合来完成,其被证明是一种相当良好的近似。换言之,认为吸收系数μa具有以下结构:ci为成分i的浓度,μai(λ)为成分i的(依赖于波长的)摩尔吸收系数,r为混浊介质内的位置,并且对所有考虑的成分i执行求和。基于来自衰减测量的数据,对于混浊介质1内的每个位置重构所述多种成分的局部组合。这通过利用以下方程来完成:
其中ΦX为波长λX处的测量的强度,Φref为参考测量中测量的强度,i定义混浊介质的成分(例如氧合血红蛋白和脱氧血红蛋白、脂质、水等等),ε为对于这些成分中的每一种的光谱消光,ci为成分i的对应浓度,s定义源位置,d定义检测位置,r定义混浊介质中的相应位置并且GX定义波长λX处的格林函数。
可以依照该方法为其确定局部浓度的成分包括血液、Hb(氧合血红蛋白或脱氧血红蛋白)、水、脂质等等。上面提到的波长特别地非常适合这个方面。因此,根据对于若干离散波长的衰减测量确定混浊介质1在第二波长范围内的吸收特性。
现在,依照本发明,基于重构的关于成分的局部浓度的信息,外插混浊介质1在第一波长范围内的吸收特性。这通过以下方式完成:如上面描述的,若干成分的局部浓度根据衰减测量获悉。此外,例如根据检查这样的成分的单独的测量等等获悉这些成分在第一波长范围内的光谱行为(其与下面描述的荧光测量有关)。因此,基于成分的局部浓度以及这些成分在第一波长范围内的吸收特性的知识,计算混浊介质1在第一波长范围内的吸收特性。换言之,对于第一波长范围的波长λ并且对于混浊介质的每个体元(体像素)计算吸收系数μa(λ,r)。结果,根据对于若干离散波长的衰减测量计算对于第一波长范围的混浊介质1的吸收特性的空间分辨图像。因此,在无需第一波长范围内的附加衰减测量的情况下获得该信息,所述附加衰减测量将需要附加的硬件和附加的测量时间。
此外,依照该实施例,执行荧光测量,其中随后利用来自所述多个源位置的光照射混浊介质1。该荧光测量中使用的光被选择成使得荧光对比剂(由优先地地积聚在混浊介质1中的病变处的荧光染料形成)受该光刺激以发射包含第一波长范围内的波长的荧光。对于所述多个源位置中的每一个,在所述多个检测位置处检测对比剂发射的荧光,作为荧光测量数据。荧光的检测通过在测量体4与检测单元7之间的光路上的适当滤光来执行。
基于荧光测量数据,控制和处理单元8重构混浊介质1中分布的对比剂的荧光发射的体图像。依照该实施例,为了重构荧光发射的体图像,将对于第一波长范围的计算的吸收特性用作重构问题的输入。换言之,通过考虑计算的吸收特性计算对于第一波长范围的格林函数。依照该实施例,通过使用结合第一波长范围处的外插吸收的齐次格林函数的波恩荧光重构来执行荧光图像的重构。可以对于荧光波段的典型荧光波长λf(例如荧光对比剂的发射率具有其最大值的波长)或者对于若干选择的荧光波长合并外插吸收以便将整个荧光光谱考虑在内。例如,以下方程用于计算荧光对比剂的浓度cdye:
其中Φf为荧光的检测的强度,γ(λf)为量子产额(每吸收光子的发射荧光光子的百分比),εsystem为依赖于波长的检测灵敏度(例如在校准测量中确定),并且Gf为相应荧光波长λf处的格林函数,总和在多个荧光波长上求取。
与其中第一波长范围内的吸收的光谱行为不被考虑或者仅被估计为与衰减测量的波长的相似的荧光图像重构相比,重构的荧光图像包括关于荧光对比剂浓度的更高的精确度。因此,提高了关于有关混浊介质内不同位置处荧光对比剂浓度的信息的图像质量。结果,与先前不考虑诸如混浊介质内发色团(血红蛋白、氧合血红蛋白、水、脂质等等)的浓度之类的因素的荧光重构相比,实现了改善的荧光图像。
第二实施例
下面将描述第二实施例。第二实施例基本上与上面描述的第一实施例相对应,并且不同之处仅在于,将荧光对比剂看作所述成分之一,其浓度根据第二波长范围内的衰减测量来确定。因此,依照第二实施例,考虑了荧光对比剂对于荧光的自吸收。这特别在高浓度荧光对比剂的情况下导致改善的荧光图像。
依照第二实施例,荧光对比剂的浓度cf(r)及其吸收系数μaf包含在方程中,即在求和中考虑了项cf(r)μaf(λ,r)。因此,已经根据衰减测量(预先)计算了荧光对比剂的浓度。然后,与第一实施例类似地计算对于第一波长范围的吸收特性,并且附加地考虑荧光对比剂吸收的已知的光谱依赖性。如关于第一实施例所述计算对于相应荧光波长(或多个波长)的格林函数。
然后,像在第一实施例中一样,根据荧光测量计算荧光对比剂的浓度。因此,依照第二实施例,对于第一波长范围的吸收特性的外插包括根据衰减测量确定荧光对比剂浓度的第一估计。根据荧光测量,再次地但是现在更精确地计算荧光对比剂的浓度。
结果,即使对于高浓度的荧光对比剂,也以更高的精确度确定对比剂的浓度。
修改
依照第二实施例的修改,如关于第二实施例所述执行所有步骤。然而,在根据荧光测量重构了荧光对比剂浓度的空间分布之后,将该信息用作再次吸收重构的输入。这意味着根据荧光测量确定的荧光对比剂的浓度在根据衰减测量确定混浊介质的所述多种成分的浓度时用作cf(r),并且重构重复进行。这些步骤然后在具有多次迭代的“反馈环”中重复。对于每次迭代,实现了染料浓度的改进的重构。在理想的情况下,重复上面描述的步骤,直到荧光对比剂的重构浓度收敛到恒定图像。
依照该修改,即使出现高浓度的荧光对比剂,使得荧光对比剂的相当数量的自吸收发生,也可以以提高的精确度根据荧光测量确定荧光对比剂的浓度的空间分布。
尽管上面结合实施例描述了在荧光测量之前执行衰减测量,但是技术人员应当理解,该顺序可以反转。
在上面描述的实施例中,检测单元包括多个不同的检测器,其相应检测位置围绕测量体分布。这允许实现快速而高效的测量。然而,检测单元也可以包括仅仅一个检测器,其相对于测量体的位置在测量期间改变,以便在所述多个不同检测位置处检测光。
此外,上面描述了测量体以具有杯状形状的接纳部分为界。然而,测量体也可以具有不同的形状,例如可以以两个平行板为界,测量期间混浊介质以压缩状态容纳在这两个平行板之间。
Claims (9)
1.用于重构混浊介质内部的荧光图像的方法,该方法包括步骤:
-在测量体(4)中容纳向其施用了荧光对比剂的混浊介质(1),荧光对比剂能够在用光照射时发射第一波长范围内的光;
-通过随后利用来自多个不同源位置的光照射混浊介质(1)并且对于每个源位置在多个检测位置检测从混浊介质(1)发出的光,在多个不同的波长λ1,...,λk下执行衰减测量;
-根据衰减测量对于所述多个不同的波长重构作为混浊介质(1)内部的位置的函数的吸收系数μa(r,λ),以便以空间分辨的方式在所述多个不同的波长λ1,...,λk下根据衰减测量确定混浊介质(1)中多种物质i的局部浓度ci;
-对于第一波长范围的典型荧光波长λf或者对于第一波长范围的若干选择的荧光波长,基于所述多种物质i的局部浓度ci和所述多种物质i在第一波长范围内的吸收系数计算作为混浊介质(1)内部的位置的函数的吸收系数;
-通过随后利用来自所述多个源位置的使得荧光对比剂发射第一波长范围内的光的照射混浊介质(1)并且对于每个源位置在所述多个检测位置检测从荧光对比剂发出的光,执行荧光测量;
-根据荧光测量使用计算的吸收系数重构混浊介质(1)内部荧光对比剂的空间分布的荧光图像。
2.依照权利要求1的方法,其中在重构吸收系数的步骤中,计算吸收系数μa(r,λ)。
3.依照权利要求2的方法,其中通过将混浊介质(1)看作多种物质i的线性组合来计算吸收系数μa(r,λ)。
4.依照权利要求3的方法,其中使用所述多种物质i的已知光谱行为计算对于第一波长范围的波长的吸收系数μa(r,λ)。
5.依照权利要求3或4中任何一项的方法,其中所述多种物质i包括荧光对比剂。
6.依照权利要求5的方法,其中,将重构的荧光图像作为输入反馈到重构作为混浊介质内部的位置的函数的吸收系数μa(r,λ)的步骤,并且迭代地重复计算吸收系数和重构荧光图像的步骤。
7.用于对混浊介质内部成像的设备;该设备包括:
-测量体(4),其适于容纳向其施用了荧光对比剂的混浊介质(1),所述荧光对比剂在照射时能够发射第一波长范围的光;
-光源单元(6),其适于随后利用来自多个不同源位置的光照射测量体(4);光源单元能够有选择地发射多个不同波长λ1,...,λk的光;
-检测单元(7),其适于在多个不同检测位置检测从测量体(4)发出的光;以及
-控制和重构单元(8),其适于控制所述设备以
-通过随后利用来自多个不同源位置的光照射测量体(4)在多个不同的波长下执行衰减测量,并且对于每个源位置在多个检测位置处检测从测量体(4)发出的光;
-通过随后利用来自所述多个源位置的使得荧光对比剂发射第一波长范围内的光的照射测量体(4)并且对于每个源位置在所述多个检测位置检测从荧光对比剂发出的光,执行荧光测量;
-根据衰减测量对于所述多个不同的波长重构作为测量体(4)内的位置的函数的吸收系数μa(r,λ),以便以空间分辨的方式在所述多个不同的波长λ1,...,λk下根据衰减测量确定混浊介质(1)中多种物质i的局部浓度ci;
-对于第一波长范围的典型荧光波长λf或者对于第一波长范围的若干选择的荧光波长,基于所述多种物质i的局部浓度ci和所述多种物质i在第一波长范围内的吸收系数计算作为测量体(4)内的位置的函数的吸收系数;
-根据荧光测量使用计算的吸收系数重构测量体(4)内荧光对比剂的空间分布的荧光图像。
8.依照权利要求7的设备,其中在该设备中提供可在测量体(4)与检测单元(7)之间的光路上引入的滤光器以用于荧光测量。
9.依照权利要求7或8中任何一项的设备,其中该设备是医疗图像采集设备。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP08153363.0 | 2008-03-27 | ||
EP08153363 | 2008-03-27 | ||
PCT/IB2009/051179 WO2009118672A1 (en) | 2008-03-27 | 2009-03-20 | Method for reconstructing a fluorescent image of the interior of a turbid medium and device for imaging the interior of a turbid medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101980656A CN101980656A (zh) | 2011-02-23 |
CN101980656B true CN101980656B (zh) | 2013-12-04 |
Family
ID=40720002
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009801110542A Expired - Fee Related CN101980656B (zh) | 2008-03-27 | 2009-03-20 | 用于重构混浊介质内部的荧光图像的方法以及用于对混浊介质内部成像的设备 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8520921B2 (zh) |
EP (1) | EP2259717A1 (zh) |
CN (1) | CN101980656B (zh) |
WO (1) | WO2009118672A1 (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103635131B (zh) * | 2011-06-28 | 2018-07-27 | 皇家飞利浦有限公司 | 对相关组织样本进行光学分析的设备 |
US9442005B2 (en) | 2014-07-30 | 2016-09-13 | Corning Optical Communications LLC | Non-contact methods of measuring insertion loss in optical fiber connectors |
US9730649B1 (en) | 2016-09-13 | 2017-08-15 | Open Water Internet Inc. | Optical imaging of diffuse medium |
US10506181B2 (en) | 2018-03-31 | 2019-12-10 | Open Water Internet Inc. | Device for optical imaging |
US10778912B2 (en) | 2018-03-31 | 2020-09-15 | Open Water Internet Inc. | System and device for optical transformation |
US10778911B2 (en) | 2018-03-31 | 2020-09-15 | Open Water Internet Inc. | Optical transformation device for imaging |
US10966612B2 (en) | 2018-06-14 | 2021-04-06 | Open Water Internet Inc. | Expanding beam optical element |
US10962929B2 (en) | 2018-09-14 | 2021-03-30 | Open Water Internet Inc. | Interference optics for optical imaging device |
US10874370B2 (en) | 2019-01-28 | 2020-12-29 | Open Water Internet Inc. | Pulse measurement in optical imaging |
US10955406B2 (en) | 2019-02-05 | 2021-03-23 | Open Water Internet Inc. | Diffuse optical imaging with multiple beams |
US11320370B2 (en) | 2019-06-26 | 2022-05-03 | Open Water Internet Inc. | Apparatus for directing optical and acoustic signals |
US11581696B2 (en) | 2019-08-14 | 2023-02-14 | Open Water Internet Inc. | Multi-channel laser |
US11622686B2 (en) | 2019-11-22 | 2023-04-11 | Open Water Internet, Inc. | Optical imaging with unshifted reference beam |
US11819318B2 (en) | 2020-04-27 | 2023-11-21 | Open Water Internet Inc. | Optical imaging from light coherence |
US11259706B2 (en) | 2020-05-19 | 2022-03-01 | Open Water Internet Inc. | Dual wavelength imaging and out of sample optical imaging |
US11559208B2 (en) | 2020-05-19 | 2023-01-24 | Open Water Internet Inc. | Imaging with scattering layer |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1200174A (zh) * | 1995-08-24 | 1998-11-25 | 普渡研究基金会 | 基于荧光寿命的人体组织及其它无规则介质成象技术和光谱技术 |
CN1429089A (zh) * | 2000-05-09 | 2003-07-09 | 成象诊断系统公司 | 利用多波长同时数据采集进行乳房成像的医用光学成像扫描仪 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7556602B2 (en) * | 2000-11-24 | 2009-07-07 | U-Systems, Inc. | Breast cancer screening with adjunctive ultrasound mammography |
US8110281B2 (en) * | 2004-07-02 | 2012-02-07 | 3Dtl, Inc. | Systems and methods for creating optical effects on media |
WO2007054846A2 (en) * | 2005-11-10 | 2007-05-18 | Philips Intellectual Property & Standards Gmbh | Absorption and scattering map reconstruction for optical fluorescence tomography |
US8239006B2 (en) * | 2006-07-06 | 2012-08-07 | The University Of Connecticut | Method and apparatus for medical imaging using near-infrared optical tomography and fluorescence tomography combined with ultrasound |
US8070682B2 (en) * | 2006-07-19 | 2011-12-06 | The University Of Connecticut | Method and apparatus for medical imaging using combined near-infrared optical tomography, fluorescent tomography and ultrasound |
JP2008115223A (ja) * | 2006-11-01 | 2008-05-22 | Nec Lighting Ltd | 蛍光体含有ガラスシート、その製造方法及び発光装置 |
JP2009018048A (ja) * | 2007-07-12 | 2009-01-29 | Fujifilm Corp | 医用画像表示装置、方法及びプログラム |
WO2009089543A2 (en) * | 2008-01-10 | 2009-07-16 | The Ohio State University Research Foundation | Fluorescence detection system |
US8112146B2 (en) * | 2008-07-24 | 2012-02-07 | Massachusetts Institute Of Technology | Three-dimensional imaging using a luminescent surface and a differentially attenuating medium |
US20100187441A1 (en) * | 2009-01-28 | 2010-07-29 | Bio-Rad Laboratories, Inc. | Transilluminator adaptor for conversion of ultraviolet radiation to visible light |
-
2009
- 2009-03-20 US US12/933,887 patent/US8520921B2/en not_active Expired - Fee Related
- 2009-03-20 WO PCT/IB2009/051179 patent/WO2009118672A1/en active Application Filing
- 2009-03-20 EP EP09725479A patent/EP2259717A1/en not_active Ceased
- 2009-03-20 CN CN2009801110542A patent/CN101980656B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1200174A (zh) * | 1995-08-24 | 1998-11-25 | 普渡研究基金会 | 基于荧光寿命的人体组织及其它无规则介质成象技术和光谱技术 |
CN1429089A (zh) * | 2000-05-09 | 2003-07-09 | 成象诊断系统公司 | 利用多波长同时数据采集进行乳房成像的医用光学成像扫描仪 |
Also Published As
Publication number | Publication date |
---|---|
CN101980656A (zh) | 2011-02-23 |
US20110026851A1 (en) | 2011-02-03 |
US8520921B2 (en) | 2013-08-27 |
EP2259717A1 (en) | 2010-12-15 |
WO2009118672A1 (en) | 2009-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101980656B (zh) | 用于重构混浊介质内部的荧光图像的方法以及用于对混浊介质内部成像的设备 | |
Darne et al. | Small animal fluorescence and bioluminescence tomography: a review of approaches, algorithms and technology update | |
Graves et al. | A submillimeter resolution fluorescence molecular imaging system for small animal imaging | |
CN102058393B (zh) | 基于反射光谱测量的皮肤生理参数与光学特性参数的测量方法 | |
CN101902952B (zh) | 用于检测混浊介质内部中不均匀性的存在的方法和用于对混浊介质内部成像的设备 | |
US5840035A (en) | Method for the spectroscopic examination of a biological tissue | |
WO2007111669A2 (en) | Combined x-ray and optical tomographic imaging system | |
JP6573732B2 (ja) | 近赤外分光法及びフォトニックミキサ装置での断層撮影法を用いる腫瘍低酸素症の継続的モニタリング | |
Ueda et al. | Time-resolved optical mammography and its preliminary clinical results | |
EP1954177B1 (en) | A device for imaging an interior of a turbid medium | |
CN101262812B (zh) | 光学层析成像装置 | |
Mastanduno et al. | Automatic and robust calibration of optical detector arrays for biomedical diffuse optical spectroscopy | |
Yu et al. | Phantom verification for a ring-scanning and prone diffuse optical imaging system | |
Alhemsi et al. | Time-resolved near-infrared spectroscopic imaging systems | |
US20100249606A1 (en) | Method and device for reconstructing a fluorescence image of the interior of a turbid medium | |
WO2009060383A2 (en) | Method for calibrating a device for imaging the interior of turbid media and method for correcting measured intensities in such a device | |
Michaelsen et al. | Using Partial Frequency Domain Data to Improve Reconstruction Accuracy in Continuous Wave Breast Tomosynthesis Guided Diffuse Optical Spectroscopy | |
Ueda et al. | Development of optical mammography based on analysis of time-resolved photon path distribution | |
Wang et al. | Flying spot remote sensing of ICG kinetics of undeformed tissues | |
Pöschinger | Non-Contact Optical Fluorescence Tomography for Small Animal Imaging: System Development and Multispectral Applications | |
WO2009077953A1 (en) | Method for reconstructing an image of an interior of a turbid medium and device for imaging the interior of turbid media | |
Brooksby et al. | Development of a system for simultaneous MRI and Nearinfrared diffuse tomography to diagnose breast cancer | |
WO2010103440A1 (en) | Irradiating an absorption contrast agent contained in a turbid medium | |
Jiang et al. | Multispectral NIR diffuse optical tomography system development | |
WO2009122342A1 (en) | Method and device for visualizing accumulation or depletion of a contrast agent in a turbid medium |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131204 Termination date: 20180320 |