CN114403840A - 一种基于近红外光谱的深层组织灌注度计算方法 - Google Patents
一种基于近红外光谱的深层组织灌注度计算方法 Download PDFInfo
- Publication number
- CN114403840A CN114403840A CN202210068601.0A CN202210068601A CN114403840A CN 114403840 A CN114403840 A CN 114403840A CN 202210068601 A CN202210068601 A CN 202210068601A CN 114403840 A CN114403840 A CN 114403840A
- Authority
- CN
- China
- Prior art keywords
- signal
- deep tissue
- perfusion
- infrared light
- index
- 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.)
- Pending
Links
- 230000010412 perfusion Effects 0.000 title claims abstract description 36
- 238000002329 infrared spectrum Methods 0.000 title claims abstract description 8
- 238000004364 calculation method Methods 0.000 title claims abstract description 7
- 238000010521 absorption reaction Methods 0.000 claims abstract description 18
- 238000001914 filtration Methods 0.000 claims abstract description 18
- 230000003287 optical effect Effects 0.000 claims description 11
- 238000000034 method Methods 0.000 claims description 7
- 239000008280 blood Substances 0.000 abstract description 15
- 210000004369 blood Anatomy 0.000 abstract description 15
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 abstract description 13
- 229910052760 oxygen Inorganic materials 0.000 abstract description 13
- 239000001301 oxygen Substances 0.000 abstract description 13
- 108010054147 Hemoglobins Proteins 0.000 description 11
- 102000001554 Hemoglobins Human genes 0.000 description 11
- 238000012544 monitoring process Methods 0.000 description 7
- 239000000523 sample Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000017531 blood circulation Effects 0.000 description 4
- 108010064719 Oxyhemoglobins Proteins 0.000 description 3
- 210000001367 artery Anatomy 0.000 description 3
- 230000008081 blood perfusion Effects 0.000 description 3
- 210000001736 capillary Anatomy 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 210000004204 blood vessel Anatomy 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 210000001061 forehead Anatomy 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 238000004497 NIR spectroscopy Methods 0.000 description 1
- 239000002250 absorbent Substances 0.000 description 1
- 230000002745 absorbent Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 210000003128 head Anatomy 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000000541 pulsatile effect Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000024883 vasodilation Effects 0.000 description 1
- 210000003462 vein Anatomy 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/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0075—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by spectroscopy, i.e. measuring spectra, e.g. Raman spectroscopy, infrared absorption spectroscopy
-
- 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/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- 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/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Cardiology (AREA)
- Physiology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Optics & Photonics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Hematology (AREA)
- Vascular Medicine (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明公开了一种基于近红外光谱的深层组织灌注度计算方法,包括以下步骤:S1:对近端红外光信号和远端红外光信号进行低通滤波,得到信号和所述近端红外光信号和远端红外光信号的波长λ相同,均采用等吸收点的波长;S2:根据所述信号和计算ΔOD(t,λ)信号;S3:对所述ΔOD(t,λ)信号做带通滤波得到信号D;S4:对所述信号做自相关得到信号H;S5:对所述信号H做N点FFT变换,并得到最大能量点索引值index;S6:对所述信号D做N点FFT变换,并在频域取所述最大能量点索引值index处的能量值,所述能量值即为所述深层组织灌注度。本发明能够有效排除血氧饱和度变化对灌注度指标的影响,获得灌注度的绝对值。
Description
技术领域
本发明涉及灌注度测试技术领域,特别涉及一种基于近红外光谱的深层组织灌注度计算方法。
背景技术
灌注,是指流经一个器官或组织体的毛细循环血流速率。为了在不同大小的器官间做个衡量,灌注指标一般被归一化为除以总血流速率。
灌注度是反映动脉血对组织的供血情况的一个指标,临床上更关注灌注指标的相对变化而非绝对值。目前,血流灌注度主要通过计算反射光强的交流成分和直流成分的比值得到。该方法主要测量的是浅层组织毛细血管的灌注度,且容易受到血氧饱和度变化的影响。
发明内容
针对上述问题,本发明旨在提供一种基于近红外光谱的深层组织灌注度计算方法,利用810nm处光波对氧合血红蛋白和还原血红蛋白吸收系数相等的特点,能有效排除血氧饱和度变化对灌注度指标的影响。
本发明的技术方案如下:
一种基于近红外光谱的深层组织灌注度计算方法,包括以下步骤:
S3:对所述ΔOD(t,λ)信号做带通滤波得到信号D;
S5:对所述信号H做N点FFT变换,并得到最大能量点索引值index;
S6:对所述信号D做N点FFT变换,并在频域取所述最大能量点索引值index处的能量值,所述能量值即为所述深层组织灌注度。
作为优选,步骤S1中,进行低通滤波时,滤波器截至频带为4Hz。
作为优选,步骤S1中,所述等吸收点的波长为810nm或805nm。
作为优选,步骤S2中,所述ΔOD(t,λ)信号通过下式进行计算:
作为优选,步骤S3中,进行带通滤波时,滤波器频带为0.5-4Hz。
本发明的有益效果是:
本发明利用等吸收点处光波对氧合血红蛋白和还原血红蛋白吸收系数相等的特点,能有效排除血氧饱和度变化对灌注度指标的影响,获得灌注度的绝对值;且利用ΔOD(t,λ)的空间差分原理能够检测深层组织的灌注度,使其不受表层血流的影响。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一个具体实施例的血氧监测探头监测的30s人体前额数据结果示意图;
图2为图1数据进行低通滤波的结果示意图;
图3为图2数据计算光密度差值ΔOD(t,λ)的结果示意图;
图4为图3数据进行带通滤波的结果示意图;
图7为对信号H做N点FFT变换的结果示意图;
图8为对信号D做N点FFT变换的结果示意图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的技术特征可以相互结合。需要指出的是,除非另有指明,本申请使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。本发明公开使用的“包括”或者“包含”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同,而不排除其他元件或者物件。
本发明提供一种基于近红外光谱的深层组织灌注度计算方法,包括以下步骤:
在一个具体的实施例中,进行低通滤波时,滤波器截至频带为4Hz;所述等吸收点的波长为810nm或805nm。
S3:对所述ΔOD(t,λ)信号做带通滤波得到信号D;
在一个具体的实施例中,进行带通滤波时,滤波器频带为0.5-4Hz。
S5:对所述信号H做N点FFT变换,并得到最大能量点索引值index。
S6:对所述信号D做N点FFT变换,并在频域取所述最大能量点索引值index处的能量值,所述能量值即为所述深层组织灌注度。
需要说明的是,本发明的步骤中,对信号做自相关、N点FFT变换等均为现有技术,具体步骤在此不再赘述。
本发明的基本原理主要是基于血氧监测探头,所述血氧监测探头由一个两波长LED灯和距离LED分别r1和r2(r2>r1)远的两个光接收器组成,LED和两个光接收器位于同一平面同一直线上。LED分时发出波长分别为λ1和λ2的光,远近两个光接收器分时分别得到两种波长的信号。需要说明的是,所述血氧监测探头可使用本公司发明的多通道组织血氧检测探头,也可采用现有技术中其他具有上述结构的血氧监测探头。
在本发明中,所述血氧监测探头的远端和近端波长均采用等吸收点的波长,利用等吸收点处光波对氧合血红蛋白和还原血红蛋白吸收系数相等的特点,有效排除血氧饱和度变化对灌注度指标的影响,获取灌注度的绝对值。
血液(动脉、静脉和毛细血管)中,主要的吸收物质为Hb和HbO2,定义吸光系数为μa(λ):
式中:εHbO2和εHb分别为氧合血红蛋白和还原血红蛋白在波长λ处的吸光系数;CHbO2和CHb分别为氧合血红蛋白和还原血红蛋白的浓度。
则修正的朗伯-比尔定律可以简写为:
OD(t,λ)=μa(λ)*r*DPFλ+G(λ) (3)
式中:OD(t,λ)为波长λ经过组织反射得到的光密度;r为光接收器与发光点的直线距离;DPFλ为差分路径因子;G(λ)为背景吸收。
远近两个光接收器做空间差分即可得到深层组织的光密度:
ΔOD(t,λ)=OD2-OD1=μa(λ)*Δr*DPFλ+ΔG(λ) (4)
动脉血管舒张和收缩会引起光通过血管的长度呈现周期性脉动变化,因此,实际光路由两部分构成:
L=L0+Lac (5)
式中:L0为没有脉搏时,光通过血管的有效光路;Lac为动脉搏动引起的光路变化量;
由此,公式(4)可重新写为:
ΔOD(t,λ)=μa(λ)*(L0+Lac)+ΔG(λ) (6)
根据公式(6)可知,μa(λ)即为经过脉搏波(PPG)调制的波幅值,从ΔOD(t,λ)中分离出交流分量即为:
ΔODac(t,λ)=μa(λ)*Lac (7)
血液中主要有氧合血红蛋白和还原两种吸收体,结合(2)、(7)式可知,ΔODac(t,λ)的幅度主要由总吸光系数μa(λ)和动脉的脉动幅度决定。脉动幅度大小反映了局部组织的灌注强弱,因此可作为局部组织灌注度的指标。总吸光系数由氧合血红蛋白和还原血红蛋白两种物质的吸光系数加权得到,而在等吸收点两种物质的权重相等,因此总的吸光系数不随氧饱和度的变化而变化。
另外,通过上述检测原理,还能使得本发明能够检测深层组织的灌注度,使其不受表层血流的影响。
在一个具体的实施例中,以某人体为例,采用本发明所述的基于近红外光谱的深层组织灌注度计算方法计算该人体的深层脑前叶组织血流灌注度,具体包括以下步骤:
(1)利用血氧监测探头监测等吸收点810nm处该人体在平静状态下,30s的人体前额数据,结果如图1所示。
(3)根据公式(1)计算得到光密度差值ΔOD(t,λ),结果如图3所示。
(4)对所述光密度差值ΔOD(t,λ)做带通滤波得到信号D,滤波器截至频带为[0.5-4Hz];结果如图4所示。
(6)对信号H做N点FFT变换并得到最大幅度点的索引index=0.9766Hz,结果如图7所示。
(7)对步骤(4)得到的信号D做N点FFT变换,并在频域取index=0.9766Hz处的能量值,结果如图8所示。从图8可以看出该处的能量值PI=1.309,四舍五入得到本实施例该人体的深层脑前叶组织血流灌注度为1.3。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (5)
2.根据权利要求1所述的基于近红外光谱的深层组织灌注度计算方法,其特征在于,步骤S1中,进行低通滤波时,滤波器截至频带为4Hz。
3.根据权利要求1所述的基于近红外光谱的深层组织灌注度计算方法,其特征在于,步骤S1中,所述等吸收点的波长为810nm或805nm。
5.根据权利要求1所述的基于近红外光谱的深层组织灌注度计算方法,其特征在于,步骤S3中,进行带通滤波时,滤波器频带为0.5-4Hz。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210068601.0A CN114403840A (zh) | 2022-01-20 | 2022-01-20 | 一种基于近红外光谱的深层组织灌注度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210068601.0A CN114403840A (zh) | 2022-01-20 | 2022-01-20 | 一种基于近红外光谱的深层组织灌注度计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114403840A true CN114403840A (zh) | 2022-04-29 |
Family
ID=81274634
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210068601.0A Pending CN114403840A (zh) | 2022-01-20 | 2022-01-20 | 一种基于近红外光谱的深层组织灌注度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114403840A (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5553615A (en) * | 1994-01-31 | 1996-09-10 | Minnesota Mining And Manufacturing Company | Method and apparatus for noninvasive prediction of hematocrit |
US20090270745A1 (en) * | 2008-02-14 | 2009-10-29 | Yoshiyuki Sankai | Blood flow measuring apparatus and brain activity measuring apparatus using the same |
CN105078441A (zh) * | 2015-09-23 | 2015-11-25 | 广州医软智能科技有限公司 | 人体微循环血流灌注检测仪及检测方法 |
CN107714050A (zh) * | 2017-10-13 | 2018-02-23 | 广东乐心医疗电子股份有限公司 | 一种三波长血氧饱和度检测方法与装置以及可穿戴设备 |
CN109009064A (zh) * | 2018-08-10 | 2018-12-18 | 重庆名希医疗器械有限公司 | 一种颅内压无创监测装置及方法 |
CN113940671A (zh) * | 2021-10-15 | 2022-01-18 | 成都云卫康医疗科技有限公司 | 一种基于近红外光谱的深层组织血氧饱和度计算方法 |
-
2022
- 2022-01-20 CN CN202210068601.0A patent/CN114403840A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5553615A (en) * | 1994-01-31 | 1996-09-10 | Minnesota Mining And Manufacturing Company | Method and apparatus for noninvasive prediction of hematocrit |
US20090270745A1 (en) * | 2008-02-14 | 2009-10-29 | Yoshiyuki Sankai | Blood flow measuring apparatus and brain activity measuring apparatus using the same |
CN105078441A (zh) * | 2015-09-23 | 2015-11-25 | 广州医软智能科技有限公司 | 人体微循环血流灌注检测仪及检测方法 |
CN107714050A (zh) * | 2017-10-13 | 2018-02-23 | 广东乐心医疗电子股份有限公司 | 一种三波长血氧饱和度检测方法与装置以及可穿戴设备 |
CN109009064A (zh) * | 2018-08-10 | 2018-12-18 | 重庆名希医疗器械有限公司 | 一种颅内压无创监测装置及方法 |
CN113940671A (zh) * | 2021-10-15 | 2022-01-18 | 成都云卫康医疗科技有限公司 | 一种基于近红外光谱的深层组织血氧饱和度计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7001337B2 (en) | Monitoring physiological parameters based on variations in a photoplethysmographic signal | |
US6709402B2 (en) | Apparatus and method for monitoring respiration with a pulse oximeter | |
US4955379A (en) | Motion artefact rejection system for pulse oximeters | |
EP1861000B1 (en) | Method and device for determining the perfusion of blood in a body member | |
US5431159A (en) | Pulse oximetry | |
Addison et al. | A novel time–frequency-based 3D Lissajous figure method and its application to the determination of oxygen saturation from the photoplethysmogram | |
IL163173A (en) | Venous pulse oximetry | |
US6658277B2 (en) | Signal processing method and device for signal-to-noise improvement | |
AU2009298937B2 (en) | System and method for photon density wave pulse oximetry and pulse hemometry | |
US20050004479A1 (en) | Locating features in a photoplethysmograph signal | |
JP2003532107A (ja) | 分光光度法により血液酸素添加を非観血的にモニターする方法 | |
KR101746307B1 (ko) | 혈액 중의 산소 포화도 측정 방법 및 시스템 | |
JP2004135854A (ja) | 反射型光電脈波検出装置および反射型オキシメータ | |
WO2011014297A1 (en) | Patient monitoring system and method utilising photon density waves in transmission mode | |
CN113940671B (zh) | 一种基于近红外光谱的深层组织血氧饱和度计算方法 | |
CN108185992B (zh) | 一种无创光学脑组织氧代谢的测量方法 | |
US20090030296A1 (en) | Predictive oximetry model and method | |
CN104739425A (zh) | 混合和中央静脉氧饱和度的光学无创检测方法 | |
CN114403840A (zh) | 一种基于近红外光谱的深层组织灌注度计算方法 | |
Wukitsch | Pulse oximetry: historical review and Ohmeda functional analysis | |
Nitzan et al. | Measurement of the variability of the skin blood volume using dynamic spectroscopy | |
WO2019026062A1 (en) | METHOD FOR MEASURING OXYGEN SATURATION OF ARTERIAL BLOOD AND APPARATUS THEREFOR | |
RU2731414C1 (ru) | Способ комплексной оценки состояния артериального русла | |
Santha et al. | Design considerations of small size reflective type pulse oximeter heads in special applications | |
Lennmarken et al. | Advances in pulse oximetry |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20220429 |
|
RJ01 | Rejection of invention patent application after publication |