CN110547785B - 微循环检测系统及检测方法 - Google Patents
微循环检测系统及检测方法 Download PDFInfo
- Publication number
- CN110547785B CN110547785B CN201910477770.8A CN201910477770A CN110547785B CN 110547785 B CN110547785 B CN 110547785B CN 201910477770 A CN201910477770 A CN 201910477770A CN 110547785 B CN110547785 B CN 110547785B
- Authority
- CN
- China
- Prior art keywords
- detection system
- array
- microcirculation
- processing unit
- skin
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 142
- 230000004089 microcirculation Effects 0.000 title claims abstract description 128
- 238000009826 distribution Methods 0.000 claims abstract description 114
- 238000012545 processing Methods 0.000 claims abstract description 99
- 230000008859 change Effects 0.000 claims abstract description 72
- 238000010792 warming Methods 0.000 claims abstract description 41
- 238000010438 heat treatment Methods 0.000 claims abstract description 37
- 230000010355 oscillation Effects 0.000 claims description 21
- 230000003595 spectral effect Effects 0.000 claims description 13
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 210000003491 skin Anatomy 0.000 description 73
- 238000013186 photoplethysmography Methods 0.000 description 59
- 210000004204 blood vessel Anatomy 0.000 description 37
- 210000001519 tissue Anatomy 0.000 description 33
- 230000001121 heart beat frequency Effects 0.000 description 31
- 238000010586 diagram Methods 0.000 description 27
- 210000004207 dermis Anatomy 0.000 description 16
- 230000033001 locomotion Effects 0.000 description 16
- 238000001228 spectrum Methods 0.000 description 16
- 230000000694 effects Effects 0.000 description 14
- 238000000034 method Methods 0.000 description 13
- 238000012544 monitoring process Methods 0.000 description 13
- 238000004364 calculation method Methods 0.000 description 12
- 230000008602 contraction Effects 0.000 description 12
- 239000008280 blood Substances 0.000 description 11
- 210000004369 blood Anatomy 0.000 description 11
- 238000001816 cooling Methods 0.000 description 11
- 238000006243 chemical reaction Methods 0.000 description 8
- 230000017531 blood circulation Effects 0.000 description 7
- 239000003814 drug Substances 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 210000005036 nerve Anatomy 0.000 description 6
- 230000002889 sympathetic effect Effects 0.000 description 6
- 230000002792 vascular Effects 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 5
- 230000035790 physiological processes and functions Effects 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 230000007423 decrease Effects 0.000 description 4
- 206010033675 panniculitis Diseases 0.000 description 4
- 230000035939 shock Effects 0.000 description 4
- 210000004304 subcutaneous tissue Anatomy 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000004087 circulation Effects 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 3
- 206010012601 diabetes mellitus Diseases 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 230000004217 heart function Effects 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 208000033808 peripheral neuropathy Diseases 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 201000001320 Atherosclerosis Diseases 0.000 description 2
- 206010053615 Thermal burn Diseases 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 230000002051 biphasic effect Effects 0.000 description 2
- 230000036770 blood supply Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 238000001647 drug administration Methods 0.000 description 2
- 210000002615 epidermis Anatomy 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000035479 physiological effects, processes and functions Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 210000005070 sphincter Anatomy 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000013271 transdermal drug delivery Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 206010019332 Heat exhaustion Diseases 0.000 description 1
- 206010019345 Heat stroke Diseases 0.000 description 1
- 206010021138 Hypovolaemic shock Diseases 0.000 description 1
- 208000010718 Multiple Organ Failure Diseases 0.000 description 1
- 206010028851 Necrosis Diseases 0.000 description 1
- 206010031264 Osteonecrosis Diseases 0.000 description 1
- 208000018262 Peripheral vascular disease Diseases 0.000 description 1
- 208000034189 Sclerosis Diseases 0.000 description 1
- 206010040030 Sensory loss Diseases 0.000 description 1
- 206010065769 Soft tissue necrosis Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 210000002565 arteriole Anatomy 0.000 description 1
- 208000011775 arteriosclerosis disease Diseases 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000008081 blood perfusion Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000005779 cell damage Effects 0.000 description 1
- 230000030833 cell death Effects 0.000 description 1
- 208000037887 cell injury Diseases 0.000 description 1
- 230000019522 cellular metabolic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 210000001035 gastrointestinal tract Anatomy 0.000 description 1
- 230000017525 heat dissipation Effects 0.000 description 1
- 238000002639 hyperbaric oxygen therapy Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 208000028867 ischemia Diseases 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003340 mental effect Effects 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 239000002207 metabolite Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 208000018731 motor weakness Diseases 0.000 description 1
- 208000029744 multiple organ dysfunction syndrome Diseases 0.000 description 1
- 230000017074 necrotic cell death Effects 0.000 description 1
- 230000035764 nutrition Effects 0.000 description 1
- 235000016709 nutrition Nutrition 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 210000000578 peripheral nerve Anatomy 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000144 pharmacologic effect Effects 0.000 description 1
- 230000037081 physical activity Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 206010040560 shock Diseases 0.000 description 1
- 230000008326 skin blood flow Effects 0.000 description 1
- 210000002460 smooth muscle Anatomy 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
- 230000007556 vascular defect Effects 0.000 description 1
- 210000001835 viscera 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/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/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/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- 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/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/746—Alarms related to a physiological condition, e.g. details of setting alarm thresholds or avoiding false alarms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Hematology (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
一种微循环检测系统,该微循环检测系统包含加热装置、阵列式传感器以及处理单元。所述加热装置用于对皮肤区域加温。所述阵列式传感器用于检测所述皮肤区域的出射光并在加温期间内的不同时间点分别输出多个亮度变化信号。所述处理单元用于根据所述多个亮度变化信号计算阵列能量分布在所述加温期间内的变化并据以判断微循环状态。
Description
技术领域
本发明说明有关一种微循环检测,更特别有关一种检测微循环随皮肤表面温度上升而变化的检测系统及其检测方法。
背景技术
目前,可携式电子装置(portable electronic device)及穿戴式电子装置(wearable electronic device)已成为生活中不可或缺的电子产品,其功能也随着人们生活型态的改变而不断地演进。
同时,身体健康成为现代忙碌的生活中人人所关心的问题,因此生理检测功能也开始逐渐被应用至可携式电子装置及穿戴式电子装置上,以符合用户的需求。
发明内容
有鉴于此,本发明说明提出一种阵列式生理检测系统及其运作方法,其至少可检测并记录用户的三个维度以上的生理特征。
本发明说明提供一种阵列式生理检测系统及其运作方法,其通过多个感测像素分别检测不同被检测组织区域的生理特征,以产生三维生理特征分布。
本发明说明提供一种微循环检测系统,包含光源、加热装置、感光阵列以及处理单元。所述光源用于发光照明皮肤区域。所述加热装置用于对所述皮肤区域加温。所述感光阵列用于检测所述皮肤区域的出射光并输出多个PPG信号。所述处理单元用于将所述多个PPG信号转换为阵列能量分布,辨识所述阵列能量分布中的环状式样及计算加温期间内的所述环状式样的振荡的频率变化。
本发明说明还提供一种微循环检测系统,包含光源、加热装置、计时器、感光阵列以及处理单元。所述光源用于发光照明皮肤区域。所述加热装置用于对所述皮肤区域加温。所述计时器用于计时所述加热装置的加温期间。所述感光阵列用于检测所述皮肤区域的出射光并输出多个PPG信号。所述处理单元用于将所述多个PPG信号转换为阵列能量分布,辨识第一时间点的所述阵列能量分布中的第一环状式样及第二时间点的所述阵列能量分布中的第二环状式样以及判断所述第一环状式样及所述第二环状式样的往复频率的频率变化是否在所述加温期间的第一时间区间及第二时间区间分别具有峰值。
本发明说明还提供一种微循环检测系统的检测方法,该微循环检测系统包含光源、加热装置、感光阵列以及处理单元。所述检测方法包含下列步骤:以所述光源照明皮肤区域;以所述加热装置对所述皮肤区域加温;以所述感光阵列检测所述皮肤区域的出射光并输出多个PPG信号;以所述处理单元将所述多个PPG信号转换为阵列能量分布;以及以所述处理单元辨识所述阵列能量分布中预定位置的能量振荡在加温期间的频率变化在所述加温期间的两个预定时间区间中的两个峰值。
本发明说明实施例的阵列式生理检测系统及其运作方法中,还可建立表示生理特征分布的三维能量分布随时间的三维能量变化,以形成四维生理检测系统。
为了让本发明说明的上述和其他目的、特征和优点能更明显,下文将配合所附图示,详细说明如下。此外,于本发明说明中,相同的构件以相同的符号表示,于此先述明。
附图说明
图1为本发明说明实施例的生理检测系统获得浅层微循环血管涨缩变化的流程图;
图2A为本发明说明实施例的生理检测系统获取的图像帧及其观测窗口的示意图;
图2B为本发明说明实施例的生理检测系统获取的多个图像帧的亮度变化的示意图;
图2C为本发明说明实施例的生理检测系统求得的微循环血管涨缩变化信号的频谱图;
图2D为本发明说明实施例的生理检测系统求得的目前心跳频率下多个像素区域的能量分布的示意图;
图3A为本发明说明实施例的生理检测系统所检测的变化量变化的示意图;
图3B为本发明说明实施例的生理检测系统所检测的平均值变化的示意图;
图4为本发明说明实施例的生理检测系统的示意图;
图5为本发明说明实施例的生理检测系统的运作方法的流程图;
图6A及6B为本发明说明实施例的生理检测所检测的525纳米光的三维能量分布的示意图;
图7A及7B为本发明说明实施例的生理检测所检测的880纳米光的三维能量分布的示意图;
图8为本发明说明另一实施例的生理检测装置的示意图;
图9为本发明说明再一实施例的生理检测系统的方框图;
图10A为本发明说明实施例的微循环检测系统的方框图;
图10B为本发明说明实施例的微循环检测系统的操作示意图;
图11为本发明说明实施例的微循环检测系统的加温示意图;
图12为本发明说明实施例的微循环检测系统检测的频率变化的示意图;
图13为本发明说明实施例的微循环检测系统的检测方法的流程图。
具体实施方式
以下说明内容包含本发明说明的实施方式,以便理解本发明说明如何应用于实际状况。须注意的是,在以下图式中,与本发明说明技术无关的部分已省略,同时为彰显组件之间的关系,图式里各组件之间的比例与真实的组件之间的比例并不一定相同。
请参照图1所示,其为本发明说明一实施例的阵列式生理检测系统获得浅层微循环血管涨缩变化的流程图。阵列式生理检测系统用于通过皮肤表面检测身体组织的浅层微循环血管涨缩变化的三维能量分布,以帮助使用者对自身健康的监测。此外,本发明说明实施例的阵列式生理检测系统可应用于可携式电子装置或穿戴式电子装置以实现随身式生理监测装置,并适合长期进行自我监测。例如,可长期监测所述三维能量分布随时间的三维能量变化,也即微循环信息随时间的变化。藉此,所得到的监测数据可搭配于医疗机构进行的短期健康检查的检测结果,以得到高信赖度的生理状态信息。
首先,阵列式生理检测系统读取多个像素区域输出的多个浅层微循环血管涨缩变化信号,如步骤101所示,例如光体积变化描述波形(PPG)信号。为了获取多个像素的浅层微循环血管涨缩变化信号,生理检测系统需要取得真皮层(dermis)的血管涨缩变化信号以表示微循环数据。举例而言,可以使用光学检测方式来达成,藉由使用特定波长的光线,使得光线能够穿透表皮层(epidermis)但却不会穿过真皮层,然后利用感光阵列来检测皮肤区域内的微循环血管涨缩变化信号;其中,感光阵列包含多个感光像素,每一个感光像素都可以产生一个微循环血管涨缩变化信号,其各种统计数值可以提供给后续应用。
举例来说,可以使用波长525纳米(nm)的光线,皮肤穿透深度低于1毫米(mm)。在不同的身体部位,可以使用不同波长的光线,以探知在真皮层活动的微循环血管变化状态。由于真皮层的深度大约介于1-3毫米,因此光线波长优选选择为无法穿透3毫米深度,例如300-940纳米。
接着,生理检测系统根据所述微循环血管涨缩变化信号建立三维能量分布,如步骤102所示;其中,所述三维能量分布指频谱能量分布。此步骤中,由于各像素所检测到的微循环血管涨缩变化信号的能量包含各种不同的频率,因此可以选定其中一个频率来进行数据分析。一实施例中,可以先利用所述等微循环血管涨缩变化信号估算出目前心跳,再以此为基础,搜集每一个像素在此心跳频率之下的振幅信号的振幅变化数值,来代表浅层微循环血管的涨缩变化。
由于浅层微循环血管的涨缩变化跟随着心脏跳动而产生,因此像素的振幅变化在心跳频率或者其倍数频率会比在其他的频率上更为明显,以利进行后续的分析。
接着,根据所述能量分布的特性参数判断生理状态,如步骤103所示。在这个步骤中,可以根据能量分布的各种特性,例如振幅变化、平均数值、心跳频率等特性来估测生理状态。在其后说明里会描述如何根据特性参数来判断生理特征。
接着,可提供生理状态警示给使用者,如步骤104所示,以便使用者可以据此调整作息与活动内容。
图2A-2D为本发明说明实施例的阵列式生理检测系统获得浅层微循环血管涨缩变化的示意图。以光学式生理检测系统为例,图2A为所获取的图像帧及其观测窗口(WOI)的示意图;其中,所述观测窗口WOI的尺寸及位置可调整。图2B为多个图像帧(例如显示为6秒内)或所述多个图像帧的观测窗口的亮度变化的示意图;其中,所述亮度变化反应微循环血管的涨缩变化。图2C为浅层微循环血管涨缩变化信号的频谱图,其将图2B的亮度变化(即浅层微循环血管涨缩变化信号)转换至频率所得,并显示有目前心跳频率。图2D为目前心跳频率下,多个像素区域分别获得的能量值所组合成的阵列式变化示意图,即振幅分布;其中,柱状高度表示相对所述目前心跳频率的频谱能量。通过图2D可以明了每一个像素区域所获得的检测结果(即能量值)会产生变化的,此变化状态便可以代表生理特征的变化(例如真皮层中微血管的分布及运作等),详述如下。必须说明的是,图2D中每个像素的幅度是一个像素的能量值或多个像素的平均能量值。
本发明说明可利用不同微循环状态,来估算运动状态。例如可将不同微循环状态分为四个状态,分别是运动前状态I、热身完毕状态II、运动中状态III、运动后冷却状态IV,举例说明如下:
请同时参照图3A及3B,图3A为本发明说明实施例的生理检测系统所检测的变化量变化的示意图;图3B为本发明说明实施例的生理检测系统所检测的平均值变化的示意图。
当使用者处于运动前状态I时,振幅信号的振幅变化(A)不高,但振幅信号的平均数值(B)很高。
当使用者处于热身完毕状态II时,振幅信号的振幅变化(A)逐渐升高,但振幅信号的平均数值(B)开始下降。当振幅信号的平均数值(B)低于热身平均阀值(例如THa1)时,代表热身完毕;或者,当振幅信号的平均数值(B)低于热身阀值(例如THa1)且振幅信号的振幅变化(A)高于热身变化阀值(例如THv1)时,代表热身完毕。
当使用者处于运动中状态III时,振幅信号的振幅变化(A)维持在高数值,但平均数值(B)降至一个相对低的水平。当振幅信号的平均数值(B)低于运动平均阀值(例如THa2)时,代表处于运动状态当中;或者,当振幅信号的平均数值(B)低于运动阀值(例如THa2)且振幅信号的振幅变化(A)高于运动变化阀值(例如THv2)时,代表处于运动状态当中。
当使用者处于运动后冷却状态IV时,振幅信号的振幅变化(A)逐渐降低,振幅信号的平均数值(B)开始升高,当平均数值(B)恢复到超过冷却平均阀值(例如THa3)时,代表冷却完毕;或者,当平均数值(B)恢复到超过冷却平均阀值(例如THa3)且振幅信号的振幅变化(A)恢复到低于冷却平均阀值(例如THv3)时,代表冷却完毕。
必须说明的是,虽然图3A-3B中显示四个微循环状态、三个变化阀值及三个平均阀值,其仅用于说明而非用于限定本发明说明,微循环状态、变化阀值及平均阀值的数量及数值根据不同应用而定。
请参照图4所示,其为本发明说明一实施例的阵列式生理检测系统400的示意图。阵列式生理检测系统400用于检测皮肤微循环的变化,并包含光源41、感光阵列43以及处理单元45。
所述光源41可为同调光源、非同调光源或部分同调光源,例如发光二极管、雷射二极管等。所述光源41用于提供光线L照射皮肤区域,且所述光线穿透至所述皮肤区域的真皮层。必须说明的是,本发明说明的阵列式生理检测系统400仅用于检测真皮层内的微循环血管变化状态而不检测真皮层之下的皮下组织中其他组织状态,故选择适当光波长方能达成所述功效。因此,所述光源41的波长选择为不穿透至所述皮肤区域的真皮层的皮下组织,例如所述光源波长选择介于300-940纳米。
其他实施例中,光源模块可具有多个光源,例如525纳米、880纳米及606纳米等不同波长的光源,以得到从人体反射光及散射光的不同结果。例如,当使用525纳米的短波长光源时,三维能量分布的结果会相应施加于人体的物理压力而呈现环状式样(arc-likepattern)。该环状式样可用于评断是否系统已被使用者适当的穿戴。
使用300至940纳米中的短波长光源可造成人体更显著的吸收变化,其推测为造成如图6A-6B的环状式样的一个原因。使用300-940纳米中的长波长光源可在高压状态下呈现三维能量分布的检测结果,故系统可在不同情况(紧密穿戴或宽松穿戴)下切换至不同光源以得到良好的三维能量分布。例如,图7A及图7B为使用长波长光照射时的三维能量分布的示意图,其中在如同图6A及6B所使用的感光阵列的尺寸下没有观测到清楚的环状式样。
所述感光阵列43优选为主动式图像感测阵列,例如CMOS图像传感器,如此可根据取样结果实时选择观测窗口WOI(如图2A所示)的尺寸及位置,例如根据图像质量或亮度分布等决定所述观测窗口WOI,而所述感光阵列43中所述观测窗口WOI以外的像素数据可不被所述感光阵列43输出。所述感光阵列43包含多个感光像素,所述多个感光像素中每一者用于连续检测经过所述皮肤区域的真皮层的出射光,以输出多个亮度信号作为PPG信号(即浅层微循环血管涨缩变化信号),如图2B所示亮度变化信号。某些实施例中,所述多个亮度信号为数字信号,也即所述感光阵列43可包含模拟数字转换器(ADC)用于进行模拟-数字转换。
所述处理单元45用于将相对所述多个感光像素的所述多个亮度变化信号(即PPG信号)转换为频域数据(如图2C所示),藉以形成微循环血管涨缩变化的三维能量分布(如图2D所示)。所述处理单元45还计算所述多个频域数据的变化量及平均值,以根据所述变化量的变化(例如图3A)及所述平均值的变化(例如图3B)判断不同微循环状态。所述处理单元45例如可为数字信号处理器(DSP)、中央处理器(CPU)、微控制器(MCU)等可用于计算感测阵列所输出的数据的装置,并无特定限制。
所述处理单元45可利用软件、硬件、韧体或其组合完成上述运算。例如,所述处理单元45可包含频域转换模块451、心跳计算模块452、变化计算模块453、平均计算模块454、比较单元455以及储存单元456。可以了解的是,图4以不同组件说明不同运算功能,然而由于该些组件均位于所述处理单元45内,故该些组件所进行的运算即为所述处理单元45所执行的运算。此外,所述处理单元45内还可包含其他运算功能,例如滤波、放大等,本发明说明中省略了其他不直接相关功能的说明。
例如,所述感光阵列43的每一感光像素均输出随时间变化的亮度信号作为PPG信号(如图2B所示),所述处理单元45用于根据多个所述PPG信号计算心跳频率。
一实施例中,所述频域转换模块451将相对每一感光像素的PPG信号(如图2B所示)转换至频域以产生频域数据(如图2C所示),所述心跳计算模块452根据相对所述多个感光像素中每一者的所述频域数据计算估测心跳频率,并将相对所述多个感光像素的多个所述估测心跳频率之中统计量最高的估测心跳频率作为所述心跳频率。也即,相对所述多个感光像素中每一者均可算出估测心跳频率,当计算出某一估测心跳频率的感光像素的数目最多时,则以所述估测心跳频率作为心跳频率。藉此,可降低噪声干扰所导致的错误而增加计算正确度。
另一实施例中,所述处理单元45计算所述感光阵列43输出的每一图像帧中(或所述观测窗口内)所有或部分的所述多个感光像素的所述多个亮度信号的亮度和,并根据相对多个图像帧的多个亮度和计算所述心跳频率。也即,本实施例中,所述处理单元45相对每一图像帧求得亮度和,多个图像帧则可求得亮度和变化以作为PPG信号,如图2B所示。本实施例中,所述心跳计算模块452可直接在时域(time domain)计算心跳频率,例如计算图2B中时间间隔THR的倒数;或者,所述频域转换模块451先将亮度和变化转换至频域以产生频域数据,如图2C所示,接着所述心跳计算模块452根据所述频域数据计算心跳频率,如图2C中频谱能量值最高者。换句话说,本发明说明中,图2B可表示单一感光像素输出的亮度变化或多个图像帧输出的亮度和变化;图2C可表示单一感光像素输出的亮度变化的频域数据或多个图像帧输出的亮度和变化的频域数据,端视不同应用而定。本发明说明中,所述时域-频率转换可以适当频域转换算法进行,例如快速傅立叶,并无特定限制。
所述心跳频率决定后,所述变化计算模块453可相对所述多个感光像素中每一者产生位于所述心跳频率的频谱能量值,以形成三维能量分布或能量集合,如图2D所示。所述变化计算模块453计算所述三维能量分布或能量集合的能量变化量以作为振幅变化数值,例如计算相邻像素能量间的差的总和、每一像素能量与平均能量的差的总和、能量集合的变化数(variance)等,只要能计算出表示所述三维能量分布或能量集合各能量成份(components)间的变化即可,并无特别限制。本实施例中,所述变化量为所述心跳频率的频谱能量变化量。
所述心跳频率决定后,所述平均计算模块454可相对所述多个感光像素中每一者产生位于所述心跳频率的频谱能量值,以形成三维能量分布或能量集合,如图2D所示。所述平均计算模块454计算所述三维能量分布或能量集合的平均值以作为振幅平均值。本实施例中,所述平均值为所述心跳频率的频谱能量平均值。
图2D中的三维能量分布可用以计算并产生代表位置以供后续评断。例如,超过阈值的分布值可用以计算质心位置、重心位置或使用的分布值的中心位置。所述代表位置是根据用户的生理状态而改变,其中所述生理状态是与使用者的肢体动作或身心变化相关。
必须说明的是,虽然上述实施例中以所述变化计算模块453及所述平均计算模块454产生所述三维能量分布或能量集合来说明,然而本发明说明并不以此为限,所述三维能量分布或能量集合可由所述处理单元45所包含的其他模块来运算,例如所述频域转换模块451、所述心跳计算模块452等,并无特定限制。
某些实施例中,除了所述变化量的变化及所述平均值的变化外,所述处理单元45还可同时配合心跳频率判断所述不同微循环状态。也即,本发明说明中,所述处理单元45可根据变化量、平均值及心跳频率的变化的不同组合来判断不同微循环状态(microcirculation state),例如前述的运动前状态、热身完毕状态、运动状态以及运动后冷却状态等,然而并不以此为限。
所述比较单元455可用于比较所述变化量与至少一个变化阈值(例如图3A的THv1-THv3)以判断不同微循环状态。所述比较单元455可用于比较所述平均值与至少一个平均阈值(例如图3B的THa1-THa3)以判断不同微循环状态。所述比较单元455可用于比较心跳频率与至少一个心跳阈值以判断不同微循环状态。所述阈值可事先储存于所述储存单元456中;其中,所述储存单元456例如可为已知内存,并无特定限制。
请参照图5所示,其为本发明说明实施例的阵列式生理检测系统的运作方法,用于通过多个感光像素检测皮肤微循环的变化。所述运作方法包含下列步骤:利用光源提供光线照射皮肤区域并穿透至所述皮肤区域的真皮层(步骤S51);利用每一感光像素连续检测经过所述皮肤区域的所述真皮层的出射光以分别输出亮度变化信号(步骤S52);转换相对每一感光像素的所述亮度变化信号为频域数据(步骤S53);计算相对所述多个感光像素的所述多个频域数据的变化量及/或平均值(步骤S54);根据所述变化量及/或所述平均值的变化判断微循环状态(步骤S55)。
请同时参照图2A-2D、3A-3B及4-5所示,本运作方法的实施方式说明如下。
步骤S51:所述光源41提供光线L照射皮肤区域并穿透至所述皮肤区域的真皮层。如前所述,所述光线L的波长选择为不穿透至皮下组织,因此多个感光像素仅检测微循环血管的数据而不检测皮下组织的数据。
步骤S52:所述感光阵列43的每一感光像素连续检测经过所述皮肤区域的真皮层的出射光以分别输出亮度变化信号,如图2B的PPG信号。因此,所述感光阵列43输出的PPG信号的数目与有效像素数目相同。
步骤S53:所述处理单元45则转换相对每一感光像素的所述亮度变化信号为频域数据,如图2C所示。因此,所产生的频域数据的数目同样与有效像素数目相同。
步骤S54:所述处理单元45接着计算相对所述多个感光像素的所述多个频域数据的变化量及/或平均值。如前所述,由于所述多个频域数据位于心跳频率或其倍频的特征较为明显。因此,计算所述变化量及/或所述平均值前,所述处理单元45先根据所述多个亮度变化信号计算心跳频率,如前述可于时域直接计算或于频域以不同方式计算,并产生位于所述心跳频率的三维能量分布或能量集合,如图2D所示。接着,所述处理单元45可根据所述三维能量分布或能量集合计算所述心跳频率的频谱能量平均值及/或所述心跳频率的频谱能量变化量;其中,所述变化量的计算已说明如前,故于此不再赘述。
步骤S55:所述处理单元45可根据所述变化量及/或所述平均值随时间的变化判断微循环状态;其中,判断的方式例如为比较所述变化量与至少一个变化阈值及/或比较所述平均值与至少一个平均阈值,例如图3A-3B所示者。
如前所述,某些实施例中,所述处理单元45还可搭配心跳频率随时间的变化来判断所述微循环状态。
最后,所述处理单元45可通过不同方式提示使用者所判断的微循环状态,例如通过图像、声音等,并无特定限制。
综上所述,本发明说明并不依据最大心跳率的百分比以及使用者自觉判断来判断运动状态,而是依据与血液流动分布相关的皮肤浅层微循环血管变化来达成;其中,所述皮肤浅层微循环血管变化的数据例如以感光阵列的多个感光像素连续检测经过真皮层的出射光所输出的多个亮度信号呈现,所述亮度信号的变化可称为光体积变化描述波形(PPG)信号。
如前例所述,当所述平均值降低为小于热身平均阈值(例如THa1)及/或所述变化量增加为大于热身变化阈值(例如THv1)时,所述处理单元45判断进入所述热身完毕状态。当所述平均值进而降低为小于运动平均阈值(例如THa2)及/或所述变化量进而增加为大于运动变化阈值(例如THv2)时,所述处理单元45判断进入所述运动状态。当所述平均值从所述运动平均阈值以下增加为大于冷却平均阈值(例如THa3)及/或从所述运动变化阈值以上减少为小于冷却变化阈值(例如THv3)时,所述处理单元45则判断进入所述运动后冷却状态。此外,不同状态与阈值的比较根据不同应用而定。
另一实施例中,本发明说明的生理检测装置,例如图4的400,还用于确认所述生理检测装置是否与肤面有良好的接触,以使得所述生理检测装置正常运作。已知的是,肤面与生理检测装置之间的相对运动会降低图像质量。因此,确认接触状态甚为重要。
例如参照图4,本实施例的生理检测装置400包含光源模块41、感光阵列43及处理单元45。此外,本实施例的生理检测装置400还包含显示器47(如图8所示)用于显示生理检测装置的检测结果,例如,显示警示信息、指示信息等。
本实施例中,光源模块41用于发出不同波长的光至皮肤下方的组织区域以检测所述组织区域的不同深度。如前所述,短波长可用于确认是否已经形成良好穿戴。例如,光源模块41发出第一波长的光,例如介于500纳米和550纳米之间,至欲量测的组织区域。
感光阵列43用于检测所述组织区域的出射光并输出多个PPG信号,每一个PPG信号如图2B所示。如前所述,本发明配置成一个像素输出一个PPG信号(如图2B所示),或者多个像素输出一个平均PPG信号(如图2B所示),例如以硬件电路计算而得。
处理单元45用于将所述多个PPG信号转换为三维能量分布,例如图2D所示、辨识所述三维能量分布中的环状式样、以及当该环状式样被确认时,例如图6A及6B所示,控制显示器47显示生理检测装置400已经就续(也即穿戴良好)的信息。产生三维能量分布的方式已说明于前。
一种非限定的实施方式中,所述环状式样包含至少一个环,该至少一个环是由三维能量分布中大于能量阈值的能量值所形成,例如图6A及6B中较浅颜色的峰值所形成。从图6A及图6B可看出一个以上的环。处理单元45还可以其他方法得到环,例如计算相邻像素的能量值的差值并找出三维能量分布中的区域极值以作为所述环。
然而,如果三维能量分布中不存在环状式样,则表示生理检测装置400没有穿戴好或足够贴紧以进行生理检测。因此,当无法确认出环状式样时,处理单元45还用于控制显示器47显示改变生理检测装置400的穿戴位置或松紧状态的信息。生理检测装置400配置成提供警示信息给使用者,直到环状式样被确认。在确认环状式样之后,生理检测装置400所检测的三维能量分布则视为包含有效数据。
如前所述,本发明说明的生理检测装置400能检测不同组织深度的浅层微循环。例如,处理单元45还用于控制光源模块41发出比所述第一波长还长的第二波长的发射光,以使得感光阵列43在所述环状式样确认之后,检测来自不同组织深度的出射光。一种非限定的实施例中,第二波长是介于850纳米与900纳米之间或介于590纳米与620纳米之间,并无特定限制。通过改变光波长及分析相对不同光波长的三维能量分布,能够得到被检测组织区域的更多细节信息。一种非限定的实施例中,处理单元45配置成控制显示器47显示改变光波长的信息以得到适合的三维能量分布。
例如,图7A及图7B为相对光源模块41发出的880纳米发射光的三维能量分布的示意图,此处使用480×480个像素的感光阵列,且每个像素尺寸为5微米×5微米。可看出图7A及图7B没有环状式样。这是因为当使用较长波长的光时,发射光穿过较多组织(包括浅层组织及较深层组织),故三维能量分布反映了较复杂的数据。所侦测的与较长波长的光相关的三维能量分布需要较复杂的处理。因此,一种非限定的实施例中,为了简化处理,处理单元45并不用于在与第二波长的光相关的三维能量分布中辨识环状式样。也即,本发明说明中,处理单元45配置成用于根据较短波长光,例如第一波长,确认生理检测装置400是否适当穿戴,但不使用较长波长光,例如第二波长,判断穿戴状况。
此外,处理单元45还用于控制显示器47显示信息,该信息表示移动生理检测装置400的方向来获得有意义的数据。图7A及图7B是使用相同波长光并在不同时间点时,感光阵列43所检测的三维能量分布。该三维能量分布随着时间反复的在图7A及图7B之间变化。图中可看出在较低Y-轴的部分(大约在Y=0到10)始终有不随时间变化的较高能量值,这可视为检测数据超过系统的可检测范围。因此,处理单元45通过显示器47通知使用者朝向正Y-轴方向(例如Y=40)移动生理检测装置400以避免检测始终高能量值的区域。
其他实施例中,处理单元45计算质心位置、重心位置或三维能量分布的中心点。如果求得的位置不位于三维能量分布的中心,处理单元45控制显示器47指引使用者移动生理检测装置400以使得所述位置靠近三维能量分布的中心。也即,处理单元45配置成控制显示器47显示移动生理检测装置400的方向至预定区域,例如具有较多血管的区域的指引。
上述实施例中,当生理检测装置400没有适当穿戴时,生理检测装置400通知用户改变位置或以不同松紧度穿戴。
其他实施例中,生理检测装置400先进行自我调整,而且如果无法通过自我调整满足需求,例如检测出环状式样,才如前所述的通知用户调整位置或松紧度。参照图8,其为本发明另一实施例的生理检测装置400的示意图,其同样包含光源模块、感光阵列43以及处理单元45。本实施例中,生理检测装置400同样包含显示器47。
本实施例的光源模块包含多个发光二极管,例如图8显示发光二极管411至416。光源模块以所述多个发光二极管中的不同组发光二极管朝向组织区域发出第一波长的发射光,其中所述第一波长介于500纳米与550纳米之间。例如,第一组发光二极管包含发光二极管411至413。必须说明的是,发光二极管的数目和感光阵列的尺寸并不限于图8所示。
感光阵列43用于检测来自组织区域的出射光并输出多个PPG信号,例如每个信号如图2B所示。必须说明的是,感光阵列43与多个发光二极管411至416的配置并不限于图8所示,只要当点量不同组的发光二极管时,感光阵列43能够检测来自不同方向的出射光即可。
处理单元45用于将所述多个PPG信号转换为三维能量分布、在点亮第一组发光二极管时获得的三维能量分布中辨识环状式样、并当在所述三维能量分布中无法确认环状式样时控制第二组发光二极管,例如发光二极管414至416,发光。环状式样已说明于上,故于此不再赘述。
本实施例与前一实施例不同之处在于,当在与第一组发光二极管相关的三维能量分布中无法确认环状式样时,处理单元45改变另外的发光二极管群组照明组织区域,但同样使用第一波长。除了改变不同位置的发光二极管发光,处理单元45还选择在感光阵列43获取的图像帧中改变感兴趣窗口(例如图2A的WOI),以获取适当的三维能量分布。如果通过自我调整,例如点量不同发光二极管或改变WOI,即能够检测到环状式样,处理器45则不控制显示器47显示信息来进行手动调整。
如前所述,如果自我调整行的通,生理检测装置400可用于检测不同组织深度的浅层微循环。也即,当在与第一组发光二极管相关的三维能量分布中确认存在环状式样时,处理器45还用于控制光源模块以比第一波长还长的第二波长发光。如前所述,第二波长选择介于850纳米与900纳米之间,或介于590纳米与620纳米之间,但并不限于此。
另一实施例中,生理检测系统包含两个生理检测装置,例如称为阵列式PPG检测器,以监控人体不同部位的浅层微循环。例如参照图9,其为本发明说明再一实施例的生理检测系统500的方框图。
生理检测系统500包含第一阵列式PPG检测器501、第二阵列式PPG检测器503、处理单元505以及显示器507。必须说明的是,虽然图9中显示处理单元505配置于第一阵列式PPG检测器501及第二阵列式PPG检测器503之外,本发明并不限于此。一种非限定的实施例中,处理单元505配置于第一阵列式PPG检测器501或第二阵列式PPG检测器503内部。
第一阵列式PPG检测器501及第二阵列式PPG检测器503包含类似于图4的感光阵列43。本实施例中,第一阵列式PPG检测器501用于产生多个第一PPG信号,且第二阵列式PPG检测器503用于产生多个第二PPG信号。感光阵列产生多个PPG信号(如图2B所示)的方式已说明如上。
例如,第一阵列式PPG检测器501包含第一光源模块及第一感光阵列。所述第一光源模块用于发出第一波长的光以照明第一组织区域。所述第一感光阵列用于检测来自所述第一组织区域的出射光并产生多个第一PPG信号。第二阵列式PPG检测器503包含第二光源模块及第二感光阵列。所述第二光源模块用于发出第二波长的光以照明第二组织区域。所述第二感光阵列用于检测来自所述第二组织区域的出射光并产生多个第二PPG信号。例如,所述第一波长介于500纳米与550纳米之间。
显示器507用于显示生理检测系统的检测结果。
处理单元505使用类似上述的方法将多个第一PPG信号及多个第二PPG信号分别转换为第一三维能量分布及第二三维能量分布。本实施例中,所述第一组织区域例如位于用户的手而所述第二组织区域例如位于用户的脚,并无特定限制,只要两个阵列式PPG检测器位于不同被检测的皮肤表面即可。例如,处理单元505比较第一三维能量分布与第二三维能量分布以决定是否手部或脚部的微循环变差了,例如长时间处于相同坐姿所造成。同理,在进行比较之前,处理单元505确认第一阵列式PPG检测器501及第二阵列式PPG检测器503是否穿戴良好。也即,处理单元505辨识第一三维能量分布及第二三维能量分布中的环状式样。当第一三维能量分布与第二三维能量分布两者分别包含环状式样时,表示生理检测系统500运作正常。
在第一三维能量分布与第二三维能量分布中均确认了环状式样之后,生理检测系统500还可持续监测第一三维能量分布与第二三维能量分布随时间的变化。
一种实施例中,处理单元505分别计算第一三维能量分布的第一平均值并计算第二三维能量分布的第二平均值。处理单元505还计算第一平均值与第二平均值的比值或差值,并监控所述比值或差值的变化。当所述比值或差值改变而超过变化阈值时,表示监测中的两个身体部位的微循环状态不同,此时处理器505配置成控制显示器507显示警示信号以通知用户活动身体。
另一实施例中,处理单元505比较第一三维能量分布及第二三维能量分布与能量阈值,且计算第一三维能量分布中能量大于所述能量阈值的第一区域及第二三维能量分布中能量大于所述能量阈值的第二区域。处理单元505监测所述第一区域与所述第二区域的比值或差值的变化,并当所述比值或差值改变而超过变化阈值时控制显示器507显示警示信号。
必须说明的是,本实施例的生理检测系统可包含两个以上的生理检测装置以监测不同的身体部位,且当所求得的多个三维能量分布之间不平衡或具有明显差异时,处理单元45控制显示器47进行提示。
其他实施例中,所述皮肤浅层微循环血管变化的数据可利用非光学式方式检测,例如多普勒(Doppler)检测,只要能符合分辨率要求即可,例如感测像素尺寸优选介于5×5μm-10×10μm,感测阵列尺寸优选介于240×240-480×480,并不限于光学检测方式。也即,不论生理检测系统是否包含光源,均包含感测阵列以及处理单元。所述感测阵列用于检测皮肤真皮层的阵列式微循环数据以同时反映不同皮肤区域状态;其中,所述感测阵列包含多个像素区域。于光学式检测中,所述多个像素区域则为感光像素;于其他检测方式中,所述多个像素区域为相对应感测像素。所述处理单元则用于根据所述阵列式微循环数据随时间的变化判断不同微循环状态;其中,所述变化例如包含所述微循环数据的变化量的变化以及平均值的变化。
监测微循环的其他特性也有助于早期发现末梢血管病变(peripheral vasculardisease)。例如图6A及图6B的环状式样随时间的往返频率可反应前毛细管括约肌(precapillary sphincter)每分钟启闭前毛细管的频率。本发明说明还提出不使用杜普勒方式检测二相血流量响应(biphasic blood flow response)的微循环检测系统及检测方法。本发明说明的微循环检测系统检测到的环状式样随时间的振荡频率与心律不同,大约5-10次/分钟。
请参照图10A所示,其为本发明说明实施例的微循环检测系统1000的方框图。类似图4的阵列式生理检测系统400,本实施例的微循环检测系统1000同样包含光源1001、感光阵列1003以及处理单元1005,其中光源1001、感光阵列1003及处理单元1005的型式分别相同于上述光源41、感光阵列43及处理单元45,故于此不再赘述。
光源1001及感光阵列1003相对皮肤表面的配置方式例如可参照图4。光源1001用于发光照明皮肤区域,其中,所述光的波长优选介于500纳米与550纳米之间,以利检测前毛细管的组织深度。感光阵列1003用于检测所述皮肤区域的出射光并输出多个PPG信号,其中,每个PPG信号例如可参照图2B。感光阵列1003包含阵列排列(例如图2A所示)的多个像素区域用于分别输出亮度变化信号以作为所述多个PPG信号的其中一者,其中所述多个像素区域的每一者包含至少一个感光像素。当一个像素区域包含多个感光像素时,感光阵列1003具有电路对一个像素区域的多个感光像素的检测信号进行加法运算并输出一个亮度变化信号总和以作为该像素区域的PPG信号。
处理单元1005同样将所述多个PPG信号转换为阵列能量分布(如图2D所示的三维能量分布)并辨识所述阵列能量分布中的环状式样,例如辨识第一时间点的所述阵列能量分布中的第一环状式样ED1(或图6A所示者)及第二时间点的所述阵列能量分布中的第二环状式样ED2(或图6B所示者),其中第一时间点不同于第二时间点。如前所述,所述阵列能量分布是所述多个PPG信号在预定频率(例如心律或其倍频)的频谱能量相对所述多个像素区域的二维空间的能量值分布。所述阵列能量分布中的环状式样会随时间振荡,例如在ED1及ED2之间反复振荡。当检测中皮肤区域的皮肤表面温度大致维持固定时,所述环状式样的振荡频率(次/分钟)大致维持固定。
此外,本实施例的微循环检测系统1000还包含加热装置1002用于对所述皮肤区域加温以及计时器1006用于计时所述加热装置1002的加温期间。计时器1006可从已知装置选择并无特定限制,只要能受到处理单元1005的控制以当所述加热装置1002开始加热时开始计时即可。处理单元1005还可在每次开始计时前重置计时器1006。
请参照图10B所示,其为本发明说明实施例的微循环检测系统1000的运作示意图。一种非限定的实施方式中,加热装置1002例如包含腔室1021及加热器1022。所述腔室1021用于容纳欲检测皮肤区域。例如,当欲检测皮肤区域位于使用者的手部时,腔室1021具有开口1023以供使用者的手通过所述开口1023伸进腔室1021中。可以了解的是,当所述腔室1021用于容纳其他身体部位(例如脚)时,所述开口1023可位于腔室1021的不同表面,以利使用者置入身体部位以供检测。
加热器1022例如可为红外线灯或电热管,其设置于腔室1021内并用于加温腔室1021内的气体。红外线灯还可通过辐射热加温欲检测皮肤区域。可以了解的是,当使用者的手置入腔室1021中时,通过加温腔室1021内的气体可均匀的加温欲检测皮肤区域。其他实施方式中,也可选择使用接触式加热直接对欲检测皮肤区域进行加温。更详言之,加热器1022的型式并无定限制,只要其能加热皮肤表面温度即可。
同时,为了记录欲检测皮肤区域的皮肤表面温度,本实施例的微循环检测系统1000还包含温度传感器1043用于量测皮肤表面温度。均匀加热使用者的手部的优点是,温度传感器1043量测不同手指(例如图10B显示第四指)的皮肤表面温度即可视为欲检测皮肤区域(例如图10B显示置放于光源1001及感光阵列1003上的第二指)的皮肤表面温度。温度传感器1043的量测温度Ts则传送至缓冲存储1051(buffer memory)纪录以供处理单元1005存取。此外,微循环检测系统1000还可另包含温度传感器1041用于量测腔室1021内的腔室温度Tc并传送至缓冲存储1051纪录以供处理单元1005存取。该缓冲存储1051可为挥发性存储并包含于处理单元1005之内或之外,并无特定限制。
参照图11所示,其显示本实施例的微循环检测系统1000的腔室温度Tc及加温皮肤温度(例如图10B显示为左手)在加温期间的示意图。某些实施例中,可选择性地使用另一温度传感器对常温皮肤温度(例如右手)进行温度纪录以确认皮肤加温过程是否正常。一般而言,加温期间可依据不同使用者设定为12至15分钟。从图11中可看出随着腔室温度Tc的上升,加温皮肤温度会被加温至较高的皮肤表面温度。
在使用加热装置1002对欲检测皮肤区域进行加温的过程中,持续以感光阵列1003检测所述皮肤区域的出射光并在不同时间点输出多个PPG信号,且每一时间点的多个PPG信号用以形成一个阵列能量分布。处理单元1005则以预定频率(例如帧率)辨识每个阵列能量分布中的环状式样,并计算所述环状式样的振荡在加温期间内的频率变化。一种非限定的实施方式中,所述环状式样的振荡是以所述阵列能量分布的环状式样中与所述多个像素区域至少一者相对应的位置的能量振幅的振荡来判定。
另一种非限定的实施方式中,所述环状式样的振荡是以第一环状式样ED1及第二环状式样ED2的式样振荡判定。例如,处理单元1005先选定并储存第一环状式样ED1及第二环状式样ED2(例如储存于帧缓冲器中),其中所述第一环状式样ED1及所述第二环状式样ED2大致彼此反相。在加温期间中,处理单元1005将每次辨识出的环状式样与储存的第一环状式样ED1及第二环状式样ED2进行相似度或相关度比对,以确认其往复变化。当在所述加温期间辨识出某些环状式样依序属于所述第一环状式样ED1(例如相似度或相关度高于阈值)及所述第二环状式样ED2,即可计算出振荡频率。
例如参照图12所示,其为本发明说明实施例的微循环检测系统1000检测的频率变化的示意图。图12显示所述环状式样的振荡频率在开始加热后的第2分钟及第11分钟分别出现峰值(即最高振荡频率),此现象通称为二相流量响应。为了检测此二相(即二峰值),处理单元1005还判断频率变化(例如每分钟往复于ED1及ED2的次数的变化)是否在所述加温期间的第一时间区间及第二时间区间分别具有峰值,其中所述第一时间区间选择为所述加温期间的第2至5分钟;所述第二时间区间选择为所述加温期间的第10至15分钟。图12中还同时显示了常温皮肤的频率变化以作为对照组。于实际运作时,微循环检测系统1000可以仅纪录加温皮肤的环状式样的频率变化而不记录常温皮肤(即未加温皮肤)的频率变化。
此外,本发明实施例的微循环检测系统1000还另包含提示装置1008用以通过图像、声音、震动、灯号、无线电波等方式提示检测结果。例如,当提示装置1008为显示器时,该显示器可用于显示不同时间点的环状式样(例如图10A的第一环状式样ED1及第二环状式样ED2)、频率变化(例如以图12的柱状图或线形图)、加温时间(例如以数字或线形图表示)、腔室温度及皮肤表面温度(例如以数字或线形图表示)至少其中之一。处理单元1005当判断所述第一时间区间及所述第二时间区至少其中一者不具有峰值时,表示使用者的微循环的血流调控可能不正常,则控制提示装置1008发出提示。
请参照图13所示,其为本发明实施例的微循环检测系统的检测方法的流程图,该检测方法例如适用于图10A及图10B的微循环检测系统1000。本实施例的检测方法包含:以光源照明皮肤区域(步骤S131);以加热装置对所述皮肤区域加温(步骤S132);以感光阵列检测所述皮肤区域的出射光并输出多个PPG信号(步骤S133);以处理单元将所述多个PPG信号转换为阵列能量分布(步骤S134);以及以所述处理单元辨识所述阵列能量分布中预定位置的能量振荡在加温期间的频率变化在所述加温期间的两个预定时间区间中的两个峰值(步骤S135)。
首先,开启光源1001照射皮肤区域并以加热装置1002对所述皮肤区域加热(步骤S131-132)。一种非限定的实施例中,光源1001可设定成当加热装置1002开启时才开始发光。光源1001、温度传感器1043及感光阵列1003可直接设置于腔室1021内部,或者先配置于用户的身体部位后再置入腔室1021内。
一种非限定的实施例中,感光阵列1003设定成当加热装置1002开始加温腔室1021时才开始输出多个PPG信号,其中PPG信号的数目是根据像素区域的数目而决定(步骤S133)。
处理单元1005则将所述多个PPG信号在预定频率的频谱能量分布于相对所述多个像素区域的二维空间以形成阵列能量分布,例如图10A中的ED1及ED2为不同时间点的阵列能量分布(步骤S134)。
处理单元1005配置为选定阵列能量分布中的预定位置,例如相对所述多个像素区域其中至少一者能量振荡超过阈值的位置,其可为环状式样的质心位置、重心位置或中心点等位置。处理单元1005计算该预定位置的振荡频率在加温期间内的频率变化,如图12所示。如前所述,处理单元1005判断加温期间的两个预定时间区间(第2至5分钟及第10至15分钟)中的两个峰值(步骤S135)。当至少一个预定时间区间不存在频率变化的峰值时,处理单元1005则控制提示装置1008进行提示。此外,也可使用显示器1008显示处理单元1005的运算结果,包括环状式样的图像、频率变化的数字或图形、记录温度以及加温时间等。处理单元1005还可通过通讯界面或网络将运算结果传送至外部装置。
本发明说明的检测方法中,处理单元1005可以如上所述的根据两个选定的环状式样之间的振荡来计算频率变化。
本实施例中,优选在能够辨识出环状式样之后才开始纪录频率变化。若处理单元1005没办法辨识出环状式样,则可以上述的改变点亮光源位置、调整穿戴松紧度等方式以获取包含环状式样的阵列能量分布后,再开始加温及检测皮肤区域。本发明说明的环状式样类似石头掉入水中时所形成的水波纹。
本发明说明可适用于经皮药物释放系统监测。经皮药物释放系统(transdermaldrug delivery system)是指药物在经皮肤给药后,以一定的速率通过皮肤,经过微循环血管吸收后,进入人体循环而产生药效的一种投药机制。好处是可以避免肝脏的首过反应以及肠胃道对于药物的破坏,进而可以达成减少投药次数、延长投药间的时间间隔、维持血液中有效血药浓度等效果,以提高疗效。
本发明说明可用于监测微循环血管对药物的吸收反应,当微循环血管的振幅变化升高同时心跳频率也增快时,可判定为经皮药物释放系统持续在作用当中。当微循环血管的振幅变化、心跳频率、振幅平均数值都回复过往的正常范围时,便可得知经皮药物释放系统已经作用完毕,可以进行后续的疗程,例如再次投药等等。换句话说,本发明说明的阵列式生理检测系统可通过三维频谱能量反映微循环的投药状态,其可显示投药效果。
糖尿病患者的四肢容易发生的有动脉粥状硬化以及末梢神经病变等疾病。动脉粥状硬化会导致组织缺血与坏死,末梢神经病变会导致运动衰弱并且感觉丧失。由于微循环血管受到交感神经支配,因此监测微循环血管的变化,可以在早期警示糖尿病患者是否发生前述疾病。
本发明说明可以监测微循环血管的变化,当糖尿病患者的微循环血管显示其振幅变化与振幅平均数值随着时间降低,便可得知血管正逐渐失去功能。换句话说,本发明说明的阵列式生理检测系统可通过三维频谱能量反映微循环的退化状态,其可显示病灶程度。
本发明说明也可在对患者施加外加刺激时,观察患者的微循环反应,例如观察患者有无末梢神经病变时,会施加外部冷热刺激,此时若微循环呈现振幅变化降低,同时心跳频率也不会升高时,便代表末梢神经不活跃,可能出现病变。换句话说,本发明说明的阵列式生理检测系统可通过三维频谱能量反映微循环的反应状态,其可显示神经活跃度。
烧烫伤病患因为局部面积失去皮肤保护,容易出现低血容量性休克、微循环血管脆性与通透性增加等状况。由于此类状况发展快,因此如果抢救不及,可能出现多器官功能障碍综合症。利用本发明说明可以监测烧烫伤病患的末梢组织循环,藉此监测病患的病程变化,以避免更危险的状况发生。换句话说,本发明说明的阵列式生理检测系统可通过三维频谱能量反映微循环的运作状态,其可显示病程变化。
高压氧治疗在临床上证明可有效改善放射线照射后的组织微循环,对于放射性骨坏死或者软组织坏死具有显着的疗效。患者在接受高压氧治疗时,可藉由本发明说明来监测治疗效果变化,当微循环血管的振幅变化增加,且心跳频率也增加时,代表微循环血管逐渐恢复活性,治疗逐渐发生效果。换句话说,本发明说明的阵列式生理检测系统可通过三维频谱能量反映微循环的恢复状态,其可显示治疗效果。
休克是一种进行性的过程,当循环系统失去支持身体代谢能力,造成身体组织或者器官的血流灌注不足时,送到身体的氧气无法让各部分组织充分利用,以至于细胞代谢产生异常,引起细胞损伤或者死亡。当病患开始发生休克状态时,将引起微循环血管扩张,此时血液积蓄于微循环内,若无法有效排除,则可能导致严重的休克状况。
通过应用本发明说明,可以在对病患进行休克缓解处置时,同步观测是否发生效果。若一直显示微循环血管的振幅信号的平均数值相当高,且振幅信号的振幅变化不高,同时心跳频率持续维持在较高频率,则代表缓解处置尚未发生效果,反之亦然。换句话说,本发明说明的阵列式生理检测系统可通过三维频谱能量反映微循环的恢复状态,其可显示缓解效果。
在运动过量的状态中典型的有热衰竭以及中暑反应。当人体处于这类状态下,皮肤血液循环会增加,此时心脏所输送出去的血液也要跟着增加。当血液不足的时候,体内的血液会重新分配,使内脏器官的血液循环减少,皮肤血液循环增加协助排汗,把体内的热散逸出去。配合在运动中使用本发明说明,可以在显示微循环血管的振幅信号的平均数值相当高,且振幅信号的振幅变化不高,同时心跳频率持续维持在较高频率时,适度提醒用户可能已经处于运动过量的状态下,而不适宜持续运动。换句话说,本发明说明的阵列式生理检测系统可通过三维频谱能量反映微循环的血液分配状态,其可显示散热效果。
微循环具有调节组织血流量、供给细胞营养、排除代谢产物等功能,而局部血量多寡可以代表相对的温度变化,通过本发明说明,可以检测末梢组织微循环温度的相对变化,当局部组织的微循环血管振幅信号的平均数值高时,代表温度升高,反之亦然。换句话说,本发明说明的阵列式生理检测系统可通过三维频谱能量反映微循环的温度状态,其可显示局部血量。
目前为止,仍没有有效的可以方便携带的末稍微循环血管检测产品可供自闭症患者/婴儿/宠物末梢组织交感神经检测。在微循环血管的动脉与小动脉的管壁由平滑肌构成,其受交感神经支配,其控制着微循环血管的开放与关闭,藉此来决定组织的血液供应量。通过本发明说明,可以藉由观察末稍微循环血管的变化趋势,来间接推测交感神经活性状态。当交感神经活跃时,微循环血管的变化趋势也会趋向活跃,反之亦然。换句话说,本发明说明的阵列式生理检测系统可通过三维频谱能量反映微循环的供血状态,其可显示交感神经的活跃状态。
本发明说明也可应用在判断心脏功能或者全身性的血管缺损或硬化。将所有的浅层微血管涨缩变化信号整合为一个单一结果后,在不同频率之下会具有不同的能量。一般而言,代表能量的信号应该出现在心跳频率的倍频,同时信号能量会维持在一个正常的区间,这个区间因人而异,但针对同一个使用者,其随时间的变化性应该不会太大。因此,如果代表能量的信号出现在心跳频率的倍频以外的区间,又例如当随着时间流逝,信号的能量脱离正常的区间,则代表用户的心脏或者血管功能有异常,需要进行进一步检查。
举例而言,当代表能量的信号出现在心跳频率的倍频以外的区间,可能代表用户的心脏功能缺损,例如瓣膜缺损。当随着时间流逝,信号的能量超过正常区间甚多,可能代表使用者的血管硬化,因此心脏需要增大其输出的功率,来将血液运送到全身。换句话说,本发明说明的阵列式生理检测系统可通过三维频谱能量反映微循环的异常状态,其可显示心脏功能。
上述说明中,振幅变化指三维频谱能量的变化,而振幅平均数值指三维频谱能量的平均值;其中,三维频谱能量的分布类似于图2D。上述说明中,有效像素数目指该观测窗口WOI范围内的像素数目。本发明说明实施例中所举出的数值仅用于说明,并非用于限定本发明。
虽然本发明说明已通过前述实例披露,但是其并非用于限定本发明说明,任何本发明说明所属技术领域中具有通常知识的技术人员,在不脱离本发明说明的精神和范围内,当可作各种的更动与修改。因此本发明说明的保护范围当视后附的权利要求所界定的范围为准。
Claims (18)
1.一种微循环检测系统,该微循环检测系统包含:
光源,用于发光照明皮肤区域,其中所述光的波长介于500纳米与550纳米之间;
加热装置,用于对所述皮肤区域加温;
感光阵列,用于检测所述皮肤区域的出射光并输出多个PPG信号;及
处理单元,用于
将所述多个PPG信号转换为阵列能量分布,
辨识所述阵列能量分布中的环状式样,及
计算加温期间内的所述环状式样的振荡的频率变化,
其中所述感光阵列包含阵列排列的多个像素区域用于分别输出亮度变化信号以作为所述多个PPG信号。
2.根据权利要求1所述的微循环检测系统,还包含显示器用于显示不同时间点的所述环状式样、所述频率变化、加温温度及皮肤表面温度至少其中之一。
3.根据权利要求1所述的微循环检测系统,还包含计时器用于计时所述加温期间。
4.根据权利要求1所述的微循环检测系统,其中所述加温期间为15分钟。
5.根据权利要求1所述的微循环检测系统,其中所述加热装置包含:
腔室,用于容纳所述皮肤区域;以及
红外线灯,设置于所述腔室内并用于加温所述腔室。
6.根据权利要求1所述的微循环检测系统,其中所述多个像素区域的每一者包含至少一个感光像素。
7.根据权利要求6所述的微循环检测系统,其中所述阵列能量分布是所述多个PPG信号在预定频率的频谱能量相对所述多个像素区域的二维空间的能量值分布。
8.根据权利要求7所述的微循环检测系统,其中所述环状式样的振荡是所述阵列能量分布的所述环状式样中与所述多个像素区域至少一者相对应的位置的能量振幅的振荡。
9.一种微循环检测系统,该微循环检测系统包含:
光源,用于发光照明皮肤区域,其中所述光的波长介于500纳米与550纳米之间;
加热装置,用于对所述皮肤区域加温;
计时器,用于计时所述加热装置的加温期间;
感光阵列,用于检测所述皮肤区域的出射光并输出多个PPG信号;及
处理单元,用于
将所述多个PPG信号转换为阵列能量分布,
辨识第一时间点的所述阵列能量分布中的第一环状式样及第二时间点的所述阵列能量分布中的第二环状式样,及
判断所述第一环状式样及所述第二环状式样的往复频率的频率变化是否在所述加温期间的第一时间区间及第二时间区间分别具有峰值,其中所述感光阵列包含阵列排列的多个像素区域用于分别输出亮度变化信号以作为所述多个PPG信号。
10.根据权利要求9所述的微循环检测系统,其中
所述第一时间区间为所述加温期间的第2至5分钟;且
所述第二时间区间为所述加温期间的第10至15分钟。
11.根据权利要求9所述的微循环检测系统,还包含提示装置,所述处理单元还用于
当判断所述第一时间区间及所述第二时间区间至少其中一者不具有所述峰值时,控制所述提示装置发出提示。
12.根据权利要求11所述的微循环检测系统,其中所述提示装置为显示器,该显示器用于显示所述加温期间中的所述频率变化。
13.根据权利要求9所述的微循环检测系统,其中所述加热装置包含:
腔室,用于容纳所述皮肤区域;以及
红外线灯,设置于所述腔室内并用于加温所述腔室。
14.根据权利要求9所述的微循环检测系统,其中所述多个像素区域的每一者包含至少一个感光像素。
15.根据权利要求14所述的微循环检测系统,其中所述阵列能量分布是所述多个PPG信号在预定频率的频谱能量相对所述多个像素区域的二维空间的能量值分布。
16.一种微循环检测系统的检测方法,该微循环检测系统包含光源、加热装置、感光阵列以及处理单元,所述检测方法包含:
以所述光源照明皮肤区域,其中所述光源发光的波长介于500纳米与550纳米之间;
以所述加热装置对所述皮肤区域加温;
以所述感光阵列检测所述皮肤区域的出射光并输出多个PPG信号;
以所述处理单元将所述多个PPG信号转换为阵列能量分布;以及
以所述处理单元辨识所述阵列能量分布中预定位置的能量振荡在加温期间的频率变化在所述加温期间的两个预定时间区间中的两个峰值,
其中所述感光阵列包含多个像素区域用于分别输出亮度变化信号以作为所述多个PPG信号。
17.根据权利要求16所述的检测方法,其中
所述阵列能量分布是所述多个PPG信号在预定频率的频谱能量相对所述多个像素区域的二维空间的能量值分布。
18.根据权利要求17所述的检测方法,其中所述预定位置是所述阵列能量分布中,相对所述多个像素区域其中至少一者能量振荡超过能量阈值的位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210836923.5A CN115251876A (zh) | 2018-06-04 | 2019-06-03 | 生理检测装置 |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/996,652 | 2018-06-04 | ||
US15/996,652 US11020013B2 (en) | 2015-01-29 | 2018-06-04 | Array physiological detection system and device |
US16/355,864 | 2019-03-18 | ||
US16/355,864 US11311203B2 (en) | 2015-01-29 | 2019-03-18 | Microcirculation detection system |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210836923.5A Division CN115251876A (zh) | 2018-06-04 | 2019-06-03 | 生理检测装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110547785A CN110547785A (zh) | 2019-12-10 |
CN110547785B true CN110547785B (zh) | 2022-08-05 |
Family
ID=68736395
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210836923.5A Pending CN115251876A (zh) | 2018-06-04 | 2019-06-03 | 生理检测装置 |
CN201910477770.8A Active CN110547785B (zh) | 2018-06-04 | 2019-06-03 | 微循环检测系统及检测方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210836923.5A Pending CN115251876A (zh) | 2018-06-04 | 2019-06-03 | 生理检测装置 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN115251876A (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105832307A (zh) * | 2015-01-29 | 2016-08-10 | 原相科技股份有限公司 | 阵列式生理检测系统及其运作方法 |
CN105877723A (zh) * | 2016-06-22 | 2016-08-24 | 中国科学院苏州生物医学工程技术研究所 | 无创连续血压测量装置 |
CN106880351A (zh) * | 2015-12-15 | 2017-06-23 | 德州仪器公司 | 减少光体积描记(ppg)信号中运动诱发的伪影 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI552719B (zh) * | 2014-12-23 | 2016-10-11 | 原相科技股份有限公司 | 三維生理偵測系統及其運作方法 |
-
2019
- 2019-06-03 CN CN202210836923.5A patent/CN115251876A/zh active Pending
- 2019-06-03 CN CN201910477770.8A patent/CN110547785B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105832307A (zh) * | 2015-01-29 | 2016-08-10 | 原相科技股份有限公司 | 阵列式生理检测系统及其运作方法 |
CN106880351A (zh) * | 2015-12-15 | 2017-06-23 | 德州仪器公司 | 减少光体积描记(ppg)信号中运动诱发的伪影 |
CN105877723A (zh) * | 2016-06-22 | 2016-08-24 | 中国科学院苏州生物医学工程技术研究所 | 无创连续血压测量装置 |
Also Published As
Publication number | Publication date |
---|---|
CN115251876A (zh) | 2022-11-01 |
CN110547785A (zh) | 2019-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6608388B2 (ja) | ユーザの心拍数の変動を検出するためのシステム及び方法 | |
US20220183577A1 (en) | Microcirculation detection system and method | |
US11213217B2 (en) | Physiological detection system and operating method thereof | |
EP3030145B1 (en) | Monitoring system and method for monitoring the hemodynamic status of a subject | |
US6261236B1 (en) | Bioresonance feedback method and apparatus | |
US20180360373A1 (en) | Wearable device and method for determining electro-dermal activity of a subject | |
Liu et al. | PCA-based multi-wavelength photoplethysmography algorithm for cuffless blood pressure measurement on elderly subjects | |
RU2724426C2 (ru) | Устройство и система мониторинга глаза субъекта | |
CN105832307B (zh) | 阵列式生理检测系统及其运作方法 | |
US11918329B2 (en) | Array physiological detection system | |
US20100262025A1 (en) | Apparatus for measurement of heart rate variability | |
Mishra et al. | Blood oxygen saturation measurement using polarization-dependent optical sectioning | |
CN114947765A (zh) | 信号检测方法和相关产品 | |
CN109793510B (zh) | 阵列式生理检测系统及生理检测装置 | |
CN110547785B (zh) | 微循环检测系统及检测方法 | |
US10827972B1 (en) | Sleep staging using an in-ear photoplethysmography (PPG) | |
US11583201B2 (en) | Device for recording the vascular response of the human spinal cord triggered by a suprasensible stimulus through the use of functional near-infrared spectroscopy | |
US10610164B2 (en) | Determining changes to autoregulation | |
Ishibashi et al. | Effects of mental task on heart rate variability during graded head-up tilt | |
US10646158B2 (en) | Method of using autonomic nerve fluctuations and system using the same | |
CN114224309A (zh) | 生理检测系统及生理检测装置 | |
US20100217138A1 (en) | Recognizing awareness in an anesthetized subject | |
KR20190081527A (ko) | 혈압 모니터링을 위한 이어폰 및 이를 이용한 혈압 모니터링 방법 | |
JP2023506542A (ja) | 対象者の酸素飽和度を特定するためのデバイス、システム、及び方法 | |
Khalil et al. | Measuring peripheral vascular reactivity with diffusive optical imaging |
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 |