OSA's Digital Library

Virtual Journal for Biomedical Optics

Virtual Journal for Biomedical Optics

| EXPLORING THE INTERFACE OF LIGHT AND BIOMEDICINE

  • Editors: Andrew Dunn and Anthony Durkin
  • Vol. 6, Iss. 8 — Aug. 26, 2011

Optics InfoBase > Virtual Journal for Biomedical Optics > Volume 6 > Issue 8 > Page 3589

« Show journal navigation

Automated imaging, identification, and counting of similar cells from digital hologram reconstructions

Mona Mihailescu, Mihaela Scarlat, Alexandru Gheorghiu, Julia Costescu, Mihai Kusko, Irina Alexandra Paun, and Eugen Scarlat  »View Author Affiliations


Applied Optics, Vol. 50, Issue 20, pp. 3589-3597 (2011)
http://dx.doi.org/10.1364/AO.50.003589


View Full Text Article

Acrobat PDF (1227 KB)


Advanced Search

Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

This paper presents our method, which simultaneously combines automatic imaging, identification, and counting with the acquisition of morphological information for at least 1000 blood cells from several three-dimensional images of the same sample. We started with seeking parameters to differentiate between red blood cells that are similar but different with respect to their development stage, i.e., mature or immature. We highlight that these cells have different diffractive patterns with complementary central intensity distribution in a given plane along the propagation axis. We use the Fresnel approximation to simulate propagation through cells modeled as spheroid-shaped phase objects and to find the cell property that has the dominant influence on this behavior. Starting with images obtained in the reconstruction step of the digital holographic microscopy technique, we developed a code for automated simultaneous individual cell image separation, identification, and counting, even when the cells are partially overlapped on a slide, and accurate measuring of their morphological features. To find the centroids of each cell, we propose a method based on analytical functions applied at threshold intervals. Our procedure separates the mature from the immature red blood cells and from the white blood cells through a decision based on gradient and radius values.

© 2011 Optical Society of America

1. Introduction

Red blood cell imaging, identification, and counting are the first steps in a rigorous analysis to establish the properties of similar cells exhibiting different development stages, i.e., mature or immature. Automated devices developed for blood cell count, such as flow cytometers and impedance counters, have made work in hospital laboratories easier, offering a fast method for obtaining accurate results. However, it is well established that pathologists can obtain useful information about cell morphology by visualizing the blood smear on a slide, an important diagnostic tool in detecting many diseases. In these procedures, since many cells are transparent in bright-field microscopy, they must use chromatographic agents. For this reason, the procedures remain time-consuming; for example, in the case of immature red blood cells, sample preparation and analysis may take several hours.

In 1865, W. H. Erb [1

E. B. Krumbhaar, “Reticulosis-increased percentage in the peripheral blood,” in Hematology: Landmark Papers of the Twentieth Century , M. A. Lichtman, L. A. Boxer, E. Henderson, and J. L. Spivak, eds. (Academic, 2000), p. 805.

] described immature red blood cells (IRBCs), also named reticulocytes in their last stage, as transitional forms of mature red blood cells (MRBCs). The reticulocyte is formed in the bone marrow and finishes its development in the peripheral circulation [2

L. M. Lowenstein, “The Mammalian reticulocyte,” in International Review of Cytology , G. H. Bourne and J. F. Danielli, eds. (Academic, 1959), Vol. 8. [CrossRef] [PubMed]

]. Reticulocyte maturation into erythrocytes is the final step of erythropoiesis, which occurs in circulating blood [3

L. Blanc, A. De Gassart, C. Géminard, P. Bette-Bobillo, and M. Vidal, “Exosome release by reticulocytes—an integral part of the red blood cell differentiation system,” Blood Cells Mol. Dis. 35, 21–26 (2005). [CrossRef] [PubMed]

]. Many biochemical and physical studies or complicated analyses have been conducted to find their properties and their role in the detection, evolution, or treatment of different diseases [4

J. G. Stephens, “Surface and fragility differences between mature and immature red cells,” J. Physiol. 99, 30–48 (1940). [PubMed]

, 5

A. K. Percy, E. Schmell, B. J. Earles, and W. J. Lennarz, “Phospholipid biosynthesis in the membranes of immature and mature red blood cells,” Biochemistry 12, 2456–2461 (1973). [CrossRef] [PubMed]

, 6

T. Lecklin, A. Tuominen, and M. Nikinmaa, “The adrenergic volume changes of immature and mature rainbow trout (Oncorhynchus mykiss) erythrocytes,” J. Exp. Biol. 203, 3025–3031 (2000). [PubMed]

]. Morphological changes in red blood cells can be induced by a variety of factors, including the lethal human malaria parasite [7

Y. M. Serebrennikova, J. Patel, W. K. Milhous, and L. H. Garcıa-Rubio, “Quantitative analysis of morphological alterations in Plasmodium falciparum infected red blood cells through theoretical interpretation of spectral measurements,” J. Therm. Biol. 265, 493–500 (2010). [CrossRef]

]; during the development of hemoglobinemia following infection with Leptospira interrorgans serovar bullum [8

J. C. Thompson and B. W. Manktelow, “Pathogenesis and red blood cell destruction in hemoglobinemic leptospirosis,” J. Comp. Pathol. 96, 529–540 (1986). [CrossRef] [PubMed]

]; during blood heating [9

S. Christel and C. Little, “Morphological changes during heating of erythrocytes from stored human blood,” J. Therm. Biol. 9, 221–228 (1984). [CrossRef]

]; and after UV irradiation [10

I. A. Kabat, W. Leyko, B. Kwiatkowski, and I. Zakrzewska, “Osmotic properties and morphological changes of submicroscopic surface structures of mammalian red blood cells subjected to UV-irradiation in vitro ,” Zentralbl. Bakteriol. Orig. B 159, 88–94 (1974). [PubMed]

].

Since IRBCs contain a residual ribosomal RNA, when they are combined with dye agents, they permit fluorescent analysis in laser-based counters or bright-field microscopy analysis on a smear. Reticulocyte maturation into erythrocytes is achieved in about 48 h for in vivo conditions, when the IRBC loses its nucleus, changes its shape from oblate spheroid to that of a typical biconcave disk, and minimizes its area and volume. The IRBC percentage in the whole blood has a prominent role in the initial classifying of anemia and in determining the efficiency of chemotherapy treatments. The behavior of the IRBC and MRBC in oxygenation/ deoxygenation cycles gives information that is useful in sickle cell anemia [11

A. J. McGoron, C. H. Joiner, M. B. Palascak, W. J. Claussen, and R. S. Franco, “Dehydration of mature and immature sickle red blood cells during fast oxygenation/deoxygenation cycles: role of KCl cotransport and extracellular calcium,” Blood 95, 2164–2168 (2000). [PubMed]

].

Mohandas et al. [12

N. Mochandas, M. R. Clark, M. S. Jacobs, and S. B. Shohet, “Analysis of factors regulating erythrocyte deformability,” J. Clin. Invest. 66, 563–573 (1980). [CrossRef]

] performed one of the first optical analyses of MRBCs using their diffraction pattern by monitoring two points using a photomultiplier. From these intensity values, they deduced the dependence of erythrocyte deformability versus shear stress. In recent years, many analyses based on optical techniques centered on MRBC behavior in different conditions, or diseases have been reported [13

M. Mir, H. Ding, Z. Wang, J. Reedy, K. Tangella, and G. Popescu, “Blood screening using diffraction phase cytometry,” J. Biomed. Opt. 15, 027016 (2010). [CrossRef] [PubMed]

, 14

S. Rancourt-Grenier, M.-T. Wei, J.-J. Bai, A. Chiou, P. P. Bareil, P.-L. Duval, and Y. Sheng, “Dynamic deformation of red blood cell in dual-trap optical tweezers,” Opt. Express 18, 10462–10472 (2010). [CrossRef] [PubMed]

, 15

C. Fang-Yen, S. Oh, Y. Park, W. Choi, S. Song, H. S. Seung, R. R. Dasari, and M. S. Feld, “Imaging voltage-dependent cell motions with heterodyne Mach–Zehnder phase microscopy,” Opt. Lett. 32, 1572–1574 (2007). [CrossRef] [PubMed]

, 16

B.-W. Yang and Z. Li, “Measuring microinteractions between coagulating red blood cells using optical tweezers,” Biomed. Opt. Express 1, 1217–1224 (2010). [CrossRef]

, 17

Y. K. Park, C. A. Best, K. Badizadegan, R. R. Dasari, M. S. Feld, T. Kuriabova, M. L. Henle, A. J. Levine, and G. Popescu, “Measurement of red blood cell mechanics during morphological changes,” Proc. Natl. Acad. Sci. USA 107, 6731–6736 (2010). [CrossRef] [PubMed]

].

In this paper, we use a holographic method combined with image processing for IRBC and MRBC identification, counting, and morphological properties measurement in a blood smear on a slide without chromatographic agents. The digital holographic microscopy (DHM) technique offers nanometric resolution along the propagation axis and is used in many laboratories for MRBC analysis [18

Y.-S. Choi and S. Lee, “Three-dimensional volumetric measurement of red blood cell motion using digital holographic microscopy,” Appl. Opt. 48, 2983–2990 (2009). [CrossRef] [PubMed]

, 19

B. Rappaz, A. Barbul, A. Hoffmann, D. Boss, R. Korenstein, C. Depeursinge, P. J. Magistretti, and P. Marquet, “Spatial analysis of erythrocyte membrane fluctuations by digital holographic microscopy,” Blood Cells Mol. Dis. 42, 228–232 (2009). [CrossRef] [PubMed]

]. The off-axis DHM technique has the major advantage of reconstructing the MRBC details positioned at various distances along the propagation axis using a single hologram. In recent years, few methods of enhancing the lateral resolution have been reported [20

M. Paturzo, F. Merola, S. Grilli, S. De Nicola, A. Finizio, and P. Ferraro, “Super-resolution in digital holography by a two-dimensional dynamic phase grating,” Opt. Express 16, 17107–17118 (2008). [CrossRef] [PubMed]

, 21

V. Mico, Z. Zalevsky, C. Ferreira, and J. García, “Super-resolution digital holographic microscopy for three-dimensional samples,” Opt. Express 16, 19260–19270 (2008). [CrossRef]

].

Using DHM, we record (on a CCD camera) the interference pattern of the reference and object beams (the last one also carries the diffraction pattern of the object). First, we record the diffraction pattern from IRBCs and MRBCs situated in the same frame. Then, we highlight that, for a given distance between the sample and the objective, their diffraction patterns have complementary central intensity distribution. Next, we develop ellipsoidal models for IRBC and MRBC shapes to use in a propagation simulation analysis in order to establish which of the cell properties has the dominant influence on this behavior.

The reconstructed images from DHM are processed in order to separate edges of the tangential cells, to distinguish IRBCs from the MRBCs, to count automatically, and to further obtain information on their morphological parameters. For accurate numerical values, we remove the background using successive filters. We propose a method based on nonlinear functions applied at threshold intervals and distance transform aiming to separately counting partially overlapping cells. The histograms obtained automatically after simultaneously processing more than six reconstructed images are sub sequently compared with those obtained using the impedance counter.

2. Experimental Procedure

Using Nikon equipment (Eclipse Ti-U, DS-Fi1 5 megapixel CCD camera, and dedicated software, NIS Elements, Nikon Corp., Japan), the experiment includes analyses in bright-field microscopy of the samples with a chromatographic agent and in differential interference contrast (DIC) configuration without a chromatographic agent.

Samples without a chromatographic agent are also studied in a DHM setup based on the Mach–Zehnder interferometer, with the investigated object placed in one arm. Identical microscope objectives are used in the reference and object arms to obtain waves with the same wavefront curvature in the CCD plane. A single recorded hologram (an intensity image) is enough to reconstruct the phase shift introduced by the transparent sample in the optical path. To scan the entire area of the sample, a computer-controlled two-dimensional (2D) motion was performed in a plane perpendicular to the propagation axis. The 40× objective that collects images of the samples allows a 0.9 μm transversal resolution. Two linear polarizers are used to adjust light intensity, and a neutral density filter in the reference arm ensures proper visibility in the holograms recorded by the CCD (Pike F421, Allied Vision Technologies, Germany, with Kodak sensor 2048×2048 pixels, pixel pitch Δx=7.4 μm, acquisition rate of 16 fps at full resolution). Starting from this pitch, we have chosen the minimum interfringe at 9Δx, which corresponds to a spatial frequency of 1/(9Δx), which is in accordance with the maximum spatial frequency established by the Nyquist–Shannon sampling theorem of 1/(2Δx). In this case, the angle between the reference and the object beams is around 0.5 deg.

