CN111493853A - 一种针对血管性皮肤病的血管参数评价方法及系统 - Google Patents
一种针对血管性皮肤病的血管参数评价方法及系统 Download PDFInfo
- Publication number
- CN111493853A CN111493853A CN202010331727.3A CN202010331727A CN111493853A CN 111493853 A CN111493853 A CN 111493853A CN 202010331727 A CN202010331727 A CN 202010331727A CN 111493853 A CN111493853 A CN 111493853A
- Authority
- CN
- China
- Prior art keywords
- light
- image
- emitted
- blood vessel
- vascular
- 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.)
- Withdrawn
Links
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 53
- 238000011156 evaluation Methods 0.000 title claims abstract description 12
- 201000010099 disease Diseases 0.000 title claims description 13
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 title claims description 13
- 239000013307 optical fiber Substances 0.000 claims abstract description 34
- 230000017531 blood circulation Effects 0.000 claims abstract description 29
- 238000001228 spectrum Methods 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000012545 processing Methods 0.000 claims abstract description 27
- 208000026349 Vascular skin disease Diseases 0.000 claims abstract description 22
- 238000003384 imaging method Methods 0.000 claims abstract description 15
- 230000002792 vascular Effects 0.000 claims abstract description 10
- 238000011002 quantification Methods 0.000 claims abstract description 6
- 208000031050 Skin vascular disease Diseases 0.000 claims abstract description 4
- 239000000835 fiber Substances 0.000 claims description 21
- 230000003287 optical effect Effects 0.000 claims description 18
- FCKYPQBAHLOOJQ-UHFFFAOYSA-N Cyclohexane-1,2-diaminetetraacetic acid Chemical compound OC(=O)CN(CC(O)=O)C1CCCCC1N(CC(O)=O)CC(O)=O FCKYPQBAHLOOJQ-UHFFFAOYSA-N 0.000 claims description 11
- 210000004088 microvessel Anatomy 0.000 claims description 7
- 230000009286 beneficial effect Effects 0.000 claims description 6
- 230000003068 static effect Effects 0.000 claims description 6
- 238000004445 quantitative analysis Methods 0.000 claims description 4
- 230000001678 irradiating effect Effects 0.000 claims 1
- 208000017520 skin disease Diseases 0.000 abstract description 2
- 206010048768 Dermatosis Diseases 0.000 abstract 1
- 230000015572 biosynthetic process Effects 0.000 abstract 1
- 238000003759 clinical diagnosis Methods 0.000 abstract 1
- 238000000605 extraction Methods 0.000 abstract 1
- 210000003491 skin Anatomy 0.000 description 25
- 238000012014 optical coherence tomography Methods 0.000 description 16
- 210000001519 tissue Anatomy 0.000 description 12
- 208000003163 Cavernous Hemangioma Diseases 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 238000002428 photodynamic therapy Methods 0.000 description 4
- 208000006787 Port-Wine Stain Diseases 0.000 description 3
- 206010043189 Telangiectasia Diseases 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 201000011066 hemangioma Diseases 0.000 description 3
- 230000001575 pathological effect Effects 0.000 description 3
- 231100000241 scar Toxicity 0.000 description 3
- 208000009056 telangiectasis Diseases 0.000 description 3
- 208000032544 Cicatrix Diseases 0.000 description 2
- 206010067193 Naevus flammeus Diseases 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000002583 angiography Methods 0.000 description 2
- 208000002026 familial multiple nevi flammei Diseases 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000037387 scars Effects 0.000 description 2
- 231100000216 vascular lesion Toxicity 0.000 description 2
- 206010071368 Psychological trauma Diseases 0.000 description 1
- 238000001574 biopsy Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000002845 discoloration Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 210000002615 epidermis Anatomy 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 238000002329 infrared spectrum Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000010412 perfusion Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 239000003229 sclerosing agent Substances 0.000 description 1
- 210000000434 stratum corneum Anatomy 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 238000012795 verification Methods 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Cardiology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Hematology (AREA)
- Physiology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
一种针对血管性皮肤病的血管参数评价方法，包括血管扫描、形成干涉光谱、数据采集、成像处理、参数量化及提取、生成血流图像及得到临床数据；实现该方法的系统由宽带红外光源、红光光源、环形器、光纤耦合器、参考臂、扫描装置、光谱采集模块及信号处理模块构成；系统结构简单，方法容易实现，得到的数据精度高，有助于医生对血管性皮肤病的临床诊断。
Description
(一)技术领域：
本发明涉及医学成像相关技术领域，尤其是一种针对血管性皮肤病的血管参数评价方法及系统。
(二)背景技术：
血管性皮肤病是皮肤美容临床常见疾病之一。最重要的是血管性皮损会使患者产生严重的心理创伤。常见的血管性皮肤疾病主要包括：鲜红斑痣、草莓状血管瘤、海绵状血管瘤、混合性血管瘤和毛细血管扩张症等。治疗方法包括手术切除法、电烙术、冷冻、注射硬化剂等,但这些方法对广泛性、较深的血管性皮损难以控制且副作用较多。解放军总医院激光医学科顾瑛教授发明的光动力疗法(photodynamic therapy，PDT)既能有效消除病变，又不遗留疤痕，疗效好，副作用小。该疗法于1991年1月在临床应用并取得成功，经数千例临床治疗验证，有效率达98％以上。治疗后的病变色可完全消褪，无疤痕，在十余年的随访中未见复发。目前，光动力疗法治疗血管性皮肤疾病的疗效明显优于国内外现有的其它各种治疗手段，是一种理想的治疗方法。
经过近20年对该项技术从基础到临床的系统研究发现，PDT治疗血管性皮肤疾病的疗效主要取决于几个关键治疗参数即：血管直径、血管深度、血管密度及表皮层厚度。获取这三个关键治疗参数可以通过病理学活检，但其属于有创检测，活检将遗留斑痕而难于被患者接受。因此，发展一种适合于血管性皮肤疾病病理特点的无创、高精度快速成像手段一直为临床所期待。
随着光学相干断层(Optical Coherence Tomography，简称OCT)成像技术的问世，OCT以其非侵入性，高分辨和快速成像的特点在诸多领域得到了成功应用。由于OCT具备了在皮肤中穿透一定深度的成像能力，为血管性皮肤病的无创检测提供了一项新的可能的光学技术，但由于OCT只能提供结构信息，不能对血管结构实现高精度成像，在获得微细血管的血流灌注信息及血流参数信息的方面却力不从心。近年来，随着光学微血管造影技术与OCT技术的突破，OCT微血管造影(Optical Coherence Tomography Angiography，简称OCT-A)实现了无须造影剂、高空间分辨率和快速成像的完美结合，有望为微血管疾病的诊断带来革命性的突破，具有非常广阔的应用前景。
(三)发明内容：
本发明的目的在于提供一种针对血管性皮肤病的血管参数评价方法及系统，它可以弥补现有技术的不足，是一种基于OCT-A技术的方法和系统，不但可以对血流进行三维成像，同时还可以对OCT-A图像进行定量描述，可实现血管的参数量化处理，为血管性皮肤病的治疗和监测提供临床指导，其系统结构简单、操作方便。
本发明的技术方案：一种针对血管性皮肤病的血管参数评价方法，其特征在于它包括以下步骤：
(1)光源发出的光经环形器后被耦合器一分为二，一束光进入参考臂，经反射后按照原光路返回成为参考光；发射到光纤准直器，经光纤准直器出射的光发射到扫描振镜；扫描振镜的出射光经物镜聚焦照射于皮肤下微血管处；光束经反射或散射形成测量光，沿着原光路返回到光纤准直器；
(2)光纤准直器出射的测量光与参考光同时发射到所述光纤耦合器发生干涉；干涉光经环形器被光谱采集模块采集，并传输到信号处理模块进行成像处理和参数量化；
(3)OCT系统采集装置在人体组织内血管的同一个位置采集不少于5帧的图像数据，且每帧图像数据中都包括该图像的振幅分量A和相位分量ω；对两个参数进行提取；
(4)根据每幅图像的振幅分量差值ΔA、相位分量差值Δω或两种差值的和ΔA+Δω，求取血流信息，生成OCTA血流图像；
(5)对步骤(4)得到的OCTA血流图像进行量化处理，得到有助于医生客观评估疾病进展和检测治疗的血管密度(Vessel Area Density，简称VAD)、血管骨架密度(VesselSkeleton Density，VSD)、血管直径指数(Vessel Diameter Index，VDI)、血管周长指数(Vessel Perimeter Index，VPI)和血管复杂指数(Vessel Complexity Index，VCI)，对这几个参数进行量化分析，有利于医生对于疾病的客观诊断。
所述步骤(4)中求取血流信息及生成OCTA血流图像的方法具体为：
①通过计算幅度图像相关性的方法提取血流信号；在组织的后向散射光场中，静态组织区域对应的散射光场是保持不变的，具有很强的相关性；而动态血流区域对应的散射光场是随时间变化的，相关性很弱；
②对血管中的同一位置进行至少两次的连续B扫描，根据连续两次B扫图像中局部区域的相关性，即可将组织中的静态信号和动态信号进行分离；
③由于在一帧图像中，无样品信号的区域内的信号幅度值随机变化，因而具有很弱的相关性，导致相关组织结构图像中没有信号的区域具有很强的背景噪声；
④对组织结构图像进行强度阈值设置，每一个强度阈值下都可以生成与其对应的二值化图像，并对二值化图像进行叠加处理，达到去除背景噪声的目的。
所述VAD为微血管血流图像中微血管所占面积与图像总面积之比：
其中，A(i,j)表示血管区域占用的像素；X(i,j)表示微血管图像总像素；
所述VSD为血管图中血管所占长度与总面积之比：
其中，S(i,j)表示血管长度占用的像素；
所述VDI为血管区域占用的像素与血管所占长度之比：
所述VPI为血管周长与OCTA图像总面积之比：
其中P(i,j)表示血管周界内的所有像素；
所述VCI的计算方式为：
一种可实现上述方法的针对血管性皮肤病的血管参数评价系统，其特征在于它包括宽带红外光源110、红光光源120、环形器130、光纤耦合器140、参考臂150、扫描装置160、光谱采集模块170及信号处理模块180；所述参考臂150是由准直器I151、透镜I152和反射镜153构成；所述扫描装置160是由准直器II161、透镜II162、物镜163、扫描振镜164和振镜驱动165构成；所述光谱采集系统170是由准直器III171、扩束器172、光栅173、透镜III174和光电传感器件线性CCD相机175构成；
所述宽带红外光源110产生一个光信号，依次通过环形器130和光纤耦合器140分成两路光信号，一束光进入参考臂150后被反射镜153反射按照原光路返回成为参考光；另一束光发射到扫描装置160的光纤准直器II161，经光纤准直器II161出射的光经过透镜II162发射到扫描振镜164；经扫描振镜出射的光经物镜163聚焦照射到皮肤下微血管处，经反射或散射的光沿着原光路返回，即依次发射到物镜163、扫描振镜164和光纤准直器II161；经光纤准直器II161出射的测量光与参考光发射到光纤耦合器140发生干涉，经光纤耦合器140出射的干涉光经环形器130发射到光谱采集模块170；所述光谱采集模块170将采集的干涉信号发送到信号处理模块180中。
所述宽带红外光源110采用中心波长为1310nm、半高全宽60nm、发射功率25mW的红外光源，照射到皮肤上功率为20mW。
所述红光光源120采用光源中心波长为655nm的红外光源。
所述光纤耦合器140采用10/90耦合器。
本发明的优越性：基于OCT-A技术，图像成像清晰度提高，在对血流进行三维成像的基础上，使用各种指标来定量描述OCT-A图像，有利于医生多角度解释OCT-A血管造影图像，帮助简化识别病理血管特征的过程，以便能够客观评估疾病进展和治疗监测；可以对于皮肤下微血管的血管直径、血管深度、血管密度等参数进行量化分析，通过参数分析极大提高了血管性皮肤病，如鲜红斑痣、草莓状血管瘤、海绵状血管瘤、混合性血管瘤和毛细血管扩张症等的诊断准确性。
(四)附图说明：
图1为本发明所涉一种针对血管性皮肤病的血管参数评价方法的工作流程示意图。
图2为本发明所涉一种针对血管性皮肤病的血管参数评价系统的结构示意图。
其中，110为宽带红外光源；120为红光光源；130为环形器；140为光纤耦合器；150为参考臂；151为准直器I；152为透镜I；153为反射镜；160为扫描装置；161为准直器II；162为透镜II；163为物镜；164为扫描振镜；165为振镜驱动；170为光谱采集模块；171为准直器III；172为扩束器；173为光栅；174为透镜III；175为光电传感器件线性CCD相机；180为信号处理模块。
(五)具体实施方式：
实施例：一种针对血管性皮肤病的血管参数评价方法，如图1所示，其特征在于它包括以下步骤：
(1)光源发出的光经环形器后被耦合器一分为二，一束光进入参考臂，经反射后按照原光路返回成为参考光；发射到光纤准直器，经光纤准直器出射的光发射到扫描振镜；扫描振镜的出射光经物镜聚焦照射于皮肤下微血管处；光束经反射或散射形成测量光，沿着原光路返回到光纤准直器；
(2)光纤准直器出射的测量光与参考光同时发射到所述光纤耦合器发生干涉；干涉光经环形器被光谱采集模块采集，并传输到信号处理模块进行成像处理和参数量化；
(3)OCT系统采集装置在人体组织内血管的同一个位置采集不少于5帧的图像数据，且每帧图像数据中都包括该图像的振幅分量A和相位分量ω；对两个参数进行提取；
(4)根据每幅图像的振幅分量差值ΔA、相位分量差值Δω或两种差值的和ΔA+Δω，求取血流信息，生成OCTA血流图像；
(5)对步骤(4)得到的OCTA血流图像进行量化处理，得到有助于医生客观评估疾病进展和检测治疗的血管密度(Vessel Area Density，简称VAD)、血管骨架密度(VesselSkeleton Density，VSD)、血管直径指数(Vessel Diameter Index，VDI)、血管周长指数(Vessel Perimeter Index，VPI)和血管复杂指数(Vessel Complexity Index，VCI)，对这几个参数进行量化分析，有利于医生对于疾病的客观诊断。
所述步骤(4)中求取血流信息及生成OCTA血流图像的方法具体为：
①通过计算幅度图像相关性的方法提取血流信号；在组织的后向散射光场中，静态组织区域对应的散射光场是保持不变的，具有很强的相关性；而动态血流区域对应的散射光场是随时间变化的，相关性很弱；
②对血管中的同一位置进行至少两次的连续B扫描，根据连续两次B扫图像中局部区域的相关性，即可将组织中的静态信号和动态信号进行分离；
③由于在一帧图像中，无样品信号的区域内的信号幅度值随机变化，因而具有很弱的相关性，导致相关组织结构图像中没有信号的区域具有很强的背景噪声；
④对组织结构图像进行强度阈值设置，每一个强度阈值下都可以生成与其对应的二值化图像，并对二值化图像进行叠加处理，达到去除背景噪声的目的。
所述VAD为微血管血流图像中微血管所占面积与图像总面积之比：
其中，A(i,j)表示血管区域占用的像素；X(i,j)表示微血管图像总像素；
所述VSD为血管图中血管所占长度与总面积之比：
其中，S(i,j)表示血管长度占用的像素；
所述VDI为血管区域占用的像素与血管所占长度之比：
所述VPI为血管周长与OCTA图像总面积之比：
其中P(i,j)表示血管周界内的所有像素；
所述VCI的计算方式为：
一种可实现上述方法的针对血管性皮肤病的血管参数评价系统，其特征在于它包括宽带红外光源110、红光光源120、环形器130、光纤耦合器140、参考臂150、扫描装置160、光谱采集模块170及信号处理模块180；所述参考臂150是由准直器I151、透镜I152和反射镜153构成；所述扫描装置160是由准直器II161、透镜II162、物镜163、扫描振镜164和振镜驱动165构成；所述光谱采集系统170是由准直器III171、扩束器172、光栅173、透镜III174和光电传感器件线性CCD相机175构成；
所述宽带红外光源110产生一个光信号，依次通过环形器130和光纤耦合器140分成两路光信号，一束光进入参考臂150后被反射镜153反射按照原光路返回成为参考光；另一束光发射到扫描装置160的光纤准直器II161，经光纤准直器II161出射的光经过透镜II162发射到扫描振镜164；经扫描振镜出射的光经物镜163聚焦照射到皮肤下微血管处，经反射或散射的光沿着原光路返回，即依次发射到物镜163、扫描振镜164和光纤准直器II161；经光纤准直器II161出射的测量光与参考光发射到光纤耦合器140发生干涉，经光纤耦合器140出射的干涉光经环形器130发射到光谱采集模块170；所述光谱采集模块170将采集的干涉信号发送到信号处理模块180中。
所述宽带红外光源110采用中心波长为1310nm、半高全宽60nm、发射功率25mW的红外光源，照射到皮肤上功率为20mW。
所述红光光源120采用光源中心波长为655nm的红外光源。
所述光纤耦合器140采用10/90耦合器。
以下结合附图来对本发明进行详细的描绘。然而应当理解，附图的提供仅为了更好地理解本发明，它们不应该理解成对本发明的限制。
如图1所示，本发明的一种针对血管性皮肤病的血管参数评价系统基于OCT血管造影实现，包括宽带红外光源110、红光光源120、环形器130、光纤耦合器140、参考臂150、扫描装置160、光谱采集模块170、信号处理模块180。其中，参考臂150包括准直器I151，透镜II152，反射镜153；扫描装置160包括准直器II161、透镜II162、物镜163、扫描振镜164、振镜驱动165。光谱采集模块170包括准直器III171、扩束器172、光栅173、透镜III174和光电传感器件线性CCD相机175。
本实施例的宽带红外光源110采用光源中心波长为1310nm，半高全宽60nm，发射功率25mW，照射到皮肤上功率为20mW；红光光源120采用光源中心波长为655nm；以此为例，不限于此。宽带红外光源110发出的红外光依次经过环形器130和光纤耦合器140被分为两束光，一束光进入参考臂150后被反射镜153反射按照原光路返回成为参考光。另一束光发射到扫描装置160的光纤准直器II161，经光纤准直器II161出射的光经过透镜II162发射到扫描振镜164，经扫描振镜出射的光经物镜163聚焦照射到皮肤下微血管处，经反射或散射的光沿着原光路返回，即依次发射到物镜163、扫描振镜164和光纤准直器II161，经光纤准直器II161出射的测量光与参考光发射到光纤耦合器140发生干涉，经光纤耦合器140出射的干涉光经环形器130发射到光谱采集模块170，光谱采集模块170将采集的干涉信号发送到信号处理模块180采用现有的插值、色散补偿、傅里叶变换等图像处理方法进行处理得到三维皮肤下微血管图像。通过量化分析血管直径、血管深度、血管密度等来诊断血管性皮肤病。
振镜驱动165和光谱采集模块170同时连接信号处理模块180，以保证扫描和采集的同步性。红光光源120发出的红光与红外宽带光同路，起到指示光的作用，红光功率较小且在返回过程中被环形器过滤掉，所以不会影响红外光谱信息。
下面，我们结合上述针对血管性皮肤病的血管参数评价系统工作流程图进行具体描述。
如图2所示，本实施方式提供了一种针对血管性皮肤病的血管参数评价系统系统工作流程，包括如下步骤：
S100、红光确定扫描的皮肤位置。调节扫描装置160的角度，使出射的指示红光垂直照射于待测皮肤区域。
S200、调试扫描装置与皮肤距离。调节扫描装置160与待测皮肤之间的距离，使由物镜出射的红外光聚焦于待测皮肤下微血管处。
S300、扫描装置采集皮肤微血管信息。宽带红外光源110发出红外光，依次经过环形器130和光纤耦合器140到达扫描装置160，经光纤准直器II161出射的光经过透镜II162发射到扫描振镜164，经扫描振镜出射的光经物镜163聚焦照射到皮肤下微血管；信号处理模块180利用振镜驱动165控制扫描振镜164对皮肤下微血管进行扫描，并接收散射和反射的光按原路返回，与参考光在光纤耦合器140中干涉，形成干涉光谱；同时，信号处理模块180控制光谱采集模块170采集干涉光谱，使扫描和采集达到一致；并且，光谱采集模块170将采集到的干涉光谱信号转化为电信号传递到信号处理模块180进行处理成像和参数量化等。
S400、对采集信息进行提取、处理、成像。由于步骤S300的OCT三维数据既包括待测皮肤下微血管的OCT三维数据，还包括皮肤表层、角质层等的三维数据，因此，我们需要从OCT三维数据中提取出皮肤下微血管。另外，步骤S130的OCT三维数据只是初步采集到的数据，我们利用图像处理技术对其进行了进一步的处理。
S500、对结果进行量化。利用软件处理图像信息，对皮肤下微血管的血管面积密度(VAD)、血管骨架密度(VSD)、血管直径指数(VDI)、血管周长指数(VPI)和血管复杂指数(VCI)等参数进行量化分析，通过参数分析极大提高了血管性皮肤病，如鲜红斑痣、草莓状血管瘤、海绵状血管瘤、混合性血管瘤和毛细血管扩张症等的诊断准确性。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种针对血管性皮肤病的血管参数评价方法，其特征在于它包括以下步骤：
(1)光源发出的光经环形器后被耦合器一分为二，一束光进入参考臂，经反射后按照原光路返回成为参考光；发射到光纤准直器，经光纤准直器出射的光发射到扫描振镜；扫描振镜的出射光经物镜聚焦照射于皮肤下微血管处；光束经反射或散射形成测量光，沿着原光路返回到光纤准直器；
(2)光纤准直器出射的测量光与参考光同时发射到所述光纤耦合器发生干涉；干涉光经环形器被光谱采集模块采集，并传输到信号处理模块进行成像处理和参数量化；
(3)OCT系统采集装置在人体组织内血管的同一个位置采集不少于5帧的图像数据，且每帧图像数据中都包括该图像的振幅分量A和相位分量ω；对两个参数进行提取；
(4)根据每幅图像的振幅分量差值ΔA、相位分量差值Δω或两种差值的和ΔA+Δω，求取血流信息，生成OCTA血流图像；
(5)对步骤(4)得到的OCTA血流图像进行量化处理，得到有助于医生客观评估疾病进展和检测治疗的血管密度VAD、血管骨架密度VSD、血管直径指数VDI、血管周长指数VPI和血管复杂指数VCI，对这几个参数进行量化分析，有利于医生对于疾病的客观诊断。
2.根据权利要求1所述一种针对血管性皮肤病的血管参数评价方法，其特征在于所述步骤(4)中求取血流信息及生成OCTA血流图像的方法具体为：
①通过计算幅度图像相关性的方法提取血流信号；在组织的后向散射光场中，静态组织区域对应的散射光场是保持不变的，具有很强的相关性；而动态血流区域对应的散射光场是随时间变化的，相关性很弱；
②对血管中的同一位置进行至少两次的连续B扫描，根据连续两次B扫图像中局部区域的相关性，即可将组织中的静态信号和动态信号进行分离；
③由于在一帧图像中，无样品信号的区域内的信号幅度值随机变化，因而具有很弱的相关性，导致相关组织结构图像中没有信号的区域具有很强的背景噪声；
④对组织结构图像进行强度阈值设置，每一个强度阈值下都可以生成与其对应的二值化图像，并对二值化图像进行叠加处理，达到去除背景噪声的目的。
3.根据权利要求1所述一种针对血管性皮肤病的血管参数评价方法，其特征在于所述VAD为微血管血流图像中微血管所占面积与图像总面积之比：
其中，A(i,j)表示血管区域占用的像素；X(i,j)表示微血管图像总像素。
4.根据权利要求1所述一种针对血管性皮肤病的血管参数评价方法，其特征在于所述VSD为血管图中血管所占长度与总面积之比：
其中，S(i,j)表示血管长度占用的像素。
5.根据权利要求1所述一种针对血管性皮肤病的血管参数评价方法，其特征在于所述VDI为血管区域占用的像素与血管所占长度之比：
6.根据权利要求1所述一种针对血管性皮肤病的血管参数评价方法，其特征在于所述VPI为血管周长与OCTA图像总面积之比：
其中P(i,j)表示血管周界内的所有像素。
7.根据权利要求1所述一种针对血管性皮肤病的血管参数评价方法，其特征在于所述VCI的计算方式为：
8.一种可实现上述方法的针对血管性皮肤病的血管参数评价系统，其特征在于它包括宽带红外光源110、红光光源120、环形器130、光纤耦合器140、参考臂150、扫描装置160、光谱采集模块170及信号处理模块180；所述参考臂150是由准直器I151、透镜I152和反射镜153构成；所述扫描装置160是由准直器II161、透镜II162、物镜163、扫描振镜164和振镜驱动165构成；所述光谱采集系统170是由准直器III171、扩束器172、光栅173、透镜III174和光电传感器件线性CCD相机175构成；
所述宽带红外光源110产生一个光信号，依次通过环形器130和光纤耦合器140分成两路光信号，一束光进入参考臂150后被反射镜153反射按照原光路返回成为参考光；另一束光发射到扫描装置160的光纤准直器II161，经光纤准直器II161出射的光经过透镜II162发射到扫描振镜164；经扫描振镜出射的光经物镜163聚焦照射到皮肤下微血管处，经反射或散射的光沿着原光路返回，即依次发射到物镜163、扫描振镜164和光纤准直器II161；经光纤准直器II161出射的测量光与参考光发射到光纤耦合器140发生干涉，经光纤耦合器140出射的干涉光经环形器130发射到光谱采集模块170；所述光谱采集模块170将采集的干涉信号发送到信号处理模块180中。
9.根据权利要求8所述一种针对血管性皮肤病的血管参数评价系统，其特征在于所述宽带红外光源110采用中心波长为1310nm、半高全宽60nm、发射功率25mW的红外光源，照射到皮肤上功率为20mW。
10.根据权利要求9所述一种针对血管性皮肤病的血管参数评价系统，其特征在于所述红光光源120采用光源中心波长为655nm的红外光源；所述光纤耦合器140采用10/90耦合器。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010331727.3A CN111493853A (zh) | 2020-04-24 | 2020-04-24 | 一种针对血管性皮肤病的血管参数评价方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010331727.3A CN111493853A (zh) | 2020-04-24 | 2020-04-24 | 一种针对血管性皮肤病的血管参数评价方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111493853A true CN111493853A (zh) | 2020-08-07 |
Family
ID=71848470
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010331727.3A Withdrawn CN111493853A (zh) | 2020-04-24 | 2020-04-24 | 一种针对血管性皮肤病的血管参数评价方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111493853A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114209278A (zh) * | 2021-12-14 | 2022-03-22 | 复旦大学 | 一种基于光学相干层析成像的深度学习皮肤病诊断系统 |
| CN114305321A (zh) * | 2022-03-15 | 2022-04-12 | 汕头大学·香港中文大学联合汕头国际眼科中心 | 一种测量视网膜血管壁厚度的方法及系统 |
| CN114469008A (zh) * | 2022-03-09 | 2022-05-13 | 浙江大学 | 基于血管近红外成像的血管穿刺难易智能分级装置 |
| CN118236034A (zh) * | 2024-03-26 | 2024-06-25 | 深圳市利孚医疗技术有限公司 | 一种基于红外透皮成像的皮下毛发检测技术 |
Citations (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1083340A (zh) * | 1992-08-31 | 1994-03-09 | 林秀列 | 电子诊脉装置及方法 |
| CN2311374Y (zh) * | 1997-09-10 | 1999-03-24 | 张一铭 | 光纤装饰品 |
| WO2002015786A1 (en) * | 2000-08-25 | 2002-02-28 | Cytometrics, Llc | System, method and computer program product for measuring blood properties form a spectral image |
| CN101996329A (zh) * | 2010-11-17 | 2011-03-30 | 沈阳东软医疗系统有限公司 | 一种对血管形变区域的检测装置和方法 |
| CN103732132A (zh) * | 2011-06-13 | 2014-04-16 | 安吉奥梅特里克斯公司 | 用于分析解剖学参数和功能参数的多功能导丝组件和系统 |
| CN204562099U (zh) * | 2014-12-29 | 2015-08-19 | 浙江大学 | 基于复数互相关的微血管光学造影及抖动补偿系统 |
| CN105574855A (zh) * | 2015-12-10 | 2016-05-11 | 南京理工大学 | 云背景下基于模板滤波和虚警抑制的红外小目标检测方法 |
| CN106137134A (zh) * | 2016-08-08 | 2016-11-23 | 浙江大学 | 多角度复合的血流成像方法及系统 |
| CN106539622A (zh) * | 2017-01-28 | 2017-03-29 | 北京欣方悦医疗科技有限公司 | 基于血流动力学分析的冠脉虚拟支架植入方法和系统 |
| CN106815853A (zh) * | 2016-12-14 | 2017-06-09 | 海纳医信（北京）软件科技有限责任公司 | 对眼底图像中视网膜血管的分割方法和装置 |
| CN107137073A (zh) * | 2017-05-27 | 2017-09-08 | 天津优视眼科技术有限公司 | 一种基于甲壁血流分析的青光眼和糖尿病等微循环系统疾病检测技术 |
| CN107307848A (zh) * | 2017-05-27 | 2017-11-03 | 天津海仁医疗技术有限公司 | 一种基于高速大范围扫描光学微造影成像的人脸识别及皮肤检测系统 |
| US20180012359A1 (en) * | 2016-07-06 | 2018-01-11 | Marinko Venci Sarunic | Systems and Methods for Automated Image Classification and Segmentation |
| CN107595250A (zh) * | 2017-09-30 | 2018-01-19 | 浙江大学 | 基于运动与图形混合对比度的血流成像方法与系统 |
| CN107822600A (zh) * | 2017-12-14 | 2018-03-23 | 天津海仁医疗技术有限公司 | 应用于血管性皮肤病检测、定位的装置和系统及工作方法 |
| CN108042125A (zh) * | 2017-05-27 | 2018-05-18 | 天津海仁医疗技术有限公司 | 一种高速内窥光学相干血流成像系统 |
| CN108309244A (zh) * | 2018-02-13 | 2018-07-24 | 天津海仁医疗技术有限公司 | 一种新型压力传感皮肤微血管扫描装置和系统及工作方法 |
| CN108478192A (zh) * | 2018-04-11 | 2018-09-04 | 西安交通大学 | 一种用于估测皮肤组织微血管深度的测量系统 |
| CN109033949A (zh) * | 2018-06-11 | 2018-12-18 | 天津海仁医疗技术有限公司 | 基于光学相干断层扫描的人脸识别方法和系统 |
| CN109171670A (zh) * | 2018-06-25 | 2019-01-11 | 天津海仁医疗技术有限公司 | 一种基于逆向主成分分析法的3d血管成像算法 |
| CN209018702U (zh) * | 2017-12-14 | 2019-06-25 | 天津海仁医疗技术有限公司 | 应用于血管性皮肤病检测、定位的装置和系统 |
| JP2019103762A (ja) * | 2017-12-14 | 2019-06-27 | 株式会社ニデック | Ｏｃｔデータ解析装置、およびｏｃｔデータ解析プログラム |
| CN110057308A (zh) * | 2019-04-30 | 2019-07-26 | 北京理工大学 | 一种形变监测方法及系统 |
| CN110146791A (zh) * | 2019-05-07 | 2019-08-20 | 国网山东省电力公司电力科学研究院 | 一种基于图像处理的电晕检测方法 |
| CN110291556A (zh) * | 2016-09-22 | 2019-09-27 | R·K·王 | 增强光学图像质量的方法和系统 |
| CN110292359A (zh) * | 2019-07-09 | 2019-10-01 | 浙江大学 | 一种无标记全光学神经调控与成像的方法与装置 |
| JP2019209136A (ja) * | 2018-05-31 | 2019-12-12 | キヤノン株式会社 | 情報処理装置、情報処理方法及びプログラム |
-
2020
- 2020-04-24 CN CN202010331727.3A patent/CN111493853A/zh not_active Withdrawn
Patent Citations (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1083340A (zh) * | 1992-08-31 | 1994-03-09 | 林秀列 | 电子诊脉装置及方法 |
| CN2311374Y (zh) * | 1997-09-10 | 1999-03-24 | 张一铭 | 光纤装饰品 |
| WO2002015786A1 (en) * | 2000-08-25 | 2002-02-28 | Cytometrics, Llc | System, method and computer program product for measuring blood properties form a spectral image |
| CN101996329A (zh) * | 2010-11-17 | 2011-03-30 | 沈阳东软医疗系统有限公司 | 一种对血管形变区域的检测装置和方法 |
| CN103732132A (zh) * | 2011-06-13 | 2014-04-16 | 安吉奥梅特里克斯公司 | 用于分析解剖学参数和功能参数的多功能导丝组件和系统 |
| CN204562099U (zh) * | 2014-12-29 | 2015-08-19 | 浙江大学 | 基于复数互相关的微血管光学造影及抖动补偿系统 |
| CN105574855A (zh) * | 2015-12-10 | 2016-05-11 | 南京理工大学 | 云背景下基于模板滤波和虚警抑制的红外小目标检测方法 |
| US20180012359A1 (en) * | 2016-07-06 | 2018-01-11 | Marinko Venci Sarunic | Systems and Methods for Automated Image Classification and Segmentation |
| CN106137134A (zh) * | 2016-08-08 | 2016-11-23 | 浙江大学 | 多角度复合的血流成像方法及系统 |
| CN110291556A (zh) * | 2016-09-22 | 2019-09-27 | R·K·王 | 增强光学图像质量的方法和系统 |
| CN106815853A (zh) * | 2016-12-14 | 2017-06-09 | 海纳医信（北京）软件科技有限责任公司 | 对眼底图像中视网膜血管的分割方法和装置 |
| CN106539622A (zh) * | 2017-01-28 | 2017-03-29 | 北京欣方悦医疗科技有限公司 | 基于血流动力学分析的冠脉虚拟支架植入方法和系统 |
| CN107137073A (zh) * | 2017-05-27 | 2017-09-08 | 天津优视眼科技术有限公司 | 一种基于甲壁血流分析的青光眼和糖尿病等微循环系统疾病检测技术 |
| CN107307848A (zh) * | 2017-05-27 | 2017-11-03 | 天津海仁医疗技术有限公司 | 一种基于高速大范围扫描光学微造影成像的人脸识别及皮肤检测系统 |
| CN108042125A (zh) * | 2017-05-27 | 2018-05-18 | 天津海仁医疗技术有限公司 | 一种高速内窥光学相干血流成像系统 |
| CN107595250A (zh) * | 2017-09-30 | 2018-01-19 | 浙江大学 | 基于运动与图形混合对比度的血流成像方法与系统 |
| CN209018702U (zh) * | 2017-12-14 | 2019-06-25 | 天津海仁医疗技术有限公司 | 应用于血管性皮肤病检测、定位的装置和系统 |
| JP2019103762A (ja) * | 2017-12-14 | 2019-06-27 | 株式会社ニデック | Ｏｃｔデータ解析装置、およびｏｃｔデータ解析プログラム |
| CN107822600A (zh) * | 2017-12-14 | 2018-03-23 | 天津海仁医疗技术有限公司 | 应用于血管性皮肤病检测、定位的装置和系统及工作方法 |
| CN108309244A (zh) * | 2018-02-13 | 2018-07-24 | 天津海仁医疗技术有限公司 | 一种新型压力传感皮肤微血管扫描装置和系统及工作方法 |
| CN108478192A (zh) * | 2018-04-11 | 2018-09-04 | 西安交通大学 | 一种用于估测皮肤组织微血管深度的测量系统 |
| JP2019209136A (ja) * | 2018-05-31 | 2019-12-12 | キヤノン株式会社 | 情報処理装置、情報処理方法及びプログラム |
| CN109033949A (zh) * | 2018-06-11 | 2018-12-18 | 天津海仁医疗技术有限公司 | 基于光学相干断层扫描的人脸识别方法和系统 |
| CN109171670A (zh) * | 2018-06-25 | 2019-01-11 | 天津海仁医疗技术有限公司 | 一种基于逆向主成分分析法的3d血管成像算法 |
| CN110057308A (zh) * | 2019-04-30 | 2019-07-26 | 北京理工大学 | 一种形变监测方法及系统 |
| CN110146791A (zh) * | 2019-05-07 | 2019-08-20 | 国网山东省电力公司电力科学研究院 | 一种基于图像处理的电晕检测方法 |
| CN110292359A (zh) * | 2019-07-09 | 2019-10-01 | 浙江大学 | 一种无标记全光学神经调控与成像的方法与装置 |
Non-Patent Citations (3)
| Title |
|---|
| ALAM, MINHAJ；THAPA, DAMBELIM, JENNIFER I.R；: "Computer-aided classification of sickle cell retinopathy using quantitative features in optical coherence tomography angiography", 《BIOMEDICAL OPTICS EXPRESS》 * |
| CHU, ZD，LIN, JASON，GAO, CHEN: "Quantitative assessment of the retinal microvasculature using optical coherence tomography angiography", 《JOURNAL OF BIOMEDICAL OPTICS》 * |
| 陈朝良: "基于谱域光学相干层析术的人体皮肤三维血流成像研究", 《CNKI博士学位论文全文库》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114209278A (zh) * | 2021-12-14 | 2022-03-22 | 复旦大学 | 一种基于光学相干层析成像的深度学习皮肤病诊断系统 |
| CN114209278B (zh) * | 2021-12-14 | 2023-08-25 | 复旦大学 | 一种基于光学相干层析成像的深度学习皮肤病诊断系统 |
| CN114469008A (zh) * | 2022-03-09 | 2022-05-13 | 浙江大学 | 基于血管近红外成像的血管穿刺难易智能分级装置 |
| CN114305321A (zh) * | 2022-03-15 | 2022-04-12 | 汕头大学·香港中文大学联合汕头国际眼科中心 | 一种测量视网膜血管壁厚度的方法及系统 |
| CN118236034A (zh) * | 2024-03-26 | 2024-06-25 | 深圳市利孚医疗技术有限公司 | 一种基于红外透皮成像的皮下毛发检测技术 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111493853A (zh) | 一种针对血管性皮肤病的血管参数评价方法及系统 | |
| Deegan et al. | Microvascular imaging of the skin | |
| van den Berg et al. | Feasibility of photoacoustic/ultrasound imaging of synovitis in finger joints using a point-of-care system | |
| US7322972B2 (en) | In vivo port wine stain, burn and melanin depth determination using a photoacoustic probe | |
| CN107595250B (zh) | 基于运动与图形混合对比度的血流成像方法与系统 | |
| JP4849755B2 (ja) | イメージング装置および試料分析方法 | |
| US8036727B2 (en) | Methods for noninvasively measuring analyte levels in a subject | |
| KR100490461B1 (ko) | 청진기 | |
| CN100477955C (zh) | 用于光电检测肿瘤的装置 | |
| WO2021003782A1 (zh) | 一种无标记全光学神经调控与成像的方法与装置 | |
| JP2008537995A (ja) | 高分解能空間変調蛍光画像化および断層撮影のための方法ならびに装置 | |
| CN112168144B (zh) | 一种用于烧伤皮肤的光学相干层析成像系统 | |
| US10746706B2 (en) | Photoacoustic physio-chemical tissue analysis | |
| JP2007503224A (ja) | 心拍動の間に血液を介して血管壁をスペクトル検査する方法およびシステム | |
| JPH11514549A (ja) | レーザー光学音響画像形成装置 | |
| CN109662735B (zh) | 皮肤血流灌注量的测量方法 | |
| CN108478192B (zh) | 一种用于估测皮肤组织微血管深度的测量系统 | |
| CN107228904B (zh) | 一种光致超声的血糖无创检测装置及方法 | |
| KR101746763B1 (ko) | 망막 또는 맥락막 내 혈관조영 광가간섭 단층촬영 장치 및 이를 이용한 질병 진단방법 | |
| Jeeva et al. | Reconstruction of optical scanned images of inhomogeneities in biological tissues by Monte Carlo simulation | |
| Wang et al. | Quantitative classification of melasma with photoacoustic microscopy: a pilot study | |
| CN108604377B (zh) | 处理光学相干断层成像扫描图 | |
| CN105249933A (zh) | 光声分子三维成像仪 | |
| CN209474591U (zh) | 一种基于近红外光吸收的手持式紧凑型水肿检测装置 | |
| EP1797818A2 (en) | Method and system for tomographic imaging using fluorescent proteins |
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 | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20200807 |