Open access peer-reviewed chapter

Integral Photography Technique for Three-Dimensional Imaging of Dusty Plasmas

Written By

Akio Sanpei

Submitted: 10 May 2019 Reviewed: 28 July 2019 Published: 28 August 2019

DOI: 10.5772/intechopen.88865

From the Edited Volume

Progress in Fine Particle Plasmas

Edited by Tetsu Mieno, Yasuaki Hayashi and Kun Xue

Chapter metrics overview

788 Chapter Downloads

View Full Metrics

Abstract

The integral photography technique has an advantage in which instantaneous three-dimensional (3D) information of objects can be estimated from a single-exposure picture obtained from a single viewing port. Recently, the technique has come into use for scientific research in diverse fields and has been applied to observe fine particles floating in plasma. The principle of integral photography technique and a design of a light-field camera for dusty plasma experiments are reported. The important parameters of the system, dependences of the size of the imaging area, and the spatial resolution on the number of lenses, pitch, and focal length of the lens array are calculated. Designed recording and reconstruction system is tested with target particles located on known positions and found that it works well in the range of dusty plasma experiment. By applying the integral photography technique to the obtained experimental image array, the 3D positions of dust particles floating in an RF plasma are identified.

Keywords

  • dusty plasma
  • integral photography
  • three-dimensional reconstruction
  • particle measurements
  • light-field
  • plenoptic camera

1. Introduction

Fine particles immersed in plasma are charged up negatively, show three-dimensional (3D) motion, and form 3D-ordered state, i.e., Coulomb crystal [1, 2, 3, 4, 5]. Diagnostic methods for 3D information about the positions of fine particles in a plasma have therefore been widely researched. Among the various dusty plasma experiments, 90° separated two CCD cameras with helping 3D computed tomographic reconstruction [6] and stereoscopic [7, 8] are widely used to determine the 3D position of each fine particle [9]. They require two or more detectors; however, the locations and numbers of observation ports are considerably restricted in many plasma experiment devices. Planar laser scanning technique can obtain the 3D information of particles with one CCD camera [10, 11], but it requires a little while to scan across the wide field of view. In-line holographic techniques [12] and two-color gradient methods [13, 14] can obtain 3D position of dust particles from a single-exposed photograph taken from a direction; however, these methods require a 12 bit or higher dynamic range sensors. It is required that a technique can acquire the 3D information of a dusty plasma with a single-exposed photograph taken from one viewing port with a conventional dynamic range sensor.

The “integral photography technique” [15] is known as a principle used in naked eye 3D display and in commercial refocus cameras. Such refocus camera is also called as “plenoptic camera” or “light-field camera.” It provides 3D imaging technologies based on a small lens array or a pinhole array to capture light rays from slightly different directions. This technique has an advantage in which instantaneous 3D information of objects can be estimated from a single-exposure picture obtained from a single viewing port.

Recently, the integral photography technique has come into use for scientific research in diverse fields. Example applications are particle tracking for velocimetry [16, 17, 18, 19], microscopy measurement [20, 21], spray imaging [22], etc. In the research filed of plasma, 3D reconstructions of positions of particles levitating in a plasma have been demonstrated using commercial light-field Lytro cameras [23], and the time evolution of dusty plasmas has been measured using a commercial light-field Raytrix camera [24]. An open-ended plenoptic camera, which is constructed with a lens array and a typical reflex CMOS camera, obtained the 3D positions of dust particles in a radio-frequency (RF) plasma [25, 26]. Dual-filter plenoptic imaging system has been applied to observe lithium pellets in a high-temperature plasma [27].

In this chapter, the principle, design, and experimental results of the integral photography technique for 3D imaging of dusty plasmas will be presented.

Advertisement

2. Principle of the integral photography analysis

