CN106296747A - 基于结构决策图的鲁棒多模型拟合方法 - Google Patents
基于结构决策图的鲁棒多模型拟合方法 Download PDFInfo
- Publication number
- CN106296747A CN106296747A CN201610592964.9A CN201610592964A CN106296747A CN 106296747 A CN106296747 A CN 106296747A CN 201610592964 A CN201610592964 A CN 201610592964A CN 106296747 A CN106296747 A CN 106296747A
- Authority
- CN
- China
- Prior art keywords
- hypothesis
- decision diagram
- data
- model
- structure decision
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
基于结构决策图的鲁棒多模型拟合方法,涉及鲁棒多模型拟合。包括以下步骤:对输入样本数据随机采样产生大量的假设;基于核密度估计和内点尺度对产生的每个假设计算权重分数;根据权重分数的大小,对所有的假设进行排序;基于连续一致性集和皮尔逊积矩相关系数对排序后的每个假设分别计算最短抵达距离;根据权重分数和最短抵达距离构建结构决策图;在结构决策图上确定所有结构对应的结构原型并计算结构数量;根据结构原型进行内点和异常点的划分,输出每个结构对应的模型参数,完成基于结构决策图的鲁棒多模型拟合。通过利用假设的一致性信息来选择结构原型。不涉及过滤或聚类过程,解决可能删除有代表性的假设和忽略较小的结构的问题。
Description
技术领域
本发明涉及鲁棒多模型拟合,特别是涉及一种基于结构决策图的鲁棒多模型拟合方法。
背景技术
在三维平面重建和运动分割等一系列计算机视觉应用中,来自多模型分布的观察数据通常包含大量的离群点。如何从这些包含离群点的数据中提取多结构信息是这些应用所面临的的一个主要挑战。包含离群点的鲁棒多模型拟合的任务包括:1)准确地估计数据中存在的结构数量;2)精确地恢复这些结构对应的模型参数。
在计算机视觉领域,大量的鲁棒多模型拟合算法已经被提出用来处理包含离群点的数据。在这些提出的算法中,一些算法通过使用数据点的倾向集(preference set)来优化一个拟合准则以实现拟合和分割数据中多结构的任务。例如,J-Linkage算法(R.Toldo,A.Fusiello.Robust multiple structures estimation with J-Linkage.EuropeanConference on Computer Vision,2008)被提出用来同时估计结构的数量及其对应的模型参数。但是该算法要求用户指定一个内点尺度参数来区分内点和异常点,而这个内点尺度参数通常需要手动设置。核拟合(kernel fitting)算法(T.-J.Chin,H.Wang,D.Suter.Robust fitting of multiple structures:The statistical learningapproach.IEEE International Conference on Computer Vision,2009)基于来自同一结构的数据点在再生核希尔伯特空间收敛于同一位置的假设实现模型拟合。该算法能在有效地移除数据中存在的大量离群点的同时找到数据中存在的多个结构。另外一些算法通过分析模型假设的一致性集(consensus set)来找到能很好地描述数据中有代表性的模型假设。比如,AKSWH算法(H.Wang,T.-J.Chin,D.Suter.Simultaneously fitting andsegmenting multiple-structure data with outliers.IEEE Transactions on PatternAnalysis and Machine Intelligence,2012)同时估计数据中结构的数量以及每个结构的尺度和参数。RCG算法(H.Liu,S.Yan.Efficient structure detection via randomconsensus graph.IEEE Conference on Computer Vision and Pattern Recognition,2012)从超图的观点提出了一个快速而鲁棒的多模型拟合算法。
AKSWH和核拟合使用“过滤和聚类”的策略在鲁棒多模型拟合上显示出很好的性能。具体来说,这两类算法首先对所有的模型假设或数据点根据它们的重要性赋予一个权重分数,然后基于一个过滤步骤去除那些不重要的假设或数据点(即实现去除离群点的目的),最后用一个聚类算法来对剩下的模型假设或数据点聚类,来获得多个结构的分割。但是,一个可能的问题就是当数据中结构不平衡时,过滤步骤可能会去除一些包含少量内点的较小结构。
发明内容
本发明的目的在于提供一种基于结构决策图的鲁棒多模型拟合方法。
本发明包括以下步骤:
A.对输入样本数据随机采样产生大量的假设;
B.基于核密度估计和内点尺度对产生的每个假设计算权重分数;
C.根据权重分数的大小,对所有的假设进行排序;
D.基于连续一致性集和皮尔逊积矩相关系数对排序后的每个假设分别计算最短抵达距离;
E.根据权重分数和最短抵达距离构建结构决策图;
F.在结构决策图上确定所有结构对应的结构原型并计算结构数量;
G.根据结构原型进行内点和异常点的划分,输出每个结构对应的模型参数,完成基于结构决策图的鲁棒多模型拟合。
在步骤A中,所述对输入样本数据随机采样产生大量的假设的具体方法可为:
A1.给定输入样本数据D={d1,d2,....,dN},其中di表示第i个样本数据;N为样本数目且N为自然数;
A2.对输入数据随机抽取p个数据点形成一个模型假设(例如,拟合一条直线,需要抽取两个数据点,即p=2;拟合一个圆,需要抽取三个数据点,即p=3;拟合一个平面,需要抽取三个数据点,即p=3;拟合一个运动目标,需要抽取八个数据点,即p=8);
A3.计算所有数据点到这个模型假设的距离,构成残差向量其中为对应于第j假设的第i个残差值;
在步骤B中,所述基于核密度估计和内点尺度对产生的每个假设计算权重分数的具体方法可为:
B1.计算每个假设Hj(j=1,...,M)的内点尺度估计νj,计算公式为:
其中为对应于第j个结构的第K个最小的绝对值残差;θ-1(·)为累积密度函数;Nj为对应于第j个假设的内点数目;K为常数,通常设置为所有数据点数目的10%,即假定数据中最小结构的内点数量为所有数据的10%;
B2.计算每个假设Hj(j=1,...,M)的核密度函数f(γ;Hj),计算公式为:
其中K为Epanechnikov核函数,即且贝塔函数Γ(·)是高斯函数且满足Γ(n)=(n-1)!;hj为基于核函数和内点噪声尺度估计νj得到的可变带宽,其计算公式为:
B3.计算核密度函数f(γ;Hj)在原点处的值f(0;Hj);
B4.计算每个假设Hj(j=1,...,M)的权重分数,计算公式为:
在步骤D中,所述基于连续一致性集和皮尔逊积矩相关系数对排序后的每个假设计算最短抵达距离的具体方法可为:
D1.根据步骤A中得到的每个假设Hj(j=1,...,M)的残差向量得到每个假设对应的连续一致性集,计算公式为:
其中,当时,否则E0为一常数阈值,通常设置为2.5(即假设数据服从高斯分布时,当E0=2.5时,大约98%的数据点被认为是内点);
D2.利用皮尔逊积矩相关系数计算任意两个假设(假定为Hj和Hk)之间的相似性,计算公式为:
其中,协方差cov(Δ(Hj),Δ(Hk))=E(Δ(Hj)-μ)E(Δ(Hk)-μ),μ为所有假设对应的连续一致性集的均值。E(.)是期望操作符;
D3.对于权重分数最高的假设其最小抵达距离MinAD的计算公式如下:
D4.对于其他的假设其最小抵达距离MinAD定义为当前的随机假设到任何其他比当前的假设具有更高的权重分数的随机假设的最短的距离,计算公式如下:
在步骤E中,所述构建结构决策图的方式可为:将所有假设的最短抵达距离MinAD按照从小到大的递增顺序排序,作为纵轴(Y轴),而对应的权重分数作为横轴(X轴),从而构成一个结构决策图。
在步骤F中,所述在结构决策图上确定所有结构对应的结构原型并计算结构数量的具体方法包括以下步骤:
F1.计算相邻两个假设点的最短抵达距离MinAD之间的差值,得到一个差值序列;
F2.根据差值序列中差值最大的元素位置k,从该元素位置k到差值序列最后一个元素位置M之间所对应的所有假设均认为是结构原型;
F3.计算结构数量为M-k+1。
本发明提出一个新颖的基于结构决策图的鲁棒多模型拟合算法来拟合和分割存在大量异常点的数据中的多个结构。本发明基于下面的观察:1)对于多结构数据,每个结构往往可以用至少一个有代表性的假设来表示,在这里把这个有代表性的假设称之为结构原型;2)这些结构原型间通常有一个相对较大的距离。这些观察启发在多模型拟合中引入结构决策图的概念。通过结构决策图,多结构数据的分割和拟合问题可以很容易地解决。此外,本发明提出的方法不需要任何迭代过程。
与现有技术相比,本发明提出基于结构决策图的方法通过利用假设的一致性信息来选择结构原型。与现有的鲁棒多模型拟合方法相比,本发明不涉及过滤或聚类过程,从而解决了这些过程可能删除有代表性的假设(或内点)和忽略较小的结构的问题。
附图说明
图1为本发明实施例的流程图。
图2为本发明与其它几种鲁棒多模型拟合方法在四直线数据集(内点高斯噪声尺度为σ=0.01和异常点比例为90%)上的对比拟合结果图。
图3为本发明实施例的两视角运动分割实例图。其中标记:乘号为离群点;方块、菱形和圆圈分别对应数据中的不同的结构。
具体实施方式
下面结合附图和实施例对本发明的方法作详细说明。
参见图1,本发明实施例包括以下步骤:
S1.对输入样本数据随机采样产生大量的假设。
具体包括:
(1)给定输入样本数据D={d1,d2,....,dN},其中di表示第i个样本数据;N为样本数目且N为自然数。
(2)对输入数据随机抽取p个数据点形成一个模型假设(例如,拟合一条直线,需要抽取两个数据点,即p=2;拟合一个圆,需要抽取三个数据点,即p=3;拟合一个平面,需要抽取三个数据点,即p=3;拟合一个运动目标,需要抽取八个数据点,即p=8)。
(3)计算所有数据点到这个模型假设的距离,构成残差向量其中为对应于第j假设的第i个残差值。
S2.基于核密度估计和内点尺度对产生的每个假设计算权重分数。
具体包括:
(1)计算每个假设Hj(j=1,...,M)的内点尺度估计νj,计算公式为:
其中为对应于第j个结构的第K个最小的绝对值残差;θ-1(·)为累积密度函数;Nj为对应于第j个假设的内点数目;K为常数,通常设置为所有数据点数目的10%,即假定数据中最小结构的内点数量为所有数据的10%。
(2)计算每个假设Hj(j=1,...,M)的核密度函数f(γ;Hj),计算公式为:
其中K为Epanechnikov核函数,即且贝塔函数Γ(·)是高斯函数且满足Γ(n)=(n-1)!;hj为基于核函数和内点噪声尺度估计νj得到的可变带宽,其计算公式为:
(3)计算核密度函数f(γ;Hj)在原点处的值f(0;Hj)。
(4)计算每个假设Hj(j=1,...,M)的权重分数,计算公式为:
S3.根据权重分数的大小,对所有的假设进行排序得到M个排序后的假设其中λ={λ1,λ2,...,λM}满足
S4.基于连续一致性集和皮尔逊积矩相关系数对排序后的每个假设分别计算最短抵达距离。
具体包括:
(1)根据步骤S1中得到的每个假设Hj(j=1,...,M)的残差向量得到每个假设对应的连续一致性集,计算公式为:
其中,当时,否则E0为一常数阈值,通常设置为2.5(即假设数据服从高斯分布时,大约98%的数据点被认为是内点)。
(2)利用皮尔逊积矩相关系数计算任意两个假设(假定为Hj和Hk)之间的相似性,计算公式为:
其中,协方差cov(Δ(Hj),Δ(Hk))=E(Δ(Hj)-μ)E(Δ(Hk)-μ)。μ为所有假设对应的连续一致性集的均值。E(.)是期望操作符。
(3)对于权重分数最高的假设其最小抵达距离MinAD的计算公式如下:
(4)对于其他的假设其最小抵达距离MinAD定义为当前的随机假设到任何其他比当前的假设具有更高的权重分数的随机假设的最短的距离,计算公式如下:
S5.将所有假设的最短抵达距离MinAD按照从小到大的递增顺序排序,作为纵轴(Y轴),而对应的权重分数作为横轴(X轴),从而构成一个结构决策图。
S6.在结构决策图上确定所有结构对应的结构原型并计算结构数量。
具体包括:
(1)计算相邻两个假设点的最短抵达距离MinAD之间的差值,得到一个差值序列。
(2)根据差值序列中差值最大的元素位置k,从该元素位置k到差值序列最后一个元素位置M之间所对应的所有假设均认为是结构原型。
(3)计算结构数量为M-k+1。
S7.根据结构原型进行内点和异常点的划分。对于第p(p=1,2,....,M-k+1)个结构原型,如果其中为对应于第i个样本数据到第p个结构原型的残差;νp为每个假设的内点尺度估计;E0为一常数阈值,通常设置为2.5,则第i个样本数据判断为该结构原型的内点,否则判断为离群点。根据内点输出每个结构对应的模型参数,完成最终模型拟合。
本发明与其它几种鲁棒多模型拟合方法在四直线数据集(内点高斯噪声尺度为σ=0.01和异常点比例为90%)上的对比拟合结果图参见图2。在图2中,(a)为初始输入数据;(b)方法1对应为R.Toldo等人提出的方法(R.Toldo,A.Fusiello.Robust multiplestructures estimation with J-Linkage.European Conference on Computer Vision,2008);(c)方法2对应为T.-J.Chin等人提出的方法(T.-J.Chin,H.Wang,D.Suter.Robustfitting of multiple structures:The statistical learning approach.IEEEInternational Conference on Computer Vision,2009);(d)方法3对应为H.Liu等人提出的方法(H.Liu,S.Yan.Efficient structure detection via random consensus graph[C].IEEE Conference on Computer Vision and Pattern Recognition,2012);(e)方法4对应为H.Wang等人提出的方法(H.Wang,T.-J.Chin,D.Suter.Simultaneously fittingand segmenting multiple-structure data with outliers.IEEE Transactions onPattern Analysis and Machine Intelligence,2012);(f)方法5对应为L.Magri等人提出的方法(L.Magri,A.Fusiello.T-Linkage:A continuous relaxation of J-Linkage formulti-model fitting[C].IEEE Conference on Computer Vision and PatternRecognition,2014);(g)方法6为本发明提出的方法。
本发明实施例的两视角运动分割实例图参见图3。
本发明可以对多结构数据进行有效的鲁棒模型拟合。相对于传统的鲁棒多模型拟合算法,本发明不仅可以自动找到结构的数量及对应的模型参数,而且还可以有效地处理多结构数据中的不平衡问题。
Claims (6)
1.基于结构决策图的鲁棒多模型拟合方法,其特征在于包括以下步骤:
A.对输入样本数据随机采样产生大量的假设;
B.基于核密度估计和内点尺度对产生的每个假设计算权重分数;
C.根据权重分数的大小,对所有的假设进行排序;
D.基于连续一致性集和皮尔逊积矩相关系数对排序后的每个假设分别计算最短抵达距离;
E.根据权重分数和最短抵达距离构建结构决策图;
F.在结构决策图上确定所有结构对应的结构原型并计算结构数量;
G.根据结构原型进行内点和异常点的划分,输出每个结构对应的模型参数,完成基于结构决策图的鲁棒多模型拟合。
2.如权利要求1所述基于结构决策图的鲁棒多模型拟合方法,其特征在于在步骤A中,所述对输入样本数据随机采样产生大量的假设的具体方法为:
A1.给定输入样本数据D={d1,d2,....,dN},其中di表示第i个样本数据;N为样本数目且N为自然数;
A2.对输入数据随机抽取p个数据点形成一个模型假设,所述假设,当拟合一条直线时,抽取两个数据点,即p=2;当拟合一个圆时,抽取三个数据点,即p=3;当拟合一个平面时,抽取三个数据点,即p=3;当拟合一个运动目标时,抽取八个数据点,即p=8;
A3.计算所有数据点到这个模型假设的距离,构成残差向量其中为对应于第j假设的第i个残差值。
3.如权利要求1所述基于结构决策图的鲁棒多模型拟合方法,其特征在于在步骤B中,所述基于核密度估计和内点尺度对产生的每个假设计算权重分数的具体方法为:
B1.计算每个假设Hj(j=1,...,M)的内点尺度估计vj,计算公式为:
其中为对应于第j个结构的第K个最小的绝对值残差;θ-1(·)为累积密度函数;Nj为对应于第j个假设的内点数目;K为常数,通常设置为所有数据点数目的10%,即假定数据中最小结构的内点数量为所有数据的10%;
B2.计算每个假设Hj(j=1,...,M)的核密度函数f(γ;Hj),计算公式为:
其中K为Epanechnikov核函数,即且贝塔函数Γ(·)是高斯函数且满足Γ(n)=(n-1)!;hj为基于核函数和内点噪声尺度估计vj得到的可变带宽,其计算公式为:
B3.计算核密度函数f(γ;Hj)在原点处的值f(0;Hj);
B4.计算每个假设Hj(j=1,...,M)的权重分数,计算公式为:
4.如权利要求1所述基于结构决策图的鲁棒多模型拟合方法,其特征在于在步骤D中,所述基于连续一致性集和皮尔逊积矩相关系数对排序后的每个假设计算最短抵达距离的具体方法为:
D1.根据步骤A中得到的每个假设Hj(j=1,...,M)的残差向量得到每个假设对应的连续一致性集,计算公式为:
其中,当时,否则E0为一常数阈值,通常设置为2.5,即假设数据服从高斯分布,当E0=2.5时,98%的数据点被认为是内点;
D2.利用皮尔逊积矩相关系数计算任意两个假设之间的相似性,所述两个假设记为Hj和Hk,计算公式为:
其中,协方差cov(Δ(Hj),Δ(Hk))=E(Δ(Hj)-μ)E(Δ(Hk)-μ),μ为所有假设对应的连续一致性集的均值,E(.)是期望操作符;
D3.对于权重分数最高的假设其最小抵达距离MinAD的计算公式如下:
D4.对于其他的假设其最小抵达距离MinAD定义为当前的随机假设到任何其他比当前的假设具有更高的权重分数的随机假设的最短的距离,计算公式如下:
5.如权利要求1所述基于结构决策图的鲁棒多模型拟合方法,其特征在于在步骤E中,所述构建结构决策图的方式为:将所有假设的最短抵达距离MinAD按照从小到大的递增顺序排序,作为纵轴,而对应的权重分数作为横轴,从而构成一个结构决策图。
6.如权利要求1所述基于结构决策图的鲁棒多模型拟合方法,其特征在于在步骤F中,所述在结构决策图上确定所有结构对应的结构原型并计算结构数量的具体方法包括以下步骤:
F1.计算相邻两个假设点的最短抵达距离MinAD之间的差值,得到一个差值序列;
F2.根据差值序列中差值最大的元素位置k,从该元素位置k到差值序列最后一个元素位置M之间所对应的所有假设均认为是结构原型;
F3.计算结构数量为M-k+1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610592964.9A CN106296747A (zh) | 2016-07-26 | 2016-07-26 | 基于结构决策图的鲁棒多模型拟合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610592964.9A CN106296747A (zh) | 2016-07-26 | 2016-07-26 | 基于结构决策图的鲁棒多模型拟合方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106296747A true CN106296747A (zh) | 2017-01-04 |
Family
ID=57652534
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610592964.9A Pending CN106296747A (zh) | 2016-07-26 | 2016-07-26 | 基于结构决策图的鲁棒多模型拟合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106296747A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107123119A (zh) * | 2017-04-27 | 2017-09-01 | 厦门大学 | 一种针对多结构数据的指导性采样方法 |
CN108921871A (zh) * | 2018-07-27 | 2018-11-30 | 福建农林大学 | 一种基于全局贪婪搜索的鲁棒模型拟合方法 |
CN110163865A (zh) * | 2019-05-28 | 2019-08-23 | 闽江学院 | 一种针对模型拟合中不平衡数据的采样方法 |
-
2016
- 2016-07-26 CN CN201610592964.9A patent/CN106296747A/zh active Pending
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107123119A (zh) * | 2017-04-27 | 2017-09-01 | 厦门大学 | 一种针对多结构数据的指导性采样方法 |
CN107123119B (zh) * | 2017-04-27 | 2019-12-17 | 厦门大学 | 一种针对多结构数据的指导性采样方法 |
CN108921871A (zh) * | 2018-07-27 | 2018-11-30 | 福建农林大学 | 一种基于全局贪婪搜索的鲁棒模型拟合方法 |
CN108921871B (zh) * | 2018-07-27 | 2022-02-01 | 福建农林大学 | 一种基于全局贪婪搜索的鲁棒模型拟合方法 |
CN110163865A (zh) * | 2019-05-28 | 2019-08-23 | 闽江学院 | 一种针对模型拟合中不平衡数据的采样方法 |
CN110163865B (zh) * | 2019-05-28 | 2021-06-01 | 闽江学院 | 一种针对模型拟合中不平衡数据的采样方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Harouni et al. | Online Persian/Arabic script classification without contextual information | |
CN110287942B (zh) | 年龄估计模型的训练方法、年龄估计方法以及对应的装置 | |
CN110852168A (zh) | 基于神经架构搜索的行人重识别模型构建方法及装置 | |
CN107633226B (zh) | 一种人体动作跟踪特征处理方法 | |
CN108229588B (zh) | 一种基于深度学习的机器学习识别方法 | |
CN110674752A (zh) | 一种基于隐马尔可夫模型的刀具磨损状态识别与预测方法 | |
CN104392253B (zh) | 一种草图数据集的交互式类别标注方法 | |
JPWO2014002554A1 (ja) | 画像処理装置、画像処理方法、及びプログラム | |
CN111027140B (zh) | 基于多视角点云数据的飞机标准件模型快速重构方法 | |
CN101477529A (zh) | 一种三维对象的检索方法和装置 | |
CN108171119B (zh) | 基于残差网络的sar图像变化检测方法 | |
CN106296747A (zh) | 基于结构决策图的鲁棒多模型拟合方法 | |
CN108764280A (zh) | 一种基于症状向量的医学数据处理方法和系统 | |
CN111862167B (zh) | 一种基于稀疏紧凑相关滤波器的快速鲁棒目标跟踪方法 | |
Gerych et al. | Recovering the propensity score from biased positive unlabeled data | |
CN108710761A (zh) | 一种基于谱聚类去除离群点的鲁棒模型拟合方法 | |
CN113516019B (zh) | 高光谱图像解混方法、装置及电子设备 | |
CN108921853B (zh) | 基于超像素和免疫稀疏谱聚类的图像分割方法 | |
Ananth et al. | A geodesic active contour level set method for image segmentation | |
CN109190505A (zh) | 基于视觉理解的图像识别方法 | |
Firouznia et al. | Adaptive chaotic sampling particle filter to handle occlusion and fast motion in visual object tracking | |
CN104200222B (zh) | 一种基于因子图模型的图片中对象识别方法 | |
CN110223378A (zh) | 基于层次聚类的二次曲面基元提取方法、系统、装置 | |
CN109165586A (zh) | 用于ai芯片的智能图像处理方法 | |
CN109165587A (zh) | 智能图像信息抽取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170104 |