CN106771169B - A kind of colloidal gold detection method and system based on image procossing - Google Patents
A kind of colloidal gold detection method and system based on image procossing Download PDFInfo
- Publication number
- CN106771169B CN106771169B CN201710037576.9A CN201710037576A CN106771169B CN 106771169 B CN106771169 B CN 106771169B CN 201710037576 A CN201710037576 A CN 201710037576A CN 106771169 B CN106771169 B CN 106771169B
- Authority
- CN
- China
- Prior art keywords
- point
- gray
- line positions
- peak
- colloidal
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/53—Immunoassay; Biospecific binding assay; Materials therefor
- G01N33/558—Immunoassay; Biospecific binding assay; Materials therefor using diffusion or migration of antigen or antibody
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
Landscapes
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Chemical & Material Sciences (AREA)
- Hematology (AREA)
- Urology & Nephrology (AREA)
- Biotechnology (AREA)
- Microbiology (AREA)
- Cell Biology (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a kind of colloidal gold detection method and system based on image procossing, it carries out gray proces by the region to be identified of the detection image to colloidal-gold detecting-card and calculates gray integration curve, all trough points of the gray integration curve are searched, and calculate the half-peak breadth of each trough point, peak center point, effective integral area；According to the half-peak breadth of maximum trough point judge the validity of the colloidal-gold detecting-card, and judged and searched C line positions and T line positions according to the position of the peak center point of the maximum trough point, to realize the qualitative analysis of colloidal-gold detecting-card, it is few to the image qualifications of acquisition, C lines, the fixation and recognition analysis of T lines are accurate, good compatibility, it is applied widely；The ratio of the T line positions and the gray integration value of the C line positions is also further calculated, and the ratio is compared with typical ratio, obtains the concentration value of sample to be tested, to realize that the quantitative analysis of colloidal-gold detecting-card, accuracy of detection are high.
Description
Technical field
The present invention relates to colloidal gold immunochromatographimethod technical field, especially a kind of colloidal gold detection side based on image procossing
The system of method and its application this method.
Background technology
In food security scene and the field of fast detection of circulation, colloidal gold immunochromatographimethod technology is easy to operate with its, sentences
Convenient advantage is read, the medicament residue detection in tissue, urine, blood is had been widely used for, such as：Clenbuterol, Lake
Dopamine, tetracycline etc..Its principle is that special antibody is first fixed on to a certain zone of nitrocellulose filter, when the drying
After sample is immersed in nitrocellulose one end, sample will be moved forward along the film, when being moved to the region for being fixed with antibody, sample
Corresponding antigen is specifically bound with the antibody in product, using immuno-gold labeling antigen-antibody, keeps the region aobvious
The certain color of color, to realize the immunodiagnosis of specificity.Testing result is the color depth of the C lines of observing response plate, T lines
Shallow, the color ratio of T lines/C lines is smaller, then test substance concentration is higher, currently, most of only lean on gross visualization to carry out interpretation,
Subjectivity is larger, causes its scope of application by larger limitation.
Invention content
The present invention is simultaneous to solve the above problems, provide a kind of colloidal gold detection method and system based on image procossing
Capacitive is good, applied widely, accuracy of detection is high.
To achieve the above object, the technical solution adopted by the present invention is：
A kind of colloidal gold detection method based on image procossing comprising following steps：
(10) detection image of colloidal-gold detecting-card is obtained, and the region to be identified of the detection image is carried out at gray scale
Reason, obtains gray level image；
(20) the gray integration curve of the gray level image is calculated；
(30) all trough points of the gray integration curve are searched, and calculate the half-peak breadth of each trough point, peak center
Point, effective integral area；
(40) the maximum trough point of effective integral area is searched, maximum trough point is obtained；
(50) if the half-peak breadth of the maximum trough point judges that the maximum trough point is invalid except preset range
Point, and it is dead card to regard as the colloidal-gold detecting-card；If the half-peak breadth of the maximum trough point within preset range,
The maximum trough point is judged for available point, and it is effectively to block to regard as the colloidal-gold detecting-card, then executes step (60)；
(60) according to the position of the peak center point of the maximum trough point judge the maximum trough point for C line positions or
Person's T line positions；If being judged as C line positions, to other significant wave valley points of T lines direction finding, if finding other significant wave valley points
The trough point is then regarded as into T line positions, no T lines are regarded as if it can not find other significant wave valley points；If being judged as T lines position
It sets, then to other significant wave valley points of C lines direction finding, the trough point is regarded as into C lines position if finding other significant wave valley points
It sets；No C lines are regarded as if it can not find other significant wave valley points, and it is dead card to regard as the colloidal-gold detecting-card.
Preferably, in the step (10), also further the detection image is stretched or is compressed to the length of reaction plate
Degree, and according to the detection zone of the reaction plate detection image is cut to obtain the area to be identified of the detection image
Domain.
Preferably, the detection image is cut to obtain the region to be identified of the detection image, refers to described
Centered on reaction plate, retain the pixel region of transverse direction M*L, longitudinal direction M*L*0.3 as the region to be identified；Wherein, M is calibration
Coefficient, L are the length of the reaction plate；The value of the M is ranging from：0.75~0.95.
Preferably, the gray integration curve that the gray level image is calculated in the step (20), further comprises：
(21) gray value [X, Y, Z] of the gray level image is obtained, wherein X is that the pixel of the long side of the reaction plate is sat
Mark, Y are the pixel coordinate of the short side of the reaction plate, and Z is the gray value at coordinate position (X, Y)；
(22) using X-coordinate as abscissa, using Y coordinate as ordinate, gray integration summation is carried out to Y coordinate, obtains ash
Spend integral curve [X, Zsum]；
(23) the gray integration curve is filtered.
Preferably, in the step (30), the half-peak breadth, peak center point, effective integral area computational methods such as
Under：
(31) peak center point is calculated：
Centered on trough point, inflection point is found to both sides respectively, calculates inflection point and trough point ordinate difference Z1, Z2, and
Abscissa difference X1, X2 is found according to the ordinate difference, it is peak center point to take the central point of X1, X2；
(32) half-peak breadth is calculated：
Straight line is done with X1, X2, the distance of trough point to straight line is peak height, and calculates half-peak breadth according to peak height；
(33) effective integral area is calculated：
Calculate the summation of the peak height of all the points between X1, X2, as the effective integral area of trough point.
Preferably, it carries out calculating peak center point in the step (31), also further to ordinate difference Z1, Z2
Calibrated, that is, ordinate difference after calibration is N*Z1, N*Z2, and according to ordinate difference N*Z1, N*Z2 after calibration into
Row finds abscissa difference X1, X2；Wherein, N is calibration factor, and the value of N is ranging from：0.75~0.95.
Preferably, also further the half-peak breadth is calibrated, that is, the preset range of the half-peak breadth after calibration is P/V
Between P*V, the half-peak breadth then judges that the maximum trough point is available point, the half-peak breadth within the scope of P/V and P*V
Then judge that the maximum trough point is Null Spot except P/V and P*V ranges；Wherein, P is half-peak breadth preset value, and V is half-peak breadth
Coefficient of alteration；The value of the V is ranging from：1.1~1.8.
Preferably, still further comprise step (70), calculated according to the T line positions and the C line positions described in
The concentration value of detection sample in colloidal-gold detecting-card；Specifically include following steps：
(71) the gray integration value of the T line positions is calculated according to the T line positions in the gray integration curve,
In, when regarding as no T lines, then the gray integration value of T lines is in terms of zero；
(72) the gray integration value of the C line positions is calculated according to the C line positions in the gray integration curve；
(73) ratio R of the T line positions and the gray integration value of the C line positions is calculated；
(74) ratio R is compared with typical ratio Rs, obtains sample to be tested concentration value.
Preferably, the step (74) are the standards of the standard sample by measuring a certain range of series concentration
Ratio R s, each typical ratio Rs corresponds to a unique standard sample concentration value, to obtain standard curve or standard
Deck watch；Then search corresponding typical ratio Rs according to the ratio R, and by the corresponding standard samples of typical ratio Rs
Product concentration value is as the sample to be tested concentration value.
Correspondingly, the present invention also provides a kind of colloidal gold detecting system based on image procossing comprising：
Image pre-processing module obtains the detection image of colloidal-gold detecting-card, and to the area to be identified of the detection image
Domain carries out gray proces, obtains gray level image；
Curve processing module, the gray integration curve for calculating the gray level image；And it is bent to search the gray integration
All trough points of line, and calculate the half-peak breadth of each trough point, peak center point, effective integral area；And search effectively product
Facet accumulates maximum trough point, obtains maximum trough point；
Effective judgement module, if the half-peak breadth of the maximum trough point judges the maximum wave except preset range
Valley point is Null Spot, and it is dead card to regard as the colloidal-gold detecting-card；If the half-peak breadth of the maximum trough point is default
Within the scope of, then the maximum trough point is judged for available point, and it is effectively to block to regard as the colloidal-gold detecting-card；
CT line locating modules carry out judging the maximum trough point according to the position of the peak center point of the maximum trough point
For C line positions or T line positions；If being judged as C line positions, to other significant wave valley points of T lines direction finding, if finding other
The trough point is then regarded as T line positions by significant wave valley point, and no T lines are regarded as if it can not find other significant wave valley points；If sentencing
Break and recognizes the trough point if finding other significant wave valley points then to other significant wave valley points of C lines direction finding for T line positions
It is set to C line positions；No C lines are regarded as if it can not find other significant wave valley points, and it is nothing to regard as the colloidal-gold detecting-card
Effect card.
Preferably, the curve processing module further comprises curve generation unit and curve computation unit, wherein：
For the curve generation unit by the gray value [X, Y, Z] of the acquisition gray level image, X is the long side of reaction plate
Pixel coordinate, Y are the pixel coordinate of the short side of reaction plate, and Z is the gray value at coordinate position (X, Y)；Using X-coordinate as horizontal seat
Mark carries out gray integration summation to Y coordinate, obtains gray integration curve [X, Zsum] using Y coordinate as ordinate；And further
The gray integration curve is filtered；
The curve computation unit carries out calculating the half-peak breadth of the trough point, peak center point, effective integral area, specifically
's：
Calculate peak center point：It is by centered on trough point, finding inflection point to both sides respectively, calculating inflection point and trough point
Ordinate difference Z1, Z2, and abscissa difference X1, X2 is found according to the ordinate difference, it is in peak to take the central point of X1, X2
Heart point；
Calculate half-peak breadth：It is by doing straight line with X1, X2, the distance of trough point to straight line is peak height, and according to peak height
Calculate half-peak breadth；
Calculate effective integral area：It is the summation by calculating the peak height of all the points between X1, X2, as trough point has
Imitate integral area.
Preferably, quantitative analysis module is still further comprised, is counted according to the T line positions and the C line positions
Calculate the concentration value of the detection sample in the colloidal-gold detecting-card；Specifically：
The gray integration value of the T line positions is calculated according to the T line positions in the gray integration curve, wherein recognize
When being set to no T lines, then the gray integration value of T lines is in terms of zero；
The gray integration value of the C line positions is calculated according to the C line positions in the gray integration curve；
Calculate the ratio R of the T line positions and the gray integration value of the C line positions；
The ratio R is compared with typical ratio Rs, obtains sample to be tested concentration value.
Preferably, the ratio R is compared with typical ratio Rs, is by measuring a certain range of series concentration
Standard sample typical ratio Rs, each typical ratio Rs corresponds to a unique standard sample concentration value, to obtain
Standard curve or standard comparison table；Then search corresponding typical ratio Rs according to the ratio R, and by the standard ratio
The corresponding standard sample concentration values of value Rs are as the sample to be tested concentration value.
The beneficial effects of the invention are as follows：
(1) present invention by obtain colloidal-gold detecting-card detection image, and to the region to be identified of the detection image into
Row gray proces obtain gray level image, calculate the gray integration curve of the gray level image, search the gray integration curve
All trough points, and the half-peak breadth of each trough point, peak center point, effective integral area are calculated, and search effective integral area
Maximum trough point obtains maximum trough point；It is carried out judging the colloidal gold detection according to the half-peak breadth of the maximum trough point
The validity of card is further judged and is searched C lines if effectively according to the position of the peak center point of the maximum trough point
Position and T line positions, to realize the qualitative analysis of colloidal-gold detecting-card, few to the image qualifications of acquisition, C lines, T lines
Automatic positioning discriminance analysis is accurate, and good compatibility is applied widely.
(2) present invention also further calculates the ratio of the T line positions and the gray integration value of the C line positions, and will
The ratio is compared with typical ratio, obtains the concentration value of sample to be tested, to realize the quantitative analysis of colloidal-gold detecting-card,
Accuracy of detection is high.
Description of the drawings
Attached drawing described herein is used to provide further understanding of the present invention, and constitutes the part of the present invention, this hair
Bright illustrative embodiments and their description are not constituted improper limitations of the present invention for explaining the present invention.In the accompanying drawings：
Fig. 1 is the detection image comparison signal of the colloidal-gold detecting-card of three kinds of samples to be tested of the specific embodiment of the invention
Figure；
Fig. 2 is the contrast schematic diagram of the corresponding gray integration curve of three kinds of detection images of Fig. 1；
The gray integration curve that Fig. 3 is Fig. 2 carries out the effect diagram after mean value smoothing is filtered；
Fig. 4 is the canonical plotting of the typical ratio of the standard sample of series concentration；
Fig. 5 is that the testing result of tri- kinds of samples to be tested of Fig. 1 compares chart.
Specific implementation mode
In order to keep technical problems, technical solutions and advantages to be solved clearer, clear, tie below
Closing accompanying drawings and embodiments, the present invention will be described in further detail.It should be appreciated that specific embodiment described herein is only used
To explain the present invention, it is not intended to limit the present invention.
As shown in Figure 1, a kind of colloidal gold detection method based on image procossing of the present invention comprising following steps：
(10) detection image of colloidal-gold detecting-card is obtained, and the region to be identified of the detection image is carried out at gray scale
Reason, obtains gray level image (as shown in Figure 1)；
(20) the gray integration curve (as shown in Figure 2) of the gray level image is calculated；
(30) all trough points of the gray integration curve are searched, and calculate the half-peak breadth of each trough point, peak center
Point, effective integral area；
(40) the maximum trough point of effective integral area is searched, maximum trough point is obtained；
(50) if the half-peak breadth of the maximum trough point judges that the maximum trough point is invalid except preset range
Point, and it is dead card to regard as the colloidal-gold detecting-card；If the half-peak breadth of the maximum trough point within preset range,
The maximum trough point is judged for available point, and it is effectively to block to regard as the colloidal-gold detecting-card, then executes step (60)；
(60) according to the position of the peak center point of the maximum trough point judge the maximum trough point for C line positions or
Person's T line positions；If being judged as C line positions, to other significant wave valley points of T lines direction finding, if finding other significant wave valley points
The trough point is then regarded as into T line positions, no T lines are regarded as if it can not find other significant wave valley points；If being judged as T lines position
It sets, then to other significant wave valley points of C lines direction finding, the trough point is regarded as into C lines position if finding other significant wave valley points
It sets；No C lines are regarded as if it can not find other significant wave valley points, and it is dead card to regard as the colloidal-gold detecting-card.
In the step (10), the detection image of colloidal-gold detecting-card is obtained, the colloidal-gold detecting-card of the acquisition
The pixel size of detection image need to include the region of entire reaction plate without limitation, detection image；Also, also further by the inspection
Altimetric image stretches or is compressed to the length L of reaction plate, and in the present embodiment, L is 500 pixels；And according to the reaction plate
Detection zone cuts the detection image to obtain the region to be identified of the detection image.In the present embodiment, the inspection
The method that altimetric image is stretched or compressed is as follows：Using StretchDIBits functions by picture in rectangular area in detection image
The color data that element uses copies in specified target rectangle.If target rectangle size is bigger than source rectangle, to color data
Row and column stretched, to be matched with target rectangle；If target rectangle size is smaller than source rectangle, by using specified light
Ranks are compressed in grid operation.In the present embodiment, carry out cutting the region to be identified, be centered on the reaction plate,
The detection image after stretching/compressing is cut, filtering may interfere with the region of calculating, obtain region to be identified.It is right in the present embodiment
The detection image is cut to obtain the region to be identified of the detection image, refers to being retained centered on the reaction plate
Lateral M*L, longitudinal direction M*L*0.3 pixel region as the region to be identified；Wherein, M is calibration factor, and L is the reaction
The length of plate；The value of the M is ranging from：0.75~0.95.In the present embodiment, M=0.9, L=400 are taken, it is lateral to retain image
For 0.9*400, longitudinal pixel region for being 0.9*400*0.3 is as region to be identified.
The gray integration curve that the gray level image is calculated in the step (20), further comprises：
(21) gray value [X, Y, Z] of the gray level image is obtained, wherein X is that the pixel of the long side of the reaction plate is sat
Mark, Y are the pixel coordinate of the short side of the reaction plate, and Z is the gray value at coordinate position (X, Y)；
(22) using X-coordinate as abscissa, using Y coordinate as ordinate, gray integration summation is carried out to Y coordinate, obtains ash
Spend integral curve [X, Zsum]；
(23) the gray integration curve is filtered (as shown in Figure 3)；It is described to be filtered in the present embodiment
It is filtered using mean value smoothing, to eliminate noise.
In the step (30), the half-peak breadth, peak center point, the computational methods of effective integral area are as follows：
(31) peak center point is calculated：
Centered on trough point, inflection point is found to both sides respectively, calculates inflection point and trough point ordinate difference Z1, Z2, and
Abscissa difference X1, X2 is found according to the ordinate difference, it is peak center point to take the central point of X1, X2；In the present embodiment, also
Further ordinate difference Z1, Z2 is calibrated, that is, the ordinate difference after calibration is N*Z1, N*Z2, and according to school
Ordinate difference N*Z1, N*Z2 after standard carries out finding abscissa difference X1, X2；Wherein, N is calibration factor, and the model of the value of N
Enclose for：0.75~0.95；In the present embodiment, N=0.8 is taken.
(32) half-peak breadth is calculated：
Straight line is done with X1, X2, the distance of trough point to straight line is peak height, and calculates half-peak breadth according to peak height；
(33) effective integral area is calculated：
Calculate the summation of the peak height of all the points between X1, X2, as the effective integral area S of trough point.
In the step (50), also further the half-peak breadth is calibrated, that is, half-peak breadth after calibration it is default
Between ranging from P/V and P*V, the half-peak breadth then judges that the maximum trough point is available point within the scope of P/V and P*V,
The half-peak breadth then judges that the maximum trough point is Null Spot except P/V and P*V ranges；Wherein, P is default for half-peak breadth
Value, V are half-peak breadth coefficient of alteration；Preferably, the value of the V is ranging from：1.1~1.8.In the present embodiment, the P values are
32, the V values are 1.6, and the preset range of the half-peak breadth is between 32/1.6~32*1.6, that is, the half-peak breadth of trough point
It is available point between 32/1.6~32*1.6.
In the present embodiment, step (70) is still further comprised, is calculated according to the T line positions and the C line positions
The concentration value of detection sample in the colloidal-gold detecting-card；Specifically include following steps：
(71) the gray integration value of the T line positions is calculated according to the T line positions in the gray integration curve,
In, when regarding as no T lines, then the gray integration value of T lines is in terms of zero；
(72) the gray integration value of the C line positions is calculated according to the C line positions in the gray integration curve；
(73) ratio R of the T line positions and the gray integration value of the C line positions is calculated；
(74) ratio R is compared with typical ratio Rs, obtains sample to be tested concentration value.
Specifically, the step (74), is the standard of the standard sample by measuring a certain range of series concentration
Ratio R s, each typical ratio Rs corresponds to a unique standard sample concentration value, to obtain standard curve (such as Fig. 4 institutes
Show) or standard comparison table；Then search corresponding typical ratio Rs according to the ratio R, and by typical ratio Rs
Corresponding standard sample concentration value is as the sample to be tested concentration value.
In the present embodiment, the standard curve can choose different computational methods according to different colloidal-gold detecting-cards, including
Linear fit, logistic fit, floating-point computing method etc., the present embodiment preferably use floating-point computing method.
It is corresponding with detection method above-mentioned, the colloidal gold detecting system based on image procossing that the present invention also provides a kind of,
It includes：
Image pre-processing module obtains the detection image of colloidal-gold detecting-card, and to the area to be identified of the detection image
Domain carries out gray proces, obtains gray level image；
Curve processing module, the gray integration curve for calculating the gray level image；And it is bent to search the gray integration
All trough points of line, and calculate the half-peak breadth of each trough point, peak center point, effective integral area；And search effectively product
Facet accumulates maximum trough point, obtains maximum trough point；
Effective judgement module, if the half-peak breadth of the maximum trough point judges the maximum wave except preset range
Valley point is Null Spot, and it is dead card to regard as the colloidal-gold detecting-card；If the half-peak breadth of the maximum trough point is default
Within the scope of, then the maximum trough point is judged for available point, and it is effectively to block to regard as the colloidal-gold detecting-card；
CT line locating modules carry out judging the maximum trough point according to the position of the peak center point of the maximum trough point
For C line positions or T line positions；If being judged as C line positions, to other significant wave valley points of T lines direction finding, if finding other
The trough point is then regarded as T line positions by significant wave valley point, and no T lines are regarded as if it can not find other significant wave valley points；If sentencing
Break and recognizes the trough point if finding other significant wave valley points then to other significant wave valley points of C lines direction finding for T line positions
It is set to C line positions；No C lines are regarded as if it can not find other significant wave valley points, and it is nothing to regard as the colloidal-gold detecting-card
Effect card.
The curve processing module further comprises curve generation unit and curve computation unit, wherein：
For the curve generation unit by the gray value [X, Y, Z] of the acquisition gray level image, X is the long side of reaction plate
Pixel coordinate, Y are the pixel coordinate of the short side of reaction plate, and Z is the gray value at coordinate position (X, Y)；Using X-coordinate as horizontal seat
Mark carries out gray integration summation to Y coordinate, obtains gray integration curve [X, Zsum] using Y coordinate as ordinate；And further
The gray integration curve is filtered；
The curve computation unit carries out calculating the half-peak breadth of the trough point, peak center point, effective integral area, specifically
's：
Calculate peak center point：It is by centered on trough point, finding inflection point to both sides respectively, calculating inflection point and trough point
Ordinate difference Z1, Z2, and abscissa difference X1, X2 is found according to the ordinate difference, it is in peak to take the central point of X1, X2
Heart point；
Calculate half-peak breadth：It is by doing straight line with X1, X2, the distance of trough point to straight line is peak height, and according to peak height
Calculate half-peak breadth；
Calculate effective integral area：It is the summation by calculating the peak height of all the points between X1, X2, as trough point has
Imitate integral area.
In the present embodiment, still further comprise quantitative analysis module, according to the T line positions and the C line positions into
Row calculates the concentration value of the detection sample in the colloidal-gold detecting-card；Specifically：
The gray integration value of the T line positions is calculated according to the T line positions in the gray integration curve, wherein recognize
When being set to no T lines, then the gray integration value of T lines is in terms of zero；
The gray integration value of the C line positions is calculated according to the C line positions in the gray integration curve；
Calculate the ratio R of the T line positions and the gray integration value of the C line positions；
The ratio R is compared with typical ratio Rs, obtains sample to be tested concentration value.
Preferably, the ratio R is compared with typical ratio Rs, is by measuring a certain range of series concentration
Standard sample typical ratio Rs, each typical ratio Rs corresponds to a unique standard sample concentration value, to obtain
Standard curve or standard comparison table；Then search corresponding typical ratio Rs according to the ratio R, and by the standard ratio
The corresponding standard sample concentration values of value Rs are as the sample to be tested concentration value.
A kind of the colloidal gold detection method and system based on image procossing of the present invention, by the glue for obtaining sample to be tested
The image of body gold detection card, and image compress/stretch, the pretreatments such as cutting processing, then pass through gray integration, filtering
Algorithm obtains the gray integration curve of image, finds out trough point and by parameter preset and calibration factor, according in half-peak breadth, peak
Heart point position, effective integral area carry out C lines, the fixation and recognition of T lines and the calculating of integral area ratio, last to be waited for according to
Calculating is compared with the series standard sample in the integral ratio of the T/C lines of sample, to realize colloidal-gold detecting-card
Qualitative and quantitative analysis, few to the image qualifications of acquisition, C lines, the automatic positioning discriminance analysis of T lines are accurate, and good compatibility is fitted
With range is wide, accuracy of detection is high.
Colloidal-gold detecting-card described in the present embodiment mainly by taking Clenbuterol colloidal-gold detecting-card as an example, is of course also apply to it
His colloidal-gold detecting-card, is not limited, only can be according to different colloidal-gold detecting-card and difference in design parameter setting
Sample is detected to be configured, operation principle is substantially similar, in the contemplated scope of the present invention, and all should be in institute of the present invention
In attached scope of the claims.
It should be noted that each embodiment in this specification is described in a progressive manner, each embodiment weight
Point explanation is all difference from other examples, and the same or similar parts between the embodiments can be referred to each other.
For system embodiments, since it is basically similar to the method embodiment, thus description it is fairly simple, related place referring to
The part of embodiment of the method illustrates.Also, herein, the terms "include", "comprise" or its any other variant meaning
Covering non-exclusive inclusion, so that the process, method, article or equipment including a series of elements includes not only that
A little elements, but also include other elements that are not explicitly listed, or further include for this process, method, article or
The intrinsic element of equipment.In the absence of more restrictions, the element limited by sentence "including a ...", is not arranged
Except there is also other identical elements in the process, method, article or apparatus that includes the element.In addition, this field is general
Logical technical staff is appreciated that realize that all or part of step of above-described embodiment can be completed by hardware, can also pass through
Program instructs the relevant hardware to complete, and the program can be stored in a kind of computer readable storage medium, above-mentioned to carry
To storage medium can be read-only memory, disk or CD etc..
The preferred embodiment of the present invention has shown and described in above description, it should be understood that the present invention is not limited to this paper institutes
The form of disclosure is not to be taken as excluding other embodiments, and can be used for other combinations, modifications, and environments, and energy
Enough in this paper invented the scope of the idea, modifications can be made through the above teachings or related fields of technology or knowledge.And people from this field
The modifications and changes that member is carried out do not depart from the spirit and scope of the present invention, then all should be in the protection of appended claims of the present invention
In range.
Claims (12)
1. a kind of colloidal gold detection method based on image procossing, which is characterized in that include the following steps：
(10) detection image of colloidal-gold detecting-card is obtained, and gray proces are carried out to the region to be identified of the detection image, is obtained
To gray level image；
(20) the gray integration curve of the gray level image is calculated；
(30) all trough points of the gray integration curve are searched, and calculates the half-peak breadth of each trough point, peak center point, have
Imitate integral area；
(40) the maximum trough point of effective integral area is searched, maximum trough point is obtained；
(50) if the half-peak breadth of the maximum trough point judges the maximum trough point for Null Spot except preset range, and
It is dead card to regard as the colloidal-gold detecting-card；If the half-peak breadth of the maximum trough point judges within preset range
The maximum trough point is available point, and it is effectively to block to regard as the colloidal-gold detecting-card, then executes step (60)；
(60) it carries out judging the maximum trough point for C line positions or T according to the position of the peak center point of the maximum trough point
Line position；It, will if finding other significant wave valley points to other significant wave valley points of T lines direction finding if being judged as C line positions
The trough point regards as T line positions, and no T lines are regarded as if it can not find other significant wave valley points；If being judged as T line positions,
To other significant wave valley points of C lines direction finding, the trough point is regarded as into C line positions if finding other significant wave valley points；If
It can not find other significant wave valley points and then regard as no C lines, and it is dead card to regard as the colloidal-gold detecting-card；
Wherein, in the step (30), the half-peak breadth, peak center point, the computational methods of effective integral area are as follows：
(31) peak center point is calculated：
Centered on trough point, respectively to both sides searching inflection point, calculating inflection point and trough point ordinate difference Z1, Z2, and according to
The ordinate difference finds abscissa difference X1, X2, and it is peak center point to take the central point of X1, X2；
(32) half-peak breadth is calculated：
Straight line is done with X1, X2, the distance of trough point to straight line is peak height, and calculates half-peak breadth according to peak height；
(33) effective integral area is calculated：
Calculate the summation of the peak height of all the points between X1, X2, as the effective integral area of trough point.
2. a kind of colloidal gold detection method based on image procossing according to claim 1, it is characterised in that：The step
Suddenly in (10), also further the detection image is stretched or is compressed to the length of reaction plate, and according to the inspection of the reaction plate
Region is surveyed the detection image is cut to obtain the region to be identified of the detection image.
3. a kind of colloidal gold detection method based on image procossing according to claim 2, it is characterised in that：To the inspection
Altimetric image is cut to obtain the region to be identified of the detection image, refers to retaining transverse direction M* centered on the reaction plate
L, the pixel region of longitudinal direction M*L*0.3 is as the region to be identified；Wherein, M is calibration factor, and L is the length of the reaction plate
Degree；The value of the M is ranging from：0.75~0.95.
4. a kind of colloidal gold detection method based on image procossing according to claim 2, it is characterised in that：The step
Suddenly the gray integration curve that the gray level image is calculated in (20), further comprises：
(21) gray value [X, Y, Z] of the gray level image is obtained, wherein X is the pixel coordinate of the long side of the reaction plate, Y
For the pixel coordinate of the short side of the reaction plate, Z is the gray value at coordinate position (X, Y)；
(22) using X-coordinate as abscissa, using Y coordinate as ordinate, gray integration summation is carried out to Y coordinate, obtains gray scale product
Component curve [X, Zsum]；
(23) the gray integration curve is filtered.
5. a kind of colloidal gold detection method based on image procossing according to claim 1, it is characterised in that：The step
Suddenly it carries out calculating peak center point in (31), also further ordinate difference Z1, Z2 is calibrated, that is, is after calibration vertical
Coordinate difference is N*Z1, N*Z2, and carries out searching abscissa difference X1, X2 according to ordinate difference N*Z1, N*Z2 after calibration；
Wherein, N is calibration factor, and the value of N is ranging from：0.75~0.95.
6. a kind of colloidal gold detection method based on image procossing according to claim 1, it is characterised in that：Also further
The half-peak breadth is calibrated, that is, the preset range of the half-peak breadth after calibration is between P/V and P*V, and the half-peak breadth is in P/
Then judge that the maximum trough point is available point within the scope of V and P*V, the half-peak breadth then judges except P/V and P*V ranges
The maximum trough point is Null Spot；Wherein, P is half-peak breadth preset value, and V is half-peak breadth coefficient of alteration；The range of the value of the V
For：1.1~1.8.
7. a kind of colloidal gold detection method based on image procossing according to any one of claims 1 to 6, feature exist
In：Step (70) is still further comprised, carries out calculating the colloidal-gold detecting-card according to the T line positions and the C line positions
In detection sample concentration value；Specifically include following steps：
(71) the gray integration value of the T line positions is calculated according to the T line positions in the gray integration curve, wherein recognize
When being set to no T lines, then the gray integration value of T lines is in terms of zero；
(72) the gray integration value of the C line positions is calculated according to the C line positions in the gray integration curve；
(73) ratio R of the T line positions and the gray integration value of the C line positions is calculated；
(74) ratio R is compared with typical ratio Rs, obtains sample to be tested concentration value.
8. a kind of colloidal gold detection method based on image procossing according to claim 7, it is characterised in that：The step
Suddenly (74) are the typical ratio Rs, each typical ratio Rs of the standard sample by measuring a certain range of series concentration
A corresponding unique standard sample concentration value, to obtain standard curve or standard comparison table；Then according to the ratio R
It carries out searching corresponding typical ratio Rs, and using the corresponding standard sample concentration values of typical ratio Rs as the sample to be tested
Concentration value.
9. a kind of colloidal gold detecting system based on image procossing, which is characterized in that including：
Image pre-processing module, obtain colloidal-gold detecting-card detection image, and to the region to be identified of the detection image into
Row gray proces, obtain gray level image；
Curve processing module, the gray integration curve for calculating the gray level image；And search the gray integration curve
All trough points, and calculate the half-peak breadth of each trough point, peak center point, effective integral area；And search effective integral face
The maximum trough point of product, obtains maximum trough point；
Effective judgement module, if the half-peak breadth of the maximum trough point judges the maximum trough point except preset range
For Null Spot, and it is dead card to regard as the colloidal-gold detecting-card；If the half-peak breadth of the maximum trough point is in preset range
Within, then the maximum trough point is judged for available point, and it is effectively to block to regard as the colloidal-gold detecting-card；
CT line locating modules carry out judging the maximum trough point for C according to the position of the peak center point of the maximum trough point
Line position or T line positions；If being judged as C line positions, to other significant wave valley points of T lines direction finding, if finding other has
The trough point is then regarded as T line positions by effect trough point, and no T lines are regarded as if it can not find other significant wave valley points；If judging
For T line positions, then to other significant wave valley points of C lines direction finding, by the trough point identification if finding other significant wave valley points
For C line positions；No C lines are regarded as if it can not find other significant wave valley points, and it is invalid to regard as the colloidal-gold detecting-card
Card；
Wherein, the curve processing module further comprises curve computation unit：Described in the curve computation unit is calculated
The half-peak breadth of trough point, peak center point, effective integral area；Specifically：
Calculate peak center point：It is by centered on trough point, finding inflection point to both sides respectively, calculating inflection point and the vertical seat of trough point
Difference Z1, Z2 is marked, and abscissa difference X1, X2 is found according to the ordinate difference, it is peak center to take the central point of X1, X2
Point；
Calculate half-peak breadth：It is by doing straight line with X1, X2, the distance of trough point to straight line is peak height, and is calculated according to peak height
Half-peak breadth；
Calculate effective integral area：It is the summation by calculating the peak height of all the points between X1, X2, as effective product of trough point
Facet is accumulated.
10. a kind of colloidal gold detecting system based on image procossing according to claim 9, it is characterised in that：The song
Line processing module further comprises curve generation unit：The gray value that the curve generation unit passes through the acquisition gray level image
[X, Y, Z], X are the pixel coordinate of the long side of reaction plate, and Y is the pixel coordinate of the short side of reaction plate, and Z is coordinate position (X, Y)
The gray value at place；Using X-coordinate as abscissa, using Y coordinate as ordinate, gray integration summation is carried out to Y coordinate, obtains gray scale
Integral curve [X, Zsum]；And further the gray integration curve is filtered.
11. a kind of colloidal gold detecting system based on image procossing according to claim 9, it is characterised in that：Also into one
Step includes quantitative analysis module, according to the T line positions and the C line positions calculate in the colloidal-gold detecting-card
Detect the concentration value of sample；Specifically：
The gray integration value of the T line positions is calculated according to the T line positions in the gray integration curve, wherein regard as
When without T lines, then the gray integration value of T lines is in terms of zero；
The gray integration value of the C line positions is calculated according to the C line positions in the gray integration curve；
Calculate the ratio R of the T line positions and the gray integration value of the C line positions；
The ratio R is compared with typical ratio Rs, obtains sample to be tested concentration value.
12. a kind of colloidal gold detecting system based on image procossing according to claim 11, it is characterised in that：It will be described
Ratio R is compared with typical ratio Rs, is the typical ratio of the standard sample by measuring a certain range of series concentration
Rs, each typical ratio Rs corresponds to a unique standard sample concentration value, to obtain standard curve or standard comparison
Table；Then it is carried out searching corresponding typical ratio Rs according to the ratio R, and the corresponding standard samples of typical ratio Rs is dense
Angle value is as the sample to be tested concentration value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710037576.9A CN106771169B (en) | 2017-01-19 | 2017-01-19 | A kind of colloidal gold detection method and system based on image procossing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710037576.9A CN106771169B (en) | 2017-01-19 | 2017-01-19 | A kind of colloidal gold detection method and system based on image procossing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106771169A CN106771169A (en) | 2017-05-31 |
| CN106771169B true CN106771169B (en) | 2018-08-10 |
Family
ID=58943370
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710037576.9A Active CN106771169B (en) | 2017-01-19 | 2017-01-19 | A kind of colloidal gold detection method and system based on image procossing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106771169B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108830314B (en) * | 2018-06-06 | 2020-12-15 | 西安纸贵互联网科技有限公司 | Semi-quantitative identification method and device for in-vitro diagnostic reagent |
| CN108805072B (en) * | 2018-06-06 | 2021-01-01 | 西安纸贵互联网科技有限公司 | Qualitative identification method and device for in-vitro diagnostic reagent |
| CN109146882A (en) * | 2018-10-26 | 2019-01-04 | 北京纳百生物科技有限公司 | A kind of food safety colloid gold test paper analytical equipment and its application |
| CN111157673A (en) * | 2018-11-08 | 2020-05-15 | 江苏达骏生物科技有限公司 | Quantitative detection method and device |
| CN113112453B (en) * | 2021-03-22 | 2022-03-22 | 深圳市华启生物科技有限公司 | Colloidal gold detection card identification method and system, electronic equipment and storage medium |
| CN113362279A (en) * | 2021-05-13 | 2021-09-07 | 北京理工大学 | Intelligent concentration detection method of immunochromatographic test paper |
| CN113450383B (en) * | 2021-06-29 | 2022-11-01 | 华中科技大学 | Quantitative analysis method, device, equipment and medium of immunochromatographic test paper |
| CN114511522A (en) * | 2022-01-21 | 2022-05-17 | 长沙协大生物科技有限公司 | Automatic reagent judgment method based on fluorescence immunoassay and colloidal gold method |
| CN116109815B (en) * | 2023-04-11 | 2023-07-28 | 深圳市易瑞生物技术股份有限公司 | Positioning method and device for test card calculation area and terminal equipment |
| CN117475246B (en) * | 2023-12-28 | 2024-03-22 | 苏州诚检生物科技有限公司 | Multithreading high-flux colloidal gold detection system and method |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7128266B2 (en) * | 2003-11-13 | 2006-10-31 | Metrologic Instruments. Inc. | Hand-supportable digital imaging-based bar code symbol reader supporting narrow-area and wide-area modes of illumination and image capture |
| JP4235046B2 (en) * | 2003-07-08 | 2009-03-04 | 新日本製鐵株式会社 | Steel plate surface inspection method, system, image processing apparatus, and computer program |
| CN102183510B (en) * | 2011-01-26 | 2012-11-14 | 上海奥普生物医药有限公司 | Colloidal gold detecting method and device based on digital image processing |
| CN103245780B (en) * | 2013-04-27 | 2015-05-27 | 北京福乐云检测科技有限公司 | Colloidal gold anti-counterfeiting tracing cloud platform |
-
2017
- 2017-01-19 CN CN201710037576.9A patent/CN106771169B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106771169A (en) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106771169B (en) | A kind of colloidal gold detection method and system based on image procossing | |
| US8687856B2 (en) | Methods, systems and computer program products for biometric identification by tissue imaging using optical coherence tomography (OCT) | |
| CN103528617B (en) | A kind of cockpit instrument identifies and detection method and device automatically | |
| Gustavsson et al. | A dynamic programming procedure for automated ultrasonic measurement of the carotid artery | |
| US8787623B2 (en) | Methods, systems and computer program products for diagnosing conditions using unique codes generated from a multidimensional image of a sample | |
| CN103900972B (en) | Multi-feature fusion-based meat freshness hyperspectral image visual detection | |
| US20160228008A1 (en) | Image diagnosis device for photographing breast by using matching of tactile image and near-infrared image and method for aquiring breast tissue image | |
| CN107845090A (en) | A kind of silicon chip detection method and silicon chip detection means | |
| CN105510364B (en) | The nondestructive detection system and its detection method of industrial part defect based on X-ray | |
| US11125687B2 (en) | Method and device for identifying fragmented red blood cells, blood cell analyzer and analysis method | |
| CN106526177A (en) | Biomarker detection system and method based on colloidal gold test strip | |
| CN107167594B (en) | Immunochromatographic test strip quantitative detection device and method | |
| CN109406529A (en) | A kind of property regulation method of AOI defect detecting system | |
| CN107678153B (en) | Optical fiber bundle image processing method and device | |
| Rahmat et al. | Automated color classification of urine dipstick image in urine examination | |
| CN103987317B (en) | Measurement equipment and method | |
| CN108320799A (en) | Image analysis and recognition method for lateral flow paper strip disease diagnosis | |
| Zhang et al. | Multispectral heterogeneity detection based on frame accumulation and deep learning | |
| WO2019076267A1 (en) | Optical fibre bundle image analysis method and apparatus | |
| CN206603773U (en) | The burn degree and area detection system measured based on blood erythrocyte concentration | |
| Neuwinger et al. | Computerized semen analysis with sperm tail detection | |
| CN105931221B (en) | A kind of multi-modal magnetic resonance image difference detecting method and device based on image group | |
| CN1903131A (en) | Preset scaling of ultrasonic apparatus type ultrasonic attenuation quantitative analysis tech. | |
| CN114219780A (en) | Deep learning-based method for detecting and segmenting nail fold microcirculation blood vessel | |
| Wang et al. | Nondestructive detecting cracks of preserved eggshell using polarization technology and cluster analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |