Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Geometrical-optics approach to increase the accuracy in LED-based photometers for point-of-care testing

Open Access Open Access

Abstract

A geometrical-optics approach is proposed to increase the accuracy in photometric measurements, using a point-of-care testing (POCT) LED-based sensor. Due to stray-light effects, the measurement accuracy depends on the dimension of the CMOS area, where the radiation is detected. We propose two image processing approaches and evaluate the influence of the sensor area. In addition, we demonstrate that with the same measurement, both absorption coefficient and refractive index can be determined, measuring the beam attenuation and the spot-size enlargement due to ray refraction.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Healthcare is becoming more consumer-focused, and requires screening diagnostic tests, that are portable and easy to operate, to provide results in real time and at the site of patient care, far from traditional laboratories [1]. The point-of-care testing (POCT) market is estimated to have grown 9.3% between 2013 and 2018, and focuses on early detection, prevention and managing of chronic illnesses, as well as on environmental monitoring. A rapid screening of infectious diseases is critical in remote or resource-limited areas, where there is not an easy access to a clinical infrastructure, and a prompt treatment may prevent infections from spreading. Malaria and dengue fever diagnosis in remote locations are now possible with new POCTs [2,3], as well as blood typing [4], and water quality monitoring [5,6]. Many POCT devices are based on colorimetric o spectrophotometric analysis of a sample liquid, such as blood or saliva, using reagents to form a colored compound: the biomarker concentration is proportional to the color intensity [7]. Mobile phones are used for colorimetric tests, by taking a picture of the sample and estimating the substance concentration from the RGB pixel values in the image [8].

First photometers for clinical chemistry date back to 1930, and nowadays, they play a fundamental role in many biochemical, pharmaceutical and microbiological experiments that involve DNA, RNA, protein isolation, and enzyme kinetics [9–11]. Spectrophotometry is routinely used to quantify the hemoglobin concentration in blood [12], and in Enzyme-Linked Immunosorbent Assay (ELISA) tests, to measure proteins and antibodies for many diseases ranging from HIV to cancer [13].

Photometric and spectrophotometric sensors measure the number of photons absorbed after that a light beam passes through the sample [14], giving information on the absorption parameter α at different wavelengths. The transmittance T is related to the absorbance A, also known as optical density (OD), according to the Beer-Lambert law

T=ΦtΦi=eτ=10A,
where Φt and Φi are the intensities transmitted through the sample and the reference liquid (blank), respectively. The optical depth τ = αh depends on the effective lightpath h travelled inside the liquid, so that it is strictly required that the monochromatic light beam is a pencil of rays with a reduced waist. Therefore, the analysis is performed over only a small amount of the liquid, and this could be a shortcoming in the case that the solution is not homogeneous.

Most bench-top spectrophotometers use lasers or white sources and monochromators to select the testing wavelength, and some optical architectures are bulky and expensive. Simplifying the optical design and reducing the overall cost is strategic for POCT devices, and different architectures have been proposed for spectrophotometric analysis of a liquid inside a cuvette, using incoherent light sources and CMOS camera [15–17]. However, the accuracy of the absorbance measurements as function of the sensor area has not been fully investigated yet.

We propose a geometrical-optics approach to increase the measurement accuracy in a low-cost optical photometer, composed of a LED, a lens and a CMOS camera. In general, the use of incoherent light hampers the precision of absorbance measurements, due to light diffraction and stray light effects. We demonstrate that a low-complexity image processing can enhance the accuracy of concentration measurements, evaluating not only the beam attenuation, but also the spot-size changes. We consider different dimensions of the sensor area where the beam intensity is detected and analyze how this parameter affects the measurement accuracy. In addition, the spot-size enlargement, due to light refraction and optical ray displacement in the liquid, allows us also to determine the refractive index.

Refractometry has a large number of biomedical and chemical applications, and it is often used to identify a substance, confirm its purity, or measure its concentration. Generally, the change of refractive index is used to quantify the concentration of a solute, such as sugar, in an aqueous solution (Brix degree). It can also be used in the determination of drug concentration in pharmaceutical industry [18].

The POCT photometer is compact and portable, and the image processing is implemented real-time, using a commercial and low-cost data processor, such as a Raspberry Pi.

In the present paper, we consider only homogeneous liquids, and we demonstrate that the beam spot-size does not change with the solute concentration. However, the same device can be also used to measure the OD of turbid liquids, such as bacterial suspensions [19]. In nephelometric and turbidimetric analysis, the light beam attenuation is due to scattering from particles dispersed in the liquid, the spot-size enlarges, and a different geometrical-optics approach is required to enhance the photometer accuracy [20].

The proposed photometric approach refers to a single wavelength, but it can be used to increase the accuracy of any colorimetric and spectrophotometric POCT devices, using an incoherent source and image sensors.

2. Methods

We use a compact and low-cost spectrophotometric device (WeLab, DNAPhone), with a commercial RGB LED source (Flora RGB Smart NeoPixel version 2, Adafruit) and a 5- Megapixel sensor (OmniVision OV5647, OmniBSI) (Raspberry Pi Camera). The device has a do-it-yourself (DIY) architecture, to enable any user to build a customized low-cost sensor, for different applications [21]. The optical schematics is shown in Fig. 1, and the LED wavelengths and intensities are reported in Table 1. The camera mounts a 2592 × 1944 CMOS image sensor (pixel size 1.4 × 1.4 μm) and 1/4” optics. The lens has focal length f = 3.6 mm and f-number f/2.9. Commercial, low-cost disposable polystyrene cuvette (2xOptical Sarstedt) with h = 10 mm width have been used. The sensor architecture has very basic optical and mechanical alignment, to reduce the fabrication costs [21]; however, the measurement accuracy can be then improved by additional image processing. Furthermore, the CMOS sensor presents different in-system programming features to adjust the image brightness to a desired range by setting exposure time and gain. Therefore, it is possible to perform an image sensor calibration both in manual and automatic modes, and this is highly recommended when the liquids to be examined may present very different colours and densities.

 figure: Fig. 1

Fig. 1 Photometer layout (drawing not in scale).

Download Full Size | PDF

For each measurement, an 8-bit raw image is acquired at the central 540 × 540 pixels, using a single wavelength illumination, and only the corresponding color in the RGB matrix is considered. In the experiments, we use only red-light source to obtain data compatible with conventional OD at 600nm (OD600) measurements; however, the other spectrophotometric acquisitions, using blue and green wavelengths, can be used for additional colorimetric measurements [16–18].

In front of the LED, there is a diffuser and a pinhole with diameter D1 = 4 mm, and we assume that the source can be approximated as a generalized incoherent disk, at a distance d from the lens, with angular aperture θ1 = atan(D1/2d) = 4.52 deg. The LED luminous flux ranges within the solid angle 2θ1 are reported in Table 1.

In the numerical simulations, we assume that the intensity distribution on the CMOS image sensor has a Gaussian profile (Fig. 2 (a) and (b))

Φ(X,Y)=4Φ1πw12e4X2+Y2w12,
where Φ1 is the average intensity and
w1=D1dd'=2d'tanθ1=0.4mm
is the full beam spot-size, that corresponds to the geometrical image of the D1 = 4 mm pinhole (magnification factor M = d’/d = 0.1). We observe that the Gaussian profile model fits quite well with the measurements, but a Lambertian intensity distribution may also be used [22]. The key parameter of the proposed approach is the beam spot-size, that is evaluated using Eq. (3) in a ray-optics model. If X2 + Y2 = w12/4, the radiant intensity of Eq. (2) is e−1 = 0.37 of the maximum. Therefore, we measure the beam spot-size as the sensor area (approximated to a circle) covered by the 8-bit pixels (value range 0 - 255) with values larger than R = 0.37*255~94.

 figure: Fig. 2

Fig. 2 False-color beam intensity distribution on the image sensor in vacuum. (a) Simulated 3D beam pattern. (b) Simulated contour plot (c) Measured contour plot.

Download Full Size | PDF

3. Results

3.1 Absorbance measurements

