CN105866035B - 一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法与技术 - Google Patents

一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法与技术 Download PDF

Info

Publication number
CN105866035B
CN105866035B CN201610197338.XA CN201610197338A CN105866035B CN 105866035 B CN105866035 B CN 105866035B CN 201610197338 A CN201610197338 A CN 201610197338A CN 105866035 B CN105866035 B CN 105866035B
Authority
CN
China
Prior art keywords
spatial frequency
modulation
tissue
large area
microstructure
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.)
Active
Application number
CN201610197338.XA
Other languages
English (en)
Other versions
CN105866035A (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.)
Wenzhou Medical University
Original Assignee
Wenzhou Medical University
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 Wenzhou Medical University filed Critical Wenzhou Medical University
Priority to CN201610197338.XA priority Critical patent/CN105866035B/zh
Publication of CN105866035A publication Critical patent/CN105866035A/zh
Application granted granted Critical
Publication of CN105866035B publication Critical patent/CN105866035B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • G01N2021/1765—Method using an image detector and processing of image signal
    • G01N2021/177—Detector of the video camera type

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明提供一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法与技术，其包括以下步骤：一、控制光源输出多个不同频率的空间调制光条纹图案到组织样品；二、通过CCD多次采集经组织样品散射后的反射光强；三、对采集到的反射光强数据进行分析与处理；四、采用标准三相移法或SSMD解调法等解调出不同频率的交流、直流分量和调制传递函数；五、获得散射结构系数SSI，通过高空间频率的SFDI调制不仅可以获得吸收系数和散射系数，而且可以获得散射光的相位函数和散射结构系数(SSI)。高空间频率SFDI在大面积组织学诊断上具有重要应用，它可以定量获取大面积的组织光学特性及散射特征的分布图，用于生物组织的客观诊断。

Description

一种基于空间频域调制大面积解析微观结构的快速无损组织 活检方法与技术
技术领域
本发明涉及一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法与技术。
背景技术
通过组织的光散射可以用于诊断组织，并已大量应用于生物医学。在许多不同成像方式中，空间频域成像SFDI目前作为一种微观结构成像方式在广泛的领域引起了极大重视。SFDI可以大面积地用于量化浑浊介质的光子性能。尤其是高空间频域成像HSFDI可以提高成像的分辨率并定量地提取光散射的相位函数，HSFDI的问题是缺少一种能够在高空间频率分析上描述光反射系数的合适模型。另外，传统的 SFDI结构不能实现光学成像***所要求的准确性和高分辨率。
发明内容
为了解决上述问题，本发明提供了一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法与技术。
本发明提供：一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法，其包括以下步骤：
控制光源输出多个不同频率的空间调制光条纹图案到组织样品；
通过CCD多次采集经组织样品散射后的反射光强；
对采集到的反射光强数据进行分析与处理；
采用标准三相移法或SSMD解调法解调出不同频率的交流、直流分量和调制传递函数；
通过获得散射结构系数SSI。
所述散射结构系数SSI通过以下步骤获得，首先根据
其中μb≡μt/p(π)以及 lΘ≡Θ/μs，散射介质的随机相位函数pSAA(θ)为获得该相位函数pSAA(θ)及散射结构系数SSI。
一种基于上述快速无损组织活检方法的高空间频域成像技术，其包括光源发射装置、透镜、分束镜以及采集装置，所述光源发射装置发射的光经过透镜照射在分束镜上，并通过分束镜反射在组织样品上，组织样品的反向散射光经过图像本身成像于采集装置处。
所述光源发射装置为数字微镜设备。
所述采集装置为相机。
本发明通过高空间频率调制实现的HSFDI不仅可以获得吸收系数和散射系数，而且可以获得散射光的相位函数和散射结构系数(SSI)。 HSFDI在大面积组织学诊断中具有重要应用，它可以定量获取大面积组织的光学特性及散射特征分布图用于生物组织的客观诊断。
附图说明
图1为本发明的原理示意图。
图2是在SFDI调制频率低(f＝0.244mm-1)，中(f＝0.569mm-1)，高(f＝1.365mm-1)频率下2％浓度的脂肪乳下计算所得的光学性质，底部显示了由低频(红色)，中频(绿色)和高频(蓝色)频率的切面图。
图3通过SFDI调制在低频(f＝0.682mm-1中间列)和高频((f＝3.412 mm-1右列)所获得的猪肉组织的光学特征图，SSI特征图也显示在图上。
图4为肾癌组织的μa，μ′s(f＝0.682mm-1)分布图和SSI(f＝1.706 to 3.412mm-1)分布图，图中虚线框为癌化显著区域。
图5为胃癌组织的μa,μ′s(f＝0.682mm-1)分布图和SSI(f＝1.706 to 3.412mm-1)分布图，图中虚线框为癌化显著区域。
图6为健康男性手掌的μa，μ′s(f＝0.682mm-1)分布图和SSI(f＝1.706 to 3.412mm-1)分布图。
具体实施方式
下面结合附图对本发明实施例作进一步说明：
光在浑浊介质(如生物组织)中的散射是一个复杂的过程，光的传播可以用玻尔兹曼辐射传播方程(RTE)精确地描述，但目前其计算耗时而难以应用。另外RTE的扩散近似方程只适用于光源距离探测器的距离是传输步长的3/4倍以上时才有效，而在光源与探测器近距离时是无效的。然而，散射介质的相位函数包含了介质形态的基本信息，当光源与探测器的距离小于单个传输步长时非扩散光对反射率影响显著。因此，HSFDI最重要的一步是将非扩散光的反射率与相位函数的确切形式联系在一起。
而徐敏在申请日2015-10-15，申请号为2015106669451申请的发明一种大面积获取组织光学参数及微观结构的光反射成像技术，则提出的相位函数分析模型包含在光源与探测器近距离时非扩散反射光辐射传输的小角度散射近似，这个模型还表达了在近距离时散射介质的相位函数与反射光的关系。
在前向光散射介质中的后向散射光包含3种类型：第一种是经历了多次小角度散射和一次完全地大角度散射的SAA光子；第二种是经历了两次完全地大角度散射的“蛇”形光；第三种是经历了两次以上大角度散射的扩散光。对于任意光源-探测器距离和在整个空间频域中，根据它们各自的范围反射光公式可以分成低空间频率和高空间频率两部分：
其中qc～2πμr是两范围的交点。低空间频率表达式中的μb/2μt是由 SAA光子中的未散射成分引起的，公式(1)可以准确描述具有低吸收到中等吸收(μa/μ′s＜1)前向光散射介质的反射光。尤其，当qltr＞＞1 时表示SAA光子反射的主要表达式ISAA(q)可以简化为：
其中后向散射系数μb≡μt/p(π)，传播距离lΘ≡Θ/μs。散射介质的随机相位函数pSAA(θ)可以如下表示：
其中
模拟体系真实的相位函数p(θ)，即pb值由后向散射到后半球的光子数目决定，散射角的均方根Θ与组织和细胞中的大结构体相关，尤其是各项异性因子g：
g＝(1-2pb)(1-Θ2/2) (6)
当μa/μs＜＜1,μt＝μa+μs≈μs和μb/μt＝pb/4π时，公式(2)可以简化为：
因此，在可测量的频率范围内用各种高空间频率测量多重反射光 ISAA(q)可以拟合出公式(7)得到pb和lΘ。
相位函数的一个新特性定义为散射结构系数(Scattering Structural Index，SSI)，该系数可以用来评估浑浊介质的形态学特征。由于在大多数生物组织中Θ＞＞pb，所以将SSI定义为：
而本发明则提供了一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法，其包括以下步骤：
控制光源输出多个不同频率的空间调制光条纹图案到组织样品；
通过CCD多次采集经组织样品散射后的反射光强；
对采集到的反射光强数据进行分析与处理；
采用标准三相移法或SSMD解调法调出不同频率的交流、直流分量和调制传递函数；
通过获得散射结构系数SSI。
空间频域成像(SFDI)是一种大范围，非接触的方法，其使用空间调制构成光图案以一维形式，不同的频率投影到组织上：
其中的S0,M0,fx和α分别为光源密度，调制深度，空间频率和空间相位，反射形式的强度如下式：
IOUT＝IDC+MAC(x)cos(2πfxx+α) (10)
浑浊介质中的反射光子密度波MAC(x)的调制转移依赖于组织的光学特性，并能通过光的传播方程进行模拟。为了能解析出MAC(x)需要利用单次解调方法，通常采用的方法是三相位解调法而另外一种的单频快照解调技术(SSMD)已经在实时成像平台上实现了。
为了获得非扩散光子对光反射的敏感性，当扩散模型是无效时， MAC(x)可通过多个高空间频率(fx＞μ′s/3)的测量获得。随后，我们利用 SAA分析模型去定量的估计散射结构值(SSI)的分布图。不同于使用低频成像的传统SFDI方案，HSFDI方案需要更高精度、更高分辨率的同步调制照明***和光学成像***。
而本发明还公开了一种基于上述快速无损组织活检方法的高空间频域成像技术，其包括光源发射装置、透镜lens、分束镜beam splitter 以及采集装置，所述光源发射装置发射的光经过透镜照射在分束镜上，并通过分束镜反射在组织样品上，组织样品的反向散射光经过图像本身成像于采集装置处。出于对像差优化的考虑，为制照明***的中心光学设备使用了薄膜分束器，来消除鬼影，DMD的图案被放大后投影到样本上，样本的反向散射光经过图像本身成像于相机上。
所述光源发射装置为数字微镜设备，其型号为数字微镜设备LightCrafterTM 4500，Texas Instruments)。采用红色通道(623nm)，DMD 设置为图案模式(pattern mode)，每幅画的跳变时间为10000μs，
所述采集装置为相机，其型号为Canon 5D Mark III相机，用于收集反射图像，曝光时间为1/10S。
DMD输出的光强并不是线性的，通过漫反射利用朗伯反射标准来纠正这种非线性。照在样本表面的是13mmx13mm照射范围，HSFDI ***的解析度为7.7μmx7.7μm，最大的调制空间频率为3.412mm-1。
通过均匀体模实验。我们发现边界效应会影响HSFDI的表现，实际上是切去边界后的空间调制图案，其长度是固定的。我们得出由于边界效应使图像质量退化，这将导致计算介质的调制传递函数不准确，从理论上说，当图像中的调制模式中的周期数减少时，伪影会随着时间的增加而增加。因此，尤其我们SFDI设备在不特别大的13mm*13mm 的范围内时，减少伪影的方式是最大限度的使用高空间频率的调制光图案，减少工件。
分别照射不同频率(0.244,0.569和1.365mm-1)的调制正弦图案于2％的脂肪乳悬乳液(10％的脂肪乳通过1:4进行稀释)上，如图2 所示。基于扩散传播方程计算的μa和μ′s的平均值在表1所示。当排除边界效应对HSFDI的影响下，计算所得的μa和μ′s值与理论值随着频率的增加而越来越接近，图2中入射光的波长为0.623μm，像素格点为 7.7μm×7.7μm。
表1，为2％浓度的脂肪乳悬乳液在不同频率下计算所得的μa和μ′s。
基于扩散光的传播方法的工作条件下，对在低，高调制频率对10mm*10mm的猪肉组织进行SFDI成像实验观察。得到如图3所描述的μa和μ′s光学分布图。在中间一列显示的结果为低频(f＝0.682mm-1)，在右边一列现实的结果为高频(f＝3.412mm-1)。可以发现，在高调制频率的空间频域成像中利用基于扩散近似传输理论的查表法重建出的光学参数出现溢出现象，这跟我们所预期的一样，即在探测器与光源点近距离时RTE的扩散近似传输理论将失效。
然后，我们通过拟合的反射光的小角度散射近似(SAA)模型(看公式7)，从三个高空间频率提取猪肉组织散射结构指数(SSI)分布图。此结果如图3所示。从HSFDI得到的非-漫散射反射光携带着相位函数的关键特征，与约化散射系数相比，SSI具有更好的结构对比度 (μ′s)。通过SFDI技术可以发现新鲜的肾癌和胃癌组织，其结果分别对应图4和图5。图中虚线框区域为医生诊断出的癌变范围。我们可以看到SSI可以对诊断区域关联的结构变换的定量表征。高的SSI值通常表示生物组织的制粒的细胞核，可被用于癌化的前期诊断。
最后，我们论证了HSDFI方法确定超过10mm*10mm男性手掌的 SSI分布图。同时，浅表周围的组织的相位函数相关的SSI有更好的结构对比度并显示出清晰的掌纹。
实施例不应视为对本发明的限制，但任何基于本发明的精神所作的改进，都应在本发明的保护范围之内。

Claims (2)

1.一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法，其特征在于：其包括以下步骤：
一、控制光源输出多个不同频率的空间调制光条纹图案到组织样品；
二、通过CCD多次采集经组织样品散射后的反射光强；
三、对采集到的反射光强数据进行分析与处理；
四、采用标准三相移法或SSMD解调法解调出不同频率的交流、直流分量和调制传递函数；
五、通过获得散射结构系数SSI。
2.根据权利要求1所述的基于空间频域调制大面积解析微观结构的快速无损组织活检方法，其特征在于：所述散射结构系数SSI通过以下步骤获得，首先根据 其中以及散射介质的随机相位函数pSAA(θ)为获得该相位函数pSAA(θ)及散射结构系数SSI。
CN201610197338.XA 2016-03-30 2016-03-30 一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法与技术 Active CN105866035B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610197338.XA CN105866035B (zh) 2016-03-30 2016-03-30 一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法与技术

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610197338.XA CN105866035B (zh) 2016-03-30 2016-03-30 一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法与技术

Publications (2)

Publication Number Publication Date
CN105866035A CN105866035A (zh) 2016-08-17
CN105866035B true CN105866035B (zh) 2019-04-05

Family

ID=56626845

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610197338.XA Active CN105866035B (zh) 2016-03-30 2016-03-30 一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法与技术