Figure 1 shows a schematic of recording and 3D reconstruction system with integral photography for dusty plasma experiment. A small lens array is placed in front of the particles levitating in a plasma to obtain an array of projected image. The rays emerging from 3D objects, i.e., scattering light rays from dusts pass through the small lens array and are captured on a sensor device. A 3D spatial point of object that is a position of a dust particle should be projected to two-dimensional (2D) image points (X (i, j), Y (i, j)) through each (i, j)th lens of the array on a detector. The 3D reconstruction is then carried out computationally by creating inverse propagating rays within a virtual system similar to the recorded one.

Figure 1.

Schematic of recording and 3D reconstruction system with the integral photography technique: (A) object, (B) lens array, (C) array of projected images on the detector, (D) virtual lens array, and (E and E′) virtual observation planes [25].

2.1 Recording system for dusty plasma experiments

In this section, it is shown how to design a recording system for dusty plasma experiments. Figure 2 shows a schematic of relationship among lens array, detector, and imaging area on the recording system. The Cartesian (x, y, z) coordinate axes are indicated in the lower right corner of the figure. Rays are projected directly onto the sensor after passing through the lens array, which is placed at z = 0. Therefore, 3D spatial point of light source, i.e., the position of a dust particle, is projected through each lens of the array to 2D image points (X (i, j), Y (i, j)) on the detector. The array of projected images on the sensor is called as “elemental image array.”

Figure 2.

The schematic of relationship among the small lens array, the detector, and the imaging area of the recording system. Rays are projected directly onto the sensor after passing through the lens array to form an elemental image array [26].

The considerable parameters of the lens array are the number of lenses, lens pitch, and focal length. In the following discussion, we deal only with convex lens array. The focal length F of a convex lens array is calculated from the lens law as follows:

F = Ll L + l , E1

where L is the distance between the object and the lens array and l is the distance between the lens array and the sensor. Figure 3 shows F dependence on L with three values of l. The most suitable value of F can be determined by changing the experimental configuration. Particles located at misaligned positions from focused z-plane create blurred spots on the sensor.

Figure 3.

Focal length F of a convex lens array plotted as a function of the distance L between the object and the lens array for various l values between the lens and a sensor [26].

In order to make an efficient recording system, the configuration of the CMOS sensor must be considered. The number of lenses in the array should be a multiple of the aspect ratio of the CMOS sensor. In addition, the lens pitch must take into account the size of the imaging area. After passing through the lens, rays are projected onto the sensor directly. If we assume that rays from a given dust particle are projected onto all elemental images, the limit of the imaging area is calculated using straight lines passing through the center of the outermost lens

x y z = x p n y 1 / 2 0 E2

and through pixels on the edges of the sensor area corresponding to the outermost lens

x y z = x H / 2 l , x H n y 2 / 2 l . E3

Here, n is the number of lenses along a given axis, the subscript “y” denotes the direction of the y-axis, H is the length of the CMOS sensor in the y-axis, and p is the lens pitch. Furthermore, the limits for the plane perpendicular to the z-axis are expressed as

y 1 = p n y 1 2 L H n y 1 p 2 l , E4
y 2 = p n y 1 2 + L l p n y 1 2 H n y n y 2 2 . E5

In the same manner, the limits for the z-axis of the imaging area are expressed as

z 1 = n y 1 2 p × 2 n y l n y 2 H n y n y 1 p , E6
z 2 = n y 1 2 p × 2 l H n y 1 p . E7

The number nx of lenses along the x-direction is calculated from ny using the aspect ratio of the CMOS sensor. Then we can design p and n from experimental requirement of L and the size of the imaging area using above equations.

An uncertainty in the reconstructed image will be attributed to the spatial resolution of an elemental image on CMOS sensor [24]. The length h of a side of a pixel on the CMOS sensor produces uncertainties in the plane perpendicular to the z-axis as

hL l E8

and along the z-axis itself as

hL 2 l × n y 1 L + l H . E9

The above equation indicates that the large ratio of distances L/l and finite size of a pixel on the CMOS sensor cause a relatively large uncertainty along the z-axis.

