CN101229068B - 晶状体三维成像超声扫描方法 - Google Patents
晶状体三维成像超声扫描方法 Download PDFInfo
- Publication number
- CN101229068B CN101229068B CN2008100206051A CN200810020605A CN101229068B CN 101229068 B CN101229068 B CN 101229068B CN 2008100206051 A CN2008100206051 A CN 2008100206051A CN 200810020605 A CN200810020605 A CN 200810020605A CN 101229068 B CN101229068 B CN 101229068B
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- ultrasonic probe
- echo
- crystalline lens
- signal
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明涉及一种人眼晶状体三维成像超声扫描方法。通过数据处理单元对由多个超声探头规则地设置在探头支架上组成的超声探头阵列进行扫描控制,逐个进行发射信号扫描,并对超声探头阵列中其余的超声探头上的回波信号进行处理得到晶状体的三维图像。本晶状体三维成像超声扫描方法不仅在于能得到三维立体图像,最大优点还在于:由于超声探头阵列对人眼进行超声扫描时,无需对受检者眼球进行麻醉,利用人眼在正常视物时,两眼的晶状体是同时做相同调节的原理,进行实时、全范围、直接地观测晶状体前后曲面屈光度调节变化,观测睫状肌如何调节晶状体的功能,从而为研究眼睛的调节机理提供了有效的技术手段。
Description
技术领域
本发明涉及一种人眼晶状体三维成像超声扫描方法。
背景技术
人眼是人身体中最重要的感觉器官,非常完善、精巧和不可思议,是生命长期进化到高级形式的必然产物。人眼的构造相当于一架摄像机或照相机。前面由角膜、晶状体、前房、后房、玻璃体所共同组成的具备镜头功能的组合,把物体发出的光线聚焦到后面的用于检测光线的视网膜上。在人感觉的外界信息中,有85%以上是通过眼睛获得的。晶状体是眼球中重要的屈光介质之一。它呈双凸透镜状,前面的曲率半径约10mm,后面的约6mm,富有弹性。晶状体的直径约9mm,厚约4~5mm,对光线有屈光作用,同时也能滤去一部分紫外线,保护视网膜,但它最重要的作用是通过睫状肌的收缩或松弛改变屈光度,使看远或看近时眼球聚光的焦点都能准确地落在视网膜上。眼睛的调节是指人眼在看近距离物体时晶状体的屈光力增强,使近处的发散光线能聚焦在视网膜上,从而看清近处物体。眼睛的调节功能需要靠角膜、晶状体与睫状体共同完成,而晶状体和睫状体的协调作用是最主要的,不过调节的关键细节至今并未完全弄清。有关调节机制的研究可以追溯到19世纪中后叶。几百年来人们进行了大量研究,提出了不同的假说和理论。影响较大的主要有经典的Helmholtz松弛理论,即调节时晶状体悬韧带处于松弛状态,晶状体依赖自身的弹性使前后极变凸,但该理论和实际观察到的许多调节时的晶状体改变不相符;Tscherning紧张理论则认为调节时晶状体悬韧带至少一部分是处于紧张状态;近年来Schachar也提出了自己的调节理论。调节机制之所以一直是研究争论的热点,主要是因为睫状体区位于虹膜后,是常规眼科检查的盲区,无法直接观察其在活体眼调节中的变化。目前,应用于眼科晶状体研究的测量手段,主要有:A型超声生物测量、B型超声生物测量、超声生物显微镜(ultrasound biomicroscopy,UBM)测量、相干光断层成像技术(optical coherence tomography,OCT)、眼前节图像分析系统、Orbscan II系统等影像技术。其中A型超声是时间振幅扫描。它的特点是:超声探头以固定位置和方向对人体扫查,显示器的纵坐标代表反射回声的强度,以波幅的高低表示,回声强者波峰高、回声弱者波峰低;横坐标代表回声声源的距离和深度。波峰的形状代表界面状况。该技术目前广泛用于眼球的角膜厚度、前房深度、晶状体厚度、玻璃体腔长度、眼轴长度测量,以及预测人工晶状体的屈光度。但是缺点也很明显,获得是一维图像,仅显示波峰,直观性差。B型超声是超声声速的扇形扫描或电子开关线性扫描形成的二维切面图像,能显示器官和病变的大小、形态,以及与周围组织的关系,根据回声高低可间接判断组织性质。但测量的精确度稍差,分辨率也不够理想。UBM使用30-50M的高频率探头,其分辨率可达到50um左右,但组织穿透力低,仅为4~5mm,因此仅适用于眼前节结构的放大显示和测量。OCT作为一种新型的成像技术,能对活体组织内部微小结构进行实时、在体、高分辨力的断层成像。但是和UBM相似,两者均使用人工方法测量前房直径、辨别房角隐窝的位置,会引起主观误差。
这些技术,比如A超测量,虽然可以实现无创、无痛,但是分辨率不高,且只有一维的图像信息;有些具有较高分辨率和实时性的,如UBM,则需要借助麻醉剂,可能会影响眼睛的神经系统,且人眼在麻醉情况下,无法进行正常的调节,并不能为解决晶状体调节机制的争论提供有效的技术支持,另外,无法三维立体成像也是其存在的不足之处。而OCT、眼前节图像分析系统、Orbscan II系统等技术能得到晶状体的一些参数,但是无法实时地观测到晶状体的调节情况。
发明内容
本发明的目的在于为克服现有的眼科晶状体研究的测量手段存在着的虽然可以测量、观察晶状体,但无法进行实时、全范围、直接观察,也无法绘制晶状体在调节过程中的三维图像,只能依靠有限的数据进行图像重建的缺陷而提供一种晶状体三维成像超声扫描方法,实现人眼正常视物的情况下,无创、无痛地多方位实时观测晶状体前后曲面屈光度调节变化,观测睫状肌如何调节晶状体的功能。
为实现上述目的,本发明公开了一种晶状体三维成像超声扫描方法,通过数据处理单元对由多个超声探头规则地设置在探头支架上组成的超声探头阵列进行扫描控制,由超声探头阵列中的超声探头逐个进行发射信号扫描,并对超声探头阵列中其余的超声探头上的回波信号进行处理得到晶状体的三维图像。
扫描方式可以有多种,其中之一可以是包括如下步骤:
1)由一个超声探头发射信号,其余探头同时接收回波信号;
2)数据处理单元读取其余探头的回波信号,确定其中较强的两个回波信号,时间靠前的定为距离发射探头较近处的前表面界面的反射波,靠后的定为距离发射探头较远处的后表面界面的反射波,根据两个回波信号的时间,确定晶状体附近反射界面的一对前后表面距离探头的距离ΔX,ΔX′,从而确定出晶状体前后表面附近的反射点一对点;
3)更换一个超声探头发射信号,重复步骤1)和2),直至超声探头阵列上所有的超声探头依次轮流发射信号一次;
4)数据处理单元根据得到的点的集合,得出晶状体的整个三维轮廓。
除了上述方式,还可以采取如下扫描方式,具体包括如下步骤:
1)由超声探头阵列上的一个超声探头发射信号,由超声探头阵列上其余的超声探头中的一个接收信号;
2)数据处理单元读取该超声探头的回波信号;
3)由超声探头阵列上其余的超声探头中的另一个超声探头接收信号,重复步骤1)和2),直至超声探头阵列上其余所有的超声探头均接收信号一次;
4)数据处理单元读取其余探头的回波信号,确定其中最强的两个回波信号,时间靠前的定为前表面反射波,靠后的定为后表面反射波,根据两个回波信号的时间,确定晶状体附近反射界面的一对前后表面距离探头的距离ΔX,ΔX′,从而确定出晶状体附近的表面一对反射点;
5)更换另一个超声探头发射信号,重复步骤1)至4),直至超声探头阵列上所有的超声探头依次轮流发射信号一次;
6)数据处理单元根据得到的点的集合,得出晶状体的整个三维轮廓。
扫描专用超声探头阵列包括探头支架及规则地设置在探头支架上的多个超声探头,支架为平面板状,所述超声探头方向一致,以均匀间隔设置在支架上,所有超声探头处在同一平面上;也可以将支架设置成为球面弧形状,所述超声探头均朝向球面中心以均匀间隔设置在支架上,所有超声探头处在同一球面上。
本晶状体三维成像超声扫描方法与现有技术相比,其有益效果不仅在于能得到三维立体图像,最大优点还在于:由于超声探头阵列对人眼进行超声扫描时,无需对受检者眼球进行麻醉,利用人眼在正常视物时,两眼的晶状体是同时做相同调节的原理,进行实时、全范围、直接观测晶状体前后曲面屈光度调节变化,观测睫状肌如何调节晶状体的功能,从而为研究眼睛的调节机理提供了有效的技术手段。
下面将结合附图及具体实施例对本发明作进一步说明。
附图说明
附图1为平面探头阵列超声波的传播示意图。
附图2为球面弧形探头阵列超声波的传播示意图。
附图3为平面探头阵列结构示意图。
附图4为球面弧形探头阵列结构示意图。
附图5为球面弧形探头结构剖视图。
附图6为数据处理单元结构框图。
具体实施方式
本晶状体三维成像超声扫描方法的总体实质是通过数据处理单元对由多个超声探头规则地设置在探头支架上组成的超声探头阵列进行扫描控制,由超声探头阵列中的超声探头逐个进行发射信号扫描,并对超声探头阵列中其余的超声探头上的回波信号进行处理得到晶状体的三维图像。
具体扫描方式可以有多种,其中之一可以是包括如下步骤:
1)由一个超声探头发射信号,其余探头同时接收回波信号;
2)数据处理单元读取其余探头的回波信号,确定其中最强的两个回波信号,时间靠前的定为前表面反射波,靠后的定为后表面反射波,根据两个回波信号的时间,确定晶状体附近反射界面的一对前后表面距离探头的距离ΔX,ΔX′,从而确定出晶状体附近的表面一对反射点;
3)更换一个超声探头发射信号,重复步骤1)和2),直至超声探头阵列上所有的超声探头依次轮流发射信号一次;
4)数据处理单元根据得到的点的集合,得出晶状体的整个三维轮廓。
即假如超声探头阵列中设有n个探头,则每一次一个探头发射,其余n-1个探头同时接收回波信号。从该n-1个强弱不同的回波信号中选取最强的两点,即峰值点,时间靠前的定为前表面反射波,靠后的定为后表面反射波。每一个探头发射,可以得到一组回波信号,每组信号可以确定晶状体的一对前后表面距离探头的距离ΔX1,ΔX2,即可以绘出一对点。整个探头阵列一轮发射完毕,得到n对晶状体前后表面点的集合。多轮发射过后,即可以根据得到的点的集合,准确的绘出晶状体的整个三维轮廓。
另外还可以采取如下扫描方式,具体包括如下步骤:
1)由超声探头阵列上的一个超声探头发射信号,由超声探头阵列上其余的超声探头中的一个接收信号;
2)数据处理单元读取该超声探头的回波信号;
3)由超声探头阵列上其余的超声探头中的另一个超声探头接收信号,重复步骤1)和2),直至超声探头阵列上其余所有的超声探头均接收信号一次;
4)数据处理单元读取其余探头的回波信号,确定其中最强的两个回波信号,时间靠前的定为前表面反射波,靠后的定为后表面反射波,根据两个回波信号的时间,确定晶状体附近反射界面的一对前后表面距离探头的距离ΔX,ΔX′,从而确定出晶状体附近的表面一对反射点;
5)更换另一个超声探头发射信号,重复步骤1)至4),直至超声探头阵列上所有的超声探头依次轮流发射信号一次;
6)数据处理单元根据得到的点的集合,得出晶状体的整个三维轮廓。
该方式即是由同一个探头多次发射,不同探头分时进行接收。假如超声探头阵列中设有n个探头,在t1时刻,1号探头发射,2号探头接收;t2时刻,1号探头发射,3号探头接收;……发射次数与接受次数相同。在完成n次发射后,接收到n个回波信号,根据方法一所述之原理确定晶状体前后表面距离探头的距离ΔX1,ΔX2,即可以绘出一对点。整个探头阵列一轮发射完毕,可以得到n2对晶状体前后表面点的集合。多轮发射过后,即可以根据得到的点集合,准确的绘出晶状体的整个三维轮廓。
扫描专用超声探头阵列包括探头支架1及规则地设置在探头支架1上的多个超声探头2,如图3所示,支架1为平面板状,所述超声探头2方向一致,以均匀间隔设置在支架1上,所有超声探头2处在同一平面上;也可以将支架1设置成为规则曲面,如为球面的一部分,如图4、5所示,所述超声探头2均朝向球面中心以均匀间隔设置在支架1上,所有超声探头2处在同一球面上。探头支架1上超声探头2的个数越多,绘出的晶状体三维图像精度就越高,但对超声探头2及超声探头2阵列制作的要求更高,根据目前的实际情况,本发明最好采用支架1为球面状,如图4、5所示,曲率直径为26mm-30mm,所得截面直径为18mm-22mm,均采用直径为1mm的眼科专用B型超声探头2,超声探头2均朝向球面中心以均匀间隔设置在支架1上,所有超声探头处在同一球面上,相邻超声探头2之间的间距为0.3-0.6mm,探头的发射频率为:20MHz-30MHz。
如图6所示,本发明中的数据处理单元包括超声驱动单元、矩阵控制单元、回波信号通道、模数转换单元及中央处理单元,回波信号经回波信号通道和模数转换单元连接至中央处理单元数据输入端,中央处理单元数据输出端与矩阵控制单元连接,矩阵控制单元输出与超声驱动单元连接,超声驱动单元驱动超声探头2阵列。
晶状体前后表面到探头距离的确定:
1.平面探头阵列
平面探头阵列超声波信号的发射和接收如图1所示,超声探头2垂直发射,1号接收探头接收到晶状体前表面反射波信号,经过的时间为t1;2号接收探头接收到晶状体后表面反射波信号,经过的时间为t2。则可以得到:
前表面回波信号:d1+d2=t1*v (1)
后表面回波信号:l1+l2=t2*v (2)
已知发射探头和接受探头之间的距离为d3和l3,且探头为垂直发射,根据勾股定理,则有:
d12+d32=d22 (3)
l12+l32=l22 (4)
2.球面弧形探头阵列
球面弧形探头阵列超声波信号的发射和接收如图2所示,发射探头和接受探头之间的距离为d3和l3,探头垂直发射,且可以得到发射探头和1号接收探头之间的夹角为θ,与2号接收探头之间的夹角为φ,那么根据余弦定理有:
d22=d12+d32-2*d1*d2*cosθ (5)
l22=l12+l32-2*l1*l2*cosφ (6)
由式(1)和式(3),可以求得平面探头阵列方式下的距离d1,d2,再根据已知的探头阵列坐标,即可获得晶状体前表面的反射点坐标;由式(2)和式(4),可以求得方式下的距离l1,l2,再根据已知的探头阵列坐标,即可获得晶状体后表面的反射点坐标。
同样,由式(1)和式(5),可以求得球面弧形探头阵列方式下的距离d1,d2,再根据已知的探头阵列坐标,即可获得晶状体前表面的反射点坐标;由式(2)和式(6),可以求得球面弧形探头阵列方式下的距离l1,l2,再根据已知的探头阵列坐标,即可获得晶状体后表面的反射点坐标。
若发射探头处于阵列边缘,虽然无法接收到典型的回波信号,但由于一个探头发射,其余探头均进行接受,在发射探头处可接收到较强的界面回波信号,同样可以获得晶状体前后表面的反射点坐标。
利用具有数据处理能力的设备如计算机根据上述算式,可很方便地计算出晶状体前后表面的反射点坐标,此处不再赘述。
Claims (4)
1.一种晶状体三维成像超声扫描方法,其特征在于通过由多个规则排列的超声探头组成的超声探头阵列对晶状体的扫描并由数据处理单元对数据后续处理得到晶状体三维成像,所述扫描包括如下步骤:
1)由一个超声探头发射信号,其余探头同时接收回波信号;
2)数据处理单元读取其余探头的回波信号,确定其中最强的两个回波信号,时间靠前的定为前表面反射波,靠后的定为后表面反射波,根据两个回波信号的时间,确定晶状体附近反射界面的一对前后表面距离探头的距离ΔX,ΔX′,从而确定出晶状体附近的表面一对反射点;
3)更换一个超声探头发射信号,重复步骤1)和2),直至超声探头阵列上所有的超声探头依次轮流发射信号一次;
4)数据处理单元根据得到的点的集合,得出晶状体的整个三维轮廓。
2.根据权利要求1所述晶状体三维成像超声扫描方法,其特征在于:所述数据处理单元包括超声驱动单元、矩阵控制单元、回波信号通道、模数转换单元及中央处理单元,回波信号经回波信号通道和模数转换单元连接至中央处理单元数据输入端,中央处理单元数据输出端与矩阵控制单元连接,矩阵控制单元输出与超声驱动单元连接,超声驱动单元驱动超声探头阵列。
3.一种晶状体三维成像超声扫描方法,其特征在于通过由多个规则排列的超声探头组成的超声探头阵列对晶状体的扫描并由数据处理单元对数据后续处理得到晶状体三维成像,所述扫描包括如下步骤:
1)由超声探头阵列上的一个超声探头发射信号,由超声探头阵列上其余的超声探头中的一个接收信号;
2)数据处理单元读取该超声探头的回波信号;
3)由超声探头阵列上其余的超声探头中的另一个超声探头接收信号,重复步骤1)和2),直至超声探头阵列上其余所有的超声探头均接收信号一次;
4)数据处理单元读取其余探头的回波信号,确定其中最强的两个回波信号,时间靠前的定为前表面反射波,靠后的定为后表面反射波,根据两个回波信号的时间,确定晶状体附近反射界面的一对前后表面距离探头的距离ΔX,ΔX′,从而确定出晶状体附近的表面一对反射点;
5)更换另一个超声探头发射信号,重复步骤1)至4),直至超声探头阵列上所有的超声探头依次轮流发射信号一次;
6)数据处理单元根据得到的点的集合,得出晶状体的整个三维轮廓。
4.根据权利要求3所述晶状体三维成像超声扫描方法,其特征在于:所述数据处理单元包括超声驱动单元、矩阵控制单元、回波信号通道、模数转换单元及中央处理单元,回波信号经回波信号通道和模数转换单元连接至中央处理单元数据输入端,中央处理单元数据输出端与矩阵控制单元连接,矩阵控制单元输出与超声驱动单元连接,超声驱动单元驱动超声探头阵列。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100206051A CN101229068B (zh) | 2008-02-06 | 2008-02-06 | 晶状体三维成像超声扫描方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100206051A CN101229068B (zh) | 2008-02-06 | 2008-02-06 | 晶状体三维成像超声扫描方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101229068A CN101229068A (zh) | 2008-07-30 |
CN101229068B true CN101229068B (zh) | 2010-06-02 |
Family
ID=39896101
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100206051A Expired - Fee Related CN101229068B (zh) | 2008-02-06 | 2008-02-06 | 晶状体三维成像超声扫描方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101229068B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102283680B (zh) * | 2011-05-18 | 2013-03-06 | 中国医学科学院生物医学工程研究所 | 一种白内障晶状体硬度超声测量系统 |
CN102426194A (zh) * | 2011-11-15 | 2012-04-25 | 北京理工大学 | 复杂曲面微缺陷的阵列超声检测技术 |
CN103006274A (zh) * | 2012-12-21 | 2013-04-03 | 深圳大学 | 一种超声检测角膜粘弹性的方法和系统 |
CN105997152A (zh) * | 2016-06-13 | 2016-10-12 | 杭州融超科技有限公司 | 一体式瞳孔测量装置、使用该装置的数据处理方法和系统 |
CN106568844B (zh) * | 2016-11-14 | 2019-09-17 | 核动力运行研究所 | 一种分布式超声阵列全捕捉检测方法 |
CN114027883B (zh) * | 2021-11-08 | 2023-06-23 | 中山大学中山眼科中心 | 一种晶状体生物参数的测量方法、装置及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6949071B1 (en) * | 1998-01-12 | 2005-09-27 | Centre National De La Recherche Scientifique | Method for exploring and displaying tissue of human or animal origin from a high frequency ultrasound probe |
US7048690B2 (en) * | 2001-03-20 | 2006-05-23 | Cornell Research Foundation, Inc. | Precision ultrasound measurement for intraocular lens placement |
CN1257698C (zh) * | 2000-03-24 | 2006-05-31 | 博士伦公司 | 测量病人眼睛轴向距离的眼睛测量系统和方法 |
-
2008
- 2008-02-06 CN CN2008100206051A patent/CN101229068B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6949071B1 (en) * | 1998-01-12 | 2005-09-27 | Centre National De La Recherche Scientifique | Method for exploring and displaying tissue of human or animal origin from a high frequency ultrasound probe |
CN1257698C (zh) * | 2000-03-24 | 2006-05-31 | 博士伦公司 | 测量病人眼睛轴向距离的眼睛测量系统和方法 |
US7048690B2 (en) * | 2001-03-20 | 2006-05-23 | Cornell Research Foundation, Inc. | Precision ultrasound measurement for intraocular lens placement |
Non-Patent Citations (2)
Title |
---|
JP平1-256939A 1989.10.13 |
JP特开平6-339464A 1994.12.13 |
Also Published As
Publication number | Publication date |
---|---|
CN101229068A (zh) | 2008-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101229068B (zh) | 晶状体三维成像超声扫描方法 | |
US8496588B2 (en) | Procedures for an ultrasonic arc scanning apparatus | |
Silverman | Focused ultrasound in ophthalmology | |
US9039623B2 (en) | Compound scanning head for an ultrasonic scanning apparatus | |
Schachar et al. | In vivo increase of the human lens equatorial diameter during accommodation | |
US8894207B2 (en) | Enhanced biometry using optical coherence tomography | |
US8317709B2 (en) | Alignment and imaging of an eye with an ultrasonic scanner | |
US20130144171A1 (en) | Alignment and imaging of an eye with an ultrasonic scanner | |
US20070121120A1 (en) | Apparatus and method for measuring scleral curvature and velocity of tissues of the eye | |
CN106963338A (zh) | 晶状体生物力学和光学特性无创在体成像系统及测量方法 | |
JPH09509334A (ja) | 眼圧および血圧を監視する装置 | |
CN103079475B (zh) | 利用凹面形态阵列的超声波图像生成方法 | |
KR102312209B1 (ko) | 초음파 음향 방사력을 이용한 안압 측정 장치 | |
CN105496349B (zh) | 哈特曼人眼色差测量系统 | |
Zhang et al. | Quantitative optical coherence elastography of the optic nerve head in vivo | |
CN109009229A (zh) | 一种巩膜力学特性在体定量测量的设备及方法 | |
CN114209275A (zh) | 一种兼容oct的光声传感器 | |
Urs | Investigation of accommodation and presbyopia using ultrasound imaging during ex vivo simulated accommodation | |
KR102470241B1 (ko) | 백내장 진단을 위한 초음파 검사 장치 | |
Thijssen | Echo-ophthalmology: physical principles and diagnostic value | |
RU2424770C1 (ru) | Способ ультразвуковой диагностики косоглазия у детей | |
CN117064326A (zh) | 一种基于oct主体的双探头双目同步眼前节成像系统 | |
CN114532974A (zh) | 一种基于睫状体的人眼实物焦距个体实测探测方法、装置 | |
RU2407424C1 (ru) | Способ ранней доклинической диагностики изменения положения в глазу искусственного хрусталика глаза | |
US20130197400A1 (en) | System and method for determining the optimum size of an ultrasonic wave generation device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100602 Termination date: 20150206 |
|
EXPY | Termination of patent right or utility model |