CN1945204A - 镜面物体表面三维轮廓测量装置和方法 - Google Patents

镜面物体表面三维轮廓测量装置和方法 Download PDF

Info

Publication number
CN1945204A
CN1945204A CN 200610117315 CN200610117315A CN1945204A CN 1945204 A CN1945204 A CN 1945204A CN 200610117315 CN200610117315 CN 200610117315 CN 200610117315 A CN200610117315 A CN 200610117315A CN 1945204 A CN1945204 A CN 1945204A
Authority
CN
China
Prior art keywords
computing machine
flat
guide rail
image
measuring
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
CN 200610117315
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.)
Shanghai University
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for Science and Technology
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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN 200610117315 priority Critical patent/CN1945204A/zh
Publication of CN1945204A publication Critical patent/CN1945204A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明提供了一种镜面物体表面三维轮廓测量装置和方法,其测量装置由LCD平板显示器、工业摄像机、计算机和导轨组成,LCD平板显示器被固定在导轨上并通过电缆与计算机相连接,工业摄像机通过视频线与计算机相连接。测量方法包括测量装置调整、数据采集和数据处理等步骤。本发明的镜面物体表面三维轮廓测量装置组成简单,测量方法操作方便,数据处理的计算复杂性较低。

Description

镜面物体表面三维轮廓测量装置和方法
技术领域
本发明涉及一种检测技术,尤其涉及一种镜面物体表面三维轮廓测量装置和测量方法。
背景技术
镜面物体表面三维轮廓测量技术有重要的应用价值。在工程中,大量零件的表面性质为镜面反射。例如:抛光模具等精加工零部件;某些表面涂镀零件,如喷涂后的汽车覆盖件等;某些玻璃及塑料制品;印刷线路板焊点等。目前,对于该类零件的三维测量主要采用两种方法。其一,采用坐标测量机或三维测量手臂等接触式测量设备,效率很低;其二,喷涂其表面,改变其性质为漫反射后用光学方法进行测量,最后再除去涂层。显然,该方法严重削弱了光学方法的非接触优点。其它的非接触测量方法包括被动式和主动式两类。其中被动式方法包括“基于光度学的测量技术”和“从镜面反射成分恢复面形技术”,这些技术的数据处理过程复杂,并且精度较低。主动测量主要是采用结构光测量技术。其中,点结构光技术一般采用激光光源,通过电流计扫描部件精确控制出射光线的方向,对目标物体进行逐点扫描,测量其梯度,测量效率很低。采用光带扫描技术部分了提高测量的效率。面结构光技术,根据面结构光的编码方式的不同又可分为图案编码和相位编码两类方法,后者一般称为条纹投射技术。这种技术采用漫射结构光源,将强度为正弦分布的结构光投射到镜面物体表面,通过条纹位相计算被测表面的梯度。这种方法可直接获得全场信息,无需扫描装置,效率很高。上述三种结构光技术都只能测得被测面的梯度分布,必须通过数值积分方能得到被测面的高度信息,计算较为复杂,有较大的累积误差。
发明内容
本发明的目的,在于提供一种通用的、非接触的、以全场方式工作的镜面物体表面三维轮廓测量装置和测量方法,可实现非接触测量,测量效率高,精度高。
为达到上述目的,本发明采用下述技术方案:
一种镜面物体表面三维轮廓测量装置,其特征在于:由LCD平板显示器、工业摄像机、计算机和导轨组成;所述的LCD平板显示器被垂直安装在导轨上,可以沿导轨平移至不同位置,并通过数据电缆与计算机相连接,显示由计算机生成的正弦条纹图案,并将这些条纹图案投射在被测量的物体表面上;所述的工业摄像机通过视频线与计算机相连接,工业摄像机拍摄被测量物体表面上的投射条纹图像,并通过视频线将条纹图像传输到计算机;所述的计算机将条纹图像转换成被测量物体表面轮廓。
一种镜面物体表面三维轮廓测量方法,通过上述的镜面物体表面三维轮廓测量装置实施,其特征在于操作步骤如下:
第一步、测量装置调整:
调整测量装置,确定工业摄像机和LCD平板显示器的参数,以及导轨的空间位置和方向;
第二步、数据采集:
①使平板显示器沿导轨移动到第一个位置,记录该位置;在计算机中生成正弦条纹图像,将这些正弦条纹图像显示在平板显示器上;用摄像机拍摄正弦条纹图像被被测物体表面所成的虚像,并输入计算机;②使平板显示器沿导轨移动到第二个位置,记录该位置;将这些正弦条纹图像显示在平板显示器上;用摄像机拍摄正弦条纹图像被被测物体表面所成的虚像,并输入计算机;
第三步、数据处理:
根据测量系统的几何参数和摄像机采集到的图像信息由计算机计算被测量表面的三维轮廓。
上述镜面物体表面三维轮廓测量方法中,第三步中所述的数据处理方法为:①对平板显示器在两个不同位置时由摄像机采集到的条纹图像进行处理,计算条纹图像的位相分布;②根据图像的像素坐标和摄像机参数确定被测表面的出射光线的方程;③根据像素对应的两个位相值和平板显示器参数确定被测表面的入射光线的方程;④计算出射光线和入射光线交点的坐标,这些坐标就是被测表面点的坐标。
本发明与现有技术相比,具有如下显而易见的突出实质性特点和显著优点:
1、整个测量系统主要由LCD平板显示器、摄像机、计算机和导轨组成,系统结构简单。可实现非接触测量,无需传统测量技术中的各类接触式测头,降低了测量成本,提高了测量效率。
2、采用可移动平板结构光源(平板显示器)替代原有技术中的固定光源,可获得多组物面位相分布,从而实现入射光线方向的直接追踪。原理简单,物面点的空间坐标可以直接计算获得,避免了原有方法中存在的复杂的数值积分过程,同时也避免了累积误差,测量精度高。
附图说明
图1是本发明镜面物体表面三维轮廓测量装置的结构示意图;
图2是本发明镜面物体表面三维轮廓测量方法的原理图;
图3是本发明的一实施例中的被测量凹面镜照片图;
图4是本发明的一实施例中,根据测量结果绘出的等高图。
具体实施方式
本发明的一个优选实例结合附图详述如下:
参见图1,本镜面物体表面三维轮廓测量装置由LCD平板显示器(1)、工业摄像机(2)、计算机(3)和导轨(4)组成。LCD平板显示器(1)被固定在导轨(4)上并通过数据电缆(5)与计算机相(3)连接,可以显示由计算机(3)生成的正弦条纹图案,并将这些条纹图案投射在被测量的物体表面上。工业摄像机(2)通过视频线(6)与计算机(3)相连接,可拍摄被测量物体表面上的投射条纹图像,并通过视频线(6)将条纹图像传输到计算机(3)。LCD平板显示器(1)被垂直安装在导轨(4)上,可以沿导轨(4)平移至不同位置。计算机内(3)设有将条纹图像转换成被测量物体表面轮廓的软件。
本镜面物体表面三维轮廓测量方法的测量原理如图2所示。点C是摄像机的入瞳中心。摄像机光轴与空间一个平面即参考面R相交于点O。以O点为原点建立坐标系,其中X轴与图面平行,Z轴与图面垂直,Y轴与摄像机光轴重合。LCD平板显示器在导轨驱动下可沿着直线OQ运动。OQ垂直于LCD平板显示器表面。LCD平板显示器在运动前后,其表面分别位于OAA和OBB。直线OQ过LCD中心OA,OB。对于摄像机每一像素(i,j),其接收的物面反射光线为直线CR。入射光线可通过LCD在不同位置时获取的相位信息来确定。LCD在两个位置OAA,OBB两个位置时,点(i,j)处条纹图像分别来自点A和B,其相位分别为φA和φB。通过φA和φB,可计算A与B的位置坐标。直线AB就是CR所对应的物面入射光线。其交点坐标即为物面点坐标。如此,遍历图像上的每一个像素点,就可以重建摄像机视野内的三维面形。
本镜面物体表面三维轮廓测量方法的测量过程包括:第一步、调整测量装置:确定工业摄像机和LCD平板显示器的参数,以及导轨的空间位置和方向。第二步、进行数据采集:①使平板显示器沿导轨移动到第一个位置,记录该位置;在计算机中生成正弦条纹图像,将这些正弦条纹图像显示在平板显示器上;用摄像机拍摄正弦条纹图像被被测物体表面所成的虚像,并输入计算机;②使平板显示器沿导轨移动到第二个位置,记录该位置;将这些正弦条纹图像显示在平板显示器上;用摄像机拍摄正弦条纹图像被被测物体表面所成的虚像,并输入计算机。第三步、进行数据处理:根据测量系统的几何参数和摄像机采集到的图像信息由计算机计算被测量表面的三维轮廓。首先,利用下列公式计算平板显示器在两个不同位置时由摄像机采集到的条纹图像的位相分布。
φ ( i , j ) = - tan - 1 [ Σ i = 0 N - 1 I k sin α k Σ k = 0 N - 1 I k cos α k ]
式中(i,j)是像素的坐标。Ik是图像的灰度值。αk=2kπ/N(k=0,1,ΛN-1)是相移量,由计算机控制。然后利用下列公式计算被测表面点的坐标(x,y,z)。
x = x R [ h ( L 1 - L 2 ) sin θ + K ( h cos θ - L 1 ) φ B + K ( L 2 - h cos θ ) φ A ( x R cos θ - h sin θ ) ( L 2 - L 1 ) + K ( φ B - φ A ) ( x R sin θ + h cos θ ) ]
y = h [ K ( x R sin θ + L 1 ) · φ B - K ( x R sin θ + L 2 ) · φ A + x R ( L 2 - L 1 ) · cos θ ( x R cos θ - h sin θ ) ( L 2 - L 1 ) + K ( φ B - φ A ) ( x R sin θ + h cos θ ) ]
z = z R [ h ( L 1 - L 2 ) sin θ + K ( h cos θ - L 1 ) φ B + K ( L 2 - h cos θ ) φ A ( x R cos θ - h sin θ ) ( L 2 - L 1 ) + K ( φ B - φ A ) ( x R sin θ + h cos θ ) ]
式中,h是摄像机入瞳中心到参考面的距离。θ是直线OQ与Y轴的夹角。L1和L2分别是LCD显示器沿着直线OQ平移到两个投射位置时,中心OA和OB与原点之间的距离。K=P/(2π)是相位转换系数(P是投射条纹的节距)。xR=Si,yR=Sj(S是摄像机放大因子)。这些参数均已在测量装置调整时被确定。
以下通过一个具体实例对本发明镜面物体表面三维轮廓测量装置和测量方法的测量原理和测量过程作进一步说明:图3是一个被测凹面镜,其直径为113mm。调整测量装置,使其满足R=0.3884mm/pixel,h=771.2766mm。在L=178mm和228mm两个位置投射两组条纹正弦条纹。求解位相,再计算被测物面点的坐标。图4是根据测量结果画出的等高图。

Claims (3)

1.一种镜面物体表面三维轮廓测量装置,其特征在于:由LCD平板显示器(1)、工业摄像机(2)、计算机(3)和导轨(4)组成;所述的LCD平板显示器(1)被垂直安装在导轨(4)上,可以沿导轨(4)平移至不同位置,并通过数据电缆(5)与计算机相(3)连接,显示由计算机(3)生成的正弦条纹图案,并将这些条纹图案投射在被测量的物体表面上;所述的工业摄像机(2)通过视频线(6)与计算机(3)相连接,工业摄像机(2)拍摄被测量物体表面上的投射条纹图像,并通过视频线(6)将条纹图像传输到计算机(3);所述的计算机(3)将条纹图像转换成被测量物体表面轮廓。
2.一种镜面物体表面三维轮廓测量方法,通过根据权利要求1所述的镜面物体表面三维轮廓测量装置实施,其特征在于操作步骤如下:
第一步、测量装置调整:调整测量装置的几何参数,这些参数包括工业摄像机(2)和导轨(4)的空间位置和方向;
第二步、数据采集:①使平板显示器(1)沿导轨移动到第一个位置,记录该位置;在计算机(3)中生成正弦条纹图像,将这些正弦条纹图像显示在平板显示器(1)上;用摄像机(2)拍摄正弦条纹图像被被测物体表面所成的虚像,并输入计算机(3);②使平板显示器(1)沿导轨(4)移动到第二个位置,记录该位置;将这些正弦条纹图像显示在平板显示器(1)上;用摄像机(2)拍摄正弦条纹图像被被测物体表面所成的虚像,并输入计算机(3);
第三步、数据处理:根据测量系统的几何参数和摄像机(2)采集到的图像信息,由计算机计算被测量表面的三维轮廓。
3.根据权利要求2所述的镜面物体表面三维轮廓测量方法,其特征在于所述的第三步中所述的数据处理方法为:①对平板显示器(1)在两个不同位置时由摄像机(2)采集到的条纹图像进行处理,计算条纹图像的位相分布;②根据图像的像素坐标和摄像机(2)参数确定被测表面的出射光线的方程;③根据像素对应的两个位相值和平板显示器(1)参数确定被测表面的入射光线的方程;④计算出射光线和入射光线交点的坐标,这些坐标就是被测表面点的坐标。
CN 200610117315 2006-10-19 2006-10-19 镜面物体表面三维轮廓测量装置和方法 Pending CN1945204A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200610117315 CN1945204A (zh) 2006-10-19 2006-10-19 镜面物体表面三维轮廓测量装置和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200610117315 CN1945204A (zh) 2006-10-19 2006-10-19 镜面物体表面三维轮廓测量装置和方法

Publications (1)

Publication Number Publication Date
CN1945204A true CN1945204A (zh) 2007-04-11

Family

ID=38044708

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200610117315 Pending CN1945204A (zh) 2006-10-19 2006-10-19 镜面物体表面三维轮廓测量装置和方法

Country Status (1)

Country Link
CN (1) CN1945204A (zh)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102901465A (zh) * 2012-10-29 2013-01-30 大连宏海新能源发展有限公司 一种测量光学镜片斜率误差的设备和方法
CN103217126A (zh) * 2013-04-24 2013-07-24 中国科学院电工研究所 一种太阳能槽式聚光器面形检测系统及方法
CN103487441A (zh) * 2013-09-24 2014-01-01 电子科技大学 一种用于硅晶片缺陷检测和面形测量的方法
CN103630541A (zh) * 2013-11-28 2014-03-12 广州市香港科大霍英东研究院 一种双光源镜面反射面测量系统
CN104101611A (zh) * 2014-06-06 2014-10-15 华南理工大学 一种类镜面物体表面光学成像装置及其成像方法
CN104111036A (zh) * 2013-04-18 2014-10-22 中国科学院沈阳自动化研究所 一种基于双目视觉的镜面物体测量装置及方法
CN104180773A (zh) * 2013-05-20 2014-12-03 沈阳同联集团高新技术有限公司 一种三维形貌测量装置
CN104266607A (zh) * 2014-09-22 2015-01-07 电子科技大学 镜面目标轮廓光学测量系统和方法
CN104279981A (zh) * 2014-10-21 2015-01-14 苏州大学 一种基于条纹反射的镜面/类镜面物体绝对面形的测量方法及装置
CN104613898A (zh) * 2013-10-12 2015-05-13 石泉 检测仪、具有该检测仪的测量系统以及测量方法
CN105091784A (zh) * 2015-06-30 2015-11-25 东莞市盟拓光电科技有限公司 用于表面为镜面或透明表面的被测物的三维成像系统
CN105764860A (zh) * 2013-11-25 2016-07-13 康宁股份有限公司 用于确定基本上圆柱形镜面反射表面的形状的方法
CN105806257A (zh) * 2016-03-12 2016-07-27 上海大学 一种高反射物体表面光场偏折术测量系统和方法
CN106216815A (zh) * 2016-09-21 2016-12-14 兰州理工大学 一种基于双屏的物体表面三维形貌测量方法
CN106546192A (zh) * 2016-10-12 2017-03-29 上海大学 一种高反射自由曲面非接触测量方法和系统
CN108627121A (zh) * 2018-05-15 2018-10-09 浙江中控太阳能技术有限公司 一种镜面面形检测装置及其检测方法
CN109855573A (zh) * 2018-12-29 2019-06-07 中国科学院长春光学精密机械与物理研究所 一种类镜面多尺度表面形貌特征检测系统及其检测方法
CN111656260A (zh) * 2018-01-24 2020-09-11 赛博光学公司 用于镜像表面的结构式光投射
CN111879258A (zh) * 2020-09-28 2020-11-03 南京理工大学 基于条纹图像转换网络FPTNet的动态高精度三维测量方法
CN113237436A (zh) * 2021-04-07 2021-08-10 上海大学 偏振光空间相移非朗伯金属物体形貌测量装置及方法
TWI735953B (zh) * 2019-09-18 2021-08-11 財團法人工業技術研究院 三維量測裝置與其操作方法

Cited By (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102901465A (zh) * 2012-10-29 2013-01-30 大连宏海新能源发展有限公司 一种测量光学镜片斜率误差的设备和方法
CN104111036A (zh) * 2013-04-18 2014-10-22 中国科学院沈阳自动化研究所 一种基于双目视觉的镜面物体测量装置及方法
CN103217126A (zh) * 2013-04-24 2013-07-24 中国科学院电工研究所 一种太阳能槽式聚光器面形检测系统及方法
CN103217126B (zh) * 2013-04-24 2016-09-14 中国科学院电工研究所 一种太阳能槽式聚光器面形检测系统及方法
CN104180773A (zh) * 2013-05-20 2014-12-03 沈阳同联集团高新技术有限公司 一种三维形貌测量装置
CN103487441A (zh) * 2013-09-24 2014-01-01 电子科技大学 一种用于硅晶片缺陷检测和面形测量的方法
CN104613898B (zh) * 2013-10-12 2018-04-13 石泉 检测仪、具有该检测仪的测量系统以及测量方法
CN104613898A (zh) * 2013-10-12 2015-05-13 石泉 检测仪、具有该检测仪的测量系统以及测量方法
CN105764860A (zh) * 2013-11-25 2016-07-13 康宁股份有限公司 用于确定基本上圆柱形镜面反射表面的形状的方法
CN105764860B (zh) * 2013-11-25 2019-01-11 康宁股份有限公司 用于确定基本上圆柱形镜面反射表面的形状的方法
CN103630541A (zh) * 2013-11-28 2014-03-12 广州市香港科大霍英东研究院 一种双光源镜面反射面测量系统
CN103630541B (zh) * 2013-11-28 2016-05-25 广州市香港科大霍英东研究院 一种双光源镜面反射面测量系统
CN104101611A (zh) * 2014-06-06 2014-10-15 华南理工大学 一种类镜面物体表面光学成像装置及其成像方法
CN104266607A (zh) * 2014-09-22 2015-01-07 电子科技大学 镜面目标轮廓光学测量系统和方法
CN104266607B (zh) * 2014-09-22 2017-03-29 电子科技大学 镜面目标轮廓光学测量方法
CN104279981A (zh) * 2014-10-21 2015-01-14 苏州大学 一种基于条纹反射的镜面/类镜面物体绝对面形的测量方法及装置
CN104279981B (zh) * 2014-10-21 2017-07-07 苏州大学 一种基于条纹反射的镜面/类镜面物体绝对面形的测量方法及装置
CN105091784A (zh) * 2015-06-30 2015-11-25 东莞市盟拓光电科技有限公司 用于表面为镜面或透明表面的被测物的三维成像系统
CN105806257A (zh) * 2016-03-12 2016-07-27 上海大学 一种高反射物体表面光场偏折术测量系统和方法
CN105806257B (zh) * 2016-03-12 2019-02-22 上海大学 一种高反射物体表面光场偏折术测量系统和方法
CN106216815A (zh) * 2016-09-21 2016-12-14 兰州理工大学 一种基于双屏的物体表面三维形貌测量方法
CN106216815B (zh) * 2016-09-21 2018-07-03 兰州理工大学 一种基于双屏的物体表面三维形貌测量方法
CN106546192A (zh) * 2016-10-12 2017-03-29 上海大学 一种高反射自由曲面非接触测量方法和系统
CN106546192B (zh) * 2016-10-12 2019-08-06 上海大学 一种高反射自由曲面非接触测量方法和系统
CN111656260A (zh) * 2018-01-24 2020-09-11 赛博光学公司 用于镜像表面的结构式光投射
CN108627121A (zh) * 2018-05-15 2018-10-09 浙江中控太阳能技术有限公司 一种镜面面形检测装置及其检测方法
CN109855573A (zh) * 2018-12-29 2019-06-07 中国科学院长春光学精密机械与物理研究所 一种类镜面多尺度表面形貌特征检测系统及其检测方法
TWI735953B (zh) * 2019-09-18 2021-08-11 財團法人工業技術研究院 三維量測裝置與其操作方法
US11248903B2 (en) 2019-09-18 2022-02-15 Industrial Technology Research Institute Three-dimension measurement device and operation method thereof
CN111879258A (zh) * 2020-09-28 2020-11-03 南京理工大学 基于条纹图像转换网络FPTNet的动态高精度三维测量方法
CN113237436A (zh) * 2021-04-07 2021-08-10 上海大学 偏振光空间相移非朗伯金属物体形貌测量装置及方法
CN113237436B (zh) * 2021-04-07 2022-08-02 上海大学 偏振光空间相移非朗伯金属物体形貌测量装置及方法

Similar Documents

Publication Publication Date Title
CN1945204A (zh) 镜面物体表面三维轮廓测量装置和方法
US8180101B2 (en) Calibration method for structure parameters of structured-light vision sensor
US10812694B2 (en) Real-time inspection guidance of triangulation scanner
Zhou et al. Complete calibration of a structured light stripe vision sensor through planar target of unknown orientations
JP3511450B2 (ja) 光学式測定装置の位置校正方法
CN203231736U (zh) 一种基于双目视觉的镜面物体测量装置
JPH09105613A (ja) 非接触型の三次元測定装置および測定方法
AU2704595A (en) Method and apparatus for transforming coordinate systems in an automated video monitor alignment system
CN102288131A (zh) 物体360°轮廓误差的自适应条纹测量装置和方法
CN101813462A (zh) 单处理器控制的三维形貌光学测量系统及测量方法
CN106168466A (zh) 全域式影像检测系统及其检测方法
CN111811433B (zh) 基于红蓝正交条纹的结构光系统标定方法及装置和应用
CN111765851B (zh) 基于双屏透射显示的镜面物体三维形貌测量方法及装置
Tutsch et al. Optical three-dimensional metrology with structured illumination
CN113888626A (zh) 惯性测量单元与线扫描3d相机在线联合标定装置及方法
CN1267699C (zh) 基于相移的光栅投影式三维轮廓测量装置及测量方法
CN105841618A (zh) 二维三维复合式测量仪及其数据融合标定方法
Yang et al. A novel projector ray-model for 3D measurement in fringe projection profilometry
CN1093935C (zh) 一种快速投影结构光的三维轮廓相位测量方法及装置
CN115855955A (zh) 一种基于多束线激光的模具表面结构缺陷检测装置和方法
CN113899321A (zh) 凹面镜辅助成像的镜面物体三维形貌测量方法及系统
JPH09218022A (ja) ワークの拡散面の輪郭決定方法
Xu et al. A robot-assisted back-imaging measurement system for transparent glass
CN111397511A (zh) 一种利用物体平移进行单目三维测量的方法及装置
CN110966937A (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
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication