CN101660992B - A method for rapidly detecting the sedimentation velocity of algae cells - Google Patents

A method for rapidly detecting the sedimentation velocity of algae cells Download PDF

Info

Publication number
CN101660992B
CN101660992B CN2009100931394A CN200910093139A CN101660992B CN 101660992 B CN101660992 B CN 101660992B CN 2009100931394 A CN2009100931394 A CN 2009100931394A CN 200910093139 A CN200910093139 A CN 200910093139A CN 101660992 B CN101660992 B CN 101660992B
Authority
CN
China
Prior art keywords
sampling volume
settling
virtual sampling
cells
settling velocity
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
CN2009100931394A
Other languages
Chinese (zh)
Other versions
CN101660992A (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN2009100931394A priority Critical patent/CN101660992B/en
Publication of CN101660992A publication Critical patent/CN101660992A/en
Application granted granted Critical
Publication of CN101660992B publication Critical patent/CN101660992B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analysing Biological Materials (AREA)

Abstract

The invention relates to a method for quickly detecting settling velocity of alga cells. The method includes the following steps: adopting volumetric solution with known concentration and corresponding particle number to calibrate the virtual sampling volume and virtual sampling volume height corresponding to a microscope objective, specifically, taking alga cell solution in exponential phase to conduct centrifugal weight suspension; continuously diluting the alga cell solution obtained in step b at equal multiplying power to obtain a plurality of samples; using an inverted microscope to conduct imaging observation on the bottom of all holes bearing the alga cell solution and selecting a settling column for detecting the settling velocity; continuously imaging the bottom of the settling column at an interval of equal time to obtain a plurality of images; adopting a digital image analysis method to calculate and obtain the number of alga cells in each image; and adopting a linear regression method to analyze the corresponding relation between the number of alga cells and settling time and using the calibrated virtual sampling volume height to calculate the settling velocity of alga cells. The method has the advantages of simple method and device, economical and practical properties, high detection accuracy, simple and convenient operation, etc.

Description

A kind of method of fast detecting settling velocity of alga cells
Technical field
The invention belongs to the detection technique field of phytoplankton settling velocity, especially relate to the method for settling velocity that a kind of inverted microscope that adopts imaging continuous time and digital image analysis detects the alga cells of 5-50um diameter.
Background technology
Settling velocity directly reflects physiological status and the nourishment of phytoplankton, so phytoplankton cell settlement speed is significant on phytoplankton ecology.Yet because the frustule size is less, settling velocity is slow, there is very big technical difficulty in the settling velocity of directly measuring frustule.In the current measuring methods, simple and efficient method measuring accuracy is limited, and its checkout equipment complexity of method that precision is high, and product-free product can not satisfy the research work demand.
Existing frustule settling velocity detection method mainly adopts settling leg to detect, and is broadly divided into two class methods according to the difference of detection position, is respectively that the settling leg fixed position is detected and the settling leg entire scan.The fixed position is the bottom of settling leg normally, the inverted microscope method that typical method has Smayda and Boleyn to adopt, this method preparation concentration known algae liquid sedimentation in settling leg, use the quantity of inverted microscope naked eyes counting settling leg bottom frustule, with the cut off of 80% frustule sedimentation as sedimentation, this method has been ignored as yet not the sedimentation cell to the influence of bulk settling speed; People such as Eppley adopt the fluorescence signal of photofluorometer monitoring settling leg bottom to change the settling velocity that detects alga cells; Bienfang has invented ' SETCOL ' method, the settling leg that this method uses a top and bottom all to be provided with sample tap is carried out sedimentation experiment, after a set time, analyze the settling velocity that frustule is asked in the distribution of frustule concentration in settling leg by hand sampling.
The entire scan method is divided into continuous sweep and interrupted scanning two classes again.Walsby and Xypolyta have invented and have a kind ofly determined the method for time from the detection alga cells sedimentation of the diverse location sample analysis of settling leg; Walsby and Holland combining laser scanning and Percoll gradient solution have been invented a kind of frustule settling velocity detection method of continuous sweep.There is following weak point in these methods: 1) manually participate in morely, can not guarantee the accuracy and the consistance of testing result; 2) the equipment complexity of most methods uses, and criteria of right and wrong instrument is unfavorable for the laboratory study application.
Summary of the invention
For overcoming prior art in the weak point that exists aspect the detection frustule settling velocity, a kind of method of fast detecting settling velocity of alga cells is proposed, this method has advantages such as square law device is simple, economical and practical, accuracy of detection is high, easy and simple to handle based on imaging continuous time of inverted microscope and the detection method of digital image analysis.The common lab that disposes the inverted microscope of imaging device just can be carried out the settling velocity detection of frustule.
The method of a kind of fast detecting settling velocity of alga cells of the present invention is characterized in that, may further comprise the steps:
The volumetric solution of a, employing concentration known and corresponding amounts of particles is demarcated micro objective corresponding virtual sampling volume and virtual sampling volume height Δ H:
Wherein the floorage S of virtual sampling volume is:
Figure G2009100931394D00021
Virtual sampling volume height Δ H is: Δ H=K 1/ C 0S; K in the formula 1Corresponding to C in the definite virtual sampling volume of linear regression method 0The endocorpuscular quantity of calibration solution of concentration;
B, get the frustule solution that is in exponential phase, use low speed centrifuge to carry out centrifugal resuspended with centrifuge tube;
C, with etc. the frustule solution that obtains of multiplying power serial dilution b step, obtain to be not less than three samples, place continuous several holes of 24 orifice plates respectively;
D, with inverted microscope imaging is carried out in each hole that is loaded with frustule solution and observed, the imaging observation place is the bottom in hole, and selecting initial cell quantity in the objective angular field is that 50-200 hole is as the settling leg of detection settling velocity;
E, to settling leg with interval equal time, settling leg bottom continuous imaging more than 3 times, is obtained 3 width of cloth altogether with epigraph;
F, based on the image that e step obtains, adopt method of digital image analysis to calculate the quantity of frustule in every width of cloth image;
The quantity of g, the frustule that calculates based on f step adopts linear regression method to analyze frustule quantity and the corresponding relation of settling time, and with a settling velocity that the virtual sampling volume height Δ H substitution of demarcation can calculate frustule in the step is:
V a=K 2ΔH/b 2
K in the formula 2And b 2Be respectively to increase frustule quantity and time t=0 moment frustule quantity (being initial cell quantity) by the unit interval in the definite virtual sampling volume of linear regression method.
Characteristics of the present invention and effect:
1, compared with the prior art, this method has good exploitativeness and higher precision, and device therefor all has configuration in general laboratory, can obtain high-precision frustule settling velocity under simple condition.
2, this method is introduced the notion of virtual sampling volume, only needs to demarcate once for a kind of microcobjective, does not need repetition in follow-up test, effectively reduces the workload of measuring.
3, this method adopts method of digital image analysis that continuous imaging is carried out in the frustule sedimentation, adopts linear regression method analysis of cells subsidence rate, has improved the robotization of testing process and the stability of testing result.
Description of drawings
Fig. 1 is the structural representation of the employed device of the inventive method.
Fig. 2 is the synoptic diagram of virtual sampling volume and settling leg.
Fig. 3 is the distribution schematic diagram of volumetric solution sample in 24 orifice plates of six kinds of concentration in the inventive method.
The method of digital image analysis FB(flow block) that Fig. 4 adopts for the present invention.
Embodiment
The method of a kind of fast detecting settling velocity of alga cells that the present invention proposes reaches embodiment in conjunction with the accompanying drawings and is described in detail as follows:
The pick-up unit that the present invention adopts as shown in Figure 1, this pick-up unit comprises light source 1, by being successively set on the inverted microscope that the light source below is formed with the condenser on the optical axis 2, article carrying platform 4, micro objective 5; Also comprise CCD camera 7, be arranged on the reflective mirror 6 between micro objective 5 and the CCD camera 7, reflective mirror 6 is used to make CCD camera 7 can receive the emergent light of micro objective 5; And multi-purpose computer 9, the data transmission interface 8 of connection multi-purpose computer 9 and CCD camera 7.This device all adopts the conventional universal product to assemble.
Among the embodiment of this device: inverted microscope is standard configuration Olympus X71, and light source 1 adopts 100W transmission light beam, and micro objective 5 differs object lens for Olympus long reach achromatism, numerical aperture 0.4, and enlargement ratio is 20X; It is 1392x 1040 that CCD camera 7 adopts the CCD camera resolution of Qimaging, and pixel size is 6.45um; Data transmission interface 8 is the IEEE-1394 interface, and this interface is finished view data and the transmission of control command between multi-purpose computer 9 and CCD camera 7.
Adopt the method for the foregoing description device fast detecting settling velocity of alga cells may further comprise the steps:
A, demarcation micro objective corresponding virtual sampling volume, determine the floorage S and the height Δ H of virtual sampling volume: said virtual sampling volume, as shown in Figure 2, right cylinder is a settling leg among the figure, cubic space is virtual sampling volume, this cubical floorage S is a CCD camera apparent field, can directly calculate according to the correlation parameter of object lens and CCD camera; It is relevant that clear degree falls in the image that causes owing to out of focus in virtual sampling volume height Δ H and the microcobjective depth of field and the imaging process;
The inventive method specifically may further comprise the steps:
1) the floorage S of the virtual sampling volume of calculating:
Figure G2009100931394D00031
In the present embodiment: S=1392 * 1040 * 6.45 * 6.45/20 2=150568um 2(1);
2) adopting diameter is the polystyrene standard spheric grain configuration volumetric solution 2mL of 5-30um (for example 10um), the volumetric solution concentration C 010 5-10 7Granules of polystyrene/mL;
3) this volumetric solution is carried out 2 multiplying powers dilutions continuously, obtain and comprise original volumetric solution totally 6 samples, each sample concentration is C i=C 0/ 2 i, wherein i=0-5 is the mark mark of sample, get the 1mL sample respectively and place continuous 6 holes of 24 orifice plates with delegation, as shown in Figure 2, C 0Sample is at A-1 position, C 1Sample is at A-2 position, C 2Sample is at A-3 position, C 3Sample is at A-4 position, C 4Sample is at A-5 position, C 5Sample is in the A-6 position;
4) use inverted microscope that the orifice plate bottom granules of polystyrene that sample respectively is housed is carried out imaging, present embodiment needs solution in the hole is no less than 5 piping and druming before imaging, move horizontally micro objective the diverse location imaging of each bottom surface, hole three times, add up the quantity of granules of polystyrene in every width of cloth image, obtain N after averaging i, N iCorresponding C iThe quantity of the granules of polystyrene of concentration granules of polystyrene solution, wherein i=0-5;
5) calculate virtual sampling volume height Δ H:
The quantity N of granules of polystyrene in each hole in theory iHas following relation with i
N i = C 0 × ΔH × S × 1 2 i - - - ( 2 )
Wherein S is the floorage of above-mentioned virtual sampling volume; In practical operation, adopt the first-order linear homing method to 6 groups of N iWith 1/2 iCarry out regretional analysis, can obtain a linear relationship
N i = K 1 × 1 2 i + b 1 - - - ( 3 )
K wherein 1C in the corresponding virtual sampling volume 0The quantity of concentration granules of polystyrene, b 1Since the error that factors such as operation are introduced, K 1, b 1It is real number; Suitable operation in granules of polystyrene solution allocation process (being that each sample concentration is accurate in the step 3)), b 1With respect to N iBe to ignore event, therefore, can make
K 1=C 0×ΔH×S (4)
With S and C 0, K 1Go into this formula, can try to achieve Δ H; The relative position relation in Δ H and hole (being right cylinder among the figure) as shown in Figure 3, Δ H is the height of virtual sampling volume (among the figure for cuboid);
B, get the chlorella that is in exponential phase to be detected (this algae mature cell diameter is about 10um) frustule solution 1mL with the 1.5mL centrifuge tube, this solution is added resuspended liquid uses low speed centrifuge to carry out centrifugal resuspended 2 times, last resuspension is a 2mL frustule solution, centrifuge speed 3000rpm wherein, centrifugation time is 1min, and resuspended liquid adopts deionized water;
C, the frustule solution that obtains with 2 multiplying power serial dilution b steps obtain three samples of 1mL, and sample concentration is respectively C a, C a/ 2, C a/ 4; Respectively get every kind of frustule solution example and place continuous three holes of 24 orifice plates, with pipettor sample well is blown and beaten gently, the piping and druming number of times is no less than 5 times;
D, with inverted microscope imaging is carried out in each hole that is loaded with frustule solution and observed numeration, the imaging observation place is the bottom in hole, and selecting initial cell quantity in the objective angular field is that 50-200 hole is as the settling leg of detection settling velocity;
E, to settling leg to settling leg bottom imaging once with interval equal time (for example every 5min) with the CCD camera, continuous imaging 10 times obtains 10 width of cloth images altogether;
F, based on 10 width of cloth images that e step obtains, adopt method of digital image analysis to calculate the quantity of frustule in every width of cloth image; Method of digital image analysis flow process in the present embodiment as shown in Figure 4, comprise: at first carry out the medium filtering noise reduction process for the original image of gathering, image adopts OSTU (big Tianjin threshold method) method to extract the image segmentation threshold value behind the noise reduction, to cut apart the background that obtains adopts pixel to adopt 0 gray scale to fill, image after the filling adopts corrosion and expanding method to remove the noise spot that produces in the image segmentation, adopts the quantity of the method counting frustule of connected component labeling at last; The embodiment of this method is the known technology of this area;
The quantity of g, 10 width of cloth image frustules that calculate based on f step adopts linear regression method to analyze frustule quantity and the corresponding relation of settling time, and with a in the step the virtual sampling volume height Δ H substitution of demarcation can calculate the settling velocity of frustule;
Theoretical derivation is the result show, for the uniform frustule solution of initial concentration, frustule quantity and time have linear relationship in the microscopic fields of view, and this point has also obtained affirmation from the regretional analysis of experimental result; By mathematical derivation as can be known, the quantity N of frustule in the virtual sampling volume of t Δ H height during the time Δ H(t) available following formulate:
N ΔH ( t ) = N t S Al × ( ΔH + V a t ) - - - ( 5 )
N in the formula tBe the total quantity of frustule in the hole, A is the floorage in hole, and l is the liquid level of frustule solution in sedimentation hole, V aSettling velocity for frustule; The physical significance of following formula clearly, wherein
Figure G2009100931394D00052
Be the highly inherent t=0 of Δ H frustule quantity constantly, Be the frustule quantity that increases in the Δ H height in the time at t; By the first-order linear regretional analysis, can obtain N Δ H(t) has following relation with time t
N ΔH(t)=K 2t+b 2(6)
Order K 2 = N t S Al × V a , b 2 = N t S Al × ΔH , Have thus
K 2/b 2=V a/ΔH (7)
Just the Δ H substitution following formula that obtains of step a can obtain V a:
V a=K 2ΔH/b 2=K 1K 2/C 0Sb 2 (8)。

Claims (1)

1.一种快速检测藻类细胞沉降速度的方法,其特征在于,包括以下步骤:1. A method for rapidly detecting algae cell settling velocity, is characterized in that, comprises the following steps: a、采用已知C0浓度及相应颗粒数量的标定溶液标定显微镜物镜对应的虚拟采样容积和虚拟采样容积高度ΔH:a. Calibrate the virtual sampling volume corresponding to the microscope objective lens and the height ΔH of the virtual sampling volume by using the calibration solution with the known C0 concentration and the corresponding number of particles: 其中虚拟采样容积的底面积S为:The bottom area S of the virtual sampling volume is:
Figure FSB00000220873200011
Figure FSB00000220873200011
虚拟采样容积高度ΔH为:ΔH=K1/C0S;式中K1对应通过线性回归方法确定的虚拟采样容积内C0浓度的标定溶液内颗粒的数量;The virtual sampling volume height ΔH is: ΔH=K 1 /C 0 S; where K 1 corresponds to the number of particles in the calibration solution of the C 0 concentration in the virtual sampling volume determined by the linear regression method; b、用离心管取处于对数生长期的藻细胞溶液,使用低速离心机进行离心重悬;b. Use a centrifuge tube to take the algae cell solution in the logarithmic growth phase, and use a low-speed centrifuge to centrifuge and resuspend; c、以等倍率连续稀释b步骤得到的藻细胞溶液,获得不小于三个样品,分别置于24孔板中的连续数个孔中;c. Serially dilute the algae cell solution obtained in step b at an equal ratio to obtain no less than three samples, which are respectively placed in several consecutive wells in a 24-well plate; d、用倒置显微镜对载有藻细胞溶液的各孔进行成像观察,成像观察位置为孔的底部,选择物镜视场内初始细胞数量为50-200个的孔作为检测沉降速度的沉降柱;d. Use an inverted microscope to image and observe each hole containing the algae cell solution. The imaging observation position is the bottom of the hole, and the hole with an initial cell number of 50-200 in the field of view of the objective lens is selected as a settling column for detecting the settling velocity; e、对沉降柱以间隔相等时间,对沉降柱底部连续成像3次以上,共获取3幅以上图像;e. Continuously image the bottom of the sedimentation column for more than 3 times at equal intervals for the sedimentation column, and obtain more than 3 images in total; f、基于e步获取的图像,采用数字图像分析方法计算得到每幅图像中藻细胞的数量;f, based on the image obtained in step e, the number of algae cells in each image is calculated by using a digital image analysis method; g、基于f步计算的藻细胞的数量,采用线性回归方法分析藻细胞数量与沉降时间的对应关系,并将a步中标定的虚拟采样容积高度ΔH代入即可计算出藻细胞的沉降速度为:g. Based on the number of algal cells calculated in step f, the linear regression method is used to analyze the corresponding relationship between the number of algal cells and the settling time, and the virtual sampling volume height ΔH calibrated in step a can be substituted to calculate the settling velocity of the algal cells as : Va=K2ΔH/b2V a =K 2 ΔH/b 2 , 式中K2和b2分别是通过线性回归方法确定的虚拟采样容积内单位时间增加藻细胞数量和时间t=0时刻藻细胞数量(即初始细胞数量)。In the formula, K 2 and b 2 are the number of algae cells increased per unit time in the virtual sampling volume and the number of algae cells at time t=0 (ie, the initial cell number) determined by the linear regression method, respectively.
CN2009100931394A 2009-09-29 2009-09-29 A method for rapidly detecting the sedimentation velocity of algae cells Expired - Fee Related CN101660992B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100931394A CN101660992B (en) 2009-09-29 2009-09-29 A method for rapidly detecting the sedimentation velocity of algae cells

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100931394A CN101660992B (en) 2009-09-29 2009-09-29 A method for rapidly detecting the sedimentation velocity of algae cells

