Simulation study of phase retrieval for hard X-ray in-line phase contrast imaging
|
|
- Jeffrey Gaines
- 6 years ago
- Views:
Transcription
1 450 Science in China Ser. G Physics, Mechanics & Astronomy 2005 Vol.48 No Simulation study of phase retrieval for hard X-ray in-line phase contrast imaging YU Bin 1,2, PENG Xiang 1,2, TIAN Jindong 1, NIU Hanben 1, DIAO Luhong 3 & LI Hua 3 1. Key Laboratory of Optoelectronic Devices and Systems of Ministry of Education, Institute of Optoelectronics, Shenzhen University, Shenzhen , China; 2. State Key Laboratory of Precision Measuring Technology and Instruments, Tianjin University, Tianjin ,China; 3. Key Laboratory of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing ,China Correspondence should be addressed to Yu Bin ( yubin@szu.edu.cn) Received March 21, 2005; revised April 2, 2005 Abstract Two algorithms for the phase retrieval of hard X-ray in-line phase contrast imaging are presented. One is referred to as Iterative Angular Spectrum Algorithm (IASA) and the other is a hybrid algorithm that combines IASA with TIE (transport of intensity equation). The calculations of the algorithms are based on free space propagation of the angular spectrum. The new approaches are demonstrated with numerical simulations. Comparisons with other phase retrieval algorithms are also performed. It is shown that the phase retrieval method combining the IASA and TIE is a promising technique for the application of hard X-ray phase contrast imaging. Keywords: X-ray optics, phase contrast imaging, phase retrieval, angular spectrum, hard X-ray. DOI: / X-ray phase contrast imaging (XPCI) is a novel method that exploits the phase shift for the incident X-ray to form an image. For light elements such as carbon, hydrogen and oxygen, the phase-shift term can be up to 1000 times greater than the absorption term in the hard X-ray energy region. So XPCI has attracted much attention in recent years. Various methods for XPCI have been proposed and demonstrated on synchrotron devices and other X-ray sources [1 13], particularly the in-line method [7 8] proposed by Wilkins et al., which simply exploits Fresnel diffraction to provide contrast. The in-line hard X-ray phase contrast imaging has become one of the promising XPCI techniques. A key problem of X-ray phase contrast imaging is how to retrieve the phase from the image intensity. Nugent et al. established a method based on the transport of intensity equation (TIE) [9 13], which can be applied successfully to retrieve the phase information in near-field region. Cheng J. et al. [4] also presented a general treatment of X-ray image
2 Simulation study of phase retrieval for hard X-ray in-line phase contrast imaging 451 formation by direct Fresnel diffraction with partially coherent hard X- rays. In this article, two algorithms for the phase retrieval of hard X-ray in-line phase contrast imaging are presented. One is referred to as Iterative Angular Spectrum Algorithm (IASA) and the other is a hybrid algorithm that combines IASA with TIE. The new approaches are demonstrated with numerical simulations. Comparisons with other phase retrieval algorithms are also performed. It is shown that the phase retrieval method combining the IASA and TIE is a promising technique for the application of hard X-ray phase contrast imaging. 1 Algorithm for phase retrieval The simplified schematic setup for in-line hard X-ray phase contrast imaging is shown in Fig. 1. We consider the imaging setup as a plane monochromatic incident wave propagating along the direction of the optic axis z. Fig. 1. Schematic set-up of hard X-ray in-line phase contrast imaging. According to angular spectrum of plane-waves diffraction theory [14], the complex wave field U ( x, y ) ρ ( x, y ) exp i φ ( x, y ) = in the image plane can be described as the result of propagation of the complex wave field (, ) ρ ( x, y ) exp i φ ( x, y ) U x y = in the object plane to the distance z, with the Fourier transforms of these two (initial and propagated) wave functions given by A 1 (α, β, 0) and A 2 (α, β, z) correspondingly. In the Fresnel approximation, the Fourier transform of the diffracted wave function is related to the Fourier transform of the initial function via the frequency transfer characteristic of the free space H(α, β, z), given by: ( αβ,, z) A H = = A exp ( ) 1 ikz 1 α β ( αβ,,0) ( αβ,,0) = (,,0exp ) { i ( α + β )} (1a) (1b) A U x y k x y dx dy
3 452 Science in China Ser. G Physics, Mechanics & Astronomy 2005 Vol.48 No ( αβ ) = ( ) { ( α + β )} 2 2 (1c) A,, z U x, y, z exp i k x y dx dy, 2π where α, β is the x and y direction cosine, k = is wave number, λ is the radiation λ wavelength. The intensity distribution in the image plane is proportional to the intensity of this field as ( ) ( ) 2 I x, y = U x, y. (2) Thus the problem of numerical reconstruction of X-ray phase contrast imaging can be described as follows: if we know the data of intensity measurements, I 2 (x 2, y 2 ), collected U x, y, in one or more planes orthogonal to the optic axis z, and the amplitude ρ 1 of ( ) how can we recover the phase φ 1 of (, ) U x y where the recovered wave functions satisfy expression (1) with high accuracy? Expressions (1) provide a symmetrical relation between the initial and diffracted wave functions in the Fresnel approximation. So we can reconstruct the phase information based on the iterative angular spectrum approach (IASA) proposed by Mellin and Nordin [15]. The IASA algorithm is based on Gerchberg-Saxton (GS) [16 18] and applied in the order (1b) (1a) (1c), resulting in the diffracted wave function; while being applied in the reversed order, they allow for reconstruction of the initial wave function from the result of diffraction. We shall denote the forward and the reversed propagation operations defined by expressions (1a), (1b), (1c) respectively. The algorithm described above can be implemented numerically by using Fast Fourier Transform (FFT) on a finite rectangular grid under periodic border conditions. Thus an iterative phase retrieval algorithm from two intensity planes in the Fresnel domain is simply obtained by using free space calculation, as shown in Fig. 2, instead of a single FFT block. Fig. 2. Free space propagation using angular spectrum calculation. Fig. 3 shows the flow chart of the iterative procedure based IASA algorithm. First, the phase object can be used as the initial value of φ 1 to start the iteration. The algorithm starts with the grabbed object plane and a random phase. The obtained wavefront is then transformed to the second plane using free space propagation calculation. Secondly, the phase at this point is maintained unchanged while the amplitude is forced to be the one grabbed at the image plane. The obtained wavefront is transformed back using free space propagation backwards where the obtained phase is again maintained unchanged and the amplitude is forced to be the one grabbed at the object plane. The algorithm is
4 Simulation study of phase retrieval for hard X-ray in-line phase contrast imaging 453 Fig. 3. IASA algorithm used for phase retrieval. running back and forth between the object and image planes until the sum-squared error (SSE) between the grabbed and the calculated amplitudes, defined by 2 2 (3) ( n) SSE = ρ ρ ρ ( n) where ρ2 expresses the result of the nth iteration for initial guess phase φ 1, reaches a given small quantity for determining the calculation accuracy, or until the number of iterations exceeds a given value N max. At this point, the calculated phase of the object and image planes has converged towards the values of the physical phase. 2 Numerical Simulation In this section we solve the phase-and-amplitude retrieval problem, which occurs in quantitative in-line hard X-ray phase contrast imaging. In the following calculations we use the simulated object-plane phase distribution shown in Fig. 4.The well-known Lena image shown here has a variety of features on different resolution scales that tend to be quite helpful for visual assessment of the reconstruction accuracy. The image had pixels and was assigned a physical size of µm. The radiation wavelength was set to λ = 0.1 nm which corresponds to hard X-rays. The phase values in Fig. 4 varied in the interval [ 1, 0] rad. These conditions are realistic for quantitative in-line hard X-ray phase contrast imaging. The object-plane intensity distribution was uniform, I 1 (x 1, y 1 ) 1, which corresponds to the case of a plane incident Fig. 4. Simulated phase distribution in the wave and a non- absorbing object. object plane.
5 454 Science in China Ser. G Physics, Mechanics & Astronomy 2005 Vol.48 No We calculated the intensity distribution, (, ) I x y, in the image planes for different propagation distances z = 0.1, 1 and 10 m (Fig. 5), which was calculated by numerically evaluating expression (1). Fig. 5. Intensity distributions in the different image plans: (a) z = 0.1 m; (b) z = 1 m; (c) z = 10 m. From Fig. 5, we can see that the image is an edge-enhanced image in near-field region. Then we used the TIE to recover the object phase φ 1 (Figs. 6 (a) (c)) from the image intensities (, ) I x y at z = 0.1, 1 and 10 m. Obviously, the reconstruction quality is very good at shorter distances, and quite bad at the longest one where only the low spatial frequencies have been reconstructed successfully. Fig. 6. Reconstructed phase with the TIE algorithm from intensity distributions in Fig. 5. Next we applied the GS and IASA algorithms for phase reconstruction from the im- I x, y at z = 0.1, 1 and 10 m, respectively. We applied the two algo- age intensities ( ) rithms with N =100 iterations between the object and the image planes, and with the random phases used as an initial approximation. As assumed the experimental intensity distributions in the object and image planes is known, we forced the intensity distribu-, I x, y in the object and image planes respec- tions to be equal to I ( x y ) and ( ) tively at each iteration of the algorithm. The results of the phase retrievals using the two new methods described above are presented in Figs. 7 8
6 Simulation study of phase retrieval for hard X-ray in-line phase contrast imaging 455 Fig. 7. Reconstructed phase with the GS algorithm from intensity distributions in Fig. 5. Fig. 8. Reconstructed phase with the IASA algorithm from intensity distributions in Fig. 5. From above images, at the two longer distances the two methods produced the better results. But the results produced by IASA algorithm were better than by GS algorithm, because the IASA method is based on angular spectrum propagation. In view of expression (1), the transfer function has no approximation, so the IASA method can be used in any field. Nevertheless, because the random of originally guess phase and iteration of the algorithm, the reconstructed results had some errors. We can improve it by combining the iterative algorithm with TIE method or adding iterative numbers. Then we applied IASA or GS algorithm with N = 100 iterations between the object and the image planes, and with the corresponding TIE-reconstructed phases used as an initial approximation (IASA+TIE algorithm or IASA+TIE algorithm). The results of the phase retrievals using the two new methods described above are presented in Figs One can see that the IASA+TIE algorithm produced the best results at shorter distances, but all the results show some improvement over the original TIE-reconstructed phase, GS-reconstructed phase and the original IASA-reconstructed phase. A quantitative estimation of the phase retrieval accuracy can be calculated with the following relative root-mean-square (RMS) error formula:
7 456 Science in China Ser. G Physics, Mechanics & Astronomy 2005 Vol.48 No Fig. 9. Reconstructed phase with the GS+TIE algorithm from intensity distributions in Fig. 5 Fig. 10. Reconstructed phase with the IASA+TIE algorithm from intensity distributions in Fig. 5 where RMS = 100% φ φ φ rec true 2 true 2 ij ij ij i, j i, j (4) true φ ij was the value of the phase at pixel (i, j) of the original simulated phase distribution and rec φ ij was the reconstructed phase value at the same pixel. The results of phase retrieval are presented in Table 1 for different values of the object-to-detector distances. These error values are in good agreement with the visual impression of the reconstruction quality. Table 1 Phase reconstruction errors RMS for varying object-to-image distance z Method/distance z = 0.1 m z = 1 m z = 10 m TIE 8.03% 15.41% 24.90% GS 35.73% 31.46% 21.92% IASA 35.06% 29.41% 19.28% GS+TIE 2.58% 5.05% 9.52% IASA+TIE 2.45% 4.24% 6.70% 1 2
8 Simulation study of phase retrieval for hard X-ray in-line phase contrast imaging 457 Fig. 11 illustrates the relationship between the SSE and the number of iterations for four iterative algorithms at z = 1 m. As only 100 iterations are used, the SSE are for IASA+TIE algorithm. It can be seen that the SSE is reduced faster when the IASA+TIE algorithm is used. We spectulate that better results could be obtained with the IASA+TIE algorithm if more iteration were used. Fig. 11. Relationship between the SSE and the number of iterations at z =1 m. 3 Conclusions We have demonstrated that IASA algorithm can be used to phase retrieval. It is also shown that the new method combining the TIE method has improved the quality of phase retrieval from hard X-ray in-line phase contrast imaging. The algorithms delivered substantially more accurate retrieved phase and amplitude compared to the results of the IASA and TIE methods. Numerical examples show the effectiveness of our methods. Acknowledgements This work was supported by the National Natural Science Foundation of China (Grant No ), and the Natural Science Foundation of Guangdong Province (Grant Nos and ), and Shenzhen Bureau of Science, Technology & Information (Grant No ). References 1. Fitzgerald, R., Phase sensitive X-ray imaging, Phys. Today, 2000, 53(7): Gao, D. C., Pogany, A., Stevenson, A.W. et al., Hard X-ray phase contrast imaging, Acta Phys. Sin., 2000, 49: Chen Min, Xiao Tiqiao, Luo yuyu et al., Phase contrast imaging with microfocus X-ray source, Acta Phys. Sin., 2004, 53:
9 458 Science in China Ser. G Physics, Mechanics & Astronomy 2005 Vol.48 No Cheng Jing, Han Shensheng, Shao Wenwen et al., Study of partial coherent X-ray phase imaging, Acta Optica Sinica, 1999, 19(5): Cheng, J., Han, S., Phase imaging with partially coherent X rays, Opt. Lett., 1999, 24(4): Yu Hong, Zhu Pinpin, Han Shensheng, A study of diffractive phase imaging and phase restoration with partially coherent X-ray, Acta Optica Sinica, 2003, 23 (4): Wilkins, S. W., Gureyev, T. E., Gao, D. et al., Phase-contrast imaging using polychromatic hard X-rays, Nature, 1996, 384: Pogany, A., Gao, D., Wilkins, S. W., Contrast and resolution in imaging with a microfocus X-ray source, Rev. Sci. Insirum., 1997, 68(7): Nugent, K. A., Gureyev, T. E., Cookson, D. F. et al., Quantitative phase imaging using hard X-rays, Phys. Rev. Lett., 1996, 77(14): Gureyev, T. E., Wilkins, S. W., On X-ray phase imaging with a point source, J. Opt. Soc. Am. A, 1998, 15(3): Gureyev, T. E., Composite techniques for phase retrieval in the Fresnel region, Optics Communications, 2003, 220: Gureyev, T. E., Pogany, A., Paganin, D. M. et al., Linear algorithms for phase retrieval in the Fresnel region, Optics Communications, 2004, 231: Gureyev, T. E., Nugent, K. A., Rapid quantitative phase imaging using the transport of intensity equation, Optics Communications, 1997, 133: Goodman, J., Introduction to Fourier Optics (New York: McGraw-Hill), 1968, Stephen, D., Mellin, Gregory, P., Nordin, Limits of scalar diffraction theory and an iterative angular spectrum algorithm for finite aperture diffractive optical element design, Opt. Express, 2001, 8(13): Gerchberg, R. W., Saxton, W. O., A practical algorithm for the determination of phase from image and diffraction plane pictures, Optik (Stuttgart), 1972, 35: Fienup, J. R., Reconstruction of an object from the modulus of its Fourier transform, Opt. Lett., 1978, 3(1): Fienup, J. R., Phase retrieval algorithms: a comparison, Appl. Opt., 1982, 21:
Fourier, Fresnel and Image CGHs of three-dimensional objects observed from many different projections
Fourier, Fresnel and Image CGHs of three-dimensional objects observed from many different projections David Abookasis and Joseph Rosen Ben-Gurion University of the Negev Department of Electrical and Computer
More informationPhase problem and the Radon transform
Phase problem and the Radon transform Andrei V. Bronnikov Bronnikov Algorithms The Netherlands The Radon transform and applications Inverse problem of phase-contrast CT Fundamental theorem Image reconstruction
More informationSupplementary Figure 1: Schematic of the nanorod-scattered wave along the +z. direction.
Supplementary Figure 1: Schematic of the nanorod-scattered wave along the +z direction. Supplementary Figure 2: The nanorod functions as a half-wave plate. The fast axis of the waveplate is parallel to
More informationPerformance analysis of ghost imaging lidar in background light environment
Research Article Vol. 5, No. 5 / October 2017 / Photonics Research 431 Performance analysis of ghost imaging lidar in background light environment CHENJIN DENG, 1,2 LONG PAN, 1,2 CHENGLONG WANG, 1,2 XIN
More informationImaging outside the box: Resolution enhancement in X-ray coherent diffraction imaging by extrapolation of diffraction patterns
Zurich Open Repository and Archive University of Zurich Main Library Strickhofstrasse 39 CH-8057 Zurich www.zora.uzh.ch Year: 2015 Imaging outside the box: Resolution enhancement in X-ray coherent diffraction
More informationLimits of computational white-light holography
Journal of Physics: Conference Series Limits of computational white-light holography To cite this article: Sebastian Mader et al 2013 J. Phys.: Conf. Ser. 415 012046 View the article online for updates
More informationDigital correlation hologram implemented on optical correlator
Digital correlation hologram implemented on optical correlator David Abookasis and Joseph Rosen Ben-Gurion University of the Negev Department of Electrical and Computer Engineering P. O. Box 653, Beer-Sheva
More informationSolution to the Twin Image Problem in Holography
1 Solution to the Twin Image Problem in Holography Tatiana Latychevskaia & Hans-Werner Fink Institute of Physics, University of Zurich, Winterthurerstrasse 190, CH-8057, Switzerland While the invention
More information1 Laboratory #4: Division-of-Wavefront Interference
1051-455-0073, Physical Optics 1 Laboratory #4: Division-of-Wavefront Interference 1.1 Theory Recent labs on optical imaging systems have used the concept of light as a ray in goemetrical optics to model
More informationUniversal Iterative Phasing Method for Near-Field and Far-Field Field Coherent Diffraction Imaging
Universal Iterative Phasing Method for Near-Field and Far-Field Field Coherent Diffraction Imaging Qun Shen (APS, ANL) Xianghui Xiao (Cornell) Imaging in different regimes near-field field phase contrast
More informationTo determine the wavelength of laser light using single slit diffraction
9 To determine the wavelength of laser light using single slit diffraction pattern 91 Apparatus: Helium-Neon laser or diode laser, a single slit with adjustable aperture width, optical detector and power
More informationFormula for the asymmetric diffraction peak profiles based on double Soller slit geometry
REVIEW OF SCIENTIFIC INSTRUMENTS VOLUME 69, NUMBER 6 JUNE 1998 Formula for the asymmetric diffraction peak profiles based on double Soller slit geometry Takashi Ida Department of Material Science, Faculty
More informationLED holographic imaging by spatial-domain diffraction computation of. textured models
LED holographic imaging by spatial-domain diffraction computation of textured models Ding-Chen Chen, Xiao-Ning Pang, Yi-Cong Ding, Yi-Gui Chen, and Jian-Wen Dong* School of Physics and Engineering, and
More informationModule 18: Diffraction-I Lecture 18: Diffraction-I
Module 18: iffraction-i Lecture 18: iffraction-i Our discussion of interference in the previous chapter considered the superposition of two waves. The discussion can be generalized to a situation where
More informationCondenser Optics for Dark Field X-Ray Microscopy
Condenser Optics for Dark Field X-Ray Microscopy S. J. Pfauntsch, A. G. Michette, C. J. Buckley Centre for X-Ray Science, Department of Physics, King s College London, Strand, London WC2R 2LS, UK Abstract.
More informationMEASUREMENT OF THE WAVELENGTH WITH APPLICATION OF A DIFFRACTION GRATING AND A SPECTROMETER
Warsaw University of Technology Faculty of Physics Physics Laboratory I P Irma Śledzińska 4 MEASUREMENT OF THE WAVELENGTH WITH APPLICATION OF A DIFFRACTION GRATING AND A SPECTROMETER 1. Fundamentals Electromagnetic
More informationMeasurement of period difference in grating pair based on analysis of grating phase shift
Measurement of period difference in grating pair based on analysis of grating phase shift Chao Guo, Lijiang Zeng State Key Laboratory of Precision Measurement Technology and Instruments Department of Precision
More informationHolographic elements for Fourier transform
Optica Applicata, Vol. XXXIV, No. 1, 2004 Holographic elements for Fourier transform EUGENIUSZ JAGOSZEWSKI, ANDRZEJ ANDRUCHÓW Institute of Physics, Wrocław University of Technology, Wybrzeże Wyspiańskiego
More informationCoherent Diffraction Imaging with Nano- and Microbeams
Diffraction Imaging with Nano- and Microbeams Why does lensless need? Mark A Pfeifer Cornell High Energy Synchrotron Source Cornell University Ithaca, NY 14850 map322@cornell.edu XLD 2011 June 28, 2011
More informationTHz Transmission Properties of Metallic Slit Array
THz Transmission Properties of Metallic Slit Array Guozhong Zhao * Department of Physics, Capital Normal University, Beijing 100048, China Beijing Key Lab for Terahertz Spectroscopy and Imaging Key Lab
More informationContrast Optimization: A faster and better technique for optimizing on MTF ABSTRACT Keywords: INTRODUCTION THEORY
Contrast Optimization: A faster and better technique for optimizing on MTF Ken Moore, Erin Elliott, Mark Nicholson, Chris Normanshire, Shawn Gay, Jade Aiona Zemax, LLC ABSTRACT Our new Contrast Optimization
More informationNarcissus of Diffractive OpticalSurfaces. Jonathan B. Cohen. ElOp E1ectrooptics Industries Ltd. P.O.B , Rehovot , Israel LINTRODUCTION
Narcissus of Diffractive OpticalSurfaces Jonathan B. Cohen ElOp E1ectrooptics Industries Ltd. P.O.B. 1 165, Rehovot 7611 1, Israel 1. ABSTRACT Narcissus is usually approximated by means of a paraxial ray
More informationDetermining Wave-Optics Mesh Parameters for Modeling Complex Systems of Simple Optics
Determining Wave-Optics Mesh Parameters for Modeling Complex Systems of Simple Optics Dr. Justin D. Mansell, Steve Coy, Liyang Xu, Anthony Seward, and Robert Praus MZA Associates Corporation Outline Introduction
More informationExperimental Observation of Invariance of Spectral Degree of Coherence. with Change in Bandwidth of Light
Experimental Observation of Invariance of Spectral Degree of Coherence with Change in Bandwidth of Light Bhaskar Kanseri* and Hem Chandra Kandpal Optical Radiation Standards, National Physical Laboratory,
More informationOptical simulations within and beyond the paraxial limit
Optical simulations within and beyond the paraxial limit Daniel Brown, Charlotte Bond and Andreas Freise University of Birmingham 1 Simulating realistic optics We need to know how to accurately calculate
More informationFRESNEL DIFFRACTION AND PARAXIAL WAVE EQUATION. A. Fresnel diffraction
19 IV. FRESNEL DIFFRACTION AND PARAXIAL WAVE EQUATION A. Fresnel diffraction Any physical optical beam is of finite transverse cross section. Beams of finite cross section may be described in terms of
More informationAdvanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude
Advanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude A. Migukin *, V. atkovnik and J. Astola Department of Signal Processing, Tampere University of Technology,
More informationOptics Vac Work MT 2008
Optics Vac Work MT 2008 1. Explain what is meant by the Fraunhofer condition for diffraction. [4] An aperture lies in the plane z = 0 and has amplitude transmission function T(y) independent of x. It is
More informationModeling of diffractive optical elements for lens design.
Juan L. Rayces and Lan Lebich OA Applied Optics, 7421 Orangewood Ave., Garden Grove, A 92641. Abstract. The use of a standard aspheric profile to describe conventional optical elements in lens design programs
More informationRay Optics I. Last time, finished EM theory Looked at complex boundary problems TIR: Snell s law complex Metal mirrors: index complex
Phys 531 Lecture 8 20 September 2005 Ray Optics I Last time, finished EM theory Looked at complex boundary problems TIR: Snell s law complex Metal mirrors: index complex Today shift gears, start applying
More informationMEASUREMENT OF SMALL TRANSVERSE BEAM SIZE USING INTERFEROMETRY
MEASUREMENT OF SMALL TRANSVERSE BEAM SIZE USING INTERFEROMETRY T. Mitsuhashi High Energy Accelerator Research Organisation, Oho, Tsukuba, Ibaraki, 35-81 Japan Abstract The principle of measurement of the
More informationChapter 35 &36 Physical Optics
Chapter 35 &36 Physical Optics Physical Optics Phase Difference & Coherence Thin Film Interference 2-Slit Interference Single Slit Interference Diffraction Patterns Diffraction Grating Diffraction & Resolution
More informationDIFFRACTION 4.1 DIFFRACTION Difference between Interference and Diffraction Classification Of Diffraction Phenomena
4.1 DIFFRACTION Suppose a light wave incident on a slit AB of sufficient width b, as shown in Figure 1. According to concept of rectilinear propagation of light the region A B on the screen should be uniformly
More informationModifications of detour phase computer-generated holograms
Modifications of detour phase computer-generated holograms Uriel Levy, Emanuel Marom, and David Mendlovic The detour phase method for the design of computer-generated holograms can be modified to achieve
More informationDistortion Correction for Conical Multiplex Holography Using Direct Object-Image Relationship
Proc. Natl. Sci. Counc. ROC(A) Vol. 25, No. 5, 2001. pp. 300-308 Distortion Correction for Conical Multiplex Holography Using Direct Object-Image Relationship YIH-SHYANG CHENG, RAY-CHENG CHANG, AND SHIH-YU
More informationReconstruction of purely absorbing, absorbing and phaseshifting, and strong phase-shifting objects from their single-shot in-line holograms
Reconstruction of purely absorbing, absorbing and phaseshifting, and strong phase-shifting objects from their single-shot in-line holograms Tatiana Latychevskaia* and Hans-Werner Fink Physics Institute,
More informationThree Dimensional Measurements by Deflectometry and Double Hilbert Transform
Three Dimensional Measurements by Deflectometry and Double Hilbert Transform Silin Na*, Sanghoon Shin**, Younghun Yu* * Department of Physics, Jeju National University, Jeju, 63243, Korea ** Kanghae Precision
More informationSecuring of Two and Three Dimensional Information Based on In-line Digital Holography
Securing of Two and Three Dimensional Information Based on In-line Digital Holography Hesham Eldeeb Computer & System Department Electronic Research Institute National Research Center, Dokki, Giza Cairo,
More informationSingle Slit Diffraction
Name: Date: PC1142 Physics II Single Slit Diffraction 5 Laboratory Worksheet Part A: Qualitative Observation of Single Slit Diffraction Pattern L = a 2y 0.20 mm 0.02 mm Data Table 1 Question A-1: Describe
More informationFiber Fourier optics
Final version printed as of 4/7/00 Accepted for publication in Optics Letters Fiber Fourier optics A. E. Siegman Ginzton Laboratory, Stanford University, Stanford, California 94305 Received The Fourier
More informationCoupling of surface roughness to the performance of computer-generated holograms
Coupling of surface roughness to the performance of computer-generated holograms Ping Zhou* and Jim Burge College of Optical Sciences, University of Arizona, Tucson, Arizona 85721, USA *Corresponding author:
More informationPhase-Contrast Imaging and Tomography at 60 kev using a Conventional X-ray Tube
Phase-Contrast Imaging and Tomography at 60 kev using a Conventional X-ray Tube T. Donath* a, F. Pfeiffer a,b, O. Bunk a, W. Groot a, M. Bednarzik a, C. Grünzweig a, E. Hempel c, S. Popescu c, M. Hoheisel
More informationEffective Medium Theory, Rough Surfaces, and Moth s Eyes
Effective Medium Theory, Rough Surfaces, and Moth s Eyes R. Steven Turley, David Allred, Anthony Willey, Joseph Muhlestein, and Zephne Larsen Brigham Young University, Provo, Utah Abstract Optics in the
More informationPhysics I : Oscillations and Waves Prof. S Bharadwaj Department of Physics & Meteorology Indian Institute of Technology, Kharagpur
Physics I : Oscillations and Waves Prof. S Bharadwaj Department of Physics & Meteorology Indian Institute of Technology, Kharagpur Lecture - 20 Diffraction - I We have been discussing interference, the
More informationChapter 38 Wave Optics (II)
Chapter 38 Wave Optics (II) Initiation: Young s ideas on light were daring and imaginative, but he did not provide rigorous mathematical theory and, more importantly, he is arrogant. Progress: Fresnel,
More informationTitle. Author(s)Yamaguchi, Kazuhiro; Sakamoto, Yuji. CitationApplied Optics, 48(34): H203-H211. Issue Date Doc URL. Rights.
Title Computer generated hologram with characteristics of Author(s)Yamaguchi, Kazuhiro; Sakamoto, Yuji CitationApplied Optics, 48(34): H203-H211 Issue Date 2009-12-01 Doc URL http://hdl.handle.net/2115/52148
More informationPhase retrieval of diffraction patterns from noncrystalline samples using the oversampling method
Phase retrieval of diffraction patterns from noncrystalline samples using the oversampling method Jianwei Miao* Stanford Synchrotron Radiation Laboratory, Stanford Linear Accelerator Center, Stanford University,
More informationPhysical & Electromagnetic Optics: Basic Principles of Diffraction
24/06/2017 Physical & Electromagnetic Optics: Basic Principles of Diffraction Optical Engineering Prof. Elias N. Glytsis School of Electrical & Computer Engineering National Technical University of Athens
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi:10.1038/nature10934 Supplementary Methods Mathematical implementation of the EST method. The EST method begins with padding each projection with zeros (that is, embedding
More informationA Single Grating-lens Focusing Two Orthogonally Polarized Beams in Opposite Direction
, pp.41-45 http://dx.doi.org/10.14257/astl.2016.140.08 A Single Grating-lens Focusing Two Orthogonally Polarized Beams in Opposite Direction Seung Dae Lee 1 1* Dept. of Electronic Engineering, Namseoul
More informationA new fractal algorithm to model discrete sequences
A new fractal algorithm to model discrete sequences Zhai Ming-Yue( 翟明岳 ) a) Heidi Kuzuma b) and James W. Rector b)c) a) School of EE Engineering North China Electric Power University Beijing 102206 China
More informationMultitemporal Geometric Distortion Correction Utilizing the Affine Transformation
Purdue University Purdue e-pubs LARS Symposia Laboratory for Applications of Remote Sensing 10-1-1973 Multitemporal Geometric Distortion Correction Utilizing the Affine Transformation R. A. Emmert Lawrence
More informationControl of Light. Emmett Ientilucci Digital Imaging and Remote Sensing Laboratory Chester F. Carlson Center for Imaging Science 8 May 2007
Control of Light Emmett Ientilucci Digital Imaging and Remote Sensing Laboratory Chester F. Carlson Center for Imaging Science 8 May 007 Spectro-radiometry Spectral Considerations Chromatic dispersion
More informationUniversity Physics (Prof. David Flory) Chapt_37 Monday, August 06, 2007
Name: Date: 1. If we increase the wavelength of the light used to form a double-slit diffraction pattern: A) the width of the central diffraction peak increases and the number of bright fringes within
More informationImage quality improvement of polygon computer generated holography
Image quality improvement of polygon computer generated holography Xiao-Ning Pang, Ding-Chen Chen, Yi-Cong Ding, Yi-Gui Chen, Shao-Ji Jiang, and Jian-Wen Dong* School of Physics and Engineering, and State
More informationOverview of techniques applicable to self-interference incoherent digital holography
J. Europ. Opt. Soc. Rap. Public. 8, 13077 (2013) www.jeos.org Overview of techniques applicable to self-interference incoherent digital holography J. Hong jisoohong@mail.usf.edu Department of Physics,
More informationOPTI-521 Graduate Report 2 Matthew Risi Tutorial: Introduction to imaging, and estimate of image quality degradation from optical surfaces
OPTI-521 Graduate Report 2 Matthew Risi Tutorial: Introduction to imaging, and estimate of image quality degradation from optical surfaces Abstract The purpose of this tutorial is to introduce the concept
More informationPhysics 202, Lecture 23
Physics 202, Lecture 23 Today s Topics Lights and Laws of Geometric Optics Nature of Light Reflection and Refraction Law of Reflection Law of Refraction Index of Reflection, Snell s Law Total Internal
More informationKeyhole coherent diffractive imaging
Keyhole coherent diffractive imaging LETTERS BRIAN ABBEY 1, KEITH A. NUGENT 1, GARTH J. WILLIAMS 1 *, JESSE N. CLARK 2, ANDREW G. PEELE 2, MARK A. PFEIFER 2, MARTIN DE JONGE 3 AND IAN MCNULTY 3 1 School
More informationDiffraction: Propagation of wave based on Huygens s principle.
Diffraction: In addition to interference, waves also exhibit another property diffraction, which is the bending of waves as they pass by some objects or through an aperture. The phenomenon of diffraction
More informationBasic optics. Geometrical optics and images Interference Diffraction Diffraction integral. we use simple models that say a lot! more rigorous approach
Basic optics Geometrical optics and images Interference Diffraction Diffraction integral we use simple models that say a lot! more rigorous approach Basic optics Geometrical optics and images Interference
More informationMulti-projector-type immersive light field display
Multi-projector-type immersive light field display Qing Zhong ( é) 1, Beishi Chen (í ì) 1, Haifeng Li (Ó ô) 1, Xu Liu ( Ê) 1, Jun Xia ( ) 2, Baoping Wang ( ) 2, and Haisong Xu (Å Ø) 1 1 State Key Laboratory
More informationSubpixel Corner Detection Using Spatial Moment 1)
Vol.31, No.5 ACTA AUTOMATICA SINICA September, 25 Subpixel Corner Detection Using Spatial Moment 1) WANG She-Yang SONG Shen-Min QIANG Wen-Yi CHEN Xing-Lin (Department of Control Engineering, Harbin Institute
More informationMatthew Schwartz Lecture 19: Diffraction and resolution
Matthew Schwartz Lecture 19: Diffraction and resolution 1 Huygens principle Diffraction refers to what happens to a wave when it hits an obstacle. The key to understanding diffraction is a very simple
More informationDiffraction and Interference of Plane Light Waves
PHY 92 Diffraction and Interference of Plane Light Waves Diffraction and Interference of Plane Light Waves Introduction In this experiment you will become familiar with diffraction patterns created when
More informationOptimization of modified volume Fresnel zone plates
Birck Nanotechnology Center Birck and NCN Publications Purdue Libraries Year 2009 Optimization of modified volume Fresnel zone plates Pornsak Srisungsitthisunti Okan Ersoy Xianfan Xu Purdue University,
More information14 Chapter. Interference and Diffraction
14 Chapter Interference and Diffraction 14.1 Superposition of Waves... 14-14.1.1 Interference Conditions for Light Sources... 14-4 14. Young s Double-Slit Experiment... 14-4 Example 14.1: Double-Slit Experiment...
More informationUsing a multipoint interferometer to measure the orbital angular momentum of light
CHAPTER 3 Using a multipoint interferometer to measure the orbital angular momentum of light Recently it was shown that the orbital angular momentum of light can be measured using a multipoint interferometer,
More informationVibration parameter measurement using the temporal digital hologram sequence and windowed Fourier transform
THEORETICAL & APPLIED MECHANICS LETTERS 1, 051008 (2011) Vibration parameter measurement using the temporal digital hologram sequence and windowed Fourier transform Chong Yang, 1, 2 1, a) and Hong Miao
More informationSupplementary Figure 1 Optimum transmissive mask design for shaping an incident light to a desired
Supplementary Figure 1 Optimum transmissive mask design for shaping an incident light to a desired tangential form. (a) The light from the sources and scatterers in the half space (1) passes through the
More informationExperimental reconstruction of a highly reflecting fiber Bragg grating by using spectral regularization and inverse scattering
3284 J. Opt. Soc. Am. A/ Vol. 24, No. 10/ October 2007 Rosenthal et al. Experimental reconstruction of a highly reflecting fiber Bragg grating by using spectral regularization and inverse scattering Amir
More informationWORCESTER POLYTECHNIC INSTITUTE
WORCESTER POLYTECHNIC INSTITUTE MECHANICAL ENGINEERING DEPARTMENT Optical Metrology and NDT ME-593L, C 2018 Introduction: Wave Optics January 2018 Wave optics: coherence Temporal coherence Review interference
More informationDesigning LED array for uniform illumination based on local search algorithm
J. Europ. Opt. Soc. Rap. Public. 9, 14014 (2014) www.jeos.org Designing LED array for uniform illumination based on local search algorithm P. Lei leipldu@163.com School of Information and Electric Engineering,
More informationAnalysis of the effect of target scattering on the operational effectiveness of seeker
2nd International Conference on Materials Science, Machinery and Energy Engineering (MSMEE 2017) Analysis of the effect of target scattering on the operational effectiveness of seeker Jun-chao Yuan1, a,
More informationFresnel and Fourier digital holography architectures: a comparison.
Fresnel and Fourier digital holography architectures: a comparison. Damien P., David S. Monaghan, Nitesh Pandey, Bryan M. Hennelly. Department of Computer Science, National University of Ireland, Maynooth,
More informationChapter 36. Diffraction. Copyright 2014 John Wiley & Sons, Inc. All rights reserved.
Chapter 36 Diffraction Copyright 36-1 Single-Slit Diffraction Learning Objectives 36.01 Describe the diffraction of light waves by a narrow opening and an edge, and also describe the resulting interference
More informationSupplementary materials of Multispectral imaging using a single bucket detector
Supplementary materials of Multispectral imaging using a single bucket detector Liheng Bian 1, Jinli Suo 1,, Guohai Situ 2, Ziwei Li 1, Jingtao Fan 1, Feng Chen 1 and Qionghai Dai 1 1 Department of Automation,
More informationImage Fusion in X-Ray Differential Phase-Contrast Imaging
Image Fusion in X-Ray Differential Phase-Contrast Imaging W. Haas a,b, M. Polyanskaya a, F. Bayer b, K. Gödel b, H. Hofmann a, J. Rieger b, A. Ritter b, T. Weber b, L. Wucherer b, J. Durst b, T. Michel
More informationMetallic Transmission Screen for Sub-wavelength Focusing
Metallic Transmission Screen for Sub-wavelength Focusing A.M.H. Wong, C.D. Sarris and G.V. leftheriades Abstract: A simple metallic transmission screen is proposed that is capable of focusing an incident
More informationA raytracing code for zone plates
A raytracing code for zone plates Alexei Erko *, Franz Schaefers, Nikolay Artemiev a BESSY GmbH, Albert-Einstein-Str.15, 12489 Berlin, Germany a Laboratoire d'optique Appliquee ENSTA Ecole Polytechnique
More informationCalculation of the diffraction efficiency on concave gratings based on Fresnel Kirchhoff s diffraction formula
Calculation of the diffraction efficiency on concave gratings based on Fresnel Kirchhoff s diffraction formula Yuanshen Huang, 1, * Ting Li, 1 Banglian Xu, 1 Ruijin Hong, 1 Chunxian Tao, 1 Jinzhong Ling,
More informationRobust color segmentation algorithms in illumination variation conditions
286 CHINESE OPTICS LETTERS / Vol. 8, No. / March 10, 2010 Robust color segmentation algorithms in illumination variation conditions Jinhui Lan ( ) and Kai Shen ( Department of Measurement and Control Technologies,
More informationIterative aperture mask design in phase space using a rank constraint
Iterative aperture mask design in phase space using a rank constraint Roarke Horstmeyer, 1,* Se Baek Oh, 2 and Ramesh Raskar 1 1 The Media Lab, Massachusetts Institute of Technology, 75 Amherst Street,
More informationFormulas of possible interest
Name: PHYS 3410/6750: Modern Optics Final Exam Thursday 15 December 2011 Prof. Bolton No books, calculators, notes, etc. Formulas of possible interest I = ɛ 0 c E 2 T = 1 2 ɛ 0cE 2 0 E γ = hν γ n = c/v
More informationFundamental Optics for DVD Pickups. The theory of the geometrical aberration and diffraction limits are introduced for
Chapter Fundamental Optics for DVD Pickups.1 Introduction to basic optics The theory of the geometrical aberration and diffraction limits are introduced for estimating the focused laser beam spot of a
More informationTwo pixel computer generated hologram using a zero twist nematic liquid crystal spatial light modulator
Two pixel computer generated hologram using a zero twist nematic liquid crystal spatial light modulator Philip M. Birch, Rupert Young, David Budgett, Chris Chatwin School of Engineering, University of
More informationCMOS compatible highly efficient grating couplers with a stair-step blaze profile
CMOS compatible highly efficient grating couplers with a stair-step blaze profile Zhou Liang( ) a), Li Zhi-Yong( ) a), Hu Ying-Tao( ) a), Xiong Kang( ) a), Fan Zhong-Chao( ) b), Han Wei-Hua( ) b), Yu Yu-De
More informationIterative procedure for in-situ EUV optical testing with an incoherent source
APS/123-QED Iterative procedure for in-situ EUV optical testing with an incoherent source Ryan Miyakawa and Patrick Naulleau Lawrence Berkeley National Laboratory, Berkeley, CA 94720 Avideh Zakhor Dept.
More informationUNIT VI OPTICS ALL THE POSSIBLE FORMULAE
58 UNIT VI OPTICS ALL THE POSSIBLE FORMULAE Relation between focal length and radius of curvature of a mirror/lens, f = R/2 Mirror formula: Magnification produced by a mirror: m = - = - Snell s law: 1
More informationNumerical Propagation Manual
Numerical Propagation Manual Numerical Propagation Manual 1 Numerical Propagation Manual Version 1.2 Juan Pablo Piedrahita Quintero, Raúl Andrés Castañeda Quintero and Jorge Iván García Sucerquia Monday,
More informationRegion segmentation and parallel processing for creating large-scale CGHs in polygon source method
Practical Holography XXIII: Materials and Applications, eds. H. I. Bjelkhagen, R. K. Kostuk, SPIE#733, 7330E(009). 1 Region segmentation and parallel processing for creating large-scale CGHs in polygon
More informationPlane Wave Imaging Using Phased Array Arno Volker 1
11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic More Info at Open Access Database www.ndt.net/?id=16409 Plane Wave Imaging Using Phased Array
More informationNONINTERFEROMETRIC TOMOGRAPHIC RECONSTRUCTION OF 3D STATIC AND DYNAMIC AMPLITUDE AND PHASE OBJECTS. Thesis. Submitted to
NONINTERFEROMETRIC TOMOGRAPHIC RECONSTRUCTION OF 3D STATIC AND DYNAMIC AMPLITUDE AND PHASE OBJECTS Thesis Submitted to The School of Engineering of the UNIVERSITY OF DAYTON In Partial Fulfillment of the
More information*
Resolution enhancement in digital holography by self-extrapolation of holograms Tatiana Latychevskaia * and Hans-Werner Fink Physics Institute, University of Zurich, Winterthurerstrasse 190, CH-8057, Switzerland
More informationSpectrographs. C. A. Griffith, Class Notes, PTYS 521, 2016 Not for distribution.
Spectrographs C A Griffith, Class Notes, PTYS 521, 2016 Not for distribution 1 Spectrographs and their characteristics A spectrograph is an instrument that disperses light into a frequency spectrum, which
More informationAnisotropic scaling of hard disk surface structures
JOURNAL OF APPLIED PHYSICS VOLUME 88, NUMBER 6 15 SEPTEMBER 2000 Anisotropic scaling of hard disk surface structures T. Karabacak a) and Y.-P. Zhao Department of Physics, Applied Physics and Astronomy,
More informationPolarizing properties of embedded symmetric trilayer stacks under conditions of frustrated total internal reflection
University of New Orleans ScholarWorks@UNO Electrical Engineering Faculty Publications Department of Electrical Engineering 3-1-2006 Polarizing properties of embedded symmetric trilayer stacks under conditions
More informationHolographic Elements in Solar Concentrator and Collection Systems
Holographic Elements in Solar Concentrator and Collection Systems Raymond K. Kostuk,2, Jose Castro, Brian Myer 2, Deming Zhang and Glenn Rosenberg 3 Electrical and Computer Engineering, Department University
More informationThree Dimensional Diffraction Imaging
Three Dimensional Diffraction Imaging Leigh Martin Research Group of Dr. John Miao, Physics and Astronomy Department, UCLA 1 Introduction Microscopy has incited many of the greatest revolutions in scientific
More informationManual Diffraction. Manual remote experiment Project e-xperimenteren+ J. Snellenburg, J.M.Mulder
Manual remote experiment Project e-xperimenteren+ J. Snellenburg, J.M.Mulder 30-01-006 Colofon Manual diffraction Manual remote experiment Project e-xperimenteren+ Stichting Digitale Universiteit Oudenoord
More information