CN103679774A - 一种多边形农田作业区域边界建模方法 - Google Patents
一种多边形农田作业区域边界建模方法 Download PDFInfo
- Publication number
- CN103679774A CN103679774A CN201410002900.XA CN201410002900A CN103679774A CN 103679774 A CN103679774 A CN 103679774A CN 201410002900 A CN201410002900 A CN 201410002900A CN 103679774 A CN103679774 A CN 103679774A
- Authority
- CN
- China
- Prior art keywords
- summit
- farmland operation
- summits
- abscissa value
- last
- 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.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种多边形农田作业区域边界建模方法,采集多边形农田作业区域边界的n个顶点,将该n个顶点放到直角坐标系的第一象限,确定横坐标值最大的顶点和横坐标值最小的顶点,建立上边界函数模型和下边界函数模型,完成边界建模。本发明的方法可以将直线型多边形农田作业区域的边界以及农田作业区域内的障碍物顶点精准地记录并绘制出来,从而为农业植保机的作业人员提供准确的农田作业区域,为农业植保机的路径和航迹规划以及计算农田作业面积提供可靠的依据。
Description
技术领域
本发明涉及一种多边形农田作业区域边界建模方法。
背景技术
随着科技的不断发展,采用现代化机械来替换人工劳动已成为各行各业的流行趋势。我国作为一个传统的农业大国,耕地面积十分广阔,然而,目前国内农田作业方面还是采用传统的人工作业方式。其中农药喷洒采取的是人工喷洒的方式,这种方式不仅效率低,而且对作业人员的身心有极大的伤害,因此,急需先进的技术来改变这一现象。农业植保机应运而生,然而,现今所有的植保机都是靠人工进行操作,对农田的喷洒区域及边界完全靠操作员的眼睛进行判断,难免产生误判、漏喷及多喷的现象。因此,如果操作员能够获得准确的农田作业区域的边界,在无人直升机进行作业时,实时的将路径和航迹在显示器上显示出来,操作员便可根据航迹判断出飞机的位置,从而很大程度的解决误判、漏喷及多喷的问题,且大大提高了农业喷洒的安全性。而这一问题的关键便在于要准确的获取农田的边界信息。因此,急需一种好的农田作业区域边界的建模方法来解决这一实际问题。
北京农业信息技术研究中心的发明专利“采集农田关键顶点测绘成图的方法”提供了一种采集农田关键顶点测绘成图的方法,通过采集四种交界顶点、两种轮廓顶点一级一种辅助顶点后,根据关键顶点间的拓扑关系实时自动校验并查看测绘结果。包括以下步骤:S1:获取GPS位置信息;S2:勾勒待测区域轮廓;S3:测绘农田关键分界顶点,并标注地块名称和注释;S4:实时校验、提示与分割农田;S5:上传测绘数据;S6:获取矢量地图。
其中,在S2中,勾勒待测区域轮廓包括:顺序标定和无序标定,所述顺序标定,以轮廓顶点的编号为顺序,依次标定待测区域轮廓;所述无序标定,是指自动标定包含全部轮廓顶点的面积最大的多边形为区域轮廓。
现有方法是以绘出农田的矢量地图为目的,而并非以画出农田作业区域边界为目的,因此,所绘出的农田轮廓包含两个部分:作业区域和非作业区域。在将此农田轮廓图用于农业植保机喷洒农药的路径和航迹规划等场合时,会出现非作业区域的多喷现象,浪费农药;在用于计算农田作业面积时,会出现作业面积的计算不准确现象;现有方法未对农田区域边界建立数学模型,从而无法为农田作业区域的路径和航迹规划提供可靠依据;没有将农田作业区域以内可能存在的障碍物顶点(如电线杆、树、信号发射塔等)标记出来,在将此农田轮廓图用于农业植保机喷洒农药的路径和航迹规划等场合时,可能导致农业植保机出现撞机或坠毁危险现象发生。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供一种多边形农田作业区域边界建模方法,将直线型多边形农田作业区域的边界以及农田作业区域内的障碍物顶点精准地记录并绘制出来,从而为农业植保机的作业人员提供准确的农田作业区域,为农业植保机的路径和航迹规划以及计算农田作业面积提供可靠的依据。
为解决上述技术问题,本发明所采用的技术方案是:一种多边形农田作业区域边界建模方法,该方法为:
1)从多边形农田作业区域边界的某一个顶点D1开始,以顺时针方向采集多边形农田作业区域边界的所有顶点,设共有n个顶点,分别为D1(x1,y1),D2(x2,y2),D3(x3,y3),...,Dn(xn,yn);建立直角坐标系,将所述n个顶点设置在直角坐标系的第一象限内;
2)比较上述n个顶点的横坐标,找出横坐标值最大的顶点和最小的顶点;若有多个顶点的横坐标值均为最小值,且该多个顶点是第一个顶点与最后一个顶点之间的顶点,则将次序最先的顶点作为横坐标值最小的顶点,若第一个顶点和最后一个顶点的横坐标值都是最小值,则将最后一个顶点作为横坐标值最小的顶点;若有多个顶点的横坐标值均为最大值,且该多个顶点是第一个顶点与最后一个顶点之间的顶点,则将次序最后的顶点作为横坐标值最大的顶点,若第一个顶点和最后一个顶点的横坐标值都是最大值,则将第一个顶点作为横坐标值最大的顶点;将横坐标值最小的顶点Di(xi,yi)和最大的顶点Dj(xj,yj)分别作为所有的顶点在直角坐标系中的左端和右端;
3)以Di(xi,yi)为起始点,Dj(xj,yj)为终点;若j>i,按顺时针方向依次求得横坐标位于区间[xi,xj]中的次序相邻的两顶点的直线方程;若j<i,则将第一个顶点D1作为最后一个顶点Dn的下一点,按顺时针方向依次求得横坐标位于区间[xi,xj]中的次序相邻的两顶点的直线方程;确定多边形农田作业区域上边界函数模型S(x):
4)以Di(xi,yi)为起始顶点,Dj(xj,yj)为终顶点,若j<i,则按逆时针方向依次求得横坐标位于区间[xi,xj]中的次序相邻的两顶点的直线方程;若j>i,则将最后一个顶点Dn作为第一个顶点D1的下一点,按逆时针方向依次求得横坐标位于区间[xi,xj]中的次序相邻的两顶点的直线方程;确定多边形农田作业区域下边界函数模型W(x)。
确定了上边界函数模型S(x)和下边界函数模型W(x),即确定了整个多边形农田作业区域的模型。
当所述多边形农田作业区域内存在障碍物时,采集多边形农田作业区域内障碍物顶点的坐标Z(x0,y0),建立以下障碍物警戒线函数方程:
一级障碍物警戒线的函数方程为(x-x0)2+(y-y0)2=a2;
二级障碍物警戒线的函数方程为(x-x0)2+(y-y0)2=b2;
三级障碍物警戒线的函数方程为(x-x0)2+(y-y0)2=c2;
其中,0<a<b<c<50。
与现有技术相比,本发明所具有的有益效果为:本发明的方法可以将直线型多边形农田作业区域的边界以及农田作业区域内的障碍物顶点精准地记录并绘制出来,从而为农业植保机的作业人员提供准确的农田作业区域,为农业植保机的路径和航迹规划以及计算农田作业面积提供可靠的依据。
附图说明
图1为本发明一实施例直线型边界模型示意图;
图2为本发明一实施例有障碍物的农田作业区域障碍物警戒线模型。
具体实施方式
以下结合附图详细说明本发明的具体实现过程。
1)假设采点工作者以第一个点开始,用GPS定位的方式以顺时针方向采集了多边形农田作业区域的n个顶点,分别为D1(x1,y1),D2(x2,y2),D3(x3,y3),...,Dn(xn,yn),如图1所示(以n=4为例)。
2)比较这n个顶点的横坐标,找出其中横坐标最小和最大的点。假设横坐标最小的点为Di(xi,yi)(1≤i≤n),横坐标最大的点为Dj(xj,yj)(1≤j≤n)(若遇到多个点的横坐标都为最小值,则把采点次序最先的点作为最小点,在此有个特例:如果D1和Dn的横坐标都相等,但是在此取Dn作为最小点;若遇到多个点的横坐标都为最大值,则把取点次序最后的点作为最大点,在此有个特例:如果Dn和D1的横坐标相等,但是在此取D1作为最大点),那么Di(xi,yi)即为左端,Dj(xj,yj)即为右端。
3)以左端Di为起始点,右端Dj为终点,按顺时针方向依次求得横坐标位于区间[xi,xj]中的采点次序相邻的两点的直线方程。(若j<i,则在顺时针方向移点的过程中,Dn出现在Dj的左边(即xn<xj),那么在计算由Dn和Dn的下一点所确定的直线方程时,Dn的下一点即为D1,计算出Dn和D1两点所确定的直线方程,若D1为此时的右端,则移点结束,否则继续移点计算直线方程直到计算出由Dj的上一个点和Dj所确定的直线方程为止。)那么这些由采点次序相邻的两点所确定的所有直线就构成了上边界。
4)以横坐标Di为起始点,Dj为终点,按逆时针方向依次求得横坐标位于区间[xi,xj]中的采点次序相邻的两点的直线方程。(若j>i,则在逆时针方向移点的过程中,D1出现在Dj的左边(即x1<xj),那么在计算由D1和D1的下一点所确定的直线方程时,D1的下一点即为Dn,计算出由D1和Dn两点所确定的直线方程,若Dn为此时的右端,则移点结束,否则继续移点计算直线方程直到计算出由Dj的上一个点和Dj所确定的直线方程为止。)那么这些由采点次序相邻的两点所确定的所有直线就构成了下边界。
如果农田中有障碍物,则由采点工作者采集到障碍物点的坐标为Z(x0,y0),假设一级警戒线为离障碍物点a米的圆,二级警戒线为离障碍物点b米的圆,三级警戒线为离障碍物点c米的圆(50>c>b>a>0),如图2所示,图2中三个圆从内至外分别为一级警戒线、二级警戒线、三级警戒线。
则一级警戒线的函数方程为(x-x0)2+(y-y0)2=a2;
二级警戒线的函数方程为(x-x0)2+(y-y0)2=b2;
三级警戒线的函数方程为(x-x0)2+(y-y0)2=c2。
Claims (2)
1.一种多边形农田作业区域边界建模方法,其特征在于,该方法为:
1)从多边形农田作业区域边界的某一个顶点D1开始,以顺时针方向采集多边形农田作业区域边界的所有顶点,设共有n个顶点,分别为D1(x1,y1),D2(x2,y2),D3(x3,y3),...,Dn(xn,yn);建立直角坐标系,将所述n个顶点设置在直角坐标系的第一象限内;
2)比较上述n个顶点的横坐标,找出横坐标值最大的顶点和最小的顶点;若有多个顶点的横坐标值均为最小值,且该多个顶点是第一个顶点与最后一个顶点之间的顶点,则将次序最先的顶点作为横坐标值最小的顶点,若第一个顶点和最后一个顶点的横坐标值都是最小值,则将最后一个顶点作为横坐标值最小的顶点;若有多个顶点的横坐标值均为最大值,且该多个顶点是第一个顶点与最后一个顶点之间的顶点,则将次序最后的顶点作为横坐标值最大的顶点,若第一个顶点和最后一个顶点的横坐标值都是最大值,则将第一个顶点作为横坐标值最大的顶点;将横坐标值最小的顶点Di和最大的顶点Dj分别作为所有的顶点在直角坐标系中的左端和右端;
3)以Di为起始点,Dj为终点;若j>i,按顺时针方向依次求得横坐标位于区间[xi,xj]中的次序相邻的两顶点的直线方程;若j<i,则将第一个顶点D1作为最后一个顶点Dn的下一点,按顺时针方向依次求得横坐标位于区间[xi,xj]中的次序相邻的两顶点的直线方程;确定多边形农田作业区域上边界函数模型S(x):
4)以Di为起始顶点,Dj为终顶点,若j<i,则按逆时针方向依次求得横坐标位于区间[xi,xj]中的次序相邻的两顶点的直线方程;若j>i,则将最后一个顶点Dn作为第一个顶点D1的下一点,按逆时针方向依次求得横坐标位于区间[xi,xj]中的次序相邻的两顶点的直线方程;确定多边形农田作业区域下边界函数模型W(x)。
2.根据权利要求1所述多边形农田作业区域边界建模方法,其特征在于,当所述多边形农田作业区域内存在障碍物时,采集多边形农田作业区域内障碍物顶点的坐标Z(x0,y0),建立以下障碍物警戒线函数方程:
一级障碍物警戒线的函数方程为(x-x0)2+(y-y0)2=a2;
二级障碍物警戒线的函数方程为(x-x0)2+(y-y0)2=b2;
三级障碍物警戒线的函数方程为(x-x0)2+(y-y0)2=c2;
其中,0<a<b<c<50。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410002900.XA CN103679774B (zh) | 2014-01-03 | 2014-01-03 | 一种多边形农田作业区域边界建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410002900.XA CN103679774B (zh) | 2014-01-03 | 2014-01-03 | 一种多边形农田作业区域边界建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103679774A true CN103679774A (zh) | 2014-03-26 |
CN103679774B CN103679774B (zh) | 2016-06-22 |
Family
ID=50317221
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410002900.XA Active CN103679774B (zh) | 2014-01-03 | 2014-01-03 | 一种多边形农田作业区域边界建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103679774B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155993A (zh) * | 2014-08-11 | 2014-11-19 | 江苏恒创软件有限公司 | 一种基于无人直升机的坡地区域农作物作业方法 |
CN104503464A (zh) * | 2014-12-30 | 2015-04-08 | 中南大学 | 基于计算机的凸多边形农田无人机喷洒作业航迹规划方法 |
CN106599565A (zh) * | 2016-12-07 | 2017-04-26 | 中国电子科技集团公司第五十四研究所 | 基于位置时‑距分层的通航飞行器越界高效监视方法 |
CN106643719A (zh) * | 2016-09-23 | 2017-05-10 | 江西洪都航空工业集团有限责任公司 | 一种智能割草车的路径规划算法 |
WO2018006216A1 (en) * | 2016-07-04 | 2018-01-11 | SZ DJI Technology Co., Ltd. | Aerial operation support and real-time management |
CN108007437A (zh) * | 2017-11-27 | 2018-05-08 | 北京航空航天大学 | 一种基于多旋翼飞行器测量农田边界与内部障碍的方法 |
CN108257403A (zh) * | 2017-12-28 | 2018-07-06 | 北京千方科技股份有限公司 | 禁行区域标识方法及装置、识别方法及装置、智能车辆 |
CN109118502A (zh) * | 2018-10-18 | 2019-01-01 | 信阳师范学院 | 基于断点分割的作业覆盖区域实时重构方法及系统 |
CN109215096A (zh) * | 2018-10-30 | 2019-01-15 | 中国农业大学 | 一种快速的播种机已播区生成和判断方法 |
CN109633542A (zh) * | 2018-11-30 | 2019-04-16 | 中交遥感农业科技江苏有限公司 | 一种用于无人机的局域性地图网络的布置方法 |
CN109949198A (zh) * | 2019-02-22 | 2019-06-28 | 中国农业机械化科学研究院 | 一种麦田边界检测装置和检测方法 |
CN110598694A (zh) * | 2019-08-14 | 2019-12-20 | 广东奥普特科技股份有限公司 | 一种快速处理感兴趣区域图像的算法 |
CN110647141A (zh) * | 2018-06-27 | 2020-01-03 | 西安合众思壮导航技术有限公司 | 避障路径的生成方法、装置和系统 |
WO2023231286A1 (zh) * | 2022-06-02 | 2023-12-07 | 丰疆智能(深圳)有限公司 | 地块面积获取方法、系统、终端及计算机可读存储介质 |
WO2024125041A1 (zh) * | 2022-12-16 | 2024-06-20 | 美智光电科技股份有限公司 | 感应范围标定方法、装置、吊顶电器、设备和存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130100130A1 (en) * | 2011-10-21 | 2013-04-25 | IntegrityWare, Inc. | Methods and Systems for Generating and Editing Surfaces |
CN103247062A (zh) * | 2013-04-28 | 2013-08-14 | 北京农业信息技术研究中心 | 采集农田关键点测绘成图的方法 |
CN103279986A (zh) * | 2013-06-17 | 2013-09-04 | 陈超东 | 三维水平地质剖面图制作方法及其用途 |
-
2014
- 2014-01-03 CN CN201410002900.XA patent/CN103679774B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130100130A1 (en) * | 2011-10-21 | 2013-04-25 | IntegrityWare, Inc. | Methods and Systems for Generating and Editing Surfaces |
CN103247062A (zh) * | 2013-04-28 | 2013-08-14 | 北京农业信息技术研究中心 | 采集农田关键点测绘成图的方法 |
CN103279986A (zh) * | 2013-06-17 | 2013-09-04 | 陈超东 | 三维水平地质剖面图制作方法及其用途 |
Non-Patent Citations (1)
Title |
---|
于国英 等: "行播作物农田图像边界提取研究", 《安徽农业科学》 * |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155993B (zh) * | 2014-08-11 | 2017-07-25 | 江苏恒创软件有限公司 | 一种基于无人直升机的坡地区域农作物作业方法 |
CN104155993A (zh) * | 2014-08-11 | 2014-11-19 | 江苏恒创软件有限公司 | 一种基于无人直升机的坡地区域农作物作业方法 |
CN104503464A (zh) * | 2014-12-30 | 2015-04-08 | 中南大学 | 基于计算机的凸多边形农田无人机喷洒作业航迹规划方法 |
US11016488B2 (en) | 2016-07-04 | 2021-05-25 | SZ DJI Technology Co., Ltd. | Aerial operation support and real-time management |
WO2018006216A1 (en) * | 2016-07-04 | 2018-01-11 | SZ DJI Technology Co., Ltd. | Aerial operation support and real-time management |
WO2018006454A1 (zh) * | 2016-07-04 | 2018-01-11 | 深圳市大疆创新科技有限公司 | 一种无人机的飞行路径规划、控制方法及系统 |
US11703865B2 (en) | 2016-07-04 | 2023-07-18 | SZ DJI Technology Co., Ltd. | Aerial operation support and real-time management |
CN106643719A (zh) * | 2016-09-23 | 2017-05-10 | 江西洪都航空工业集团有限责任公司 | 一种智能割草车的路径规划算法 |
CN106643719B (zh) * | 2016-09-23 | 2020-04-14 | 江西洪都航空工业集团有限责任公司 | 一种智能割草车的路径规划算法 |
CN106599565B (zh) * | 2016-12-07 | 2018-11-06 | 中国电子科技集团公司第五十四研究所 | 基于位置时-距分层的通航飞行器越界高效监视方法 |
CN106599565A (zh) * | 2016-12-07 | 2017-04-26 | 中国电子科技集团公司第五十四研究所 | 基于位置时‑距分层的通航飞行器越界高效监视方法 |
CN108007437A (zh) * | 2017-11-27 | 2018-05-08 | 北京航空航天大学 | 一种基于多旋翼飞行器测量农田边界与内部障碍的方法 |
CN108257403A (zh) * | 2017-12-28 | 2018-07-06 | 北京千方科技股份有限公司 | 禁行区域标识方法及装置、识别方法及装置、智能车辆 |
CN110647141A (zh) * | 2018-06-27 | 2020-01-03 | 西安合众思壮导航技术有限公司 | 避障路径的生成方法、装置和系统 |
CN109118502B (zh) * | 2018-10-18 | 2021-11-02 | 信阳师范学院 | 基于断点分割的作业覆盖区域实时重构方法及系统 |
CN109118502A (zh) * | 2018-10-18 | 2019-01-01 | 信阳师范学院 | 基于断点分割的作业覆盖区域实时重构方法及系统 |
CN109215096A (zh) * | 2018-10-30 | 2019-01-15 | 中国农业大学 | 一种快速的播种机已播区生成和判断方法 |
CN109215096B (zh) * | 2018-10-30 | 2022-09-30 | 中国农业大学 | 一种快速的播种机已播区生成和判断方法 |
CN109633542A (zh) * | 2018-11-30 | 2019-04-16 | 中交遥感农业科技江苏有限公司 | 一种用于无人机的局域性地图网络的布置方法 |
CN109949198A (zh) * | 2019-02-22 | 2019-06-28 | 中国农业机械化科学研究院 | 一种麦田边界检测装置和检测方法 |
CN109949198B (zh) * | 2019-02-22 | 2020-12-22 | 中国农业机械化科学研究院 | 一种麦田边界检测装置和检测方法 |
CN110598694A (zh) * | 2019-08-14 | 2019-12-20 | 广东奥普特科技股份有限公司 | 一种快速处理感兴趣区域图像的算法 |
WO2023231286A1 (zh) * | 2022-06-02 | 2023-12-07 | 丰疆智能(深圳)有限公司 | 地块面积获取方法、系统、终端及计算机可读存储介质 |
WO2024125041A1 (zh) * | 2022-12-16 | 2024-06-20 | 美智光电科技股份有限公司 | 感应范围标定方法、装置、吊顶电器、设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103679774B (zh) | 2016-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103679774A (zh) | 一种多边形农田作业区域边界建模方法 | |
CN107368094A (zh) | 一种无人机植保作业航线规划方法及装置 | |
CN105629973B (zh) | 一种基于rtk技术的农机直线行驶的方法及装置 | |
CN106643719B (zh) | 一种智能割草车的路径规划算法 | |
CN103499973B (zh) | 一种主-从机协同作业农业机械智能导航系统 | |
CN103950540B (zh) | 一种基于无线传感器网络的植保无人机喷施作业方法 | |
CN107272726A (zh) | 基于无人机植保作业的作业区域确定方法及装置 | |
CN109240284A (zh) | 一种无人驾驶农机的自主路径规划方法及装置 | |
CN112099556A (zh) | 农业无人飞行器的控制方法、地面控制端及存储介质 | |
CN107289950A (zh) | 植保无人机作业航线规划方法及植保无人机 | |
CN103988824A (zh) | 一种基于双目视觉技术的自动对靶喷施系统 | |
CN105557672A (zh) | 一种果树靶标探测系统 | |
CN103679775A (zh) | 一种直线和曲线相结合的农田作业区域边界建模方法 | |
CN102800083A (zh) | 基于双目视觉网格划分匹配算法的农作物喷雾定位方法 | |
CN105197243A (zh) | 一种农用无人机机载变量施药系统及方法 | |
CN206524983U (zh) | 一种基于uwb的仓库叉车实时高精度定位系统 | |
CN112154447A (zh) | 地表特征识别方法、设备、无人机及计算机可读存储介质 | |
CN206057975U (zh) | 一种避障装置及植保无人机 | |
CN106547276A (zh) | 自动喷洒回字形路径规划方法及植保机喷洒作业方法 | |
CN110619649A (zh) | 作业面积确定方法及装置、终端 | |
CN205180165U (zh) | 北斗导航定位农药喷洒控制系统 | |
CN204104581U (zh) | 一种基于双目视觉技术的自动对靶喷施系统 | |
CN203860304U (zh) | 一种自动对靶喷施系统 | |
CN105045280A (zh) | 一种基于卫星系统的自走式喷施机导航系统以及方法 | |
CN106774434A (zh) | 应用于无人机植保的自动避障方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |