Coupon reader
11023773 · 2021-06-01
Assignee
Inventors
Cpc classification
G06F2218/10
PHYSICS
G06F2218/18
PHYSICS
International classification
Abstract
A method of reading a coupon channel that displays a test section pattern after being exposed to a target substance, the method uses a device having a computer readable memory, digital camera, logic assembly and user interface; providing a pixel target intensity profile; placing the coupon in the device and exposing the coupon channel to a test fluid mixture; automatically using the digital camera to take a digital image of the coupon channel test section after the exposure. The improvement in the method includes finding the contiguous set of pixels from the test section of the coupon channel that best matches the intensity profile of the target pattern representation and determining if this best match set of pixels exceeds a similarity threshold and in response to a best match set of pixels passing the similarity threshold, automatically providing a human perceptible indication that the target substance has been detected.
Claims
1. A method of reading a coupon channel that displays a test section pattern after being exposed to a target substance, said method includes providing a device having a computer readable memory, a digital camera, a logic assembly and a user interface; providing a pixel target intensity profile representation, having m pixels and a shape, and representative of an exposed and at least partially developed coupon channel control section intensity profile; placing said pixel target intensity profile into said memory; placing said coupon in said device and exposing said coupon channel to a test fluid mixture; automatically using said digital camera to take a digital image of said coupon channel test section after said exposure; and a) wherein the improvement in said method includes finding the contiguous set of m pixels from said test section of said coupon channel that best matches said shape of said pixel target intensity profile and determining if this best match set of pixels exceeds a similarity threshold to said shape of said pixel target intensity profile and, in response to a best match contiguous set of pixels passing said similarity threshold, automatically providing a human perceptible indication that the target substance has been detected.
2. The method of claim 1, wherein said coupon channel is hosted on a coupon having a single channel only.
3. The method of claim 1, wherein said coupon channel is hosted on a coupon having multiple channels.
4. The method of claim 1, wherein said contiguous set of pixels that best matches the intensity profile of said target pattern representation is the set of contiguous pixels that has the best least squares match.
5. The method of claim 1, wherein linear bias is subtracted out of each contiguous set of pixels to obtain the best match.
6. The method of claim 1, wherein said finding said contiguous set of m pixels from said test section of said coupon channel that best matches said shape of said pixel target intensity profile includes creating a fitting equation that includes a term that is said target pixel intensity profile times a coefficient and for each set of m contiguous pixels in said test section computing the value of said coefficient that yields the lowest sum of the squares of the difference between each pixel of said contiguous set of m pixel values and each corresponding pixel of said fitting equation, and picking said contiguous set of m pixel values yielding a highest value for said coefficient.
7. The method of claim 6, wherein said fitting equation has another term that is a second coefficient times pixel position, and wherein said second coefficient is set to minimize the effect of a changing pixel magnitude of each set of m pixels, such as may be caused by an uneven illumination of said test section.
8. The method of claim 6, wherein said fitting equation has another term that is a constant, designed to subtract out any constant difference between each set of m pixels and said target pixel intensity profile.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) Exemplary embodiments are illustrated in referenced drawings. It is intended that the embodiments and figures disclosed herein are to be considered illustrative rather than restrictive.
(2)
(3)
(4)
(5)
(6)
DETAILED DESCRIPTION
(7) Referring to
(8) In one embodiment, means are provided to introduce a target pattern representation (an image of a developed coupon) into a computer readable memory 14 of reader 10. This data entry may be performed by placing a developed coupon or a control section into reader 10 and using user interface 12 to command reader 10 to use a digital camera 16 to take a digital photograph of the developed coupon and store it in the computer readable memory 14, properly labeled as a digitized image of a target pattern. In another preferred embodiment, reader 10 is provided with digitized target images already stored. Otherwise digitized target images may be introduced into reader 10 by way of the USB port (not shown) or may be downloaded by way of a cellular network transceiver 18.
(9) After an exposed coupon is introduced into reader 10, the microcontroller 24 controls the digital camera 16 to repeatedly form images of each coupon channel and effectively compare each image with the digitized target image stored in computer readable memory 14.
(10) Coupon manufacturers typically identify a specific time period that should elapse between exposing the coupon to the substance to be tested and reading the results on the coupon. In some cases, results may be visible before this predetermined development time has elapsed and time may be wasted by continuing to wait for the entire pre-determined time period. The ability of the coupon reader 10 to determine when an indicator pattern has been detected and to alert any nearby people that detection has occurred helps to eliminate this wasted time. It will be appreciated that this ability is particularly important when the coupon reader is being used to detect the presence of toxic or harmful compounds in a public safety context, where time may be of the essence.
(11) Many other algorithms, including least squares detection and various linear algorithms are used in alternative embodiments.
(12) Prior to channel analysis, all embodiments require that the coupon image is recognized in the digital image being analyzed. Immediately following coupon image capture, image recognition routines analyze the captured image and identify channel locations based on the detection of characteristic coupon features. In embodiments, the user has entered the coupon type, and this information facilitates the recognition of the channel locations. The processing described below relates to each channel analyzed.
(13) Known optical features of the coupon are pre-loaded into the device before use. Image recognition tools identify the coordinates of these characteristic features for each different coupon type. The pixel coordinates of these features provide reference points and allow each coupon image to be overlaid by a coordinate map which can then be used to locate the assay channels.
(14) Upon positive recognition and satisfaction of location criteria, coupon analysis is allowed to continue. Upon detection of an abnormal situation, such as the coupon being inserted 90 or 180 degrees from the correct orientation, the user is alerted to the type of error. Once the error is remedied and location criteria are satisfied, sample analysis procedures may proceed.
(15) In one embodiment, each time a reader 10 processes a coupon, the results are time, position and device code stamped, with position being provided by a GPS reader 20, time being provided by a clock 22, and device being provided by device firmware, all controlled by microcontroller 24 and stored in computer readable memory 14. In a preferred embodiment, this information is uploaded periodically via cell network link 28 using cell network transceiver 18, to a remote station 30. In one embodiment, remote station 30 includes a display and may be commanded to display a map of an area, showing the locations, times, results and device codes of each coupon reading.
(16) Referring to
(17) Referring now to
(18) In one embodiment, the light diffusing assembly comprises side walls 36 which form a rectangular cross-section approximately 8 cm on a side. The side walls 36 meet the coupon surface 44 at angles that maximize uniformity of illumination over the coupon surface and compensate for the light source's off-center location. In a preferred embodiment, the included angle 46 between the side walls 36 and the coupon surface 44 is:
(19) Face 1=70 degrees
(20) Face 2=70 degrees
(21) Face 3=70 degrees
(22) Face 4=70 degrees
(23) In one embodiment, the distance between the coupon surface 44 and the negative axicon lens 40 is 15 cm.
(24) Diffuse reflection may be obtained from the reflective interior surfaces 38 by covering them with a flat or gloss white paint, or by constructing the structure using a white polymer. However, commercially available white polymers may absorb a significant amount of the optical power impinging on them, and this may need to be taken into consideration.
(25) Optionally, the reflective interior surfaces 38 may be covered with a colored coating rather than a white coating. A colored coating may be desirable in circumstances where a specific light spectrum distribution is preferred for illuminating certain types of coupon. In a further option, the reflective interior surfaces may be covered by a hybrid coating that provides a mixture of specular and diffuse reflection.
(26) Axicons are special lenses that have at least one active surface that is conic in shape. They will typically take a point source of light and transform it into either a line or circle of light. A negative axicon lens is used for coupon reader 10. Typical commercially available axicon lenses are designed to form a specific geometric shape, and for that reason are considered “positive axicons.” The negative axicon, in contrast, takes a point of light and transforms it into an annular fan of light. Instead of a centrally located solid cone of material as in a positive axicon, the negative axicon has a conic pocket into which source light expands.
(27) Refraction at the conic top surface and at the optionally planar lower surface cause rays from a point source to be distributed over a range of angles determined by the conic profile, the distance of the source from the lens, the lens thickness, and the refractive index of the lens material. Angular emission of the lens can be tuned so that no light hits the coupon surface 44 without first reflecting off the reflective interior surfaces 38. That is, in a preferred embodiment, no light passing through the negative axicon lens 40 falls directly on the coupon surface 44, particularly at the spot immediately below the lens where a hotspot would normally occur. In this way, all light impinging on the coupon surface 44 is diffused light re-emitted or reflected by the reflective interior surfaces 38 and the uniformity of illumination over the coupon is greatly improved, even if the coupon is in close proximity to the light source and lens.
(28) It is not necessary that the conic profile of the negative axicon lens 40 be a simple cone. By adjusting the surface profile of the lens, the distribution of emitted optical power can be modified in an angular sense. Since the objective of the negative axicon lens and light diffusing assembly is to provide diffused illumination, minor errors in the surface profiles are not critical. In one embodiment, the lenses are manufactured on a CNC lathe, followed by vapor polishing to remove macroscopic tooling marks.
(29) In one embodiment, a negative axicon lens with planar output surface for supplying optical power to the reflective interior surfaces is constructed of clear polycarbonate. In one embodiment, the negative axicon lens has the following dimensions: the overall height of the lens may be 3.5 mm, the diameter may be 10 mm, and the maximum depth of the cone may be 2.5 mm.
(30) In one embodiment, the cone profile is described as:
Radius=1.98543z−0.18142z.sup.2+0.04361z.sup.3 (where z is the distance in mm from the cone's vertex).
(31) In one embodiment, nonuniformity of illumination across the coupon surface is compensated for by color-based ratiometric analysis of pixel data. Most color cameras provide at least three color sub-pixels per color pixel; typically red, green, and blue sub-pixels. It has been found that light intensity variations over a coupon due to geometric lighting issues are very similar in the three color sub-pixels. Most lateral flow immunoassay coupons use colloidal gold as the line colorant and its peak absorption is typically in the wavelength range of 540 nm to 580 nm, which is principally apparent in the camera pixels sensitive to green and red. It has been discovered that by using the color pixels to form “ratiometric” or, termed slightly differently, “color-ratio” pixels, that is a pixel having an intensity corresponding to the ratio of two sub-pixels, an unexpected and remarkable reduction in background nonuniformity is seen. For the coupons having the common colloidal gold colorant, the red and green sub-pixels are typically most useful in this regard but creating color-ratio pixels using other color sub-pixels may be preferred for other colorants.
(32)
(33) In one embodiment, a new variable is created when the green pixel intensities are divided by the corresponding red pixel intensities, to create green-to-red color-ratio pixels and the resultant curve is normalized to an arbitrary value of 1.0 (or 100%) at some arbitrary pixel location distant from the test and control lines. The resulting color-ratio pixel values are also shown in
(34) The fluid dynamic and chemical processes that create the test and control lines are typically identical or very similar. For most coupons, even those with multiple identification channels, all test and control lines produce pixel signal responses that are similar in shape and substantively differ only in overall intensity. Therefore, the relative concentration of a targeted material may be estimated by simply comparing the maximum deviation of the color-ratio pixel signal from background and comparing it to the maximum deviation of the fully-formed color-ratio pixel control line. However, this may not be the most suitable method for test lines that have a small peak intensity because there may be too high of a signal-to-noise ratio to determine the peak value accurately. In addition, the test and target lines are of finite width and vary in intensity over the line's width. Therefore, the quantity of target material in the sample is more appropriately proportional to the integrated total response over the entire line width. Using the integrated response over the entire line width will also serve to smooth out pixel signal noise. For example, if a signal is averaged over 100 pixels, random background noise as a fraction of the mean signal value will typically be reduced by a factor of 10 times compared to the noise associated with a single pixel measurement process.
(35) In one embodiment, the following mathematical strategy and corresponding algorithm automatically provide a measure of the test and control line integrated areas, and the degree of match between the shapes or “intensity profiles” of the two lines. This technique yields the pixel (or stated slightly differently, “coupon”) location of the test section intensity profile that best matches the reference pattern. That is, the set of contiguous pixels in the test section having an intensity profile that most closely matches that of the reference pattern. The control and test line locations may indicate whether the assay coupon was correctly assembled or not, or if it was subsequently mishandled by the user.
(36) The coupon typically includes a control line, which develops even in the absence of the target substance to indicate that the coupon is functioning properly and that the correct procedures were used to expose and develop the coupon. In one embodiment, the control line is used as the source of reference information that permits the coupon reader to determine whether a test line is present or absent. In other embodiments, a different source of reference information may be used, such as a pre-printed line on the coupon, data from a different coupon, or an abstract representation of a developed control line, showing the expected result when the tested for substance is present.
(37) Referring to
Y.sub.f(i)=a+b×i+c×Y*(i)
(38) In this equation, “a” is the average baseline value, “b” is the rate of baseline change per pixel, and “c(p)” is a measure of the goodness-of-fit of experimental data over a set of pixels of equal width to the reference signature, Y*(i), and starting at pixel index “p”. As noted, in the example of
η=Σ.sub.i=.sup.p+m[Y(i)−Y.sub.f(i−p)].sup.2
(39) Where p=the starting pixel position and m (as noted above) is the number of values in the Y*(i) set of values and i=0, is the starting pixel position of Y.sub.f(i). The values of a(p), b(p) and c(p) yielding the minimum value of η are found (in the manner described below) repeatedly, for each value of “p” over the test section of the coupon. In the example of
(40) Skilled persons will recognize that solving for the c(p) that minimizes η, for each value of p, and then thresholding the largest c(p), is more likely to yield a detection in the instance in which a section of the coupon not only has some intensity (that is, contrast or difference relative to a background noise level, which may be shown as darkness against a white background), but also where the portion tested has a similar shape of intensity profile to Y*(i). As shown in
(41) Also, the value of c(p).sub.max is not directly affected by position of p, within the pixel linear array, so that c(p).sub.max may occur at some other location than the value of p that would be expected based on information from the coupon manufacturer. This has a benefit, because coupons sometimes have flaws in the position of the prospective test line, and it is accordingly advantageous to “look” for a test line in a wider area than only where it would be expected based on the intended position. Moreover, when c(p).sub.max is found, the pixel position “p” may yield valuable information concerning the state of the coupon that has been tested.
(42) The procedure determining the values of the parameters “a(p)”, “b(p)” and “c(p)” to minimize are well-known in statistical mathematics. The procedure is performed over the test pixel section to find a value of “c(p)” for each possible pixel starting position “p”. When this procedure is performed, it is found that the best-fit value of “c” is given by the following equation:
C=K.sub.1×S.sub.iyK.sub.2×S.sub.y−K.sub.3×S.sub.ys
(43) The values K.sub.1, K.sub.2, and K.sub.3 are constants that are derived based on the reference signature Y*(i), and are easily calculated by someone knowledgeable in statistical mathematics. The other equation factors are calculated for each sectional pixel array of Y(p) through Y(p+m) to be examined as follows:
S.sub.iy=Σ.sub.i=p.sup.p+mi×Y(i) . . . S.sub.y=Σ.sub.i=p.sup.p+mY(i) . . . S.sub.ys=Σ.sub.i=p.sup.p+mY(i)×Y*(i−p)
(44) The three summations are easily and quickly calculated and the constants K.sub.1, K.sub.2, and K.sub.3 are only calculated once: After that they may be stored as look-up values for the channel. If the pixel data presented for comparison is the reference pixel array Y*.sub.i itself, the value of “c(p)” will be exactly 1.0 as discussed previously.
(45) As previously discussed, the test line, or pattern, in most cases is essentially identical in shape (that is, in cross-sectional variation in intensity) to the control line, but less intense overall. The value of “c” will have a maximum value when the signal data most closely matches the control line shape Y*(i). This allows identifying the location of each test line center point while the value of “c” represents the relative integral size of the discovered test line as compared to the control line.
(46) This method is favored for several reasons:
(47) Small variations in the baseline are common due to lighting variations, shadows, and residual color in the coupon substrate. This method automatically estimates and corrects for a non-constant baseline. The value of c(p) is an unbiased best estimate of the baseline-corrected integral value.
(48) Due to the use of all data points associated with a line, signal-to-noise is improved and a wide dynamic range is realized of as much as 100 to 1. That is, a test line with a value of c(p).sub.max on the order of 0.01 may be detected.
(49) The value of c(p).sub.max reflects the target substance concentration, which allows development of a calibration curve that quantitatively correlates the targeted material's concentration with the value of c(p).sub.max.
(50) The method allows a quantitative determination that a test line's response is statistically above background noise.
(51) The method provides the lines' peak value location, which is useful in identification of manufacturing or user errors.
(52) While the method can be applied to raw data as well as ratiometric data, the integral value from ratiometric data provides a measure of signal strength that is independent of absolute light intensity or light exposure time as long as the color channels all change proportionately the same when these factors are varied. This is a reasonable assumption for changes in exposure time or temperature, since both the control line and test line will see the same environmental shifts.
(53) In a preferred embodiment, for a multichannel coupon the three constants K.sub.1, K.sub.2, and K.sub.3 are empirically determined for each individual channel by running one or more assays and calculating these least-squares constants for each channel's control line. This may be of particular value if there are differences in channel construction or if assay reagents are significantly different among the channels.
(54) In addition, in a multi-channel coupon the channel control lines may be compared one against another and these relative “c” values stored in the device's memory. During later use, the control lines can be compared to each other to see if any have deviated from the expected relative “c” values. The device may then issue a warning to the user that there may be a problem with one or more channels.
(55) In a preferred embodiment, digital images of the coupon are repeatedly taken over a period of time, starting shortly after addition of the sample and well before the time period suggested by the manufacturer has passed. These images are processed sequentially to determine if a test line has developed. As soon as detection positive result is found, a human perceivable signal is given. Depending on the coupon type, it is possible to gather useful information as early as 1 minute after development has been initiated. Some explosives coupons have a peak color at 2 minutes after application of a test fluid, and noticeably fade after 3 minutes. Accordingly, the time after exposure at which a first digital image is formed may be as little as 10 seconds or even 1 second. The techniques described herein have been shown to reduce the time to detection for bio-coupons to as little as 5 minutes, a third of the time period that the manufacturer recommends waiting, before evaluating an exposed coupon.
(56) Accordingly, in embodiments the taking of digital imagery begins well before the manufacturer's development time has passed. In various embodiments, the taking of digital imagery begins at least 4, 3, 2 and 1 minute prior to the passing of the recommended time to develop, after exposure to said test fluid mixture. Described in different terms, in various embodiments a first digital image is formed no more than 1, 2, 3 and 4 minutes after the coupon channel is exposed to said test fluid mixture.
(57) In a further inventive feature, digital imagery is also repeatedly taken of the control section of each channel, beginning shortly after exposure to a test fluid mixture. This permits early detection of a coupon that has been ineffectively exposed to the test fluid mixture. As soon as ineffective exposure is determined, the test can be terminated and a new test begun, thereby greatly decreasing the time before an accurate reading is determined.
(58) In a further embodiment, the sequence of images are integrated together, to enhance the signal to noise ratio and, in some instances, provide earlier detection than would otherwise be possible. In one variant, time integration is performed on the raw pixels, but in another least squares fitting as described above is applied to the data from each digital image and the multiple results are then integrated over time.
(59) The invention may be embodied in other specific forms besides and beyond those described herein. The foregoing embodiments are therefore to be considered in all respects illustrative rather than limiting, and the scope of the invention is defined and limited only by the appended claims and their equivalents, rather than by the foregoing description.