Publications (2)

Publication Number Publication Date
CN101660992A CN101660992A (en) 2010-03-03
CN101660992B true CN101660992B (en) 2011-01-12

Family

ID=41789133

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100931394A Expired - Fee Related CN101660992B (en) 2009-09-29 2009-09-29 A method for rapidly detecting the sedimentation velocity of algae cells

Country Status (1)

Country Link
CN (1) CN101660992B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105987866A (en) * 2015-02-09 2016-10-05 中国科学院上海药物研究所 Heterogeneous liquid settlement automatic-monitoring method and device
CN109085098A (en) * 2017-06-14 2018-12-25 北京精微视觉科技有限公司 A kind of detection method and detection device of cell settlement progress
CN108709836B (en) * 2018-07-26 2024-03-19 宁夏大学 Aerosol detection method and system
FR3091924B1 (en) * 2019-01-23 2021-10-15 Erba Diagnostics Ltd Method for optimizing the concentration of elements of interest for visual measurements on a biological sample
MX2022004703A (en) * 2019-10-22 2022-05-11 S D Sight Diagnostics Ltd ACCOUNTING FOR ERRORS IN OPTICAL MEASUREMENTS.
CN112033868B (en) * 2020-09-01 2025-01-17 湖南伊鸿健康科技有限公司 Method for accelerating cell precipitation, metal plate and microscope
CN112691710A (en) * 2020-12-03 2021-04-23 中山大学 Micro-fluidic system
CN113074657B (en) * 2021-03-16 2022-02-01 华中科技大学 Three-coordinate phase mapping table calibration method based on virtual plane sampling
CN113092346B (en) * 2021-04-06 2022-02-11 中国水利水电科学研究院 Algae cell counting detection system and detection method thereof
CN116153388B (en) * 2023-04-19 2023-09-19 北京蓝晶微生物科技有限公司 Quantitative relation model and application, cell state monitoring method, device and system

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1196485A (en) * 1998-01-21 1998-10-21 天津中医学院 Method for analysing sedimentation curve and red blood cell settling

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1196485A (en) * 1998-01-21 1998-10-21 天津中医学院 Method for analysing sedimentation curve and red blood cell settling

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A. E. Walsby.The form resistance of chitan fibres attached to the cells.《European Journal of Phycology》.1977,第12卷215-223. *
JP昭62-231140A 1987.10.09
JP特开平8-136447A 1996.05.31

Also Published As

Publication number Publication date
CN101660992A (en) 2010-03-03

Similar Documents

Publication Publication Date Title
CN101660992B (en) A method for rapidly detecting the sedimentation velocity of algae cells
Bisgaard Velocity fields around spheres and bubbles investigated by laser-Doppler anemometry
CN108444897B (en) Label-free microfluidic cytometer based on light sheet illumination and sheath flow technology and method
CN101842751A (en) Tracking and characterizing particles with holographic video microscopy
CN103267543B (en) Microscopic particles image micro-flow measurement instrument and method
CN103822868A (en) Device and method for detecting grain size of phytoplankton in seawater
CN103063626A (en) Light path auto-correction cell laser excitation detecting device and detecting method thereof
CN101122555A (en) High concentration super fine granule measuring device and method based on backward photon related spectrum
CN104155426B (en) A quantitative monitoring system and method for soil splash erosion dynamics
CN107703027A (en) A kind of method based on quantum dot three-dimensional tracer determination cytoplasm viscosity
CN103592103A (en) Mini-channel liquid-solid two-phase flow parameter measurement device and method based on laser extinction method
CN104390895B (en) Microimaging-based method for measuring particle size by utilizing image gray scale
CN111198188A (en) A method for determining the concentration of surfactant solution based on image acquisition
Harvey et al. Optical coherence tomography velocimetry in controlled shear flow
KR20040010141A (en) Collecting and measuring method of yellow sand particles
JP6264229B2 (en) Bubble diameter distribution measuring method and bubble diameter distribution measuring apparatus
Singh et al. Analysis of cellular structure by light scattering measurements in a new cytometer design based on a liquid-core waveguide
EP1729109A1 (en) Method and equipment for the determination of surfactant concentrations in aqueous solutions by determining contact angle.
CN209342750U (en) A kind of experimental facilities using Laser Doppler speed
CN201876402U (en) Underwater optical measuring device for field measurement of grain diameter of phytoplankton
RU115484U1 (en) DEVICE FOR DETERMINING THE DISTRIBUTION OF PARTICLES OF DISPERSED MATERIALS BY DIMENSIONS
CN105891153B (en) Pulsed sieve-plate extraction column hydraulic performance measurement method based on four probe of optical fiber
CN205157378U (en) Large -scale portable suspended solid section high -resolution is measured and slidable sampling device
CN112033867B (en) On-site counting method for phytoplankton in seawater
Agrawal et al. Grain size distribution and sediment flux structure in a river profile, measured with a LISST-SL Instrument

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

Termination date: 20120929