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
images
interpolation
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

The present invention relates to a thin layer chromatography quantitative analysis method based on image processing technology, which belongs to the field of digital image processing and chemical analysis. The present invention is characterized in that a digital camera or a video camera is used to shoot a thin layer plate under the illumination condition of different wave length, the distortion of scenes of images which are shot are corrected, the filter of the noise of the images is preprocessed, the characteristic of uniform spot distribution when samples are spotted and the equipartition principle of a longitudinal distance are used to determine control spots, an interpolation method is used to construct an image background according to the control spots, background images which are constructed are eliminated from the thin layer images of which the scene distortion is corrected, the brightness of prime spots of each of the background images which are constructed is normalized, thin layer images from which the background images are eliminated are separated, and a gray scale in each region is integrated according to a separation result. The method can correctly describe the background information of the images which are provided, separate the target images from the images, make the illumination brightness of each spot corrected, and have the advantages of very high image separation precision, very high precision of a sample quantitative analysis result, very good repeatability and easy popularization.

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
CN100347541C (en) Thin layer chromatography quantitative analysis method based on image processing technology
US11721018B2 (en) System and method for calculating focus variation for a digital microscope
EP3988985B1 (en) Fast auto-focus in microscopic imaging
US20150256739A1 (en) Fast Auto-Focus in Imaging
JP6787076B2 (en) Color reaction detection system, color reaction detection method and program
CN107462741B (en) Moving object speed and acceleration measuring device
US9253449B2 (en) Mosaic picture generation
CN106030653A (en) Image processing system and method for generating high dynamic range image
CN106408616B (en) The inconsistent bearing calibration of perspective view background in a kind of CT imaging
CN1411046A (en) Pattern evaluating apparatus, pattern evaluating method and program
JP2023534175A (en) Neural network analysis of LFA specimens
CN108548824A (en) PVC facial mask detection method and device
JP2010276540A (en) Living tissue surface analyzer, living tissue surface analysis program, and living tissue surface analysis method
US6535624B1 (en) Gel electrophoresis image combining for improved dynamic range
US8705698B2 (en) X-ray analyzer and mapping method for an X-ray analysis
CN1472527A (en) Multiple-marknig-object biological chip signal analyzing systems
CN110308117B (en) Diffraction distance calibration method and system and image reconstruction method and system
CN105427315B (en) Digital instrument image position testing method and device
CN115616018B (en) Positioning method and device for scanning electron microscope, electronic equipment and storage medium
KR20070091236A (en) Defective particle measuring apparatus and defective particle measuring method
JP2010276492A (en) Method and apparatus for measuring distribution of specular reflection of point image
Mascalchi et al. Which elements to build co-localization workflows? from metrology to analysis
JP2000019064A (en) Film evaluation method and film evaluation apparatus
CN103646390A (en) Image processing system and method based on multi-level image positioning
CN111060437A (en) Detection equipment and detection method for 3D printing sand box

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