Certified phlebotomists collected blood from patients in a hospital, respecting all the appropriate hygienic and legal regulations. The blood was then stored in coated containers at 2–3 °C, following the standard procedures of in-vitro preservation. Each container was marked uniquely, and all patient personal data were removed. The content of each container was then split into three identical subsamples, for subsequent comparison: one for automatic counts, the second for chromatographic investigations on a slide, and the last for investigations without a chromatographic agent on a slide. A complete blood count for routine analysis was then performed on each first subsample, using the impedance counter in the hospital hematology laboratory (Multisizer 4 Coulter Counter from Beckman Coulter, Inc., USA). The blood from the second and third subsamples was smeared on a microscope slide. A cover slip was placed on top of the blood drop for samples without a chromatographic agent to avoid the cells suffering morphological changes. In this case, the recorded images in the DHM setup and DIC configuration were taken in only few minutes to avoid alterations and drying.

First, we study the blood smear on the microscopic slide using diluted dye agents in bright-field microscopy and the same region in DHM, to visualize which cells are IRBCs or MRBCs. Figure 1 shows the same region of one sample visualized in bright-field microscopy and DHM. In Fig. 1a, the IRBC appears discolored and slightly greater in diameter. In Fig. 1b, the same cell has a different diffractive pattern than neighboring cells, and in Fig. 1c, reconstructed in DHM, it appears as being mountain- shaped, as opposed to the crater-shaped MRBCs.

To avoid the time-consuming procedure with chromatographic agents, we use the third subsample (without a chromatographic agent) in the DHM setup [see Figs. 2a, 2c]. We present images from two regions with at least one IRBC in each, many MRBCs, and one white blood cell [the largest one in Figs. 2c, 2d] being differentiated by its dimension. One can see that the diffractive patterns are different for these three types of investigated cells. Only in the first studies, we observe the same slide in DIC configuration [see Fig. 2b, 2d] to confirm our identification.

Figure 3 presents the holograms recorded when the distance between the sample and the objective is modified (with an actuator step resolution of 0.5 μm). One can observe the differences in the diffraction pattern of these two cell types at different distances between the sample and the microscope objective, which differentiate the IRBC from the MRBC in the preprocessing step. In order to better observe these differences, we record the same region with the reference beam blocked.

The images from Fig. 4 were recorded at different distances between the sample and the objective lens, but without superposing the reference beam onto the object beam. Here, the central maxima and minima in MRBC and IRBC diffraction patterns are more visible and are complementary in the studied interval. The reference distance of z=0 μm is chosen in the point where the reconstructed object amplitude has minimum diameter.

3. Diffraction Pattern from Modeled IRBCs and MRBCs

Starting from these experimental observations concerning MRBC and IRBC diffractive pattern behavior, we create models for their shape and refractive index. The phase shift introduced by them is proportional to their thickness z(x,y), thus with the length traveled by the incident monochromatic plane wave within the medium of the sample in every point (x,y),
ΔΦ(x,y)= 2πλ[(n(x,y)− nS)z(x,y)],
(1)
where n is the refractive index of the cell, ns is the refractive index of the surrounding medium, and λ is the wavelength of the incident plane wave. In the simulation scheme, we consider the cells to be ellipsoids. At every point, the cell thickness is linked with the coordinates (x,y) in the projection plane through the equation of the centered ellipsoid in the Cartesian coordinate system:
x2 a2+ y2 b2+ z2 c2=1,
(2)
where a and b are the equatorial radii and c is the polar radius. IRBCs are simple oblate spheroids with circular projection aI= bI; their measured values were in the range 7.9–9.4 μm. MRBCs are biconcave oblate spheroids with circular projection aM= bM; their measured values were in the range 5.8–7.7 μm. The subscript letters M and I stand for MRBC and IRBC, respectively.

For the concave shape (which is cropped), different models were considered: semisphere, oblate spheroid, cylindrical or truncated cone, in accordance with experimental images. Similar ellipsoidal models have been introduced in the literature for the study of red blood cell deformability in ektacytometry [22

G. J. Streekstra1, J. G. G. Dobbe, and A. G. Hoekstra, “Quantification of the fraction poorly deformable red blood cells using ektacytometry,” Opt. Express 18, 14173–14182 (2010). [CrossRef]

] or for interpreting their spectral response [7

Y. M. Serebrennikova, J. Patel, W. K. Milhous, and L. H. Garcıa-Rubio, “Quantitative analysis of morphological alterations in Plasmodium falciparum infected red blood cells through theoretical interpretation of spectral measurements,” J. Therm. Biol. 265, 493–500 (2010). [CrossRef]

].

All these equations become matrices, implemented in MATLAB, which digitally describe cell transmission functions. The phase objects are considered as diffractive elements in a circular aperture of radius Ri with the transmittance function
t(x,y)=exp⁡(iπ·ΔΦ(x,y)),
(3)
where the phase function, ΔΦ(x,y), is determined with respect to Eqs. (1, 2) for different values of the parameters aM, bM, cM, aI, bI, and cI and the refractive indices n M and n I.

To calculate the Fresnel diffraction pattern, we have considered a monochromatic plane wave of wavelength λ, normally illuminating the structure of the samples with the transmittance function described by Eq. (3) situated in the input plane. The diffracted field in the output plane at the distance z from the sample, U( x0, y0,z), can be expressed using the 2D convolution operator [23

J. W. Goodman, Introduction to Fourier Optics (McGraw-Hill, 1968).

]:
U( x0, y0,z)=t(x,y,0)⊗h( x0, y0,x,y,z),
(4)
where
h( x0, y0,x,y,z)= exp⁡(ik r01) iλ r01
(5)
is the impulse response function, x0, y0 are the coordinates in the output plane, k=2π/λ, and r01 is the distance between any two points from the input plane to the output plane. This equation is a direct consequence of the Rayleigh–Sommerfeld theory. An extended version of the theoretical considerations can be found elsewhere [24

M. Mihailescu, “Natural quasi-periodic binary structure with focusing property in near-field diffraction pattern,” Opt. Express 18, 12526–12536 (2010). [CrossRef] [PubMed]

]. In Fresnel approximations, the impulse response function from Eq. (5) becomes
h(x,y,z)= exp⁡(ikz) iλzexp⁡{ ik 2z[ x2+ y2]}.
(6)

Consequently, the complex amplitude of the diffracted field is given by:
U( x 0, y 0,z)= exp⁡(ikz) iλzexp⁡[ ik 2z( x 0 2+ y 0 2)]⋅ ∫ −∞ +∞ ∫ −∞ +∞ U ′(x,y,0)⁢exp⁡[−i2π( x 0x+ y 0y)]dxdy,
(7)
where
U ′(x,y,0)=t(x,y,0)exp⁡[ ik 2z( x2+ y2)]
(8)
and U( x0, y0,z) is the 2D Fourier transform of the envelope U ′(x,y,0) at the spatial frequencies fx= x0/λz, fy= y0/λz up to multiplicative amplitude and phase factors.

For the numerical evaluation, these functions are sampled at the Nyquist frequency and Eqs. (7, 8) are implemented in discrete forms using 2D matrices and the fast Fourier transform routines in MATLAB. The values for equatorial radii were measured experimentally. The values for the refractive index depend on the hemoglobin concentration of the cell [13

M. Mir, H. Ding, Z. Wang, J. Reedy, K. Tangella, and G. Popescu, “Blood screening using diffraction phase cytometry,” J. Biomed. Opt. 15, 027016 (2010). [CrossRef] [PubMed]

]:
n= ns+β·MCHC,
(9)
where β is the refractive increment of hemoglobin ( 0.002 dl/g), MCHC is the concentration of dry protein expressed in the g/dl, and n s is the refractive index of the surrounding medium. For the refractive index of an MRBC, we consider a mean value of nM=1.418 [25

B. Rappaz, A. Barbul, Y. Emery, R. Korenstein, C. Depeursinge, P. J. Magistretti, and P. Marquet, “Comparative study of human erythrocytes by digital holographic microscopy, confocal microscopy, and impedance volume analyzer,” Cytometry Part A 73A, 895–903 (2008). [CrossRef]

], and the thickness values for all MRBCs are calculated from phase differences measured experimentally. For IRBCs, the refractive index is smaller [2

L. M. Lowenstein, “The Mammalian reticulocyte,” in International Review of Cytology , G. H. Bourne and J. F. Danielli, eds. (Academic, 1959), Vol. 8. [CrossRef] [PubMed]

, 26

W. Groner, N. Mohandas, and M. Bessis, “New optical technique for measuring erythrocyte deformability with the ektacytometer,” Clin. Chem. 26, 1435–1442 (1980). [PubMed]

, 27

F. M. Gaffney, “Experimental hemolytic anemia with particular reference to the corpuscular hemoglobin concentrations of the erythrocytes,” Br. J. Haematol. 3, 311–319 (1957). [CrossRef] [PubMed]

] and the hemoglobin concentration is about 20.5% lower than for MRBCs [28

G. d’Onofrio, R. Chirillo, G. Zini, G. Caenaro, M. Tommasi, and G. Micciulli, “Simultaneous measurement of reticulocyte and red blood cell indices in healthy subjects and patients with microcytic and macrocytic anemia,” Blood 85, 818–823 (1995). [PubMed]

]. Based on linear dependence, the IRBC mean refractive index is calculated at nI=1.396, assuming that the surrounding medium is homogeneous.

Simulations performed in MATLAB, in Fresnel approximation, confirm the diffraction pattern observed experimentally (Fig. 4). As in the experimental case, we also consider five distances, and the simulated diffraction patterns after propagation can be seen in Fig. 5 for aI= bI=4.5 μm, aM= bM=3.5 μm, cI=2 μm, and the height of the MRBC after concave cropping 3.3 μm. The intensity values and corresponding gray levels are arbitrary in all figures, normalized at the maximum value in each. For this reason, some backgrounds appear brighter than others.

To study the influence of cell parameters on their diffractive patterns, we scan different characteristic values for IRBCs and MRBCs at two given distances, zi and z i+1. Figure 6 illustrates the central intensity dependence on the equatorial and polar radii of the spheroid. The increment i, i+1 designates all distances where we found in the central intensity of the diffracted patterns maximum values for MRBC and minimum values for IRBC or vice versa.

4. Automated Procedure for Cell Separation, Identification, and Counting

In this procedure, the starting points are holograms of the transparent cells in plasma, recorded in hospital conditions in a DHM configuration. An algorithm for object image reconstruction and aberration compensation is used [29

F. Montfort, F. Charrière, T. Colomb, E. Cuche, P. Marquet, and C. Depeursinge, “Purely numerical compensation for microscope objective phase curvature in digital holographic microscopy: influence of digital phase masks position,” J. Opt. Soc. Am. A 23, 2944–2953 (2006). [CrossRef]

]. The sample images corresponding to the phase shift introduced by each cell (samples without a chromatographic agent) and their three-dimensional (3D) representations are given in Fig. 7. We apply a low-pass filter in the Fourier plane to remove small, unwanted details. In Fig. 7a, one can see that the entire sample plane is elevated in the back right corner. To adjust this, we use morphological structuring elements, which scan all pixels with their center and probe the neighboring area using a binary block. Its dimension can be scanned until the proper tilt is achieved [the background values are similar in all regions from Fig. 7c.] In Fig. 7b, the removed background is shown in a resized manner only for representation reasons (the indexing syntax is applied to view only one out of eight pixels in each direction; otherwise the surface plot would be too dense). In this sample, a white blood cell (of the largest radius), two IRBCs (the mountain shapes), and many MRBCS (crater shapes) are visible.

Some typical profiles are subtracted from these reconstructions (see Fig. 8). Their shapes prove that our choice to model the cells as oblate ellipsoids was appropriate. All these reconstructions are obtained from the hologram recorded in the in-focus plane, corrections at nanometric scale being achieved using the reconstruction algorithm (to find the minimum cell radius in amplitude reconstructions).

For automated cell–edge separation, identification, and volume calculation for IRBCs and MRBCs, we built a MATLAB code. It can read many images simultaneously and in the end exhibits the histograms, which include information extracted from these images for at least 1000 cells. To design this code, after background corrections presented in the previous paragraph, we use “the distance transform” which provides a metric or calculates the separation points in the image. The aim is to separate all cell images and detect their edges, centroids, radius, eccentricity, and volume. The resultant output images are shown in Fig. 9 (the input images are those of the reconstructed holograms). Notice that partial cell images are excluded (look at the edge for both samples). The distance transform individually identifies cells when they are far apart, but not when some cells are partially overlapped or adhere to each other (see the right region of Fig. 9).

To remove this shortcoming, we propose the following scheme based on several steps. Initially, we enhance the contrast using nonlinear functions applied to the image from Fig. 10a in threshold intervals. We apply an exponential decay for low-value pixels and a logarithmic function for large-value pixels (the commercial programs for contrast enhancement did not give good results). The threshold value is computed statistically for each image based on its pixel weights [Fig. 10b]. However, these functions enhance unwanted details as well, which are further removed through convolution with a Gaussian filter applied on 5×5 neighboring pixels [Fig. 10c].

Images obtained by applying the distance transform after contrast enhancement and Gaussian filter convolution are shown in Fig. 11. The improvement is less spectacular for spaced-out cells [see Fig. 11a], but highly visible for closely packed cells [compare the central right side of the Fig. 10b with the right side of Fig. 11b, where four cells are separately identified]. Comparing the border regions from Figs. 9b, 11b one can see that the incompletely pictured cells are excluded. As the whole image enlarges, the same cells are complete and visible enough to consider these cells in the counting process.

Because the contrast enhancement modifies the values of the pixels, to compute the morphological parameters, we consider the initial images only with Fourier filters for small details and tilt adjustment. Next, finding the properties of the individual cells in our image implies determining the coordinates of the centroids and the contained pixels based on the images with distance transform and edge detection using the “canny” operator from MATLAB. Based on this information and on the profile shapes, we can establish which cells are MRBCs and which are IRBCs by calculating a gradient from their centroids. For each cell, we add up the differences between the height of its centroid and the height of every other constituent pixel. These calculations are performed on the input images in the projected areas delimited after edge detection. For immature cells, which have profiles as plotted in Fig. 8b, the slope will be positive, summing up to a very large number. For mature cells [profiles as in Fig. 8a], the slope will be negative for the first points, starting from their centroid, then positive. Assuming that there is at least one immature cell in the image (our original images have 2048×2048 pixels and contain approximately 200 cells), the cell with the largest gradient number will certainly be an IRBC (or a white blood cell, but white blood cells will be excluded by the scheme, because they have greater radii). To identify the other immature cells, we chose, based on empirical results, those cells whose gradient numbers amount to at least one-third of the maximum. To go from phase shift to thickness values, we consider Eq. (1) and ns= n plasma=1.3515 [30

M. Mir, Z. Wang, K. Tangella, and G. Popescu, “Phase cytometry: blood on a CD-ROM,” Opt. Express 17, 2579–2585 (2009). [CrossRef] [PubMed]

].

In the same iterations, while calculating the gradient for each cell, we also determine its volume by adding up the heights of all the pixels inside the cell projection. Using this parameter, the IRBCs are well-differentiated by the MRBCs. At a 95% confidence limit, the percentage of reticulocytes is 1.65%±0.05% [see Fig. 12a] for a healthy person and 1.22%±0.05% for a suspected patient [see Fig. 12b]. These results were consistent with other measurement techniques: 1.72%±0.05% and 1.18±0.05% (impedance counter) and 1.70%±0.07% and 1.23%±0.07% (bright-field microscopy) performed on the first and second subsamples from each blood sample.

5. Discussion and Conclusions

Besides being used for automated counts, correct blood smear investigation is important in a wide range of diseases. Chromatographic agents are used to visualize and identify different cell types, but some procedures are time-consuming. This is also the case for the classical method for finding the percentage of IRBCs in the total number of red blood cells. For this reason, methods to simultaneously image, and identify/count are sought by scientists.

We experimentally revealed that at different distance intervals between the sample and the objective lens, IRBCs and MRBCs exhibit complementary central maxima and minima in their diffractive patterns. To simulate the propagation in Fresnel approximation, we model IRBCs and MRBCs as oblate spheroids with experimentally established dimensions. We use a mean value for the refractive index, which takes into account hemoglobin concentration. One can observe the succession of the central maxima and minima, which is similar to that seen in the experimental diffraction pattern. The dependencies of the central intensity with regard to the spheroid parameters reveal that the transversal dimensions have a more dominant influence on the diffractive pattern than the longitudinal ones. Between the mean radius of MRBCs and the mean radius of IRBCs there is only one central intensity maximum in their respective diffractive patterns. The equatorial and polar radii do influence the central intensity in the same manner: for both we have a maximum at MRBC parameters and a minimum for IRBCs at a given zi and vice versa at z i+1. We find few values for the distance zi within the studied interval. Between these values are transition zones, but we never find maximum central intensity for both IRBC and MRBC radii. These preliminary data will be used as preprocessing observations in our microfluidic applications, where every cell has its own distance to the focal plane.

Next, we developed a MATLAB code for combining the advantage of a less time-consuming cell identification-count procedure with the acquisition of other morphological parameters, for discriminating similar cells. The input images are obtained in the re construction step of the DHM technique. Since the distance transform did not separate some adjacent cells, we proposed a scheme based on the nonlinear functions applied at threshold intervals for contrast enhancement and Gaussian filter convolution. As a result, our MATLAB code recognizes each centroid, identifies cell type, counts IRBCs and MRBCs from at least six images (over 1000 cells), and builds the corresponding histograms automatically. White blood cells are separated for their greater radii. Our measurements are in good agreement with those obtained in bright-field microscopy or using the impedance counter.

The developed code can also exhibit histograms after radius, projection area, projection eccentricity, or other morphological parameters, using the values given by the “regionprops” MATLAB operator. The eccentricity parameter is important because a high value indicates that the cells have abnormal shapes— a known signal for many diseases. The eccentricity is greater than 0.7 for oval RBCs and even greater than 1 for pear-shaped, comma-shaped or sickle cells. For the latter types, careful observations are necessary.

Our results prove that DHM combined with image processing is a suitable instrument for simultaneously performing noncontact, label-free individual RBC morphological investigation and real time imaging, identification, and counting. Important in our studies is the fact that the measurements were performed on different slides prepared in hospital conditions with many tangential cells and the procedure separates IRBCs from MRBCs.

Acknowledgments

