CN107203761B - 基于高分辨率卫星影像的道路宽度估计方法 - Google Patents
基于高分辨率卫星影像的道路宽度估计方法 Download PDFInfo
- Publication number
- CN107203761B CN107203761B CN201710454604.7A CN201710454604A CN107203761B CN 107203761 B CN107203761 B CN 107203761B CN 201710454604 A CN201710454604 A CN 201710454604A CN 107203761 B CN107203761 B CN 107203761B
- Authority
- CN
- China
- Prior art keywords
- road
- width
- pixel
- road width
- descriptor
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/182—Network patterns, e.g. roads or rivers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
Abstract
本发明公开了一种基于高分辨率卫星影像的道路宽度估计方法,根据高分辨率卫星影像中道路的特征,定义了一种新颖的道路宽度描述符,通过统计固定大小的道路块内的像素分布并用直方图的形式表示,可以描述道路的宽度特征;在道路宽度描述符的基础上,结合卷积神经网络,提出了基于卷积神经网络的道路宽度估计方法;针对道路宽度估计结果中存在的宽度不连续问题,根据道路宽度的连续性,定义了一种基于空间一致性的道路宽度估计能量函数,通过最小化该能量函数可以实现卫星影像中道路宽度类别的合理估计。
Description
技术领域
本发明涉及遥感应用领域,特别涉及一种基于高分辨率卫星影像的道路宽度估计方法。
背景技术
卫星遥感技术具有高效性、实时性以及信息多元化等特点,其作为一种先进的对地观测方法在城市交通领域中发挥了重要的作用。而且随着遥感技术的发展,分辨率的极大提高,让卫星图像能够更清楚的表达地物目标的特征,为计算机从卫星图像中提取道路信息提供了必要的条件,使得基于高分辨率卫星影像的道路宽度估计成为可能。
从上个世纪七十年代开始,众多学者在道路信息提取领域进行了大量的研究工作并取得了丰硕的成果,特别是对道路的现状特征提取更是提出了很多具有实用价值的方法。但是对于道路宽度估计,相关的研究成果却并不多见。
发明内容
本发明的目的在于提供一种基于高分辨率卫星影像的道路宽度估计方法。
为实现上述目的,本发明采用以下技术方案:
基于高分辨率卫星影像的道路宽度估计方法,包括以下步骤:
S1、截取不同宽度等级的道路块样本,将道路块样本转换为道路宽度描述符,利用不同宽度等级道路块样本转换得到的宽度描述符训练卷积神经网络以得到道路宽度类别预测模型;
S2、提取待进行宽度估计的卫星图像的道路中心线,沿道路中心线进行滑窗,以道路中心线上每个像素点为中心截取固定大小的道路块并转换为道路宽度描述符,利用S1所获得的道路宽度类别预测模型对该步骤所获得的道路宽度描述符进行预测,得到道路块属于每个宽度类别的概率;
S3、根据道路宽度的连续性及空间一致性定义道路宽度估计能量函数,对道路中心线上的每个像素点寻求一个最优的宽度类别使其对应的道路宽度估计能量函数取值最小,则各像素点所对应的最优的宽度类别的宽度范围即以各像素点为中心的道路块的宽度范围;
其中,道路宽度估计能量函数记为式中,L={l1…ln},其表示以像素p为中心的道路块的宽度类别集合,像素p为位于道路中心线上的像素点;D(P)=1-P(p,li),P(p,li)为以像素p为中心的道路块属于宽度类别li的概率,P(p,li)由训练好的道路宽度类别预测模型给出;S(L)为平滑项,表示对以像素p为中心的道路块宽度进行估计时,来自相邻道路宽度类别的影响;参数ε为平衡因子,调整平滑项对宽度估计结果的影响程度。
进一步地,步骤S1具体如下:
S11、将卫星图像中的道路宽度分为若干个等级,针对每种道路宽度图像等级,截取大量属于其宽度范围的道路块样本,所截取的道路块样本大小相同;
S12、将所截取的道路块样本转换为道路宽度描述符;
S121、以截取的道路块样本的中心为圆心,构建a个同心圆,a个同心圆的半径分别为道路块内切圆半径的1/a、2/a……7/a及1倍,将a个同心圆标记为Cn(n=1,2……a);
S122、根据道路块中的像素值分布,将道路块中的像素值分为b个等长的像素值区间,记为Im(m=1,2……b),令max和min分别为道路块中最大像素值和最小像素值,则每个像素值区间的长度为(max-min)/b;统计每个圆范围内落在各像素值区间里的像素个数,将第n个圆内属于第m个区间范围内的像素个数记为Nn,m,并将每个圆范围内落在各像素值区间内的像素个数归一化到0-255范围内;
S123、将第n个圆内的像素值分布作为描述符的第n行,获得道路宽度描述符;
S13、将道路宽度描述符制作为训练样本,所述训练样本包括训练集和验证集,利用卷积神经网络对训练样本进行学习,获得道路宽度类别预测模型。
进一步地,步骤S2采用与步骤S12相同的方式将其截取的道路块转化为道路宽度描述符。
进一步地,步骤S1与S2中,还包括采用图像放大算法将所获得的道路宽度描述符放大若干倍的步骤。
进一步地,所述图像放大算法为最近邻插值算法。
进一步地,步骤S3中S(L)的计算具体为:
其中,Np表示道路中心线上与像素p相邻的若干个像素的集合;li和lq分别表示以像素p和q为中心的道路块的宽度类别;||p-q||为像素p与像素q之间的距离;当li=lq时,符号函数sign(li,lq)取值为0,否则sign(li,lq)取值为1。
采用上述技术方案后,本发明能够通过道路宽度描述对道路的宽度特征进行很好描述并对不同宽度等级的道路进行区分,在此基础之上引入卷积神经网络得到道路宽度类别预测模型,该模型可以给出道路块属于每种宽度类别的概率,最后结合空间一致性可以对道路的宽度进行合理的估计,为基于高分辨率卫星影像的道路宽度估计提供了一种新颖的解决方案。
附图说明
图1为本发明流程图;
图2为本发明的道路宽度描述符转换结果示例图,图(a)-图(c)为原道路块,图(d)-图(f)分别为图(a)-图(c)转换得到的道路宽度描述符;
图3为本发明在不同图像上进行道路宽度估计的效果图,图3(a),(c),(e)是Pleiades-1A测试图像、Gaofen2测试图像和航拍测试图像的道路宽度估计结果的局部图,其中局部放大图为宽度估计错误区域,图3(b),(d),(f)分别为对应的道路宽度类别真值图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
如图1所示的是本发明的流程图,主要包括以下三个步骤:
S1、截取不同宽度等级的道路块样本,将道路块样本转换为道路宽度描述符,利用不同宽度等级道路块样本转换得到的宽度描述符训练卷积神经网络以得到道路宽度类别预测模型。
具体为:
S11、将卫星图像中的道路宽度分为若干个等级,针对每种道路宽度图像等级,截取大量属于其宽度范围的道路块样本,所截取的道路块样本大小相同;
S12、将所截取的道路块样本转换为道路宽度描述符;
S121、以截取的道路块样本的中心为圆心,建立一个a*a大小的二维直方图来统计道路块中像素的分布。具体为构建a个同心圆(本实施例中a=8),a个同心圆的半径分别为道路块内切圆半径的1/a、2/a……7/a及1倍,将a个同心圆标记为Cn(n=1,2……a);
S122、根据道路块中的像素值分布,将道路块中的像素值分为b个等长的像素值区间,记为Im(m=1,2……b),令max和min分别为道路块中最大像素值和最小像素值,则每个像素值区间的长度为(max-min)/b,本实施例中b=8;统计每个圆范围内落在各像素值区间里的像素个数,将第n个圆内属于第m个区间范围内的像素个数记为Nn,m,并将每个圆范围内落在各像素值区间内的像素个数归一化到0-255范围内;
S123、将第n个圆内的像素值分布作为描述符的第n行(即将各个圆范围内归一化后的像素值分布作为描述符的一行,第一个圆内的像素值分布作为描述符的第一行,第二个圆内的像素值分布作为描述符的第二行,依此类推),获得道路宽度描述符;
则道路宽度描述符中第n行、第m列的像素值为
其中,Nn,m为圆Cn内属于区间Im的像素个数,K为圆Cn内的总像素个数。
S124、图像放大算法将所获得的道路宽度描述符放大若干倍。
为了便于利用卷积神经网络对这些道路宽度描述符进行训练,需要将这些道路宽度描述符从8×8放大为32×32,为了不改变道宽度描述符对道路块宽度信息的描述,所使用的图像放大算法为最近邻插值算法。
图2是道路宽度描述符转换实例。其中图(a)与图(b)为相同宽度等级的道路块,图(c)为较窄的道路块,图(d)-图(f)分别为图(a)-图(c)转换得到的道路宽度描述符。从图中可以看出,相同宽度等级的道路块转换得到的宽度描述符非常相似,而不同宽度等级的道路块的宽度描述符之间具有明显差异。
S13、将道路宽度描述符制作为训练样本,所述训练样本包括训练集和验证集,利用Cifar-10卷积神经网络对训练样本进行学习,获得道路宽度类别预测模型。
S2、提取待进行宽度估计的卫星图像的道路中心线,沿道路中心线进行滑窗,以道路中心线上每个像素点为中心截取固定大小的道路块并转换为道路宽度描述符,利用S1所获得的道路宽度类别预测模型对该步骤所获得的道路宽度描述符进行预测,得到道路块属于每个宽度类别的概率;
S3、根据道路宽度的连续性,结合空间一致性定义了一种道路宽度估计能量函数,通过最小化该能量函数可以实现对道路宽度的合理估计。
对于一个位于道路中心线C上的像素p,用L={l1…ln}表示以p像素为中心的道路块的宽度类别集合,则宽度估计能量函数可以记为:
其中,D(P)为数据项,由概率图P决定;S(L)为平滑项,由类别集合L决定。参数ε为平衡因子,该参数可以决定平滑项对整个宽度估计结果的影响程度,ε为经验值,其取值范围为[0.1,1]。
数据项中D(P)的定义如下所示:
D(P)=1-P(p,li)
其中,P(p,li)为以像素p为中心的道路块属于宽度类别li的概率,表示道路自身的宽度信息,该概率值由训练好的卷积神经网络模型给出。
平滑项S(li,lq)表示在对以像素p为中心的道路块宽度进行估计时,来自相邻道路块宽度类别的影响,定义如下:
其中,Np表示道路中心线上与像素p相邻的像素集合,是一个一维的范围,通常取为与像素p最近的12个像素。li和lq分别表示以像素p和q为中心的道路块的宽度类别。||p-q||为像素p与像素q之间的距离。当li=lq时,符号函数sign(li,lq)取值为0,否则sign(li,lq)取值为1。
目标是,对于每一个位于道路中心线C上的像素p,找到一个最优的宽度类别lp使得能量函数取值最小,这是一个经典的马尔科夫随机场(MRF)问题,可以使用图割(GraphCuts)算法解决该问题,最后道路中心线上的每个像素都获得一个道路宽度类别,表示以该像素为中心的道路块的宽度类别,该宽度类别的宽度范围即为道路块的宽度范围。
图3是在不同图像上进行宽度估计的效果图,呈现了本发明应用于高分辨率卫星影像上的道路宽度估计效果。本发明在Pleiades-1A、Gaofen2的卫星图像以及航拍图像上进行了测试,这些测试图像的分辨率均为0.5m,Pleiades-1A测试图像大小为3200×3200,Gaofen2测试图像和航拍测试图像的大小均为6400×4800。其中图3(a),(c),(e)是分别本发明对Pleiades-1A测试图像、Gaofen2测试图像和航拍测试图像的道路宽度估计结果的局部图,其中局部放大图为宽度估计错误区域,图3(b),(d),(f)分别为对应的道路宽度类别真值图。根据测试图像中所存在道路的宽度信息,将Pleiades-1A测试图像、Gaofen2测试图像和航拍测试图像中的道路分别分为3个、4个和3个宽度类别,并用不同颜色来对不同宽度类别(A、B、C及D)进行标记,则三幅测试图像中的宽度范围如下表所示:
将算法估计结果与真值图相比较可以发现,对于三种不同传感器记录的数据,算法都能够正确的估计大部分道路的宽度。经统计,算法对Pleiades-1A测试图像、Gaofen2测试图像和航拍测试图像的道路宽度估计准确率分别为85.2%、90.8%、92.1%。
以上,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (6)
1.基于高分辨率卫星影像的道路宽度估计方法,其特征在于,包括以下步骤:
S1、截取不同宽度等级的道路块样本,将道路块样本转换为道路宽度描述符,利用不同宽度等级道路块样本转换得到的宽度描述符训练卷积神经网络以得到道路宽度类别预测模型;
S2、提取待进行宽度估计的卫星图像的道路中心线,沿道路中心线进行滑窗,以道路中心线上每个像素点为中心截取固定大小的道路块并转换为道路宽度描述符,利用S1所获得的道路宽度类别预测模型对该步骤所获得的道路宽度描述符进行预测,得到道路块属于每个宽度类别的概率;
S3、根据道路宽度的连续性及空间一致性定义道路宽度估计能量函数,对道路中心线上的每个像素点寻求一个最优的宽度类别使其对应的道路宽度估计能量函数取值最小,则各像素点所对应的最优的宽度类别的宽度范围即以各像素点为中心的道路块的宽度范围;
其中,道路宽度估计能量函数记为L={l1…ln},其表示以像素p为中心的道路块的宽度类别集合,像素p为位于道路中心线上的像素点;D(P)=1-P(p,li),P(p,li)为以像素p为中心的道路块属于某一宽度类别li的概率,P(p,li)由训练好的道路宽度类别预测模型给出;S(L)为平滑项,表示对以像素p为中心的道路块宽度进行估计时,来自相邻道路宽度类别的影响;参数ε为平衡因子,调整平滑项对宽度估计结果的影响程度。
2.根据权利要求1所述的基于高分辨率卫星影像的道路宽度估计方法,其特征在于,步骤S1具体如下:
S11、将卫星图像中的道路宽度分为若干个等级,针对每种道路宽度图像等级,截取大量属于其宽度范围的道路块样本,所截取的道路块样本大小相同;
S12、将所截取的道路块样本转换为道路宽度描述符;
S121、以截取的道路块样本的中心为圆心,构建a个同心圆,a个同心圆的半径分别为道路块内切圆半径的1/a、2/a……7/a及1倍,将a个同心圆标记为Cn,其中n=1,2……a;
S122、根据道路块中的像素值分布,将道路块中的像素值分为b个等长的像素值区间,记为Im,其中m=1,2……b,令max和min分别为道路块中最大像素值和最小像素值,则每个像素值区间的长度为(max-min)/b;统计每个圆范围内落在各像素值区间里的像素个数,将第n个圆内属于第m个区间范围内的像素个数记为Nn,m,并将每个圆范围内落在各像素值区间内的像素个数归一化到0-255范围内;
S123、将第n个圆内的像素值分布作为描述符的第n行,获得道路宽度描述符;
S13、将道路宽度描述符制作为训练样本,所述训练样本包括训练集和验证集,利用卷积神经网络对训练样本进行学习,获得道路宽度类别预测模型。
3.根据权利要求2所述的基于高分辨率卫星影像的道路宽度估计方法,其特征在于:步骤S2采用与步骤S12相同的方式将其截取的道路块转化为道路宽度描述符。
4.根据权利要求3所述的基于高分辨率卫星影像的道路宽度估计方法,其特征在于:步骤S1与S2中,还包括采用图像放大算法将所获得的道路宽度描述符放大若干倍的步骤。
5.根据权利要求4所述的基于高分辨率卫星影像的道路宽度估计方法,其特征在于:所述图像放大算法为最近邻插值算法。
6.根据权利要求1所述的基于高分辨率卫星影像的道路宽度估计方法,其特征在于,步骤S3中S(L)的计算具体为:
其中,Np表示道路中心线上与像素p相邻的若干个像素的集合;li和lq分别表示以像素p和q为中心的道路块的宽度类别;‖p-q‖为像素p与像素q之间的距离;当li=lq时,符号函数sign(li,lq)取值为0,否则sign(li,lq)取值为1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710454604.7A CN107203761B (zh) | 2017-06-15 | 2017-06-15 | 基于高分辨率卫星影像的道路宽度估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710454604.7A CN107203761B (zh) | 2017-06-15 | 2017-06-15 | 基于高分辨率卫星影像的道路宽度估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107203761A CN107203761A (zh) | 2017-09-26 |
CN107203761B true CN107203761B (zh) | 2019-09-17 |
Family
ID=59908428
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710454604.7A Active CN107203761B (zh) | 2017-06-15 | 2017-06-15 | 基于高分辨率卫星影像的道路宽度估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107203761B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109558801B (zh) * | 2018-10-31 | 2020-08-07 | 厦门大学 | 道路网络提取方法、介质、计算机设备及系统 |
CN111121797B (zh) * | 2018-11-01 | 2021-11-12 | 百度在线网络技术(北京)有限公司 | 道路筛选方法、装置、服务器及存储介质 |
CN110348383B (zh) * | 2019-07-11 | 2020-07-31 | 重庆市地理信息和遥感应用中心(重庆市测绘产品质量检验测试中心) | 一种基于卷积神经网络回归的道路中心线和双线提取方法 |
CN113139454B (zh) * | 2021-04-19 | 2024-04-23 | 国交空间信息技术(北京)有限公司 | 基于单张图像的道路宽度提取方法及装置 |
CN113901342B (zh) * | 2021-09-30 | 2022-07-26 | 北京百度网讯科技有限公司 | 一种路网数据处理方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106446914A (zh) * | 2016-09-28 | 2017-02-22 | 天津工业大学 | 基于超像素和卷积神经网络的道路检测 |
CN106845547A (zh) * | 2017-01-23 | 2017-06-13 | 重庆邮电大学 | 一种基于摄像头的智能汽车定位与道路标识识别系统及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9514391B2 (en) * | 2015-04-20 | 2016-12-06 | Xerox Corporation | Fisher vectors meet neural networks: a hybrid visual classification architecture |
-
2017
- 2017-06-15 CN CN201710454604.7A patent/CN107203761B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106446914A (zh) * | 2016-09-28 | 2017-02-22 | 天津工业大学 | 基于超像素和卷积神经网络的道路检测 |
CN106845547A (zh) * | 2017-01-23 | 2017-06-13 | 重庆邮电大学 | 一种基于摄像头的智能汽车定位与道路标识识别系统及方法 |
Non-Patent Citations (4)
Title |
---|
Road network extraction via deep learning and line integral convolution;Peikang Li 等;《2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS)》;20160715;第1599-1602页 |
Simultaneous extraction of roads and buildings in remote sensing imagery with convolutional neural networks;Rasha Alshehhi 等;《ISPRS Journal of Photogrammetry and Remote Sensing》;20170609;第130卷;第139-149页 |
一种新的高分辨率遥感图像道路提取方法;王振国 等;《测绘与空间地理信息》;20101231;第33卷(第6期);第149-155页 |
基于深度卷积神经网络的场景自适应道路分割算法;王海 等;《电子与信息学报》;20170228;第39卷(第2期);第263-269页 |
Also Published As
Publication number | Publication date |
---|---|
CN107203761A (zh) | 2017-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107203761B (zh) | 基于高分辨率卫星影像的道路宽度估计方法 | |
Brenning | Benchmarking classifiers to optimally integrate terrain analysis and multispectral remote sensing in automatic rock glacier detection | |
Hamedianfar et al. | Improving detailed rule-based feature extraction of urban areas from WorldView-2 image and lidar data | |
Stow et al. | Object‐based classification of residential land use within Accra, Ghana based on QuickBird satellite data | |
Huang et al. | A comparative study of spatial approaches for urban mapping using hyperspectral ROSIS images over Pavia City, northern Italy | |
Jia et al. | Object-oriented feature selection of high spatial resolution images using an improved Relief algorithm | |
CN108830870B (zh) | 基于多尺度结构学习的卫星影像高精度农田边界提取方法 | |
Pham et al. | Urban growth and change analysis using remote sensing and spatial metrics from 1975 to 2003 for Hanoi, Vietnam | |
Cao et al. | A new change-detection method in high-resolution remote sensing images based on a conditional random field model | |
Xu et al. | A supervoxel approach to the segmentation of individual trees from LiDAR point clouds | |
Wang et al. | Segmentation of multispectral high-resolution satellite imagery based on integrated feature distributions | |
Zhang et al. | A neighbourhood-constrained k-means approach to classify very high spatial resolution hyperspectral imagery | |
CN103578110A (zh) | 基于灰度共生矩阵的多波段高分辨率遥感影像分割方法 | |
Peng et al. | Object-based change detection from satellite imagery by segmentation optimization and multi-features fusion | |
CN112069985B (zh) | 基于深度学习的高分辨率大田图像稻穗检测与计数方法 | |
Xiao et al. | Segmentation of multispectral high-resolution satellite imagery using log Gabor filters | |
Su et al. | Image segmentation using mean shift for extracting croplands from high-resolution remote sensing imagery | |
Su et al. | Classification of MODIS images combining surface temperature and texture features using the Support Vector Machine method for estimation of the extent of sea ice in the frozen Bohai Bay, China | |
Chen et al. | Plant leaf segmentation for estimating phenotypic traits | |
CN105354824A (zh) | 基于区域提取的双参数恒虚警检测方法 | |
Zhang et al. | Feature extraction for high-resolution imagery based on human visual perception | |
Li et al. | A label propagation method using spatial-spectral consistency for hyperspectral image classification | |
Chattaraj et al. | Semi-automated object-based landform classification modelling in a part of the Deccan Plateau of central India | |
Marvaniya et al. | Small, sparse, but substantial: techniques for segmenting small agricultural fields using sparse ground data | |
Paul et al. | Estimation of shade tree density in tea garden using remote sensing images and deep convolutional neural network |
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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Cang Yu Inventor after: Wang Cheng Inventor after: Wang Xiaofang Inventor after: Xia Zhichao Inventor after: Li Peikang Inventor after: Li Jun Inventor before: Cang Yu Inventor before: Wang Cheng Inventor before: Xia Zhichao Inventor before: Li Peikang Inventor before: Li Jun |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |