CN104808679B - General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction - Google Patents

General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction Download PDF

Info

Publication number
CN104808679B
CN104808679B CN201510089009.9A CN201510089009A CN104808679B CN 104808679 B CN104808679 B CN 104808679B CN 201510089009 A CN201510089009 A CN 201510089009A CN 104808679 B CN104808679 B CN 104808679B
Authority
CN
China
Prior art keywords
file
direction angle
interest
target direction
information
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.)
Expired - Fee Related
Application number
CN201510089009.9A
Other languages
Chinese (zh)
Other versions
CN104808679A (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.)
Civil Aviation University of China
Original Assignee
Civil Aviation University of China
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 Civil Aviation University of China filed Critical Civil Aviation University of China
Priority to CN201510089009.9A priority Critical patent/CN104808679B/en
Publication of CN104808679A publication Critical patent/CN104808679A/en
Application granted granted Critical
Publication of CN104808679B publication Critical patent/CN104808679B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

A kind of General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction.Its mapping relations and generation XML configuration file for including setting up navigational information assembling file and geography information;Record flight path information;The movement tendency of aircraft is predicted based on polynomial prediction method;Screening POI points of interest;Reversely search the steps such as navigational information assembling file.The General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction that the present invention is provided have advantages below:1. the mapping relations of navigational information assembling file and geography information are established, the brand-new file hierarchy being easily managed is formd.2. the Forecasting Methodology of flight path is employed, the intelligent Matching of file is realized, solves the problems, such as that traditional ff is difficult.3. computing is simple, highly reliable.4. more convenient, intelligent, three-dimensional navigational intelligence service can be provided.

Description

General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction
Technical field
The invention belongs to blank pipe technical field of automation, more particularly to a kind of General Aviation based on Flight Trajectory Prediction AIP (Aeronautial Information Publication, navigational information compilation) file intelligent Matching method.
Background technology
General Aviation refers to that the Civil Aviation Activities beyond public air transport, events in operation are engaged in using civil aircraft Cover agriculture, woods, herd, fishing, industry, building, scientific research, traffic, multiple industries such as amusement.The particular content of General Aviation covers ten thousand As, with maneuverability, rapidly and efficiently the features such as, in the last few years, with the development of civil aviaton of China cause and opening for low altitude airspace Put, the industry size of General Aviation expands day by day, flight demand is gradually vigorous, therefore with very big market prospects.
Navigational intelligence service is the important component in operational service.Accurately and timely navigational intelligence is provided to unit Service is the basis of safe flight.General Aviation operation has very strong flexibility, and airborne equipment is relatively easy, therefore, logical It is very high with being required navigational intelligence service convenience, ease for use in aviation.Navigational intelligence service include to pilot provide chart, Navigational information collects and informs the multinomial content such as notice to navigator.Chart and navigational information assembling file refer to need in flight course Using files such as regulations handbook, airport and the en-route datas for arriving, quantity is big, and content is numerous and diverse.Traditional navigational information ways of distribution It is to provide a large amount of papery version files to unit, uses extremely inconvenience, therefore be unsuitable for General Aviation.
At present, there are many manufacturers have developed both at home and abroad and used the EFB in intelligent movable equipment (Electronic Flight Bag system, EFB), the EFB is by chart and navigational information assembling file electronics Change, improve the convenience for using.But Railway Project is there is also simultaneously, and such as these manufacturers are generally planemaker, such as Boeing, sky Visitor etc., with very strong monopoly, therefore, equipment price is high, and maintenance cost is huge.Due to General Aviation operating cost compared with It is low, thus do not possess the condition for being equipped with EFB extensively.In addition, EFB is not carried out further to navigational information file Treatment, still needing to be browsed by file from level to level can just find the file for in-flight needing.
In view of the foregoing, General Aviation is badly in need of a kind of intellectual search method of electronics navigational information assembling file to solve Certainly complicated ground ff problem in traditional approach and existing EFB, so as to ensure the flight safety of General Aviation.
The content of the invention
In order to solve the above problems, it is easily managed it is an object of the invention to provide one kind, realizes intelligent of file Match somebody with somebody, broken away from traditional heap file and searched difficult problem, using the teaching of the invention it is possible to provide the navigational intelligence service of intelligent three-dimensional based on flying The General Aviation AIP file intelligent Matching methods of row trajectory predictions.
In order to achieve the above object, intelligent of the General Aviation AIP files based on Flight Trajectory Prediction that the present invention is provided Method of completing the square includes the following steps for carrying out in order:
Step (1) sets up navigational information assembling file with the mapping relations of geography information and the S01 of generation XML configuration file Stage:It is according to letter including the catalogue including course line, guidance station, Aerial Ports File in initial navigational information assembling file It is tactic, these files, and self-defined two XML configuration files are arranged successively to describe navigational information assembling file with ground Mapping relations between reason information;
Step (2) records the S02 stages of flight path information:In flight course, current newest 5 of aircraft real time record Track points information, including speed, longitude and latitude, deflection, and store that data in local database;
The S03 stages that step (3) is predicted based on polynomial prediction method to the movement tendency of aircraft:If ti(i= 1,2,3 ..., n) moment target direction angle s (ti) can be represented with one group of orthogonal polynomial function of time, further according to ti(i=1,2, 3 ..., n) the measured value x at moment target direction angleiUsing least square method to coefficient ajLeast-squares estimation is carried out to obtainMost Afterwards according to coefficientTo (tn+ dt) moment target direction angle is predicted;
Step (4) screens the S04 stages of POI points of interest:Target direction angle according to prediction defines a deflection nargin θ, the geography information in traversal XML configuration file, using the mapping relations between geography information and navigational information assembling file, will These points with navigational information assembling file attribute are defined as POI points of interest, afterwards the latitude and longitude coordinates meter according to aircraft The declinate α between POI points of interest and the target direction angle of prediction is calculated, and whether judges declinate α less than deflection nargin θ, if Declinate α then continues to calculate the distance between the POI points of interest and aircraft within deflection nargin θ, and the minimum point of distance is just It is matching degree highest point;
Step (5) reversely searches the S05 stages of navigational information assembling file:Navigated by water by parsing XML configuration file Mapping relations between data compilation file and geography information, i.e., can be reverse according to the geography information of the POI points of interest of matching Corresponding navigational information assembling file is found, so as to realize the intelligent Matching of navigational information assembling file.
In step (3), the described method being predicted to the movement tendency of aircraft based on polynomial prediction method Including the following steps for performing in order:
1) t is seti(i=1,2,3 ..., n) moment target direction angle s (ti) an orthogonal polynomial function of time table can be used Show, wherein ti=iT, T are the sampling interval of track:
2) according to ti(i=1,2,3 ..., n) the measured value x at moment target direction angleiUsing least square method to coefficient aj Least-squares estimation is carried out, can be solved by following formula:
A=(PT·P)-1PTX (2)
In formula:
Therefore the optimal estimation of target direction angle s (t) is:
Thus can be using formula (3) to (tn+ dt) the target direction angle at moment is predicted, i.e.,:
In above-mentioned formula, n=5, n represents historical track sampled point, xiIt is 5 prescriptions of database output in the S02 stages To angular data, the deflection of t is exactlyThen the polynomial function at target direction angle is:
Recycle 5 newest prescriptions to angular data, coefficient a is determined using formula (2)iLeast-squares estimationSo Just being capable of basisTo tnThe movement locus at moment is predicted, i.e.,:
In step (4), the method for described screening POI points of interest includes the following steps for performing in order:
1) the target direction angle according to prediction defines a deflection nargin θ, and direction angle nargin θ must be at 0 ° -180 ° Between, according to the motor-driven parameter of navigation aircraft, deflection nargin θ is typically taken as 60 °;
2) XML configuration file is parsed, navigational information assembling file and geography information is obtained, is had according to geographical information definition The POI points of interest of navigational information assembling file attribute;
3) latitude and longitude coordinates (x according to aircraft0,y0), the latitude and longitude coordinates (x of POI points of interest1,y1), and prediction Target direction angle φ calculate declinate α between POI points of interest and the target direction angle of prediction, specific formula is as follows:
4) whether declinate α is judged less than deflection nargin θ, if declinate calculates the POI within deflection nargin θ Then these POI points of interest are ranked up by the distance between point of interest and aircraft by distance, and the minimum point of distance is exactly With degree highest point.
The General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction that the present invention is provided have following excellent Point:1. the mapping relations of navigational information assembling file and geography information are established, the brand-new file layers being easily managed are formd It is secondary.2. the Forecasting Methodology of flight path is employed, the intelligent Matching of file is realized, traditional ff difficulty is solved Problem.3. computing is simple, highly reliable.4. more convenient, intelligent, three-dimensional navigational intelligence service can be provided.
Brief description of the drawings
The General Aviation AIP file intelligent Matching method flow diagrams based on Flight Trajectory Prediction that Fig. 1 is provided for the present invention.
Fig. 2 is the method flow being predicted to the movement tendency of aircraft based on polynomial prediction method in the present invention Figure.
Fig. 3 is the method flow diagram of screening POI points of interest in the present invention.
Fig. 4 is the schematic diagram that target direction angle in the present invention according to prediction defines deflection nargin θ.
Specific embodiment
The General Aviation AIP files based on Flight Trajectory Prediction for being provided the present invention with instantiation below in conjunction with the accompanying drawings Intelligent Matching method is described in detail.
As shown in figure 1, the General Aviation AIP file intelligent Matching method bags based on Flight Trajectory Prediction that the present invention is provided Include the following steps for performing in order:
Step (1) sets up navigational information assembling file with the mapping relations of geography information and the S01 of generation XML configuration file Stage:It is according to letter including the catalogue including course line, guidance station, Aerial Ports File in initial navigational information assembling file It is tactic, arrange these files successively, and self-defined two XML (Extensible Markup Language, it is expansible Markup language) configuration file describes the mapping relations between navigational information assembling file and geography information;
Step (2) records the S02 stages of flight path information:In flight course, current newest 5 of aircraft real time record Track points information, including speed, longitude and latitude, deflection etc., and stores that data in local database;
The S03 stages that step (3) is predicted based on polynomial prediction method to the movement tendency of aircraft:If ti(i= 1,2,3 ..., n) moment target direction angle s (ti) can be represented with one group of orthogonal polynomial function of time, further according to ti(i=1,2, 3 ..., n) the measured value x at moment target direction angleiUsing least square method to coefficient ajLeast-squares estimation is carried out to obtainMost Basis afterwardsTo (tn+ dt) moment target direction angle is predicted;
As shown in Fig. 2 the described method bag being predicted to the movement tendency of aircraft based on polynomial prediction method Include the following steps for performing in order:
1) t is seti(i=1,2,3 ..., n) moment target direction angle s (ti) an orthogonal polynomial function of time table can be used Show, wherein ti=iT, T are the sampling interval of track:
2) according to ti(i=1,2,3 ..., n) the measured value x at moment target direction angleiUsing least square method to coefficient aj Least-squares estimation is carried out, can be solved by following formula:
A=(PT·P)-1PTX (2)
In formula:
Therefore the optimal estimation of target direction angle s (t) is:
Thus can be using formula (3) to (tn+ dt) the target direction angle at moment is predicted, i.e.,:
In above-mentioned formula, n=5 (n represents historical track sampled point), xiIt is 5 prescriptions of database output in the S02 stages To angular data, the deflection of t is exactlyThen the polynomial function at target direction angle is:
Recycle 5 newest prescriptions to angular data, coefficient a is determined using formula (2)iLeast-squares estimationThus Being capable of basisTo tnThe movement locus at moment is predicted, i.e.,:
Step (4) screens the S04 stages of POI points of interest:Target direction angle according to prediction defines a deflection nargin θ, the geography information in traversal XML configuration file, using the mapping relations between geography information and navigational information assembling file, will These points with navigational information assembling file attribute are defined as POI points of interest, are calculated according to aircraft latitude and longitude coordinates afterwards Declinate α between POI points of interest and the target direction angle of prediction, and whether declinate α is judged less than deflection nargin θ, if poor Angle α then continues to calculate the distance between the POI points of interest and aircraft within deflection nargin θ, and the minimum point of distance is exactly Matching degree highest point;
As shown in figure 3, the method for described screening POI points of interest includes the following steps for performing in order:
1) as shown in figure 4, the target direction angle according to prediction defines a deflection nargin θ, nargin θ is necessary at direction angle Between 0 ° -180 °, and can not be excessive or too small, the effect of screening POI points of interest is otherwise just lost, according to navigation aircraft Motor-driven parameter, is typically taken as 60 ° by deflection nargin θ;
2) XML configuration file is parsed, navigational information assembling file and geography information is obtained, is had according to geographical information definition The POI points of interest of navigational information assembling file attribute;
3) latitude and longitude coordinates (x according to oneself aircraft0,y0), the latitude and longitude coordinates (x of POI points of interest1,y1), and The target direction angle φ of prediction calculates the declinate α between POI points of interest and the target direction angle of prediction, and specific formula is as follows:
4) whether declinate α is judged less than deflection nargin θ, if declinate calculates the POI within deflection nargin θ Then these POI points of interest are ranked up by the distance between point of interest and aircraft by distance, and the minimum point of distance is exactly With degree highest point.
Step (5) reversely searches the S05 stages of navigational information assembling file:Navigated by water by parsing XML configuration file Mapping relations between data compilation file and geography information, i.e., can be reverse according to the geography information of the POI points of interest of matching Corresponding navigational information assembling file is found, so as to realize the intelligent Matching of navigational information assembling file.

Claims (3)

1. a kind of General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction, it is characterised in that:It is included by suitable The following steps that sequence is performed:
Step (1) sets up navigational information assembling file with the mapping relations of geography information and the S01 ranks of generation XML configuration file Section:It is suitable according to letter including the catalogue including course line, guidance station, Aerial Ports File in initial navigational information assembling file Sequence arrangement, these files, and self-defined two XML configuration files are arranged successively to describe navigational information assembling file with geography Mapping relations between information;
Step (2) records the S02 stages of flight path information:In flight course, the current 5 newest flight paths of aircraft real time record Point information, including speed, longitude and latitude, deflection, and store that data in local database;
The S03 stages that step (3) is predicted based on polynomial prediction method to the movement tendency of aircraft:If ti, i=1,2, 3 ..., n, moment target direction angle s (ti) can be represented with one group of orthogonal polynomial function of time, further according to ti, i=1,2,3 ..., N, the measured value x at moment target direction angleiUsing least square method to coefficient ajLeast-squares estimation is carried out to obtainLast root According to coefficientTo (tn+ dt) moment target direction angle is predicted;
Step (4) screens the S04 stages of POI points of interest:Target direction angle according to prediction defines a deflection nargin θ, time The geography information gone through in XML configuration file, using the mapping relations between geography information and navigational information assembling file, by these Point with navigational information assembling file attribute is defined as POI points of interest, and the latitude and longitude coordinates according to aircraft calculate POI afterwards Declinate α between point of interest and the target direction angle of prediction, and whether declinate α is judged less than deflection nargin θ, if declinate α Then continue to calculate the distance between the POI points of interest and aircraft within deflection nargin θ, the minimum point of distance is exactly to match Degree highest point;
Step (5) reversely searches the S05 stages of navigational information assembling file:Navigational information is obtained by parsing XML configuration file Mapping relations between assembling file and geography information, i.e., can reversely search according to the geography information of the POI points of interest of matching To corresponding navigational information assembling file, so as to realize the intelligent Matching of navigational information assembling file.
2. the General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction according to claim 1, its feature It is:In step (3), described is included based on polynomial prediction method to the method that the movement tendency of aircraft is predicted The following steps for performing in order:
1) t is seti, i=1,2,3 ..., n, moment target direction angle s (ti) can be represented with an orthogonal polynomial function of time, its Middle ti=iT, T are the sampling interval of track:
2) according to ti, i=1,2,3 ..., n, the measured value x at moment target direction angleiUsing least square method to coefficient ajCarry out Least-squares estimation, can be solved by following formula:
A=(PT·P)-1PTX (2)
In formula:
Therefore target direction angle s (ti) optimal estimation be:
Thus can be using formula (3) to (tn+ dt) the target direction angle at moment is predicted, i.e.,:
In above-mentioned formula, n=5, n represents historical track sampled point, xiIt is 5 groups of deflection numbers of database output in the S02 stages According to the deflection of t is exactlyThen the polynomial function at target direction angle is:
Recycle 5 newest prescriptions to angular data, coefficient a is determined using formula (2)iLeast-squares estimationThis makes it possible to According toTo tnThe target direction angle at moment is predicted, i.e.,:
3. the General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction according to claim 1, its feature It is:In step (4), the method for described screening POI points of interest includes the following steps for performing in order:
1) the target direction angle according to prediction defines a deflection nargin θ, direction angle nargin θ must between 0 ° -180 °, According to the motor-driven parameter of navigation aircraft, deflection nargin θ is taken as 60 °;
2) XML configuration file is parsed, navigational information assembling file and geography information is obtained, navigation is had according to geographical information definition The POI points of interest of data compilation file attribute;
3) latitude and longitude coordinates (x according to aircraft0,y0), the latitude and longitude coordinates (x of POI points of interest1,y1), and the target predicted DeflectionThe declinate α between POI points of interest and the target direction angle of prediction is calculated, specific formula is as follows:
4) judge that whether declinate α, less than deflection nargin θ, if declinate is within deflection nargin θ, calculates the POI interest Then these POI points of interest are ranked up by the distance between point and aircraft by distance, and the minimum point of distance is exactly matching degree Highest point.
CN201510089009.9A 2015-02-27 2015-02-27 General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction Expired - Fee Related CN104808679B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510089009.9A CN104808679B (en) 2015-02-27 2015-02-27 General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510089009.9A CN104808679B (en) 2015-02-27 2015-02-27 General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction

Publications (2)

Publication Number Publication Date
CN104808679A CN104808679A (en) 2015-07-29
CN104808679B true CN104808679B (en) 2017-06-23

Family

ID=53693596

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510089009.9A Expired - Fee Related CN104808679B (en) 2015-02-27 2015-02-27 General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction

Country Status (1)

Country Link
CN (1) CN104808679B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106529705A (en) * 2016-10-25 2017-03-22 合肥飞友网络科技有限公司 Low-frequency latitude and longitude data supplement optimization flight locus method
CN106933246B (en) * 2017-03-29 2020-04-14 厦门大学 Multi-unmanned aerial vehicle complex task planning method
CN108663683B (en) * 2017-04-01 2020-07-14 湖北华中光电科技有限公司 Distance wave gate algorithm
CN112286049A (en) * 2019-07-27 2021-01-29 华为技术有限公司 Motion trajectory prediction method and device
CN110608737B (en) * 2019-09-16 2023-03-10 中国民航大学 Land-based radio navigation station selection method based on total system error estimation
CN112911497B (en) * 2019-12-02 2022-03-29 南京航空航天大学 Real-time prediction method and system for track of cooperative unmanned aerial vehicle

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6775659B2 (en) * 1998-08-26 2004-08-10 Symtec Limited Methods and devices for mapping data files
CN202221566U (en) * 2011-07-08 2012-05-16 中国民航科学技术研究院 Flight programming system and verification platform of performance-based navigation
CN102867073A (en) * 2011-07-08 2013-01-09 中国民航科学技术研究院 Flight program design system for performance-based navigation, verification platform and verification method
CN103973263A (en) * 2014-05-16 2014-08-06 中国科学院国家天文台 Novel approximation filter method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6775659B2 (en) * 1998-08-26 2004-08-10 Symtec Limited Methods and devices for mapping data files
CN202221566U (en) * 2011-07-08 2012-05-16 中国民航科学技术研究院 Flight programming system and verification platform of performance-based navigation
CN102867073A (en) * 2011-07-08 2013-01-09 中国民航科学技术研究院 Flight program design system for performance-based navigation, verification platform and verification method
CN103973263A (en) * 2014-05-16 2014-08-06 中国科学院国家天文台 Novel approximation filter method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
基于Android的通航航行资料电子包的设计与实现;吴恩森等;《计算机光盘软件与应用》;20141231(第24期);第93-95页 *
基于混杂系统和位置空间离散的冲突探测算法;石磊等;《系统工程与电子技术》;20140430;第36卷(第4期);第624-629页 *
通用航空综合运行支持流程设计;吴仁彪等;《中国民航大学学报》;20140630;第32卷(第3期);第1-5、19页 *
通用航空飞行计划航迹估算实现方法;吴仁彪等;《中国民航大学学报》;20140228;第32卷(第1期);第1-4、9页 *

Also Published As

Publication number Publication date
CN104808679A (en) 2015-07-29

Similar Documents

Publication Publication Date Title
CN104808679B (en) General Aviation AIP file intelligent Matching methods based on Flight Trajectory Prediction
Liu et al. Predicting aircraft trajectories: A deep generative convolutional recurrent neural networks approach
CN103557867B (en) The collaborative path planning method of a kind of many UAV of three-dimensional based on sparse A* search
CN103336863B (en) The flight intent recognition methods of flight path observed data of flying based on radar
CN101692315B (en) Method for analyzing high precision 4D flight trajectory of airplane based on real-time radar data
AU2011261838B2 (en) Determining emergency landing sites for aircraft
Pang et al. A recurrent neural network approach for aircraft trajectory prediction with weather features from sherlock
CN104615849B (en) A kind of flight plan assessment system and implementation method suitable for General Aviation
Patron et al. Flight trajectories optimization under the influence of winds using genetic algorithms
CN105701090A (en) Simulator for flight management system
CN104240541A (en) 4D track generating method
CN102692621A (en) ADS-B (automatic dependent surveillance broadcast) and radar combined system error estimation method
Eng Path planning, guidance and control for a uav forced landing
Holzäpfel et al. Assessment of wake-vortex encounter probabilities for crosswind departure scenarios
US10460608B2 (en) System and method for determining uncertainty in a predicted flight path for an aerial vehicle
Alligier et al. Energy rate prediction using an equivalent thrust setting profile
Ayhan et al. Preflight contingency planning approach for fixed wing UAVs with engine failure in the presence of winds
CN107909106A (en) A kind of detection method of aircraft flight environment
Schopferer et al. Trajectory risk modelling and planning for unmanned cargo aircraft
CN110135451A (en) A kind of track clustering method arriving line-segment sets distance based on point
Jiang et al. Research on method of trajectory prediction in aircraft flight based on aircraft performance and historical track data
Sibin et al. Research and modelling on performance database of flight management system
Qin et al. Research on dynamic route planning under adverse weather
Donkels Trajectory risk evaluation for autonomous low-flying air transport
Soares Flight Data Monitoring and its Application on Algorithms for Precursor Detection

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170623

Termination date: 20190227

CF01 Termination of patent right due to non-payment of annual fee