CN101877034A - Land utilization automatic division method - Google Patents
Land utilization automatic division method Download PDFInfo
- Publication number
- CN101877034A CN101877034A CN200910272814XA CN200910272814A CN101877034A CN 101877034 A CN101877034 A CN 101877034A CN 200910272814X A CN200910272814X A CN 200910272814XA CN 200910272814 A CN200910272814 A CN 200910272814A CN 101877034 A CN101877034 A CN 101877034A
- Authority
- CN
- China
- Prior art keywords
- soil
- particle
- land
- land utilization
- utilizes
- 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
Abstract
The invention relates to a land utilization automatic division method, and belongs to the field of land utilization planning. The automatic division method comprises the following steps: extracting and integrating the basic data of land utilization division; constructing a land utilization automatic division model based on a particle swarm optimization algorithm, establishing a mapping relationship between an unknown question and a particle swarm by taking the data as model input data and taking a land utilization patch as a data processing unit, finally solving a land utilization optimized layout result; and combining patches of adjacent same type land, and extracting the boundary of the combined patches so as to obtain a land utilization division boundary, and finally generating a land utilization division chart. The method has the advantages of well simulating the intelligent behavior of human being in the process of land utilization decision to improve the rationality and suitability of land utilization division so as to overcome the defect that the rationality of land-use zoning is poor caused by low efficiency of the conventional land utilization efficiency and hard simulation of the intelligent behavior of human being.
Description
Technical field
The present invention relates to a kind of soil and utilize the auto-partition method, refer to be applied to the soil especially and utilize the soil in the space optimization configuration to utilize the auto-partition method, belong to the land use planning field based on particle swarm optimization algorithm.
Background technology
It is the major technique means of land resource purposes control that the soil utilizes subregion, also is the key content that the soil utilizes the space optimization configuration.The soil utilizes subregion research to start from the regional soil that 20 th America start and utilizes the comprehensive sub-areas in the world, begun in mid-term in 20th century subsequently that the soil of land resource reasonable disposition and layout utilizes subregion research between all departments to solve, the classification subregion system of the territory planning of carrying out as Korea S, Japan.In the last few years, the soil utilized subregion research to begin to pay attention to the research of partitioning technique means and subregion theory, and subregion accuracy and practicality have also obtained certain raising.The soil of China utilizes subregion research beginning later, first representative achievement is the national present status of land utilization zoning of finishing the sixties in last century, the national potential of land resource zoning thought that the eighties, Chinese Academy of Sciences's geographical synthesis board of review proposed in the land resource map process in the establishment whole nation in 1: 100 ten thousand is divided into 9 REGION OF WATER INJECTION OILFIELD with national land resource; The soil utilizes the research of partitioning technique method to obtain paying close attention to more widely afterwards, proposed to utilize partition method than soils such as Y-factor method Y, cellular automatons based on improvement cluster analysis, space overlay analysis, planisphere method, mould, the proposition of these methods has partly overcome traditional soil and has utilized subregion based on qualitative, partition boundaries is fuzzy, subregion is random greatly, the subregion work efficiency is low, the subregion result utilizes the not strong defective of directive significance to the soil, and brought into play vital role in preceding two-wheeled overall plan for land use compilation process.Yet because the soil utilizes the complex nature of the problem, it is a land ecology economy composite system that relates to nature, society, economy that the soil utilizes system, the natural quality that the soil utilizes system can effectively be analyzed and simulate to simple measurement geography method and geographical computing system, the influence to the soil utilization at aspects such as the social economy of soil utilization and people's subjective desires then is difficult to obtain analytical effect preferably, particle swarm optimization algorithm also is the simulation that originates from simple social system, is the process that the simulation flock of birds is looked for food at first.
The tradition soil utilizes partition method, and geographical and geographical computing method are difficult to simulate the soil and utilize the complex behavior of people in the system soil to be utilized the influence of system based on metering, rationality is not high as a result to cause the soil to utilize subregion, and land resource is rationally utilized and the not strong defective of purposes control directiveness.
Summary of the invention
For addressing the above problem, the invention provides a kind of soil and utilize the auto-partition method, this method is that a kind of soil based on particle cluster algorithm utilizes the auto-partition method, utilize the newest research results of correlation techniques such as artificial intelligence and biological computation, come the anthropomorphic dummy to utilize intelligent behavior in the system in the soil, to improve rationality and the operability that the soil utilizes subregion.This method comprises the steps:
(1) the extraction soil utilizes the basic data of subregion and integrates, and the data after the integration are land use combination and each figure spot information wherein;
(2) make up the soil based on particle swarm optimization algorithm and utilize the auto-partition model, data with step (1) are model input data, with the land-use map spot is data processing unit, sets up the mapping relations wait to ask between problem and the population, finally finds the solution the soil utilization result that optimizes distribution.
(3) the soil utilization that step (2) is obtained is optimized distribution, and merges by the figure spot to contiguous identical land use pattern and handles, and extracts figure spot boundary line after merging and promptly obtains the soil and utilize zone boundary, and finally generate the soil and utilize block plan.
Described in the step (2) foundation to wait to ask the mapping relations between problem and the population as follows: the land-use map spot is abstract to be the particle in the particle swarm optimization algorithm, the position of each figure spot and land use pattern correspond to particle position (x, y) and kind (i), the soil utilizes the comprehensive benefit function to be the population fitness function.
Finding the solution described in the step (2) draws optimize distribution result's step of soil utilization and comprises:
(a) scale of population is set, the maximum flying speed (υ of particle
Xmax, υ
Ymax) and quicken weight (c
1, c
2), maximum iteration time, error amount, and each particle position of initialization, speed and attribute;
(b), calculate the fitness of each particle according to fitness function;
(c) calculate individual optimal value (p
b) and global optimum (p
g), and calculate the individual optimal location and the global optimum position of particle;
(d) calculate inertia weight, upgrade the flying speed of each particle, the individual optimal location and the global optimum position that draw according to step (c), the more current location of new particle;
(e) loop iteration: when satisfying Rule of judgment, search finishes, otherwise the value that step (d) is upgraded is as the input value of (b), continues to search for to obtain global optimum and global optimum position.
The current location of the described more new particle in the step (c) adopts the arest neighbors rule, and present bit is equipped with when a plurality of, and the arest neighbors rule adopts the area method that is dominant to determine and the most contiguous figure spot of current location.
Rule of judgment described in the step (d) is: the absolute value of the difference of twice global optimum in front and back is not more than error amount, and the number of times of loop iteration is not more than maximum iteration time.
The invention has the advantages that: particle swarm optimization algorithm can make full use of the computation capability of computing machine on the one hand, effectively improves soil subregion efficient; Simultaneously, utilize population intelligent optimization analog capability preferably simulating human utilize intelligent behavior in the decision process in the soil, utilize the rationality and the applicability of subregion to improve the soil, thereby solve the deficiency that current soil utilizes the subregion inefficiency, the intelligent behavior that is difficult to the anthropomorphic dummy causes subregion rationality difference.
Description of drawings
Fig. 1 model process flow diagram of the present invention.
Fig. 2 particle position of the present invention upgrades arest neighbors rule synoptic diagram.
The soil of Fig. 3 embodiment 1 utilizes subregion present situation figure.
The analog result figure of Fig. 4 embodiment 1.
Embodiment
Model flow process of the present invention such as Fig. 1 show.
This soil utilizes the auto-partition method to comprise the steps:
Step 1. extraction soil utilizes the basic data of subregion and integrates, data after the integration are land use combination and each figure spot information wherein, the extraction basic data is meant from other databases or other system obtains the data that the soil utilizes section post to need, utilize suitability evaluation data, statistical yearbook data etc. as present status of land utilization data, soil, integral data is meant that these data are data of multiple source, different-format, to unify standardization processing to it, form a unified basic database.
Step 2. makes up the soil based on particle swarm optimization algorithm and utilizes the auto-partition model, data with step (1) are model input data, with the land-use map spot is data processing unit, sets up the mapping relations wait to ask between problem and the population, finally finds the solution the soil utilization result that optimizes distribution.The land-use map spot is abstract to be particle in the particle swarm optimization algorithm, and the position of each figure spot and land use pattern correspond to particle position (x, y) and kind (i), the soil utilizes the comprehensive benefit function to be the population fitness function.
Step 3. is provided with the population scale, the maximum flying speed (υ of particle
Xmax, υ
Ymax) and quicken weight (c
1, c
2), maximum iteration time, error amount, each particle position of initialization, speed and attribute;
In the formula: υ
XijAnd υ
YijBe respectively the x of particle and the speed on the y direction of principal axis.
Step 4. is calculated the fitness of each particle according to fitness function.
The fitness function of particle is as follows:
F=(C
k,S
k,Z
k)
In the formula: n is a land-use map spot sum; c
IkThe land use that is i figure spot changes the needed expense of k kind land use pattern into; s
IkIt is the suitability evaluation index of the land use k of i land-use map spot; T
iIt is the set that the adjacent map spot of i land-use map spot is constituted; n
IjCommon edge number for land-use map spot i and figure spot j; l
IjhThe length of side for the h bar common edge of figure spot i and figure spot j; a
iLand area for figure spot i; A
1kAnd A
2kBe respectively required minimum of k kind land use pattern and maximum area; If i land unit is x during land use pattern among the k
Ik=1, otherwise x
Ik=0, C
kBe soil change expense target; S
kIt is the land suitability target of k kind land use; Z
kFor figure spot compactedness is the shape target; F is the multiple goal composite evaluation function.
Step 5. is calculated individual optimal value p
bWith global optimum p
g, and calculate the individual optimal location and the global optimum position of particle, suppose the fitness function determined among the f (X), the optimum position formula of the particle after so each the renewal is:
Hypothetical particle group scale is s, global optimum p
g, then:
P
g(t)∈{P
1b(t),P
1b(t),........,P
sb(t)}=min{f(P
1b(t),f(P
2b(t)),.......,f(P
sb(t)))}
Step 6. is calculated inertia weight, upgrades the flying speed of each particle, draws individual optimal location and global optimum position according to step 5, upgrades the current location of each particle.
The function of inertia weight is: ω (t)=ω
Max-t (ω
Max-ω
Min)/I
Max
In the formula: t is iterations t=1,2 ..., I
Max, I wherein
MaxBe maximum iteration time.ω
MaxBe maximum inertia weight, ω
MinBe minimum inertia weight.
The flying speed renewal function is:
Wherein: " i " expression particle i, the j dimension of " j " expression particle, w is an inertia weight, t represents iterations, c
1, c
2The expression accelerator coefficient, c
1Regulate the step-length of particle self flying speed, c
1Then regulate particle and fly to the step-length of overall desired positions, r
1~U (0,1), r
2~U (0,1) is separate random function, and a moment speed also is influence coefficient to the influence degree of present speed on the expression particle.
In order to prevent the unlimited increase of particle flying speed, present the blast disordered state, need to increase by one group of constraint condition:
Particle maximum flying speed υ
MaxLimited the flying speed υ of particle
i, υ
MaxDetermined particle search accuracy in the space, when value is too big, then particle is crossed optimum solution easily; Work as υ
MaxValue is too little, and then the particle flying speed is slow, is absorbed in the Local Search space easily and can't obtains globally optimal solution.
The current location arest neighbors rule of new particle more, if present bit is equipped with when a plurality of, the arest neighbors rule adopts the area method that is dominant.This rule proposes at proximity relations how to determine scrambling figure spot in the model, is the criterion that is used for judging the next position of particle flight when upgrading particle position.When the current location of known particle and speed, be the center of circle with the figure spot center of current location correspondence, the present speed vector is that radius is done circle, then it next certain figure spot that possible position should be and this circle intersects constantly.When having a plurality of intersection graph spot, at first calculate the area that each intersection graph spot falls into circle and account for the ratio of this figure spot area; the center at figure spot place of getting the ratio maximum is as its next position; as shown in Figure 2; when the current location of known particle and speed; its next possible position should be to be the center of circle (A) with the current location; velocity is on the circle of radius; the dashed circle shown in the accompanying drawing 2; then its next position may be a figure spot 114; 111; 204 or 112; according to the arest neighbors rule of the method that is dominant based on area; getting the center that the area that falls into circle in these figure spots accounts for that figure spot place of figure spot area ratio maximum is its next position, is A1 position among the figure.
Step 7. loop iteration, when satisfying Rule of judgment, the absolute value of the difference of twice global optimum is not more than error amount promptly, and when the number of times of loop iteration is not more than maximum iteration time, search finishes, otherwise the value that will upgrade continues search and obtains global optimum and global optimum position as the input value of step 4.
Step 8. pair step 2) the soil utilization that obtains is optimized distribution, and merges by the figure spot to contiguous identical land use pattern and handles, and extracts figure spot boundary line after merging and promptly obtains the soil and utilize zone boundary, and finally generate the soil and utilize block plan.
Embodiment 1:
1. extract the basic data of a certain small towns soil utilization and integrate, the data after the integration are land use combination and each figure spot information wherein, as shown in Figure 3.
2. use particle swarm optimization algorithm to above-mentioned data-optimized, carry out modeling and select optimum soil to utilize layout.
3. 20 of population scales are set, maximum flying speed of particle (3,4) and acceleration weight (1.44,1.44), maximum iteration time is 200 times, error amount is 10
-6, the position of primary correspondence is determined by the barycentric coordinates of figure spot, the corresponding ten kinds of land use patterns of attribute, and scope is value between (0,10), and initial velocity is that maximum flying speed multiply by the speed that the random number between (0,1) produces at random.
4. according to fitness function, calculate the fitness of each particle.
F=(C
k,S
k,Z
k)
Wherein:
5. calculate individual optimal value p
bWith global optimum p
g, and calculate the individual optimal location and the global optimum position of particle.
6. calculating inertia weight upgrades the flying speed of each particle, draws individual optimal location and global optimum position according to step 5, upgrades the current location of each particle.
Inertia weight w regards the linear function of iterations as, its from 0.98 to 0.48 linearity is set reduces, and maximum inertia weight is 0.98, and minimum inertia weight is 0.48, and the function of inertia weight is:
7. loop iteration.When satisfying Rule of judgment, promptly the absolute value of the difference of twice global optimum in front and back is not more than 10
-6, and the number of times of loop iteration is when being not more than maximum iteration time 200, search finishes, otherwise the value that will upgrade is as the input value of step 4, continues to search for to obtain global optimum and global optimum position.
8. the extraction of optimization result's processing, zone boundary is with definite.To the soil utilization of finishing the loop iteration result figure that optimizes distribution, be close to the merging of similar figure spot and handle, extract different land used categories subareas boundary line, make the soil and utilize subregion figure as a result, as shown in Figure 4.
Claims (6)
1. a soil utilizes the auto-partition method, it is characterized in that this method comprises the steps:
(1) the extraction soil utilizes the basic data of subregion and integrates, and the data after the integration are land use combination and each figure spot information wherein;
(2) make up the soil based on particle swarm optimization algorithm and utilize the auto-partition model, data with step (1) are model input data, with the land-use map spot is data processing unit, sets up the mapping relations wait to ask between problem and the population, finally finds the solution the soil utilization result that optimizes distribution.
(3) the soil utilization that step (2) is obtained is optimized distribution, and merges by the figure spot to contiguous identical land use pattern and handles, and extracts figure spot boundary line after merging and promptly obtains the soil and utilize zone boundary, and finally generate the soil and utilize block plan.
2. soil according to claim 1 utilizes the auto-partition method, it is as follows to it is characterized in that foundation in the step (2) waits to ask the mapping relations between problem and the population: the land-use map spot is abstract to be the particle in the particle swarm optimization algorithm, the position of each figure spot and land use pattern correspond to particle position (x, y) and kind (i), the soil utilizes the comprehensive benefit function to be the population fitness function.
3. soil according to claim 1 utilizes the auto-partition method, it is characterized in that finding the solution in the step (2) drawing optimize distribution result's step of soil utilization and comprising:
(a) scale of population is set, the maximum flying speed (υ of particle
Xmax, υ
Ymax) and quicken weight (c
1, c
2), maximum iteration time, error amount, and each particle position of initialization, speed and attribute;
(b), calculate the fitness of each particle according to fitness function;
(c) calculate individual optimal value (p
b) and global optimum (p
g), and calculate the individual optimal location and the global optimum position of particle;
(d) calculate inertia weight, upgrade the flying speed of each particle, the individual optimal location and the global optimum position that draw according to step (c), the more current location of new particle;
(e) loop iteration: when satisfying Rule of judgment, search finishes, otherwise the value that step (d) is upgraded is as the input value of (b), continues to search for to obtain global optimum and global optimum position.
4. soil according to claim 3 utilizes the auto-partition method, it is characterized in that the current location of the more new particle in the step (c) adopts the arest neighbors rule.
5. soil according to claim 4 utilizes the auto-partition method, it is characterized in that: present bit is equipped with when a plurality of, and the arest neighbors rule adopts the area method that is dominant.
6. soil according to claim 3 utilizes the auto-partition method, and it is characterized in that Rule of judgment is in the step (e): the absolute value of the difference of twice global optimum in front and back is not more than error amount, and the number of times of loop iteration is not more than maximum iteration time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910272814XA CN101877034B (en) | 2009-11-19 | 2009-11-19 | Land utilization automatic division method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910272814XA CN101877034B (en) | 2009-11-19 | 2009-11-19 | Land utilization automatic division method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101877034A true CN101877034A (en) | 2010-11-03 |
CN101877034B CN101877034B (en) | 2012-02-29 |
Family
ID=43019589
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910272814XA Expired - Fee Related CN101877034B (en) | 2009-11-19 | 2009-11-19 | Land utilization automatic division method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101877034B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102521443A (en) * | 2011-12-06 | 2012-06-27 | 河海大学 | Logistics node facility layout optimization method based on computer vision |
CN102567523A (en) * | 2011-12-29 | 2012-07-11 | 武汉大学 | Method for adaptively designing space sampling scheme of geographic elements with aggregated distribution characteristic |
CN103136603A (en) * | 2013-03-26 | 2013-06-05 | 武汉大学 | Intelligent land utilization layout optimal allocation method |
CN103150614A (en) * | 2013-03-26 | 2013-06-12 | 武汉大学 | Automatic configuration method for land utilization space |
CN103150613A (en) * | 2013-03-26 | 2013-06-12 | 武汉大学 | Intelligent optimization method for land utilization layout |
CN106679850A (en) * | 2017-01-03 | 2017-05-17 | 深圳市蜂联科技有限公司 | Method of improving temperature measuring accuracy during the thermal starting process of air box |
CN112801838A (en) * | 2020-12-25 | 2021-05-14 | 广州市城市规划勘测设计研究院 | Urban wetland ecological unit dividing method, device and storage medium thereof |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6608931B2 (en) * | 2001-07-11 | 2003-08-19 | Science Applications International Corporation | Method for selecting representative endmember components from spectral data |
CN1343949A (en) * | 2001-10-26 | 2002-04-10 | 中国科学院合肥智能机械研究所 | Geographic information computer system for intelligent planning of forestry and grassland |
CN1987906B (en) * | 2005-12-22 | 2010-10-13 | 中国煤炭地质总局航测遥感局 | Method for dynamicaly predicting land use change |
CN101256590A (en) * | 2008-04-03 | 2008-09-03 | 北京艺龙天地文化传播有限公司 | Simulation system for three-dimensional on-line virtual reality of environment combining with WebGis and method thereof |
CN101504663B (en) * | 2009-03-17 | 2010-09-08 | 北京大学 | Swarm intelligence based spatial data copy self-adapting distribution method |
-
2009
- 2009-11-19 CN CN200910272814XA patent/CN101877034B/en not_active Expired - Fee Related
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102521443A (en) * | 2011-12-06 | 2012-06-27 | 河海大学 | Logistics node facility layout optimization method based on computer vision |
CN102521443B (en) * | 2011-12-06 | 2013-05-08 | 河海大学 | Logistics node facility layout optimization method based on computer vision |
CN102567523A (en) * | 2011-12-29 | 2012-07-11 | 武汉大学 | Method for adaptively designing space sampling scheme of geographic elements with aggregated distribution characteristic |
CN102567523B (en) * | 2011-12-29 | 2014-02-26 | 武汉大学 | Method for adaptively designing space sampling scheme of geographic elements with aggregated distribution characteristic |
CN103136603A (en) * | 2013-03-26 | 2013-06-05 | 武汉大学 | Intelligent land utilization layout optimal allocation method |
CN103150614A (en) * | 2013-03-26 | 2013-06-12 | 武汉大学 | Automatic configuration method for land utilization space |
CN103150613A (en) * | 2013-03-26 | 2013-06-12 | 武汉大学 | Intelligent optimization method for land utilization layout |
CN103150614B (en) * | 2013-03-26 | 2016-04-13 | 武汉大学 | A kind of Automatic configuration method for land utilization space |
CN103136603B (en) * | 2013-03-26 | 2017-02-22 | 武汉大学 | Intelligent land utilization layout optimal allocation method |
CN103150613B (en) * | 2013-03-26 | 2017-10-03 | 武汉大学 | A kind of intelligent optimization method for land utilization layout |
CN106679850A (en) * | 2017-01-03 | 2017-05-17 | 深圳市蜂联科技有限公司 | Method of improving temperature measuring accuracy during the thermal starting process of air box |
CN112801838A (en) * | 2020-12-25 | 2021-05-14 | 广州市城市规划勘测设计研究院 | Urban wetland ecological unit dividing method, device and storage medium thereof |
Also Published As
Publication number | Publication date |
---|---|
CN101877034B (en) | 2012-02-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101877034B (en) | Land utilization automatic division method | |
CN110458048A (en) | Take population distribution Spatio-temporal Evolution and the cognition of town pattern feature into account | |
CN106529674B (en) | Multiple no-manned plane cooperates with mine to target assignment method | |
CN109493599A (en) | A kind of Short-time Traffic Flow Forecasting Methods based on production confrontation network | |
CN102507586B (en) | Remote sensing monitoring method for carbon emission | |
CN103150614B (en) | A kind of Automatic configuration method for land utilization space | |
CN1987906A (en) | Method for dynamicaly predicting land use change | |
CN106022614A (en) | Data mining method of neural network based on nearest neighbor clustering | |
CN101916337A (en) | Method for dynamically predicting potential productivity of paddy rice based on geographical information system | |
CN104200387A (en) | Land use zoning method giving consideration to geographic entity semantic similarity | |
CN101419623B (en) | Geographical simulation optimizing system | |
Fang et al. | Optimization transmission theory and technical pathways that describe multiscale urban agglomeration spaces | |
CN109920248A (en) | A kind of public transport arrival time prediction technique based on GRU neural network | |
CN102184328A (en) | Method for optimizing land use evolution CA model transformation rules | |
CN104732091A (en) | Cellular automaton river bed evolution prediction method based on natural selection ant colony algorithm | |
CN107748940A (en) | A kind of energy conservation potential Quantitative prediction methods | |
CN103839113A (en) | Microscopic simulation method based on house selecting models of house-renting selectors and house selecting models of house-purchasing selectors | |
Yang et al. | Spatial distribution characteristics and influencing factors of agricultural specialized villages in Guangdong Province, China | |
CN101807218B (en) | Heterogeneous network-based land pattern succession simulation system | |
Cao et al. | Analysis of spatiotemporal changes in cultural heritage protected cities and their influencing factors: Evidence from China | |
CN113743659A (en) | Urban layout prediction method based on component method and Markov cellular automaton and application | |
CN108537581A (en) | Based on the GMDH energy consumption Time Series Forecasting Methods selectively combined and device | |
CN117114176A (en) | Land utilization change prediction method and system based on data analysis and machine learning | |
CN105488346A (en) | Spatial prediction analogy method of large-scale land cover change | |
Huang et al. | A Novel Hybrid Model for PM2. 5 Concentration Forecasting Based on Secondary Decomposition Ensemble and Weight Combination Optimization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120229 Termination date: 20151119 |
|
EXPY | Termination of patent right or utility model |