Country Status (1)

Country Link
CN (1) CN105866035B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11867625B2 (en) 2019-02-03 2024-01-09 Bar Ilan University System and method for imaging via scattering medium
CN111265187B (zh) * 2020-01-21 2023-04-07 温州医科大学 一种干燥综合症无创病理诊断装置及方法
CN114018819B (zh) * 2021-03-22 2022-11-11 北京航空航天大学 基于空间频域成像的光学特性测量方法、装置及***
CN116990237B (zh) * 2023-09-26 2023-12-12 国网江苏省电力有限公司电力科学研究院 一种适用于窄带宽的增强式大景深太赫兹收发装置

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102883658A (zh) * 2009-11-19 2013-01-16 调节成像公司 用于使用结构化照明经由单元件检测来分析浑浊介质的方法和设备

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008530536A (ja) * 2005-02-09 2008-08-07 インライト ソリューションズ，インク． 検体の非侵襲性を測定する方法および装置
US9968261B2 (en) * 2013-01-28 2018-05-15 The General Hospital Corporation Apparatus and method for providing diffuse spectroscopy co-registered with optical frequency domain imaging
US9927350B2 (en) * 2013-10-17 2018-03-27 Trustees Of Boston University Thermal property microscopy with frequency domain thermoreflectance and uses thereof
CN105245761B (zh) * 2015-09-30 2018-08-24 温州医科大学 单次快照多频解调方法
CN105300928B (zh) * 2015-10-15 2018-08-24 温州医科大学 一种大面积获取组织光学参数及微观结构的光反射成像技术

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102883658A (zh) * 2009-11-19 2013-01-16 调节成像公司 用于使用结构化照明经由单元件检测来分析浑浊介质的方法和设备

Also Published As

Publication number Publication date
CN105866035A (zh) 2016-08-17

Similar Documents

Publication Publication Date Title
US20230363647A1 (en) Structured-light imaging systems and methods for determining sub-diffuse scattering parameters
CN105866035B (zh) 一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法与技术
Rehbinder et al. Ex vivo Mueller polarimetric imaging of the uterine cervix: a first statistical evaluation
Nandy et al. Characterizing optical properties and spatial heterogeneity of human ovarian tissue using spatial frequency domain imaging
Lue et al. Portable optical fiber probe-based spectroscopic scanner for rapid cancer diagnosis: a new tool for intraoperative margin assessment
Foschum et al. Fully automated spatially resolved reflectance spectrometer for the determination of the absorption and scattering in turbid media
JP6192297B2 (ja) 被検体情報取得装置、表示制御方法、およびプログラム
JP2018529929A (ja) 生検標本蛍光イメージング装置及び方法
US20090005676A1 (en) Method and Device For Reconstructing a Three-Dimensional Fluorescence Optical Tomography Image by Double Measurement
Lapointe et al. A multi-view time-domain non-contact diffuse optical tomography scanner with dual wavelength detection for intrinsic and fluorescence small animal imaging
US10485425B2 (en) Apparatus and methods for structured light scatteroscopy
JP5984547B2 (ja) 被検体情報取得装置およびその制御方法
Gallwas et al. Detection of cervical intraepithelial neoplasia by using optical coherence tomography in combination with microscopy
US7692160B2 (en) Method and system of optical imaging for target detection in a scattering medium
Gordon et al. Quantitative phase and polarization imaging through an optical fiber applied to detection of early esophageal tumorigenesis
JP5847454B2 (ja) 被検体情報取得装置、表示制御方法およびプログラム
McClatchy III et al. Light scattering measured with spatial frequency domain imaging can predict stromal versus epithelial proportions in surgically resected breast tissue
Stier et al. Imaging sub-diffuse optical properties of cancerous and normal skin tissue using machine learning-aided spatial frequency domain imaging
CN112535465B (zh) 一种基于片层光的三维血流速度成像方法及装置
WO2017166082A1 (zh) 一种基于空间频域调制大面积解析微观结构的快速无损组织活检方法与技术
JP6072206B2 (ja) 被検体情報取得装置および表示方法
CN112449683A (zh) 用于确定表面下荧光物体的深度和浓度的装置和方法
Petusseau et al. Subsurface fluorescence time-of-flight imaging using a large-format single-photon avalanche diode sensor for tumor depth assessment
JP2016185429A (ja) 情報処理装置および情報処理方法
Wang et al. A hybrid reconstruction algorithm for fluorescence tomography using Kirchhoff approximation and finite element method

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