2.2 Reconstruction of 3D position of light source

To extract the positions of projected particles (X (i, j), Y (i, j)) from the elemental image array, a subtraction technique is applied together with the color and profile thresholds. Background image is subtracted from obtained experimental image to reduce background noises such as the reflection of the laser light from the wall of the vacuum vessel. Then, the color space of the subtracted image is converted from RGB to “Lab” scale [28]. To obtain bright pixels colored with irradiating laser from the image, the threshold values of “L” and “a,” which depend on the experimental configurations, were considered. For example, by analyzing the dust experiment shown in the Section 3.3 (see Figure 9), pixels with “L” values greater than 110 out of 255 and “a” values less than 110 out of 255 are used. Filtering with the geometrical features of the luminance distribution patch on pixels was also applied to detect the positions of dust particles from noise. The luminance centroid of each luminance distribution patch was estimated, and it was treated as the position of the particle.

Using the observed elemental image array, the 3D image of particle distribution is reconstructed in the computer. The light path arriving at the point (pz , qz ) on a virtual observation plane located on z from a point (X (i, j), Y (i, j)) corresponding to (i, j)th lens on the elemental image array is calculated with ray tracing according to geometrical information, such as distance from lens array to detector and focal length of lenses. If the thin lens approximation can be applied, (pz , qz ) on giving z is easily calculated as

p z = X i j X center i j × z l + X center i j E10
q z = Y i j Y center i j × z l + Y center i j E11

where X center i j Y center i j is the center of the (i, j)th lens.

Luminosity I(x, y, z) of 3D light-field is estimated as a summation of light intensity I′(X (i, j), Y (i, j)) over all (i, j)th lenses [25, 29, 30] as

I x y z = i j I X i j Y i j G i j X i j Y i j cos 2 α / r 2 i j G i j X i j Y i j , E12

where r is the distance between (X (i, j), Y (i, j)) and the position (x, y, z) and α is the angle between the optical axis and the incident ray. G (i, j) is a function indicating whether (i, j)th lenslet exists in a field of view or not defined as

G i j = 1 , i f X Y lenslet 0 , otherwise E13

A target light source should locate on where I(x, y, z) shows the extremal value with respect to z, i.e., at convergent points of the rays. To detect the positions of dust particles, therefore, the virtual observation plane, indicated as E and E′ in Figure 1, is scanned along z-axis.

Advertisement

3. Experimental setup and results

3.1 Estimation of measurement parameters for a dusty plasma experiment

In order to determine the parameters of multi-convex lens array for a dusty plasma experiment, we adjusted the side of the imaging area and L to be approximately 5 and 80 mm, respectively. The commercial, standard reflex camera (D810, Nikon) was applied as the CMOS sensor. The sensor has dimensions of 36 × 24 mm2 and h = 1/205 mm. Due to the geometrical limit of the camera D810, the distance l must be more than 65 mm. Rays are projected onto the CMOS sensor directly after passing through the lens array. Subsequently, F is estimated as 35 mm from Eq. (1). The number n of lenses is inversely proportional to p and to the size of a lens. Because of the increased number of tracing rays, an increase in n increases the resolution; however, the number of pixels illuminated by each lens decreases. For a lens array including (nx , ny ) = 9 × 6 lenses, p is led as 2.2 mm and the D810 camera has 818 × 818 square pixels for each lens. The imaging area for the plane perpendicular to the z-axis lies between y 1 = −2.5 mm and y 2 = 2.423 mm, as estimated from Eqs. (4) and (5). The area of the measurable plane is a function of the position z of virtual plane, and it is maximized to be approximately 24 mm2 at z = 80 mm. Regarding the z-axis, the imaging area lies between z 1 = 143 mm and z 2 = 55 mm, as estimated from Eqs. (6) and (7). Table 1 shows the parameters of the designed multi-convex lens array and those of the recording system.

Parameter Value
F 35 mm (convex)
nx  × ny 9 × 6
Size of a lens 2.2 × 2.2 mm2
Pixels per lens 818 × 818 pixels
l ≥65 mm
Working area on xy plane at z = 80 mm ∼24 mm2
Working range of z from lens array 55–143 mm

Table 1.

Parameters of multi-convex lens array and the recording system for dusty plasma experiment with the CMOS camera D810.

A picture of the designed lens array using acrylic plastic is shown in Figure 4. The rim around the periphery facilitates the holding of the array. Figure 5 shows a sample image obtained with the designed lens array system. The object appears as a single green circle with ∼1 mm of diameter which is located on z = 80 mm.

Figure 4.

Photograph of the designed lens array. It includes 9 × 6 convex spherical lenses with rectangular boundaries. All lenses have the same F value of 35 mm [26].

Figure 5.

Sample image obtained with the designed lens array system. Single green circles are recorded as 9 × 6 circles on the CMOS sensor.

3.2 Reconstruction of known target light sources

Developed recording and reconstruction system has been tested using target light sources located on known positions. In this test experiment, the typical exposure time of the camera is 1/400 s. The elemental image array obtained with the system was stored in a computer as 10 bits of data. In addition to the position (X (i, j), Y (i, j)) of the particles identified in each elemental images, the intensity I′(X (i, j), Y (i, j)) for each particle is recorded as well. Then I(x, y, z) on a virtual observation plane located on z is calculated according to Eq. (12). The position of the maxima of the brightness of the light-field I(x, y, z) determines the z, which represents the depth of particles. After calibrating the apparent pixel sizes on the images to the real scales of the dust particle cloud, the absolute (x, y, z) coordinates can be determined.

In Figure 6, the 3D positions of the known target and reconstructed particles are marked by open squares and closed circles, respectively. In this test experiment, the optical axis is set along the z-axis, the multi-convex lens array is located at z = 0, and l is set to be 77 mm. Mechanical setting errors cause 30 μm of the error bars for the target data points. The error bars of the reconstructed data points are determined according to the pixel dimensions of the recording system. From Eqs. (8) and (9), the uncertainties for the x- (y-) and z-directions are estimated as ∼6 and 108 μm, respectively. The relative error between the positions of the target and reconstructed data fits into the known range. Therefore it is concluded that the developed recording and reconstruction system works well in the range of dusty plasma experiment.

Figure 6.

Top view of reconstruction result for the test data. The target and reconstructed data are indicated by open squares and closed circles, respectively [26].

3.3 Apply to dusty plasma experiment

Finally, the developed system is applied to a dusty plasma comprising monodiverse polymer spheres (diameter = 6.5 μm) floating in a horizontal, parallel-plate RF plasma. Figures 7 and 8 show a photograph and schematic of the experimental setup. A piezoelectric vibrator is contained in an RF electrode as the injector of fine particles into a plasma. A grounded counter electrode is positioned at the upper side of the 13.56-MHz-powered electrode at the distance of 14 mm. Fine particles levitate in the plasma generated between the electrodes. A solid-state laser, which radiates light of 532 nm in wavelength, 4 mm in diameter, and ∼10 mW in radiation power, was used in our experiment to observe fine particles in the plasma using scattered laser light.

Figure 7.

A photograph of experimental setup. Plasma is generated inside of an octagonal pillar shape chamber shown in the center of the figure. Designed lens array is located between the chamber, and the camera is shown in right-hand side of the figure.

Figure 8.

Schematic of experimental setup (not to scale). (A) 13.56 MHz RF source, (B) dust source, (C and C′) powered and grounded electrode respectively, (D) illuminating laser (10 mW at 532 nm), (E) dust particles in plasma, (F) lens array, (G) CMOS sensor.

The lens array and the CMOS sensor were located at a side port of the chamber with a distance of l = 65 mm to obtain the elemental image array. An enlarged experimentally obtained elemental image array from above experimental device is shown in Figure 9. This figure shows roughly 3 × 2 elements out of the 9 × 6 elements captured by the CMOS camera. Each image element is recorded with approximately 818 × 818 square pixels on the sensor. Scattered light from the dust particles levitating in the RF plasma appears as green dots, and slightly different elemental images result from parallax differences. As shown in Figure 9, only five dust particles floated in the plasma in this experiment. They oscillated vertically in the field of view and did not form any ordered array. The 3D positions of particles determined from Figure 9 are shown in Table 2 and Figure 10. Note that the z-axis is along the optical axis and value of z is the distance from the lens array. Figure 10(a) shows a bird’s-eye view of the reconstructed distribution. Green dots indicate the 3D positions of the levitating dusts, and cross symbols indicate the projected position on the xy plane with z = 90 mm. The xy distribution of dust particles in Figure 10(c) agrees with the observed configuration of dust particles in Figure 9. From the reconstructed image, we can recognize that particles are randomly distributed between z = 88.4 and 89.8 mm, which are the distance from the lens array. The mean distance among particles is estimated approximately 780 μm.

Figure 9.

Roughly 3 × 2 elements are enlarged from experimentally obtained array of 54 elemental images after subtraction of a background image. The green dots indicate the presence of dust particles. Each image element is recorded with approximately 818 × 818 square pixels [26].

x y z
−0.4226 0.7535 88.542
−0.2322 0.2773 88.950
0.1896 0.4043 89.676
0.3075 0.3272 88.860
0.6975 0.1367 89.540

Table 2.

x, y, and z coordinates in mm of each particle, as determined from Figure 9.

FIgure 10.

(a) Bird’s-eye view. Cross symbols indicate projected position of the x−y plane with z = 90 mm. (b) Projection of the x−z plane. (c) Projection of the x−y plane. Each panel shows the reconstructed positions of the levitating particles, as obtained from Figure 9. Dust are randomly distributed between z = 88.4 and 89.8 mm, which are the distance from the lens array [26].

Advertisement

4. Expected future of the integral photography technique for plasma measurement

The integral photography technique has great potential of versatile applications for plasma measurement. With the help of Mie-scattering ellipsometry technique [31, 32], it would bring information about the size of particles in addition to six-dimensional information about position and velocity. Combined with intrinsic fluorescence spectroscopy [33], specification of dust’s materials will be available not only for standard polymer but also for unusual target such as microorganisms [34, 35]. Moreover, deconvolution techniques [36, 37] will extend the integral photography to determine 3D distribution of spatially continuous light sources [38]. 3D information of bremsstrahlung emissivity distribution should be obtained with pinhole ultraviolet or soft X-ray detector [39, 40, 41], instead of lenslet array for visible light.

It is still unclear how many dust particles can be counted by the system using the integral photography technique because many parameters trade off against each other. Well-designed optical and recording systems are required to identify the 3D positions for a large number of particles. The defocusing effect of objects is another considerable problem. The effect makes the positions of objects on sensor difficult to identify, and uncertainties in the reconstructed image may increase. In order to avoid such a problem, some commercial light-field cameras mount different F lens arrays simultaneously. Moreover, subpixel analyses, modern particle detection, and interpolation algorithms would enable the achievement of enhanced accuracy [42, 43].

Advertisement

5. Conclusions

The integral photography technique is useful for 3D observation of dusty plasmas. This technique has an advantage in which instantaneous 3D information of objects can be estimated from a single-exposure picture obtained from a single viewing port. The principle of integral photography technique and its analytical method has been explained in detail. A design of a light-field camera for dusty plasma experiments has been reported. The important parameters of the system, dependences of the size of the imaging area, and the spatial resolution on the number of lenses, pitch, and focal length of the lens array are calculated. Then, the recording and reconstruction system has been tested with target particles located on known positions and found that it works well in the range of dusty plasma experiment. By applying the integral photography technique to the obtained experimental image array, the 3D positions of dust particles floating in an RF plasma are identified.

Advertisement

Acknowledgments

The author appreciates Prof. Y. Hayashi and Prof. Y. Awatsuji of Kyoto Institute of Technology for fruitful suggestions from the perspective of researches for the dusty plasma and the integral photography technique, respectively. The author also thanks Prof. S. Masamune of Chubu University and Prof. H. Himura of Kyoto Institute of Technology for the comments on this study. Finally, the author would like to thank Mr. K. Tokunaga of Kyoto Institute of Technology for experimental assistance. This research is partly supported by JSPS KAKENHI Grant Numbers 15K05364, 18K18750, and 24244094.

References

  1. 1. Hayashi Y, Tachibana K. Observation of Coulomb-crystal formation from carbon particles grown in a methane plasma. Japanese Journal of Applied Physics. 1994;33:L804. DOI: 10.1143/JJAP.33.L804
  2. 2. Chu JH, I L. Direct observation of Coulomb crystals and liquids in strongly coupled rf dusty plasmas. Physical Review Letters. 1994;72:4009. DOI: 10.1103/PhysRevLett. 72.4009
  3. 3. Thomas H, Morfill GE, Demmel V. Plasma crystal: Coulomb crystallization in a dusty plasma. Physical Review Letters. 1994;73:652. DOI: 10.1103/PhysRevLett.73.652
  4. 4. Piel A. Plasma crystals—Structure and dynamics. Plasma Fusion Reserch. 2009;4:013. DOI: 10.1585/pfr.4.013
  5. 5. Hayashi Y. Structural transitions of 3-dimensional Coulomb crystals in a dusty plasma. Physica Scripta. 2001;T89:112. DOI: 10.1238/Physica.Topical.089a00112/meta
  6. 6. Kak AC, Slaney M. Principles of Computerized Tomographic Imaging: Classics in Applied Mathematics. Vol. 33. Philadelphia, PA: Society for Industrial and Applied Mathematics; 2001. p. 49
  7. 7. Thomas E, Williams JD, Silver J. Application of stereoscopic particle image velocimetry to studies of transport in a dusty (complex) plasma. Physics of Plasmas. 2004;11:L37. DOI: 10.1063/1.1755705
  8. 8. Käding S, Melzer A. Three-dimensional stereoscopy of Yukawa (Coulomb) balls in dusty plasmas. Physics of Plasmas. 2006;13:090701. DOI: 10.1063/1.2354149
  9. 9. Hayashi Y. Radial ordering of fine particles in plasma under microgravity condition. Japanese Journal of Applied Physics. 2005;44:1436. DOI: 10.1143/JJAP.44.1436
  10. 10. Thomas HM, Morfill GE, Fortov VE, Ivlev AV, Molotkov VI, Lipaev AM, et al. Complex plasma laboratory PK-3 plus on the International Space Station. New Journal of Physics. 2008;10:033036. DOI: 10.1088/1367-2630/10/3/033036
  11. 11. Samsonov D, Elsaesser A, Edwards A, Thomas HM, Morfill GE. High speed laser tomography system. The Review of Scientific Instruments. 2008;79:035102. DOI: 10.1063/1.2885683
  12. 12. Kroll M, Harms S, Block D, Piel A. Digital in-line holography of dusty plasmas. Physics of Plasmas. 2008;15:063703. DOI: 10.1063/1.2932109
  13. 13. Goldbeck DD. Analyse dynamischer Volumenprozesse in komplexen Plasmen Dr. Thesis. [in German]. Fakultät für Physik, LMU München; 2003. Available from: https://edoc.ub.uni-muenchen.de/1111/
  14. 14. Annaratone BM, Antonova T, Goldbeck DD, Thomas HM, Morfill GE. Complex-plasma manipulation by radiofrequency biasing. Plasma Physics and Controlled Fusion. 2004;46:B495. Available from: https://iopscience.iop.org/article/10.1088/0741-3335/46/12B/041/meta
  15. 15. Lippmann G. PHOTOGRAPHY—Reversible prints. Integral photographs. Comptes Rendus de l’Académie des Sciences. 1908;146:446. Available from: https://people.csail.mit.edu/fredo/PUBLI/Lippmann.pdf
  16. 16. Lynch K, Fahringer T, Thurow B. Three-dimensional particle image velocimetry using a plenoptic camera. In: Proceeding of 50th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition; 2012. p. 1056. Available from: https://arc.aiaa.org/doi/pdf/10.2514/6.2012-1056
  17. 17. Cenedese A, Cenedese C, Furia F, Marchetti M, Moroni M, Shindler L. 3D particle reconstruction using light field imaging. In: Proceeding of 16th Int Symp on Applications of Laser Techniques to Fluid Mechanics; 2012. Available from: http://ltces.dem.ist.utl.pt/lxlaser/lxlaser2012/upload/253_paper_qoisoj.pdf
  18. 18. Shi S, Wang J, Ding J, Zhao Z, New TH. Parametric study on light field volumetric particle image velocimetry. Flow Measurement and Instrumentation. 2016;49:70-88. DOI: 10.1016/j.flowmeasinst.2016.05.006
  19. 19. Hall EM, Thurow BS, Guildenbecher DR. Comparison of three-dimensional particle tracking and sizing using plenoptic imaging and digital in-line holography. Applied Optics. 2016;55:6410-6420. DOI: 10.1364/AO.55.006410
  20. 20. Pégard NC, Liu H-Y, Antipa N, Gerlock M, Adesnik H, Waller L. Compressive light-field microscopy for 3D neural activity recording. Optica. 2016;3:517-524. DOI: 10.1364/OPTICA.3.000517
  21. 21. Kim J, Jeong Y, Kim H, Lee C-K, Lee B, Hong J, et al. F-number matching method in light field microscopy using an elastic micro lens array. Optics Letters. 2016;41:2751-2754. DOI: 10.1364/OL.41.002751
  22. 22. Nonn T, Jaunet V, Hellman S. Spray droplet size and velocity measurement using light-field velocimetry. In: Proceeding of ICLASS 2012, 12th Triennial International Conference on Liquid Atomization and Spray Systems; 2012. Available from: http://www.ilasseurope.org/ICLASS/iclass2012_Heidelberg/Contributions/Paper-pdfs/Contribution1219_b.pdf
  23. 23. Hartmann P, Donkó I, Donkó Z. Single exposure three-dimensional imaging of dusty plasma clusters. The Review of Scientific Instruments. 2013;84:023501. DOI: 10.1063/1.4789770
  24. 24. Jambor M, Nosenko V, Zhdanov SK, Thomas HM. Plasma crystal dynamics measured with a 3D plenoptic camera. The Review of Scientific Instruments. 2016;87:033505. DOI: 10.1063/1.4943269
  25. 25. Sanpei A, Takao N, Kato Y, Hayashi Y. Initial result of three-dimensional reconstruction of dusty plasma through integral photography technique. IEEE Transactions on Plasma Science. 2016;44:508. Available from: https://ieeexplore.ieee.org/abstract/document/7336542
  26. 26. Sanpei A, Tokunaga K, Hayashi Y. Design of an open-ended plenoptic camera for three-dimensional imaging of dusty plasmas. Japanese Journal of Applied Physics. 2017;56:080305. DOI: 10.7567/JJAP.56.080305/meta
  27. 27. Sun Z, Baldwin JK, Wei X, Wang Z, Jiansheng H, Maingi R, et al. Initial results and designs of dual-filter and plenoptic imaging for high-temperature plasmas. The Review of Scientific Instruments. 2018;89(10):E112. DOI: 10.1063/1.5036633
  28. 28. Hunter RS. Photoelectric color difference meter. Journal of the Optical Society of America. 1958;48:985-995. DOI: 10.1364/JOSA.48.000985
  29. 29. Ng R, Levoy M, Bredif M, Duval G, Horowitz M, Hanrahan P. Light field photography with a hand-held plenoptic camera. Stanford CTSR, Tech. Rep. No. 2005-02; 2005
  30. 30. Lumsdaine A, Georgiev T. The focused plenoptic camera. In: IEEE Int. Conf. Computational Photography; 2009. p. 1. DOI: 10.1109/ICCPHOT.2009.5559008
  31. 31. Hayashi Y, Kawano M, Sanpei A, Masuzaki S. Mie-scattering ellipsometry system for analysis of dust formation process in large plasma device. IEEE Transactions on Plasma Science. 2016;44:1032. Available from: https://ieeexplore.ieee.org/document/7475947
  32. 32. Hayashi Y, Sanpei A. Mie-scattering ellipsometry. In: Ellipsometry—Principles and Techniques for Materials Characterization. Rijeka, Croatia: Faustino Wahaia, IntechOpen; 2017. DOI: 10.5772/intechopen.70278
  33. 33. Lakowicz JR. Principles of Fluorescence Spectroscopy. 3rd ed. Boston, MA: Springer; 1999. DOI: 10.1007/978-0-387-46312-4
  34. 34. Sanpei A, Kigami T, Kanaya H, Hayashi Y, Sampei M. Levitation of microorganisms in the sheath of an RF plasma. IEEE Transactions on Plasma Science. 2018;46:718. Available from: https://ieeexplore.ieee.org/document/8063380
  35. 35. Sanpei A, Kigami T, Hayashi Y, Himura H, Masamune S, Sampei M. First observation of crystallike configuration of microorganisms in an RF plasma. IEEE Transactions on Plasma Science. 2019;47:3074. Available from: https://ieeexplore.ieee.org/document/8726398
  36. 36. Richardson WH. Bayesian-based iterative method of image restoration. Journal of the Optical Society of America. 1972;62:55. DOI: 10.1364/JOSA.62.000055
  37. 37. Lucy LB. An iterative technique for the rectification of observed distributions. Astronomical Journal. 1974;79:745. DOI: 10.1086/111605
  38. 38. Bolan J, Johnson KC, Thurow BS. Preliminary investigation of three-dimensional flame measurements with a plenoptic camera. In: Proceeding in 30th AIAA Aerodynamic Measurement Technology and Ground Testing Conference; 2014. DOI: 10.2514/6.2014-2520
  39. 39. Onchi T, Fujisawa A, Sanpei A. A prototype diagnostics system to detect ultraviolet emission for plasma turbulence. The Review of Scientific Instruments. 2014;85:113502. DOI: 10.1063/1.4900660
  40. 40. Onchi T, Ikezoe R, Oki K, Sanpei A, Himura H, Masamune S. Tangential soft-x ray imaging for three-dimensional structural studies in a reversed field pinch. The Review of Scientific Instruments. 2010;81:073502. DOI: 10.1063/1.3455216
  41. 41. Nishimura K, Sanpei A, Tanaka H, Ishi G, Kodera R, Ueba R, et al. 2-D electron temperature diagnostic using soft-x ray imaging technique. The Review of Scientific Instruments. 2014;85:33502. DOI: 10.1063/1.4867076
  42. 42. Feng Y, Goree J, Liu B. Accurate particle position measurement from images. The Review of Scientific Instruments. 2007;78:053704. DOI: 10.1063/1.2735920
  43. 43. Rogers SS, Waigh TA, Zhao X, Lu JR. Precise particle tracking against a complicated background: Polynomial fitting with Gaussian weight. Physical Biology. 2007;4:220. DOI: 10.1088/1478-3975/4/3/008

Written By

Akio Sanpei

Submitted: 10 May 2019 Reviewed: 28 July 2019 Published: 28 August 2019