CN109668543A - 基于激光雷达的倾斜度测量方法 - Google Patents
基于激光雷达的倾斜度测量方法 Download PDFInfo
- Publication number
- CN109668543A CN109668543A CN201910057998.1A CN201910057998A CN109668543A CN 109668543 A CN109668543 A CN 109668543A CN 201910057998 A CN201910057998 A CN 201910057998A CN 109668543 A CN109668543 A CN 109668543A
- Authority
- CN
- China
- Prior art keywords
- laser radar
- under
- coordinate system
- polar
- measured
- 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
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C9/00—Measuring inclination, e.g. by clinometers, by levels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
Abstract
本发明提出了一种基于激光雷达的倾斜度测量方法,首先,将激光雷达正对待测量平面设置,激光雷达沿水平方向测量待测量平面获取激光雷达极坐标系下的角度信息θ0,以激光雷达为极点、水平方向为极轴重新构建极坐标系W;然后,改变激光雷达的照射角度,使激光雷达沿垂直方向旋转扫描获取待测量平面在此垂直方向上的点的激光雷达极坐标系下极坐标,获得对应极坐标系W下的极坐标(ρi,βi);最后,将极坐标W下的点转换为直角坐标并将各直角坐标下的点连接成线,该线的斜率即为墙面倾斜度。本发明通过连线的斜率求出墙面倾斜度,计算更加精确,同时,设备简单,对环境要求不高。
Description
技术领域
本发明属于计量检测技术领域,具体为一种基于激光雷达的倾斜度测量方法。
背景技术
建筑工程项目是一个系统化复杂的施工过程,是从设计到施工和管理的一个系统化施工流程和施工管理模式,是一项涉及到多个学科和因素的过程,在施工中针对各个施工环节进行综合控制是确保施工管理和制度合理有效利用的基础。然而随着目前各种自然因素的影响,建筑物结构稳定性与安全也受到严重考验,因此对建筑物墙体倾斜度进行测量对于保障建筑结构整体性与建筑安全有着重要作用与意义。
全站仪也被称之为全站型电子速测仪,这种测量措施与测量方法是在传统的测量中人们已经提到过速测法的应用环节,主要是从仪器站同时测定某一点的平面位置和高层的施工应用方法,在有些时候,这种测试措施与方法也被称之为速测术、速测法。而速测仪则是这种测量方式中不可缺少的一部分,是根据速测法原理而设计的一种新型的测量模式和测量设备。最初的建筑倾斜度测量过程中常用的速测仪的距离是通过光学方法来实现的,我们也称之为光学速测仪,这种速测方式是通过带有视距丝的经纬仪以及相关需要测定的平面位置为基础记性确定,然后再根据高层和三角测量方法来进行深入总结和分析的过程。这种倾斜度测量的方法比较复杂,而且测量起来也稍显麻烦,精度也不高。
现在,激光三角法也被应用到倾斜度测量中。通过激光测距和墙面的高度差来求得墙面的倾斜度。可这种方法对于距离稍远的墙体来说就存在着标定的问题。对于远距离的标定,方法会复杂些,对于整个倾斜度测量的耗时会比较久。对于一些需要快速得到倾斜度的应用场景,这些方法都稍显麻烦。
发明内容
本发明的目的在于提出了一种基于激光雷达的倾斜度测量方法。
实现本发明的技术解决方案为:一种基于激光雷达的倾斜度测量方法,具体步骤为:
步骤1、将激光雷达正对待测量平面设置,激光雷达沿水平方向测量待测量平面获取激光雷达极坐标系下的角度信息θ0,以激光雷达为极点、水平方向为极轴重新构建极坐标系W,构建的极坐标系W下的角度坐标与激光雷达极坐标系下的角度坐标之间的关系为β=θ-θ0;
步骤2、改变激光雷达的照射角度,使激光雷达沿垂直方向旋转扫描获取待测量平面在此垂直方向上的点的激光雷达极坐标系下极坐标(ρi,θi),获得对应极坐标系W下的极坐标(ρi,βi),其中βi=θi-θ0;
步骤3、将极坐标W下的点转换为直角坐标并将各直角坐标下的点连接成线,该线的斜率即为墙面倾斜度。
优选地,将极坐标变化为直角坐标的公式为:
xi=ρi*cosβi
yi=ρi*sinβi
式中,ρi为激光雷达到待测量平面被测点的距离。
优选地,步骤3中各直角坐标连接成的线为直线或曲线,曲线的斜率通过均值法获得。
本发明与现有技术相比,其显著优点为:(1)本发明使用激光雷达扫描,获取墙面信息速度快;(2)本发明将墙面上接受垂直扫描的区域与激光雷达形成的极坐标式的位置关系转化为直角坐标系下的坐标关系,并用连线还原墙面,更加直观;(3)本发明通过连线的斜率求出墙面倾斜度,计算更加精确,同时,设备简单,对环境要求不高。
下面结合附图对本发明做进一步详细的描述。
附图说明
图1是基于激光雷达的倾斜度测量方法的实物结构图。
图2是激光雷达工作原理示意图。
图3是基于激光雷达的倾斜度测量方法的雷达扫描信息图。
图4是基于激光雷达的倾斜度测量方法的直角坐标系下的信息图。
具体实施方式
一种基于激光雷达的倾斜度测量方法,具体步骤为:
步骤1、将激光雷达正对待测量平面设置,激光雷达沿水平方向测量待测量平面获取激光雷达极坐标系下的角度信息θ0,以激光雷达为极点、水平方向为极轴重新构建极坐标系W,即此时起始角度为0,根据获得的θ0,则此时构构建的极坐标系W下的角度坐标与激光雷达极坐标系下的角度坐标之间的关系为β=θ-θ0;
步骤2、改变激光雷达的照射角度,使激光雷达沿垂直方向旋转扫描获取待测量平面在此垂直方向上的点的激光雷达极坐标系下极坐标(ρi,θi),获得对应极坐标系W下的极坐标(ρi,βi),其中βi=θi-θ0;
步骤3、将极坐标W下的点变化为直角坐标并将各直角坐标下的点连接成线,该线的斜率即为墙面倾斜度。
进一步的实施例中,将极坐标变化为直角坐标的公式为:
xi=ρi*cosβi
yi=ρi*sinβi
式中,ρi为激光雷达到待测量平面被测点的距离。
进一步的实施例中,步骤3中各直角坐标连接成的线为直线或曲线,曲线的斜率通过均值法获得。
实施例
如图1所示,实现本实施例的硬件结构主要分为两个模块:雷达扫描模块和后端处理模块。雷达扫描模块是由激光雷达和机架组成。将激光雷达正对待测量平面设置,激光雷达沿水平方向开始测量待测量平面获取待测量平面在垂直方向上的极坐标信息。后端处理模块使用ARM开发板和LINUX系统完成,主要是将激光雷达获取的距离信号和极坐标角度转换为直角坐标系下的二维信息。
如图2所示,通过激光雷达获取墙面信息,所用的激光雷达为RPLIDAR A1,采用三角测距法来获取墙面信息点。通过激光雷达扫描待测平面可获得如图3所示的位置信息关系,具体步骤为:
步骤1、对激光雷达进行标定,激光雷达沿水平方向测量待测量平面获取激光雷达极坐标系下的角度信息θ0,以激光雷达为极点、水平方向为极轴重新构建极坐标系W,构建的极坐标系W下的角度坐标与激光雷达极坐标系下的角度坐标之间的关系为β=θ-θ0。
步骤2、改变激光雷达的照射角度,使激光雷达沿垂直方向旋转扫描获取待测量平面在此垂直方向上的点的激光雷达极坐标系下极坐标(ρi,θi),获得对应极坐标系W下的极坐标(ρi,βi),其中βi=θi-θ0;其中,将极坐标变化为直角坐标的公式为:
xi=ρi*cosβi
yi=ρi*sinβi
式中,ρi为激光雷达到待测量平面被测点的距离。
步骤3、如图4所示,将极坐标W下的点变化为直角坐标并将各直角坐标下的点连接成线,该线的斜率即为墙面倾斜度。获得的直角坐标下的墙面被扫描部分的关系曲线不一定是直线,由于扫描误差或者墙面本身的不平滑,可能得到的是不平滑的曲线。通过均值法求出这个曲线对应的直线的斜率,即为所测墙面的倾斜率。直线斜率即为墙面倾斜度。
Claims (3)
1.一种基于激光雷达的倾斜度测量方法,其特征在于,具体步骤为:
步骤1、将激光雷达正对待测量平面设置,激光雷达沿水平方向测量待测量平面获取激光雷达极坐标系下的角度信息θ0,以激光雷达为极点、水平方向为极轴重新构建极坐标系W,构建的极坐标系W下的角度坐标与激光雷达极坐标系下的角度坐标之间的关系为β=θ-θ0;
步骤2、改变激光雷达的照射角度,使激光雷达沿垂直方向旋转扫描获取待测量平面在此垂直方向上的点的激光雷达极坐标系下极坐标(ρi,θi),获得对应极坐标系W下的极坐标(ρi,βi),其中βi=θi-θ0;
步骤3、将极坐标W下的点转换为直角坐标并将各直角坐标下的点连接成线,该线的斜率即为墙面倾斜度。
2.根据权利要求1所述的基于激光雷达的倾斜度测量方法,其特征在于,将极坐标变化为直角坐标的公式为:
xi=ρi*cosβi
yi=ρi*sinβi
式中,ρi为激光雷达到待测量平面被测点的距离。
3.根据权利要求1所述的基于激光雷达的倾斜度测量方法,其特征在于,步骤3中各直角坐标连接成的线为直线或曲线,曲线的斜率通过均值法获得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910057998.1A CN109668543A (zh) | 2019-01-22 | 2019-01-22 | 基于激光雷达的倾斜度测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910057998.1A CN109668543A (zh) | 2019-01-22 | 2019-01-22 | 基于激光雷达的倾斜度测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109668543A true CN109668543A (zh) | 2019-04-23 |
Family
ID=66149660
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910057998.1A Pending CN109668543A (zh) | 2019-01-22 | 2019-01-22 | 基于激光雷达的倾斜度测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109668543A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111060919A (zh) * | 2019-12-30 | 2020-04-24 | 广东博智林机器人有限公司 | 基于智能机器人的测量方法和智能机器人 |
CN111158012A (zh) * | 2019-12-26 | 2020-05-15 | 大庆思特传媒科技有限公司 | 弧形墙面互动识别系统和识别方法 |
CN111208493A (zh) * | 2020-01-08 | 2020-05-29 | 同济大学 | 一种车载激光雷达在整车坐标系下的快速标定方法 |
CN113064145A (zh) * | 2021-03-24 | 2021-07-02 | 盎锐(上海)信息科技有限公司 | 基于激光雷达的水平标定方法、系统及激光雷达 |
CN113311422A (zh) * | 2020-02-27 | 2021-08-27 | 富士通株式会社 | 坐标转换方法、装置和数据处理设备 |
CN113701687A (zh) * | 2021-09-17 | 2021-11-26 | 中国联合网络通信集团有限公司 | 空间曲线长度测量系统及测量方法 |
Citations (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4311384A (en) * | 1975-08-04 | 1982-01-19 | Raytheon Company | Laser scanning system |
US20030095249A1 (en) * | 2000-04-11 | 2003-05-22 | Pierre Ferdinand | Inclinometer with bragg grating |
CN103017734A (zh) * | 2012-12-11 | 2013-04-03 | 湖北省电力公司检修分公司 | 一种基于激光雷达的杆塔倾斜率测量方法 |
EP2719475A1 (en) * | 2012-10-12 | 2014-04-16 | Nivora IP B.V. | Measurement system and method for measuring an angle |
CN103925872A (zh) * | 2013-12-23 | 2014-07-16 | 中国神华能源股份有限公司 | 一种用于获取目标物空间分布的激光扫描测量装置和方法 |
CN104950313A (zh) * | 2015-06-11 | 2015-09-30 | 同济大学 | 一种路面提取及道路坡度识别方法 |
CN105044039A (zh) * | 2015-07-24 | 2015-11-11 | 北京环境特性研究所 | 一种根据激光雷达数据自动反演水平能见度的方法 |
CN105241395A (zh) * | 2015-08-07 | 2016-01-13 | 上海卫星工程研究所 | 用于空间超大尺寸结构的面形与距离高精度测量方法 |
CN105333861A (zh) * | 2015-12-02 | 2016-02-17 | 中国测绘科学研究院 | 基于激光点云的杆塔倾斜检测的方法及装置 |
US20160076877A1 (en) * | 2014-09-12 | 2016-03-17 | Mitutoyo Corporation | Position measuring device |
US20160131479A1 (en) * | 2014-11-12 | 2016-05-12 | Kabushiki Kaisha Topcon | Tilt Detecting System And Tilt Detecting Method |
CN106019297A (zh) * | 2016-07-29 | 2016-10-12 | 成都希德电子信息技术有限公司 | 一种2d激光雷达扫描系统及其装置 |
CN106020180A (zh) * | 2015-03-26 | 2016-10-12 | 株式会社日立制作所 | 倾斜度检测方法、倾斜度检测装置和倾斜度检测用器材 |
CN106500663A (zh) * | 2017-01-06 | 2017-03-15 | 中国人民解放军火箭军工程大学 | 一种倾斜度激光测量方法及系统 |
CN106646508A (zh) * | 2016-11-24 | 2017-05-10 | 中国科学院自动化研究所 | 面向斜坡区域的基于多线激光雷达的斜坡角度估计方法 |
CN106885531A (zh) * | 2017-04-20 | 2017-06-23 | 河北科技大学 | 基于二维激光雷达的车斗形容器三维扫描系统标定方法 |
CN107101579A (zh) * | 2017-04-26 | 2017-08-29 | 长沙迪迈数码科技股份有限公司 | 一种采空区点云坐标修正方法 |
CN107388992A (zh) * | 2017-07-26 | 2017-11-24 | 中国电建集团西北勘测设计研究院有限公司 | 一种基于三维激光扫描的高耸塔筒垂直度检测方法 |
CN107607929A (zh) * | 2017-09-20 | 2018-01-19 | 云南电网有限责任公司电力科学研究院 | 一种基于激光点云数据的测量杆塔倾斜角度的方法及装置 |
CN108389250A (zh) * | 2018-03-08 | 2018-08-10 | 武汉大学 | 基于点云数据快速生成建筑物断面图的方法 |
CN108445505A (zh) * | 2018-03-29 | 2018-08-24 | 南京航空航天大学 | 线环境下基于激光雷达的特征显著性检测方法 |
EP3382335A1 (en) * | 2017-03-31 | 2018-10-03 | Pliant Holding B.V. | Measuring attitude of constructions |
CN109239721A (zh) * | 2018-09-04 | 2019-01-18 | 广东电网有限责任公司 | 一种输电线路杆塔倾斜度自动测量系统及其测量方法 |
-
2019
- 2019-01-22 CN CN201910057998.1A patent/CN109668543A/zh active Pending
Patent Citations (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4311384A (en) * | 1975-08-04 | 1982-01-19 | Raytheon Company | Laser scanning system |
US20030095249A1 (en) * | 2000-04-11 | 2003-05-22 | Pierre Ferdinand | Inclinometer with bragg grating |
EP2719475A1 (en) * | 2012-10-12 | 2014-04-16 | Nivora IP B.V. | Measurement system and method for measuring an angle |
CN103017734A (zh) * | 2012-12-11 | 2013-04-03 | 湖北省电力公司检修分公司 | 一种基于激光雷达的杆塔倾斜率测量方法 |
CN103925872A (zh) * | 2013-12-23 | 2014-07-16 | 中国神华能源股份有限公司 | 一种用于获取目标物空间分布的激光扫描测量装置和方法 |
US20160076877A1 (en) * | 2014-09-12 | 2016-03-17 | Mitutoyo Corporation | Position measuring device |
US20160131479A1 (en) * | 2014-11-12 | 2016-05-12 | Kabushiki Kaisha Topcon | Tilt Detecting System And Tilt Detecting Method |
CN106020180A (zh) * | 2015-03-26 | 2016-10-12 | 株式会社日立制作所 | 倾斜度检测方法、倾斜度检测装置和倾斜度检测用器材 |
CN104950313A (zh) * | 2015-06-11 | 2015-09-30 | 同济大学 | 一种路面提取及道路坡度识别方法 |
CN105044039A (zh) * | 2015-07-24 | 2015-11-11 | 北京环境特性研究所 | 一种根据激光雷达数据自动反演水平能见度的方法 |
CN105241395A (zh) * | 2015-08-07 | 2016-01-13 | 上海卫星工程研究所 | 用于空间超大尺寸结构的面形与距离高精度测量方法 |
CN105333861A (zh) * | 2015-12-02 | 2016-02-17 | 中国测绘科学研究院 | 基于激光点云的杆塔倾斜检测的方法及装置 |
CN106019297A (zh) * | 2016-07-29 | 2016-10-12 | 成都希德电子信息技术有限公司 | 一种2d激光雷达扫描系统及其装置 |
CN106646508A (zh) * | 2016-11-24 | 2017-05-10 | 中国科学院自动化研究所 | 面向斜坡区域的基于多线激光雷达的斜坡角度估计方法 |
CN106500663A (zh) * | 2017-01-06 | 2017-03-15 | 中国人民解放军火箭军工程大学 | 一种倾斜度激光测量方法及系统 |
EP3382335A1 (en) * | 2017-03-31 | 2018-10-03 | Pliant Holding B.V. | Measuring attitude of constructions |
CN106885531A (zh) * | 2017-04-20 | 2017-06-23 | 河北科技大学 | 基于二维激光雷达的车斗形容器三维扫描系统标定方法 |
CN107101579A (zh) * | 2017-04-26 | 2017-08-29 | 长沙迪迈数码科技股份有限公司 | 一种采空区点云坐标修正方法 |
CN107388992A (zh) * | 2017-07-26 | 2017-11-24 | 中国电建集团西北勘测设计研究院有限公司 | 一种基于三维激光扫描的高耸塔筒垂直度检测方法 |
CN107607929A (zh) * | 2017-09-20 | 2018-01-19 | 云南电网有限责任公司电力科学研究院 | 一种基于激光点云数据的测量杆塔倾斜角度的方法及装置 |
CN108389250A (zh) * | 2018-03-08 | 2018-08-10 | 武汉大学 | 基于点云数据快速生成建筑物断面图的方法 |
CN108445505A (zh) * | 2018-03-29 | 2018-08-24 | 南京航空航天大学 | 线环境下基于激光雷达的特征显著性检测方法 |
CN109239721A (zh) * | 2018-09-04 | 2019-01-18 | 广东电网有限责任公司 | 一种输电线路杆塔倾斜度自动测量系统及其测量方法 |
Non-Patent Citations (1)
Title |
---|
宋晶 等: "《一种基于车载三维激光扫描技术的电杆倾斜巡查方法》", 《贵州大学学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111158012A (zh) * | 2019-12-26 | 2020-05-15 | 大庆思特传媒科技有限公司 | 弧形墙面互动识别系统和识别方法 |
CN111060919A (zh) * | 2019-12-30 | 2020-04-24 | 广东博智林机器人有限公司 | 基于智能机器人的测量方法和智能机器人 |
CN111208493A (zh) * | 2020-01-08 | 2020-05-29 | 同济大学 | 一种车载激光雷达在整车坐标系下的快速标定方法 |
CN113311422A (zh) * | 2020-02-27 | 2021-08-27 | 富士通株式会社 | 坐标转换方法、装置和数据处理设备 |
CN113064145A (zh) * | 2021-03-24 | 2021-07-02 | 盎锐(上海)信息科技有限公司 | 基于激光雷达的水平标定方法、系统及激光雷达 |
CN113064145B (zh) * | 2021-03-24 | 2024-03-29 | 盎锐(杭州)信息科技有限公司 | 基于激光雷达的水平标定方法、系统及激光雷达 |
CN113701687A (zh) * | 2021-09-17 | 2021-11-26 | 中国联合网络通信集团有限公司 | 空间曲线长度测量系统及测量方法 |
CN113701687B (zh) * | 2021-09-17 | 2023-07-21 | 中国联合网络通信集团有限公司 | 空间曲线长度测量系统及测量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109668543A (zh) | 基于激光雷达的倾斜度测量方法 | |
CN107121123B (zh) | 卫星精度单机测量方法 | |
CN108051835B (zh) | 一种基于双天线的倾斜测量装置及测量与放样方法 | |
CN103292748B (zh) | 一种基于激光测量的多基板拼合检测方法 | |
CN102589488B (zh) | 光学标尺检测平面度和倾斜度的方法 | |
JP5538929B2 (ja) | 三次元位置計測及び墨出しシステムとその使用方法 | |
CN105606129B (zh) | 一种辅助飞机惯导成品组件安装的测量校准方法 | |
CN107036557B (zh) | 一种二维测角系统和方法 | |
CN102322896B (zh) | 隧道tsp几何参数全方位测量装置与方法 | |
CZ2014489A3 (cs) | Způsob měření, zpracování a využívání dat digitálního modelu terénu pro objektivní hodnocení geometrických parametrů měřených objektů a měřící zařízení pro provádění uvedeného způsobu | |
CN110017823B (zh) | 基于全站仪免置平设站的轨道定位坐标计算方法及装置 | |
CN106017319A (zh) | 一种基于高精度点位测量的三维激光扫描数据坐标转换工具及方法 | |
KR20120051981A (ko) | Gnss를 이용한 건물의 수직도 측정 방법 및 시스템 | |
CN113267156B (zh) | 一种利用惯导测量水泥地坪平整度的方法及测量系统 | |
CN103292795A (zh) | 一种双棱镜精确测量的方法及测量杆 | |
CN114858140B (zh) | 基于标靶装置的深埋隧洞结构面点云坐标变换方法及设备 | |
US10690499B2 (en) | Analysis system, analysis method, and storage medium in which analysis program is stored | |
WO2018099051A1 (zh) | 一种rtk倾斜补偿测量精度验证的方法 | |
CN107917693B (zh) | 一种基于光学测距的倾斜测量装置及测量方法 | |
CN106767926A (zh) | 一种投线仪的数字化校准系统及方法 | |
CN107167109B (zh) | 两点激光测量平面空间偏角的方法及装置 | |
CN105180827A (zh) | 一种检测起重机主梁上拱曲线的装置 | |
CN103822580B (zh) | 超长框架变形与姿态的多点实时测量系统与方法 | |
CN105466455B (zh) | 一种测角精度标定系统及方法 | |
CN105043280B (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190423 |
|
RJ01 | Rejection of invention patent application after publication |