CN101630007B - Method for modifying navigation results of satellite remote sensing images under fault condition of K mirror - Google Patents
Method for modifying navigation results of satellite remote sensing images under fault condition of K mirror Download PDFInfo
- Publication number
- CN101630007B CN101630007B CN2009100853179A CN200910085317A CN101630007B CN 101630007 B CN101630007 B CN 101630007B CN 2009100853179 A CN2009100853179 A CN 2009100853179A CN 200910085317 A CN200910085317 A CN 200910085317A CN 101630007 B CN101630007 B CN 101630007B
- Authority
- CN
- China
- Prior art keywords
- data
- remote sensing
- sensing images
- mirror
- satellite remote
- 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
Links
Images
Landscapes
- Image Input (AREA)
- Image Processing (AREA)
- Processing Or Creating Images (AREA)
Abstract
The invention provides a method for modifying the navigation results of satellite remote sensing images under the fault condition of a K mirror, which used for the condition that a satellite scans by means of multivariate and parallelism. The method comprise the following steps of: obtaining the satellite remote sensing images and image navigation information; obtaining a mapping relationship among pixel row, line coordinate and a modified image navigation result by an unmodified image navigation result; and modifying the navigation results of the satellite remote sensing images by the mapping relationship. The invention can exactly modify an image navigation error caused by the fault of the K mirror, and has a significant applying value in processing remote satellite data.
Description
Technical field
The present invention relates to a kind of method that the navigation results of satellite remote sensing images is revised, relate in particular under the situation that a kind of K mirror in remote sensing satellite breaks down, utilize uncorrected image-guidance result to obtain pixel row, column coordinate and revise mapping relations between the image-guidance result of back, method by these mapping relations are revised the navigation results enforcement of satellite remote sensing images belongs to the remote sensing technology field.
Background technology
Satellite remote sensing images is meant that remote sensing satellite utilizes remote sensing means such as visible light, ultraviolet ray, infrared ray or microwave to survey ground object, the target image that the remotely-sensed data that is obtained is generated after treatment.It is the significant data source that people obtain terrestrial information, is subjected to remotely-sensed data user's great attention.
In spacer remote sensing,, adopt the scanning mirror scheme usually in order to enlarge the observation scope of satellite remote sensing images.45 degree rotation sweep mirror-mode are exactly a kind of common optical mechaical scanning pattern, because its scanning mirror size is little, good stability, fabric width are big, are widely used in the spaceborne and airborne optical spectrum imagers.Along with deepening continuously that remotely-sensed data is used, also more and more higher to the detection sensitivity requirement of remote sensor.In order to pursue higher resolution, the detector of current remote sensing instrument adopts polynary linear array detector usually, go up the intermediate-resolution optical spectrum imagers (MERSI) that carry as China's polar orbiting meteorological satellite wind and cloud No. three (FY-3), adopted the also imaging scheme of line scanning of 10 yuan and 40 yuan.But when the remote sensing instrument that adopts polynary linear array detector changed in the scanning view angle, the picture on the focal plane did not allow to rotate, otherwise can make the visual field of each pixel correspondence can't registration, caused pattern distortion.
Can't registration problems in order to solve the outer visual field of axle that 45 degree rotation sweep catoptrons cause, people use K illuminating apparatus structure to come the phenomenon (being called for short as revolving) of removal of images rotation.45 degree rotation sweep catoptrons and K mirror constitute optical image racemization system jointly.The K mirror is an all band reflecting system, is made up of 3 plane mirrors, when its corner remains a half of 45 degree rotation sweep catoptron corners, as not rotating, reaches the purpose of eliminating the picture rotation.Because of its thrihedral reflector is " K " word arrangement, so be named as the K mirror.
When application K mirror elimination picture revolves,, break down easily because the K mirror is a moving component.When the K mirror breaks down, need adopt the mode of Flame Image Process to carry out pin on ground and revolve processing.At Granted publication number is in the Chinese invention patent of CN1169092C, discloses a kind of method of handling the picture rotation that 45 degree scanning mirrors cause, comprises following steps: to the computing machine input as image rotating; Computing machine carries out imaging relations and calculates; Derive the imaging process of whole optics remote sensing instrument from the most basic reflection law; Acquisition is as the mapping relations between image rotating coordinate and the racemization image coordinate; Utilize coordinate Mapping relation and racemization algorithm that entire image is carried out racemization; The output of racemization image; As input, circulation is carried out with the image of above-mentioned racemization; When the K mirror broke down, what also comprise the K mirror angle obtained and imports step in this computing machine.Closing quick and Guo Qiang is published in " applicating atmosphere journal " paper " application of optical image racemization system in FY-3 MERSI framing " of the 19th the 4th phase of volume, by studying the rotary reflection characteristic of 45 degree rotation sweep catoptrons and K mirror, the object-image relation of 45 degree rotation sweep catoptrons and the reflection characteristic matrix of K mirror have been provided, and analyzed source of error, in image geo-location result, eliminated error outside the axle.In addition, be published in " ocean journal " paper " " HY-1A " satellite ocean water colour scanner is as the image rotating treatment for correcting " of the 30th the 1st phase of volume at Ceng Tao etc., based on the optical reflection vector theory, analyzed the track while scan feature that 45 degree scanning mirrors under the parallel scan mode of multiunit detector add the K mirror, utilize the correlativity between the image, at first automatically, seek the accurate angle of K mirror stall, according to the stall angle of K mirror and polynary and sweep and theoretically picture is revolved image carried out treatment for correcting, obtained available correcting image then.
But, implements above-mentioned image pin and revolve operation and will change the row, column coordinate of pixel in image, and the result of satellite remote sensing images navigation is one to one with the row, column coordinate of pixel, unique warp, the latitude information of each pixel correspondence.Therefore, be index with warp, latitude at this moment, satellite remote sensing images is waited the longitude and latitude projection, will produce mistake.As shown in Figure 1, typical satellite remote sensing images will produce dislocation, fracture, drop-out, striped or the like abnormality, and this is owing to the navigation results (i.e. Dui Ying warp, latitude information) of each pixel not being revised and caused.
Summary of the invention
Technical matters to be solved by this invention is to provide a kind of method of the K of correction mirror navigation results of satellite remote sensing images under fault condition.Image-guidance result before this method utilization correction carries out after the image pin revolves processing on ground, by pixel row, column coordinate and the mapping relations revised between the image-guidance result of back navigation results is revised.
For realizing above-mentioned purpose, the present invention adopts following technical scheme:
A kind of method of revising K mirror navigation results of satellite remote sensing images under fault condition is used in satellite and carries out it is characterized in that under the situation that polynary and line scanning and K mirror break down:
(1) once reads remote sensing image data double and that line scanning obtained;
(2) in described remote sensing image data, with the longitude of each row pixel as component, form one group of data, the most last partial data of the last part data of single pass and back single pass carries out curve fitting before selecting in these group data, and carries out the resampling of partial data in these group data;
(3) repeating step (2) finishes until each row processes pixel;
(4) in described remote sensing image data, with the latitude of each row pixel as component, form one group of data, the most last partial data of the last part data of single pass and back single pass carries out curve fitting before selecting in these group data, and carries out the resampling of partial data in these group data;
(5) repeating step (4) finishes until each row processes pixel;
(6) utilize the data after resampling that navigation results of satellite remote sensing images is revised.
Wherein, in described step (1), described remote sensing image data obtains in the following way: at first carry out once and line scanning, produce the multirow data; For each line data, utilize sequential circuit repeatedly to sample, the data of each sampling are as the data of each row.
In described step (2) and (4), the data that are used to carry out curve fitting are no less than half of these group data.The most last partial data of the last part data of single pass and back single pass before described partial data comprises, the last part data and the most last partial data that resample are no less than 1/5th of single pass fetched data respectively, 1/3rd of no more than single pass fetched data.
Described method is when implementing, and the zone of reasonableness of definition longitude is [180,180], and the zone of reasonableness of latitude is [90,90].
During near handling date line data, handle again after each longitude data added 180 respectively.
Navigation results of satellite remote sensing images modification method provided by the present invention utilizes uncorrected image-guidance result, obtains the mapping relations between pixel row, column coordinate and the correction back image-guidance result; Utilize these mapping relations that the image-guidance result is revised.Utilize the present invention, can accurately revise because the caused image-guidance error of K mirror fault has important use and is worth in the work of treatment of remote sensing satellite data.
Description of drawings
The present invention is further illustrated below in conjunction with the drawings and specific embodiments.
Fig. 1 the satellite remote sensing images synoptic diagram of inconsistent phenomenon occurs for before revising;
Fig. 2 reads the synoptic diagram of twice scan-data in this satellite remote sensing images navigation modification method at every turn;
Fig. 3 forms the longitude (latitude) of each row pixel the synoptic diagram of one group of parameter in this satellite remote sensing images navigation modification method as component;
Fig. 4 adopts the synoptic diagram of curve-fitting method in this satellite remote sensing images navigation modification method;
Fig. 5 is in this satellite remote sensing images navigation modification method, and longitude (latitude) schematic diagram data afterwards resamples;
Fig. 6 be revise after, eliminated the satellite remote sensing images synoptic diagram of inconsistent phenomenon.
Embodiment
There are relation one to one in the row, column coordinate of each pixel in the satellite remote sensing images and geographical longitude and latitude, elevation information.The data such as scan mode of track via satellite, attitude, satellite earth observation instrument can obtain pairing geographical longitude and latitude of each pixel and height in the satellite remote sensing images by calculating, and this computation process just is called image-guidance and calculates.
Before address, the detector that carries on the remote sensing satellite adopts polynary linear array detector to carry out imaging usually at present, the intermediate-resolution optical spectrum imagers (MERSI) on China's polar orbiting meteorological satellite wind and cloud No. three (FY-3) for example, 1 kilometer resolution channels of this instrument adopts 10 yuan of detectors, scan by 45 degree scanning mirrors, in light path, utilize the K mirror to carry out the light path racemization.Under the situation that the K mirror breaks down, the remote sensing images that form behind the satellite scanning imagery picture can occur and revolve.If the K mirror does not break down, this phenomenon will can not take place.Revolve in order to overcome picture, carry out racemization usually earlier and handle, carry out image-guidance again and calculate.
Deviation occurs for fear of the image-guidance result calculated, the present invention proposes a kind of method of the K of correction mirror navigation results of satellite remote sensing images under fault condition.Below, be example with No. three (FY-3) 10 yuan of detectors that weather satellite was adopted of above-mentioned wind and cloud, the concrete implementation step of this modification method is launched detailed explanation.
Because the intermediate-resolution optical spectrum imagers (MERSI) in No. three (FY-3) weather satellite of wind and cloud is 10 yuan of 1 kilometer passage employings and sweep scheme, this instrument carries out single pass will produce 10 line data.For each row of data, utilize sequential circuit to sample 2048 times.The data of each sampling are as the data of each row.Thus, the image scanning of intermediate-resolution optical spectrum imagers once will produce the satellite image of 10 * 2048 pixel sizes.The image of this 10 * 2048 size is called " single pass ".
This modification method at first reads uncorrected satellite remote sensing images navigation information.As shown in Figure 2, computing machine reads the data of twice scanning at every turn, promptly has warp, the latitude information of 20 * 2048 pixels altogether.With the example that is treated to of longitude, the longitude of each row pixel as component, is formed one group of data, two parts data that single pass and back single pass are obtained before comprising in these group data.As shown in Figure 3, be designated as lon1 herein, lon2, lon3 ..., lon20.Because the size of the image that twice scanning obtained is 20 * 2048, therefore such data have 2048 groups.
The calculating of the relation of next navigating is so that obtain pixel row, column coordinate and revise the mapping relations between the image-guidance result afterwards from uncorrected navigation results.
Particularly, each group data of gathering for above-mentioned steps, at first the most last partial data of selecting the last part data of last scanning and scanning next time adopts curve-fitting method shown in Figure 4 to carry out match, obtains curve fitting equation y=f (x).Should select data as much as possible to carry out curve fitting in each group data, therefore the data that carry out curve fitting should be no less than half of every group of data, can comprise the total data of this group the most for a long time.When carrying out curve fitting, can select Chebyshev's curve fitting or least square curve fit etc.According to the experience of practical operation, preferably adopt Chebyshev's curve-fitting method.
After having carried out curve fitting, next carry out the resampling of partial data.The process that resamples is such: for each group data, the most last partial data of the last part data of single pass and back single pass resamples according to the curvilinear equation that aforementioned curve fitting obtained before selecting.The data area that resamples is unsuitable excessive or too small.Generally speaking, above-mentioned last part data and the most last partial data should not be less than 1/5th of single pass fetched data respectively, also should not be more than 1/3rd of single pass fetched data.Self-evident, the data bulk that carries out curve fitting and resample should round numbers.As shown in Figure 5, for 10 yuan and sweep for the scheme, can select last three pixel lon8 of last scanning, lon9, first three pixel lon11 of lon10 and scanning next time, lon12, lon13 resamples, and just according to y=f (x), recomputates the 8th in these 20 pixels, 9,10,11,12, the longitude of 13 pixels (counting of pixel is since 1), obtain new longitude newlon8, newlon9, newlon10 like this, newlon11, newlon12, newlon13, these new longitude are formed one group of new parameter, lon1 with original unchanged longitude, lon2 ..., lon7, newlon8, newlon9, newlon10, newlon11, newlon12, newlon13, lon14, lon15,, lon20.This group parameter is exactly to revise the longitude of this row pixel afterwards.
Processing procedure for latitude is similarly, has not just given unnecessary details at this.
Because the data that twice scanning obtained have 2048 groups, therefore such correction is carried out 2048 groups altogether.The correction of these 2048 groups of data can once walk abreast and finish, and therefore once can revise the image block of 20 * 2048 sizes.Circulate with this, dispose up to entire image.Because 1 kilometer resolution channels image size of the intermediate-resolution optical spectrum imagers of No. three (FY-3) weather satellite of wind and cloud is 2000 * 2048, therefore circulates and to finish the processing of entire image for 100 times.
For 40 yuan and sweep for the scheme, can select last 10 pixel lon31 of last scanning, lon32,, lon39, lon40 and the most preceding 10 the pixel lon41 that next time scan, lon42,, lon50 resamples, and with the longitude of revised numerical value as this row pixel.
For the polynary of other and sweep for the scheme, handling thinking also is similarly, promptly selects the most last partial data of the last part data of last scanning and scanning next time to carry out curve fitting and resample, and has not just given unnecessary details one by one at this.
Need to prove that there is a problem in above-mentioned treatment step: when the longitude information variation was violent, resampling was inaccurate.According to research, longitude information changes violent zone and is limited to two places: the first meridian and date line.For example chosen longitude data variation behind radix point, in the time of then near date line, the process that longitude changes is-179.0 ,-179.1 ... ,-179.9,179.9,179.8 ..., 179.0.If these group data are directly carried out match,, will make curve fitting produce very large error because-179.9,179.9 these variations to adjacent data are very violent.For this reason, the present invention has carried out following processing: for whole tellurian point, the zone of reasonableness of definition longitude is [180,180].In this range of definition, the saltus step of longitude information can not appear near the first meridian.When near the data the processing date line, each data in [180,180] are added 180 separately, thereby they are transformed in [0,360] scope, can make matched curve become very level and smooth like this, mistake can not appear.This matched curve is used for resampling, can obtains correction result more accurately.So far, longitude disposes.The processing of latitude information is the same with longitude.Because the zone of reasonableness of definition latitude is [90,90], so the phenomenon of saltus step can not appear in latitude, do not need to handle once more.
After resampling, the numerical value after resampling is substituted the original numerical value that does not resample, the image-guidance result after just can obtaining to revise.Shown in Figure 6 be revise after, eliminated the satellite remote sensing images synoptic diagram of inconsistent phenomenon.Before correction, phenomenon of rupture appears in the many places, river, and image has tangible striped, and after using this method and revising, phenomenon of rupture disappears, and striped has not existed yet.Therefore these remote sensing images navigation results modification method can be brought into play good effect in the processing procedure of remote sensing satellite data.
It is emphasized that these remote sensing images navigation results modification method handled be the original navigation data of satellite remote sensing images, be different from operations such as using business software is carried out the rotation of simplicial graph picture, moved, stretching fully.The present invention by method corrections such as curve fitting and resamplings the longitude and latitude attribute of pixel, make that the longitude and latitude attribute of pixel is correct, utilize the longitude and latitude of pixel to carry out projection then, so makeover process new mistake can not occur.
More than the method for correction of the present invention K mirror navigation results of satellite remote sensing images under fault condition is had been described in detail.For one of ordinary skill in the art, any conspicuous change of under the prerequisite that does not deviate from connotation of the present invention it being done all will constitute to infringement of patent right of the present invention, with corresponding legal responsibilities.
Claims (8)
1. method of revising K mirror navigation results of satellite remote sensing images under fault condition is used in satellite and carries out it is characterized in that under the situation that polynary and line scanning and K mirror break down:
(1) once reads remote sensing image data double and that line scanning obtained;
(2) in described remote sensing image data, with the longitude of each row pixel as component, form one group of data, the most last partial data of the last part data of single pass and back single pass carries out curve fitting before selecting in these group data, and carries out the resampling of partial data in these group data;
(3) repeating step (2) finishes until each row processes pixel;
(4) in described remote sensing image data, with the latitude of each row pixel as component, form one group of data, the most last partial data of the last part data of single pass and back single pass carries out curve fitting before selecting in these group data, and carries out the resampling of partial data in these group data;
(5) repeating step (4) finishes until each row processes pixel;
(6) utilize the data after resampling that navigation results of satellite remote sensing images is revised.
2. the method for correction K mirror navigation results of satellite remote sensing images under fault condition as claimed in claim 1 is characterized in that:
In the described step (1), described remote sensing image data obtains in the following way: at first carry out once and line scanning, produce the multirow data; For each line data, utilize sequential circuit repeatedly to sample, the data of each sampling are as the data of each row.
3. the method for correction K mirror navigation results of satellite remote sensing images under fault condition as claimed in claim 1 is characterized in that:
In described step (2) and (4), the data that are used to carry out curve fitting are no less than half of these group data.
4. the method for correction K mirror navigation results of satellite remote sensing images under fault condition as claimed in claim 1 is characterized in that:
In described step (2) and (4), described curve fitting adopts Chebyshev's curve-fitting method.
5. the method for correction K mirror navigation results of satellite remote sensing images under fault condition as claimed in claim 1 is characterized in that:
In described step (2) and (4), described partial data is resampled according to the curvilinear equation that curve fitting obtained.
6. the method for correction K mirror navigation results of satellite remote sensing images under fault condition as claimed in claim 5 is characterized in that:
The most last partial data of the last part data of single pass and back single pass before described partial data comprises, the described last part data and the most last described partial data that are used to resample are no less than 1/5th of single pass fetched data respectively, 1/3rd of no more than single pass fetched data.
7. the method for correction K mirror navigation results of satellite remote sensing images under fault condition as claimed in claim 1 is characterized in that:
Described method is when implementing, and the zone of reasonableness of definition longitude is [180,180], and the zone of reasonableness of latitude is [90,90].
8. the method for correction K mirror navigation results of satellite remote sensing images under fault condition as claimed in claim 7 is characterized in that:
During near handling date line data, handle again after each longitude data added 180 respectively.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100853179A CN101630007B (en) | 2009-05-26 | 2009-05-26 | Method for modifying navigation results of satellite remote sensing images under fault condition of K mirror |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100853179A CN101630007B (en) | 2009-05-26 | 2009-05-26 | Method for modifying navigation results of satellite remote sensing images under fault condition of K mirror |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101630007A CN101630007A (en) | 2010-01-20 |
| CN101630007B true CN101630007B (en) | 2011-12-28 |
Family
ID=41575181
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2009100853179A Expired - Fee Related CN101630007B (en) | 2009-05-26 | 2009-05-26 | Method for modifying navigation results of satellite remote sensing images under fault condition of K mirror |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101630007B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112130180B (en) * | 2020-09-18 | 2023-10-27 | 上海钧正网络科技有限公司 | Method, device and server for detecting faults of satellite positioning equipment of shared equipment |
| CN114812603B (en) * | 2022-03-28 | 2024-09-17 | 国家卫星气象中心(国家空间天气监测预警中心) | Method for navigating and positioning and checking static satellite optical image by using ground laser |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101246249A (en) * | 2007-02-13 | 2008-08-20 | 中国科学院上海技术物理研究所 | Three-dimensional regulating mechanism for fixing offsetting image rotation K mirror of spacing remote sensing instrument |
-
2009
- 2009-05-26 CN CN2009100853179A patent/CN101630007B/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101246249A (en) * | 2007-02-13 | 2008-08-20 | 中国科学院上海技术物理研究所 | Three-dimensional regulating mechanism for fixing offsetting image rotation K mirror of spacing remote sensing instrument |
Non-Patent Citations (2)
| Title |
|---|
| 关敏等.光学像消旋系统在FY-3 MERSI图像定位中的应用.《应用气象学报》.2008,第19卷(第4期),正文第420页到第428页. * |
| 张强等.消除FY-4辐射计图像旋转的方法研究.《红外技术》.2007,第29卷(第6期),正文第365页到第369页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101630007A (en) | 2010-01-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Tang et al. | Triple linear-array image geometry model of ZiYuan-3 surveying satellite and its validation | |
| CN107204037B (en) | Three-dimensional image generation method based on active and passive three-dimensional imaging system | |
| CN106767714B (en) | Improve the equivalent mismatch model multistage Calibration Method of satellite image positioning accuracy | |
| CN107144293A (en) | A kind of geometric calibration method of video satellite area array cameras | |
| CN109696182A (en) | A kind of spaceborne push-broom type optical sensor elements of interior orientation calibrating method | |
| CN105913435B (en) | A kind of multiscale morphology image matching method and system suitable for big region | |
| CN106679676B (en) | A kind of monoscopic multifunctional optical sensor and implementation method | |
| CN106887016B (en) | Automatic relative registration method for GF-4 satellite sequence images | |
| CN112432705B (en) | Multispectral imaging system and method based on dynamic visual axis adjustment principle | |
| CN114972545B (en) | On-orbit data rapid preprocessing method for hyperspectral satellite | |
| CN104361563B (en) | GPS-based (global positioning system based) geometric precision correction method of hyperspectral remote sensing images | |
| Liu et al. | Robust radiometric normalization of multitemporal satellite images via block adjustment without master images | |
| CN108447100B (en) | Method for calibrating eccentricity vector and visual axis eccentricity angle of airborne three-linear array CCD camera | |
| Murtiyoso et al. | Technical considerations in Low-Cost heritage documentation | |
| CN101630007B (en) | Method for modifying navigation results of satellite remote sensing images under fault condition of K mirror | |
| CN110926501B (en) | Automatic calibration method and system for optical measurement equipment and terminal equipment | |
| Qu et al. | Retrieval of 30-m-resolution leaf area index from China HJ-1 CCD data and MODIS products through a dynamic Bayesian network | |
| Cheng et al. | New on-orbit geometric interior parameters self-calibration approach based on three-view stereoscopic images from high-resolution multi-TDI-CCD optical satellites | |
| Barreau et al. | Photogrammetry based study of ceramics fragments | |
| CN111199557A (en) | Quantitative analysis method and system for decay of remote sensor | |
| CN118154420B (en) | Multi-feature-based global mosaic generation method and system for stationary meteorological satellite | |
| CN106558036B (en) | A kind of spaceborne super-resolution imaging design method | |
| CN104143179B (en) | Method for enhancing moving target through multi-linear-array time difference scanning expansion sampling | |
| CN111161186B (en) | Push-broom type remote sensor channel registration method and device | |
| CN109029496A (en) | One kind being suitable for large area array optical camera single game eradiation calibrating method |
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 | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20111228 Termination date: 20140526 |