[go: up one dir, main page]

CN101510012A - Composite phase mask plate - Google Patents

Composite phase mask plate Download PDF

Info

Publication number
CN101510012A
CN101510012A CNA2009100971796A CN200910097179A CN101510012A CN 101510012 A CN101510012 A CN 101510012A CN A2009100971796 A CNA2009100971796 A CN A2009100971796A CN 200910097179 A CN200910097179 A CN 200910097179A CN 101510012 A CN101510012 A CN 101510012A
Authority
CN
China
Prior art keywords
function
sin
phase mask
sinusoidal
question mark
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
Application number
CNA2009100971796A
Other languages
Chinese (zh)
Other versions
CN101510012B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN2009100971796A priority Critical patent/CN101510012B/en
Publication of CN101510012A publication Critical patent/CN101510012A/en
Application granted granted Critical
Publication of CN101510012B publication Critical patent/CN101510012B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Preparing Plates And Mask In Photomechanical Process (AREA)

Abstract

The invention discloses a composite phase mask, the phase distribution function of which contains two different functions, namely, a sinusoidal type function and a third power type function. The function expression of the phase distribution function Theta(x, y) is as follows: Theta(x, y) is equal to Alpha question mark (x<3> plus y<3>) plus Beta question mark (sin(Omega question mark x) plus sin(Omega question mark y)), and in the expression, Beta question mark (sin(Omega question mark x) plus sin(Omega question mark y)) is a sinusoidal type function, and Alpha question mark (x<3> plus y<3>) is a third power type function; wherein Alpha is the amplitude of the third power type function, Beta is the amplitude of the sinusoidal type function, Omega expresses the angular frequency of the sinusoidal type function, and x, y is space coordinate with normalized aperture plane. The invention also discloses an imaging system which adopts the composite phase mask compounded by the sinusoidal type function and the third power type function and has smaller numerical value of the Fisher information, and the phase mask is not sensitive to defocusing, thus having better expansion capability of depth of focus.

Description

A kind of composite phase mask plate
Technical field
The present invention relates to a kind of phase mask plate, relate in particular to two kinds of compound composite phase mask plates of function of a kind of cube type and sinusoidal pattern
Background technology
As everyone knows, any optical system all has certain depth of field, has only the object within the field depth can be by blur-free imaging, and the object that exceeds field depth will thicken through the formed image of optical system, thereby causes losing of details, texture.
As far back as nineteen sixty, W.T.Welford has just inquired into the use ring aperture, and promptly the central obscuration mode is come the possibility of the expansion optical system depth of field.1987, Ting-Chung Poon and MasoudMotamedi are on the basis of using ring aperture, rely on the digital image restoration algorithm to verify depth of field expansion effect practically, this following period of time from the end of the eighties to the beginning of the nineties, people such as Ojeda-Castaneda Jorge use classical apodization also system's this problem of depth of field continuation to be launched extensive and deep research.Studies show that, change the good result that aperture shape can be created the big depth of field, but this effect is with the logical light quantity of sacrificial system with to reduce resolution be that cost is reached.In addition, from the mid-80 so far, also have many non-mainstream methods to be suggested the depth of field that is used for the expansion optical system, use Fresnel lens, Zeev Zalevsky and Shai Ben-Yaish to use birefringece crystal as people such as Guy Indebetouw.Light, e light characteristic have obtained the effect of depth of field continuation too.
Though said method has all played positive effect to the depth of field of expansion optical system, but effect all is more limited comparatively speaking, real breakthrough starts from nineteen ninety-five, the E.R.Dowski of Univ Colorado-Boulder USA is in traditional imaging system, the ambiguity function that uses in the field of radar and out of focus optical transfer function combined proposed wavefront coded imaging technique, and derived classical cube PHASE DISTRIBUTION function mask plate, after this researcher has proposed logarithmic, exponential type, polynomial type PHASE DISTRIBUTION function mask plate again successively.
Can determine the physical arrangement of mask plate according to the PHASE DISTRIBUTION function, as in Chinese invention patent 200410018159.2, disclosing the method for a kind of discrete continuous mixed-phase type phase-plate and realization super-resolution thereof, the phase function of this phase-plate comprises two parts: a part is the continuous distribution phase-plate, and another part is the stepped ramp type phase-plate.For example the PhD dissertation of the Wanli Chi of U.S. Rochester optical centre " Computational Imaging System for Extended Depth of Field " chapter 4 " Lens Design and Fabrication " has promptly been introduced the method for determining mask plate face shape according to the PHASE DISTRIBUTION function all-sidely again.
Wavefront coded imaging technique can enlarge the depth of field of optical system effectively by image restoration is combined with optical design under the prerequisite of inexcessive sacrificial light systemic resolution and logical light quantity.Researcher has been done a large amount of theoretical analysises and experimental study, has affirmed the powerful power that wavefront coded imaging technique had, and has proved and verified that the depth of field improves 10 times effect.This light/number hybrid imaging technology has been broken the intrinsic restriction of conventional optical systems when extended depth-of-field, is a kind of imaging technique that practical application foreground is arranged very much, is subjected to paying close attention to widely.
Summary of the invention
The invention provides a kind of composite phase mask plate with better depth of field scalability.
A kind of composite phase mask plate, the PHASE DISTRIBUTION function of described composite phase mask plate comprise sinusoidal pattern and two kinds of functions of cube type, the PHASE DISTRIBUTION function # (x, functional expression y) is as follows:
θ(x,y)=α·(x 3+y 3)+β·(sin(ω·x)+sin(ω·y))
In the formula, β (sin (ω x)+sin (ω y)) is the sinusoidal pattern function, α (x 3+ y 3) be the cube type function; Wherein, α is the amplitude of cube type function, and β is the amplitude of sinusoidal pattern function, and ω represents the angular frequency of sinusoidal pattern function, and x, y are the normalized volume coordinate of aperture plane.
Wherein, as preferably, cube type function amplitude α value is between 90 to 300, and sinusoidal pattern function amplitude β value is between scope 50 to 260.The angular frequency value of sinusoidal pattern function in scope at 0.3 π between 0.6 π because not obvious when the improvement in performance of the too small composite phase plate of ω value, the excessive phase-plate that then causes of ω value is blocked up, has increased difficulty of processing.
In the phase-plate of various forms PHASE DISTRIBUTION function, cube type phase-plate is classical type, can realize bigger depth of field expansion.Superpose on the basis of the cube type function composite phase mask plate plate of sinusoidal component of the present invention has bigger depth of field extended capability.
A kind of imaging system is made up of imaging lens, composite phase mask plate, image detector, graphics processing unit.
The PHASE DISTRIBUTION function of described composite phase mask plate comprises sinusoidal pattern and two kinds of functions of cube type, and functional expression is as follows:
θ(x,y)=α·(x 3+y 3)+β·(sin(ω·x)+sin(ω·y))
In the formula, α is the amplitude of cube type function, and β is the amplitude of sinusoidal pattern function, and ω represents the angular frequency of sinusoidal pattern function, and x, y are the normalized volume coordinate of aperture plane.
In imaging system, imageable target by imaging lens and composite phase mask plate after, on image detector, form fuzzy intermediary image, graphics processing unit carries out deconvolution processing to intermediary image then, finally obtains focusing on distinct image.
Imaging lens can adopt existing various imaging lens.Image detector can adopt CCD etc.
A kind of imaging system of the present invention has increased a composite phase mask plate in traditional optical imaging system, be placed on the pupil.After light path is passed through this composite phase mask plate, optical information is encoded, be that optical system its modulating transfer function value in bigger out of focus scope changes very little (promptly insensitive to out of focus), the transport function of simultaneity factor does not have zero point or nearly zero point in passband.This means that when there was out of focus in system, the information that exceeds the primal system field depth was not lost, just be encoded in a kind of known mode.Being imaged on the CCD is a width of cloth blurred picture, with digital information processing system to its decoding of deconvoluting, thereby recover sharp keen distinct image.
In wavefront coded field, Fisher information is one of important indicator of estimating depth of field extended capability.According to the definition of Fisher information, in the ideal case, no matter the out of focus parameter is got any value, and Fisher information all should equal zero.Therefore, the numerical value of Fisher information is more little, just means that phase-plate is insensitive more to out of focus, and depth of field extended capability is good more.
The present invention adopts the Fisher information curve of imaging system of the compound composite phase mask plate of sinusoidal pattern and two kinds of functions of cube type little, has depth of field extended capability preferably.
Description of drawings
Fig. 1 is the structural representation of a kind of imaging system of the present invention;
Fig. 2 is the comparison synoptic diagram of Fisher information curve and the Fisher information curve of the imaging system that adopts cube type phase mask plate that adopts the imaging system of composite phase mask plate of the present invention.
Embodiment
A kind of composite phase mask plate, the PHASE DISTRIBUTION function of composite phase mask plate comprise sinusoidal pattern and two kinds of functions of cube type, the PHASE DISTRIBUTION function # (x, functional expression y) is as follows:
θ(x,y)=α·(x 3+y 3)+β·(sin(ω·x)+sin(ω·y))
In the formula, β (sin (ω x)+sin (ω y)) is the sinusoidal pattern function, α (x 3+ y 3) be the cube type function; Wherein, α is the amplitude of cube type function, and β is the amplitude of sinusoidal pattern function, and ω represents the angular frequency of sinusoidal pattern function, and x, y are the normalized volume coordinate of aperture plane.
The face shape of a kind of composite phase mask plate of the present invention and the face shape of exponential type phase mask plate are relatively referring to accompanying drawing 2, the face shape of a kind of composite phase mask plate of the present invention dots in accompanying drawing 2, and the face shape of exponential type phase mask plate is represented with solid line in accompanying drawing 2.
Referring to Fig. 1, a kind of imaging system of the present invention forms 5 by imaging lens 2, composite phase mask plate 3, image detector 4, graphics processing unit.
After imageable target 1 is passed through imaging lens 2 and composite phase mask plate 3, form fuzzy intermediary image on image detector 4, graphics processing unit 5 carries out deconvolution processing then, finally obtains focusing on distinct image 6.
The PHASE DISTRIBUTION function of composite phase mask plate comprises sinusoidal pattern and two kinds of functions of cube type, the PHASE DISTRIBUTION function # (x, functional expression y) is as follows:
θ(x,y)=α·(x 3+y 3)+β·(sin(ω·x)+sin(ω·y))
In the formula, β (sin (ω x)+sin (ω y)) is the sinusoidal pattern function, α (x 3+ y 3) be the cube type function; Wherein, α is the amplitude of cube type function, and β is the amplitude of sinusoidal pattern function, and ω represents the angular frequency of sinusoidal pattern function, and x, y are the normalized volume coordinate of aperture plane.
In wavefront coded field, Fisher information is one of important indicator of estimating depth of field extended capability.According to the definition of Fisher information, in the ideal case, no matter the out of focus parameter is got any value, and Fisher information all should equal zero.Therefore, the numerical value of Fisher information is more little, just means that phase-plate is insensitive more to out of focus, and depth of field extended capability is good more.
The present invention adopts the Fisher information curve ratio of Fisher information curve and the imaging system that adopts cube type phase mask plate of imaging system of composite phase mask plate referring to accompanying drawing 2.
During data after α, β, each parameter of ω adopt in the table 1 nine groups to optimize respectively, the Fisher information curve of imaging system that adopts cube type phase mask plate is shown in the solid line in the accompanying drawing 2, and the Fisher information curve of imaging system that adopts the compound composite phase mask plate of two kinds of functions of sinusoidal pattern of the present invention and cube type is shown in the dotted line in the accompanying drawing 2.
Tables of data after table 1 α, β, nine groups of optimizations of each parameter of ω
Figure A200910097179D00071
By contrast adopt the PHASE DISTRIBUTION function be cube type mask plate imaging system Fisher information curve and adopt the Fisher information curve of PHASE DISTRIBUTION function for the imaging system of the compound composite phase mask plate of sinusoidal pattern of the present invention and two kinds of functions of cube type, the Fisher information curve of imaging system of the present invention is starkly lower than the former as can be seen, therefore has better depth of field extended capability.

Claims (2)

1、一种复合型相位掩模板,其特征在于:所述的复合型相位掩模板的相位分布函数包含正弦型和三次方型两种函数,相位分布函数θ(x,y)函数式如下:1. A composite phase mask, characterized in that: the phase distribution function of the composite phase mask includes two functions of sinusoidal and cubic, and the phase distribution function θ (x, y) is as follows: θ(x,y)=α·(x3+y3)+β·(sin(ω·x)+sin(ω·y))θ(x,y)=α·(x 3 +y 3 )+β·(sin(ω·x)+sin(ω·y)) 式中,β·(sin(ω·x)+sin(ω·y))为正弦型函数,α·(x3+y3)为三次方型函数;其中,α为三次方型函数的幅值,β为正弦型函数的幅值,ω表示正弦型函数的角频率,x,y为孔径平面归一化的空间坐标。In the formula, β·(sin(ω·x)+sin(ω·y)) is a sinusoidal function, α·(x 3 +y 3 ) is a cubic function; where α is the magnitude of the cubic function value, β is the amplitude of the sinusoidal function, ω is the angular frequency of the sinusoidal function, and x, y are the normalized space coordinates of the aperture plane. 2、一种成像系统,由成像镜头、复合型相位掩模板、图像探测器、图像处理单元组成,其特征在于:所述的复合型相位掩模板的相位分布函数包含正弦型和三次方型两种函数,函数式如下:2. An imaging system is composed of an imaging lens, a composite phase mask, an image detector, and an image processing unit. It is characterized in that: the phase distribution function of the composite phase mask includes both sinusoidal and cubic. A function, the function formula is as follows: θ(x,y)=α·(x3+y3)+β·(sin(ω·x)+sin(ω·y))θ(x,y)=α·(x 3 +y 3 )+β·(sin(ω·x)+sin(ω·y)) 式中,α为三次方型函数的幅值,β为正弦型函数的幅值,ω表示正弦型函数的角频率,x,y为孔径平面归一化的空间坐标;In the formula, α is the amplitude of the cubic function, β is the amplitude of the sinusoidal function, ω is the angular frequency of the sinusoidal function, and x, y are the space coordinates normalized on the aperture plane; 在成像系统中,成像目标通过成像镜头和复合型相位掩模板后,在图像探测器上形成模糊的中间像,然后图像处理单元对中间像进行去卷积处理,最终得到聚焦清晰的图像。In the imaging system, after the imaging target passes through the imaging lens and the composite phase mask, a blurred intermediate image is formed on the image detector, and then the image processing unit deconvolutes the intermediate image to finally obtain a sharply focused image.
CN2009100971796A 2009-03-26 2009-03-26 Composite phase mask plate Expired - Fee Related CN101510012B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100971796A CN101510012B (en) 2009-03-26 2009-03-26 Composite phase mask plate

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100971796A CN101510012B (en) 2009-03-26 2009-03-26 Composite phase mask plate

Publications (2)

Publication Number Publication Date
CN101510012A true CN101510012A (en) 2009-08-19
CN101510012B CN101510012B (en) 2010-08-11

Family

ID=41002457

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100971796A Expired - Fee Related CN101510012B (en) 2009-03-26 2009-03-26 Composite phase mask plate

Country Status (1)

Country Link
CN (1) CN101510012B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011099003A1 (en) * 2010-02-09 2011-08-18 Xceed Imaging Ltd. Imaging with extended depth of focus for use with polycromatic light
US8192022B2 (en) 2004-08-16 2012-06-05 Xceed Imaging Ltd. Optical method and system for extended depth of focus
CN103729824A (en) * 2013-12-17 2014-04-16 北京智谷睿拓技术服务有限公司 Information interactive method and information interactive system
WO2015090038A1 (en) * 2013-12-17 2015-06-25 Beijing Zhigu Rui Tuo Tech Co., Ltd Information exchange
US10795168B2 (en) 2017-08-31 2020-10-06 Metalenz, Inc. Transmissive metasurface lens integration
CN113031260A (en) * 2021-03-12 2021-06-25 中国人民解放军国防科技大学 Wavefront coding imaging system based on arcsine type phase mask plate
CN113614602A (en) * 2019-03-28 2021-11-05 索尼集团公司 Optical system, endoscope, and medical image processing system
WO2023016149A1 (en) * 2021-08-13 2023-02-16 Oppo广东移动通信有限公司 Optical lens, image acquisition device, and electronic device
US11906698B2 (en) 2017-05-24 2024-02-20 The Trustees Of Columbia University In The City Of New York Broadband achromatic flat optical components by dispersion-engineered dielectric metasurfaces
US11927769B2 (en) 2022-03-31 2024-03-12 Metalenz, Inc. Polarization sorting metasurface microlens array device
US11978752B2 (en) 2019-07-26 2024-05-07 Metalenz, Inc. Aperture-metasurface and hybrid refractive-metasurface imaging systems

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA3089373A1 (en) 2018-01-24 2019-08-01 President And Fellows Of Harvard College Polarization state generation with a metasurface
EP3799626A4 (en) 2018-07-02 2022-03-30 Metalenz, Inc. META SURFACES FOR LASER SPECKLE REDUCTION

Cited By (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8192022B2 (en) 2004-08-16 2012-06-05 Xceed Imaging Ltd. Optical method and system for extended depth of focus
US11802997B2 (en) 2010-02-09 2023-10-31 Brien Holden Vision Institute Limited Optical apparatus with structure for liquid invariant performance
US11802998B2 (en) 2010-02-09 2023-10-31 Brien Holden Vision Institute Limited Imaging system with optimized extended depth of focus
US8531783B2 (en) 2010-02-09 2013-09-10 Xceed Imaging Ltd. Imaging method and system for imaging with extended depth of focus
US12313916B2 (en) 2010-02-09 2025-05-27 Brien Holden Vision Institute Limited Optical apparatus with structure for liquid invariant performance
US8913331B2 (en) 2010-02-09 2014-12-16 Brien Holden Vision Institute Imaging method and system with optimized extended depth of focus
US8955968B2 (en) 2010-02-09 2015-02-17 Brien Holden Vision Institute Imaging with extended depth of focus for use with polychromatic light
US12282175B2 (en) 2010-02-09 2025-04-22 Brien Holden Vision Institute Limited Multi-focal lens
US9134543B2 (en) 2010-02-09 2015-09-15 Brien Holden Vision Institute Imaging system with optimized extended depth of focus
US12259565B2 (en) 2010-02-09 2025-03-25 Brien Holden Vision Institute Limited Imaging system with optimized extended depth of focus
US9429768B2 (en) 2010-02-09 2016-08-30 Brien Holden Vision Institute Imaging method and system with optimized extended depth of focus
US9500875B2 (en) 2010-02-09 2016-11-22 Brien Holden Vision Institute Imaging with extended depth of focus for use with polycromatic light
WO2011099003A1 (en) * 2010-02-09 2011-08-18 Xceed Imaging Ltd. Imaging with extended depth of focus for use with polycromatic light
US10031334B2 (en) 2010-02-09 2018-07-24 Brien Holden Vision Institute Optical apparatus with structure for liquid invariant performance
US10078159B2 (en) 2010-02-09 2018-09-18 Brien Holden Vision Institute Multi-focal lens
US8169716B2 (en) 2010-02-09 2012-05-01 Xceed Imaging, Ltd. Optical apparatus with structure for liquid invariant performance
US11079517B2 (en) 2010-02-09 2021-08-03 Brien Holden Vision Institute Limited Optical apparatus with structure for liquid invariant performance
US9239471B2 (en) 2010-02-09 2016-01-19 Brien Holden Vision Institute Multi-focal lens
US10175392B2 (en) 2010-02-09 2019-01-08 Brien Holden Vision Institute Imaging system with optimized extended depth of focus
US11199651B2 (en) 2010-02-09 2021-12-14 Brien Holden Vision Institute Limited Imaging system with optimized extended depth of focus
US9836857B2 (en) 2013-12-17 2017-12-05 Beijing Zhigu Rui Tuo Tech Co., Ltd. System, device, and method for information exchange
CN103729824A (en) * 2013-12-17 2014-04-16 北京智谷睿拓技术服务有限公司 Information interactive method and information interactive system
WO2015090038A1 (en) * 2013-12-17 2015-06-25 Beijing Zhigu Rui Tuo Tech Co., Ltd Information exchange
US11906698B2 (en) 2017-05-24 2024-02-20 The Trustees Of Columbia University In The City Of New York Broadband achromatic flat optical components by dispersion-engineered dielectric metasurfaces
US10795168B2 (en) 2017-08-31 2020-10-06 Metalenz, Inc. Transmissive metasurface lens integration
US11988844B2 (en) 2017-08-31 2024-05-21 Metalenz, Inc. Transmissive metasurface lens integration
US11579456B2 (en) 2017-08-31 2023-02-14 Metalenz, Inc. Transmissive metasurface lens integration
CN113614602A (en) * 2019-03-28 2021-11-05 索尼集团公司 Optical system, endoscope, and medical image processing system
US12114833B2 (en) 2019-03-28 2024-10-15 Sony Group Corporation Optical system, endoscope, and medical image processing system
US11978752B2 (en) 2019-07-26 2024-05-07 Metalenz, Inc. Aperture-metasurface and hybrid refractive-metasurface imaging systems
CN113031260A (en) * 2021-03-12 2021-06-25 中国人民解放军国防科技大学 Wavefront coding imaging system based on arcsine type phase mask plate
WO2023016149A1 (en) * 2021-08-13 2023-02-16 Oppo广东移动通信有限公司 Optical lens, image acquisition device, and electronic device
US11927769B2 (en) 2022-03-31 2024-03-12 Metalenz, Inc. Polarization sorting metasurface microlens array device

Also Published As

Publication number Publication date
CN101510012B (en) 2010-08-11

Similar Documents

Publication Publication Date Title
CN101510012B (en) Composite phase mask plate
CN101510011A (en) Composite phase mask plate
CN110392266B (en) Light field video coding method based on pseudo video sequence, terminal equipment and storage medium
CN100476504C (en) A phase mask and an imaging system using the phase mask
CN100520496C (en) Multi-matrix depth of field wavefront-coding system and method
CN102129096B (en) Phase plate and field depth expanding imaging system applying same
WO2022199471A1 (en) Quick response code generation method and apparatus, quick response code decoding method and apparatus, and device
CN103247020A (en) Fisheye image spread method based on radial characteristics
JP2017519250A (en) Imaging system, method and application
JP6033673B2 (en) Imaging device
JP6151641B2 (en) Imaging apparatus, imaging system, imaging method, and image processing method
CN103901615B (en) Little recessed imaging optical system
Chen et al. Light field based digital refocusing using a DSLR camera with a pinhole array mask
Tang et al. Utilizing optical aberrations for extended-depth-of-field panoramas
CN104184936A (en) Image focusing processing method and system based on light field camera
Goldlücke et al. Plenoptic Cameras.
CN103885176B (en) Phase mask plate and wavefront coding imaging system capable of adjusting quality of intermediate coding image
CN116088173B (en) Scattering imaging system, optical axis non-invasive searching method and device thereof, and storage medium
CN203882012U (en) Phase mask and wavefront coding system
CN105807550A (en) Inversion type super-speed image forming method
CN106464808A (en) Image processing apparatus, image pickup apparatus, image processing method, image processing program, and storage medium
Popovic et al. State-of-the-art multi-camera systems
Gao et al. Wave-front coding technology to extend depth of field in remote sensing optical system
Silveira Pupil Function Engineered for Improved SLAM Feature Localization
Horisaki et al. High-quality image acquisition using a compact compound-eye camera

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

Termination date: 20120326