CN102141617B - Microwave staring imaging correlation method - Google Patents

Microwave staring imaging correlation method Download PDF

Info

Publication number
CN102141617B
CN102141617B CN201110000698A CN201110000698A CN102141617B CN 102141617 B CN102141617 B CN 102141617B CN 201110000698 A CN201110000698 A CN 201110000698A CN 201110000698 A CN201110000698 A CN 201110000698A CN 102141617 B CN102141617 B CN 102141617B
Authority
CN
China
Prior art keywords
field
radiation field
random
imaging
space
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
CN201110000698A
Other languages
Chinese (zh)
Other versions
CN102141617A (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.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
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 University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN201110000698A priority Critical patent/CN102141617B/en
Publication of CN102141617A publication Critical patent/CN102141617A/en
Application granted granted Critical
Publication of CN102141617B publication Critical patent/CN102141617B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a microwave staring imaging correlation method which comprises the following stepsof irradiating a target area object through a random radiation field and forming the distribution of the radiation field on an object surface; receiving a scattering field through a multi-channel or a single channel and carrying out electromagnetic field analysis to the imaging physical process of the target area object to obtain a time and space two-dimensional radiation field function corresponding to the object surface within the whole imaging area scope; counting the relevant characteristics of a time and space two-dimensional radiation field according to the time and space two-dimensional radiation field function; and carrying out correlated processing on the scattering field and the random radiation field according to the statistic characteristics to obtain the imaging of the target area object through inversion. The microwave staring imaging correlation method can be combined with the scattering field and the time and space two-dimensional radiation field to realize the high-resolution target imaging.

Description

Microwave staring imaging correlating method
Technical field
The invention belongs to Radar Technology, spacer remote sensing and precise guidance field, specifically, the present invention relates to the association process technology in a kind of microwave staring imaging.
Background technology
Traditional microwave is stared radar imagery and is mainly adopted real aperture technique.Real aperture imaging radar serves as typical case's representative with the beam scanning radar; Its imaging processing technology dependence is carried out the relevant matches realization with the signal that scattering receives signal and emissive source; Because in this process, the signal of emissive source also shall not be applied to the spatial discrimination of target, thereby the spatial resolution of target only depends on the steric configuration of emissive source array; Promptly only by the beam angle decision of emissive source array on the traditional sense, the illuminated target detail in the wave beam then can't be realized differentiating.
Imaging processing technology as wherein a kind of focal plane imaging radar of microwave staring imaging then realizes through the locus corresponding relation of establishing between each receiving array of focal plane and the target.Because the serious diffraction of focusing anteena in the focal plane imaging radar, make to have comprised irreversible fuzzification process in the above-mentioned corresponding relation, and this process is by limited the causing of spatial frequency spectrum coverage of this diffraction system transport function.Thereby the spatial resolution of the target of focal plane imaging radar is by the diffraction limit decision that is directly proportional with the focusing anteena pore size, and is also irrelevant with the signal form of emissive source.
As indicated above, traditional microwave is stared radar imagery and is generally adopted real aperture imaging radar or focal plane imaging radar.No matter be to serve as the real aperture imaging radar of typical case's representative with the beam scanning radar; Or focal plane imaging radar; Its radiation field all need satisfy the total correlation characteristic on time and space; Its spatial discrimination mainly passes through the spatial spread realization of radar array, and irrelevant with radiation field, thereby its imaging processing process all only needs directly handle scattered field.
Therefore, the microwave staring imaging in existing traditional real aperture technology is only handled to received signal, thereby causes the imaging of low resolution.
Summary of the invention
The object of the invention is intended to one of solve the aforementioned problems in the prior at least.
For this reason, embodiments of the invention propose a kind of scattered field and space-time bidimensional high-resolution microwave staring imaging correlating method of radiation field at random of uniting.
According to an aspect of the present invention, the embodiment of the invention has proposed a kind of microwave staring imaging correlating method, and said correlating method may further comprise the steps: through radiation field at random the target area object is shone, on object plane, form radiation field distribution; Carry out scattered field through hyperchannel or single channel and receive, and the imaging physical process of target area object is carried out electromagnetic field analysis, to obtain the space-time bidimensional radiation field function in the corresponding whole imaging region scope of object plane; Add up the space-time bidimensional correlation properties of radiation field at random according to said space-time bidimensional radiation field function; And according to said statistical property, to scattered field and at random radiation field carry out association process obtains the target area object with inverting imaging.
The further embodiment according to the present invention, the space-time bidimensional at random the statistical property of radiation field satisfy the single order field intensity when irrelevant, to scattered field and at random radiation field carry out the field intensity relevant treatment.
The further embodiment according to the present invention, the space-time bidimensional at random the statistical property of radiation field satisfy high-order when irrelevant, to scattered field and at random radiation field carry out the relevant treatment of high-order field amount.
The further embodiment according to the present invention, the space-time bidimensional at random the statistical property of radiation field satisfy single order or high order part when relevant, through the mode that radiation field is at random directly inverted scattered field is handled.Wherein, the space-time bidimensional at random the statistical property of radiation field satisfy the relevant and gaze duration of single order or high order part when too short, utilize the mode of separating optimum solution to scattered field and at random radiation field handle.
The present invention unites scattered field and space-time bidimensional radiation field at random; Change with space diversity at random in time through radiation field at random; According to the space-time bidimensional difference of the statistical dependence characteristic of radiation field at random; And adopt different processing modes, realized that traditional microwave is stared the high-resolution imaging in the radar beam under the same steric configuration.
Aspect that the present invention adds and advantage part in the following description provide, and part will become obviously from the following description, or recognize through practice of the present invention.
Description of drawings
Above-mentioned and/or additional aspect of the present invention and advantage are from obviously with easily understanding becoming the description of embodiment below in conjunction with accompanying drawing, wherein:
Fig. 1 is the microwave staring imaging correlating method process flow diagram of the embodiment of the invention;
Fig. 2 is the imaging scene synoptic diagram of the radiation source at random of the embodiment of the invention;
Fig. 3 (a) and Fig. 3 (b) are respectively the time domain waveform and the frequency-domain waveform of desirable white Gaussian noise;
Fig. 4 is the radiation style of the radiation source at random of the embodiment of the invention at certain time point;
Fig. 5 (a) and Fig. 5 (b) are respectively the cross correlation function between autocorrelation function and of certain point on the radiation field at 2;
Fig. 6 is original imageable target scatter intensity distribution;
When Fig. 7 is the gaze duration long enough, the embodiment of the invention carry out the imaging results of field intensity relevant treatment based on scattered field and space-time bidimensional radiation field;
Fig. 8 is for gaze duration is sampled as 2000, signal to noise ratio (S/N ratio) is under the 32dB situation, the embodiment of the invention based on the target imaging result who directly inverts;
Fig. 9 is for gaze duration is sampled as 500, signal to noise ratio (S/N ratio) is under the 32dB situation, the embodiment of the invention based on the target imaging result who directly inverts; And
Figure 10 is for gaze duration is sampled as 500, signal to noise ratio (S/N ratio) is under the 32dB situation, the embodiment of the invention based on the target imaging result who separates optimization problem.
Embodiment
Describe embodiments of the invention below in detail, the example of said embodiment is shown in the drawings, and wherein identical from start to finish or similar label is represented identical or similar elements or the element with identical or similar functions.Be exemplary through the embodiment that is described with reference to the drawings below, only be used to explain the present invention, and can not be interpreted as limitation of the present invention.
With reference now to Fig. 1,, Fig. 1 is the microwave staring imaging correlating method process flow diagram of the embodiment of the invention.This figure has provided and the present invention is based on the space-time bidimensional association process technology in the high-resolution microwave staring imaging of radiation field at random, below in conjunction with the electromagnetic field ultimate principle, detailed has provided the processing details in each concrete steps among the present invention.Specifically describe as follows:
At first, through radiation field at random the irradiation of target area object is formed radiation field distribution (step 101) on object plane.
Carrying out scattered field through hyperchannel or single channel then receives.Above-mentioned imaging physical process to the target area object is carried out electromagnetic field analysis, obtains the interior space-time bidimensional of the corresponding whole imaging region scope of the object plane distribution function (step 102) of radiation field at random.
Suppose that special actinal surface field coverage is S '; The physical centre of bore face is a true origin; Receive the locus and be
Figure BDA0000042696060000041
scattered field on
Figure BDA0000042696060000042
actinal surface field more arbitrarily, this radiation signal form in the t moment is:
Figure BDA0000042696060000043
wherein is the coordinate vector of this point.
Figure BDA0000042696060000045
arranged in addition is the backscattering coefficient of
Figure BDA0000042696060000046
position on the imaging region; The scope of imaging region is S, then the space-time bidimensional at random radiation field can be expressed as:
Figure BDA0000042696060000047
Scattered field can be expressed as:
Then, add up the space-time bidimensional correlation properties (step 103) of radiation field at random according to above-mentioned space-time bidimensional radiation field function.
And according to the difference of above-mentioned space-time bidimensional radiation field statistics correlation properties at random, to scattered field and at random radiation field carry out association process, obtain the imaging (step 104) of target area object with inverting.
Specifically; In the real aperture radar imaging of tradition; The statistical property of radiation field satisfies the total correlation on any rank; Then fully different with it among the present invention, the space-time bidimensional among the present invention statistical property of radiation field at random only two kinds of situation possibly occur: irrelevant and part correlation according to the difference of correlativity.In addition, different according to the object of space-time bidimensional radiation field statistics correlation properties at random, can be divided into first-order characteristics and high-order characteristic again.Therefore, the imaging processing among the present invention is divided into following three kinds of patterns (easy for writing, two kinds of situation that will be roughly the same, single order part correlation and high order part be relevant merges into a pattern) here:
1, when the space-time bidimensional at random the statistical property of radiation field satisfy single order (being field intensity) when irrelevant, corresponding association process algorithm need to scattered field and space-time bidimensional at random radiation field carry out field intensity relevant treatment (<,>is the inner product operation of space-time bidimensional) here:
Figure BDA0000042696060000049
2, when the space-time bidimensional at random the statistical property of radiation field satisfy high-order (like second order time be intensity) when irrelevant, corresponding association process algorithm need to scattered field and space-time bidimensional at random radiation field carry out the relevant treatment of high-order field amount.Satisfy second order when irrelevant like the statistical property of space-time bidimensional radiation field, inversion algorithm can be expressed as: (here, <,>is the inner product operation of space-time bidimensional):
3, when the space-time bidimensional at random the statistical property of radiation field satisfy single order or high order part when relevant; Corresponding association process algorithm can directly utilize the method that radiation field is at random inverted that scattered field is handled, also available method of separating optimization problem to scattered field and at random radiation field handle.
The method of directly radiation field at random being inverted is mainly used under the sufficient situation of gaze duration, under this situation, adopts the method for directly inverting can obtain higher imaging space resolution.In addition, this algorithm is very simple, is specially adapted to the scene that real-time is had relatively high expectations.
And when gaze duration was too short, this moment was owing to receive the deficiency of signal sampling, and there is serious distortion in the imaging results that the method that can cause directly inverting obtains, thereby the preferred method of separating optimization problem that adopts is handled.
With the single order is example; Directly utilize method processing procedure that radiation field is at random inverted (wherein as
Figure BDA0000042696060000052
not during full rank,
Figure BDA0000042696060000053
is the generalized inverse of
Figure BDA0000042696060000054
) as follows:
Figure BDA0000042696060000055
The step of separating optimization problem is for finding the solution following optimized problem:
Figure BDA0000042696060000056
Wherein, λ >=0 is chosen and is looked signal to noise ratio (S/N ratio) and decide; When noise was non-existent, the problem of minimum L1 norm promptly only need be found the solution in λ=0.
Below in conjunction with concrete imaging scene; Through accompanying drawing and embodiment principle of the present invention is described in detail, because the imaging scene face that the present invention relates to is wider, thereby; Be example only below with a comparatively typical instantiation; Provide embodiment of the present invention, still, should not limit practical application of the present invention and protection domain with this.
For example, concrete imaging scene is as shown in Figure 2, sets up true origin in the physical centre of bore face, and the radiation source at random on the bore face is made up of 4 point sources, and its spatial arrangement is as shown in Figure 2, and 4 emissive sources constitute a square, and neighbor distance is D=10m.Also have a single-point receiver in addition, its coordinate is (0,0,0), and the sampling rate of receiver is f s=1GHz.
Imaging plane is to be positioned at and emission array and receiver vertical range Z 0On the plane of=300m, be the x-y plane among the figure, imageable target is as shown in Figure 6, is divided into 40 * 40 resolution element, its coordinate range be (0m, 0m)~(D Obj=120m, D Obj=120m); Each some emissive source on the bore face all adopts desirable white Gaussian noise as transmitting; Its time domain waveform and frequency-domain waveform respectively shown in 3 (a) and Fig. 3 (b), the radiation field amplitude distribution on the imaging region that obtains (at certain time point of confirming) as shown in Figure 4, the cross correlation function on the radiation field between the autocorrelation function of certain point and any 2 is respectively shown in Fig. 5 (a) and Fig. 5 (b); Can know by Fig. 5; When the gaze duration long enough, for desirable white Gaussian noise, the propagation of electromagnetic field can not influence the time and space incoherence of signal.
The abbreviation of (1) by formula can be by the expression formula of computes received field
In this instance the expression formula of received field having been done following discretize handles.At first imaging region is divided into m * n=P (40 * 40) part, is about to it and is split into 1600 discrete subregions.Backscattering coefficient in each subregion is identical, is the distinguishable unit of minimum.
Setting the SF that receives signal time point is 1GHz; Gaze duration is sampled as Q; Because the signal of each emission all is desirable white Gaussian noise; Therefore on time dimension, the radiation field that front and back two sub-samplings obtain on imaging region is incoherent, can obtain the matrix form that gaze duration is sampled as Q to be:
According to the length of gaze duration, can obtain the different space-time bidimensionals statistical dependence characteristic of radiation field at random below, adopt the picture of multiple tupe inverting target area object again.
When space-time bidimensional when the statistical property single order of radiation field (being field intensity) is irrelevant at random, to scattered field and space-time bidimensional at random radiation field carry out the field intensity relevant treatment, its imaging inversion result when Fig. 7 has demonstrated the gaze duration long enough.
When space-time bidimensional when the statistical property of radiation field satisfies the single order part correlation at random, can find the solution through the method for directly inverting, can make the target information that originally is coupled separate.Fig. 8 has provided gaze duration and has been sampled as 2000, and signal to noise ratio (S/N ratio) is the inversion result of imaging region under the situation of 32dB.
When the solving method of direct inverse can not well be finally inversed by target image, the gaze duration that provided shown in Figure 9 was sampled as 500, and signal to noise ratio (S/N ratio) is to utilize the inversion chart picture of the target that the method for direct inverse obtains under the 32dB situation.Can find out that caused imaging results to have serious distortion owing to receive the deficiency of signal sampling this moment.
Thereby in this case, preferably scattered field and radiation field are handled then and can well be recovered target image with the method for separating optimization problem.Figure 10 has provided gaze duration and has been sampled as 500, signal to noise ratio (S/N ratio) be under the 32dB situation based on separating the corresponding imaging results of optimization problem method, can see that though there are several faults among the result, the global feature of target is but obviously visible.
The real aperture technique of tradition is in order to eliminate the influence of the part that transmits in the scatter echo; Just directly scattered field is handled; The methods such as matched filtering that adopt, the association process technology among the present invention are then united scattered field and space-time bidimensional radiation field at random, through the random distribution nature on the radiation field space at random; Then change through this space distribution diversity in time, realized that traditional microwave is stared the spatial discrimination in the radar beam under the same steric configuration.Say from imaging results; The present invention is through to having contained target with space-time bidimensional the be directly proportional scattered field of the cognizable space distribution information that changes of radiation field abundant information degree at random in the wave beam; With known space-time bidimensional at random radiation field carry out association process and obtain high resolving power target imaging information, realized the imaging of the high resolution target in the same wave beam of traditional real aperture imaging under the same steric configuration condition.
And the present invention uses foregoing formation method, adopts the form of numerical simulation, has verified correctness of the present invention, realizability and reliability.
Although illustrated and described embodiments of the invention; For those of ordinary skill in the art; Be appreciated that under the situation that does not break away from principle of the present invention and spirit and can carry out multiple variation, modification, replacement and modification that scope of the present invention is accompanying claims and be equal to and limit to these embodiment.

Claims (5)

1. microwave staring imaging correlating method is characterized in that said correlating method may further comprise the steps:
Through radiation field at random the target area object is shone, on object plane, form radiation field distribution;
Carry out scattered field through hyperchannel or single channel and receive, and the imaging physical process of target area object is carried out electromagnetic field analysis, to obtain the space-time bidimensional radiation field function in the corresponding whole imaging region scope of object plane;
Add up the space-time bidimensional correlation properties of radiation field at random according to said space-time bidimensional radiation field function; And
According to said statistical dependence characteristic, to scattered field and at random radiation field carry out association process obtains the target area object with inverting imaging.
2. microwave staring imaging correlating method as claimed in claim 1 is characterized in that, the space-time bidimensional at random the statistical dependence characteristic of radiation field satisfy the single order field intensity when irrelevant, to scattered field and at random radiation field carry out the field intensity relevant treatment.
3. microwave staring imaging correlating method as claimed in claim 1 is characterized in that, the space-time bidimensional at random the statistical dependence characteristic of radiation field satisfy high-order when irrelevant, to scattered field and at random radiation field carry out the relevant treatment of high-order field amount.
4. microwave staring imaging correlating method as claimed in claim 1; It is characterized in that; The space-time bidimensional at random the statistical dependence characteristic of radiation field satisfy single order or high order part when relevant, through the mode that radiation field is at random directly inverted scattered field is handled.
5. microwave staring imaging correlating method as claimed in claim 1; It is characterized in that; The space-time bidimensional at random the statistical dependence characteristic of radiation field satisfy the relevant and gaze duration of single order or high order part when too short, utilize the mode of separating optimum solution to scattered field and at random radiation field handle.
CN201110000698A 2011-01-04 2011-01-04 Microwave staring imaging correlation method Expired - Fee Related CN102141617B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110000698A CN102141617B (en) 2011-01-04 2011-01-04 Microwave staring imaging correlation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110000698A CN102141617B (en) 2011-01-04 2011-01-04 Microwave staring imaging correlation method

Publications (2)

Publication Number Publication Date
CN102141617A CN102141617A (en) 2011-08-03
CN102141617B true CN102141617B (en) 2012-08-29

Family

ID=44409275

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110000698A Expired - Fee Related CN102141617B (en) 2011-01-04 2011-01-04 Microwave staring imaging correlation method

Country Status (1)

Country Link
CN (1) CN102141617B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103513246B (en) * 2013-10-09 2015-11-04 中国科学院空间科学与应用研究中心 A kind of system and method for sub-wavelength imaging
CN103837873B (en) * 2014-03-14 2016-08-17 中国科学技术大学 A kind of microwave based on floating platform closely spaced array antenna stares relevance imaging system and formation method
CN104199028B (en) * 2014-09-03 2016-08-24 西安电子科技大学 The microwave relevance imaging method of emission array is rotated based on radar
CN104597444B (en) * 2015-01-30 2017-05-03 中国科学技术大学 Microwave gaze high resolution imaging method based on intensity association
CN104569975B (en) * 2015-02-09 2017-04-19 中国科学技术大学 Far-field measurement method of time-space random radiation field in microwave staring correlated imaging system
CN104714229A (en) * 2015-02-09 2015-06-17 中国科学技术大学 Microwave gazing correlated imaging treatment method convenient in extracting of object contour
CN104931963B (en) * 2015-06-08 2017-08-29 中国科学技术大学 A kind of method that moving-target microwave stares relevance imaging
CN107678028B (en) * 2017-11-22 2020-05-12 中国科学技术大学 Microwave staring correlated imaging method under low signal-to-noise ratio condition
CN108828593B (en) * 2018-04-26 2020-06-02 电子科技大学 Random radiation radar correlation imaging method
CN111257871B (en) * 2020-03-09 2023-06-16 中国科学技术大学 Single-antenna radiation source design method for microwave staring correlated imaging

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6996292B1 (en) * 2002-04-18 2006-02-07 Sandia Corporation Staring 2-D hadamard transform spectral imager
CN101776749A (en) * 2010-01-14 2010-07-14 北京理工大学 Vehicle-mounted staring imaging laser radar synchronous triggering system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6996292B1 (en) * 2002-04-18 2006-02-07 Sandia Corporation Staring 2-D hadamard transform spectral imager
CN101776749A (en) * 2010-01-14 2010-07-14 北京理工大学 Vehicle-mounted staring imaging laser radar synchronous triggering system

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Hai Deng等.Radar Target Imaging Using Distributed Radar Networks.《IEEE》.2009,全文. *
Lan zhao等.Study on Mechanism of Radiometer Staring Imaging.《APMC2005 Proceedings》.2005,全文. *
刘兆军等.面阵凝视型成像空间应用技术.《红外与激光工程》.2006,第35卷(第5期),541-545. *
韩兴斌等.分布式多通道雷达成像技术.《电子与信息学报》.2007,第29卷(第10期),2354-2358. *

Also Published As

Publication number Publication date
CN102141617A (en) 2011-08-03

Similar Documents

Publication Publication Date Title
CN102141617B (en) Microwave staring imaging correlation method
KR102572612B1 (en) Method and system for radar pulse collision avoidance
KR102196734B1 (en) Apparatus and Method for Jamming in Synthetic Aperture Radar
Colone et al. Polarimetric passive coherent location
JP2015179089A (en) Submillimeter radar using phase information
EP1851568B1 (en) Method for detecting in bistatic mode by passively operating on non-co-operating radioelectric transmissions
US10634790B2 (en) Digital passband processing of wideband modulated optical signals for enhanced imaging
Lu et al. An improved sliding matched filter method for interrupted sampling repeater jamming suppression based on jamming reconstruction
Colone et al. Non‐coherent adaptive detection in passive radar exploiting polarimetric and frequency diversity
KR101426226B1 (en) Signal processing method of radar
JP2023536150A (en) Multi-frequency radar array system and sensor fusion for viewing around corners during autonomous driving
JP2009236720A (en) Moving target detector
Zhu et al. Measurement and ray-tracing simulation for millimeter-wave automotive radar
Giusti et al. Electronic countermeasure for OFDM‐based imaging passive radars
O'connor et al. Waveform design considerations for modulated pulse lidar
Bączyk et al. The impact of reference channel SNR on targets detection by passive radars using DVB-T signals
Kim et al. GPR image enhancement based on frequency shifting and histogram dissimilarity
Kirk et al. Mitigation of clutter modulation in cognitive radar for spectrum sharing applications
Ma et al. Microwave staring correlated imaging and resolution analysis
Guo et al. Anti‐interrupted sampling repeater jamming method for random pulse repetition interval and intra‐pulse frequency agile radar
RU2472176C1 (en) Method for passive detection of air objects
KR20200054657A (en) Automotive pulse radar device
Colone et al. Exploiting polarimetric diversity in FM-based PCL
Abdullah et al. Target prediction in forward scattering radar
Immoreev et al. Ultrawideband radar special features & terminology

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: 20120829

Termination date: 20210104

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