CN109101949A - 一种基于彩色视频信号频域分析的人脸活体检测方法 - Google Patents

一种基于彩色视频信号频域分析的人脸活体检测方法 Download PDF

Info

Publication number
CN109101949A
CN109101949A CN201810996588.9A CN201810996588A CN109101949A CN 109101949 A CN109101949 A CN 109101949A CN 201810996588 A CN201810996588 A CN 201810996588A CN 109101949 A CN109101949 A CN 109101949A
Authority
CN
China
Prior art keywords
face
video
signal
component analysis
independent component
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
Application number
CN201810996588.9A
Other languages
English (en)
Inventor
夏祥红
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangzhou Honghuang Intelligent Technology Co Ltd
Original Assignee
Guangzhou Honghuang Intelligent Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangzhou Honghuang Intelligent Technology Co Ltd filed Critical Guangzhou Honghuang Intelligent Technology Co Ltd
Priority to CN201810996588.9A priority Critical patent/CN109101949A/zh
Publication of CN109101949A publication Critical patent/CN109101949A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/16Human faces, e.g. facial parts, sketches or expressions
    • G06V40/161Detection; Localisation; Normalisation
    • G06V40/165Detection; Localisation; Normalisation using facial parts and geometric relationships
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2135Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/40Spoof detection, e.g. liveness detection
    • G06V40/45Detection of the body part being alive

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Data Mining & Analysis (AREA)
  • Multimedia (AREA)
  • Human Computer Interaction (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Evolutionary Biology (AREA)
  • Geometry (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Image Analysis (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

一种基于彩色视频信号频域分析的人脸活体检测方法,采用以下步骤,步骤1:图像流采集装置采集人脸视频,将该人脸视频发送到处理模块;步骤2:处理模块对人脸视频中的人脸进行检测识别;步骤3:处理模块选取人脸图像的上额区域作为目标区域;步骤4:获取目标区域内的RGB颜色均值,得到三个通道的均值序列x1(t)、x2(t)、x3(t),t表示时刻;步骤5:对RGB信号进行独立成分分析,得到三个独立成分分析的分量为s1(t),s2(t),s3(t)。使用方法为用户非接触式、无需完成任何动作即可完成人脸活体检测,用户体验好。

Description

一种基于彩色视频信号频域分析的人脸活体检测方法
技术领域
本发明涉及人脸识别领域,具体涉及一种基于彩色视频信号频域分析的人脸活体检测方法。
背景技术
人脸识别技术是目前应用最广泛的生物特征识别技术之一,它具有方便易用、用户友好性、非接触式等优点。随着各行各业的广泛应用,使用安全性便成为普遍关注的问题。在人脸识别系统中,非法分子很容易采用照片、面具假冒伪造人脸来实施攻击和欺诈行为。因此,人脸活体检测技术,逐渐受到学术界和工业界的重视,其目的就是区分真实和虚假人脸。
目前的人脸活体检测技术主要有以下几种:
1、基于特征学习的方法:由于同一设备采集的真实人脸和伪造人脸相比存在纹理细节的细微差异,表面形状的差异、以及局部高光差异等不同。这类方法对真实人脸(正样本)和伪造人脸(负样本)分别进行特征提取,通过设计分类器进行模型训练,从而实现真假人脸的辨别。这类算法的缺点是,同类样本的类内差别可能会很大,对分类器分类性能有影响,即选取的真实人脸样本可能是不同的人,不同的角度、表情、环境,很难学习到区分人脸真实与否的特征。
2、基于人机交互的方法:被检测者需按系统要求做出相应动作,如点头、眨眼等,通过分析人脸的动作模式来区分真假人脸。这类方法的缺点是,对用户要求过高,使得用户体验不佳,且认证时间较长。
3、采用多模式融合判别方法:如采用红外图像分析人脸、深度摄像头获取三维信息、联合语音识别等方法,这类方法需要额外的设备辅助,成本较高。为解决以上活体检测方法的缺点,本文发明提出了一种彩色视频信号频域分析的方法验证人脸的真实性。该方法的主要依据是,真实人的心脏搏动会引起人脸血管血液流量的周期性变化,而这种细微变化直接反应在人脸RGB信号上,通过对人脸彩色图像序列的信号进行傅里叶频域分析,发现在心率频率附近具有较高的能量,以此来判断人脸真实与否。本发明方法的优点是,直接利用获取的彩色视频,无需额外设备辅助,无需被检测者的主动配合,计算简单快速,能够很好的解决基于人脸识别的身份认证中伪造人脸的欺骗问题。
发明内容
本发明针对现有技术的不足,提出一种基于彩色视频信号频域分析的人脸活体检测方法,具体技术方案如下:一种基于彩色视频信号频域分析的人脸活体检测方法,其特征在于:
采用以下步骤,
步骤1:图像流采集装置采集人脸视频,将该人脸视频发送到处理模块;
步骤2:处理模块对人脸视频中的人脸进行检测识别;
步骤3:处理模块选取人脸图像的上额区域作为目标区域;
步骤4:获取目标区域内的RGB颜色均值,得到三个通道的均值序列x1(t)、x2(t)、x3(t),t表示时刻;
步骤5:对RGB信号进行独立成分分析,得到三个独立成分分析的分量为
s1(t),s2(t),s3(t);
步骤6:提取独立成分分析的第二个分量s2(t),对s2(t)进行傅里叶频域分析,获取信号的能量谱E(ω),即
E(ω)=|F(s2(t))|
其中F表示傅里叶变换,ω为频率;
步骤7:通过统计1.07Hz频率附近的幅值占整个能量谱的比例来判断是否为真实人脸,即计算
步骤8:设定阈值h,判断Ks是否大于h,如果是,则进入下一步骤,否则,进入步骤;
步骤9:判断为真实人脸;
步骤10:判断为伪造人脸。
进一步地:所述步骤2具体为:
处理模块对人脸视频中逐帧进行人脸检测,如果有其中某一帧没有检测到人脸,则可以采取上一帧的结果,如果连续的至少三帧没有检测结果,则判断为人脸从视频中消失,当某一帧检测出多张人脸时,采用与上一帧位置最接近的矩形区域为当前检测结果。
进一步地:所述步骤5具体为,采用独立成分分析对RGB颜色信号进行独立分解,令列向量x(t)=[x1(t),x2(t),x3(t)]T为RGB观测信号,s(t)=[s1(t),s2(t),s3(t)]T为假设的三个独立的源信号,A为3x3的包含混合参数的矩阵,则观测信号可由独立源信号的线性组合而成,即x(t)=As(t);
独立成分分析的目的是寻找一个分离矩阵W,作为混合矩阵A的逆矩阵估计,从而得到源信号的估计满足
采用特征矩阵对角化联合估计算法求解独立成分分析问题,得到三个独立成分分析分量s1(t),s2(t),s3(t)。
本发明的有益效果为:第一,无需对硬件设备、摄像头进行改造,现有大部分摄像头即可满足要求,降低用户成本;
第二,使用方法为用户非接触式、无需完成任何动作即可完成人脸活体检测,用户体验好。
第三,利用本发明,能够解决目前利用虚假人脸,具体如人脸照片、人脸视频或者面具等进行违法行为的问题,且算法简单、快速、高效。
附图说明
图1为本发明的工作流程图。
具体实施方式
下面结合附图对本发明的较佳实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
如图1所示:一种基于彩色视频信号频域分析的人脸活体检测方法,采用以下步骤,
步骤1:图像流采集装置采集人脸视频,将该人脸视频发送到处理模块,图像流采集装置包括一个分辨率大于30万像素、帧率为30帧/秒的摄像头,一台连接摄像头的PC机,在稳定的可见光环境下,采用摄像头实时获取或录制一段人脸正面的彩色视频,录制期间受检者与摄像头保持半米的距离,上额部没有遮挡且头部位置和面部表情基本保持不变,可有轻微表情变化,如正常眨眼动作;
步骤2:处理模块对人脸视频中的人脸进行检测识别,对采集的彩色视频逐帧进行人脸检测,检测算法采用Viola-Jones算法,该算法采用Haar-like特征和Adaboost级联分类器实现人脸的检测。为防止整个过程的漏检和误检,采取以下措施:
处理模块对人脸视频中逐帧进行人脸检测,如果有其中某一帧没有检测到人脸,则可以采取上一帧的结果,如果连续的至少三帧没有检测结果,则判断为人脸从视频中消失,当某一帧检测出多张人脸时,采用与上一帧位置最接近的矩形区域为当前检测结果,该处可采用Opencv实现,获得检测结果为人脸矩形框左上角坐标以及矩形宽高;
步骤3:考虑到诸多影响脸部颜色变化的因素,如眨眼、鼻子对脸部光线遮挡、人脸边缘背景等,处理模块选取人脸图像的上额区域作为目标区域-ROI;
步骤4:在每一帧获取目标区域内的RGB三个通道的像素值,且分别计算三个通道在区域内的颜色均值,得到三个通道的均值序列x1(t)、x2(t)、x3(t),t表示时刻;
步骤5:RGB颜色传感器实际上获取的是不同反射光源的混合信号,因此RGB三通道之间具有很大程度的相关性。独立成分分析-ICA,常用于盲源信号分离,即从混合信号中分离出相互独立的信号分量;
对RGB信号进行独立成分分析,该独立成分分析可用ICA作为缩写替代词,采用独立成分分析对RGB颜色信号进行独立分解,令列向量x(t)=[x1(t),x2(t),x3(t)]T为RGB观测信号,s(t)=[s1(t),s2(t),s3(t)]T为假设的三个独立的源信号,A为3x3的包含混合参数的矩阵,则观测信号可由独立源信号的线性组合而成,即x(t)=As(t);
独立成分分析的目的是寻找一个分离矩阵W,作为混合矩阵A的逆矩阵估计,从而得到源信号的估计满足
采用特征矩阵对角化联合估计算法求解独立成分分析问题,得到三个独立成分分析分量s1(t),s2(t),s3(t)。
步骤6:由于人的心脏搏动使人脸额头部的血管血液容积产生周期性变化,可反映在反射光信号中,通过大量实验验证,独立成分分析分量中第二个分量s2(t)更多的包含真实人的心率信息,因此提取独立成分分析的第二个分量s2(t),对s2(t)进行傅里叶频域分析,获取信号的能量谱E(ω),即
E(ω)=|F(s2(t))|
其中F表示傅里叶变换,ω为频率。
步骤7:通过实验发现,真实人脸独立成分分析的第二分量信号的幅频响应曲线中在频率为1.07Hz附近有明显的峰值,即信号在此频率附近能量较大,这是由于信号中隐含了心率信号,符合正常人的的心率65bpm;而对于虚假人脸图像,其独立成分分析的第二分量的能量谱则不具有此特性。通过统计1.07Hz频率附近的幅值占整个能量谱的比例来判断是否为真实人脸,即计算
步骤8:设定阈值h,判断Ks是否大于h,如果是,则进入下一步骤,否则,进入步骤;
步骤9:判断为真实人脸;
步骤10:判断为伪造人脸。

Claims (3)

1.一种基于彩色视频信号频域分析的人脸活体检测方法,其特征在于:
采用以下步骤,
步骤1:图像流采集装置采集人脸视频,将该人脸视频发送到处理模块;
步骤2:处理模块对人脸视频中的人脸进行检测识别;
步骤3:处理模块选取人脸图像的上额区域作为目标区域;
步骤4:获取目标区域内的RGB颜色均值,得到三个通道的均值序列x1(t)、x2(t)、x3(t),t表示时刻;
步骤5:对RGB信号进行独立成分分析,得到三个独立成分分析的分量为s1(t),s2(t),s3(t);
步骤6:提取独立成分分析的第二个分量s2(t),对s2(t)进行傅里叶频域分析,获取信号的能量谱E(ω),即
其中F表示傅里叶变换,ω为频率;
步骤7:通过统计1.07Hz频率附近的幅值占整个能量谱的比例来判断是否为真实人脸,即计算
步骤8:设定阈值h,判断Ks是否大于h,如果是,则进入下一步骤,否则,进入步骤;
步骤9:判断为真实人脸;
步骤10:判断为伪造人脸。
2.一种基于彩色视频信号频域分析的人脸活体检测方法,其特征在于:所述步骤2具体为:
处理模块对人脸视频中逐帧进行人脸检测,如果有其中某一帧没有检测到人脸,则可以采取上一帧的结果,如果连续的至少三帧没有检测结果,则判断为人脸从视频中消失,当某一帧检测出多张人脸时,采用与上一帧位置最接近的矩形区域为当前检测结果。
3.一种基于彩色视频信号频域分析的人脸活体检测方法,其特征在于:所述步骤5具体为,采用独立成分分析对RGB颜色信号进行独立分解,令列向量x(t)=[x1(t),x2(t),x3(t)]T为RGB观测信号,s(t)=[s1(t),s2(t),s3(t)]T为假设的三个独立的源信号,A为3x3的包含混合参数的矩阵,则观测信号可由独立源信号的线性组合而成,即x(t)=As(t);
独立成分分析的目的是寻找一个分离矩阵W,作为混合矩阵A的逆矩阵估计,从而得到源信号的估计满足
采用特征矩阵对角化联合估计算法求解独立成分分析问题,得到三个独立成分分析分量s1(t),s2(t),s3(t)。
CN201810996588.9A 2018-08-29 2018-08-29 一种基于彩色视频信号频域分析的人脸活体检测方法 Pending CN109101949A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810996588.9A CN109101949A (zh) 2018-08-29 2018-08-29 一种基于彩色视频信号频域分析的人脸活体检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810996588.9A CN109101949A (zh) 2018-08-29 2018-08-29 一种基于彩色视频信号频域分析的人脸活体检测方法

Publications (1)

Publication Number Publication Date
CN109101949A true CN109101949A (zh) 2018-12-28

Family

ID=64864157

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810996588.9A Pending CN109101949A (zh) 2018-08-29 2018-08-29 一种基于彩色视频信号频域分析的人脸活体检测方法

Country Status (1)

Country Link
CN (1) CN109101949A (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110163174A (zh) * 2019-05-27 2019-08-23 成都科睿埃科技有限公司 一种基于单目摄像头的活体人脸检测方法
CN110269600A (zh) * 2019-08-06 2019-09-24 合肥工业大学 基于多元经验模态分解与联合盲源分离的非接触式视频心率检测方法
CN110279406A (zh) * 2019-05-06 2019-09-27 苏宁金融服务(上海)有限公司 一种基于摄像头的无接触式的脉率测量方法及装置
CN110298312A (zh) * 2019-06-28 2019-10-01 北京旷视科技有限公司 活体检测方法、装置、电子设备及计算机可读存储介质
CN110384491A (zh) * 2019-08-21 2019-10-29 河南科技大学 一种基于普通摄像头的心率检测方法
CN111523438A (zh) * 2020-04-20 2020-08-11 支付宝实验室(新加坡)有限公司 一种活体识别方法、终端设备和电子设备
CN112488137A (zh) * 2019-09-11 2021-03-12 广州虎牙科技有限公司 样本获取方法、装置、电子设备及机器可读存储介质
CN112507922A (zh) * 2020-12-16 2021-03-16 平安银行股份有限公司 人脸活体检测方法、装置、电子设备及存储介质
US11126829B2 (en) 2019-11-18 2021-09-21 Wistron Corporation Live facial recognition system and method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101196997A (zh) * 2007-12-29 2008-06-11 北京中星微电子有限公司 图像中最大人脸的跟踪装置和方法
CN104545864A (zh) * 2014-12-25 2015-04-29 中国科学院深圳先进技术研究院 心理调节方法和装置
CN104808776A (zh) * 2014-01-24 2015-07-29 北京奇虎科技有限公司 检测头戴式智能设备持续附着在人体上的装置和方法
CN106063702A (zh) * 2016-05-23 2016-11-02 南昌大学 一种基于人脸视频图像的心率检测系统及检测方法
CN106156757A (zh) * 2016-08-02 2016-11-23 中国银联股份有限公司 结合活体检测技术的人脸识别方法及人脸识别系统
CN106845395A (zh) * 2017-01-19 2017-06-13 北京飞搜科技有限公司 一种基于人脸识别进行活体检测的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101196997A (zh) * 2007-12-29 2008-06-11 北京中星微电子有限公司 图像中最大人脸的跟踪装置和方法
CN104808776A (zh) * 2014-01-24 2015-07-29 北京奇虎科技有限公司 检测头戴式智能设备持续附着在人体上的装置和方法
CN104545864A (zh) * 2014-12-25 2015-04-29 中国科学院深圳先进技术研究院 心理调节方法和装置
CN106063702A (zh) * 2016-05-23 2016-11-02 南昌大学 一种基于人脸视频图像的心率检测系统及检测方法
CN106156757A (zh) * 2016-08-02 2016-11-23 中国银联股份有限公司 结合活体检测技术的人脸识别方法及人脸识别系统
CN106845395A (zh) * 2017-01-19 2017-06-13 北京飞搜科技有限公司 一种基于人脸识别进行活体检测的方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110279406A (zh) * 2019-05-06 2019-09-27 苏宁金融服务(上海)有限公司 一种基于摄像头的无接触式的脉率测量方法及装置
CN110163174A (zh) * 2019-05-27 2019-08-23 成都科睿埃科技有限公司 一种基于单目摄像头的活体人脸检测方法
CN110298312A (zh) * 2019-06-28 2019-10-01 北京旷视科技有限公司 活体检测方法、装置、电子设备及计算机可读存储介质
CN110269600A (zh) * 2019-08-06 2019-09-24 合肥工业大学 基于多元经验模态分解与联合盲源分离的非接触式视频心率检测方法
CN110269600B (zh) * 2019-08-06 2021-12-21 合肥工业大学 基于多元经验模态分解与联合盲源分离的非接触式视频心率检测方法
CN110384491A (zh) * 2019-08-21 2019-10-29 河南科技大学 一种基于普通摄像头的心率检测方法
CN112488137A (zh) * 2019-09-11 2021-03-12 广州虎牙科技有限公司 样本获取方法、装置、电子设备及机器可读存储介质
US11126829B2 (en) 2019-11-18 2021-09-21 Wistron Corporation Live facial recognition system and method
TWI743593B (zh) * 2019-11-18 2021-10-21 緯創資通股份有限公司 活體臉部辨識系統與方法
CN111523438A (zh) * 2020-04-20 2020-08-11 支付宝实验室(新加坡)有限公司 一种活体识别方法、终端设备和电子设备
CN111523438B (zh) * 2020-04-20 2024-02-23 支付宝实验室(新加坡)有限公司 一种活体识别方法、终端设备和电子设备
CN112507922A (zh) * 2020-12-16 2021-03-16 平安银行股份有限公司 人脸活体检测方法、装置、电子设备及存储介质
CN112507922B (zh) * 2020-12-16 2023-11-07 平安银行股份有限公司 人脸活体检测方法、装置、电子设备及存储介质

Similar Documents

Publication Publication Date Title
CN109101949A (zh) 一种基于彩色视频信号频域分析的人脸活体检测方法
CN104915649B (zh) 一种应用于人脸识别的活体检测方法
CN108921100B (zh) 一种基于可见光图像与红外图像融合的人脸识别方法及系统
CN105574518B (zh) 人脸活体检测的方法和装置
Steiner et al. Reliable face anti-spoofing using multispectral SWIR imaging
Gu et al. Feature points extraction from faces
CN110837784B (zh) 一种基于人体头部特征的考场偷窥作弊检测系统
CN104361326A (zh) 一种判别活体人脸的方法
CN107403142B (zh) 一种微表情的检测方法
WO2019127262A1 (zh) 基于云端的人脸活体检测方法、电子设备和程序产品
CN105740781B (zh) 一种三维人脸活体检测的方法和装置
CN106650669A (zh) 一种鉴别仿冒照片欺骗的人脸识别方法
CN105740779B (zh) 人脸活体检测的方法和装置
CN106874830B (zh) 一种基于rgb-d相机和人脸识别的视障人士辅助方法
CN107480586B (zh) 基于人脸特征点位移的生物识别照片仿冒攻击检测方法
CN108985134A (zh) 基于双目摄像机的人脸活体检测及刷脸交易方法及系统
US8965068B2 (en) Apparatus and method for discriminating disguised face
CN106599660A (zh) 终端安全验证方法及装置
CN104850842B (zh) 移动终端虹膜识别的人机交互方法
KR20170082074A (ko) 관상학적 특징 정보를 이용한 얼굴인식 장치 및 방법
CN108446690A (zh) 一种基于多视角动态特征的人脸活体检测方法
CN109409343A (zh) 一种基于活体检测的人脸识别方法
CN107862298B (zh) 一种基于红外摄像装置下眨眼的活体检测方法
CN106599880A (zh) 一种面向无人监考的同人判别方法
CN112329727A (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

Address after: 511458 room 1011, No.26, Jinlong Road, Nansha District, Guangzhou City, Guangdong Province (only for office use)

Applicant after: Guangzhou yunconghonghuang Intelligent Technology Co., Ltd

Address before: 511457 Room 1011, 26 Jinlong Road, Nansha District, Guangzhou City, Guangdong Province

Applicant before: GUANGZHOU HONGHUANG INTELLIGENT TECHNOLOGY Co.,Ltd.

CB02 Change of applicant information
RJ01 Rejection of invention patent application after publication

Application publication date: 20181228

RJ01 Rejection of invention patent application after publication