CN114190893A - 一种断肢再植及游离皮瓣术后血运评估装置及方法 - Google Patents
一种断肢再植及游离皮瓣术后血运评估装置及方法 Download PDFInfo
- Publication number
- CN114190893A CN114190893A CN202111500149.2A CN202111500149A CN114190893A CN 114190893 A CN114190893 A CN 114190893A CN 202111500149 A CN202111500149 A CN 202111500149A CN 114190893 A CN114190893 A CN 114190893A
- Authority
- CN
- China
- Prior art keywords
- indocyanine green
- time
- replantation
- skin flap
- relative
- 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
- 230000017531 blood circulation Effects 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 25
- 210000004369 blood Anatomy 0.000 claims abstract description 24
- 239000008280 blood Substances 0.000 claims abstract description 24
- 239000000049 pigment Substances 0.000 claims abstract description 17
- 238000011156 evaluation Methods 0.000 claims abstract description 15
- 238000003384 imaging method Methods 0.000 claims abstract description 13
- 230000002503 metabolic effect Effects 0.000 claims abstract description 12
- MOFVSTNWEDAEEK-UHFFFAOYSA-M indocyanine green Chemical compound [Na+].[O-]S(=O)(=O)CCCCN1C2=CC=C3C=CC=CC3=C2C(C)(C)C1=CC=CC=CC=CC1=[N+](CCCCS([O-])(=O)=O)C2=CC=C(C=CC=C3)C3=C2C1(C)C MOFVSTNWEDAEEK-UHFFFAOYSA-M 0.000 claims description 84
- 229960004657 indocyanine green Drugs 0.000 claims description 84
- 239000003814 drug Substances 0.000 claims description 45
- 229940079593 drug Drugs 0.000 claims description 27
- 238000002073 fluorescence micrograph Methods 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 16
- 230000004060 metabolic process Effects 0.000 claims description 16
- 230000010412 perfusion Effects 0.000 claims description 16
- 238000013186 photoplethysmography Methods 0.000 claims description 14
- 238000000799 fluorescence microscopy Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 11
- 230000008034 disappearance Effects 0.000 claims description 10
- 238000013139 quantization Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 8
- 210000003462 vein Anatomy 0.000 claims description 8
- 238000012632 fluorescent imaging Methods 0.000 claims description 6
- 230000005622 photoelectricity Effects 0.000 claims description 6
- 125000004122 cyclic group Chemical group 0.000 claims description 5
- 230000002980 postoperative effect Effects 0.000 claims description 5
- 238000001356 surgical procedure Methods 0.000 claims description 5
- 230000002093 peripheral effect Effects 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000002266 amputation Methods 0.000 claims description 2
- 238000001228 spectrum Methods 0.000 abstract description 5
- 230000002411 adverse Effects 0.000 abstract description 2
- 210000003414 extremity Anatomy 0.000 description 30
- 210000001519 tissue Anatomy 0.000 description 22
- 239000000243 solution Substances 0.000 description 8
- 108010054147 Hemoglobins Proteins 0.000 description 6
- 102000001554 Hemoglobins Human genes 0.000 description 6
- 238000002835 absorbance Methods 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 4
- 210000004204 blood vessel Anatomy 0.000 description 4
- 238000010968 computed tomography angiography Methods 0.000 description 4
- 210000004185 liver Anatomy 0.000 description 4
- 238000010992 reflux Methods 0.000 description 4
- 238000002583 angiography Methods 0.000 description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 3
- 230000008081 blood perfusion Effects 0.000 description 3
- 230000036770 blood supply Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000004087 circulation Effects 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 231100001095 no nephrotoxicity Toxicity 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 229910052760 oxygen Inorganic materials 0.000 description 3
- 239000001301 oxygen Substances 0.000 description 3
- 230000002792 vascular Effects 0.000 description 3
- 238000000862 absorption spectrum Methods 0.000 description 2
- 210000001367 artery Anatomy 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- GNBHRKFJIUUOQI-UHFFFAOYSA-N fluorescein Chemical compound O1C(=O)C2=CC=CC=C2C21C1=CC=C(O)C=C1OC1=CC(O)=CC=C21 GNBHRKFJIUUOQI-UHFFFAOYSA-N 0.000 description 2
- 238000013534 fluorescein angiography Methods 0.000 description 2
- 230000000004 hemodynamic effect Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000002285 radioactive effect Effects 0.000 description 2
- 230000008439 repair process Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 239000002904 solvent Substances 0.000 description 2
- 230000004083 survival effect Effects 0.000 description 2
- 201000002282 venous insufficiency Diseases 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 206010027654 Allergic conditions Diseases 0.000 description 1
- 108010017384 Blood Proteins Proteins 0.000 description 1
- 102000004506 Blood Proteins Human genes 0.000 description 1
- 208000005189 Embolism Diseases 0.000 description 1
- 206010014513 Embolism arterial Diseases 0.000 description 1
- 241000282412 Homo Species 0.000 description 1
- 102000036675 Myoglobin Human genes 0.000 description 1
- 108010062374 Myoglobin Proteins 0.000 description 1
- 208000001647 Renal Insufficiency Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 210000003363 arteriovenous anastomosis Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000941 bile Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 210000000624 ear auricle Anatomy 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000295 emission spectrum Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000002189 fluorescence spectrum Methods 0.000 description 1
- 108010036302 hemoglobin AS Proteins 0.000 description 1
- 239000007943 implant Substances 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 210000000936 intestine Anatomy 0.000 description 1
- 238000007918 intramuscular administration Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 230000000302 ischemic effect Effects 0.000 description 1
- 201000006370 kidney failure Diseases 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 238000011542 limb amputation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000004089 microcirculation Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002496 oximetry Methods 0.000 description 1
- 206010033675 panniculitis Diseases 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 210000004304 subcutaneous tissue Anatomy 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 210000003371 toe Anatomy 0.000 description 1
- 238000002054 transplantation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0071—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by measuring fluorescence emission
-
- 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
-
- 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/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
- A61B5/14551—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 for measuring blood gases
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/44—Detecting, measuring or recording for evaluating the integumentary system, e.g. skin, hair or nails
- A61B5/441—Skin evaluation, e.g. for skin disorder diagnosis
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Cardiology (AREA)
- Physiology (AREA)
- Dermatology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Optics & Photonics (AREA)
- Vascular Medicine (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明公开了一种断肢再植及游离皮瓣术后血运评估装置及方法，包括光电容积脉搏波信号采集装置、荧光图像采集装置和工作站。与现有技术相比，本发明以PPG信号中提取的ICG色素浓度谱作为ICG色素在人体内的正常代谢曲线，通过比较ICGA采集的每一像素点处的ICG代谢曲线与人体内正常ICG代谢曲线的差异性，评估断指(肢)再植及游离皮瓣术后血运情况，避免或减小了不利因素对ICGA评估皮瓣血运结果准确性的影响，直接生成断指(肢)再植及游离皮瓣术后血运评估参数，并以图像化显示，具有直观性和客观性，避免外科医生依靠临床经验评估断指(肢)再植及游离皮瓣术后血运情况的主观性。
Description
技术领域
本发明涉及医疗器械领域，特别是一种断肢再植及游离皮瓣术后血运评估装置及方法。
背景技术
断肢(指、趾)再植，是指将完全或不完全离断的肢体在光学放大镜的助视下，重新接回原位，恢复血液循环，使之成活并恢复一定功能的高、精、细手术。皮瓣移植是修复重建外科最基本、最常用的术式之一，是指将本身带有血供的皮肤和皮下组织转移到其他部位，以达到对皮瓣受区部位的修复、再造目的。断指(肢)再植及游离皮瓣术后24-72小时是吻合血管出现血管危象的高发期，发生危象应第一时间进行干预，会提高再植指及游离皮瓣的术后成活率。
临床上医生常常通过观察再植肢体或皮瓣的皮肤颜色、皮肤温度、指腹张力、毛细血管回流时间测定等指标判断再植指及游离皮瓣的血运情况。从颜色上讲，指体色泽红润为正常，指体由红润变为苍白或苍白变为浅灰色，说明断指(肢)缺血，系有动脉危象发生。指体鲜红变为暗红，继而变为暗紫色，说明静脉回流受阻，系有静脉危象发生；从温度上讲，移植组织与健侧组织的皮肤温差在0.5-2.0℃以内呈平行变化，说明动静脉吻合口通畅，移植组织血循环良好。皮肤温差突然相差3℃以上时大多是动脉栓塞，应立即手术探查。皮肤温差逐渐增大，一般24-48小时后达到3℃以上时大多是静脉栓塞；从指腹张力上讲，正常手指指腹饱满而富有弹性，如再植指供血不足，指腹张力立即降低，如静脉回流不畅，则指腹张力增高，指腹张力增高到一定程度时，皮肤会出现张力性水泡；从毛细血管回流时间测定上讲，轻轻按压患指的指腹或指甲，此时被压的皮肤或指甲呈苍白色，一旦移开压迫后，受压区在1-2秒内由苍白转为红润，则说明毛细血管回流正常，而指体苍白且测不出毛细血管回流充盈现象，则指体供血障碍，说明发生动脉危象。
在临床工作中，目前常用的术前评估组织血运(包括血管定位)的医学影像学方法包括各种多普勒仪器，CT血管造影检查(computed tomographic angiography，CTA)，核磁血管造影检查(magnetic resonance angiography，MRA)，组织血氧探测仪，荧光素造影等。便携式超声多普勒(hand-held Doppler，HHD)：是一种安全、便携的仪器，是临床上最常用的术前穿支血管位置探查设备。然而，HHD对检查者具有较强的依赖性，检查结果的准确性与检查者的临床经验等具有较大相关性。彩色超声多普勒(color duplex sonography，CDS)：与HHD相同，也是利用血管回声原理进行血管成像。与HHD相比，CDS可以提供更为详细的与穿支血管的起源、口径、肌内行程和血流动力学有关的信息，同样具有无创、无辐射、经济快捷等优点。CDS的缺点包括探查所需时间较长，对检查者技术的依赖性较大，无法提供穿支血管在深层组织内的走行情况信息。激光多普勒血流监测仪(laser-Dopplerflowmetry，LDF)：是一种能够实时监测组织内微循环血流灌注的仪器，在手术过程中可用于评估皮瓣内的血流情况；然而LDF机身较重，且对检查对象的细微运动较敏感，探测穿支血管的能力较差。CTA及MRA：是较常用的术前穿支血管定位评估的两种方法，在血管可视化方面具有明显优势。然而，两种方法所使用的设备价格昂贵，CTA具有辐射性，为有创检查，对放射科医生技术具有一定依赖性，图像制作及解读需要花费较长时间等；MRA禁用于肾功能不全、过敏体质及体内有金属植入物的患者。经皮血氧饱和度监测仪(transcutaneousoxygen monitoring)：术后应用血氧饱和度监测仪可以帮助外科医生发现皮瓣灌注不良的部位，并及时地对灌注差的皮瓣采取干预措施，从而提高皮瓣的成活率。该设备探头需要与被测部位进行接触探测，操作耗时，且空间分辨率较低。荧光素造影技术(fluoresceinangiography)：可用于术中对组织灌注情况进行评估，该技术的主要缺点为所用荧光素的半衰期较长，清除时间需要12-18小时，因此，术中不能重复多次使用。
吲哚菁绿(indocyanine green，ICG)是一种经FDA批准的可用于临床的水溶性染料。ICG通过外周静脉或中心静脉团注给药进入人体后，98％的吲哚菁绿可以和血液中的血浆蛋白结合，极少渗透至组织间隙，从而可以保留在血管腔内，有较高的特异性，是评估组织灌注的理想染料。ICG主要被肝脏摄取，经胆汁代谢，无肝肠循环，无肾毒性。ICG在体内的代谢分为两个过程：90％的吲哚菁绿在第一阶段成指数规律清除，其余10％在第二阶段排除，总时间约为66min，在同一手术过程中，可以重复使用。血液中吲哚菁绿吸收光谱范围为650～850nm，最大吸收峰为805nm，经760nm的近红外光激发后产生荧光，荧光发射光谱在820～830nm左右，荧光产量随着吲哚菁绿药物浓度的增加而线性增加，直到浓度超过80mg/L时，荧光产量就会下降。小于700nm的光可被血红蛋白、肌红蛋白等组织吸收，而大于900nm的红外光可被水等液体吸收，都不易透过组织，在吲哚菁绿的吸收发射光谱范围内的光，不易被组织吸收，穿透力强，更适合体内应用。
吲哚菁绿血管造影技术(indocyanine green angiography，ICGA)是一种血管成像方法，在术中或术后可以用于评估浅表组织的血流情况，提供定量的、实时动态的、可视化的组织血运评估结果。ICGA利用吲哚菁绿近红外光吸收和荧光发射特性，通过近红外光激发血液中的吲哚菁绿，使其产生荧光，用成像器件捕捉并呈现荧光图像，以荧光强度的静态指标或者荧光强度随时间变化(time-density curve，TDC)的动态指标提供组织血流的可视化及定量评估参数，该技术作为一种辅助手段，可以用于评估断指(肢)再植指或皮瓣移植过程中吻合血管的通畅情况及组织瓣的血液灌注状态。ICGA是目前唯一可在术中定量评估组织灌注的新技术，与传统方式相比具有无创、无辐射、无肾毒性、可多次重复注射、可定量等优势，具有巨大的临床应用前景，但就目前研究进展而言仍存在局限性。
从光电容积脉搏波(photoplethysmography，PPG)信号中提取ICG色素浓度谱计算ICG色素在人体内的代谢规律已成功应用到肝脏外科以无创评估肝脏储备功能。该方法利用光电指夹、鼻夹或耳夹从人体健侧(指端、鼻翼、耳垂)检测动脉脉搏波信号，根据朗伯-比尔定律由血红蛋白浓度计算ICG色素代谢曲线。再植指或游离皮瓣因吻合血管通畅程度不同而导致患处ICG代谢规律与健侧ICG代谢规律有所差异，利用这种差异可以提高ICGA评估断指(肢)再植和游离皮瓣术后血运的准确度。
依据外科医生的临床经验来评估断指(肢)再植及游离皮瓣术后血运情况，这种做法具有很强的主观性，外科医生会根据皮瓣的颜色、温度、指压反应和毛细血管回流情况等评估皮瓣的血运，如果单独通过上述观测指标，要想准确评估组织的灌注、预测皮瓣的转归，通常非常困难。依靠现代医学影像学技术评估组织血液灌注对检查者的临床经验要求较高，部分设备比较昂贵、笨重，且具有放射性，只能在术前使用，不能在术中及时对组织血运情况进行评估。目前影响ICGA准确性的因素较多，如温度、光线、ICG浓度、术中血流动力学改变、设备稳定性等，使ICGA结果的假阳性和假阴性仍广泛存在。
发明内容
本发明的目的是要解决现有断肢再植及游离皮瓣术后血运评估依靠医生临床经验判断具有不准确、主观性强、不能实时观察的缺点，提供一种断肢再植及游离皮瓣术后血运评估装置及方法，可实时、定量、准确地评估断指(肢)再植及游离皮瓣术后血运情况。
为达到上述目的，本发明是按照以下技术方案实施的：
一种断肢再植及游离皮瓣术后血运评估装置，包括光电容积脉搏波信号采集装置、荧光图像采集装置和工作站；所述光电容积脉搏波信号采集装置采集健侧携带有吲哚菁绿药物浓度信息的光电容积脉搏波，并将数字量化后的光电容积脉搏波信号发送给工作站；荧光图像采集装置采集患处血液中的吲哚菁绿产生的荧光图像，并将荧光图像发送到工作站；工作站对接收的数字量化后的光电容积脉搏波信号和荧光图像进行处理，计算人体内吲哚菁绿药物正常代谢曲线，并结合人体内吲哚菁绿药物正常代谢曲线计算断肢再植及游离皮瓣术后血运评估参数。
进一步地，所述光电容积脉搏波信号采集装置包括光电指夹和与光电指夹连接的光电容积脉搏波信号数字量化模块；所述光电指夹夹持在再植肢解剖对称指上，光电指夹将透过手指携带有吲哚菁绿药物浓度信息的光信号转化为电信号以便光电容积脉搏波信号数字量化模块处理。
进一步地，所述荧光图像采集装置为为带有近红外光源的荧光成像相机，近红外激发光源激发患处血液中的吲哚菁绿药物产生荧光，荧光成像相机对产生的荧光信号进行成像，所述近红外光源中心波长为740-805nm。
进一步地，所述工作站为移动电脑或台式电脑。
进一步地，所述近红外光源中心波长为760nm。
另外，利用上述断肢再植及游离皮瓣术后血运评估装置可以进行断肢再植及游离皮瓣术后血运评估，包括以下步骤：
将光电指夹夹持在再植肢解剖对称指上，将患处置于带有近红外光源的荧光成像相机的成像视野中，并开启近红外光源，按照医生指示将一定浓度和剂量的吲哚菁绿药物通过外周静脉或中心静脉团注给药进入人体，同时带有近红外光源的荧光成像相机和光电指夹开始采集荧光图像和光电容积脉搏波信号，并将数据发送到工作站，连续采集3-15min后，停止采集，工作站对接收的数据进行处理，计算人体内吲哚菁绿药物正常代谢参数，并结合人体内吲哚菁绿药物正常代谢曲线计算断肢再植及游离皮瓣术后血运评估参数，断肢再植及游离皮瓣术后血运评估参数由根据荧光强度随时间变化曲线计算的参数和吲哚菁绿药物正常代谢曲线计算的参数相比计算得到，即：
相对转运时间rITT＝(ITT(x，y)-ITT)/ITT，if ITT(x，y)<ITT ITT(x，y)＝ITT；
相对达峰时间rTTP＝(TTP(x，y)-TTP)/TTP，if TTP(x，y)<TTP TTP(x，y)＝TTP；
相对平均循环时间rMTT＝(MTT(x，y)-MTT)/MTT，if MTT(x，y)<MTT MTT(x，y)＝MTT；
相对再转运时间rIRTT＝(IRTT(x，y)-IRTT)/IRTT，if IRTT(x，y)<IRTT IRTT(x，y)＝IRTT；
相对灌注速率rPR＝|(PR(x，y)-PR)|/PR，if PR(x，y)>PR PR(x，y)＝PR；
相对血浆消失率rK＝|(K(x，y)-K)|/K，if K(x，y)>K K(x，y)＝K；
其中：将吲哚菁绿药物开始注入时刻到ICG色素到达光电指夹之间的时间t0记为转运时间ITT；
将吲哚菁绿药物开始出现到第一循环曲线浓度最大值之间的时间t0-t1记为达峰时间TTP；
将第一循环曲线浓度最大值之后第5-10s的数据进行指数拟合与横轴形成封闭区间，区间面积重心对应的时刻t2记为平均循环时间MTT；
将吲哚菁绿药物开始注入时刻到吲哚菁绿药物第二循环曲线开始时刻之间的时间t3记为再转运时间IRTT；
将吲哚菁绿药物浓度从10％上升到90％的速率记为灌注速率PR；
将3-6min内的吲哚菁绿药物浓度数据进行线性拟合，直线斜率的绝对值记为血浆消失率K；
(x，y)为当前坐标，每个参数均为二维矩阵，矩阵大小和图像大小相同；
将相对转运时间、相对达峰时间、相对平均循环时间、相对再转运时间、相对灌注速率和相对血浆消失率参数进行归一化处理，即将矩阵元素值映射到0-1之间，进行伪彩色处理，即矩阵元素值越大，颜色越深，矩阵元素值越小，颜色越浅；并将结果显示在工作站的显示器上。
与现有技术相比，本发明具有以下有益效果：
1.本发明以PPG信号中提取的ICG色素浓度谱作为ICG色素在人体内的正常代谢曲线，通过比较ICGA采集的每一像素点处的ICG代谢曲线与人体内正常ICG代谢曲线的差异性，评估断指(肢)再植及游离皮瓣术后血运情况，避免或减小了不利因素对ICGA评估皮瓣血运结果准确性的影响，然后直接生成断指(肢)再植及游离皮瓣术后血运评估参数，并以图像化显示，具有直观性和客观性，避免外科医生依靠临床经验评估断指(肢)再植及游离皮瓣术后血运情况的主观性；
2.本发明公开的一种断肢再植及游离皮瓣术后血运评估装置及方法，使用吲哚菁绿作为血管指示剂，具有无肾毒性，无辐射性，半衰期短，可以在术前、术中、术后多次使用。
附图说明
图1为本发明的断肢再植及游离皮瓣术后血运评估装置的布置示意图。
图2为人体内ICG色素代谢曲线。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步的详细说明。此处所描述的具体实施例仅用于解释本发明，并不用于限定发明。
实施例1
如图1所示，本实施例提供了首先提供了一种断肢再植及游离皮瓣术后血运评估装置，包括光电容积脉搏波信号采集装置、荧光图像采集装置和工作站5；光电容积脉搏波信号采集装置包括光电指夹2(可以使用市购的DDG-3300K探头)和与光电指夹连接的光电容积脉搏波采集盒1(可以直接在市面上购买)；所述光电指夹2夹持在再植肢解剖对称指上(如图1所示，左手中指为再植指，光电指夹则夹持在右手中指)，光电指夹2将透过手指携带有吲哚菁绿药物浓度信息的光信号转化为电信号传输至光电容积脉搏波采集盒1，光电容积脉搏波采集盒1将数字量化后的光电容积脉搏波信号发送给工作站；荧光图像采集装置为带有近红外光源的荧光成像相机4，近红外激发光源激发患处3血液中的吲哚菁绿药物产生荧光，荧光成像相机4对产生的荧光信号进行成像，所述近红外光源中心波长为740-805nm(本实施例中优选760nm)，荧光成像相机4的分辨率和成像视野相关，要求在合适的工作距离下可以分辨0.4mm直径以上的毛细血管，本实施例中荧光成像相机4分辨率为500万像素(2448x2048)。荧光成像相机4成像帧率要求不低于5帧/秒，本实施例中优选为5帧/秒，评估时将患处域置于荧光成像相机的成像视野中，调整荧光成像相机4以获得合适的成像视野和良好的成像效果，荧光成像相机4采集患处3血液中的吲哚菁绿产生的荧光图像，并将荧光图像发送到工作站5；工作站5为移动电脑或台式电脑，工作站5对接收的数字量化后的光电容积脉搏波信号和荧光图像进行处理，基于朗伯-比尔定律从光电容积脉搏波信号中计算吲哚菁绿药物正常代谢曲线，并结合人体内吲哚菁绿药物正常代谢曲线计算断肢再植及游离皮瓣术后血运评估参数。
需要说明的是：朗伯-比尔定律描述的是当用波长为λ的光照射某一溶液时，入射光强I0与出射光强I之比定义为该溶液中的溶质对该波长的光的吸光度。公式为：
Aλ为吸光度，ελ为介质的吸光系数，只与溶质种类和入射光波长有关，与溶液浓度、厚度无关，C为溶质浓度，D为光程。
当溶液中含有多种溶质时，朗伯-比尔定律符合叠加原理：
脉搏色素浓度谱基本思路是：血液可视为一种血红蛋白与血浆混合的溶液，血浆充当溶剂，血红蛋白充当溶质。当吲哚菁绿作为另一种溶质注入血液时，血浆中就含有了两种溶剂。因为血红蛋白和ICG色素对不同波长的近红外光有不同的吸光系数，就会引起透射光强的变化。根据血红蛋白和ICG色素吸光系数的差异，用不同波长的近红外光照射血液，检测透射光强就可以推出物质的浓度比。以人体指端建立朗伯-比尔定律模型：
用单波长的光照射人体指端时，出射光强同时受搏动的动脉血和固定厚度的组织的影响。前者随动脉脉搏的变化而改变，记为交流成分AC；后者的透射光强保持恒定，记为直流成分DC。由于交流成分占直流成分的比例远小于1，若动脉搏动前透射光强为I，搏动时透射光强为I-ΔI，那么该波长下吸光度为：
在805nm和940nm两种波长光照射下，吸光度之比为：
则ICG色素浓度谱计算公式为：
其中CHb为血红蛋白浓度，由血常规得到。
实施例2
本实施例利用上述断肢再植及游离皮瓣术后血运评估装置，进而提供了一种断肢再植及游离皮瓣术后血运评估方法，包括以下步骤：
将光电指夹夹持在再植肢解剖对称指上，将患处置于带有近红外光源的荧光成像相机的成像视野中，并开启近红外光源，按照医生指示将一定浓度和剂量的吲哚菁绿药物通过外周静脉或中心静脉团注给药进入人体，同时带有近红外光源的荧光成像相机和光电指夹开始采集荧光图像和光电容积脉搏波信号，并将数据发送到工作站，连续采集3-15min后，停止采集，工作站对接收的数据进行处理，计算人体内吲哚菁绿药物正常代谢曲线，如图2所示。
血液流动在人体内是一个不断循环的过程，因此注入血液中的ICG色素在未完全均匀地分布在血液之前，从肢端检测的ICG色素浓度会有数个随血液循环变化的峰值过程，只不过峰值逐渐减小。通常ICG色素代谢曲线包括色素稀释的第一循环曲线和浓度按指数规律衰减的第二循环曲线。
根据图2，将吲哚菁绿药物开始注入时刻到ICG色素到达光电指夹之间的时间t0记为转运时间ITT；
将吲哚菁绿药物开始出现到第一循环曲线浓度最大值之间的时间t0-t1记为达峰时间TTP；
将第一循环曲线浓度最大值之后第5-10s的数据进行指数拟合与横轴形成封闭区间，区间面积重心对应的时刻t2记为平均循环时间MTT；
将吲哚菁绿药物开始注入时刻到吲哚菁绿药物第二循环曲线开始时刻之间的时间t3记为再转运时间IRTT；
将吲哚菁绿药物浓度从10％上升到90％的速率记为灌注速率PR；
将3-6min内的吲哚菁绿药物浓度数据进行线性拟合，直线斜率的绝对值记为血浆消失率K。
然后，结合人体内吲哚菁绿药物正常代谢曲线计算断肢再植及游离皮瓣术后血运评估参数，断肢再植及游离皮瓣术后血运评估参数由根据荧光强度随时间变化曲线计算的参数和吲哚菁绿药物正常代谢曲线计算的参数相比计算得到，即：
相对转运时间rITT＝(ITT(x，y)-ITT)/ITT，if ITT(x，y)<ITT ITT(x，y)＝ITT；
相对达峰时间rTTP＝(TTP(x，y)-TTP)/TTP，if TTP(x，y)<TTP TTP(x，y)＝TTP；
相对平均循环时间rMTT＝(MTT(x，y)-MTT)/MTT，if MTT(x，y)<MTT MTT(x，y)＝MTT；
相对再转运时间rIRTT＝(IRTT(x，y)-IRTT)/IRTT，if IRTT(x，y)<IRTT IRTT(x，y)＝IRTT；
相对灌注速率rPR＝|(PR(x，y)-PR)|/PR，if PR(x，y)>PR PR(x，y)＝PR；
相对血浆消失率rK＝|(K(x，y)-K)|/K，if K(x，y)>K K(x，y)＝K；
其中：(x，y)为当前坐标，每个参数均为二维矩阵，矩阵大小和图像大小相同；
相对转运时间、相对达峰时间、相对平均循环时间、相对再转运时间、相对灌注速率和相对血浆消失率均为二维矩阵，且矩阵中某元素值越大，表示该元素坐标所记录的患处一定大小区域的血液灌注状态越差。
将相对转运时间、相对达峰时间、相对平均循环时间、相对再转运时间、相对灌注速率和相对血浆消失率参数进行归一化处理，即将矩阵元素值映射到0-1之间，进行伪彩色处理，即矩阵元素值越大，颜色越深，矩阵元素值越小，颜色越浅；并将结果显示在工作站的显示器上。
通过本实施例的方法，可以为临床医生提供了客观的、实时动态的、可视化的组织血运评估结果，提高了ICGA结果的准确性。
本发明的技术方案不限于上述具体实施例的限制，凡是根据本发明的技术方案做出的技术变形，均落入本发明的保护范围之内。
Claims (6)
1.一种断肢再植及游离皮瓣术后血运评估装置，其特征在于：包括光电容积脉搏波信号采集装置、荧光图像采集装置和工作站；所述光电容积脉搏波信号采集装置采集健侧携带有吲哚菁绿药物浓度信息的光电容积脉搏波，并将数字量化后的光电容积脉搏波信号发送给工作站；荧光图像采集装置采集患处血液中的吲哚菁绿产生的荧光图像，并将荧光图像发送到工作站；工作站对接收的数字量化后的光电容积脉搏波信号和荧光图像进行处理，计算人体内吲哚菁绿药物正常代谢曲线，并结合人体内吲哚菁绿药物正常代谢曲线计算断肢再植及游离皮瓣术后血运评估参数。
2.根据权利要求1所述的断肢再植及游离皮瓣术后血运评估装置，其特征在于：所述光电容积脉搏波信号采集装置包括光电指夹和与光电指夹连接的光电容积脉搏波信号数字量化模块；所述光电指夹夹持在再植肢解剖对称指上，光电指夹将透过手指携带有吲哚菁绿药物浓度信息的光信号转化为电信号以便光电容积脉搏波信号数字量化模块处理。
3.根据权利要求1所述的断肢再植及游离皮瓣术后血运评估装置，其特征在于：所述荧光图像采集装置为为带有近红外光源的荧光成像相机，近红外激发光源激发患处血液中的吲哚菁绿药物产生荧光，荧光成像相机对产生的荧光信号进行成像，所述近红外光源中心波长为740-805nm。
4.根据权利要求1所述的断肢再植及游离皮瓣术后血运评估装置，其特征在于：所述工作站为移动电脑或台式电脑。
5.根据权利要求3所述的断肢再植及游离皮瓣术后血运评估装置，其特征在于：所述近红外光源中心波长为760nm。
6.一种断肢再植及游离皮瓣术后血运评估方法，其特征在于，包括以下步骤：
将光电指夹夹持在再植肢解剖对称指上，将患处置于带有近红外光源的荧光成像相机的成像视野中，并开启近红外光源，按照医生指示将一定浓度和剂量的吲哚菁绿药物通过外周静脉或中心静脉团注给药进入人体，同时带有近红外光源的荧光成像相机和光电指夹开始采集荧光图像和光电容积脉搏波信号，并将数据发送到工作站，连续采集3-15min后，停止采集，工作站对接收的数据进行处理，计算人体内吲哚菁绿药物正常代谢参数，并结合人体内吲哚菁绿药物正常代谢曲线计算断肢再植及游离皮瓣术后血运评估参数，断肢再植及游离皮瓣术后血运评估参数由根据荧光强度随时间变化曲线计算的参数和吲哚菁绿药物正常代谢曲线计算的参数相比计算得到，即：
相对转运时间rITT＝(ITT(x，y)-ITT)/ITT，if ITT(x，y)<ITT ITT(x，y)＝ITT；
相对达峰时间rTTP＝(TTP(x，y)-TTP)/TTP，if TTP(x，y)<TTP TTP(x，y)＝TTP；
相对平均循环时间rMTT＝(MTT(x，y)-MTT)/MTT，if MTT(x，y)<MTT MTT(x，y)＝MTT；
相对再转运时间rIRTT＝(IRTT(x，y)-IRTT)/IRTT，if IRTT(x，y)<IRTT IRTT(x，y)＝IRTT；
相对灌注速率rPR＝|(PR(x，y)-PR)|/PR，if PR(x，y)>PR PR(x，y)＝PR；
相对血浆消失率rK＝|(K(x，y)-K)|/K，if K(x，y)>K K(x，y)＝K；
其中：将吲哚菁绿药物开始注入时刻到ICG色素到达光电指夹之间的时间t0记为转运时间ITT；
将吲哚菁绿药物开始出现到第一循环曲线浓度最大值之间的时间t0-t1记为达峰时间TTP；
将第一循环曲线浓度最大值之后第5-10s的数据进行指数拟合与横轴形成封闭区间，区间面积重心对应的时刻t2记为平均循环时间MTT；
将吲哚菁绿药物开始注入时刻到吲哚菁绿药物第二循环曲线开始时刻之间的时间t3记为再转运时间IRTT；
将吲哚菁绿药物浓度从10％上升到90％的速率记为灌注速率PR；
将3-6min内的吲哚菁绿药物浓度数据进行线性拟合，直线斜率的绝对值记为血浆消失率K；
(x，y)为当前坐标，每个参数均为二维矩阵，矩阵大小和图像大小相同；
将相对转运时间、相对达峰时间、相对平均循环时间、相对再转运时间、相对灌注速率和相对血浆消失率参数进行归一化处理，即将矩阵元素值映射到0-1之间，进行伪彩色处理，即矩阵元素值越大，颜色越深，矩阵元素值越小，颜色越浅；并将结果显示在工作站的显示器上。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111500149.2A CN114190893B (zh) | 2021-12-09 | 2021-12-09 | 一种断肢再植及游离皮瓣术后血运评估装置及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111500149.2A CN114190893B (zh) | 2021-12-09 | 2021-12-09 | 一种断肢再植及游离皮瓣术后血运评估装置及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114190893A true CN114190893A (zh) | 2022-03-18 |
| CN114190893B CN114190893B (zh) | 2024-01-05 |
Family
ID=80651694
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111500149.2A Active CN114190893B (zh) | 2021-12-09 | 2021-12-09 | 一种断肢再植及游离皮瓣术后血运评估装置及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114190893B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117137450A (zh) * | 2023-08-30 | 2023-12-01 | 哈尔滨海鸿基业科技发展有限公司 | 一种基于皮瓣血运评估的皮瓣移植术成像方法和系统 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080249382A1 (en) * | 2007-04-04 | 2008-10-09 | Lg Electronics Inc. | Blood pressure monitoring apparatus and method |
| US20110071403A1 (en) * | 2009-09-21 | 2011-03-24 | Board Of Regents Of The University Of Texas System | Functional near-infrared fluorescence lymphatic mapping for diagnosing, accessing, monitoring and directing therapy of lymphatic disorders |
| US20140371583A1 (en) * | 2013-06-14 | 2014-12-18 | Novadaq Technologies Inc. | Quantification of absolute blood flow in tissue using fluorescence mediated photoplethysmography |
| US20150182137A1 (en) * | 2013-06-14 | 2015-07-02 | Novadaq Technologies Inc. | Quantification of absolute blood flow in tissue using fluorescence-mediated photoplethysmography |
| CN109770946A (zh) * | 2019-03-22 | 2019-05-21 | 苏州大学附属第一医院 | 一种微型探头应用于血管吻合术中评估血流动力学的方法 |
| CN112516311A (zh) * | 2020-12-23 | 2021-03-19 | 四川大学华西医院 | 靶向氧化石墨烯复合物用于治疗耐药骨肉瘤的实验方法 |
| WO2021074097A1 (en) * | 2019-10-14 | 2021-04-22 | Koninklijke Philips N.V. | Perfusion angiography combined with photoplethysmography imaging for peripheral vascular disease assessment |
| CN113261920A (zh) * | 2021-05-10 | 2021-08-17 | 吉林大学第一医院 | 一种肢体组织微循环血液灌注状态监测评估系统 |
-
2021
- 2021-12-09 CN CN202111500149.2A patent/CN114190893B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080249382A1 (en) * | 2007-04-04 | 2008-10-09 | Lg Electronics Inc. | Blood pressure monitoring apparatus and method |
| US20110071403A1 (en) * | 2009-09-21 | 2011-03-24 | Board Of Regents Of The University Of Texas System | Functional near-infrared fluorescence lymphatic mapping for diagnosing, accessing, monitoring and directing therapy of lymphatic disorders |
| US20140371583A1 (en) * | 2013-06-14 | 2014-12-18 | Novadaq Technologies Inc. | Quantification of absolute blood flow in tissue using fluorescence mediated photoplethysmography |
| US20150182137A1 (en) * | 2013-06-14 | 2015-07-02 | Novadaq Technologies Inc. | Quantification of absolute blood flow in tissue using fluorescence-mediated photoplethysmography |
| CN109770946A (zh) * | 2019-03-22 | 2019-05-21 | 苏州大学附属第一医院 | 一种微型探头应用于血管吻合术中评估血流动力学的方法 |
| WO2021074097A1 (en) * | 2019-10-14 | 2021-04-22 | Koninklijke Philips N.V. | Perfusion angiography combined with photoplethysmography imaging for peripheral vascular disease assessment |
| CN112516311A (zh) * | 2020-12-23 | 2021-03-19 | 四川大学华西医院 | 靶向氧化石墨烯复合物用于治疗耐药骨肉瘤的实验方法 |
| CN113261920A (zh) * | 2021-05-10 | 2021-08-17 | 吉林大学第一医院 | 一种肢体组织微循环血液灌注状态监测评估系统 |
Non-Patent Citations (2)
| Title |
|---|
| GUANGDA LIU: ""SVR-EEMD: An Improved EEMD Method Based on Support Vector Regression Extension in PPG Signal Denoising"", 《COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE》 * |
| 刘颂阳等: ""扩展的卡尔曼滤波在近红外光谱提取脑血流信号中的研究"", 《光谱学与光谱分析》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117137450A (zh) * | 2023-08-30 | 2023-12-01 | 哈尔滨海鸿基业科技发展有限公司 | 一种基于皮瓣血运评估的皮瓣移植术成像方法和系统 |
| CN117137450B (zh) * | 2023-08-30 | 2024-05-10 | 哈尔滨海鸿基业科技发展有限公司 | 一种基于皮瓣血运评估的皮瓣移植术成像方法和系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114190893B (zh) | 2024-01-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10463286B2 (en) | Determination of tissue oxygenation in vivo | |
| US11701016B2 (en) | Facilitating assessment of blood flow and tissue perfusion using fluorescence-mediated photoplethysmography | |
| CN103347439B (zh) | 动态光学层析成像系统 | |
| US20100210931A1 (en) | Method for performing qualitative and quantitative analysis of wounds using spatially structured illumination | |
| Khalil et al. | Dynamic diffuse optical tomography imaging of peripheral arterial disease | |
| EP2665411B1 (en) | A method and system for visualization of cardiovascular pulsation waves | |
| WO2014027293A2 (en) | Estimations of equivalent inner diameter of arterioles | |
| KR20190133026A (ko) | 당뇨병성 혈액순환 합병증을 평가하기 위한 시스템 및 방법 | |
| CN109924987A (zh) | 反射式血氧仪的标定方法、系统和可读存储介质 | |
| CN109009064A (zh) | 一种颅内压无创监测装置及方法 | |
| US20200294228A1 (en) | Non-Contact Multispectral Imaging for Blood Oxygen Level and Perfusion Measurement and Related Systems and Computer Program Products | |
| WO2022237171A1 (zh) | 一种基于ccd的双波长红外血氧检测系统 | |
| CN114190893B (zh) | 一种断肢再植及游离皮瓣术后血运评估装置及方法 | |
| Wang et al. | Non-invasive measurement of haemoglobin based on dynamic spectrum method | |
| CN104856691A (zh) | 一种基于眼底照相机的视网膜血氧饱和度测量的方法 | |
| WO2009008858A1 (en) | System for repetitive measurements of cardiac output in freely moving individuals | |
| Bhat et al. | The biophysical parameter measurements from PPG signal | |
| Qiu et al. | ICG clearance test based on photoacoustic imaging for assessment of human liver function reserve: An initial clinical study | |
| CN212913171U (zh) | 一种基于吲哚菁绿的荧光照相机 | |
| CN112107305A (zh) | 一种有效肝脏血流量的检测方法及检测装置、存储介质 | |
| KR101032479B1 (ko) | 수학적 모델에 기반한 말초조직 관류 정도 측정장치 및 측정방법 | |
| CN215078399U (zh) | 一种肢体组织微循环血液灌注状态监测评估系统 | |
| Zhang et al. | Quantitative evaluation of human's hepatic functional reserve by indicator pharmacokinetics method | |
| Esenaliev et al. | Noninvasive optoacoustic system for rapid diagnostics and management of circulatory shock | |
| Cai et al. | Quantitative evaluation of graded hindlimb ischemia based on pharmacokinetic modelling and hemodynamic analysis of indocyanine green |
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 |