CN104504665B - 一种基于特征点的义齿修复表面快速变形方法 - Google Patents
一种基于特征点的义齿修复表面快速变形方法 Download PDFInfo
- Publication number
- CN104504665B CN104504665B CN201410841521.XA CN201410841521A CN104504665B CN 104504665 B CN104504665 B CN 104504665B CN 201410841521 A CN201410841521 A CN 201410841521A CN 104504665 B CN104504665 B CN 104504665B
- Authority
- CN
- China
- Prior art keywords
- deformation
- point
- subpoint
- summit
- characteristic
- 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
Abstract
本发明公开了一种基于特征点的义齿修复表面快速变形方法,包括:S1、读取待操作的修复体模型,同时输入特征点集合;S2、将特征点集合中的所有特征点投影到变形曲面,获得每个特征点对应的投影点;S3、计算每个投影点对应的变形区域中的每个顶点与该投影点之间的拓扑距离,进而筛选获得该投影点对应的所有变形点;S4、根据密度函数计算变形点的变形向量;S5、根据计算出的变形向量对变形点进行变形操作。本发明可以实现义齿修复表面的快速变形,计算量低、效率高而且准确度高,可广泛应用于义齿修复体表面的变形领域中。
Description
技术领域
本发明涉及义齿修复的计算机图形构造领域,特别是涉及一种基于特征点的义齿修复表面快速变形方法。
背景技术
牙齿修复体曲面局部变形是口腔修复CAD/CAM系统的核心技术之一,其功能是针对患者口内特定的约束条件修改修复体的表面,该变形方法的效率直接影响到修复体的精确度。目前主要的变形方法有以下几种:一、基于细分变形空间的局部变形方法,该方法首先构造修复体模型的细分变形空间后,将数据映射到空间中,然后通过细分变形空间间接改变修复体的局部形态,该方法细分空间的界定并不客观,变形范围也难以确定,导致变形效果比较不稳定,而且这种方法计算量也比较大,效率比较低;二、基于B样条曲线插值修复体模型的变形方法,这种方法通过B样条曲线差值修复体模型后,通过改变B样条的控制点实现局部变形,该方法计算量较大,导致实时性较差,不适用于实时性要求较高的系统;三、基于曲线驱动的局部变形方法,这种方法的计算量相对于前两种方法会较低,效率有所提升,但是提升的效率不足以满足良好的交互性要求,该方法对处于凹凸处的变形点的变形向量的计算存在较大偏差,导致整个变形方法精确度较差。总的来说,目前的变形方法中,普遍存在计算量大、效率低以及准确度低的问题。
发明内容
为了解决上述技术问题,本发明提供一种基于特征点的义齿修复表面快速变形方法。
本发明解决其技术问题所采用的技术方案是:
一种基于特征点的义齿修复表面快速变形方法,包括:
S1、读取待操作的修复体模型,同时输入特征点集合;
S2、将特征点集合中的所有特征点投影到变形曲面,获得每个特征点对应的投影点;
S3、计算每个投影点对应的变形区域中的每个顶点与该投影点之间的拓扑距离,进而筛选获得该投影点对应的所有变形点;
S4、根据密度函数计算变形点的变形向量;
S5、根据计算出的变形向量对变形点进行变形操作。
进一步,所述步骤S2,包括:
S21、读取特征点集合中的第一个特征点;
S22、遍历修复体模型上的所有顶点后,获得与第一个特征点的距离平方值最小的顶点作为第一个特征点的投影点;
S23、依次读取特征点集合中的其它特征点,以第一个特征点的投影点为中心点,依次遍历中心点的N环邻域点,并将这些顶点中与读取的特征点之间的距离平方值最小的顶点作为该特征点的投影点;
其中,N=1,2,3,…,N0;N0为预设的正整数。
进一步,所述步骤S3,包括:
S31、在每个投影点对应的变形区域中,以该投影点为中心,依次遍历该投影点的N环邻域点,并将所遍历的每个顶点的路径信息均存储到FIFO中;
S32、依次从FIFO中读取一个顶点的路径信息后,计算该顶点到其对应的投影点的拓扑距离,并判断该拓扑距离是否小于或等于变形半径,若是,将该顶点作为投影点的变形点,否则,将该顶点剔除;
S33、重复执行步骤S32,直到FIFO读取完毕后,获得该投影点对应的所有变形点;
其中,所述变形半径是指预设的投影点对应的变形区域的半径,N=1,2,3,…,N0;N0为预设的正整数。
进一步,所述FIFO的每个数据单元包括四个参数:起点、终点、投影点与终点的最短路径值以及起点到终点的长度。
进一步,所述步骤S4,包括:
S41、对每个变形点,根据下式计算其密度函数:
F(x)=(x2-1)2(x∈[0,1])
上式中,F(x)表示变形点的密度函数,x=L/R,L表示该变形点到其对应的投影点的拓扑距离,R表示变形半径;
S42、根据下式计算变形点的变形向量:
Vec(pt)=F(x)×{Cnew(vt)-Cold(vt)}
上式中,Vec(pt)表示变形点的变形向量,Cnew(vt)表示该变形点对应的投影点,Cold(vt)表示该投影点对应的特征点。
进一步,所述步骤S5,其具体为:
根据计算出的变形向量,结合下式对变形点进行变形操作:
pt'=pt+Vec(pt)
上式中,pt代表变形点,pt'代表对变形点进行变形操作后得到的顶点,Vec(pt)代表变形点的变形向量。
本发明的有益效果是:本发明提供了一种基于特征点的义齿修复表面快速变形方法,包括:S1、读取待操作的修复体模型,同时输入特征点集合;S2、将特征点集合中的所有特征点投影到变形曲面,获得每个特征点对应的投影点;S3、计算每个投影点对应的变形区域中的每个顶点与该投影点之间的拓扑距离,进而筛选获得该投影点对应的所有变形点;S4、根据密度函数计算变形点的变形向量;S5、根据计算出的变形向量对变形点进行变形操作。本发明直接获取用户输入的特征点集合,没有拟合等数据计算过程,而且通过将特征点投影到变形曲面后,根据密度函数计算出变形点的变形向量,从而对变形点进行变形操作,实现义齿修复表面的快速变形,计算量低、效率高而且准确度高。
而且,本方法基于顶点的N环邻域来获得特征点的投影点以及投影点对应的变形区域中的变形点,大大地降低了计算量,显著地提高了义齿修复表面变形的效率,且准确度较高。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是本发明的一种基于特征点的义齿修复表面快速变形方法的示意图;
图2是N环邻域示意图;
图3是本发明的一具体实施例中查找特征点的投影点的原理示意图;
图4是本发明的一具体实施例中的一投影点与其N环邻域点之间的拓扑结构图;
图5是本发明的一具体实施例中的FIFO模型的示意图。
具体实施方式
为了便于下文的描述,首先给出以下术语解释:
N环邻域点:拓扑模型上,计算顶点B到顶点A所需经过的最短边数,若这个最短边数的值为1,则顶点B是顶点A的1环邻域点,如果这个最短边数的值为N,则B是A的N环邻域点。换言之,顶点的N环邻域点是指与顶点的最短拓扑距离为N的顶点,N为不小于1的自然数。例如,图2中,顶点B是顶点A的4环邻域点,实际上,顶点B到顶点A的拓扑距离有多个,根据N环邻域点的定义选取最短的拓扑距离。
FIFO:First Input First Output的缩写,先入先出队列。
参照图1,本发明提供了一种基于特征点的义齿修复表面快速变形方法,包括:
S1、读取待操作的修复体模型,同时输入特征点集合;
S2、将特征点集合中的所有特征点投影到变形曲面,获得每个特征点对应的投影点;
S3、计算每个投影点对应的变形区域中的每个顶点与该投影点之间的拓扑距离,进而筛选获得该投影点对应的所有变形点;
S4、根据密度函数计算变形点的变形向量;
S5、根据计算出的变形向量对变形点进行变形操作。
进一步作为优选的实施方式,所述步骤S2,包括:
S21、读取特征点集合中的第一个特征点;
S22、遍历修复体模型上的所有顶点后,获得与第一个特征点的距离平方值最小的顶点作为第一个特征点的投影点;
S23、依次读取特征点集合中的其它特征点,以第一个特征点的投影点为中心点,依次遍历中心点的N环邻域点,并将这些顶点中与读取的特征点之间的距离平方值最小的顶点作为该特征点的投影点;
其中,N=1,2,3,…,N0;N0为预设的正整数。
进一步作为优选的实施方式,所述步骤S3,包括:
S31、在每个投影点对应的变形区域中,以该投影点为中心,依次遍历该投影点的N环邻域点,并将所遍历的每个顶点的路径信息均存储到FIFO中;
S32、依次从FIFO中读取一个顶点的路径信息后,计算该顶点到其对应的投影点的拓扑距离,并判断该拓扑距离是否小于或等于变形半径,若是,将该顶点作为投影点的变形点,否则,将该顶点剔除;
S33、重复执行步骤S32,直到FIFO读取完毕后,获得该投影点对应的所有变形点;
其中,所述变形半径是指预设的投影点对应的变形区域的半径,N=1,2,3,…,N0;N0为预设的正整数。
进一步作为优选的实施方式,所述FIFO的每个数据单元包括四个参数:起点、终点、投影点与终点的最短路径值以及起点到终点的长度。
进一步作为优选的实施方式,所述步骤S4,包括:
S41、对每个变形点,根据下式计算其密度函数:
F(x)=(x2-1)2(x∈[0,1])
上式中,F(x)表示变形点的密度函数,x=L/R,L表示该变形点到其对应的投影点的拓扑距离,R表示变形半径;
S42、根据下式计算变形点的变形向量:
Vec(pt)=F(x)×{Cnew(vt)-Cold(vt)}
上式中,Vec(pt)表示变形点的变形向量,Cnew(vt)表示该变形点对应的投影点,Cold(vt)表示该投影点对应的特征点。
进一步作为优选的实施方式,所述步骤S5,其具体为:
根据计算出的变形向量,结合下式对变形点进行变形操作:
pt'=pt+Vec(pt)
上式中,pt代表变形点,pt'代表对变形点进行变形操作后得到的顶点,Vec(pt)代表变形点的变形向量。
下面结合具体实施例对本发明做详细说明。
参照图1,一种基于特征点的义齿修复表面快速变形方法,包括:
S1、读取待操作的修复体模型,同时输入特征点集合;本实施例中,特征点集合来源于CAD/CAE系统使用者手动对修复体模型进行调整,或者修复体模型与邻牙适配时产生;特征点集合可以位于修复体上,也可以不在修复体上;本发明中的修复体即指义齿修复预备体。
S2、将特征点集合中的所有特征点投影到变形曲面,获得每个特征点对应的投影点;
S3、计算每个投影点对应的变形区域中的每个顶点与该投影点之间的拓扑距离,进而筛选获得该投影点对应的所有变形点;每个投影点对应的变形区域是指以该投影点为中心,以预设半径R的球体与义齿修复体模型相交的区域;
S4、根据密度函数计算变形点的变形向量;
S5、根据计算出的变形向量对变形点进行变形操作,其具体如下:
根据计算出的变形向量,结合下式对变形点进行变形操作:
pt'=pt+Vec(pt)
上式中,pt代表变形点,pt'代表对变形点进行变形操作后得到的顶点,Vec(pt)代表变形点的变形向量。
参照图3,步骤S2,具体包括S21~S23:
S21、读取特征点集合中的第一个特征点T0;
S22、遍历修复体模型上的所有顶点后,获得与第一个特征点T0的距离平方值最小的顶点P0作为第一个特征点T0的投影点;
S23、依次读取特征点集合中的其它特征点Tx,以第一个特征点T0的投影点P0为中心点,依次遍历中心点P0的N环邻域点,并将这些顶点中与读取的特征点Tx之间的距离平方值最小的顶点Px作为该特征点的投影点;Tx指T1,T2,T3,T4…;
其中,N=1,2,3,…,N0;N0为预设的正整数。
图3中,T0,T1,T2,T3,T4…均表示特征点,图3中第一圈加粗线上的顶点表示P0的1环邻域点,第二圈加粗线上的顶点表示P0的2环邻域点。
依次遍历中心点P0的N环邻域点的过程如下:依次遍历中心点P0的所有1环邻域点后,遍历中心点的2环邻域点,遍历完后再遍历中心点的3环邻域点,…最后遍历中心点的N0环邻域点。
步骤S21~S23中获得特征点的投影点的方法,只需要将第一个特征点T0与修复体模型上的所有顶点进行对比,而其他特征点T1,T2,T3,…在拓扑结构上有一定的聚集关系,所以查找到第一个特征点T0的投影点P0后只需以P0作为中心点向外遍历中心点的N环邻域搜索其他特征点对应的投影点,无需遍历修复体模型的所有顶点,大大减少了计算量,提升了获取特征点的投影点的效率。
步骤S3,具体包括S31~S33:
S31、在每个投影点对应的变形区域中,以该投影点为中心,依次遍历该投影点的N环邻域点,即依次遍历的1环邻域点,2环邻域点,…,N0环邻域点,遍历过程中并将所遍历的每个顶点的路径信息均存储到FIFO中;这里,每个顶点的路径信息指该顶点到投影点的最短路径;
S32、依次从FIFO中读取一个顶点的路径信息后,计算该顶点到其对应的投影点的拓扑距离L,并判断该拓扑距离L是否小于或等于变形半径R,若是,将该顶点作为投影点的变形点,否则,将该顶点剔除;换句话说,即剔除掉L>R的顶点;
S33、重复执行步骤S32,直到FIFO读取完毕后,获得该投影点对应的所有变形点;
其中,所述变形半径是指预设的投影点对应的变形区域的半径,N=1,2,3,…,N0;N0为预设的正整数。
本实施例中,每个投影点均独立使用一个FIFO进行数据读写。图4中所示的拓扑结构中,A为投影点,B、C、D、E、F、G分别为投影点A在其对应的变形区域中的1环邻域点或2环邻域点,执行步骤S31对图4中的顶点进行遍历后,获得的FIFO模型如图5所示,该模型中,FIFO的每个数据单元包括四个参数:起点、终点、投影点与终点的最短路径值以及起点到终点的长度,即FIFO的数据单元为一个结构体,定义为:
其中A代表边的起点,B代表边的终点,id代表投影特征点与B的最短路径值,dis代表该边的长度。
由图4和图5可以看出,寻找最短路径的过程其实就是构建FIFO的过程,拓扑距离L的计算过程其实就是基于FIFO数据单元中的dis计算得出的。
图4中顶点A和顶点G之间的拓扑距离L,参照图5中的FIFO的数据信息,计算步骤如下:
1)找出终点G,列出起点,这里只有B;
2)找出终点B,列出起点,这里有且只有A;
3)|AG|=|AB|+|BG|
在计算过程中,若有多条路径,只保留最小值。
步骤S4,具体包括S41~S42:
S41、对每个变形点,根据下式计算其密度函数:
F(x)=(x2-1)2(x∈[0,1])
上式中,F(x)表示变形点的密度函数,x=L/R,L表示该变形点到其对应的投影点的拓扑距离,R表示变形半径,所以x的取值范围为[0,1];
S42、根据下式计算变形点的变形向量:
Vec(pt)=F(x)×{Cnew(vt)-Cold(vt)}
上式中,Vec(pt)表示变形点的变形向量,Cnew(vt)表示该变形点对应的投影点,Cold(vt)表示该投影点对应的特征点。
计算获得密度函数后,首先是计算获得投影点的变形向量Vec(vt)=Cnew(vt)-Cold(vt),然后根据变形点的变形向量Vec(pt)=F(x)×Vec(vt)来获得变形点的变形向量,即得到步骤S42中所示的公式,其中,Vec(vt)表示投影点的变形向量。
本方法直接获取用户输入的特征点集合,没有拟合等数据计算过程,且基于顶点的N环邻域来获得特征点的投影点以及投影点对应的变形区域中的变形点,大大地降低了计算量,显著地提高了义齿修复表面变形的效率,且准确度较高。
而且本方法采用邻域分段计算拓扑距离,解决了单纯使用空间距离的局限性。
以上是对本发明的较佳实施进行了具体说明,但本发明创造并不限于所述实施例,熟悉本领域的技术人员在不违背本发明精神的前提下还可做出种种的等同变形或替换,这些等同的变型或替换均包含在本申请权利要求所限定的范围内。
Claims (4)
1.一种基于特征点的义齿修复表面快速变形方法,其特征在于,包括:
S1、读取待操作的修复体模型,同时输入特征点集合;
S2、将特征点集合中的所有特征点投影到变形曲面,获得每个特征点对应的投影点;
S3、计算每个投影点对应的变形区域中的每个顶点与该投影点之间的拓扑距离,进而筛选获得该投影点对应的所有变形点;
S4、根据密度函数计算变形点的变形向量;
S5、根据计算出的变形向量对变形点进行变形操作;
所述步骤S2,包括:
S21、读取特征点集合中的第一个特征点;
S22、遍历修复体模型上的所有顶点后,获得与第一个特征点的距离平方值最小的顶点作为第一个特征点的投影点;
S23、依次读取特征点集合中的其它特征点,以第一个特征点的投影点为中心点,依次遍历中心点的N环邻域点,并将这些顶点中与读取的特征点之间的距离平方值最小的顶点作为该特征点的投影点;
其中,N=1,2,3,…,N0;N0为预设的正整数;
所述步骤S3,包括:
S31、在每个投影点对应的变形区域中,以该投影点为中心,依次遍历该投影点的N环邻域点,并将所遍历的每个顶点的路径信息均存储到FIFO中;
S32、依次从FIFO中读取一个顶点的路径信息后,计算该顶点到其对应的投影点的拓扑距离,并判断该拓扑距离是否小于或等于变形半径,若是,将该顶点作为投影点的变形点,否则,将该顶点剔除;
S33、重复执行步骤S32,直到FIFO读取完毕后,获得该投影点对应的所有变形点;
其中,所述变形半径是指预设的投影点对应的变形区域的半径,N=1,2,3,…,N0;N0为预设的正整数。
2.根据权利要求1所述的一种基于特征点的义齿修复表面快速变形方法,其特征在于,所述FIFO的每个数据单元包括四个参数:起点、终点、投影点与终点 的最短路径值以及起点到终点的长度。
3.根据权利要求1所述的一种基于特征点的义齿修复表面快速变形方法,其特征在于,所述步骤S4,包括:
S41、对每个变形点,根据下式计算其密度函数:
F(x)=(x2-1)2,x∈[0,1]
上式中,F(x)表示变形点的密度函数,x=L/R,L表示该变形点到其对应的投影点的拓扑距离,R表示变形半径;
S42、根据下式计算变形点的变形向量:
Vec(pt)=F(x)·{Cnew(vt)-Cold(vt)}
上式中,Vec(pt)表示变形点的变形向量,Cnew(vt)表示该变形点对应的投影点的坐标,Cold(vt)表示该投影点对应的特征点的坐标。
4.根据权利要求3所述的一种基于特征点的义齿修复表面快速变形方法,其特征在于,所述步骤S5,其具体为:
根据计算出的变形向量,结合下式对变形点进行变形操作:
pt'=pt+Vec(pt)
上式中,pt代表变形点的坐标,pt'代表对变形点进行变形操作后得到的顶点的坐标,Vec(pt)代表变形点的变形向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410841521.XA CN104504665B (zh) | 2014-12-26 | 2014-12-26 | 一种基于特征点的义齿修复表面快速变形方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410841521.XA CN104504665B (zh) | 2014-12-26 | 2014-12-26 | 一种基于特征点的义齿修复表面快速变形方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104504665A CN104504665A (zh) | 2015-04-08 |
CN104504665B true CN104504665B (zh) | 2017-11-28 |
Family
ID=52946059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410841521.XA Active CN104504665B (zh) | 2014-12-26 | 2014-12-26 | 一种基于特征点的义齿修复表面快速变形方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104504665B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109272466A (zh) * | 2018-09-19 | 2019-01-25 | 维沃移动通信有限公司 | 一种牙齿美化方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103237519A (zh) * | 2010-10-01 | 2013-08-07 | 3形状股份有限公司 | 假牙的建模和制造 |
CN103871097A (zh) * | 2014-02-26 | 2014-06-18 | 南京航空航天大学 | 基于牙齿预备体的数据柔性融合方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8083522B2 (en) * | 2008-10-29 | 2011-12-27 | Inpronto Inc. | Method for tooth implants |
-
2014
- 2014-12-26 CN CN201410841521.XA patent/CN104504665B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103237519A (zh) * | 2010-10-01 | 2013-08-07 | 3形状股份有限公司 | 假牙的建模和制造 |
CN103871097A (zh) * | 2014-02-26 | 2014-06-18 | 南京航空航天大学 | 基于牙齿预备体的数据柔性融合方法 |
Non-Patent Citations (1)
Title |
---|
全冠和全冠桥数字化设计关键技术研究与应用;闫国栋;《中国博士学位论文全文数据库 信息科技辑》;20141215(第12期);第44页第3.1节,第46页第3.2.1节,第50页第3.2.2.3节,第50-51页第3.2.3节,第56-57页第3.3.3节 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109272466A (zh) * | 2018-09-19 | 2019-01-25 | 维沃移动通信有限公司 | 一种牙齿美化方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104504665A (zh) | 2015-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106327535B (zh) | Cbct的牙根与口内扫描的牙冠的融合方法 | |
CN103198494B (zh) | 一种基于建筑外轮廓线及激光点云的建筑立面提取方法 | |
CN104318600B (zh) | 一种使用Bézier曲线实现角色行走轨迹动画的方法 | |
CN108711194B (zh) | 一种基于三次b样条插值的三维网格模型拼接方法 | |
Zhang et al. | Automatic 3D tooth segmentation using convolutional neural networks in harmonic parameter space | |
CN101964116A (zh) | 紧凑的插中 | |
CN108629838B (zh) | 三维牙冠网格模型的孔洞修补方法及装置 | |
CN105488849B (zh) | 一种基于混合水平集的三维牙齿建模方法 | |
CN105213043B (zh) | 一种数字化牙冠贴面设计方法 | |
CN102184567A (zh) | 基于球b样条曲线的三维血管模型构造方法 | |
CN111991106B (zh) | 一种牙套切割线自动生成方法及应用 | |
CN107689254B (zh) | 一种全冠修复体外表面数字化生成方法 | |
CN110021000B (zh) | 基于图层变形的发际线修复方法及装置 | |
CN106204688A (zh) | 一种牙弓曲线生成方法 | |
CN109086495A (zh) | Bim出精装施工图的方法 | |
CN110910433A (zh) | 一种基于深度学习的点云匹配方法 | |
CN112221145A (zh) | 游戏脸模生成方法及装置、存储介质及电子设备 | |
Chen et al. | Automatic surface repairing, defeaturing and meshing algorithms based on an extended B-rep | |
CN104504665B (zh) | 一种基于特征点的义齿修复表面快速变形方法 | |
CN113689538A (zh) | 一种视频生成方法、装置、电子设备及存储介质 | |
CN113298828B (zh) | 一种基于卷积神经网络的颌骨自动分割方法 | |
CN107622529A (zh) | 一种基于形态学的三维牙齿模型自动分割方法 | |
CN105741287B (zh) | 牙齿三维网格数据的分割方法及装置 | |
CN109783965B (zh) | 一种结构网格自动分块加密方法 | |
CN103700129A (zh) | 一种任意人头与任意人体的3d结合方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |