CN111436905A - 一种光学相干层析血管造影方法及装置 - Google Patents
一种光学相干层析血管造影方法及装置 Download PDFInfo
- Publication number
- CN111436905A CN111436905A CN202010259880.XA CN202010259880A CN111436905A CN 111436905 A CN111436905 A CN 111436905A CN 202010259880 A CN202010259880 A CN 202010259880A CN 111436905 A CN111436905 A CN 111436905A
- Authority
- CN
- China
- Prior art keywords
- blood flow
- scan
- image
- layering
- images
- 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
- 238000012014 optical coherence tomography Methods 0.000 title claims abstract description 57
- 238000002583 angiography Methods 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000017531 blood circulation Effects 0.000 claims abstract description 98
- 210000001525 retina Anatomy 0.000 claims abstract description 29
- 238000010586 diagram Methods 0.000 claims abstract description 28
- 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 abstract description 24
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 5
- 239000006185 dispersion Substances 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 6
- 238000012800 visualization Methods 0.000 claims description 2
- 238000003745 diagnosis Methods 0.000 abstract description 9
- 201000010099 disease Diseases 0.000 abstract description 9
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 9
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 239000010410 layer Substances 0.000 description 25
- 238000003384 imaging method Methods 0.000 description 7
- 210000004204 blood vessel Anatomy 0.000 description 6
- 210000001508 eye Anatomy 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 230000002207 retinal effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- LRSYFEZBIMVWRY-VWMHFEHESA-N (2s)-2,5-diaminopentanoic acid;2-phenylacetic acid Chemical compound NCCC[C@H](N)C(O)=O.OC(=O)CC1=CC=CC=C1 LRSYFEZBIMVWRY-VWMHFEHESA-N 0.000 description 2
- 210000003161 choroid Anatomy 0.000 description 2
- 239000002872 contrast media Substances 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000010412 Glaucoma Diseases 0.000 description 1
- 102100023057 Neurofilament light polypeptide Human genes 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 210000000601 blood cell Anatomy 0.000 description 1
- 210000005252 bulbus oculi Anatomy 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 108010090677 neurofilament protein L Proteins 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 208000019553 vascular disease Diseases 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/0062—Arrangements for scanning
- A61B5/0066—Optical coherence imaging
-
- 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/0073—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
-
- 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)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Physiology (AREA)
- Hematology (AREA)
- Cardiology (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本发明提供了一种光学相干层析血管造影方法及装置，其中方法包括：采集视网膜的OCT原始数据，获取B‑scan图像并对每组图像配准；根据每组内在同一坐标点上的最大图像值与最小图像值的差值的平方计算所述坐标点的血流信息，进而获取当前组对应的血流信号图；将每组图像合成高清图像并进行分层；根据血流信号图及分层结果获取各层之间的OCTA。本发明对噪声部分影响较小，而血流信号会得到增强，增加了血流信号与噪声的区分度，使得到的血管造影图质量更高，毛细血管细节呈现更丰富，从而有助于医生对患者眼底疾病的诊断。同时，本发明中的算法复杂度低，计算量小，减小了计算机负担，提升了血流信号的计算速度，从而也提升了工作效率。
Description
技术领域
本发明涉及光学相干层析成像领域，尤其涉及一种光学相干层析血管造影方法及装置。
背景技术
OCT(Optical Coherence Tomography，光学相干层析成像，或称光相干断层成像)是一种高灵敏度、高分辨率、高速度、无入侵的成像方式，已被广泛应用于眼底疾病的诊断中，对眼科疾病的检测和治疗有着重要意义。作为一种断层扫描成像方式，OCT利用光的相干性对眼底扫描成像，每次扫描被称为一个A-scan，相邻连续的多个扫描组合在一起称为一个B-scan，B-scan也即通常所看到的OCT截面图，是医学诊断中OCT最主要的成像方式。
在OCT领域中，OCTA(OCT Angiography，光学相干层析血管造影)是近几年发展起来的新技术，视网膜OCTA图可以直观的显示视网膜的血流形状、血流分布等情况，对于眼底疾病的诊断有着重要作用。在静止的眼球里，眼底唯一运动的物质是血管里流动的血细胞，故OCTA可以利用血管内流动的红细胞作为造影剂，测量运动细胞引起的OCT信号差异，提供血流信息，快速无创地重建视网膜脉络膜血管的三维结构，从而实现血管网络的可视化。此方法不需要注射荧光造影剂，无入侵无接触，不会对人体造成损坏，成像效果不弱于荧光剂造影，因此迅速地被应用于临床中各种眼底血管性疾病的检查，例如糖尿病眼底病变、青光眼等眼底疾病均可通过视网膜OCTA图进行筛查和诊断。
在OCTA中，血管造影时需要对同一位置重复拍摄多次，计算多次拍摄的B-scan间差异来得到当前位置的血流信号。由于OCTA图不是直接拍摄得到，而是对多次拍摄的数据进行重构，所以过程中引入的噪声较大，导致图像不够清晰。在现有技术中，虽然有一些方法已经能够实现血管造影成像，并对血管造影图像进行分析，然而发明人在实现本发明的过程中发现，这些方法并未侧重图像的质量，其缺点是血管造影图像信噪比较低，噪声较大，图像清晰度较差，无法分辨视网膜中的毛细血管，不利于医生的诊断。
发明内容
本发明提供一种光学相干层析血管造影方法及装置，以解决血管造影中血流信号较弱、噪声较大，无法获取较为清晰的血管图像的问题。
根据本发明实施例的第一方面，提供一种光学相干层析血管造影方法，所述方法包括：
采集视网膜的光学相干层析OCT原始数据，根据所述原始数据获取B-scan图像并对每组B-scan图像进行配准，其中，对应所述视网膜同一位置的多张B-scan图像为同一组；
对于每组B-scan图像：根据组内在同一坐标点上的最大图像值与最小图像值的差值的平方，计算所述坐标点的血流信息，以及，根据各坐标点的血流信息获取当前组对应的B-scan血流信号图；
将每组B-scan图像合成一张高清图像，对合成后的高清图像进行分层；
根据获取的B-scan血流信号图及所述分层的结果，获取各层之间的光学相干层析血管造影OCTA。
可选的，根据所述原始数据获取B-scan图像，包括：
对采集到的所述原始数据进行指定操作以获取所述B-scan图像，其中所述指定操作包括以下一种或多种操作：色散补偿、加窗、傅里叶变换、取阈值。
可选的，对每组B-scan图像进行配准，包括：
根据傅里叶梅林变换计算同一组内B-scan图像的偏移，根据所述偏移对所述同一组B-scan图像进行配准。
可选的，根据组内在同一坐标点上的最大图像值与最小图像值的差值的平方，计算所述坐标点的血流信息，包括：
根据
D(x,z)＝(MAX{I1(x,z),...,In(x,z)}-MIN{I1(x,z),...,In(x,z)})2
计算坐标点(x,z)的血流信息D(x,z)，其中n为每组B-scan图像中的图像数量，I(x,z)为坐标点(x,z)位置处的图像值。
可选的，对合成后的高清图像进行分层，包括：
计算合成后高清图像的垂直梯度，通过使用Dijikstra算法对梯度图计算最短路径并不断限制分层区域的方式获取分层线以得到分层结果。
可选的，根据获取的B-scan血流信号图及所述分层的结果，获取各层之间的光学相干层析血管造影OCTA，包括：
将每张B-scan血流信号图任意层与层之间的血流信息分别在纵向上累加求和，得到血管造影图，再对所述血管造影图进行归一化处理得到对应层的OCTA。
根据本发明实施例的第二方面，提供一种光学相干层析血管造影装置，所述装置包括：
图像获取模块，用于采集视网膜的光学相干层析OCT原始数据，根据所述原始数据获取B-scan图像并对每组B-scan图像进行配准，其中，对应所述视网膜同一位置的多张B-scan图像为同一组；
计算模块，用于对于每组B-scan图像，根据组内在同一坐标点上的最大图像值与最小图像值的差值的平方，计算所述坐标点的血流信息，以及，根据各坐标点的血流信息获取当前组对应的B-scan血流信号图；
分层模块，用于将每组B-scan图像合成一张高清图像，对合成后的高清图像进行分层；
造影模块，用于根据获取的B-scan血流信号图及所述分层的结果，获取各层之间的光学相干层析血管造影OCTA。
可选的，所述图像获取模块在根据所述原始数据获取B-scan图像时，具体用于：
对采集到的所述原始数据进行指定操作以获取所述B-scan图像，其中所述指定操作包括以下一种或多种操作：色散补偿、加窗、傅里叶变换、取阈值。
可选的，所述图像获取模块在对每组B-scan图像进行配准时，具体用于：
根据傅里叶梅林变换计算同一组内B-scan图像的偏移，根据所述偏移对所述同一组B-scan图像进行配准。
可选的，所述计算模块在根据组内在同一坐标点上的最大图像值与最小图像值的差值的平方，计算所述坐标点的血流信息时，具体用于：
根据
D(x,z)＝(MAX{I1(x,z),...,In(x,z)}-MIN{I1(x,z),...,In(x,z)})2
计算坐标点(x,z)的血流信息D(x,z)，其中n为每组B-scan图像中的图像数量，I(x,z)为坐标点(x,z)位置处的图像值。
可选的，所述分层模块在对合成后的高清图像进行分层时，具体用于：
计算合成后高清图像的垂直梯度，通过使用Dijikstra算法对梯度图计算最短路径并不断限制分层区域的方式获取分层线以得到分层结果。
可选的，所述造影模块具体用于：
将每张B-scan血流信号图任意层与层之间的血流信息分别在纵向上累加求和，得到血管造影图，再对所述血管造影图进行归一化处理得到对应层的OCTA。
本发明实施例提供的技术方案可以包括以下有益效果：
发明人在实现本发明的过程中发现，OCT血管造影技术最关键的部分是血流信号的计算，如果采用去相关和求方差的方法，虽然将噪声进行平均处理可减小其带来的影响，但同时会对真正的血流信号也进行了平均处理，使得噪声和真正的血流信号之间区分不明显。而由于噪声波动比血流信号的波动要小，若计算在同一坐标点上的最大图像值与最小图像值的差值，则对噪声部分的大小影响较小，真正的血流信号会得到很大程度的增强，使得噪声相对减弱，增加了血流信号与噪声的区分度；再对所得差值进行平方，则可进一步增加了信噪比，使血流信号更加清晰，进而使得到的血管造影图质量更高，图像中的噪点明显减少，毛细血管细节呈现更丰富，从而有助于医生对患者眼底疾病的诊断。
同时，对比去相关和求方差的方法，在本发明方案中算法复杂度降低，计算量大大减小，减小了计算机负担，很大程度上提升了血流信号的计算速度，拍摄完毕几秒钟内就能获得高清的血管造影图像，从而也提升了工作效率。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本发明。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。此外，这些介绍并不构成对实施例的限定，附图中具有相同参考数字标号的元件表示为类似的元件，除非有特别申明，附图中的图不构成比例限制。
图1是根据本发明一示例性实施例示出的一种光学相干层析血管造影方法的流程图；
图2是根据本发明一示例性实施例示出的扫频OCT系统示意图；
图3是根据本发明一示例性实施例示出的一种光学相干层析血管造影方法的流程图；
图4是根据本发明一示例性实施例示出的B-scan血流信号图；
图5是根据本发明一示例性实施例示出的B-scan分层示意图；
图6是根据本发明一示例性实施例示出的ILM-IPL层之间血管造影结果；
图7是根据本发明一示例性实施例示出的IPL-OPL层之间血管造影结果；
图8是根据本发明一示例性实施例示出的一种光学相干层析血管造影装置的示意图。
具体实施方式
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本发明的一些方面相一致的装置和方法的例子。
图1是根据本发明一示例性实施例示出的一种光学相干层析血管造影方法的流程图。所述方法可用于由OCT设备、计算机等组成的系统。
参见图1所示，所述方法可以包括：
步骤S101，采集视网膜的光学相干层析OCT原始数据，根据所述原始数据获取B-scan图像并对每组B-scan图像进行配准，其中，对应所述视网膜同一位置的多张B-scan图像为同一组。
例如扫描视网膜上的M个位置，则可以得到M组，若每个位置扫描N次，则每组中包含N张B-scan图像，一共可获取到M·N张B-scan图像。
对应所述视网膜同一位置的多张B-scan图像为同一组，或者换句话说，同一组B-scan图像包括了对所述视网膜同一位置拍摄得到的多张B-scan图像。
对于具体如何采集OCT原始数据以及如何根据所述原始数据获取视网膜的B-scan图像，本实施例并不进行限制，本领域技术人员可以根据不同需求\不同场景而自行选择、设计，可以在此处使用的这些选择和设计都没有背离本发明的精神和保护范围。
作为示例，本实施例中所采用的扫频OCT系统可参见图2所示，该系统可以包括扫频光源、偏振控制器、耦合器、衰减器、振镜、平衡探测器、数据采集卡、电脑、模拟输出卡等。该系统采用的是200k的扫频光源，采集人眼视网膜数据，血管造影数据范围为视网膜中央凹区域3mm×3mm。
作为示例，在本实施例或本发明其他某些实施例中，根据所述原始数据获取B-scan图像，具体可以包括：
对采集到的所述原始数据进行指定操作以获取所述B-scan图像，其中所述指定操作可以包括以下一种或多种操作：色散补偿、加窗、傅里叶变换、取阈值。
作为示例，在本实施例或本发明其他某些实施例中，对每组B-scan图像进行配准，具体可以包括：
根据傅里叶梅林变换计算同一组内B-scan图像的偏移，根据所述偏移对所述同一组B-scan图像进行配准。
通过配准可以以减小人眼抖动产生的影响。
步骤S102，对于每组B-scan图像：根据组内在同一坐标点上的最大图像值与最小图像值的差值的平方，计算所述坐标点的血流信息，以及，根据各坐标点的血流信息获取当前组对应的B-scan血流信号图。
作为示例，在本实施例或本发明其他某些实施例中，根据组内在同一坐标点上的最大图像值与最小图像值的差值的平方，计算所述坐标点的血流信息，具体可以包括：
根据
D(x,z)＝(MAX{I1(x,z),...,In(x,z)}-MIN{I1(x,z),...,In(x,z)})2
计算坐标点(x,z)的血流信息D(x,z)，其中n为每组B-scan图像中的图像数量，I(x,z)为坐标点(x,z)位置处的图像值。
血流信息具体可以为血流信号强度。图像值具体可以为图像强度值，例如灰度值。
若有M组B-scan图像，则通过本步骤处理后可以得到连续的M张B-scan血流信号图。
步骤S103，将每组B-scan图像合成一张高清图像，对合成后的高清图像进行分层。
对于具体如何进行分层，本实施例并不进行限制，本领域技术人员可以根据不同需求\不同场景而自行选择、设计，可以在此处使用的这些选择和设计都没有背离本发明的精神和保护范围。
作为示例，在本实施例或本发明其他某些实施例中，对合成后的高清图像进行分层，具体可以包括：
计算合成后高清图像的垂直梯度，通过使用Dijikstra算法对梯度图计算最短路径并不断限制分层区域的方式获取分层线以得到分层结果。
换句话说，先计算合成后高清图像的垂直梯度，再采用Dijikstra算法对梯度图计算最短路径对视网膜进行分层，根据得到的结果不断限制分层区域，最终可以得到多条分层线。
步骤S104，根据获取的B-scan血流信号图及所述分层的结果，获取各层之间的光学相干层析血管造影OCTA。
作为示例，在本实施例或本发明其他某些实施例中，根据获取的B-scan血流信号图及所述分层的结果，获取各层之间的光学相干层析血管造影OCTA，具体可以包括：
将每张B-scan血流信号图任意层与层之间的血流信息分别在纵向上累加求和，得到血管造影图，再对所述血管造影图进行归一化处理得到对应层的OCTA。
在本实施例中，由于噪声波动比血流信号的波动要小，若计算同一坐标点位置处的图像强度组内最大值与组内最小值的差值，则对噪声部分的大小影响较小，真正的血流信号会得到很大程度的增强，使得噪声相对减弱，增加了血流信号与噪声的区分度，再对所得差值进行平方，则可进一步增加了信噪比，使血流信号更加清晰，进而使得到的血管造影图质量更高，图像中的噪点明显减少，毛细血管细节呈现更丰富，从而有助于医生对患者眼底疾病的诊断。
同时，本实施例中算法复杂度降低，计算量大大减小，减小了计算机负担，很大程度上提升了血流信号的计算速度，拍摄完毕几秒钟内就能获得高清的血管造影图像，从而也提升了工作效率。
下面再结合具体应用场景对本发明方案作进一步描述。当然以下应用场景为示例性的，在实际应用中，也可以适用于其它应用场景。
图3是根据本发明一示例性实施例示出的一种光学相干层析血管造影方法的流程图。参见图3所示，所述方法可以包括：
步骤S301，拍摄得到人眼的原始OCT数据。
例如本实施例中所采用的OCT系统可参见图2所示，该系统包括了扫频光源、偏振控制器、耦合器、衰减器、振镜、平衡探测器、数据采集卡、电脑、模拟输出卡等。该系统采用的是200k的扫频光源，采集人眼视网膜数据，血管造影数据范围为视网膜中央凹区域3mm×3mm。
步骤S302，将原始数据处理得到视网膜B-scan图像。
例如采集得到的数据经过色散补偿、加窗、傅里叶变换、取阈值等一系列操作后，可得到304组B-scan视网膜图片，每组4张(即同一个位置重复扫描的次数为4次，重复更多次可取得更好的效果)，共304×4张B-scan视网膜图片，同时，每张B-scan包含304个A-scan，即图像宽为304个像素(当然也可以采用其它尺寸)。
步骤S303，对拍摄的同一位置的B-scan图像进行配准。
例如可以采用傅里叶梅林变换，计算相同位置B-scan发生的偏移，将相同位置的4张B-scan进行配准对齐以减小人眼抖动产生的影响。
步骤S304，用配准后的图像计算血流信号。
根据
D(x,z)＝(MAX{I1(x,z),...,In(x,z)}-MIN{I1(x,z),...,In(x,z)})2
计算坐标点(x,z)位置处的血流信息D(x,z)，其中n为每组B-scan图像中的图像数量，I(x,z)为坐标点(x,z)位置处的图像强度。图像强度具体可以为灰度值。n也即对同一个位置重复扫描的次数。
同一位置的4张B-scan图像计算得到的B-scan血流信号图可以参见图4所示。经过此步骤操作后可得到304张连续的B-scan血流信号图。
步骤S305，对配准后的图像合成高清并进行分层操作。
将配准后的B-scan图像每4张合成一张高清图像，对合成后的高清图像进行分层操作。先计算合成后图片的垂直梯度，再采用Dijikstra算法对梯度图计算最短路径对视网膜进行分层。根据得到的结果不断限制分层区域，最终可得到的八条分层线，如图5所示的B-scan分层示意图，从上到下分别对应着ILM、NFL、IPL、INL、OPL、ONL、RPE和CHOROID层。
步骤S306，利用分层数据和计算得到的血流信号重构出任意层之间的OCTA图。
本步骤对步骤S304中得到的304张连续B-scan血流信号图进行进一步处理，根据步骤S305中得到的分层结果，将每张宽为304像素的B-scan血流信号图任意层与层之间的血流信息分别在纵向上进行累加求和，可以得到一张304×304大小的血管造影图，再对其进行归一化处理得到对应层的视网膜OCTA图，得到的OCTA图可以准确地反映对应层间的二维血流信息和血管分布情况。
作为示例可参见图6和图7所示。其中，图6为ILM-IPL层之间的血管造影结果，图7为IPL-OPL层之间的血管造影结果，从图中可以清晰的看到视网膜毛细血管结构，图中的噪声基本对成像效果没有影响，完成如图304×304大小的OCTA图计算时间只需要5秒。
在本实施例中，由于噪声波动比血流信号的波动要小，若计算同一坐标点位置处的最大图像值与最小图像值的差值，则对噪声部分的大小影响较小，真正的血流信号会得到很大程度的增强，使得噪声相对减弱，增加了血流信号与噪声的区分度；再对所得差值进行平方，则可进一步增加了信噪比，使血流信号更加清晰，进而使得到的血管造影图质量更高，图像中的噪点明显减少，毛细血管细节呈现更丰富，从而有助于医生对患者眼底疾病的诊断。
同时，本实施例中算法复杂度降低，计算量大大减小，减小了计算机负担，很大程度上提升了血流信号的计算速度，拍摄完毕几秒钟内就能获得高清的血管造影图像，从而也提升了工作效率。
下述为本发明装置实施例，可以用于执行本发明方法实施例。对于本发明装置实施例中未披露的细节，请参照本发明方法实施例。
图8是根据本发明一示例性实施例示出的一种光学相干层析血管造影装置的示意图。所述装置可用于由OCT设备、计算机等组成的系统。
参见图8所示，所述装置可以包括：
图像获取模块801，用于采集视网膜的光学相干层析OCT原始数据，根据所述原始数据获取B-scan图像并对每组B-scan图像进行配准，其中，对应所述视网膜同一位置的多张B-scan图像为同一组。
作为示例，在本实施例或本发明其他某些实施例中，所述图像获取模块在根据所述原始数据获取B-scan图像时，具体用于：
对采集到的所述原始数据进行指定操作以获取所述B-scan图像，其中所述指定操作包括以下一种或多种操作：色散补偿、加窗、傅里叶变换、取阈值。
作为示例，在本实施例或本发明其他某些实施例中，所述图像获取模块在对每组B-scan图像进行配准时，具体用于：
根据傅里叶梅林变换计算同一组内B-scan图像的偏移，根据所述偏移对所述同一组B-scan图像进行配准。
计算模块802，用于对每组B-scan图像根据组内在同一坐标点上的最大图像值与最小图像值的差值的平方，计算所述坐标点的血流信息，以及，根据各坐标点的血流信息获取当前组对应的B-scan血流信号图。
作为示例，在本实施例或本发明其他某些实施例中，所述计算模块在根据组内在同一坐标点上的最大图像值与最小图像值的差值的平方，计算所述坐标点的血流信息时，具体用于：
根据
D(x,z)＝(MAX{I1(x,z),...,In(x,z)}-MIN{I1(x,z),...,In(x,z)})2
计算坐标点(x,z)的血流信息D(x,z)，其中n为每组B-scan图像中的图像数量，I(x,z)为坐标点(x,z)位置处的图像值。
分层模块803，用于将每组B-scan图像合成一张高清图像，对合成后的高清图像进行分层。
作为示例，在本实施例或本发明其他某些实施例中，所述分层模块在对合成后的高清图像进行分层时，具体用于：
计算合成后高清图像的垂直梯度，通过使用Dijikstra算法对梯度图计算最短路径并不断限制分层区域的方式获取分层线以得到分层结果。
造影模块804，用于根据获取的B-scan血流信号图及所述分层的结果，获取各层之间的光学相干层析血管造影OCTA。
作为示例，在本实施例或本发明其他某些实施例中，所述造影模块具体用于：
将每张B-scan血流信号图任意层与层之间的血流信息分别在纵向上累加求和，得到血管造影图，再对所述血管造影图进行归一化处理得到对应层的OCTA。
在本实施例中，由于噪声波动比血流信号的波动要小，若计算同一坐标点位置处的最大图像值与最小图像值的差值，则对噪声部分的大小影响较小，真正的血流信号会得到很大程度的增强，使得噪声相对减弱，增加了血流信号与噪声的区分度；再对所得差值进行平方，则可进一步增加了信噪比，使血流信号更加清晰，进而使得到的血管造影图质量更高，图像中的噪点明显减少，毛细血管细节呈现更丰富，从而有助于医生对患者眼底疾病的诊断。
同时，本实施例中算法复杂度降低，计算量大大减小，减小了计算机负担，很大程度上提升了血流信号的计算速度，拍摄完毕几秒钟内就能获得高清的血管造影图像，从而也提升了工作效率。
关于上述实施例中的装置，其中各个单元\模块执行操作的具体方式已经在相关方法的实施例中进行了详细描述，此处不再赘述。
本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本发明的真正范围和精神由所附的权利要求指出。
应当理解的是，本发明并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求来限制。
Claims (12)
1.一种光学相干层析血管造影方法，其特征在于，所述方法包括：
采集视网膜的光学相干层析OCT原始数据，根据所述原始数据获取B-scan图像并对每组B-scan图像进行配准，其中，对应所述视网膜同一位置的多张B-scan图像为同一组；
对于每组B-scan图像：根据组内在同一坐标点上的最大图像值与最小图像值的差值的平方，计算所述坐标点的血流信息，以及，根据各坐标点的血流信息获取当前组对应的B-scan血流信号图；
将每组B-scan图像合成一张高清图像，对合成后的高清图像进行分层；
根据获取的B-scan血流信号图及所述分层的结果，获取各层之间的光学相干层析血管造影OCTA。
2.根据权利要求1所述的方法，其特征在于，根据所述原始数据获取B-scan图像，包括：
对采集到的所述原始数据进行指定操作以获取所述B-scan图像，其中所述指定操作包括以下一种或多种操作：色散补偿、加窗、傅里叶变换、取阈值。
3.根据权利要求1所述的方法，其特征在于，对每组B-scan图像进行配准，包括：
根据傅里叶梅林变换计算同一组内B-scan图像的偏移，根据所述偏移对所述同一组B-scan图像进行配准。
4.根据权利要求1所述的方法，其特征在于，根据组内在同一坐标点上的最大图像值与最小图像值的差值的平方，计算所述坐标点的血流信息，包括：
根据
D(x,z)＝(MAX{I1(x,z),...,In(x,z)}-MIN{I1(x,z),...,In(x,z)})2
计算坐标点(x,z)的血流信息D(x,z)，其中n为每组B-scan图像中的图像数量，I(x,z)为坐标点(x,z)位置处的图像值。
5.根据权利要求1所述的方法，其特征在于，对合成后的高清图像进行分层，包括：
计算合成后高清图像的垂直梯度，通过使用Dijikstra算法对梯度图计算最短路径并不断限制分层区域的方式获取分层线以得到分层结果。
6.根据权利要求1所述的方法，其特征在于，根据获取的B-scan血流信号图及所述分层的结果，获取各层之间的光学相干层析血管造影OCTA，包括：
将每张B-scan血流信号图任意层与层之间的血流信息分别在纵向上累加求和，得到血管造影图，再对所述血管造影图进行归一化处理得到对应层的OCTA。
7.一种光学相干层析血管造影装置，其特征在于，所述装置包括：
图像获取模块，用于采集视网膜的光学相干层析OCT原始数据，根据所述原始数据获取B-scan图像并对每组B-scan图像进行配准，其中，对应所述视网膜同一位置的多张B-scan图像为同一组；
计算模块，用于对于每组B-scan图像，根据组内在同一坐标点上的最大图像值与最小图像值的差值的平方，计算所述坐标点的血流信息，以及，根据各坐标点的血流信息获取当前组对应的B-scan血流信号图；
分层模块，用于将每组B-scan图像合成一张高清图像，对合成后的高清图像进行分层；
造影模块，用于根据获取的B-scan血流信号图及所述分层的结果，获取各层之间的光学相干层析血管造影OCTA。
8.根据权利要求7所述的装置，其特征在于，所述图像获取模块在根据所述原始数据获取B-scan图像时，具体用于：
对采集到的所述原始数据进行指定操作以获取所述B-scan图像，其中所述指定操作包括以下一种或多种操作：色散补偿、加窗、傅里叶变换、取阈值。
9.根据权利要求7所述的装置，其特征在于，所述图像获取模块在对每组B-scan图像进行配准时，具体用于：
根据傅里叶梅林变换计算同一组内B-scan图像的偏移，根据所述偏移对所述同一组B-scan图像进行配准。
10.根据权利要求7所述的装置，其特征在于，所述计算模块在根据组内在同一坐标点上的最大图像值与最小图像值的差值的平方，计算所述坐标点的血流信息时，具体用于：
根据
D(x,z)＝(MAX{I1(x,z),...,In(x,z)}-MIN{I1(x,z),...,In(x,z)})2
计算坐标点(x,z)的血流信息D(x,z)，其中n为每组B-scan图像中的图像数量，I(x,z)为坐标点(x,z)位置处的图像值。
11.根据权利要求7所述的装置，其特征在于，所述分层模块在对合成后的高清图像进行分层时，具体用于：
计算合成后高清图像的垂直梯度，通过使用Dijikstra算法对梯度图计算最短路径并不断限制分层区域的方式获取分层线以得到分层结果。
12.根据权利要求7所述的装置，其特征在于，所述造影模块具体用于：
将每张B-scan血流信号图任意层与层之间的血流信息分别在纵向上累加求和，得到血管造影图，再对所述血管造影图进行归一化处理得到对应层的OCTA。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010259880.XA CN111436905B (zh) | 2020-04-03 | 2020-04-03 | 一种光学相干层析血管造影方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010259880.XA CN111436905B (zh) | 2020-04-03 | 2020-04-03 | 一种光学相干层析血管造影方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111436905A true CN111436905A (zh) | 2020-07-24 |
| CN111436905B CN111436905B (zh) | 2024-03-19 |
Family
ID=71649896
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010259880.XA Active CN111436905B (zh) | 2020-04-03 | 2020-04-03 | 一种光学相干层析血管造影方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111436905B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112450874A (zh) * | 2020-11-20 | 2021-03-09 | 爱博图湃(北京)医疗科技有限公司 | 一种泪液分布检测方法及装置 |
| CN113643184A (zh) * | 2021-10-18 | 2021-11-12 | 广东唯仁医疗科技有限公司 | 基于光学相干断层扫描的眼底血管展示方法、系统及介质 |
| WO2023045991A1 (zh) * | 2021-09-22 | 2023-03-30 | 图湃(北京)医疗科技有限公司 | 光学相干层析血管造影方法、装置、电子设备及存储介质 |
| CN118351080A (zh) * | 2024-04-18 | 2024-07-16 | 山东探微医疗技术有限公司 | 无分叉同一毛细血管识别方法及装置 |
Citations (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1342441A (zh) * | 2000-07-21 | 2002-04-03 | Ge医疗系统环球技术有限公司 | 图像处理方法和设备、记录介质以及成像设备 |
| US20070263227A1 (en) * | 2006-05-12 | 2007-11-15 | The General Hospital Corporation | Processes, arrangements and systems for providing a fiber layer thickness map based on optical coherence tomography images |
| US20090253984A1 (en) * | 2008-04-02 | 2009-10-08 | Kabushiki Kaisha Toshiba | Medical imaging method and medical imaging system |
| JP2010213760A (ja) * | 2009-03-13 | 2010-09-30 | Toshiba Corp | 画像処理装置および画像処理プログラム |
| US20110137175A1 (en) * | 2009-12-04 | 2011-06-09 | Hossack John A | Tracked ultrasound vessel imaging |
| CN102802505A (zh) * | 2010-03-02 | 2012-11-28 | 佳能株式会社 | 图像处理设备、控制方法和光学相干断层成像系统 |
| CN102835947A (zh) * | 2012-09-04 | 2012-12-26 | 上海师范大学 | 基于散斑相关度的oct图像分析方法 |
| CN102908163A (zh) * | 2011-08-02 | 2013-02-06 | 深圳迈瑞生物医疗电子股份有限公司 | 一种帧相关方法及其装置、以及超声成像系统 |
| WO2013099638A1 (ja) * | 2011-12-28 | 2013-07-04 | 富士フイルム株式会社 | 撮像素子及び撮像装置 |
| US20130345491A1 (en) * | 2011-03-09 | 2013-12-26 | A School Corporation Kansai University | Image data processing device and transcranial magnetic stimulation apparatus |
| CN103857335A (zh) * | 2011-09-26 | 2014-06-11 | 约翰斯.霍普金斯大学 | 激光散斑图像的各向异性处理 |
| US20150213629A1 (en) * | 2012-06-29 | 2015-07-30 | Consiglio Nazionale Delle Ricerche | Method of processing optical coherence tomography images |
| CN105389810A (zh) * | 2015-10-28 | 2016-03-09 | 清华大学 | 血管内斑块的识别系统及方法 |
| US20170209037A1 (en) * | 2016-01-21 | 2017-07-27 | Canon Kabushiki Kaisha | Image display method, image display apparatus, and storage medium |
| CN107835661A (zh) * | 2015-08-05 | 2018-03-23 | 深圳迈瑞生物医疗电子股份有限公司 | 超声图像处理系统和方法及其装置、超声诊断装置 |
| CN107862724A (zh) * | 2017-12-01 | 2018-03-30 | 中国医学科学院生物医学工程研究所 | 一种改进的微血管血流成像方法 |
| CN108378824A (zh) * | 2018-03-16 | 2018-08-10 | 清华大学 | 具有阵列式环形扫描内窥探头的光学相干层析成像系统 |
| US20180308233A1 (en) * | 2015-10-30 | 2018-10-25 | Canon Kabushiki Kaisha | Radiation imaging system, information processing apparatus for irradiation image, image processing method for radiation image, and program |
| US20190038150A1 (en) * | 2016-02-04 | 2019-02-07 | The Catholic University Of Korea Industry-Academic Cooperation Foundation | Method for analyzing blood flow by using medical image |
| CN110119699A (zh) * | 2019-04-28 | 2019-08-13 | 清华大学深圳研究生院 | 内指纹提取方法、装置、系统及存储介质和指纹识别方法 |
| JP2019187798A (ja) * | 2018-04-25 | 2019-10-31 | キヤノン株式会社 | 画像処理装置およびその制御方法 |
-
2020
- 2020-04-03 CN CN202010259880.XA patent/CN111436905B/zh active Active
Patent Citations (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1342441A (zh) * | 2000-07-21 | 2002-04-03 | Ge医疗系统环球技术有限公司 | 图像处理方法和设备、记录介质以及成像设备 |
| US20070263227A1 (en) * | 2006-05-12 | 2007-11-15 | The General Hospital Corporation | Processes, arrangements and systems for providing a fiber layer thickness map based on optical coherence tomography images |
| US20090253984A1 (en) * | 2008-04-02 | 2009-10-08 | Kabushiki Kaisha Toshiba | Medical imaging method and medical imaging system |
| JP2010213760A (ja) * | 2009-03-13 | 2010-09-30 | Toshiba Corp | 画像処理装置および画像処理プログラム |
| US20110137175A1 (en) * | 2009-12-04 | 2011-06-09 | Hossack John A | Tracked ultrasound vessel imaging |
| CN102802505A (zh) * | 2010-03-02 | 2012-11-28 | 佳能株式会社 | 图像处理设备、控制方法和光学相干断层成像系统 |
| US20130345491A1 (en) * | 2011-03-09 | 2013-12-26 | A School Corporation Kansai University | Image data processing device and transcranial magnetic stimulation apparatus |
| CN102908163A (zh) * | 2011-08-02 | 2013-02-06 | 深圳迈瑞生物医疗电子股份有限公司 | 一种帧相关方法及其装置、以及超声成像系统 |
| CN103857335A (zh) * | 2011-09-26 | 2014-06-11 | 约翰斯.霍普金斯大学 | 激光散斑图像的各向异性处理 |
| WO2013099638A1 (ja) * | 2011-12-28 | 2013-07-04 | 富士フイルム株式会社 | 撮像素子及び撮像装置 |
| US20150213629A1 (en) * | 2012-06-29 | 2015-07-30 | Consiglio Nazionale Delle Ricerche | Method of processing optical coherence tomography images |
| CN102835947A (zh) * | 2012-09-04 | 2012-12-26 | 上海师范大学 | 基于散斑相关度的oct图像分析方法 |
| CN107835661A (zh) * | 2015-08-05 | 2018-03-23 | 深圳迈瑞生物医疗电子股份有限公司 | 超声图像处理系统和方法及其装置、超声诊断装置 |
| CN105389810A (zh) * | 2015-10-28 | 2016-03-09 | 清华大学 | 血管内斑块的识别系统及方法 |
| US20180308233A1 (en) * | 2015-10-30 | 2018-10-25 | Canon Kabushiki Kaisha | Radiation imaging system, information processing apparatus for irradiation image, image processing method for radiation image, and program |
| US20170209037A1 (en) * | 2016-01-21 | 2017-07-27 | Canon Kabushiki Kaisha | Image display method, image display apparatus, and storage medium |
| US20190038150A1 (en) * | 2016-02-04 | 2019-02-07 | The Catholic University Of Korea Industry-Academic Cooperation Foundation | Method for analyzing blood flow by using medical image |
| CN107862724A (zh) * | 2017-12-01 | 2018-03-30 | 中国医学科学院生物医学工程研究所 | 一种改进的微血管血流成像方法 |
| CN108378824A (zh) * | 2018-03-16 | 2018-08-10 | 清华大学 | 具有阵列式环形扫描内窥探头的光学相干层析成像系统 |
| JP2019187798A (ja) * | 2018-04-25 | 2019-10-31 | キヤノン株式会社 | 画像処理装置およびその制御方法 |
| CN110119699A (zh) * | 2019-04-28 | 2019-08-13 | 清华大学深圳研究生院 | 内指纹提取方法、装置、系统及存储介质和指纹识别方法 |
Non-Patent Citations (2)
| Title |
|---|
| 张国梅: "光学相干断层扫描血管成像检查技术" * |
| 张国梅: "光学相干断层扫描血管成像检查技术", 《实用防盲技术》, 15 August 2017 (2017-08-15) * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112450874A (zh) * | 2020-11-20 | 2021-03-09 | 爱博图湃(北京)医疗科技有限公司 | 一种泪液分布检测方法及装置 |
| CN112450874B (zh) * | 2020-11-20 | 2023-12-08 | 爱博诺德(北京)医疗科技股份有限公司 | 一种泪液分布检测方法及装置 |
| WO2023045991A1 (zh) * | 2021-09-22 | 2023-03-30 | 图湃(北京)医疗科技有限公司 | 光学相干层析血管造影方法、装置、电子设备及存储介质 |
| CN113643184A (zh) * | 2021-10-18 | 2021-11-12 | 广东唯仁医疗科技有限公司 | 基于光学相干断层扫描的眼底血管展示方法、系统及介质 |
| CN113643184B (zh) * | 2021-10-18 | 2022-02-18 | 广东唯仁医疗科技有限公司 | 基于光学相干断层扫描的眼底血管展示方法、系统及介质 |
| CN118351080A (zh) * | 2024-04-18 | 2024-07-16 | 山东探微医疗技术有限公司 | 无分叉同一毛细血管识别方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111436905B (zh) | 2024-03-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101977645B1 (ko) | 안구영상 분석방법 | |
| KR102507711B1 (ko) | 의료용 화상 처리 장치, 의료용 화상 처리 방법, 및 컴퓨터 판독가능 매체 | |
| CN111436905B (zh) | 一种光学相干层析血管造影方法及装置 | |
| US10398302B2 (en) | Enhanced vessel characterization in optical coherence tomograogphy angiography | |
| US11935241B2 (en) | Image processing apparatus, image processing method and computer-readable medium for improving image quality | |
| Hosari et al. | OCT angiography: Measurement of retinal macular microvasculature with spectralis II OCT angiography–reliability and reproducibility | |
| Schwartz et al. | Phase-variance optical coherence tomography: a technique for noninvasive angiography | |
| KR20210041046A (ko) | 의료용 화상 처리 장치, 의료용 화상 처리 방법, 컴퓨터 판독가능 매체, 및 학습 완료 모델 | |
| WO2020036182A1 (ja) | 医用画像処理装置、医用画像処理方法及びプログラム | |
| US20210224957A1 (en) | Medical image processing apparatus, medical image processing method and computer-readable medium | |
| EP1085834B1 (en) | Imaging and analyzing movement of individual erythrocytes in blood vessels | |
| Beenakker et al. | Automated retinal topographic maps measured with magnetic resonance imaging | |
| JP7362403B2 (ja) | 画像処理装置および画像処理方法 | |
| CN113557714A (zh) | 医学图像处理装置、医学图像处理方法和程序 | |
| Cheng et al. | Robust three-dimensional registration on optical coherence tomography angiography for speckle reduction and visualization | |
| CN114748032A (zh) | 一种基于oct血管成像技术的运动噪声补偿方法 | |
| JP7344847B2 (ja) | 画像処理装置、画像処理方法、及びプログラム | |
| CN111242850B (zh) | 一种广域眼底光学相干血流成像分辨率提升方法 | |
| Cagini et al. | Improvement of retinal OCT angiograms by Sampling Kantorovich algorithm in the assessment of retinal and choroidal perfusion | |
| JP7446730B2 (ja) | 画像処理装置、画像処理方法及びプログラム | |
| RU2804593C1 (ru) | Способ исследования объемного микроциркуляторного кровотока в области макулы и диска зрительного нерва | |
| Li | Computational Methods for Enhancements of Optical Coherence Tomography | |
| Wang et al. | Multi-parametric Retinal Microvascular Functional Perfusion Imaging Based on Dynamic Fundus Fluorescence Angiography | |
| Massamba et al. | Evaluation of Imaging in Patients with Posterior Uveitis and Comparative Analysis of Fluorescein and Indocyanine Green Angiography with Spectralis® HRA and Optos. Sci Set J of Ophthalmology & Eye Care, 2 (1), 01-07 | |
| CN118429214A (zh) | 一种清晰光学相干断层扫描图像的合成方法及其系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Country or region after: China Address after: Room 213, Block A, Building 1, No. 29 Shengshengyuan Road, Huilongguan Town, Changping District, Beijing, 102206 (Changping Demonstration Park) Applicant after: Tupai (Beijing) Medical Technology Co.,Ltd. Address before: Room 213, Block A, Building 1, No. 29 Shengshengyuan Road, Huilongguan Town, Changping District, Beijing, 102206 (Changping Demonstration Park) Applicant before: BEIJING TOPI IMAGING TECHNOLOGY Co.,Ltd. Country or region before: China |
|
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Country or region after: China Address after: 102206 101-f401, floor 1-6, building 3, yard 9, shengshengyuan Road, Life Science Park, Changping District, Beijing Applicant after: Tupai (Beijing) Medical Technology Co.,Ltd. Address before: Room 213, Block A, Building 1, No. 29 Shengshengyuan Road, Huilongguan Town, Changping District, Beijing, 102206 (Changping Demonstration Park) Applicant before: Tupai (Beijing) Medical Technology Co.,Ltd. Country or region before: China |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |