Feature issue of digital holography and 3D imaging (DH) introduction.

The OSA Topical Meeting "Digital Holography and 3D Imaging (DH)," was held in Seattle, Washington, July 13-17, 2014. Feature issues based on the DH meeting series have been released by Applied Optics (AO) since 2007. This year Optics Express (OE) and AO jointly decided to have one such feature issue in each journal. The DH meeting will continue in the future, as expected, and the next meeting is scheduled to be held on 24 - 28 May 2015, in Shanghai Institute of Optics and Fine Mechanics, Shanghai, China.

Calculation of the scalar diffraction field from curved surfaces by decomposing the three-dimensional field into a sum of Gaussian beams.

We present a local Gaussian beam decomposition method for calculating the scalar diffraction field due to a two-dimensional field specified on a curved surface. We write the three-dimensional field as a sum of Gaussian beams that propagate toward different directions and whose waist positions are taken at discrete points on the curved surface. The discrete positions of the beam waists are obtained by sampling the curved surface such that transversal components of the positions form a regular grid. The modulated Gaussian window functions corresponding to Gaussian beams are placed on the transversal planes that pass through the discrete beam-waist position. The coefficients of the Gaussian beams are found by solving the linear system of equations where the columns of the system matrix represent the field patterns that the Gaussian beams produce on the given curved surface. As a result of using local beams in the expansion, we end up with sparse system matrices. The sparsity of the system matrices provides important advantages in terms of computational complexity and memory allocation while solving the system of linear equations.

Optical reconstruction of transparent objects with phase-only SLMs.

Three approaches for visualization of transparent micro-objects from holographic data using phase-only SLMs are described. The objects are silicon micro-lenses captured in the near infrared by means of digital holographic microscopy and a simulated weakly refracting 3D object with size in the micrometer range. In the first method, profilometric/tomographic data are retrieved from captured holograms and converted into a 3D point cloud which allows for computer generation of multi-view phase holograms using Rayleigh-Sommerfeld formulation. In the second method, the microlens is computationally placed in front of a textured object to simulate the image of the textured data as seen through the lens. In the third method, direct optical reconstruction of the micrometer object through a digital lens by modifying the phase with the Gerchberg-Saxton algorithm is achieved.

Integral imaging based 3D display of holographic data.

We propose a method and present applications of this method that converts a diffraction pattern into an elemental image set in order to display them on an integral imaging based display setup. We generate elemental images based on diffraction calculations as an alternative to commonly used ray tracing methods. Ray tracing methods do not accommodate the interference and diffraction phenomena. Our proposed method enables us to obtain elemental images from a holographic recording of a 3D object/scene. The diffraction pattern can be either numerically generated data or digitally acquired optical data. The method shows the connection between a hologram (diffraction pattern) and an elemental image set of the same 3D object. We showed three examples, one of which is the digitally captured optical diffraction tomography data of an epithelium cell. We obtained optical reconstructions with our integral imaging display setup where we used a digital lenslet array. We also obtained numerical reconstructions, again by using the diffraction calculations, for comparison. The digital and optical reconstruction results are in good agreement.

Scalar diffraction field calculation from curved surfaces via Gaussian beam decomposition.

We introduce a local signal decomposition method for the analysis of three-dimensional (3D) diffraction fields involving curved surfaces. We decompose a given field on a two-dimensional curved surface into a sum of properly shifted and modulated Gaussian-shaped elementary signals. Then we write the 3D diffraction field as a sum of Gaussian beams, each of which corresponds to a modulated Gaussian window function on the curved surface. The Gaussian beams are propagated according to a derived approximate expression that is based on the Rayleigh-Sommerfeld diffraction model. We assume that the given curved surface is smooth enough that the Gaussian window functions on it can be treated as written on planar patches. For the surfaces that satisfy this assumption, the simulation results show that the proposed method produces quite accurate 3D field solutions.

Full-complex amplitude modulation with binary spatial light modulators.

Imperfections and nonrobust behavior of practical multilevel spatial light modulators (SLMs) degrade the performance of many proposed full-complex amplitude modulation schemes. We consider the use of more robust binary SLMs for this purpose. We propose a generic method, by which, out of K binary (or 1 bit) SLMs of size M×N, we effectively create a new 2(K)-level (or K bit) SLM of size M×N. The method is a generalization of the well-known concepts of bit plane representation and decomposition for ordinary gray scale digital images and relies on forming a properly weighted superposition of binary SLMs. When K is sufficiently large, the effective SLM can be regarded as a full-complex one. Our method is as efficient as possible from an information theoretical perspective. A 4f system is discussed as a possible optical implementation. This 4f system also provides a means for eliminating the undesirable higher diffraction orders. The components of the 4f system can easily be customized for different production technologies.

Integral imaging using phase-only LCoS spatial light modulators as Fresnel lenslet arrays.

We present a digital integral imaging system. A Fresnel lenslet array pattern is written on a phase-only LCoS spatial light modulator device (SLM) to replace the regular analog lenslet array in a conventional integral imaging system. We theoretically analyze the capture part of the proposed system based on Fresnel wave propagation formulation. Because of pixelation and quantization of the lenslet array pattern, higher diffraction orders and multiple focal points emerge. Because of the multiple focal planes introduced by the discrete lenslets, multiple image planes are observed. The use of discrete lenslet arrays also causes some other artifacts on the recorded elemental images. The results reduce to those available in the literature when the effects introduced by the discrete nature of the lenslets are omitted. We performed simulations of the capture part. It is possible to obtain the elemental images with an acceptable visual quality. We also constructed an optical integral imaging system with both capture and display parts using the proposed discrete Fresnel lenslet array written on a SLM. Optical results when self-luminous objects, such as an LED array, are used indicate that the proposed system yields satisfactory results.

Circular holographic video display system.

A circular holographic video display system reconstructs holographic video. Phase-only spatial light modulators are tiled in a circular configuration in order to increase the field of view. A beam-splitter is used to align the active area of the SLMs side by side without any gap. With the help of this configuration observers can see 3D ghost-like image floating in space and can move and rotate around the object. The 3D reconstructions can be observed binocularly. Experimental results are satisfactory.

Synthesis of three-dimensional light fields with binary spatial light modulators.

Computation of a binary spatial light modulator (SLM) pattern that generates a desired light field is a challenging quantization problem for which several algorithms have been proposed, mainly for far-field or Fourier plane reconstructions. We study this problem assuming that the desired light field is synthesized within a volumetric region in the non-far-field range after free space propagation from the SLM plane. We use Fresnel and Rayleigh-Sommerfeld scalar diffraction theories for propagation of light. We show that, when the desired field is confined to a sufficiently narrow region of space, the ideal gray-level complex-valued SLM pattern generating it becomes sufficiently low pass (oversampled) so it can be successfully halftoned into a binary SLM pattern by solving two decoupled real-valued constrained halftoning problems. Our simulation results indicate that, when the synthesis region is considered, the binary SLM is indistinguishable from a lower resolution full complex gray-level SLM. In our approach, free space propagation related computations are done only once at the beginning, and the rest of the computation time is spent on carrying out standard image halftoning.

Exact solution for scalar diffraction between tilted and translated planes using impulse functions over a surface.

The diffraction relation between a plane and another plane that is both tilted and translated with respect to the first one is revisited. The derivation of the result becomes easier when the impulse function over a surface is used as a tool. Such an approach converts the original 2D problem to an intermediate 3D problem and thus allows utilization of easy-to-interpret Fourier transform properties due to rotation and translation. An exact solution for the scalar monochromatic propagating waves case when the propagation direction is restricted to be in the forward direction is presented.

Lasers: the first fifty years.

This year marks the 50th anniversary of the invention of the laser. The Optical Society of America is publishing this feature issue to celebrate this auspicious birthday.

Graphics processing unit accelerated computation of digital holograms.

An approximation for fast digital hologram generation is implemented on a central processing unit (CPU), a graphics processing unit (GPU), and a multi-GPU computational platform. The computational performance of the method on each platform is measured and compared. The computational speed on the GPU platform is much faster than on a CPU, and the algorithm could be further accelerated on a multi-GPU platform. In addition, the accuracy of the algorithm for single- and double-precision arithmetic is evaluated. The quality of the reconstruction from the algorithm using single-precision arithmetic is comparable with the quality from the double-precision arithmetic, and thus the implementation using single-precision arithmetic on a multi-GPU platform can be used for holographic video displays.

Real-time phase-only color holographic video display system using LED illumination.

A real-time full-color phase-only holographic display system generates holograms of 3D objects. The system includes a 3D object formed by voxels, an internet-based transmission capability that transmits the object information to the server, a real-time hologram generation unit, and a holographic display unit with incoherent illumination. The server calculates three phase holograms for RGB components using multiple GPUs. The resultant phase holograms are saved into an RGB bitmap image and loaded to the phase-only spatial light modulators (SLMs). SLMs are illuminated uniformly by LEDs, and reconstructed waves are aligned and overlapped by using high precision optics and stages. Experimental results are satisfactory.

Hologram synthesis for photorealistic reconstruction.

Computation of diffraction patterns, and thus holograms, of scenes with photorealistic properties is a highly complicated and demanding process. An algorithm, based primarily on computer graphics methods, for computing full-parallax diffraction patterns of complicated surfaces with realistic texture and reflectivity properties is proposed and tested. The algorithm is implemented on single-CPU, multiple-CPU and GPU platforms. An alternative algorithm, which implements reduced occlusion diffraction patterns for much faster but somewhat lower quality results, is also developed and tested. The algorithms allow GPU-aided calculations and easy parallelization. Both numerical and optical reconstructions are conducted. The results indicate that the presented algorithms compute diffraction patterns that provide successful photorealistic reconstructions; the computation times are acceptable especially on the GPU implementations.

Exact analysis of the effects of sampling of the scalar diffraction field.

If the sampled diffraction pattern due to a planar object is used to reconstruct the object pattern by backpropagation, the obtained pattern is no longer the same as the original. The effect of such sampling on the reconstruction is analyzed. The formulation uses the plane-wave expansion, and therefore the provided solution is exact for wave propagation in media where scalar wave propagation is valid. In contrast to the sampling effects under the Fresnel approximation, the exact solution indicates that there are no modulated replicas of the original object in the reconstructed pattern. Rather, the distortion is in the form of modulated, translated, and dispersed versions of the original.