• Photonics Research
  • Vol. 11, Issue 12, 2242 (2023)
Liming Yang1、2、†, Ruihai Wang2、†, Qianhao Zhao2, Pengming Song2, Shaowei Jiang2, Tianbo Wang2, Xiaopeng Shao1、3、*, Chengfei Guo2, Rishikesh Pandey2, and Guoan Zheng2、4、*
Author Affiliations
  • 1School of Optoelectronic Engineering, Xidian University, Xi’an 710071, China
  • 2Department of Biomedical Engineering, University of Connecticut, Storrs, Connecticut 06269, USA
  • 3e-mail: xpshao@xidian.edu.cn
  • 4e-mail: guoan.zheng@uconn.edu
  • show less
    DOI: 10.1364/PRJ.504378 Cite this Article Set citation alerts
    Liming Yang, Ruihai Wang, Qianhao Zhao, Pengming Song, Shaowei Jiang, Tianbo Wang, Xiaopeng Shao, Chengfei Guo, Rishikesh Pandey, Guoan Zheng. Lensless polarimetric coded ptychography for high-resolution, high-throughput gigapixel birefringence imaging on a chip[J]. Photonics Research, 2023, 11(12): 2242 Copy Citation Text show less

    Abstract

    Polarimetric imaging provides valuable insights into the polarization state of light interacting with a sample. It can infer crucial birefringence properties of specimens without using labels, thereby facilitating the diagnosis of diseases such as cancer and osteoarthritis. In this study, we present a novel polarimetric coded ptychography (pol-CP) approach that enables high-resolution, high-throughput gigapixel birefringence imaging on a chip. Our platform deviates from traditional lens-based systems by employing an integrated polarimetric coded sensor for lensless coherent diffraction imaging. Utilizing Jones calculus, we quantitatively determine the birefringence retardance and orientation information of biospecimens from the recovered images. Our portable pol-CP prototype can resolve the 435 nm linewidth on the resolution target, and the imaging field of view for a single acquisition is limited only by the detector size of 41 mm2. The prototype allows for the acquisition of gigapixel birefringence images with a 180 mm2 field of view in 3.5 min, a performance that rivals high-end whole slide scanner but at a small fraction of the cost. To demonstrate its biomedical applications, we perform high-throughput imaging of malaria-infected blood smears, locating parasites using birefringence contrast. We also generate birefringence maps of label-free thyroid smears to identify thyroid follicles. Notably, the recovered birefringence maps emphasize the same regions as autofluorescence images, underscoring the potential for rapid on-site evaluation of label-free biopsies. Our approach provides a turnkey and portable solution for lensless polarimetric analysis on a chip, with promising applications in disease diagnosis, crystal screening, and label-free chemical imaging, particularly in resource-constrained environments.

    1. INTRODUCTION

    Polarimetric imaging has attracted considerable interest in the field of biomedicine, owing to its ability to unveil structural and chemical information that often remains undetected by conventional imaging techniques [1,2]. A notable benefit of polarimetric imaging is its application in label-free biospecimen characterization. By analyzing the polarization properties of light interacting with biospecimens, microscopic-level information about the composition and structure can be obtained without requiring labeling agents. For example, collagen fibers in connective tissues can alter the polarization state of light passing through them. Measuring this alteration provides information about collagen fiber orientation and density, which are associated with various diseases such as cancer and fibrosis [3,4]. Another advantage of polarimetric imaging is image contrast enhancement. By selectively filtering out light waves with certain polarization states, it is possible to enhance image contrast and improve the visibility of certain structures or components based on the birefringence maps. For example, birefringence maps have been used to locate and identify malaria hemozoin crystals generated from the hemoglobin metabolism process of malaria parasites [5]. Similarly, birefringence maps have shown clinical value in screening and imaging meiotic spindles [6], monosodium urate [7], squamous cell carcinoma [8,9], articular cartilage [10], cerebral amyloid [11], normal versus cancer cells [12], among others. These diverse applications demonstrate the potential of polarimetric imaging in enhancing diagnostic accuracy and facilitating the early detection of various medical conditions.

    Conventional lens-based polarization light microscopy utilizes a polarizer-analyzer pair to measure the changes in the polarization state of light induced by the sample. Despite widespread use in different application settings, there are several limitations of these lens-based systems. One significant drawback is the trade-off between the spatial resolution and the imaging field of view. One can obtain high resolution with a small field of view or low resolution with a large field of view, but not both [13,14]. To obtain a high-resolution whole slide image of the sample, integration of the polarization optics with a whole slide imaging system is necessary for data acquisition. However, this integration remains elusive as whole slide scanners are often dedicated and costly tools [15,16]. Additionally, high numerical aperture (NA) objective lenses imply a small depth of field, typically on the micrometer scale, which presents a challenge when tracking the axial topography variations of biospecimen. If the sample is not positioned within the objective lens’s small depth of field, the image quality will degrade, necessitating rescanning and causing workflow delays [15]. Further, the specialized and nonportable nature of whole slide scanner creates obstacles that prevent its widespread use in different application settings, particularly in resource-constrained environments.

    To address these limitations, researchers have developed alternative methods for polarimetric imaging. One notable advancement is in the domain of conventional real-space ptychography, where researchers have integrated a vectorial imaging framework to enable the characterization of anisotropic media [1719]. This development enables mapping of the four parameters of the Jones matrix and can be used to quantify optical material properties. Another example is the lens-based Fourier ptychographic microscopy (FPM) [20] implemented with a pixelated polarization camera [21]. In this case, the achieved spatial resolution is no longer constrained by the NA of the objective lens but determined by the largest incident angle. Consequently, a half-pitch resolution of 0.55 μm has been demonstrated over a field of view of 3.78 mm2 in a polarization-sensitive FPM system, corresponding to a moderate space bandwidth product of 12.5 megapixels [21]. Similarly, vectorial FPM has been demonstrated by placing a generator polarizer at the illumination path and an analyzer at the detection path [22]. With four generator-analyzer configurations, the object’s Jones matrix can be recovered from the low-resolution raw measurements.

    While these previous demonstrations provide new insights for polarimetric imaging, the achieved imaging throughput is still orders of magnitude lower than that of a whole slide scanner [15,23]. Moreover, real-space ptychography and FPM rely on coded illumination, and their imaging models are contingent upon how the incident beam enters the sample [14,2426]. As a result, thick specimens frequently pose a challenge to these modalities, necessitating an accurate model of the light-sample interaction process [2729]. This same challenge also applies to other lensless ptychographic implementations that utilize speckle or multiangle illumination [3035]. The quality of the recovered images often falls short when compared with that of conventional bright-field polarization microscopy.

    Another alternative to conventional lens-based polarization light microscopy is to perform polarimetric imaging based on in-line holographic acquisitions [3640]. For example, multiheight holographic measurements have been implemented with a polarizer–analyzer pair for lensless polarimetric imaging, demonstrating potentials for biomedical applications [36]. However, the traditional multiheight phase-retrieval technique may not provide sufficient diversity measurements for recovering complex objects with intensity and phase properties [41]. Specifically, recovering biospecimens containing 2π phase wraps, such as bacterial colonies or cytology smears, is a challenging task for the in-line holographic implementation [41,42]. This difficulty arises because the object’s phase information cannot be effectively converted into intensity variations for detection, i.e., the phase transfer function approaches zero at certain frequency bands [24,43,44]. Moreover, existing holographic imaging implementations still face constraints in spatial resolution, field of view, and imaging throughput when compared with robotic microscope platforms.

    In this work, we present a novel polarimetric coded ptychography (pol-CP) technique that enables high-resolution, high-throughput gigapixel birefringence imaging on a chip. Instead of using conventional lens-based polarization optics, our platform adopts the concept of coded ptychography for high-resolution, high-throughput coherent diffraction imaging with a disorder-engineered surface [42,4450]. In our implementation, we illuminate the specimen with circularly polarized laser light and attach an ultrathin polarizer film to a blood-coated sensor for image acquisition. The integrated coded sensor is then mounted on a compact, customized stage, enabling control over both rotational and translational motions. By analyzing four different polarization states with Jones calculus, we retrieve the quantitative birefringence retardance and orientation information of the specimen. Our cost-effective and handheld pol-CP prototype can resolve the 435 nm linewidth on the resolution target, with the field of view of a single acquisition limited only by the size of the sensor, i.e., approximately 41 mm2. Utilizing this prototype platform, we demonstrate that gigapixel birefringence images with a 180 mm2 field of view can be acquired in approximately 3.5 min, a performance that rivals high-end whole slide scanner but at a small fraction of the cost. In contrast with conventional lens-based polarization microscopy, our platform enables image refocusing after data acquisition, effectively addressing the autofocusing issue present in traditional whole slide scanners. Furthermore, the reported technique models only the complex wavefront exiting the sample, not its entry, rendering sample thickness inconsequential in the image formation process. Unlike conventional real-space ptychography, FPM, and other common ptychographic implementations, which necessitate thin sample assumptions, our platform eliminates this constraint. Given these advantages, we anticipate that this approach will find widespread applications in disease diagnosis, crystal screening, and label-free chemical imaging, particularly in resource-constrained environments.

    2. MATERIALS AND METHODS

    A. Polarimetric Coded Ptychography (pol-CP) on a Chip

    Figures 1(a) and 1(b) illustrate the schematic of our pol-CP prototype setup. This prototype consists of three main components (from top to bottom): the illumination light source; the integrated polarimetric coded sensor; and the motion system for programmable translation and rotation control. As shown in Fig. 1(a), a fiber-coupled, 5 mW 405 nm laser diode serves as our light source. Positioned upstream of the sample, a left-handed circular polarizer converts the linearly polarized light from the laser into circularly polarized light for sample illumination. The use of laser diode in pol-CP offers superior spatial and temporal coherence compared with conventional LED-based lensless setups. It also ensures adequate illumination intensity for a short exposure time of 1  ms for image acquisition. This allows the coded sensor to operate in continuous motion without occurrence of motion blur issues.

    Schematic diagram and operation of the lensless polarimetric coded ptychography (pol-CP) platform. (a) A circular polarizer transforms linearly polarized light from the laser diode into circularly polarized light, which then interacts with the specimen. The exit wave is recorded by the integrated polarimetric coded image sensor. (b) Illustration of changes in the polarization state within the pol-CP system. (c) A 2 μL sample of goat blood is smeared onto the sensor’s coverglass and fixed with ethyl alcohol, forming a thin and dense blood-cell layer that acts as a high-performance scattering lens with a theoretically unlimited field of view. (d) An ultrathin polarizing film is adhered to the blood-cell layer, forming an integrated polarimetric coded image sensor. (e) The pol-CP prototype features a polarimetric coded sensor mounted on programmable stages, allowing control over both rotational and translational motions.

    Figure 1.Schematic diagram and operation of the lensless polarimetric coded ptychography (pol-CP) platform. (a) A circular polarizer transforms linearly polarized light from the laser diode into circularly polarized light, which then interacts with the specimen. The exit wave is recorded by the integrated polarimetric coded image sensor. (b) Illustration of changes in the polarization state within the pol-CP system. (c) A 2 μL sample of goat blood is smeared onto the sensor’s coverglass and fixed with ethyl alcohol, forming a thin and dense blood-cell layer that acts as a high-performance scattering lens with a theoretically unlimited field of view. (d) An ultrathin polarizing film is adhered to the blood-cell layer, forming an integrated polarimetric coded image sensor. (e) The pol-CP prototype features a polarimetric coded sensor mounted on programmable stages, allowing control over both rotational and translational motions.

    The imaging system of pol-CP features an integrated polarimetric coded image sensor. In our implementation, we smeared a drop of goat blood on the sensor’s coverglass followed by cell fixation with ethyl alcohol [Fig. 1(c)]. We select goat blood because of its 2–3 μm red blood cell size, which is the smallest among all animals. The smearing process naturally creates a thin, uniform, and dense monolayer of blood cells on the sensor’s coverglass for light wave modulation. This blood-cell layer serves as the effective ptychographic probe beam within the pol-CP setup [24]. Alternative solutions include disorder-engineered surface, microbeads, and metasurface [42,44,49,50]. We note that the fabrication process of this goat-blood layer can be completed in 5  min with minimal cost; notably, it does not require the use of complex lithographic tools typically required for the construction of metasurfaces or disorder-engineered surfaces [42,44]. In our pol-CP prototype, we used a drop of polydimethylsiloxane (PDMS) liquid to attach an 80  μm polarizing film to the coded sensor, as shown in Fig. 1(d). The use of PDMS ensures a secure attachment between the polarizing film and the blood-cell layer, thus forming an integrated polarimetric coded sensor capable of high-resolution ptychographic imaging. As shown in Fig. 1(b), when the specimen exhibits optical anisotropy, the transmitted light becomes elliptically polarized and is subsequently collected by the integrated coded sensor.

    The motion system of our prototype device is shown in the bottom panel of Fig. 1(a). We used stepper motors to build a compact motorized stage capable of managing translational and rotational motion controls. Figure 1(f) shows the entire pol-CP prototype platform. The operation of the entire platform can be viewed in Visualization 1.

    B. Imaging Model and Reconstruction

    Figure 2 shows the data acquisition and reconstruction process of pol-CP. In the data acquisition process, the integrated polarimetric coded sensor is rotated to four specific angles (α=0°, 45°, 90°, 135°). At each angle, a set of lensless diffraction data is obtained by translating the coded sensor to different lateral positions (xi,yi). The four sets of measurements, labelled as I1α,  ,  INα in Fig. 2(a), are subsequently utilized to recover four high-resolution images corresponding to different polarization states, as displayed in Fig. 2(b). From these recovered intensity images, the optic axis orientation and retardance maps can be derived using Jones calculus, as demonstrated in Figs. 2(c)–2(e).

    Image acquisition and reconstruction process of pol-CP. (a) The captured four sets of raw images with the integrated coded sensor at four different orientations. (b) The recovered four images with the 0° orientation as the reference. (c), (d) The optic axis orientation and retardance maps are derived using Jones calculus. (e) A pseudo-colored birefringence map by encoding the optic axis orientation as color hue and the retardance information as intensity.

    Figure 2.Image acquisition and reconstruction process of pol-CP. (a) The captured four sets of raw images with the integrated coded sensor at four different orientations. (b) The recovered four images with the 0° orientation as the reference. (c), (d) The optic axis orientation and retardance maps are derived using Jones calculus. (e) A pseudo-colored birefringence map by encoding the optic axis orientation as color hue and the retardance information as intensity.

    The forward imaging model of pol-CP is discussed as follows. Initially, the input left-handed circularly polarized light interacts with the sample, producing a wave field that travels a distance d1 between the sample and the integrated coded sensor. This wave field is then modulated by the linear polarizing film affixed to the integrated coded sensor. The imaging process above can be written as Eoutα(x,y)=Jpolα·PSFfree(d1)*(Jsample(x,y)·Ein),where · represents matrix multiplication, and * represents convolution operation applied element-wise to each component of the matrix. In Eq. (1), Ein=t/2[1,i]T, representing the Jones vector of the input circularly polarized light with intensity t. The term PSFfree(d1) represents the convolution kernel for free space propagation over a distance d1. Within this equation, Jsample(x,y) and Jpolα correspond to the Jones matrices of the birefringent sample and the linear polarizer orientated at angle α, respectively. Specifically, the Jones matrix of the sample Jsample(x,y) can be expressed as Jsample(x,y)=[cos(δ2)+isin(δ2)cos(2θ)isin(δ2)sin(2θ)isin(δ2)sin(2θ)cos(δ2)isin(δ2)cos(2θ)].Here, the symbol denotes element-wise multiplication. For brevity, the spatial coordinates (x,y) are omitted in the retardance magnitude δ(x,y) and the slow axis orientation θ(x,y) in Eq. (2). On a related note, the Jones matrix of the linear polarizer, represented as Jpolα, is given by Jpolα=[cos2(α)cos(α)sin(α)cos(α)sin(α)sin2(α)].

    In this equation, α signifies the orientation angle of the linear polarizing film on the coded sensor. By rotating the entire coded sensor assembly, we can change α to different angles in the image acquisition process. We also note that α remains constant in our implementation, implying that Jpolα is spatially invariant across different (x,y) positions. This stands in contrast to commercially available polarized sensors, such as Sony IMX 250MZR, which incorporates polarizers with four different angles α over different sensor pixels.

    In the second step of the imaging process, the light field emerging from Eq. (1) is further modulated by the blood-cell layer at the coded surface plane. The resultant light field then travels a distance d2 between the coded surface and the pixel array. The captured intensity image can be expressed as Iiα(xxi,yyi)=|PSFfree(d2)*CS(xxi,yyi)Eoutα(x,y)|M2.Here, Iiα represents the ith captured image corresponding to the sensor positional shift (xi,yi). The superscript α of Iiα represents the orientation angle of the integrated coded sensor. CS(x,y) represents the transmission profile of the coded surface, and M represents M-fold downsampling by the pixel array. The multiplication and convolution operations are carried out element-wise on the respective components of the vector. Finally, the captured image is obtained by taking the square magnitude of the resulting value in Eq. (4).

    The reconstruction process of pol-CP aims to recover the retardance magnitude δ(x,y) and the slow axis orientation θ(x,y) in Eq. (2), leveraging the four sets of intensity measurements Iiα at orientation angles α=0°, 45°, 90°, 135°. In this process, we first recover the lateral positional shift (xi,yi) directly from the captured raw images [51]. We then recover the scalar wavefront at the coded surface plane Wα(x,y)=|Eoutα(x,y)|eiρα(x,y),where eiρα(x,y) represents the phase term of the scalar exit wavefront at the coded surface plane. The comprehensive recovery process for Wα(x,y) is elaborated in Appendix A. This process commences with an initialization phase that employs digital inline holographic reconstruction, followed by the coded ptychographic phase-retrieval process [52]. Upon obtaining Wα(x,y), we digitally propagate it back to the sample plane and obtain the scalar sample wavefront Sα(x,y) as follows: Sα(x,y)=Wα(x,y)*PSFfree(d1).

    Subsequently, we rotate and align Sα(x,y) with respect to S0(x,y), obtaining Rα(x,y) in Appendix A. For image rotation, we implemented a three-step shearing process in the Fourier domain [53]. For image registration, we estimate the residual transverse shift (Δx^,Δy^) and rotation angle (Δσ^) via a closed-form equation, as detailed in Appendix A. With the equations above, we arrive at the subsequent expressions for Rα(x,y): |R0(x,y)|2=t2(1sin(δ(x,y))sin(2θ(x,y))),|R45(x,y)|2=t2(1+sin(δ(x,y))cos(2θ(x,y))),|R90(x,y)|2=t2(1+sin(δ(x,y))sin(2θ(x,y))),|R135(x,y)|2=t2(1sin(δ(x,y))cos(2θ(x,y))).

    With Eqs. (7)–(10), we can recover the spatial distribution of retardance magnitude δ(x,y) and the optic axis of the birefringent sample θ(x,y) as follows: δ(x,y)=arcsinp(x,y)2+q(x,y)2,θ(x,y)=12arctan(p(x,y)q(x,y)).

    Here, p(x,y) and q(x,y) are intermediate variables, defined as p(x,y)=|R90(x,y)|2|R0(x,y)|2|R90(x,y)|2+|R0(x,y)|2,q(x,y)=|R45(x,y)|2|R135(x,y)|2|R45(x,y)|2+|R135(x,y)|2.

    To improve visualization of the reconstructed image, a pseudo-colored birefringence map can be obtained by encoding the optic axis orientation θ(x,y) as color hue and the retardance information δ(x,y) as intensity, as shown in Fig. 2(e).

    3. RESULTS

    A. Imaging Performance Characterization

    We assessed and quantified the imaging performance of the pol-CP prototype using three types of objects: a resolution target (Fig. 3); a quantitative phase target (Fig. 4); and a potato starch granules sample (Fig. 5). A calibration experiment was designed to determine the transmission profile of the blood-coded layer on the integrated sensor. For this experiment, we selected an object devoid of slow-varying phase, typically a blood smear slide, and acquired 1500 images by translating the coded sensor across various lateral positions. These images assisted in recovery of the object and the blood-coded layer on the sensor. Once the blood-cell layer is obtained, it becomes feasible to use it for subsequent experiments, reducing the number of required images to 250 to cover one field of view (41 mm2) of the detector.

    Image quality quantification using a USAF resolution target. (a) The captured raw image of the resolution target. (b) The recovered image using the pol-CP platform. The recovered intensity (c1) and phase (c2) images of the blood-cell layer on the integrated coded sensor. The goat blood cells are smallest among all animals.

    Figure 3.Image quality quantification using a USAF resolution target. (a) The captured raw image of the resolution target. (b) The recovered image using the pol-CP platform. The recovered intensity (c1) and phase (c2) images of the blood-cell layer on the integrated coded sensor. The goat blood cells are smallest among all animals.

    Validating the pol-CP prototype using a quantitative phase target. The recovered (a) amplitude, (b) phase, and (c) color-coded birefringence map.

    Figure 4.Validating the pol-CP prototype using a quantitative phase target. The recovered (a) amplitude, (b) phase, and (c) color-coded birefringence map.

    Validation of the pol-CP prototype by imaging potato starch granules. (a) The captured raw diffraction measurement of the potato starch granules. (b) The reconstructed optic axis orientation using the pol-CP prototype. (c)–(e) Zoomed-in views of the highlighted regions of (b). The line traces show the measurements using the pol-CP prototype (black dotted line) and the conventional lens-based polarized light microscopy with a 10×, 0.45 NA objective (red dotted line).

    Figure 5.Validation of the pol-CP prototype by imaging potato starch granules. (a) The captured raw diffraction measurement of the potato starch granules. (b) The reconstructed optic axis orientation using the pol-CP prototype. (c)–(e) Zoomed-in views of the highlighted regions of (b). The line traces show the measurements using the pol-CP prototype (black dotted line) and the conventional lens-based polarized light microscopy with a 10×, 0.45 NA objective (red dotted line).

    Figure 3(a) shows the captured raw image of the resolution target; Fig. 3(b) shows the recovered image using the pol-CP prototype. The 435 nm linewidth of group 10, element 2 on the target is clearly visible in the recovered image. This level of detail is particularly remarkable considering the pixel size of the detector in the captured raw image is 1.85 μm. It is worth mentioning that the achieved resolution in our demonstration of Fig. 3(b) ranks among the highest in lensless polarimetric imaging studies. The achieved resolution is also higher than the lens-based FPM demonstrations [21,22] while enabling orders of magnitude larger field of view. In Fig. 3(c), we present the recovered intensity and phase profiles of the blood-cell layer on the sensor, demonstrating the presence of a densely packed monolayer of goat blood cells. This blood-cell layer plays a crucial role in intensity and phase modulation, contributing to the successful high-resolution recovery of birefringent objects [44].

    To validate the quantitative nature of the pol-CP platform, we extended our evaluation of the pol-CP prototype by imaging a quantitative phase target. Figures 4(a) and 4(b) show the recovered amplitude and phase of the target. Figure 4(c) presents the recovered birefringence map of the phase target, confirming the absence of retardation within the sample.

    In Fig. 5, we validated the pol-CP prototype with potato starch granules, which exhibit strong birefringence due to the high radial orientation of amylopectin crystallites. Figures 5(a) and 5(b) show the raw image and the reconstructed optic axis orientation of the granules using the pol-CP prototype. Figures 5(c)–5(e) provide magnified views of the three highlighted regions in Fig. 5(b). To validate the accuracy of our pol-CP prototype, we also reconstructed the optic axis orientation using the gold standard lens-based polarized light microscopy (ground truth). In this experiment, a 10×, 0.45 NA objective was used to acquire images. The agreement between our lensless prototype (red curves) and the ground truth (black curves) in the right panels of Fig. 5(c) validates the effectiveness of the pol-CP approach.

    B. High-Resolution Gigapixel Birefringence Imaging over a Large Field-of-View

    Conventional lens-based polarized light microscopy [54] has a limited space-bandwidth product that restricts its ability to deliver both high resolution and a large field of view [13]. Our lensless pol-CP prototype offers a solution to this limitation, thus obtaining the best of both worlds. The imaging field of view of a single acquisition is limited only by the size of the image sensor. The lateral translation of the coded sensor naturally expands the field of view, enabling polarimetric whole slide imaging in high throughput. In our implementation, we typically acquire 250 images in 9  s for one field of view of the detector. By translating and rotating the coded sensor to different lateral positions and orientations, we can obtain gigapixel birefringence images with a 180 mm2 field of view in approximately 3.5 min, a performance that rivals high-end whole slide scanner but at a small fraction of the cost [15]. As an example, Figs. 6(a) and 6(b) show the raw image and recovered gigapixel birefringence map of a whole slide potato starch granule sample. Figures 6(c)–6(f) show the magnified views of four different regions, where we compare the raw measurements with the recovered amplitude, optic axis orientation, retardance, and color-coded birefringence map.

    High-resolution, large field-of-view gigapixel birefringence imaging via pol-CP. (a) The captured raw image of the whole-slide starch grain sample. (b) The recovered birefringence map. (c)–(f) The magnified views of different regions in (b). (c1)–(f1) The captured raw images. (c2)–(f2) The recovered amplitude. (c3)–(f3) The recovered optic axis orientation. (c4)–(f4) The recovered retardance. (c5)–(f5) The birefringence maps where the color scale is used to encode the optic axis orientation and the gray scale is used to encode the retardance information.

    Figure 6.High-resolution, large field-of-view gigapixel birefringence imaging via pol-CP. (a) The captured raw image of the whole-slide starch grain sample. (b) The recovered birefringence map. (c)–(f) The magnified views of different regions in (b). (c1)–(f1) The captured raw images. (c2)–(f2) The recovered amplitude. (c3)–(f3) The recovered optic axis orientation. (c4)–(f4) The recovered retardance. (c5)–(f5) The birefringence maps where the color scale is used to encode the optic axis orientation and the gray scale is used to encode the retardance information.

    C. Imaging with Quantitative Polarimetric and Phase Contrast

    Achieving good image contrast in microscopy is crucial for the study of specimen structure and enhancing diagnostic capabilities. While various methods exist to enhance image contrast in microscopy, such as fluorescence microscopy, darkfield microscopy, and differential interference contrast (DIC) microscopy, these techniques often lack the ability to provide quantitative information about the specimens without the need for external staining. For instance, fluorescence microscopy requires the application of external dyes or tags, which can potentially produce toxic byproducts or lead to photobleaching. On the other hand, darkfield and DIC microscopy can enhance image contrast, but their outputs are predominantly qualitative in nature. Unlike other common microscopy modalities, the lensless pol-CP platform offers a unique solution by generating the quantitative phase contrast and polarimetric contrast of the specimens at a high imaging throughput. The quantitative phase information is directly related to the morphology of the sample while the quantitative polarimetric contrast provides insights into the intrinsic chemical composition of the specimen.

    This capability is demonstrated through the analysis of a corn stem sample in Fig. 7. In Figs. 7(a)–7(d), we show the raw image, recovered amplitude, recovered optic axis orientation, and quantitative birefringence map of the corn stem. The plant cells in this specimen are composed of parallelly arranged cellulose structures, thereby exhibiting the expected birefringence property. In Figs. 7(e)–7(g), we further show the zoomed-in views of three highlighted regions of Fig. 7(a). Figures 7(e1)–7(g1) show the captured raw diffraction patterns, and Figs. 7(e2)–7(g2) show the recovered amplitude and phase images of the wavefront R0(x,y). We note that the recovered phase here is intrinsically tied to the polarization state. Figures 7(e3)–7(g3) show the zoomed-in views of the recovered optic axis orientation, retardance, and color-coded birefringence map. We can see that the orientations of the cell match well with the optic axis displayed in color. The retardance of the cell is much higher than the background due to the different optical and structural properties. The birefringence map, which encodes optic axis orientation and retardance information, significantly improves the contrast and makes it possible to visualize fine details of the cell that would otherwise be obscured by the background in amplitude maps. The ability to generate phase contrast and polarimetric contrast in a single imaging technique holds great value for a wide range of applications, providing researchers with new insights into the structural and chemical properties of specimens without the need for external staining. This capability opens up new possibilities for studying a wide range of samples, from biological specimens to materials science applications, with enhanced accuracy and precision.

    High contrast background-free imaging of a corn stem sample. (a) Captured raw image of corn stem. (b)–(d) The recovered amplitude, optic axis orientation, and birefringence map of (a). (e1)–(g3) Enlarged raw images of the highlighted regions in (a). (e1)–(g1) Raw images. (e2)–(g2) The recovered amplitude and phase of the wavefront at α=0°. (e3)–(g3) The reconstructed optic axis orientation, retardance, and birefringence map.

    Figure 7.High contrast background-free imaging of a corn stem sample. (a) Captured raw image of corn stem. (b)–(d) The recovered amplitude, optic axis orientation, and birefringence map of (a). (e1)–(g3) Enlarged raw images of the highlighted regions in (a). (e1)–(g1) Raw images. (e2)–(g2) The recovered amplitude and phase of the wavefront at α=0°. (e3)–(g3) The reconstructed optic axis orientation, retardance, and birefringence map.

    D. Tracking Malaria-Infected Blood Cells

    The lensless pol-CP platform can be a valuable tool for malaria diagnosis in resource-limited settings. When a malaria parasite infects a red blood cell (RBC), changes occur in its hemoglobin content, affecting the birefringence property of the cell. As the parasite grows, it crystallizes hemoglobin into hemozoin, which exhibits high birefringence. Leveraging polarimetric imaging, the pol-CP platform can effectively track and detect malaria-infected RBCs by detecting the hemozoin crystals.

    In Fig. 8, we analyze the capability of the pol-CP platform for detecting hemozoin crystals on a blood smear sample infected with Plasmodium falciparum, a particularly harmful type of malaria parasite. Figure 8(a) shows the recovered high-resolution, large-field-of-view phase map of the blood smear slide. Three zoomed-in views of Fig. 8(a) are provided in Figs. 8(b)–8(d). The raw images and reconstructed birefringence maps presented in the first two images are overlaid in the last column of Figs. 8(b)–8(d). The presence of hemozoin crystals inside infected cells allows for clear identification, enabling the visualization of the structure and birefringence of malaria-infected blood cells in high resolution. This capability of the pol-CP platform provides a promising avenue for a wide range of biomedicine-related applications, benefiting from its high-sensitivity birefringence imaging capability as well as its ability to deliver high-resolution and large field-of-view imaging.

    Detecting malaria-infected blood cells using the lensless pol-CP prototype. (a) The recovered phase image of the sample. (b)–(d) The zoomed-in views of the raw diffraction measurements (left column), quantitative birefringence maps (middle column), and overlaid images of phase and birefringence maps (right column).

    Figure 8.Detecting malaria-infected blood cells using the lensless pol-CP prototype. (a) The recovered phase image of the sample. (b)–(d) The zoomed-in views of the raw diffraction measurements (left column), quantitative birefringence maps (middle column), and overlaid images of phase and birefringence maps (right column).

    E. Polarimetric Imaging of Crystals

    Crystal imaging plays a vital role across various scientific and industrial domains, such as materials science, pharmaceuticals, and geology. The ability to visualize and analyze crystals with quantitative polarimetric contrast offers profound insights into their physical attributes, chemical composition, and multifaceted applications. One prominent example lies in the medical field’s analysis of urine samples. Urine samples often contain diverse crystal types, including calcium phosphate, calcium oxalate, uric acid, and cystine. These are influenced by factors such as an individual’s diet, hydration level, and underlying health conditions. The presence and classification of these crystals can be symptomatic of specific diseases or disorders. For example, calcium oxalate crystals may indicate a susceptibility to kidney stones, uric acid crystals might suggest gout or metabolic disorders, and cystine crystals could reveal a rare genetic condition known as cystinuria, leading to cystine stones in the kidneys, bladder, and urethra. The ability to accurately image and analyze these crystals over a large field of view and with polarimetric contrast can provide valuable diagnostic information. It enables healthcare providers to identify the specific type of crystals present, assess their concentration, and understand their formation conditions. This information can guide treatment decisions, such as dietary modifications, medication adjustments, or further diagnostic testing.

    To demonstrate the capability in crystal structure visualization, we imaged a sediment slide prepared from urine samples and identified calcium phosphate crystals, as shown in Fig. 9. The captured raw images of the crystals are shown in Fig. 9(a); the recovered intensity images at four distinct polarization angles are shown in Fig. 9(b). Figure 9(c) reveals the recovered birefringence maps of the calcium phosphate crystals. Another key advantage of the reported platform resides in its ability to perform digital refocusing after the data have been acquired. We substantiated this aspect using a calcium phosphate crystal sample, as shown in Fig. 9(b3). Visualization 2 further shows the digital propagation process of the recovered wavefront.

    Imaging of calcium phosphate crystal prepared from urine samples. (a1)–(a3) The captured raw images. (b1)–(b3) The recovered amplitude images at four polarization angles. (c1)–(c3) The recovered birefringence maps.

    Figure 9.Imaging of calcium phosphate crystal prepared from urine samples. (a1)–(a3) The captured raw images. (b1)–(b3) The recovered amplitude images at four polarization angles. (c1)–(c3) The recovered birefringence maps.

    F. Label-Free Chemical Imaging of Cytological Smears

    Tissue samples are naturally transparent under a conventional light microscope. Staining is commonly employed to enhance visibility by adding color to specific structures or features within a sample. While staining can aid in differentiation between different cell types or tissues, it is time-consuming for rapid on-site diagnostics and can potentially damage delicate structures. In contrast, our platform can provide label-free, intrinsic contrast by utilizing the polarization properties of biological specimens. To highlight this, we imaged an unstained thyroid smear using our prototype device, as shown in Fig. 10. Figures 10(a)–10(d) show the captured raw image, recovered amplitude, recovered birefringence map, and DAPI autofluorescence map of the unstained thyroid smear, respectively. In particular, Fig. 10(c) shows differential contrast arising from calcium oxalate crystals [55]. We note that these crystals are present within the colloid of normal thyroid follicles; they also contribute to the birefringence properties, as shown in purple color in Fig. 10(c). Interestingly, the autofluorescence image captured using a regular fluorescence microscope with a DAPI filter cube highlights the same areas. This is consistent with the autofluorescence emission signal in a similar range from calcium oxalates [56]. The results demonstrated here can be used to locate follicles on unstained cytology smears for rapid on-site evaluation.

    Imaging of an unstained thyroid sample using the lensless pol-CP prototype. (a) Captured raw image. (b) Recovered amplitude. (c) Recovered birefringence map, where the calcium oxalate crystals are present within the colloid of thyroid follicles and contribute to the purple color in the birefringence map. (d) Captured autofluorescence image captured using a regular fluorescence microscope with a DAPI filter cube. The bright region indicates the autofluorescence signals from calcium oxalates. The results demonstrated here can be used to locate follicles on unstained cytology smears for rapid on-site evaluation.

    Figure 10.Imaging of an unstained thyroid sample using the lensless pol-CP prototype. (a) Captured raw image. (b) Recovered amplitude. (c) Recovered birefringence map, where the calcium oxalate crystals are present within the colloid of thyroid follicles and contribute to the purple color in the birefringence map. (d) Captured autofluorescence image captured using a regular fluorescence microscope with a DAPI filter cube. The bright region indicates the autofluorescence signals from calcium oxalates. The results demonstrated here can be used to locate follicles on unstained cytology smears for rapid on-site evaluation.

    4. DISCUSSION AND CONCLUSION

    In this study, we introduced polarimetric coded ptychography (pol-CP) as a solution to the existing challenges of conventional lens-based polarization microscopy, overcoming the trade-off between spatial resolution and imaging field of view as well as the constraints in cost, spatial resolution, field of view, and imaging throughput compared with robotic microscope platforms. Our work bridges the gap between the current limitations of conventional polarimetric imaging and the need for high-throughput, high-resolution imaging systems in biomedicine and beyond. Through the integration of a coded sensor and an ultrathin polarizer film, we achieved the acquisition of gigapixel birefringence images with a 180 mm2 field of view in a time frame comparable with regular whole slide scanners. Notably, this was achieved in a cost-effective and portable manner without using any lens, rendering the platform highly suitable for field and resource-limited applications.

    In practical applications, the pol-CP platform exhibited potential in malaria diagnosis by effectively tracking and detecting malaria-infected red blood cells. The birefringence changes associated with the infection process, particularly the formation of hemozoin crystals by the parasite, were effectively identified and visualized using our platform. Moreover, we demonstrated pol-CP’s versatility in the realm of crystal imaging, particularly in the medical diagnosis of conditions indicated by the presence of specific crystals in urine samples. The platform’s ability to provide detailed structural information through birefringence mapping offers insights into the type, concentration, and formation conditions of these crystals, thereby aiding in more accurate and timely treatment decisions. Our platform’s potential was also highlighted in the label-free chemical imaging of unstained thyroid smears. Leveraging the inherent polarization properties of biological specimens, our platform offered intrinsic contrast, mitigating the need for time-consuming and potentially damaging staining processes. This ability can significantly enhance the efficiency and preservation of sample integrity in diagnostic processes, providing advancements for on-site rapid evaluations.

    Compared with conventional real-space ptychography, pol-CP employs a coded detection scheme. In our system, the coded surface serves as an “effective” probe beam for the diffracted object light waves, analogous to the structured illumination probe beam in conventional ptychography. This coded detection scheme offers several advantages, such as micrometer-scale scanning steps between adjacent acquisitions, thereby maximizing imaging throughput at the full camera bandwidth. Furthermore, it eliminates the “thin sample” constraint commonly encountered in conventional real-space ptychography. Our polarimetric coded surface is a precalibrated, integral component of the system, ensuring consistent performance across all subsequent experiments.

    When contrasted with FPM, pol-CP employs an extended laser beam that covers the entire sample slide without the photon budget issue. The millisecond short exposure time and the small scanning step size eliminate the motion blur issue, enabling continuous image acquisition at the camera’s full frame rate. In contrast, FPM often uses an LED array for illumination, which becomes a bottleneck for darkfield images requiring longer exposure times. In addition, FPM is susceptible to spatially variant pupil aberrations, particularly when using low-NA objective lenses. This results in compromised image quality toward the edge of the field of view. The reported pol-CP does not suffer from this shift-variant challenge, generally yielding better image quality.

    The pol-CP prototype platform also excels in its ability to refocus images after data acquisition, effectively addressing the autofocusing issue associated with traditional whole slide scanners. This ability enhances image quality, streamlines workflow, and minimizes the need for rescanning due to autofocusing errors. The elimination of the need to model the complex wavefront entering the sample, a constraint present in other methods such as conventional real-space ptychography and FPM, increases flexibility, particularly in working with thick specimens [24].

    The ultimate resolution of pol-CP is fundamentally constrained by the system’s maximum collection angle. In an ideal scenario where a point-object is placed proximate to the coded surface’s center and the detector’s pixel size is smaller than half the wavelength, the maximum collection angle could approach 90°, equating to a numerical aperture (NA) of 1 in air. However, real-world limitations such as the detector pixel’s sensitivity to obliquely incident light and its large size relative to the wavelength introduce constraints. These factors essentially act as low-pass filters and subsampling mechanisms, thereby limiting the achievable resolution. In our previous work, the highest NA achieved was approximately 0.8 [44]. In the current study, the NA is 0.61, possibly due to the additional distance introduced by the polarizing film between the object and the coded surface. The goat-blood layer may also be suboptical compared with the fabricated disorder-engineered surface. However, it is worth noting that the achievable NA already surpasses those demonstrated in FPM using objective lenses. To further improve the resolution, a future direction is to perform oil-immersion imaging by placing oil in-between the object and the coded surface. Another strategy is to employ an angle-varied plane wave for sample illumination and then perform aperture synthesizing as in FPM. With large-angle illumination, it is possible to achieve an NA of close to 2 while retaining the field of view of the detector.

    It is also worth noting that the use of circularly polarized light in the current prototype could introduce ambiguity in distinguishing between retardance values of 0 and π. We opted for circularly polarized illumination primarily for its ease in facilitating birefringence recovery across individual analysis channels, thereby obviating the need for intricate polarization calibration procedures. Importantly, in the realm of biological samples, retardance values seldom exceed π [57,58]. When regions with π-retardance do occur, they are typically bordered by areas with retardance values either less than or greater than π. This spatial context aids in resolving any ambiguity, akin to the phase unwrapping techniques that rely on surrounding phase values to accurately determine the correct phase wrap. Looking ahead, one avenue for future research could involve the adoption of elliptically polarized illumination, which would necessitate a comprehensive calibration procedure to fully characterize the polarization state.

    In conclusion, the pol-CP technique exhibits promise for the future of polarimetric imaging in the biomedical field. Its high-throughput imaging, high spatial resolution, low cost, portability, and flexibility with sample thickness anticipate its valuable role in disease diagnosis, sample screening, and label-free chemical imaging. We look forward to the continued development and application of this technique, bearing in mind that further work is required for its refinement, standardization, and exploration of utility in a broader range of biomedical applications.

    Acknowledgment

    Acknowledgment. X. Shao acknowledges the support of the National Natural Science Foundation of China.

    APPENDIX A: RECONSTRUCTION PROCEDURES OF POL-CP

    The reconstruction procedures of pol-CP contain three major steps: 1) iterative ptychographic phase retrieval; 2) image rotation; 3) translation and rotation refinement.

    Step 1. Ptychographic phase retrieval for object wavefront reconstructionSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElement

    Step 2. Image rotation via Fourier transform

    A general 2D clockwise rotation of the point (x,y) onto (x,y) by angle β is performed as follows: (xy)=(cosβsinβsinβcosβ)·(xy)=(1tanβ201)·(10sinβ1)·(1tanβ201)·(xy)=Mx·My·Mx·(xy),where Mx and My are the shearing matrices along the x and y directions, respectively. System.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElement

    Step 3. Translation and rotation refinement

    Accurate subpixel registration plays a pivotal role in the successful extraction of subsequent polarization information. To ensure consistency with the previous context, we define the 0° polarization as the reference and recover the positional shift and rotation angle by minimizing a cost function as follows: (Δx^,Δy^,Δσ^)=arg min(Δx,Δy,Δσ)(x,y)ω(TΔx,Δy,Δσ{h0(x,y)}hα(x,y))2,where ω denotes tracked object region. Δx,Δy, and Δσ represent the compensated shift in x direction, y direction, and rotation angle of the tracked object region, respectively. T{·} is a transform operation, which can be expressed as TΔx,Δy,Δσ{h0(x,y)}=h0(xcos(Δσ)ysin(Δσ)+Δx,ycos(Δσ)+xsin(Δσ)+Δy).

    For small values of Δσ, we assume sin(Δσ)Δσ and cos(Δσ)1. Therefore, we can obtain the following expression: TΔx,Δy,Δσ{h0(x,y)}h0(xyΔσ+Δx,y+xΔσ+Δy)h0(x,y)+(ΔxyΔσ)gx(x,y)+(Δy+xΔσ)gy(x,y),where gx(x,y)=h0(x,y)/x and gy(x,y)=h0(x,y)/y. In the third and fourth lines of the above equation, we adopt the first three terms of the Taylor series expansion as an approximation. With Eq. (A17), the Eq. (A15) can be rewritten as (Δx^,Δy^,Δσ^)=arg min(Δx,Δy,Δσ)(x,y)ω(h0(x,y)+(ΔxyΔσ)gx(x,y)+(Δy+xΔσ)gy(x,y)hα(x,y))2=arg min(Δx,Δy,Δσ)Eα(Δx,Δy,Δσ).

    By differentiating Eα(Δx,Δy,Δσ) with respect to Δx,Δy, and Δσ and setting the derivatives equal to zero, we can obtain the (Δx^,Δy^,Δσ^) to minimize Eq. (15). The estimated registration vector can be expressed as (Δx^,Δy^,Δσ^)T=M1K.Here, M and K can be expressed as follows: M=[gx2(x,y)gx(x,y)gy(x,y)g(x,y)¯gx(x,y)gx(x,y)gy(x,y)gy2(x,y)g(x,y)¯gy(x,y)g(x,y)¯gx(x,y)g(x,y)¯gy(x,y)g(x,y)¯2],K=[f(x,y)gx(x,y)f(x,y)gy(x,y)f(x,y)g(x,y)¯].

    The terms g(x,y)¯ and f(x,y) in the above equations can be expressed as g(x,y)¯=xgy(x,y)ygx(x,y),f(x,y)=hα(x,y)h0(x,y).

    Finally, the registered images Rα(x,y) can be expressed as Rα(x,y)=TΔx,Δy,Δσ{hα(x,y)}.

    References

    [1] C. He, H. He, J. Chang, B. Chen, H. Ma, M. J. Booth. Polarisation optics for biomedical and clinical applications: a review. Light Sci. Appl., 10, 194(2021).

    [2] S. Alali, A. Vitkin. Polarized light imaging in biomedicine: emerging Mueller matrix methodologies for bulk tissue assessment. J. Biomed. Opt., 20, 061104(2015).

    [3] R. Hristu, S. G. Stanciu, D. E. Tranca, G. A. Stanciu. Improved quantification of collagen anisotropy with polarization-resolved second harmonic generation microscopy. J. Biophoton., 10, 1171-1179(2017).

    [4] A. Keikhosravi, M. Shribak, M. W. Conklin, Y. Liu, B. Li, A. Loeffler, R. M. Levenson, K. W. Eliceiri. Real-time polarization microscopy of fibrillar collagen in histopathology. Sci. Rep., 11, 19063(2021).

    [5] C. W. Pirnstill, G. L. Coté. Malaria diagnosis using a mobile phone polarized microscope. Sci. Rep., 5, 13368(2015).

    [6] D. Keefe, L. Liu, W. Wang, C. Silva. Imaging meiotic spindles by polarization light microscopy: principles and applications to IVF. Reprod. Biomed. Online, 7, 24-29(2003).

    [7] N. W. McGill, P. A. Dieppe. Evidence for a promoter of urate crystal formation in gouty synovial fluid. Ann. Rheum. Dis., 50, 558-561(1991).

    [8] R. Sharma, S. Rehani, M. Mehendiratta, P. Kardam, M. Kumra, Y. Mathias, J. Yadav, K. Sahay. Architectural analysis of picrosirius red stained collagen in oral epithelial dysplasia and oral squamous cell carcinoma using polarization microscopy. J. Clin. Diagn. Res., 9, EC13-EC16(2015).

    [9] O. Warszawik-Hendzel, M. Olszewska, M. Maj, A. Rakowska, J. Czuwara, L. Rudnicka. Non-invasive diagnostic techniques in the diagnosis of squamous cell carcinoma. J. Dermatol. Case Rep., 9, 89-97(2015).

    [10] J. Rieppo, J. Hallikainen, J. S. Jurvelin, I. Kiviranta, H. J. Helminen, M. M. Hyttinen. Practical considerations in the use of polarized light microscopy in the analysis of the collagen network in articular cartilage. Microsc. Res. Tech., 71, 279-287(2008).

    [11] L.-W. Jin, K. A. Claborn, M. Kurimoto, M. A. Geday, I. Maezawa, F. Sohraby, M. Estrada, W. Kaminksy, B. Kahr. Imaging linear birefringence and dichroism in cerebral amyloid pathologies. Proc. Natl. Acad. Sci. USA, 100, 15294-15298(2003).

    [12] K. Park, T. D. Yang, D. Seo, M. G. Hyeon, T. Kong, B.-M. Kim, Y. Choi, W. Choi, Y. Choi. Jones matrix microscopy for living eukaryotic cells. ACS Photon., 8, 3042-3050(2021).

    [13] J. Park, D. J. Brady, G. Zheng, L. Tian, L. Gao. Review of bio-optical imaging systems with a high space-bandwidth product. Adv. Photon., 3, 044001(2021).

    [14] G. Zheng, C. Shen, S. Jiang, P. Song, C. Yang. Concept, implementations and applications of Fourier ptychography. Nat. Rev. Phys., 3, 207-223(2021).

    [15] Z. Bian, C. Guo, S. Jiang, J. Zhu, R. Wang, P. Song, Z. Zhang, K. Hoshino, G. Zheng. Autofocusing technologies for whole slide imaging and automated microscopy. J. Biophoton., 13, e202000227(2020).

    [16] N. Farahani, A. V. Parwani, L. Pantanowitz. Whole slide imaging in pathology: advantages, limitations, and emerging perspectives. Pathol. Lab. Med. Int., 7, 23-33(2015).

    [17] P. Ferrand, M. Allain, V. Chamard. Ptychography in anisotropic media. Opt. Lett., 40, 5144-5147(2015).

    [18] P. Ferrand, A. Baroni, M. Allain, V. Chamard. Quantitative imaging of anisotropic material properties with vectorial ptychography. Opt. Lett., 43, 763-766(2018).

    [19] A. Baroni, V. Chamard, P. Ferrand. Extending quantitative phase imaging to polarization-sensitive materials. Phys. Rev. Appl., 13, 054028(2020).

    [20] G. Zheng, R. Horstmeyer, C. Yang. Wide-field, high-resolution Fourier ptychographic microscopy. Nat. Photonics, 7, 739-745(2013).

    [21] S. Song, J. Kim, S. Hur, J. Song, C. Joo. Large-area, high-resolution birefringence imaging with polarization-sensitive Fourier ptychographic microscopy. ACS Photon., 8, 158-165(2021).

    [22] X. Dai, S. Xu, X. Yang, K. C. Zhou, C. Glass, P. C. Konda, R. Horstmeyer. Quantitative Jones matrix imaging using vectorial Fourier ptychography. Biomed. Opt. Express, 13, 1457-1470(2022).

    [23] C. Guo, Z. Bian, S. Jiang, M. Murphy, J. Zhu, R. Wang, P. Song, X. Shao, Y. Zhang, G. Zheng. OpenWSI: a low-cost, high-throughput whole slide imaging system via single-frame autofocusing and open-source hardware. Opt. Lett., 45, 260-263(2020).

    [24] T. Wang, S. Jiang, P. Song, R. Wang, L. Yang, T. Zhang, G. Zheng. Optical ptychography for biomedical imaging: recent progress and future directions. Biomed. Opt. Express, 14, 489-532(2023).

    [25] S. Dong, R. Horstmeyer, R. Shiradkar, K. Guo, X. Ou, Z. Bian, H. Xin, G. Zheng. Aperture-scanning Fourier ptychography for 3D refocusing and super-resolution macroscopic imaging. Opt. Express, 22, 13586-13599(2014).

    [26] H. Zhang, S. Jiang, J. Liao, J. Deng, J. Liu, Y. Zhang, G. Zheng. Near-field Fourier ptychography: super-resolution phase retrieval via speckle illumination. Opt. Express, 27, 7498-7512(2019).

    [27] L. Tian, L. Waller. 3D intensity and phase imaging from light field measurements in an LED array microscope. Optica, 2, 104-111(2015).

    [28] R. Horstmeyer, J. Chung, X. Ou, G. Zheng, C. Yang. Diffraction tomography with Fourier ptychography. Optica, 3, 827-835(2016).

    [29] A. M. Maiden, M. J. Humphry, J. Rodenburg. Ptychographic transmission microscopy in three dimensions using a multi-slice approach. J. Opt. Soc. Am. A, 29, 1606-1614(2012).

    [30] Z. Zhang, Y. Zhou, S. Jiang, K. Guo, K. Hoshino, J. Zhong, J. Suo, Q. Dai, G. Zheng. Invited article: mask-modulated lensless imaging with multi-angle illuminations. APL Photon., 3, 060803(2018).

    [31] C. Lu, Y. Zhou, Y. Guo, S. Jiang, Z. Zhang, G. Zheng, J. Zhong. Mask-modulated lensless imaging via translated structured illumination. Opt. Express, 29, 12491-12501(2021).

    [32] J. Kim, S. Song, B. Kim, M. Park, S. J. Oh, D. Kim, B. Cense, Y.-M. Huh, J. Y. Lee, C. Joo. Ptychographic lens-less polarization microscopy. arXiv preprint(2022).

    [33] H. Zhang, Z. Bian, S. Jiang, J. Liu, P. Song, G. Zheng. Field-portable quantitative lensless microscopy based on translated speckle illumination and sub-sampled ptychographic phase retrieval. Opt. Lett., 44, 1976-1979(2019).

    [34] P. Li, A. Maiden. Lensless LED matrix ptychographic microscope: problems and solutions. Appl. Opt., 57, 1800-1806(2018).

    [35] M. Stockmar, P. Cloetens, I. Zanette, B. Enders, M. Dierolf, F. Pfeiffer, P. Thibault. Near-field ptychography: phase retrieval for inline holography using a structured illumination. Sci. Rep., 3, 1927(2013).

    [36] Y. Zhang, S. Y. C. Lee, Y. Zhang, D. Furst, J. Fitzgerald, A. Ozcan. Wide-field imaging of birefringent synovial fluid crystals using lens-free polarized microscopy for gout diagnosis. Sci. Rep., 6, 28793(2016).

    [37] M. J. Lopera, C. Trujillo. Linear diattenuation imaging of biological samples with digital lensless holographic microscopy. Appl. Opt., 61, B77-B82(2022).

    [38] Y. Zhou, B. Xiong, X. Li, Q. Dai, X. Cao. Lensless imaging of plant samples using the cross-polarized light. Opt. Express, 28, 31611-31623(2020).

    [39] T. Liu, K. de Haan, B. Bai, Y. Rivenson, Y. Luo, H. Wang, D. Karalli, H. Fu, Y. Zhang, J. FitzGerald. Deep learning-based holographic polarization microscopy. ACS Photon., 7, 3023-3034(2020).

    [40] X. Liu, Y. Yang, L. Han, C.-S. Guo. Fiber-based lensless polarization holography for measuring Jones matrix parameters of polarization-sensitive materials. Opt. Express, 25, 7288-7299(2017).

    [41] C. Guo, S. Jiang, P. Song, T. Wang, X. Shao, Z. Zhang, G. Zheng. Quantitative multi-height phase retrieval via a coded image sensor. Biomed. Opt. Express, 12, 7173-7184(2021).

    [42] S. Jiang, C. Guo, P. Song, T. Wang, R. Wang, T. Zhang, Q. Wu, R. Pandey, G. Zheng. High-throughput digital pathology via a handheld, multiplexed, and AI-powered ptychographic whole slide scanner. Lab Chip, 22, 2657-2670(2022).

    [43] C. Zuo, J. Li, J. Sun, Y. Fan, J. Zhang, L. Lu, R. Zhang, B. Wang, L. Huang, Q. Chen. Transport of intensity equation: a tutorial. Opt. Laser Eng., 106187(2020).

    [44] S. Jiang, C. Guo, P. Song, N. Zhou, Z. Bian, J. Zhu, R. Wang, P. Dong, Z. Zhang, J. Liao, J. Yao, B. Feng, M. Murphy, G. Zheng. Resolution-enhanced parallel coded ptychography for high-throughput optical imaging. ACS Photon., 8, 3261-3271(2021).

    [45] S. Jiang, C. Guo, Z. Bian, R. Wang, J. Zhu, P. Song, P. Hu, D. Hu, Z. Zhang, K. Hoshino. Ptychographic sensor for large-scale lensless microbial monitoring with high spatiotemporal resolution. Biosens. Bioelectron., 196, 113699(2022).

    [46] S. Jiang, C. Guo, T. Wang, J. Liu, P. Song, T. Zhang, R. Wang, B. Feng, G. Zheng. Blood-coated sensor for high-throughput ptychographic cytometry on a Blu-ray disc. ACS Sens., 7, 1058-1067(2022).

    [47] C. Guo, S. Jiang, L. Yang, P. Song, A. Pirhanov, R. Wang, T. Wang, X. Shao, Q. Wu, Y. K. Cho. Depth-multiplexed ptychographic microscopy for high-throughput imaging of stacked bio-specimens on a chip. Biosens. Bioelectron., 224, 115049(2023).

    [48] S. Jiang, P. Song, T. Wang, L. Yang, R. Wang, C. Guo, B. Feng, A. Maiden, G. Zheng. Spatial- and Fourier-domain ptychography for high-throughput bio-imaging. Nat. Protoc., 18, 2051-2083(2023).

    [49] S. Jiang, J. Zhu, P. Song, C. Guo, Z. Bian, R. Wang, Y. Huang, S. Wang, H. Zhang, G. Zheng. Wide-field, high-resolution lensless on-chip microscopy via near-field blind ptychographic modulation. Lab Chip, 20, 1058-1065(2020).

    [50] P. Song, R. Wang, J. Zhu, T. Wang, Z. Bian, Z. Zhang, K. Hoshino, M. Murphy, S. Jiang, C. Guo. Super-resolved multispectral lensless microscopy via angle-tilted, wavelength-multiplexed ptychographic modulation. Opt. Lett., 45, 3486-3489(2020).

    [51] T. Wang, P. Song, S. Jiang, R. Wang, L. Yang, C. Guo, Z. Zhang, G. Zheng. Remote referencing strategy for high-resolution coded ptychographic imaging. Opt. Lett., 48, 485-488(2023).

    [52] A. Maiden, D. Johnson, P. Li. Further improvements to the ptychographical iterative engine. Optica, 4, 736-745(2017).

    [53] P. Song, C. Guo, S. Jiang, T. Wang, P. Hu, D. Hu, Z. Zhang, B. Feng, G. Zheng. Optofluidic ptychography on a chip. Lab Chip, 21, 4549-4556(2021).

    [54] R. A. Gatter. The Compensated Polarized Light Microscope in Clinical Rheumatology, 253-255(1974).

    [55] K. S. Wong, J. S. Lewis, S. Gottipati, R. D. Chernock. Utility of birefringent crystal identification by polarized light microscopy in distinguishing thyroid from parathyroid tissue on intraoperative frozen sections. Am. J. Surg. Pathol., 38, 1212-1219(2014).

    [56] E. Estève, D. Buob, F. Jamme, C. Jouanneau, S. Kascakova, J.-P. Haymann, E. Letavernier, L. Galmiche, P. Ronco, M. Daudon. Detection and localization of calcium oxalate in kidney using synchrotron deep ultraviolet fluorescence microscopy. J. Synchrotron Radiat., 29, 214-223(2022).

    [57] Y. Wang, H. He, J. Chang, C. He, S. Liu, M. Li, N. Zeng, J. Wu, H. Ma. Mueller matrix microscope: a quantitative tool to facilitate detections and fibrosis scorings of liver cirrhosis and cancer tissues. J. Biomed. Opt., 21, 071112(2016).

    [58] Y. Dong, J. Qi, H. He, C. He, S. Liu, J. Wu, D. S. Elson, H. Ma. Quantitatively characterizing the microstructural features of breast ductal carcinoma tissues in different progression stages by Mueller matrix microscope. Biomed. Opt. Express, 8, 3643-3655(2017).

    Liming Yang, Ruihai Wang, Qianhao Zhao, Pengming Song, Shaowei Jiang, Tianbo Wang, Xiaopeng Shao, Chengfei Guo, Rishikesh Pandey, Guoan Zheng. Lensless polarimetric coded ptychography for high-resolution, high-throughput gigapixel birefringence imaging on a chip[J]. Photonics Research, 2023, 11(12): 2242
    Download Citation