Figure 2 (c) shows the false-color contour plot of the red-light beam focused on the image sensor (without any cuvette in the lightpath), that is compared with the numerical model of Fig. 2 (b). The raw images of the beam transmitted in vacuum (without cuvette) (Fig. 3 (a)), through a cuvette filled by pure water (Fig. 3 (b)), and cuvettes filled with a food-dye blue-colored water are reported in Fig. 3 (c)-(f), for serial dilution ratios (c) 1:32000, (d) 1:16000, (e) 1:8000, (f) 1:4000, respectively. The axial beam profiles and their normalized values are shown in Fig. 4 (a) and (b), respectively. We observe that the presence of water in the cuvette reduces the maximum intensity and enlarges the beam spot size, due to the ray displacement; this effect will be exploited in the following section to measure the liquid refractive index.

 figure: Fig. 3

Fig. 3 Raw images for dye concentration measurements. (a) without cuvette; (b) pure water (blank); dye-colored water with dilution ratios (c) 1:32000; (d) 1:16000; (e) 1:8000; (f) 1:4000.

Download Full Size | PDF

 figure: Fig. 4

Fig. 4 (a) Axial beam profiles. (b) Axial beam profiles normalized to their own maxima. Data have been smoothed with a moving-average filter only for plot rendering.

Download Full Size | PDF

We also remark that the light beam is attenuated in colored water, but the beam profile does not change, and the average spot size of the corresponding normalized intensity remains practically the same (see Fig. 4(b)). In fact, for small incident angles θθ1, the absorption parameter only slightly affects the ray displacement related to the Snell’s law [23]. Therefore, we can evaluate the intensities transmitted through the sample Φt and through pure water Φi (blank) on the same sensor area, that should be suitably selected. The smallest area that yields accurate OD measurements corresponds to the beam spot-size w1 = 2d’ tanθ1 = 0.4 mm, that is the image of the pinhole of diameter D1. Since some beam misalignments are possible, the circular area under observation has a diameter of 350 pixels (i.e., sensor area corresponding to the beam spot-size, with 0.49 mm diameter). We gradually increase the area dimension to reach an extended circular region with 540-pixel diameter (full area with 0.76 mm diameter). We have verified that larger sensor areas include too much noise due to stray-light effect and therefore they are discarded. Therefore, we restrict our analysis to a circular region on the CMOS region, with a diameter that ranges from N = 350 to N = 540 pixels. The intensity of the pixels outside the selected circular area are all set equal to zero.

According to the Beer-Lambert law of Eq. (1), the absorbance parameter A is evaluated as the ratio between the light intensities transmitted through the sample Φt and through the reference liquid Φi.

We consider two different methods to evaluate the absorbance parameter A, measuring the pixel intensities Φt(n,m) and Φi(n,m) on the CMOS image sensor:

  • I) A is calculated as the ratio between the averaged intensities evaluated over the circular sensor area
    A=log10ΦiΦt=log10n,m=1NΦi(n,m)n,m=1NΦdark(n,m)n,m=1NΦt(n,m)n,m=1NΦdark(n,m),

    and

  • II) using a pixel-by-pixel approach, A’ is computed by summing the ratios between the incident Φi(n,m) and transmitted Φt(n,m) intensities for each pixel
    A'=log10ΦiΦt=log10(1N2n,m=1NΦi(n,m)Φdark(n,n)Φt(n,m)Φdark(n,m)).

    In both cases, we found that the measurement accuracy is increased if we subtract the intensity Φdark(n,m) measured when the LED is off, that is related to the background dark noise.

The absorbance measurements are compared with those obtained with the commercial spectrophotometer BioPhotometer basic (Eppendorf), that is taken as a reference.

Figure 5 shows the OD measurements obtained with the two approaches of Eq. (4) and (5), considering two sensor areas of diameter N = 540 and N = 350 pixels. We observe that both approaches give satisfactory results in the linear range (A, A’ = 0.1 ÷ 1), and that in both cases the accuracy increases considering a smaller CMOS area. In addition, the average error is less than 2% using Eq. (5) in the linear absorbance range.

 figure: Fig. 5

Fig. 5 Absorbance measurements (blue line) and corresponding error (red line) as a function of the CMOS area, where the intensities are measured. Absorbance A is evaluated using Eq. (4) and A’ using Eq. (5). The reference absorbance is plotted with black line.

Download Full Size | PDF

Therefore, we can conclude that the optimal sensor area equates the beam spot-size (diameter N = 350 pixels), i.e. only pixels with normalized values larger than 0.37 should be included in the evaluation of the absorbance parameter.

In addition, the image processing based on the pixel-by-pixel approach of Eq. (5) yields an accuracy of 2%, that is quite remarkable for a low-cost LED-based device.

We have performed a large set of measurements, also using green and blue LEDs, and other coloring dyes, obtaining similar results.

3.2 Refractive index measurements

Using a single measurement, we are able to determine both the absorbance and the liquid refractive index. The refractive index n can be evaluated by measuring the enlargement of the beam spot size, due to the light refraction in the liquid. As it is shown in Fig. 1, the presence of liquid inside the cuvette creates a ray displacement δ, that is related to the liquid refractive index n by the Snell’s law [24]

n=tanθ1tanδ1+tan2δ1+tan2θ1.

For sake of simplicity, we do not consider the refraction effects on the plastic cuvette, that are quite negligible, and we put the air refractive index equal to 1. The ray displacement is evaluated as a function of the enlarged spot-size w2 = 2d’ tanθ2

tanδ=D12(dh)tanθ22h.

Figure 6(a) shows the measured spot-sizes w2 for four different liquids, evaluated as the averaged area where the pixel values (ranging between 1 and 255) are larger than R = 94.

 figure: Fig. 6

Fig. 6 (a) Contour plots of the beam intensity at R = 94 level used to measure the spot-sizes w1 and w2. (b) Refractive index n versus the spot-size variation (w2-w1)/w1. Reference curve obtained with Eqs. (6) and (7) is plotted in full blue line. The linear fitting curve of Eq. (8) is plotted in black dotted line for n0 = 1.4.

Download Full Size | PDF

From these measurements, the corresponding refractive index is evaluated using Eqs. (6) and (7), and the results are reported in Fig. 6(b) and in Table 2.

Tables Icon

Table 2. Spot size vs refractive index

Since θ1 = 4.52 deg, we can use the small-angle approximation and evaluate the refractive index as n~tanδ/tanθ1; in this case, the spot-size variation becomes

Δw=w2w1w1=hdhn1nhdh(n01n0+nn0n02)
that can be evaluated as a linear approximation about the point n0. For the central value n0 = 1.4, the spot-size variation corresponds to −28.08 + 33.43/ R.I.U. (refractive index unit). Using the proposed optical architecture, we are able to analyze a broad range that goes from n = 1.2 to 1.6, but with a limited accuracy that depends on the slope of Eq. (8). If the refractive index range is reduced, the linear approximation furnishes more accurate values.

The refractive index interval and measurement accuracy can be varied by changing the distance d of the LED source from the lens, i.e. the magnification parameter M (the cuvette width h = 10 mm is a standard parameter). A spot-size increase of at least 2 pixels can be detected, that corresponds to about Δw = 1% variation, with an average accuracy Δn = n02 Δw (d-h)/h = 0.03 R.I.U.

4. Summary

Commercial spectrophotometers used in laboratories employ expensive lasers or monochromators and photodetectors. To reduce the costs, portable POCT for colorimetric analysis mount LED sources and CMOS image sensors; however, their accuracy is limited by diffraction and stray-light effects. We propose a geometrical optics approach to evaluate their accuracy and investigate the influence of the dimension of the image sensor area, where the transmitted intensity is detected. We refer to two different real-time image processing strategies, measuring the average intensity on the sensor area, or evaluating the absorbance by a pixel-by-pixel approach. Furthermore, with the same intensity measurements, we can also evaluate the liquid refractive index, measuring the beam spot-size enlargement and identify the liquid.

Funding

MBSMART (A0114E0030, A0114-2017-14829, POR FESR LAZIO 2014 – 2020).

Acknowledgments

The Authors are grateful to Dr. Alessandro Candiani and DNAPhone for their support and assistance. We are also thankful to anonymous Reviewers whose comments helped us to improve the technical quality of the paper.

Disclosures

The authors declare that there are no conflicts of interest related to this article.

References

1. P. Luppa and R. Junker, Point-of-care Testing: Principles and Clinical Applications (Springer, 2018).

2. S. E. McBirney, D. Chen, A. Scholtz, H. Ameri, and A. M. Armani, “Rapid diagnostic for point-of-care malaria screening,” ACS Sens. 3(7), 1264–1270 (2018). [CrossRef]  

3. A. Thiha and F. Ibrahim, “A colorimetric enzyme-linked immunosorbent assay (ELISA) detection platform for a point-of-care dengue detection system on a lab-on-compact-disc,” Sensors (Basel) 15(5), 11431–11441 (2015). [PubMed]   [CrossRef]  

4. A. Thiha, J. Fernandes, S. Pimenta, F. Soares, and G. Minas, “A complete blood typing device for automatic agglutination detection based on absorption spectrophotometry,” IEEE Trans. Instrum. Meas. 64(1), 112–119 (2015). [CrossRef]  

5. S. Sumriddetchkajorn, K. Chaitavon, and Y. Intaravanne, “Mobile-platform based colorimeter for monitoring chlorine concentration in water,” Sens. Actuators B Chem. 191, 561–566 (2014). [CrossRef]  

6. C. Zhao, G. Zhong, D. E. Kim, J. Liu, and X. Liu, “A portable lab-on-a-chip system for gold-nanoparticle-based colorimetric detection of metal ions in water,” Biomicrofluidics 8(5), 052107 (2014). [PubMed]   [CrossRef]  

7. R. S. Khan, Z. Khurshid, and F. Y. I. Asiri, “Advancing point-of-care (PoC) testing using human saliva as liquid biopsy,” in Proceedings Diagnostics (2017).

8. L. Shen, J. A. Hagen, and I. Papautsky, “Point-of-care colorimetric detection with a smartphone,” Lab Chip 12(21), 4240–4243 (2012). [CrossRef]  

9. E. R. Holiday, “Spectrophotometry of proteins,” J. Biochem. 30(10), 1795–1803 (1936). [PubMed]   [CrossRef]  

10. W. S. Hoffman, “Recent advances in photometric clinical chemistry,” Amer. J. clin. Path. 12(9), 449–457 (1942).

11. A. T. Trumbo, E. Schultz, M. G. Borland, and M. E. Pugh, “Applied spectrophotometry: analysis of a biochemical mixture,” Biochem. Mol. Biol. Educ. 41(4), 242–250 (2013). [CrossRef]  

12. E. J. van Kampen and W. G. Zijlstra, “Determination of hemoglobin and its derivatives,” Adv. Clin. Chem. 8, 141–187 (1966). [CrossRef]  

13. R. M. Lequin, “Enzyme immunoassay (EIA)/enzyme-linked immunosorbent assay (ELISA),” Clin. Chem. 51(12), 2415–2418 (2005). [CrossRef]  

14. J. A. Räty, K. -E. Peiponen, and T. Asakura, UV-visible Reflection Spectroscopy of Liquids (Springer-Verlag Berlin Heidelberg 2004).

15. J.-S. Kim, A. H. Kim, H. B. Oh, B. J. Goh, E. S. Lee, J. S. Kim, G. I. Jung, J. Y. Baek, and J. H. Jun, “Simple LED spectrophotometer for analysis of color information,” Biomed. Mater. Eng. 26(s1), S1773–S1780 (2015). [CrossRef]  

16. K. Long, H. Yu, and B. Cunningham, “Smartphone instrument for portable enzyme-linked immunosorbent assays,” Biomed. Opt. Express 5(11), 3792–3806 (2014). [PubMed]   [CrossRef]  

17. C. Zhang, J. P. Kim, M. Creer, J. Yang, and Z. Liu, “A smartphone-based chloridometer for point-of-care diagnostics of cystic fibrosis,” Biosens. Bioelectron. 97, 164–168 (2017). [CrossRef]  

18. M. D. Green, H. Nettey, O. Villalva Rojas, C. Pamanivong, L. Khounsaknalath, M. Grande Ortiz, P. N. Newton, F. M. Fernández, L. Vongsack, and O. Manolin, “Use of refractometry and colorimetry as field methods to rapidly assess antimalarial drug quality,” J. Pharm. Biomed. Anal. 43(1), 105–110 (2007). [CrossRef]  

19. E. Berrocal, D. L. Sedarsky, M. E. Paciaroni, I. V. Meglinski, and M. A. Linne, “Laser light scattering in turbid media Part I: Experimental and simulated results for the spatial intensity distribution,” Opt. Express 15(17), 10649–10665 (2007). [CrossRef]  

20. M. Lucidi, M. Marsan, F. Pudda, M. Pirolo, E. Frangipani, P. Visca, and G. Cincotti, “A geometrical-optics approach to measure bacteria concentrations using a in LED-based photometer,” to be submitted to Biomed. Opt. Express.

21. https://www.we-lab.it/en/welabmaker

22. I. Moreno and C.-C. Sun, “Modeling the radiation pattern of LEDs,” Opt. Express 16(3), 1808–1819 (2008). [CrossRef]  

23. S. A. Kovalenko, “Descartes-Snell law of refraction with absorption,” Semicond. Phys. Quant. Electron. Optoelectron. 4(3) 214–218 (2001)

24. S. Nemoto, “Measurement of the refractive index of liquid using laser beam displacement,” Appl. Opt. 31(31), 6690–6694 (1992). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1
Fig. 1 Photometer layout (drawing not in scale).
Fig. 2
Fig. 2 False-color beam intensity distribution on the image sensor in vacuum. (a) Simulated 3D beam pattern. (b) Simulated contour plot (c) Measured contour plot.
Fig. 3
Fig. 3 Raw images for dye concentration measurements. (a) without cuvette; (b) pure water (blank); dye-colored water with dilution ratios (c) 1:32000; (d) 1:16000; (e) 1:8000; (f) 1:4000.
Fig. 4
Fig. 4 (a) Axial beam profiles. (b) Axial beam profiles normalized to their own maxima. Data have been smoothed with a moving-average filter only for plot rendering.
Fig. 5
Fig. 5 Absorbance measurements (blue line) and corresponding error (red line) as a function of the CMOS area, where the intensities are measured. Absorbance A is evaluated using Eq. (4) and A’ using Eq. (5). The reference absorbance is plotted with black line.
Fig. 6
Fig. 6 (a) Contour plots of the beam intensity at R = 94 level used to measure the spot-sizes w1 and w2. (b) Refractive index n versus the spot-size variation (w2-w1)/w1. Reference curve obtained with Eqs. (6) and (7) is plotted in full blue line. The linear fitting curve of Eq. (8) is plotted in black dotted line for n0 = 1.4.

Tables (2)

Tables Icon

Table 2 Spot size vs refractive index

Equations (8)

Equations on this page are rendered with MathJax. Learn more.

T = Φ t Φ i = e τ = 10 A ,
Φ ( X , Y ) = 4 Φ 1 π w 1 2 e 4 X 2 + Y 2 w 1 2 ,
w 1 = D 1 d d ' = 2 d ' tan θ 1 = 0.4 m m
A = log 10 Φ i Φ t = log 10 n , m = 1 N Φ i ( n , m ) n , m = 1 N Φ d a r k ( n , m ) n , m = 1 N Φ t ( n , m ) n , m = 1 N Φ d a r k ( n , m ) ,
A ' = log 10 Φ i Φ t = log 10 ( 1 N 2 n , m = 1 N Φ i ( n , m ) Φ d a r k ( n , n ) Φ t ( n , m ) Φ d a r k ( n , m ) ) .
n = tan θ 1 tan δ 1 + tan 2 δ 1 + tan 2 θ 1 .
tan δ = D 1 2 ( d h ) tan θ 2 2 h .
Δ w = w 2 w 1 w 1 = h d h n 1 n h d h ( n 0 1 n 0 + n n 0 n 0 2 )
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.