The research presented in this paper is supported by the Sectorial Operational Program for Human Resource Development financed by the European Social Fund and by the Romanian Government under contract POSDRU/89/1.5/S/63700. The equipment used in these experiments were acquired using funds from the Romanian contract 4/CP/I/2007-2009 Capacities.

References and links

1.

E. B. Krumbhaar, “Reticulosis-increased percentage in the peripheral blood,” in Hematology: Landmark Papers of the Twentieth Century , M. A. Lichtman, L. A. Boxer, E. Henderson, and J. L. Spivak, eds. (Academic, 2000), p. 805.

2.

L. M. Lowenstein, “The Mammalian reticulocyte,” in International Review of Cytology , G. H. Bourne and J. F. Danielli, eds. (Academic, 1959), Vol. 8. [CrossRef] [PubMed]

3.

L. Blanc, A. De Gassart, C. Géminard, P. Bette-Bobillo, and M. Vidal, “Exosome release by reticulocytes—an integral part of the red blood cell differentiation system,” Blood Cells Mol. Dis. 35, 21–26 (2005). [CrossRef] [PubMed]

4.

J. G. Stephens, “Surface and fragility differences between mature and immature red cells,” J. Physiol. 99, 30–48 (1940). [PubMed]

5.

A. K. Percy, E. Schmell, B. J. Earles, and W. J. Lennarz, “Phospholipid biosynthesis in the membranes of immature and mature red blood cells,” Biochemistry 12, 2456–2461 (1973). [CrossRef] [PubMed]

6.

T. Lecklin, A. Tuominen, and M. Nikinmaa, “The adrenergic volume changes of immature and mature rainbow trout (Oncorhynchus mykiss) erythrocytes,” J. Exp. Biol. 203, 3025–3031 (2000). [PubMed]

7.

Y. M. Serebrennikova, J. Patel, W. K. Milhous, and L. H. Garcıa-Rubio, “Quantitative analysis of morphological alterations in Plasmodium falciparum infected red blood cells through theoretical interpretation of spectral measurements,” J. Therm. Biol. 265, 493–500 (2010). [CrossRef]

8.

J. C. Thompson and B. W. Manktelow, “Pathogenesis and red blood cell destruction in hemoglobinemic leptospirosis,” J. Comp. Pathol. 96, 529–540 (1986). [CrossRef] [PubMed]

9.

S. Christel and C. Little, “Morphological changes during heating of erythrocytes from stored human blood,” J. Therm. Biol. 9, 221–228 (1984). [CrossRef]

10.

I. A. Kabat, W. Leyko, B. Kwiatkowski, and I. Zakrzewska, “Osmotic properties and morphological changes of submicroscopic surface structures of mammalian red blood cells subjected to UV-irradiation in vitro ,” Zentralbl. Bakteriol. Orig. B 159, 88–94 (1974). [PubMed]

11.

A. J. McGoron, C. H. Joiner, M. B. Palascak, W. J. Claussen, and R. S. Franco, “Dehydration of mature and immature sickle red blood cells during fast oxygenation/deoxygenation cycles: role of KCl cotransport and extracellular calcium,” Blood 95, 2164–2168 (2000). [PubMed]

12.

N. Mochandas, M. R. Clark, M. S. Jacobs, and S. B. Shohet, “Analysis of factors regulating erythrocyte deformability,” J. Clin. Invest. 66, 563–573 (1980). [CrossRef]

13.

M. Mir, H. Ding, Z. Wang, J. Reedy, K. Tangella, and G. Popescu, “Blood screening using diffraction phase cytometry,” J. Biomed. Opt. 15, 027016 (2010). [CrossRef] [PubMed]

14.

S. Rancourt-Grenier, M.-T. Wei, J.-J. Bai, A. Chiou, P. P. Bareil, P.-L. Duval, and Y. Sheng, “Dynamic deformation of red blood cell in dual-trap optical tweezers,” Opt. Express 18, 10462–10472 (2010). [CrossRef] [PubMed]

15.

C. Fang-Yen, S. Oh, Y. Park, W. Choi, S. Song, H. S. Seung, R. R. Dasari, and M. S. Feld, “Imaging voltage-dependent cell motions with heterodyne Mach–Zehnder phase microscopy,” Opt. Lett. 32, 1572–1574 (2007). [CrossRef] [PubMed]

16.

B.-W. Yang and Z. Li, “Measuring microinteractions between coagulating red blood cells using optical tweezers,” Biomed. Opt. Express 1, 1217–1224 (2010). [CrossRef]

17.

Y. K. Park, C. A. Best, K. Badizadegan, R. R. Dasari, M. S. Feld, T. Kuriabova, M. L. Henle, A. J. Levine, and G. Popescu, “Measurement of red blood cell mechanics during morphological changes,” Proc. Natl. Acad. Sci. USA 107, 6731–6736 (2010). [CrossRef] [PubMed]

18.

Y.-S. Choi and S. Lee, “Three-dimensional volumetric measurement of red blood cell motion using digital holographic microscopy,” Appl. Opt. 48, 2983–2990 (2009). [CrossRef] [PubMed]

19.

B. Rappaz, A. Barbul, A. Hoffmann, D. Boss, R. Korenstein, C. Depeursinge, P. J. Magistretti, and P. Marquet, “Spatial analysis of erythrocyte membrane fluctuations by digital holographic microscopy,” Blood Cells Mol. Dis. 42, 228–232 (2009). [CrossRef] [PubMed]

20.

M. Paturzo, F. Merola, S. Grilli, S. De Nicola, A. Finizio, and P. Ferraro, “Super-resolution in digital holography by a two-dimensional dynamic phase grating,” Opt. Express 16, 17107–17118 (2008). [CrossRef] [PubMed]

21.

V. Mico, Z. Zalevsky, C. Ferreira, and J. García, “Super-resolution digital holographic microscopy for three-dimensional samples,” Opt. Express 16, 19260–19270 (2008). [CrossRef]

22.

G. J. Streekstra1, J. G. G. Dobbe, and A. G. Hoekstra, “Quantification of the fraction poorly deformable red blood cells using ektacytometry,” Opt. Express 18, 14173–14182 (2010). [CrossRef]

23.

J. W. Goodman, Introduction to Fourier Optics (McGraw-Hill, 1968).

24.

M. Mihailescu, “Natural quasi-periodic binary structure with focusing property in near-field diffraction pattern,” Opt. Express 18, 12526–12536 (2010). [CrossRef] [PubMed]

25.

B. Rappaz, A. Barbul, Y. Emery, R. Korenstein, C. Depeursinge, P. J. Magistretti, and P. Marquet, “Comparative study of human erythrocytes by digital holographic microscopy, confocal microscopy, and impedance volume analyzer,” Cytometry Part A 73A, 895–903 (2008). [CrossRef]

26.

W. Groner, N. Mohandas, and M. Bessis, “New optical technique for measuring erythrocyte deformability with the ektacytometer,” Clin. Chem. 26, 1435–1442 (1980). [PubMed]

27.

F. M. Gaffney, “Experimental hemolytic anemia with particular reference to the corpuscular hemoglobin concentrations of the erythrocytes,” Br. J. Haematol. 3, 311–319 (1957). [CrossRef] [PubMed]

28.

G. d’Onofrio, R. Chirillo, G. Zini, G. Caenaro, M. Tommasi, and G. Micciulli, “Simultaneous measurement of reticulocyte and red blood cell indices in healthy subjects and patients with microcytic and macrocytic anemia,” Blood 85, 818–823 (1995). [PubMed]

29.

F. Montfort, F. Charrière, T. Colomb, E. Cuche, P. Marquet, and C. Depeursinge, “Purely numerical compensation for microscope objective phase curvature in digital holographic microscopy: influence of digital phase masks position,” J. Opt. Soc. Am. A 23, 2944–2953 (2006). [CrossRef]

30.

M. Mir, Z. Wang, K. Tangella, and G. Popescu, “Phase cytometry: blood on a CD-ROM,” Opt. Express 17, 2579–2585 (2009). [CrossRef] [PubMed]

Fig. 1 (a) Classical bright-field microscopy on a slide using chromatographic substances to differentiate the cells; (b) hologram of the same slide in the region of interest recorded using off-axis DHM setup; (c) 3D image in false colors of the reconstructed cells from the previous hologram. The IRBCs are highlighted in all images.
Fig. 2 IRBCs and MRBCs observed in classical DIC microscopy without chromatographic substances. (a) Sample 3, (c) Sample 13; and in DHM (b) Sample 3, (d) Sample 13. The IRBCs are highlighted in all images.
Fig. 3 Holograms from regions with at least one IRBC recorded at different distances between the microscope objective and sample: (a) Sample 27 at z1, (b) Sample 27 at z2, (c) Sample 19 at z3, (d) Sample 19 at z4. The IRBCs are highlighted in all images.
Fig. 4 Diffraction pattern (with reference beam obstructed in DHM setup) from an MRBC (top) and an IRBC (bottom) recorded at different distances between the microscope objective and Sample 27 in the same region as in Fig. 3 at (a)  z=9 μm, (b)  z=6 μm, (c)  z=3 μm, (d)  z=0 μm, (e)  z=−6 μm.
Fig. 5 Simulated diffraction pattern from MRBC (top) and IRBC (bottom) at the same distances between sample and objective as in Fig. 4.
Fig. 6 Central intensity dependence on (a) equatorial radius at zi, (b) equatorial radius at z i+1, (c) polar radius at zi, (d) polar radius at z i+1.
Fig. 7 (a) Initial 3D phase reconstruction, (b) background correction, (c) 3D phase reconstruction with tilt adjustment.
Fig. 8 Phase shift profiles for (a) MRBC, (b) IRBC.
Fig. 9 Distance transform applied on the images obtained from the hologram reconstruction of (a) Sample 3, (b) Sample 5.
Fig. 10 (a) Initial image. (b) Same image with contrast enhancement using nonlinear functions. (c) Convolution of the previous image with a Gaussian filter.
Fig. 11 Distance transform applied on the image for (a) Sample 3 and (b) Sample 5.
Fig. 12 Histograms built after volume for (a) healthy person and (b) suspected patient.

OCIS Codes
(090.0090) Holography : Holography
(090.1970) Holography : Diffractive optics
(170.0180) Medical optics and biotechnology : Microscopy
(090.1995) Holography : Digital holography

ToC Category:
Holography

History
Original Manuscript: February 28, 2011
Manuscript Accepted: May 4, 2011
Published: July 7, 2011

Virtual Issues
Vol. 6, Iss. 8 Virtual Journal for Biomedical Optics
July 6, 2011 Spotlight on Optics

Citation
Mona Mihailescu, Mihaela Scarlat, Alexandru Gheorghiu, Julia Costescu, Mihai Kusko, Irina Alexandra Paun, and Eugen Scarlat, "Automated imaging, identification, and counting of similar cells from digital hologram reconstructions," Appl. Opt. 50, 3589-3597 (2011)
http://www.opticsinfobase.org/vjbo/abstract.cfm?URI=ao-50-20-3589


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. E. B. Krumbhaar, “Reticulosis-increased percentage in the peripheral blood,” in Hematology: Landmark Papers of the Twentieth Century, M.A.Lichtman, L.A.Boxer, E.Henderson, and J.L.Spivak, eds. (Academic, 2000), p. 805.
  2. L. M. Lowenstein, “The Mammalian reticulocyte,” in International Review of Cytology, G.H.Bourne and J.F.Danielli, eds. (Academic, 1959), Vol.  8. [CrossRef] [PubMed]
  3. L. Blanc, A. De Gassart, C. Géminard, P. Bette-Bobillo, and M. Vidal, “Exosome release by reticulocytes—an integral part of the red blood cell differentiation system,” Blood Cells Mol. Dis. 35, 21–26 (2005). [CrossRef] [PubMed]
  4. J. G. Stephens, “Surface and fragility differences between mature and immature red cells,” J. Physiol. 99, 30–48 (1940). [PubMed]
  5. A. K. Percy, E. Schmell, B. J. Earles, and W. J. Lennarz, “Phospholipid biosynthesis in the membranes of immature and mature red blood cells,” Biochemistry 12, 2456–2461 (1973). [CrossRef] [PubMed]
  6. T. Lecklin, A. Tuominen, and M. Nikinmaa, “The adrenergic volume changes of immature and mature rainbow trout (Oncorhynchus mykiss) erythrocytes,” J. Exp. Biol. 203, 3025–3031(2000). [PubMed]
  7. Y. M. Serebrennikova, J. Patel, W. K. Milhous, and L. H. Garc?a-Rubio, “Quantitative analysis of morphological alterations in Plasmodium falciparum infected red blood cells through theoretical interpretation of spectral measurements,” J. Therm. Biol. 265, 493–500 (2010). [CrossRef]
  8. J. C. Thompson and B. W. Manktelow, “Pathogenesis and red blood cell destruction in hemoglobinemic leptospirosis,” J. Comp. Pathol. 96, 529–540 (1986). [CrossRef] [PubMed]
  9. S. Christel and C. Little, “Morphological changes during heating of erythrocytes from stored human blood,” J. Therm. Biol. 9, 221–228 (1984). [CrossRef]
  10. I. A. Kabat, W. Leyko, B. Kwiatkowski, and I. Zakrzewska, “Osmotic properties and morphological changes of submicroscopic surface structures of mammalian red blood cells subjected to UV-irradiation in vitro,” Zentralbl. Bakteriol. Orig. B 159, 88–94 (1974). [PubMed]
  11. A. J. McGoron, C. H. Joiner, M. B. Palascak, W. J. Claussen, and R. S. Franco, “Dehydration of mature and immature sickle red blood cells during fast oxygenation/deoxygenation cycles: role of KCl cotransport and extracellular calcium,” Blood 95, 2164–2168 (2000). [PubMed]
  12. N. Mochandas, M. R. Clark, M. S. Jacobs, and S. B. Shohet, “Analysis of factors regulating erythrocyte deformability,” J. Clin. Invest. 66, 563–573 (1980). [CrossRef]
  13. M. Mir, H. Ding, Z. Wang, J. Reedy, K. Tangella, and G. Popescu, “Blood screening using diffraction phase cytometry,” J. Biomed. Opt. 15, 027016 (2010). [CrossRef] [PubMed]
  14. S. Rancourt-Grenier, M.-T. Wei, J.-J. Bai, A. Chiou, P. P. Bareil, P.-L. Duval, and Y. Sheng, “Dynamic deformation of red blood cell in dual-trap optical tweezers,” Opt. Express 18, 10462–10472 (2010). [CrossRef] [PubMed]
  15. C. Fang-Yen, S. Oh, Y. Park, W. Choi, S. Song, H. S. Seung, R. R. Dasari, and M. S. Feld, “Imaging voltage-dependent cell motions with heterodyne Mach–Zehnder phase microscopy,” Opt. Lett. 32, 1572–1574 (2007). [CrossRef] [PubMed]
  16. B.-W. Yang and Z. Li, “Measuring microinteractions between coagulating red blood cells using optical tweezers,” Biomed. Opt. Express 1, 1217–1224 (2010). [CrossRef]
  17. Y. K. Park, C. A. Best, K. Badizadegan, R. R. Dasari, M. S. Feld, T. Kuriabova, M. L. Henle, A. J. Levine, and G. Popescu, “Measurement of red blood cell mechanics during morphological changes,” Proc. Natl. Acad. Sci. USA 107, 6731–6736 (2010). [CrossRef] [PubMed]
  18. Y.-S. Choi and S. Lee, “Three-dimensional volumetric measurement of red blood cell motion using digital holographic microscopy,” Appl. Opt. 48, 2983–2990 (2009). [CrossRef] [PubMed]
  19. B. Rappaz, A. Barbul, A. Hoffmann, D. Boss, R. Korenstein, C. Depeursinge, P. J. Magistretti, and P. Marquet, “Spatial analysis of erythrocyte membrane fluctuations by digital holographic microscopy,” Blood Cells Mol. Dis. 42, 228–232 (2009). [CrossRef] [PubMed]
  20. M. Paturzo, F. Merola, S. Grilli, S. De Nicola, A. Finizio, and P. Ferraro, “Super-resolution in digital holography by a two-dimensional dynamic phase grating,” Opt. Express 16, 17107–17118 (2008). [CrossRef] [PubMed]
  21. V. Mico, Z. Zalevsky, C. Ferreira, and J. García, “Super-resolution digital holographic microscopy for three-dimensional samples,” Opt. Express 16, 19260–19270 (2008). [CrossRef]
  22. G. J. Streekstra1, J. G. G. Dobbe, and A. G. Hoekstra, “Quantification of the fraction poorly deformable red blood cells using ektacytometry,” Opt. Express 18, 14173–14182 (2010). [CrossRef]
  23. J. W. Goodman, Introduction to Fourier Optics (McGraw-Hill, 1968).
  24. M. Mihailescu, “Natural quasi-periodic binary structure with focusing property in near-field diffraction pattern,” Opt. Express 18, 12526–12536 (2010). [CrossRef] [PubMed]
  25. B. Rappaz, A. Barbul, Y. Emery, R. Korenstein, C. Depeursinge, P. J. Magistretti, and P. Marquet, “Comparative study of human erythrocytes by digital holographic microscopy, confocal microscopy, and impedance volume analyzer,” Cytometry Part A 73A, 895–903 (2008). [CrossRef]
  26. W. Groner, N. Mohandas, and M. Bessis, “New optical technique for measuring erythrocyte deformability with the ektacytometer,” Clin. Chem. 26, 1435–1442 (1980). [PubMed]
  27. F. M. Gaffney, “Experimental hemolytic anemia with particular reference to the corpuscular hemoglobin concentrations of the erythrocytes,” Br. J. Haematol. 3, 311–319 (1957). [CrossRef] [PubMed]
  28. G. d’Onofrio, R. Chirillo, G. Zini, G. Caenaro, M. Tommasi, and G. Micciulli, “Simultaneous measurement of reticulocyte and red blood cell indices in healthy subjects and patients with microcytic and macrocytic anemia,” Blood 85, 818–823 (1995). [PubMed]
  29. F. Montfort, F. Charrière, T. Colomb, E. Cuche, P. Marquet, and C. Depeursinge, “Purely numerical compensation for microscope objective phase curvature in digital holographic microscopy: influence of digital phase masks position,” J. Opt. Soc. Am. A 23, 2944–2953 (2006). [CrossRef]
  30. M. Mir, Z. Wang, K. Tangella, and G. Popescu, “Phase cytometry: blood on a CD-ROM,” Opt. Express 17, 2579–2585 (2009). [CrossRef] [PubMed]

Cited By

 Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited