CN112998683A - 一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统与方法 - Google Patents
一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统与方法 Download PDFInfo
- Publication number
- CN112998683A CN112998683A CN202110201906.XA CN202110201906A CN112998683A CN 112998683 A CN112998683 A CN 112998683A CN 202110201906 A CN202110201906 A CN 202110201906A CN 112998683 A CN112998683 A CN 112998683A
- Authority
- CN
- China
- Prior art keywords
- module
- light source
- light
- representing
- image acquisition
- 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
- 206010025282 Lymphoedema Diseases 0.000 title claims abstract description 24
- 208000002502 lymphedema Diseases 0.000 title claims abstract description 24
- 238000005516 engineering process Methods 0.000 title claims abstract description 21
- 206010006187 Breast cancer Diseases 0.000 title claims abstract description 19
- 208000026310 Breast neoplasm Diseases 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000012634 optical imaging Methods 0.000 title claims abstract description 19
- 210000001364 upper extremity Anatomy 0.000 title claims abstract description 18
- 230000003287 optical effect Effects 0.000 claims abstract description 21
- 238000001514 detection method Methods 0.000 claims abstract description 16
- 230000008081 blood perfusion Effects 0.000 claims abstract description 11
- 238000005286 illumination Methods 0.000 claims abstract description 9
- 230000000694 effects Effects 0.000 claims abstract description 8
- 238000010521 absorption reaction Methods 0.000 claims description 27
- 238000003384 imaging method Methods 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 102000001554 Hemoglobins Human genes 0.000 claims description 9
- 108010054147 Hemoglobins Proteins 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000003595 spectral effect Effects 0.000 claims description 9
- 230000017531 blood circulation Effects 0.000 claims description 8
- 230000004089 microcirculation Effects 0.000 claims description 7
- 230000008338 local blood flow Effects 0.000 claims description 6
- 238000009792 diffusion process Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 3
- 235000013405 beer Nutrition 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 229910052760 oxygen Inorganic materials 0.000 claims description 3
- 239000001301 oxygen Substances 0.000 claims description 3
- 230000002980 postoperative effect Effects 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 abstract description 2
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 210000001519 tissue Anatomy 0.000 description 19
- 238000005259 measurement Methods 0.000 description 9
- 206010030113 Oedema Diseases 0.000 description 8
- 239000000126 substance Substances 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 208000004085 Breast Cancer Lymphedema Diseases 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000001926 lymphatic effect Effects 0.000 description 3
- 206010033675 panniculitis Diseases 0.000 description 3
- 230000010412 perfusion Effects 0.000 description 3
- 102000004169 proteins and genes Human genes 0.000 description 3
- 108090000623 proteins and genes Proteins 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 210000004304 subcutaneous tissue Anatomy 0.000 description 3
- 238000001356 surgical procedure Methods 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 208000035475 disorder Diseases 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000010992 reflux Methods 0.000 description 2
- 239000006096 absorbing agent Substances 0.000 description 1
- 238000000862 absorption spectrum Methods 0.000 description 1
- 210000002048 axillary vein Anatomy 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002512 chemotherapy Methods 0.000 description 1
- 230000004087 circulation Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 210000003722 extracellular fluid Anatomy 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 210000004880 lymph fluid Anatomy 0.000 description 1
- 210000001165 lymph node Anatomy 0.000 description 1
- 229920002521 macromolecule Polymers 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000008961 swelling Effects 0.000 description 1
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/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4869—Determining body composition
- A61B5/4875—Hydration status, fluid retention of the body
- A61B5/4878—Evaluating oedema
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Physiology (AREA)
- Hematology (AREA)
- Cardiology (AREA)
- Radiology & Medical Imaging (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供了一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统与方法,包括光源模块,图像采集模块和上位机;所述图像采集模块连接所述光源模块和所述上位机;所述光源模块包括激光、白光和两个单波长LED光源;利用激光照明,使用图像采集模块拍摄图像,通过上位机获取图像R通道的数据,并基于激光散斑效应来计算血流灌注情况。其中本发明的有益效果是:在乳腺癌术后上肢淋巴水肿检测中,避免了与患者检测部位的直接接触,并且光学检测原理受外界环境影响小,所以具有灵敏度高、可靠性高、抗干扰能力强等特点。
Description
技术领域
本发明涉及激光测量和光学吸收技术及医疗领域,特别涉及一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统与方法。
背景技术
乳腺癌相关淋巴水肿(Breast cancer-related lymphedema,BCRL)是乳腺癌根治术后最常见的并发症。据调查水肿发生率高达16.6%。淋巴水肿一旦发生,将会长期困扰患者,很难在短期内治愈,且对患者的危害随着肿胀进行逐渐加重,因此前期有效的检测和评估,将会降低术后上肢淋巴水肿地发生率,减轻水肿程度,促进患者康复。
已知造成水肿的原因主要是由于手术对腋下淋巴结的切除或者手术前后的化疗,造成患者淋巴循环障碍,使得淋巴液中的蛋白大分子因回流不畅聚集在细胞外液。也有研究认为乳腺癌腋窝淋巴清扫改变了腋静脉的血流动力学,使之容易发生淋巴回流障碍,导致术后淋巴水肿。目前临床上对淋巴水肿的评估和测量所采用的技术及效果也各有不同,由于水肿一般分为三个阶段,目前还没有一种评估方法能完全覆盖这三个阶段。
激光散斑血流成像技术也称激光散斑衬比成像(Laser speckle contrastimaging,LSCI)技术,该技术利用血管中红细胞运动产生的后向动态散斑对比度值来获取血流的速度信息,通过成像方式能够获得全场的二维高分辨率血流分布图像,无需机械扫描或注入外源物质就可实现长时间的血流监测与测量。再结合CCD相机以及高性能的并行运算设备,便可达到毫秒量级的时间分辨率和微米量级的空间分辨率。
比尔兰伯特定律的相关实验表明,人体皮肤的吸收系数是和光的波长以及皮下组织的含量相关的,可见光(400-780nm)在皮肤内经历强烈的吸收和散射。光线投射到人体时,会发生如下三种反应:反射,散射和吸收,生物体内的主要光吸收成分是水,蛋白质和脂肪等,水总共占生物体重量的约70%,水对近红外光具有很强的吸收性,在生物体中除了水以外的典型的光吸收体有血红蛋白,血红蛋白有被氧化的状态与未被氧化的状态,这两种状态的吸收光谱也是不相同的。人体对光的吸收强度直接取决于血红蛋白,水和脂肪等物质占据组织的体积分数,人体皮肤表层内的不同物质对波长不同的光的吸收系数也不同。
激光作为诊断和治疗的一种手段在生物医学中的应用研究层出不穷,本发明中利用了多模态光学成像的思想对患者局部血流灌注进行了研究,利用了人体内不同物质与光的相互作用探究了淋巴水肿患者体内的水含量,通过与患者的皮肤电阻抗和臂维等传统评估技术获得的数据的对比,也验证了本发明的有效性,这在指导临床上对水肿有效而直接的评估有着重要的意义。
但目前存在光源单一、具有一定局限性的问题。
发明内容
为了解决上述技术问题,本发明中披露了一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统与方法,本发明的技术方案是这样实施的:
一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统,包括光源模块,图像采集模块和上位机;所述图像采集模块连接所述光源模块和所述上位机;所述光源模块包括激光光源、白光光源和两个单波长LED光源;
利用激光光源照明,使用图像采集模块拍摄图像,通过上位机获取图像RGB三通道的数据,基于激光散斑效应来计算血流灌注情况,
其中IA代表激光光源照明下经过吸收和散射的反射光的光强度,K表示衬比度值,其平方与血流灌注成反比;
利用光源模块中的两个单波长LED光源单独照明,使用图像采集模块拍摄图像,通过上位机获取图像RGB三通道的数据,结合衬比度值数据,基于下述最优化目标计算白光校准矩阵S:
其中代表两个单波长LED光源单独照明时图像采集模块记录的反射光的光强度,代表入射光(单波长LED光源单独照明)强度,F代表图像RGB三通道的滤光片光谱特性,G为两个单波长LED光源分别对应的光谱矩阵,a和b为数值参数,为l1秩的平方,为l2秩的平方,算子定义如下:
将两个单波长LED光源单独照明并分别计算得到的白光校准矩阵S1和S2用于修正白光照明的图像数据:
其中,I表示记录到的白光反射光强数据,Iwhite表示修正后的白光数据;
利用光源模块中的白光光源照明,使用图像采集模块拍摄图像,通过上位机获取图像RGB三通道的数据,结合修正矩阵获得修正后的白光数据Iwhite,使用比尔兰伯特定律实现组织含量的计算,比尔兰伯特定律计算组织含量的基本定义如下:
其中Iwhite表示修正后的白光反射图像数据,Lw0表示入射白光光强,L表示光程,μλ表示吸收与散射的结合系数;
所述上位机对采集到的图像进行处理,其中处理方式包括数字图像处理和信号处理进行局部血流微循环的测量,对于经过吸收后返回的光学图像,处理方法是利用人体内组织的光谱吸收特性结合比尔兰伯特定律,进行组织含量的测量,在的图片中提取R,G,B三个通道的像素数据,利用如下的公式:
其中,IR,IG,IB,分别表示反射光在三个通道的光强度数值,Iwhite表示白光的发射光强度, 1, 1,θ1分别表示图像RGB三通道的光电响应系数,分别表示组织中的水,不含氧血红蛋白,含氧血红蛋白和脂肪的含量,伪,纬,未,尾分别表示上述四种组织在相应波长下的吸收系数,L1表示白光的光程。
优选地,所述光源模块还包括散光装置。
优选地,所述激光模块采用波长为830nm的激光器。
优选地,所述白光模块的波段为400-780nm。
优选地,所述两个单波长LED光源其中一个波长为450nm,另一个波长为590nm。
优选地,所述图像采集模块包括图像传感器和成像镜头。
优选地,所述图像采集模块将采集到的散斑和经过吸收后返回的光学图像传送至上位机,且将上位机控制光源模块的指令发送至光源模块中的的控制中枢。
优选地,所述指令包括了激光模块和白光模块的切换与闭合,采集时间的设置,分析算法的选择,采集区域的选择与调整。
一种基于多模态光学影像技术的检测方法,具体步骤如下:
光源模块分别照射患者的健侧和患侧选定部位,并利用激光模块和白光模块进行成像;
利用激光模块和图像采集模块对产生的激光散斑进行成像并采集,然后利用白光模块和图像采集模块对患者部位的常规照射光进行吸收成像;
图像采集模块对采集到的散斑图像和经过吸收后返回的光学图像发送至上位机,上位机对光学图像进行处理,获得血流和组织含量的信息。
实施本发明的技术方案可解决现有技术中乳腺癌相关淋巴水肿检测中检测效果不好的技术问题;实施本发明的技术方案,通过数字图像处理的相关计算方法分析散斑图像,得到患者照射区域的时间连续的相对位移,再对相对位移信号进行滤波等处理,就可以完成血流灌注量信号的测量,可实现无接触检测、无损检测的技术效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一种实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
其中相同的零部件用相同的附图标记表示。需要说明的是,下面描述中使用的词语“前”、“后”、“左”、“右”、“上”和“下”指的是附图中的方向,词语“底面”和“顶面”、“内”和“外”分别指的是朝向或远离特定部件几何中心的方向。
图1为检测系统的结构简图;
图2为检测方法的流程简图;
图3为使用具体检测系统的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例
在一种具体的实施例中,如图1、图2和图3所示,一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统,包括光源模块,图像采集模块和上位机;所述图像采集模块连接所述光源模块和所述上位机;所述光源模块包括激光模块和白光模块;所述光源模块还包括散光装置。
本实施例中图像采集模块包括图像传感器和成像镜头。
散光片模块用于对集中的光束进行散射,使其均匀的照射人体体表;本实施例中采用的散光模块为三个散光片。
本实施例中的图像采集模块为Basler相机,激光模块为830nm激光机,白光模块,2个单波长的LED光源,波长为LW450和LW595/590。
使用830nm波长的激光一方面的原因是要测量患者的局部血流微循环,另一方面是利用其较白光强的穿透性,有助于白光数据校准。
所述激光模块,用于照射人体健侧(健康侧,以下简称健侧)和患侧(水肿侧,以下简称患侧)部位,产生的激光散斑效应来计算血流灌注情况,利用激光光源照明,使用图像采集模块拍摄图像,通过上位机获取图像R通道的数据,基于激光散斑效应来计算血流灌注情况,
其中IA代表激光光源照明下经过吸收和散射的反射光的光强度,K表示衬比度值,其平方与血流灌注成反比。
利用光源模块中的两个单波长LED光源单独照明,使用图像采集模块拍摄图像,通过上位机获取图像RGB三通道的数据,结合衬比度值数据,基于下述最优化目标计算白光校准矩阵S:
其中代表两个LED单独照明时相机记录的反射光的光强度,代表入射光(两个LED单独照明)强度,F代表相机RGB三个通道的滤光片光谱特性,G为两个LED光源分别对应的光谱矩阵,a和b为数值参数,为l1秩的平方,为l2秩的平方,算子定义如下:
将两个LED单独照明并分别计算得到的白光校准矩阵S1和S2用于修正白光照明的图像数据:
其中,I表示记录到的白光反射光强数据,Iwhite表示修正后的白光数据。
利用光源模块中的白光光源照明,使用图像采集模块拍摄图像,通过上位机获取图像RGB三通道的数据,结合修正矩阵获得修正后的白光数据Iwhite,使用比尔兰伯特定律实现组织含量的计算,比尔兰伯特定律计算组织含量的基本定义如下:
其中Iwhite表示修正后的白光反射图像数据,Iw0表示入射白光光强,L表示光程,μλ表示吸收与散射的结合系数;
图像采集模块,用于对产生的散斑和组织吸收后的图像进行采集。图像采集模块将采集到的散斑和经过吸收后返回的光学图像传送至上位机,且将上位机控制光源模块的指令发送至光源模块中的的控制中枢。指令包括了激光模块和白光模块的切换与闭合,采集时间的设置,分析算法的选择,采集区域的选择与调整。
上位机与图像采集模块连接,用于对采集到的图像进行计算处理,获得被测部位的局部血流微循环和组织含量等数值和图像信息。其中处理方式包括数字图像处理和信号处理进行局部血流微循环的测量,对于经过吸收后返回的光学图像,处理方法是利用人体内组织的光谱吸收特性结合比尔兰伯特定律,进行组织含量的测量,对修正后的白光数据提取R,G,B三个通道的像素数据,利用如下的公式:
其中,IR,IG,IB,分别表示反射光在三个通道的光强度数值,Iwhite表示白光的发射光强度,φ1,θ1分别表示相机RGB三通道的光电响应系数,CLipid分别表示组织中的水,不含氧血红蛋白,含氧血红蛋白和脂肪的含量,α,γ,δ,β分别表示上述四种组织在相应波长下的吸收系数,L1表示白光的光程。
本实施例的检测方法具体步骤如下:
光源模块分别照射患者的健侧和患侧选定部位,并利用激光模块和白光模块进行成像。
利用激光模块和图像采集模块对产生的激光散斑进行成像并采集,然后利用白光模块和图像采集模块对患者部位的常规照射光进行吸收成像。
图像采集模块对采集到的散斑图像和经过吸收后返回的光学图像发送至上位机,上位机对光学图像进行处理,获得血流和组织含量的信息。
本实施例中的局部血流微循环灌注测量基于激光散斑原理,通过数字图像处理的相关计算方法分析散斑图像,得到患者照射区域的时间连续的相对位移,再对相对位移信号进行滤波等处理,就可以完成血流灌注量信号的测量,具有无接触检测、无损检测等特点。
本实施例中的皮下组织含量的测量基于比尔兰伯特定律,通过提取拍摄图片RGB三通道的数据,获取健侧和患侧部位所采集图像各像素点的数值,结合比尔兰伯特方程的变形式和不同物质的吸收系数,完成组织含量的计算。
本实施例避免了与患者检测部位的直接接触,并且光学检测原理受外界环境影响小,所以具有灵敏度高、可靠性高、抗干扰能力强等特点。
本发明与以往发明不同,利用同一台设备既可以检测患者健侧和患侧局部血流灌注情况,又可以利用多波长产生的数据对被测部位的组织含量进行计算,得到评估水肿最直接而有效的数据,而且可以根据测量部位和测量情形的不同,选择不同的器件组合。
本实施例针对传统水肿评估技术的缺陷,提供了一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统与方法,该系统不仅可以完成局部血流微循环灌注量的测量,而且实现了对皮下组织含量(主要包括水,脂肪,蛋白质)的非接触式测量。
需要指出的是,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统,其特征在于:包括光源模块,图像采集模块和上位机;所述图像采集模块连接所述光源模块和所述上位机;所述光源模块包括激光光源、白光光源和两个单波长LED光源;
利用激光光源照明,使用图像采集模块拍摄图像,通过上位机获取图像RGB三通道的数据,基于激光散斑效应来计算血流灌注情况,
其中IA代表激光光源照明下经过吸收和散射的反射光的光强度,K表示衬比度值,其平方与血流灌注成反比;
利用光源模块中的两个单波长LED光源单独照明,使用图像采集模块拍摄图像,通过上位机获取图像RGB三通道的数据,结合衬比度值数据,基于下述最优化目标计算白光校准矩阵S:
其中代表两个单波长LED光源单独照明时图像采集模块记录的反射光的光强度,代表入射光(单波长LED光源单独照明)强度,F代表图像RGB三通道的滤光片光谱特性,G为两个单波长LED光源分别对应的光谱矩阵,a和b为数值参数,为l1秩的平方,为l2秩的平方,算子定义如下:
将两个单波长LED光源单独照明并分别计算得到的白光校准矩阵S1和S2用于修正白光照明的图像数据:
其中,I表示记录到的白光反射光强数据,Iwhite表示修正后的白光数据;
利用光源模块中的白光光源照明,使用图像采集模块拍摄图像,通过上位机获取图像RGB三通道的数据,结合修正矩阵获得修正后的白光数据Iwhite,使用比尔兰伯特定律实现组织含量的计算,比尔兰伯特定律计算组织含量的基本定义如下:
其中Iwhite表示修正后的白光反射图像数据,Iw0表示入射白光光强,L表示光程,μλ表示吸收与散射的结合系数;
所述上位机对采集到的图像进行处理,其中处理方式包括数字图像处理和信号处理进行局部血流微循环的测量,对于经过吸收后返回的光学图像,处理方法是利用人体内组织的光谱吸收特性结合比尔兰伯特定律,进行组织含量的测量,在的图片中提取R,G,B三个通道的像素数据,利用如下的公式:
2.根据权利要求1所述的一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统,其特征在于:所述光源模块还包括散光装置。
3.根据权利要求2所述的一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统,其特征在于:所述激光模块采用波长为830nm的激光器。
4.根据权利要求3所述的一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统,其特征在于:所述白光模块的波段为400-780nm。
5.根据权利要求4所述的一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统,其特征在于:所述两个单波长LED光源其中一个波长为450nm,另一个波长为590nm。
6.根据权利要求5所述的一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统,其特征在于:所述图像采集模块包括图像传感器和成像镜头。
7.根据权利要求6所述的一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统,其特征在于:所述图像采集模块将采集到的散斑和经过吸收后返回的光学图像传送至上位机,且将上位机控制光源模块的指令发送至光源模块中的的控制中枢。
8.根据权利要求7所述的一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统,其特征在于:所述指令包括了激光模块和白光模块的切换与闭合,采集时间的设置,分析算法的选择,采集区域的选择与调整。
9.一种基于多模态光学影像技术的检测方法,其特征在于:采用如权利要求8所述的检测系统,具体步骤如下:
光源模块分别照射患者的健侧和患侧选定部位,并利用激光模块和白光模块进行成像;
利用激光模块和图像采集模块对产生的激光散斑进行成像并采集,然后利用白光模块和图像采集模块对患者部位的常规照射光进行吸收成像;
图像采集模块对采集到的散斑图像和经过吸收后返回的光学图像发送至上位机,上位机对光学图像进行处理,获得血流和组织含量的信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110201906.XA CN112998683B (zh) | 2021-02-23 | 2021-02-23 | 一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统与方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110201906.XA CN112998683B (zh) | 2021-02-23 | 2021-02-23 | 一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统与方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112998683A true CN112998683A (zh) | 2021-06-22 |
CN112998683B CN112998683B (zh) | 2024-02-20 |
Family
ID=76407799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110201906.XA Active CN112998683B (zh) | 2021-02-23 | 2021-02-23 | 一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统与方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112998683B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113974574A (zh) * | 2021-12-15 | 2022-01-28 | 潍坊医学院 | 基于多模态光学的成像系统及成像方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140357990A1 (en) * | 2013-05-29 | 2014-12-04 | University Of Washington Through Its Center For Commercialization | Methods for Laser Speckle Contrast Imaging of Blood Perfusion |
CN104887216A (zh) * | 2015-06-10 | 2015-09-09 | 上海大学 | 多光束相干人体皮肤灌注成像系统和方法 |
CN107122593A (zh) * | 2017-04-06 | 2017-09-01 | 复旦大学 | 基于深度扫描与信息分析的上肢淋巴水肿监测系统 |
CN110780437A (zh) * | 2019-05-10 | 2020-02-11 | 田允娴 | 一种用于人体皮肤上微观位置的重复显微镜成像系统 |
CN111407635A (zh) * | 2020-03-18 | 2020-07-14 | 深圳大学总医院 | 淋巴水肿治疗仪 |
CN111476143A (zh) * | 2020-04-03 | 2020-07-31 | 华中科技大学苏州脑空间信息研究院 | 获取多通道图像、生物多参数以及身份识别的装置 |
CN111481171A (zh) * | 2020-04-03 | 2020-08-04 | 上海交通大学 | 一种用于脑外科手术的多模态监测系统及方法 |
CN112150469A (zh) * | 2020-09-18 | 2020-12-29 | 上海交通大学 | 一种基于无监督领域自适应的激光散斑衬比图像分割方法 |
-
2021
- 2021-02-23 CN CN202110201906.XA patent/CN112998683B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140357990A1 (en) * | 2013-05-29 | 2014-12-04 | University Of Washington Through Its Center For Commercialization | Methods for Laser Speckle Contrast Imaging of Blood Perfusion |
CN104887216A (zh) * | 2015-06-10 | 2015-09-09 | 上海大学 | 多光束相干人体皮肤灌注成像系统和方法 |
CN107122593A (zh) * | 2017-04-06 | 2017-09-01 | 复旦大学 | 基于深度扫描与信息分析的上肢淋巴水肿监测系统 |
CN110780437A (zh) * | 2019-05-10 | 2020-02-11 | 田允娴 | 一种用于人体皮肤上微观位置的重复显微镜成像系统 |
CN111407635A (zh) * | 2020-03-18 | 2020-07-14 | 深圳大学总医院 | 淋巴水肿治疗仪 |
CN111476143A (zh) * | 2020-04-03 | 2020-07-31 | 华中科技大学苏州脑空间信息研究院 | 获取多通道图像、生物多参数以及身份识别的装置 |
CN111481171A (zh) * | 2020-04-03 | 2020-08-04 | 上海交通大学 | 一种用于脑外科手术的多模态监测系统及方法 |
CN112150469A (zh) * | 2020-09-18 | 2020-12-29 | 上海交通大学 | 一种基于无监督领域自适应的激光散斑衬比图像分割方法 |
Non-Patent Citations (3)
Title |
---|
PEI-AN HU: ""Clinical and imaging features of Kaposiform Hemangioendothelioma"", 《THE BRITISH INSTITUTE OF RADIOLOGY》 * |
刘维江: ""基于激光散斑成像的乳腺癌相关淋巴水肿患者上肢血流灌注检测*"", 《工业控制计算机》, vol. 33, no. 6 * |
刘谦: ""利用激光散斑成像监测光动力治疗的血管损伤效应"", 《中国激光》, vol. 32, no. 6 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113974574A (zh) * | 2021-12-15 | 2022-01-28 | 潍坊医学院 | 基于多模态光学的成像系统及成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112998683B (zh) | 2024-02-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2731409C (en) | Quantitative multi-spectral opto-acoustic tomography (msot) of tissue biomarkers | |
US6735458B2 (en) | Self-calibrating optical imaging system | |
US8788003B2 (en) | Monitoring blood constituent levels in biological tissue | |
Jayachandran et al. | Critical review of noninvasive optical technologies for wound imaging | |
US20100210931A1 (en) | Method for performing qualitative and quantitative analysis of wounds using spatially structured illumination | |
US20210052164A1 (en) | Sensor for tissue measurements | |
CN104997519B (zh) | 基于眼底相机的双波长视网膜血管血氧测量系统 | |
CA2754602A1 (en) | Imaging method | |
US11445947B2 (en) | Methods and apparatus for measuring blood oxygenation of tissue | |
CN112998683B (zh) | 一种基于多模态光学影像技术的乳腺癌术后上肢淋巴水肿的检测系统与方法 | |
KR20180057447A (ko) | 다수 개의 근적외선 영상센서가 구비된 스캐너 및 이를 통해 촬영한 다수 개의 시간에 따른 2차원 생체영상을 시간에 따른 당뇨성 상처 확인용 3차원 생체영상으로 변환하는 방법 | |
EP3383245A1 (en) | Hemoglobin measurement from a single vessel | |
Dremin et al. | Hyperspectral imaging of diabetes mellitus skin complications | |
Johansson et al. | Scanning, non-contact, hybrid broadband diffuse optical spectroscopy and diffuse correlation spectroscopy system | |
KR20190133107A (ko) | 혈관 검출 장치 및 방법 | |
US20110208063A1 (en) | NON-CONTACT FREQUENCY DOMAIN NEAR INFRARED ABSORPTION (fNIR) DEVICE FOR ASSESSING TISSUE DAMAGE | |
Huang et al. | The image analysis of skin tissue irradiated with difference wavelengths of LED sources | |
Chiba et al. | Advanced multispectral image-processing endoscopy system for visualizing two-dimensional hemoglobin saturation and relative hemoglobin concentration | |
Rovira | Particular aspects of the use of videopolarimetric technology for dermatological study | |
Neidrauer | Assessment of chronic wounds using in vivo diffuse near infrared spectroscopy | |
Alsadi | Diffuse correlation spectroscopy for estimation of coagulation thickness: a phantom study | |
Khalil et al. | Dynamic contact-free continuous-wave diffuse optical tomography system for the detection of vascular dynamics within the foot | |
Kanawade et al. | Monitoring of epithelial capillary density | |
Opilski et al. | Optical diagnostic system for visualization subcutaneous blood vessels | |
CN113974574A (zh) | 基于多模态光学的成像系统及成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |