[go: up one dir, main page]

CN100347541C - Thin layer chromatography quantitative analysis method based on image processing technology - Google Patents

Thin layer chromatography quantitative analysis method based on image processing technology Download PDF

Info

Publication number
CN100347541C
CN100347541C CNB200510115774XA CN200510115774A CN100347541C CN 100347541 C CN100347541 C CN 100347541C CN B200510115774X A CNB200510115774X A CN B200510115774XA CN 200510115774 A CN200510115774 A CN 200510115774A CN 100347541 C CN100347541 C CN 100347541C
Authority
CN
China
Prior art keywords
image
thin layer
gray
interpolation
pixel
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
CNB200510115774XA
Other languages
Chinese (zh)
Other versions
CN1766607A (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 CNB200510115774XA priority Critical patent/CN100347541C/en
Publication of CN1766607A publication Critical patent/CN1766607A/en
Application granted granted Critical
Publication of CN100347541C publication Critical patent/CN100347541C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/90Plate chromatography, e.g. thin layer or paper chromatography
    • G01N30/95Detectors specially adapted therefor; Signal analysis

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Image Analysis (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明属于数字图像处理和化学分析领域,其特征在于利用数码相机或摄像机在不同波长的光照条件下对薄层板进行拍摄;对所拍摄的图像进行镜头畸变校正、噪声虑除预处理;利用点样时点均匀分布的特点以及纵向距离等分原则确定控制点;根据控制点再利用插值方法构造出图像背景;将构造出的背景图像从镜头畸变校正后的薄层图像中去除,并根据所构造的背景图像对每一像素点进行亮度归一化;对去掉背景图像后的薄层图像进行分割;根据分割结果对每个区域内的灰度进行积分。本方不仅能够准确的描述所给图像背景信息,分割出目标图像,而且还可以对每一点的光照强度进行修正,图像分割以及样品定量分析结果精度都很高,重复性也非常好,易于推广。

Figure 200510115774

The invention belongs to the field of digital image processing and chemical analysis, and is characterized in that a digital camera or video camera is used to shoot thin-layer plates under illumination conditions of different wavelengths; lens distortion correction and noise removal preprocessing are performed on the captured images; The control point is determined by the characteristics of uniform distribution of points and the principle of equal division of longitudinal distance; according to the control point, the image background is constructed by interpolation method; the constructed background image is removed from the thin-layer image after lens distortion correction, and according to The brightness of each pixel is normalized for the constructed background image; the thin-layer image after removing the background image is segmented; the gray level in each region is integrated according to the segmentation result. Our method can not only accurately describe the background information of the given image and segment the target image, but also correct the light intensity of each point. The results of image segmentation and quantitative analysis of samples are very accurate and repeatable, and easy to popularize .

Figure 200510115774

Description

Thin layer chromatography quantitative analysis method based on image processing techniques
Technical field
The invention belongs to computer digital image process field and chemical analysis field.
Background technology
In fields such as chemical analysis, biomedicine, medical research, clinical medicine, food processing, detection to certain micro substance composition is a very necessary job with evaluation, wherein thin-layered chromatography or thin layer chromatography (Thin LayerChromatography is called for short TLC) are widely used.So-called thin-layered chromatography is taped against adsorbent (for example silica gel or aluminium oxide) exactly on the planar object and forms skim as glass plate, plastic sheet or metal forming, with solvent sample is launched and the method for separating then, it is a kind of chromatographic process easy, economic, micro-, applied widely.Thin-layered chromatography mainly is made up of three steps: the expansion of sample behind the point sample of sample, the point sample, the back measurement to sample of expansion.The point sample of sample has become very proven technique with expansion, and the sample after launching is measured, in the middle of particularly accurate quantitative measurment technology but still is in and constantly explores.At present, thin layer chromatography quantitative analysis has all adopted the thin-layer chromatogram scanner method basically.Though thin-layer chromatogram scanner is widely used.But thin-layer chromatogram scanner exists following shortcoming:
(a) measuring speed is slow, measures a thin layer plate and will spend dozens of minutes, is unfavorable for the detection of a large amount of samples;
(b) precision is lower, particularly big, irregular blotches for area, and error is bigger;
(c) can't directly take pictures or record a video, reduced confidence level;
(d) not directly perceived, because of seeing the location situation of scanning, easily cause the location inaccurate, thereby produce than mistake;
(e) complicated operation is because therefore the optical texture more complicated of thin-layer chromatogram scanner adjusts very difficulty;
(f) cost height, because the optical texture more complicated of thin-layer chromatogram scanner, thereby it is many to make operation, thereby the cost height;
Recent years, there is the scholar to propose to utilize flat image reading apparatus method that the sample on the thin layer plate is carried out quantitative test in the world and measures.Someone utilizes the thin layer plate after flat image reading apparatus will dye to scan, and obtains handling after the digital picture again.They with resolution be the scanner of 100DPI with thin layer image scanning in computing machine, utilize computer programming to carry out quantitative test then.Their conclusion is that the method not only can be carried out qualitative analysis, also can carry out quantitative test.For flat image reading apparatus, owing to utilize the linear light sources motion scan, its scanning ray distributes more even, thereby the effect that its thin layer plate after to dyeing obtains when handling is relatively good.But flat image reading apparatus can not scan the thin layer chromatography board that needs special light sources to carry out imaging, as fluorescence irradiation, ultraviolet ray irradiation etc.In addition, this disposal route is time-consuming long, also can't be integrated with other instruments.
Solve the problem that flat image reading apparatus brings, people will expect video camera naturally.In fact, before flat image reading apparatus is used, utilizes video camera imaging that thin-layer chromatography is carried out quantitative test and just be studied, this is because the image taking speed of video camera is fast, easy for installation, and is easy to and the integrated use of Other Instruments.Yet after having adopted video camera to carry out thin-layer chromatographic analysis, the conclusion that people provide is but not very good.Reason is to compare with traditional thin-layer chromatogram scanner, and the correctness of video camera imaging is subjected to the restriction of illumination uniformity, thereby the result is not very accurate, and repeatability is bad.Except the reason that they point out, during video camera imaging, light source is indefinite, aperture is indefinite, film speed is indefinite has also brought difficulty to research, and just repeatability is not strong.These problems make in the quantitative test of video camera imaging method in thin-layer chromatography to can not get using preferably just.
Summary of the invention
From top analysis, if the shortcoming of uneven illumination in the video camera imaging process can be remedied accordingly, video camera imaging just has great vitality in thin-layer chromatography: image taking speed is fast, analysis precision is high, visual result, be easy to the storage, with low cost.The phenomenon of uneven illumination is intrinsic in the video camera imaging process, and for this reason, the means of what use is made of Digital Image Processing remedy the defective that this phenomenon is brought, key point of the present invention that Here it is.
In addition, the present invention also aims to provide a kind of a kind of image partition method based on background estimating.
Ultimate principle of the present invention is: when sample launching on the thin layer plate and drying after, just this thin layer plate is placed in the lamp box with multi-wavelength light photograph, utilize digital camera or the video camera formation digital picture of taking pictures then, at last estimate the intensity of illumination that every bit on the thin layer plate, according to the actual content that the intensity of illumination of estimating splits sample and the integration of gray scale is determined sample per sample by computing machine.Concrete steps comprise utilizes video camera under the illumination condition of different wave length thin layer plate to be taken; Captured image is carried out lens distortion calibration, noise considers and pre-service such as removes; Utilize vertical pixel grey scale integral projection method to determine the reference mark; According to resulting reference mark, utilize interpolation method to construct image background again; To remove the thin layer image of background image behind lens distortion calibration that construct, according to the background image of being constructed each pixel is carried out intensity of illumination normalization simultaneously; The thin layer image that removes behind the background image is cut apart; According to segmentation result the gray scale in each sample area is carried out integration; Linear characteristics of gray integration and sample quality per sample, contrast standard product again calculate the actual content of each sample.
The present invention is characterised in that, successively by following step:
The thin layer plate that step 1. pair is loaded with sample utilizes digital camera or video camera to take under selected wavelength illumination condition, is stored in computing machine or the DSP disposal system after forming digital picture;
Step 2. computing machine carries out pre-service according to the following steps to captured image:
Step 2.1. lens distortion calibration: take the standard graticule with digital camera or video camera, concern one to one according to the grid intersection point again, the image rectification that step 2.1 is obtained becomes the standard graticule, notes the correction parameter of each point, is provided with doing to proofread and correct when the thin layer plate image is taken in the back using;
Step 2.2. consider to remove noise, adopts any method in the following filtering method to carry out: smothing filtering or low-pass filtering or Wiener filtering, thus obtain the digital picture of thin layer plate;
Be along perpendicular to equally distributed these characteristics of expansion direction when step 3. is utilized point sample, after the pre-service that step 2 obtains, find the horizontal ordinate at reference mark the image, according to the principle of fore-and-aft distance five equilibrium, determine the ordinate at reference mark again, thus controlled point;
The structural setting image is come at the reference mark that step 4. utilizes any method in following each interpolation method to obtain according to step 3: described interpolation method is Newton interpolation or the piecewise low-order interpolation or the spline interpolation of Lagrange's interpolation or inequality;
Step 5. deducts the background image that is constructed from the thin layer plate digital picture that step 2 obtains, if f is (x, y) be pixel (x, the gray scale of y) locating, and g (x in the thin layer plate digital picture behind the lens distortion calibration, y) be pixel (x in the background image of estimating, y) gray scale of locating, then subtract each other the pixel that obtains image (x, the gray-scale value of y) locating is (g (x, y)-f (x, y)); According to the background image of being constructed every bit is carried out normalization again, obtain the thin layer plate digital picture of homogeneous background, (x after the gray scale correction of y) locating is corresponding point in this image
Figure C20051011577400051
Perhaps
Figure C20051011577400052
Wherein, I MaxBe the highest gray scale of the background image of being constructed, I MinBe minimal gray, I bFor same point in the described background image (x, the gray-scale value of y) locating, promptly every bit (x, y) actual light intensity and bias light are inversely proportional to;
Step 6. utilizes the global threshold method to be partitioned into target image, and described global threshold gets 2 or 3;
Interior each pixel in each zone adds up with gray-scale value in the image that step 7. pair step 6 obtains, and promptly carries out gray integration;
Step 8. is determined quantitative relationship between each sample spot according to the size of the gray integration value that step 7 obtains by the linear ratio relation.Concrete computing method are: the gray integration value of establishing a certain sample is S, and the integrated value of standard items is T, and the quality of standard items is G, and then the quality of sample is
Figure C20051011577400061
Adopt this method to handle and to describe the image background information of giving accurately, be partitioned into target image, but also can carry out normalization to the intensity of illumination of every bit, segmentation result reaches very high accuracy rate, thereby makes the result of quantitative test reach higher precision.
The invention will be further described below in conjunction with accompanying drawing.
Description of drawings:
Fig. 1 is a hardware composition frame chart of the present invention;
Fig. 2 is an overall flow block diagram of the present invention;
Fig. 3 is a thin layer plate image collecting device of the present invention;
Fig. 4 is lens distortion calibration shooting figure of the present invention;
Fig. 5 is the former figure of thin layer plate image taking that the present invention collects;
Fig. 6 is that the 3 D stereo of the thin layer plate image gathered of the present invention shows;
Fig. 7 is a filtering template of the present invention;
Fig. 8 is the synoptic diagram that the present invention is used to select the reference mark, takes figure based on Fig. 5;
Fig. 9 is a 3-D display of utilizing institute of the present invention construct image background;
Figure 10 is the 3-D display result of the thin layer plate image after the removal uneven illumination background that obtains of the present invention;
Figure 11 is a sample segmentation result of the present invention.
Embodiment:
Fig. 1 has provided hardware composition frame chart of the present invention.The thin layer plate 1 that is loaded with sample is placed in the lamp box 2, by video camera or digital camera 3 it is taken pictures then, and the image of formation is sent in computing machine or the DSP disposal system 4, the result of processing be sent to that printer 6 is printed or miscellaneous equipment 7 in standby.Monitor 5 among Fig. 1 is in order to show captured image and program operation process.
Fig. 2 has provided workflow block diagram of the present invention, comprises taking thin layer plate image (first step); Photographic images is done image pre-service such as the abnormal correction of lens distortion, noise filtering (second step); Utilize the vertical pixel grey scale integral projection of image to determine the horizontal ordinate at reference mark, utilize the ordinate of obtaining the reference mark apart from five equilibrium then, thus controlled point (the 3rd step); Produce background image (the 4th step) according to above-mentioned reference mark interpolation; From original image, deduct the background that interpolation generates, and each pixel gray scale is carried out correcting process (the 5th step) according to background image; Utilize connected domain or overall thresholding split image (the 6th step); Calculate each regional gray integration (the 7th step) according to segmentation result, according to the gray integration content of calculation sample (the 8th step) as a result.
The first step of the present invention is to take the thin layer plate image, and Fig. 3 has provided concrete shooting process.The thin layer plate a that is loaded with sample is placed on and can selects illumination wavelength (wavelength is 254nm, 312nm, 365nm, and white light) among the lamp box b, then by the digital camera that optical filter is housed or video camera c to thin layer plate take pictures form digital picture and be stored in computing machine or the DSP disposal system in.Fig. 5 has provided a captured thin layer plate image, and wherein f is the sample that will measure.In order to further specify the influence that the analysis of uneven background specified rate brings, the 3 D stereo display effect that Fig. 6 has provided Fig. 5 (for clarity sake, pixel grey scale among Fig. 6 is by anti-phase), g is the position of sample in the three-dimensional plot among the figure, the influence of background is fairly obvious as can be seen.
After obtaining the digital picture of thin layer plate, in order to proofread and correct the influence that causes by lens distortion, and remove noise, must carry out image pre-service such as lens distortion calibration, noise filtering (second step) to resulting thin layer plate digital picture.The antidote of lens distortion is with digital camera or video camera.Shooting standard graticule obtains fault image, shown in Fig. 4 (a), concern one to one according to the grid intersection point then, corresponding relation as e and d, distorted image correction is become the standard graticule, shown in Fig. 4 (b), note correction parameter, when taking the thin layer plate image afterwards, just proofread and correct with this correction parameter.Noise filtering then adopts smothing filtering, or low-pass filtering, or filtering method such as Wiener filtering.For example adopt 3 * 3 smothing filtering template to carry out smothing filtering, Fig. 7 has provided this 3 * 3 template, when specifically calculating K 0 = 1 9 Σ i = 0 8 K i , K wherein iIt is the gray-scale value of each point in the template.
In order to estimate the uneven illumination background that the thin layer plate image, the 3rd step of the present invention and the 4th step search out the horizontal ordinate at reference mark exactly from pretreated thin layer plate digital picture, utilize the ordinate of obtaining the reference mark apart from halving method then, thereby controlled point, this was the 3rd step, then utilize Lagrange's interpolation according to asking reference mark again, or inequality and Newton's interpolation formula, or difference and equidistant knot interpolation, or Hermite interpolation, or piecewise low-order interpolation, or interpolation method such as spline interpolation constructs image background figure, and this was the 5th step.Fig. 8 has provided the synoptic diagram of seeking the reference mark from pretreated thin layer plate digital picture.Because sample launches along equidirectional, thereby we are in point sample, can carry out equidistant point sample along direction perpendicular to expansion direction, when measuring, equidistantly divide earlier according to the point sample port number, find the horizontal ordinate of control vertex, shown in the perpendicular line h among Fig. 8 (a), the position that these perpendicular line passed through is exactly an image background, gets from these straight lines and a little just can obtain the reference mark of structural setting image in the future.After obtaining the horizontal ordinate at reference mark, by the mode that is not more than the horizontal ordinate spacing image is carried out vertically equidistant the division again, thereby determine the ordinate at reference mark, Fig. 8 (b) has provided the last synoptic diagram at reference mark, and each point of crossing i is exactly the reference mark that we will look for.Control vertex has been arranged, utilized curved surface structure method above-mentioned just can construct image background.For example, utilize Catmull-Rom spline surface construct image background 3 times, formula is:
P ( u , v ) = Σ i = 0 n Σ j = 0 m P ij B i , 3 ( u ) B j , 3 ( v ) , 0 ≤ u ≤ 1,0 ≤ v ≤ 1
P in the following formula Ij(i=0,1 ..., n; J=0,1 ..., m) being the reference mark of finding in the step 3, u, v are interpolation parameter, the size of its value affects the density of interpolation point, B I, 3(u) be harmonic function, their expression is as follows:
B 0,3(u)=-τ+2τu 2-τu 3
B 1,3(u)=1+(τ-3)u 2+(τ-2)u 3
B 2,3(u)=τu+(3-2τ)u 2+(τ-2)u 3
B 3,3(u)=-τu 2+τu 3
τ is the dynamics control coefrficient, 0≤τ≤1, and general value is 0.2.
Fig. 9 is that the 3 D stereo of the background image j that produces according to the determined reference mark of Fig. 8 and by Catmull-Rom spline surface interpolation shows.For clarity sake, the pixel grey scale among Fig. 9 is by anti-phase.
Estimate uneven background image after, the thin layer plate digital picture behind lens distortion calibration deducts the thin layer plate digital picture that the background image that is constructed has just obtained homogeneous background, Here it is of the present invention the 5th the step.Concrete subtractive method is, if f is (x, y) be pixel (x, the gray scale of y) locating, and g (x in the thin layer plate digital picture behind the lens distortion calibration, y) be pixel (x in the background image of estimating, y) gray scale of locating, then subtract each other the pixel that obtains image (x, the gray-scale value of y) locating is (g (x, y)-f (x, y)).Figure 10 has provided the 3-D display of the thin layer plate image that removes non-homogeneous background, wherein k representative sample point.
In addition, the background of estimating has been arranged, we have also just obtained the relative light intensity value of each pixel, thereby can carry out the light intensity correction.The method of revising is: the highest gray scale of establishing estimation background j is I Max, minimum is I Min, (x, gray-scale value y) are I to a certain pixel among the j bIt is, corresponding among Figure 10 that (x, the gray-scale value of y) locating are I, and then corresponding among Figure 10 (x, the gray-scale value of y) locating is corrected for
Figure C20051011577400091
Perhaps Just the actual light intensity of every bit and bias light are inversely proportional to.
The thin layer plate image that removes non-homogeneous background has been arranged, utilized the global threshold method to be partitioned into target image (the 6th step), and remove isolated point.Global threshold is 0 in theory, but in fact because the existence of noise generally gets 2 or 3.Figure 11 is exactly last segmentation result, and 1 is exactly the sample area that is partitioned into.
According to top segmentation result, we just can carry out gray integration to respective regions in Figure 10, just the gray-scale value of each pixel in the zone are added up, and this was the 7th step.
At last, just can determine quantitative relationship between each sample spot by linear scaling according to the size of gray integration value, this is final step of the present invention, i.e. the 8th step.Concrete computing method are: the gray integration value of establishing a certain sample is S, and the integrated value of standard items is T, and the quality of standard items is G, and then the quality of sample is
Figure C20051011577400093

Claims (1)

1. based on the thin layer chromatography quantitative analysis method of image processing techniques, it is characterized in that:
The thin layer plate that step 1. pair is loaded with sample utilizes digital camera or video camera to take under selected wavelength illumination condition, is stored in computing machine or the digital signal processor processes system after forming digital picture;
Step 2. computing machine carries out pre-service according to the following steps to captured image:
Step 2.1. lens distortion calibration: take the standard graticule with digital camera or video camera, concern one to one according to the grid intersection point again, become the standard graticule taking the resulting image rectification of standard graticule with digital camera or video camera in the step 2.1, note the correction parameter of each point, be provided with doing to proofread and correct when the thin layer plate image is taken in the back using;
Step 2.2. filtering noise, adopt any method in the following filtering method to carry out: smothing filtering or low-pass filtering or Wiener filtering, thus obtain the digital picture of thin layer plate;
Be along perpendicular to equally distributed these characteristics of expansion direction when step 3. is utilized point sample, after the pre-service that step 2 obtains, find the reference mark horizontal ordinate of using for interpolation the image, again according to the principle of fore-and-aft distance five equilibrium, determine the reference mark ordinate that the confession interpolation is used, thereby obtain the reference mark that the confession interpolation is used;
The structural setting image is come at the reference mark that step 4. utilizes any method in following each interpolation method to obtain according to step 3: described interpolation method is Newton interpolation or the piecewise low-order interpolation or the spline interpolation of Lagrange's interpolation or inequality;
Step 5. deducts the background image that is constructed from the thin layer plate digital picture that step 2 obtains, if f is (x, y) be pixel (x, the gray scale of y) locating, and g (x in the thin layer plate digital picture behind the lens distortion calibration, y) be pixel (x in the background image of estimating, y) gray scale of locating, then subtract each other the pixel that obtains image (x, the gray-scale value of y) locating is (g (x, y)-f (x, y)); According to the background image of being constructed each pixel is carried out normalization again, obtain the thin layer plate digital picture of homogeneous background, (x after the gray scale correction of y) locating is respective pixel in this image
Figure C2005101157740002C1
Perhaps
Figure C2005101157740002C2
Wherein, I is pixel (x, gray-scale value y), the I that subtracts each other resultant image MaxBe the highest gray scale of the background image of being constructed, I MinBe the minimal gray of the background image of being constructed, I bFor same point in the described background image (x, the gray-scale value of y) locating, promptly every bit (x, y) actual light intensity and bias light are inversely proportional to;
Step 6. utilizes the global threshold method to be partitioned into target image, and described global threshold gets 2 or 3;
Interior each pixel in each zone adds up with gray-scale value in the image that step 7. pair step 6 obtains, and promptly carries out gray integration;
Step 8. is determined quantitative relationship between each sample spot according to the size of the gray integration value that step 7 obtains by the linear ratio relation; Concrete computing method are: the gray integration value of establishing a certain sample is S, and the integrated value of standard items is T, and the quality of standard items is G, and then the quality of sample is
CNB200510115774XA 2005-11-11 2005-11-11 Thin layer chromatography quantitative analysis method based on image processing technology Expired - Fee Related CN100347541C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB200510115774XA CN100347541C (en) 2005-11-11 2005-11-11 Thin layer chromatography quantitative analysis method based on image processing technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB200510115774XA CN100347541C (en) 2005-11-11 2005-11-11 Thin layer chromatography quantitative analysis method based on image processing technology

Publications (2)

Publication Number Publication Date
CN1766607A CN1766607A (en) 2006-05-03
CN100347541C true CN100347541C (en) 2007-11-07

Family

ID=36742614

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB200510115774XA Expired - Fee Related CN100347541C (en) 2005-11-11 2005-11-11 Thin layer chromatography quantitative analysis method based on image processing technology

Country Status (1)

Country Link
CN (1) CN100347541C (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663735A (en) * 2012-03-15 2012-09-12 天津大学 Quantitative assessment method for image heterogeneity characteristics

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101847202B (en) * 2009-03-23 2012-11-28 凯迈(洛阳)电子有限公司 Correction algorithm for image distortion of optical fingerprint collector
CN102881007B (en) * 2012-08-15 2016-01-20 百正药业股份有限公司 The image processing method of compound planar separation result and system thereof
CN104502519B (en) * 2014-12-23 2016-05-25 厦门海荭兴科技股份有限公司 A kind of thin-layer chromatography Fast Determination of Pesticide Residue method based on image processing
CN104792915B (en) * 2015-04-08 2017-01-11 广州浩意万医药科技有限公司 Method for determining nystose content of morinda officinalis roots
CN106546692A (en) * 2016-10-25 2017-03-29 烟台大学 A kind of method of imaging method detection organophosphorus pesticide residue and application thereof
CN106546693B (en) * 2016-10-25 2018-05-08 烟台大学 A kind of method of portable imaging method detection sulfa drug residue and application thereof
JP7298993B2 (en) 2018-04-09 2023-06-27 浜松ホトニクス株式会社 Specimen observation device and specimen observation method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3924948A (en) * 1973-12-17 1975-12-09 Kontes Glass Co Densitometer for use in quantitative thin layer chromatographic analysis
US6195449B1 (en) * 1997-05-18 2001-02-27 Robert Bogden Method and apparatus for analyzing data files derived from emission spectra from fluorophore tagged nucleotides
CN1369848A (en) * 2001-02-13 2002-09-18 科学与工业研究会 Method of chromatogram fingerprint atlas, single medicine and preparation standardization

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3924948A (en) * 1973-12-17 1975-12-09 Kontes Glass Co Densitometer for use in quantitative thin layer chromatographic analysis
US6195449B1 (en) * 1997-05-18 2001-02-27 Robert Bogden Method and apparatus for analyzing data files derived from emission spectra from fluorophore tagged nucleotides
CN1369848A (en) * 2001-02-13 2002-09-18 科学与工业研究会 Method of chromatogram fingerprint atlas, single medicine and preparation standardization

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
凝胶色谱图象计算机分析技术用于DNA检测的实验研究 张思祥等.生命科学仪器,第1卷第2期 2003 *
田基黄薄层色谱分类的图像分析方法研究 许建新等.福建电脑,第2期 2005 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663735A (en) * 2012-03-15 2012-09-12 天津大学 Quantitative assessment method for image heterogeneity characteristics
CN102663735B (en) * 2012-03-15 2014-08-20 天津大学 Quantitative assessment method for image heterogeneity characteristics

Also Published As

Publication number Publication date
CN1766607A (en) 2006-05-03

Similar Documents

Publication Publication Date Title
US11721018B2 (en) System and method for calculating focus variation for a digital microscope
US7835561B2 (en) Method for image processing and reconstruction of images for optical tomography
CN100557424C (en) Image inspection method and device
CN101178369B (en) Digital device for industrial radiographic inspection film
CN106408616B (en) A Correction Method for Inconsistency of Projection Image Background in CT Imaging
KR20140045331A (en) Fast auto-focus in microscopic imaging
US20130010175A1 (en) Method and device for analyzing the optical quality of a transparent substrate
JP2023534175A (en) Neural network analysis of LFA specimens
JP6908624B2 (en) Systems and methods for real-time volume control
CN107403177A (en) Brightness measurement method based on industrial camera
CN100347541C (en) Thin layer chromatography quantitative analysis method based on image processing technology
CN111784684A (en) Laser-assisted transparent product internal defect depth setting detection method and device
CN109325927A (en) Brightness compensation method for photogrammetry images of industrial cameras
CN1411046A (en) Pattern evaluating apparatus, pattern evaluating method and program
JP2010276540A (en) Living tissue surface analyzer, living tissue surface analysis program, and living tissue surface analysis method
KR100926019B1 (en) Defective particle measuring apparatus and defective particle measuring method
CN110095472B (en) HDRI-based high-reflectivity metal surface defect detection method and system
JP2005293334A (en) Template matching device
CN110308117B (en) Diffraction distance calibration method and system and image reconstruction method and system
Mascalchi et al. Which elements to build co-localization workflows? from metrology to analysis
JP2000019064A (en) Film evaluation method and film evaluation device
CN116878430A (en) Pearl layer thickness measuring method based on X-ray micro-focusing detection equipment
CN111060437A (en) Detection equipment and detection method for 3D printing sand box
CN101354357A (en) A method for analyzing micro-density images of tree rings
MUSTAFA et al. FLUDI: Using digital images for measuring fish length

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

Termination date: 20131111