Design of paraxial diffractive elements with the CAD system DIGIOPT
|
|
- Ethelbert Barber
- 5 years ago
- Views:
Transcription
1 Design of paraxial diffractive elements with the CAD system DIGIOPT Harald Aagedal, Thomas Beth, Heiko Schwarzer, Stephan Teiwes Institute of Algorithms and Cognitive Systems, University of Karlsruhe Am Fasanengarten 5, D Karlsruhe, Germany Tel.: , Fax: aagedal@ira.uka.de Abstract User-friendly computer-aided design tools for the implementation of diffractive optical elements may on the one hand influence the acceptance of diffractive optics in industry and on the other help to simplify research work in this area. In the present paper we present the software system DIGIOPT which was designed to fulfil our demands on the design of paraxial diffractive elements for optical information processing applications. Key words: DIGIOPT, DE design and simulation, optical information processing. 1. Introduction Diffractive optics is a field of growing importance not only in science but also for industrial applications such as information and material processing. In order to make diffractive optics more attractive to industry easy-to-use design and simulation tools with a high amount of functionality and flexibility have to be developed. An important aspect of diffractive optics is the modelling of inverse wave propagation under consideration of certain boundary conditions such as optical function, material constraints and structure size. On the basis of Fourier optics, methods to design paraxial diffractive elements for applications in information processing, communication technology, laser physics, etc., are developed. In this context, the computer-aided design (CAD) of diffractive elements has several advantages compared to optically produced diffractive elements, simple and flexible modelling of digital signal waves, automated diffractive element design by efficient linear and non-linear optimisation techniques, test by simulation of the optical function of the designed optical system and direct implementation of digital elements, e.g. e-beam lithography or addressable spatial light modulators (SLMs) In principle, commercially available standard software can be used to design diffractive elements [1]. However, the required flexibility and functionality are difficult to achieve with standard software. The design and coding of diffractive elements such as filter, multifacet or transformation elements using standard software, e.g. Mathematica, require a great amount of programming which can only be done by experts knowing the recent design and coding techniques. At present, there are hardly any software packages combining developed theoretical concepts in a structured way under a comfortable user interface, thus making diffractive optics available to experts as well as non-experts of the field. There are many criteria that have to be fulfilled by a CAD system for the design and simulation of diffractive elements. These can be summarised into the following keywords: Functionality. The most important feature of a software system is of course that it is capable of doing what the user wants it to do. In diffractive optics this means that a variety of standard low level operators as well as high level algorithms for the design of diffractive elements should be contained in the system. Physically correct simulation of diffractive elements and wave propagation should also be carried out by the system so that entire optical systems can be simulated correctly.
2 Flexibility. Flexibility of a program for diffractive optics is the ability to combine all kinds of operators in every imaginable way. This is of great importance for the scientific utilisation of the program. Flexibility enables the user to test new ideas and perform rapid prototyping of diffractive elements and entire optical systems. User-friendliness. User-friendliness should be a property of every end-user application. This criterion is however very hard to fulfil and normally only found in application software with a broad market such as word processing systems, spreadsheets, etc. This is unfortunately not the case of systems for the design of diffractive elements. Communication. The software system should be capable of communicating with production devices such as e-beam machines, laser printers, etc.. For scientific documentation it should also be possible to acquire bitmaps, graphs, tables, etc. which can be pasted into standard text processing systems. A program that fulfils these criteria can be used as a design and simulation tool for a scientific and industrial purpose as well as in the education of non-experts of diffractive optics. Thus, such a program could function as a link between research and industry, making new theoretical concepts of diffractive elements available to industry and perhaps even further the acceptance of diffractive optics for industrial utilisation. In the present paper we will describe the software system DIGIOPT for the design and simulation of optical systems based on paraxial diffractive elements which is a good attempt in this respect. 2. The Concept of DIGIOPT The software system DIGIOPT was designed and implemented at the Institute of Algorithms and Cognitive Systems, University of Karlsruhe with the above criteria in mind. It runs under Microsoft Windows on IBM compatible PCs. DIGIOPT is strongly windows oriented. A wavefront or signal is represented by a complex-valued two dimensional array of which either the amplitude, phase, real or imaginary distribution is shown in a window on the DIGIOPT desktop (Figure 1). A variety of elementary preprocessing operations can be performed on each window representing the signal. These include unitary operators for digital preprocessing of 2D signals such as scaling, translation, rotation, embedding, Fourier, Fresnel and Kirchhoff transformation, etc., binary operators such as addition, subtraction, multiplication, division, convolution, correlation, etc. of signals, functions for the generation of primitive wavefronts such as plane, spherical, cylindrical and gaussian waves, etc., as well as parsing of analytical expressions, Figure 1 The windows-oriented concept of DIGIOPT. The amplitude, phase, real and imaginary distribution of a Fresnel transform of the upper left image are shown in the upper row. A complex histogram and a profile of the Fresnel transform are displayed in the lower row. functions for the computation of complex histograms and profiles of signal distributions and functions determining the quality of diffractive elements such as diffraction efficiency and signal-to-noise ratio with respect to the generated signal wave.
3 One of the goals of DIGIOPT was to gather all relevant algorithms for the design of paraxial diffractive elements in one userfriendly program, including recent work concerning the iterative Fourier transform algorithm [2]. The high level features may be summarised into coding operators, a general form of the iterative Fourier transform algorithms, operators for the generation of multifacet elements, 2D display of wave propagation of one dimensional signals and an animation facility (Figure 2). DIGIOPT also contains a simple programming language. This is an extremely helpful tool Figure 2 Animation tool whenever operators are to be concatenated or signals combined in some way or another. The programmer may also use loops as in the following program giving the animation sequence shown in Figure 2. var field_object; var ani_a("test.ani"); var real_d,card_cycles; // Creates an animation window on the DigiOpt desktop field_object:=embedextractfield(loadfield("marilyn.de"),128,128); // Loads and embeds the Marilyn image (64x64) // into a zero matrix of size 128x128 real_d:=15.0; card_cycles:=50; // Distance in mm // Number of iterations repeat card_cycles times real_d:=real_d*1.05; AddToAnimation(ani_a,FRT(field_Object,real_d)); // Adds Fresnel transforms of field_object at distance // real_d to animation ani_a end; The sequence is a collection of Fresnel transforms from distances between 15 mm and 172 mm of the image shown in Window 1 of Figure 5. As depicted in Figure 3, the DIGIOPT programming language may also be used to generate wavefronts defined by analytical expressions. Figure 3 A one-line DIGIOPT program creating the distribution of Window 1. As mentioned in the introduction, a software system should be capable of communicating with production devices such as e- beam machines and laser printers. Through the MS Windows environment DIGIOPT can produce diffractive element distributions in a variety of data formats such as PostScript, GIF and BMP. At present, DIGIOPT does not support GDS2, but this will change in near future. For scientific documentation, bitmaps, graphs, tables, etc. can easily be exported from DIGIOPT and pasted into standard text processing or slide presentation systems. Likewise, online documention of research work as well as real-time simulation presentations with DIGIOPT running on a laptop connected to an overhead display are possible.
4 DIGIOPT was not only designed to simulate an optical bench, but also to be a component of opto-electronical hybrid systems (Figure 4) [3]. In such systems SLMs can be controlled by DIGIOPT. Moreover, the output of the optical system can be captured by a CCD camera and input into DIGIOPT. This makes it possible to perform real-time optical tests of diffractive elements [4]. Input CRT Projection lens DIGIOPT SLM Beam splitter Beam expander Fourierlens Laser 3. Design Example Design of diffractive elements can be a complicated task. In order to give an impression on how a non-trivial design problem can easily be tackled by the DIGIOPT system we design a non-periodic diffractive phase element generating a speckle-free reconstruction of a specified pattern in the Fourier plane. There are several steps which have to be performed in order to yield a diffractive element with satisfactory properties. We will try to focus on how to use DIGIOPT and not on the theoretical background. We have also chosen to use elementary operations which may be performed by higher level operations in one step so that the flexibility of DIGIOPT will become clear. Generating a non-periodic diffractive Fourier element, that is, a diffractive element with a continuous impulse response, is much more difficult than generating a periodic diffractive element with discrete impulse response. This is due to the sinc-interpolation of the digital impulse response which again is caused by the finite size of the diffractive element. For a discrete impulse response this sinc-interpolation does not do any harm, it simply changes the shape of the reconstruction points. In the nonperiodic case, however, speckles in the continuous impulse response will occur if the phase values of the reconstruction points are not carefully chosen. The interpolation problem with unfavourably chosen phase distributions can easily be demonstrated by DIGIOPT as shown in Figure 5. Window 2 was generated by superimposing random Output phase values onto Window 1. To sinc-interpolate Window 2 we simply compute its Fourier transform (Window 3), embeds this into a zero matrix (Window 4) and compute the inverse Fourier transform (Window 5). This is the sinc-interpolation of the signal with the random phase, containing a great amount of speckles. We will now design a non-periodic diffractive phase element with a version of the iterative Fourier transform algorithm (IFTA). Iterative Fourier transform algorithms were first applied to phase retrieval problems in electron microscopy by Gercherg and Saxton [5]. These algorithms are based on an iterative technique which makes the transmission function of the diffractive element fulfil certain constraints in the Fourier plane as well as in the reconstruction plane by allowing some parameters of freedom in the reconstruction plane. In the case of intensity signals, these parameters are freedom of phase and scalation in a signal window in which the desired image is located and freedom of phase and amplitude outside this signal window [2]. CCD Figure 4 DIGIOPT as a controlling component of an optoelectronical hybrid system. Figure 5 Speckles occur after sinc-interpolation if the phase values are not carefully chosen.
5 3.1 Finding an object independent band-limited phase distribution The success of the iterative Fourier transform algorithms concerning the quality of the generated system responses depends on the initial distribution of either the diffractive element or of the signal, depending on whether the iteration starts in the Fourier domain or in the spatial domain. In our case, the iteration begins in the spatial domain. Hence, the initial phase distribution of the object has to be carefully chosen. On the one hand, the starting phase should function as a diffusor, that is, spread the energy into the entire spectrum of the complex object distribution. This can easily be done by superimposing a random phase distribution onto the object. On the other hand, the phase distribution has to be smooth to avoid the speckle effects shown in Figure 5. These initial speckles will not be removed by the iterative Fourier transform. In this sense, smoothness means that the phase values of two neighbouring pixels do not differ by values close to π. This is related to the existence of so-called 2π-jumps [6,7,8]. A random phase distribution is thus not a good choice. In our example we will use an algorithm proposed by Bräuer [7] to find an initial phase distribution. The Bräuer phase distribution was superimposed onto the plane wave in Window 2, giving Window 3, by selection the point Bräuer Start Phase of the dialog box in Figure 7. As shown in Window 4 of Figure 6, this diffusor spreads out the energy into the entire spectrum (Window 4). In order to achieve a smooth object-independent band-limited phase Window 3 has to be interpolated. This may be done by embedding Window 4 into a zero matrix (Window 5), computing the inverse Fourier transform (Window 6) and setting all amplitudes to one. Thus, the phase distribution of Window 8 is band-limited and evenly distributed in the central part of the spectrum. Figure 6 The Modify Phase dialog (left) was used to find an initial phase distribution without 2π-jumps (Window 8). 3.2 Finding an initial object dependent band-limited phase distribution The object independent phase distribution found in Section 3.1 may now be superimposed onto the object. One should, however, be sure to have a band-limited object, that is, the spectrum should not exceed the central one-quarter of the spectral window. This may be achieved by simply sinc-interpolating the object as shown in Windows 1 5 of Figure 7. The pointwise multiplication of the object independent phase distribution and the interpolated object (Window 13) is now almost band-limited (Window 14). A simple iterative Fourier transform algorithm may now be applied to improve the phase distribution of Window 13 even further. This can be done by selecting the point Band Limited Phase (IFTA) of the Modify Phase dialog (Figure 6).
6 Figure 7 Before superimposing the object independent band-limited phase distribution onto the object, the object should be sinc-interpolated (Windows 1 5) to ensure band-limitation. The pointwise multiplication (Window 13) of Window 5 and 12 is now almost band-limited, as shown in Window 14. By applying a simple IFTA-version the phase distribution of Window 13 may be further improved (Window 15 and 16). 3.3 Generation of the diffractive element After having performed the steps in Sections 3.1 and 3.2 we have superimposed a phase distribution onto the image which the diffractive element is to generate as an intensity distribution in the Fourier plane. As already mentioned in Section 0, the iterative Fourier transform needs parameters of freedom in order to achieve diffractive elements with satisfactory properties. For intensity signals these are freedom of phase and scalation in a window containing the desired signal and freedom of phase and amplitude outside the window. The signal window may be selected by simply using the mouse or defining the coordinates (Window 21 in Figure 8). The size of the signal window is crucial for the quality of the diffractive element. If there is not enough space outside the signal window, the freedom of amplitude and phase will normally be too small to obtain an element with a satisfactory signal-to-noise-ratio [2]. Hence, the signal (Figure 8, Window 17) was embedded into a zero matrix (Window 21) and the signal window of Window 21 was selected to hold the signal with a margin of 4 pixels. The diffractive element designed was to be illuminated by a plane wave. Because of the speckle reduction we embedded a plane wave, i.e. a field I(x) with I(x)=1, into a zero matrix of the double size in order to simulate the finite size of the diffractive element. We could of course have chosen any other field (also a complex-valued one) as the illumination field of the diffractive element.
7 To start the iteration of the diffractive element the signal field (Window 21) and the illumination wave (Window 22) have to be specified in the Diffractive Element Iteration dialog box. If no illumination wave is specified, a plane wave (unembedded) is assumed. The type of the diffractive element may either be set to Amplitude, Phase or Complex, whereas Figure 8 The dialog box for a general iterative Fourier transform algorithm. Window 21 selected as signal field, whereas Window 22 defines the illumination wave. the signal type is either Intensity or Complex. In our case we choose a diffractive phase element and an intensity signal. The Animation group makes it possible to save the result of every iteration in two animation sequences one for the diffractive elements and one for their impulse responses. This lets the user examine the development of the diffractive element during the iteration process. The Freedom iterations group (Figure 8) should be described more carefully. In this group the number of iteration cycles has to be specified. The iteration process is divided into three stages, Phase freedom iteration, Amplitude-Phase-Scale freedom iteration and Quantization. During the phase freedom iteration only the phase values of the pixels in the signal window may develop freely. The noise outside the signal window is suppressed. This stage is normally used to increase the diffraction efficiency of the diffractive element. It may, however, destroy the band limitation of the phase distribution superimposed onto the intensity signal. Thus, in the case of speckle reduction no Phase freedom iteration should be performed. If a complex signal type is selected, this stage is disabled. The only constraint of this iterative Fourier transform algorithm in the space domain during the Amplitude-Phase-Scale freedom iteration is the replacement of the intensity in the signal window, that is, the freedom parameters are amplitude and phase outside the signal window and phase freedom in the signal window (only in the case of intensity signals). The constraints in the Fourier domain depend on whether Soft or Hard coding is selected in the Coding group. Hard coding forces the transmission function of every iteration to lie in the set of legal transmission functions defined by the DOE Type group. If soft coding was selected, the transmission function is slowly forced towards this set during the
8 Diff.Eff. SNR Figure 10 During the iteration the diffraction efficiency and SNR value of the current diffractive element distribution are appended to a log file displayed on the DIGIOPT desktop (left). These values may be pasted directly into standard windows software and displayed as graphs, tables, etc. (right). iteration stage. This leads only to small changes in the phase relation of neighbouring pixels after each iteration, thus speckles may be avoided. In our example, we performed 200 iterations with amplitude, phase and scale freedom. The last stage, Quantization, enables the quantization of the diffractive element. If a non-zero value is specified, a subdialog is displayed asking for the number of quantization levels and coding type. The coding type is either hard or soft as in the Amplitude-Phase-Scale freedom iteration. In our example, we decided to perform 200 iterations of soft 4-level quantization. The Log File group lets the user create a log file displayed on the DIGIOPT desktop during the iteration process. The log file enables the examination of the development of the transmission function of the diffractive element. Figure 10 shows a log file window and how this can be pasted into standard software and displayed as a graph. Figure 9 shows the results of the applied iterative Fourier transform algorithm. The upper left image is the continuous diffractive phase element that was created after the Amplitude-Phase-Scale freedom iteration. Its sinc-interpolated reconstruction is shown to the right. The lower row shows the 4-level quantized diffractive element obtained after the Quantization stage and its sinc-interpolated reconstruction. The reconstruction of the continuous element has a good signal-to-noise ratio and is practically speckle-free. Unfortunately, the computation of quantized non-periodic diffractive elements with few quantization levels is a very difficult task. Although the reconstruction of the 4-level quantized diffractive element does not contain any speckles, the signal-to-noise ratio is poor. This is easily seen in Figure Conclusion We have stated criteria for the goodness of software systems for the design and simulation of paraxial diffractive elements. The DIGIOPT system was described and demonstrated by generating a non-periodic diffractive phase element giving a continuous speckle-free reconstruction of a desired pattern which is a complicated task. The computation of the diffractive element was divided into three stages; (i) finding an object independent band-limited phase distribution, (ii) finding an initial object dependent band-limited phase distribution and (iii) generation of the diffractive element. Although not necessary, many low level operators were used instead of fewer high level ones to show the flexibility and functionality of DIGIOPT. The example demonstrated the ease of use of DIGIOPT even in solving non-trivial design problems. Figure 9 The upper row shows the computed continuous DPE with its sinc-interpolated reconstruction. The 4-level quantized DPE is shown in the lower row.
9 The programming language of DIGIOPT was demonstrated by two small examples, one producing an animation sequence containing the diffraction pattern of an image at different distances behind the image and one creating a wavefront window from an analytical expression. The programming language offers many interesting features for the advanced user. Considering the non-trivial design example, the variety of features of DIGIOPT and the ease of use, DIGIOPT seems to be a system fulfilling the criteria stated for the goodness of software systems for the design and simulation of diffractive optical systems. Thus, DIGIOPT could possibly function as a link between research and industry in the field of optical signal processing based on paraxial diffractive elements. 5. References [1] D. C. O Shea, J. W. Beltic, and Menelaous Poutous, Binary-mask generation for diffractive optical elements using microcomputers, Appl. Optics, Vol. 32, No. 14, May [2] F. Wyrowski and O. Bryngdahl, Digital holography as part of diffractive optics, Rep. Prog. in Physics 54, , [3] H. Aagedal, T. Beth, F. Reichel, H. Schwarzer and S. Teiwes, Opto-electronical Hybrid System for Rapid Prototyping of Paraxial Diffractive Elements, Presented at Spring Topical Meeting, Spatial Light Modulators and Applications, Opt.Soc.Am., Salt Lake City, [4] H. Aagedal, T. Beth, H. Schwarzer and S. Teiwes, Modern concepts for computer-aided design in diffractive optics, OSA Proceedings of the International Optical Design Conference, G.W. Forbes, ed., Vol. 22, , [5] R.W. Gerchberg and W.O. Saxton, A practical algorithm for the determination of phase from image and diffraction plane pictures, Optik 35, , [6] M.S. Scivier and M.A. Fiddy, Phase ambiguities and the zeros of multidimensional band-limited functions, J.Opt.Soc.Am. A 2, , [7] R. Bräuer, F. Wyrowski and O. Bryngdahl, Diffusors in digital holography, J.Opt.Soc.Am. A 8, , [8] S. Teiwes, H. Schwarzer, H. Aagedal, M. Schmid. T. Beth, F. Reichel and F. Wyrowski, Speckle-free Signal Reconstruction with Non-periodic Diffractive Elements on Liquid-crystal SLM, Presented at Spring Topical Meeting, Spatial Light Modulators and Applications, Opt.Soc.Am., Salt Lake City, 1995.
A Comparison of the Iterative Fourier Transform Method and. Evolutionary Algorithms for the Design of Diffractive Optical.
A Comparison of the Iterative Fourier Transform Method and Evolutionary Algorithms for the Design of Diffractive Optical Elements Philip Birch, Rupert Young, Maria Farsari, David Budgett, John Richardson,
More informationHOLOEYE Photonics. HOLOEYE Photonics AG. HOLOEYE Corporation
HOLOEYE Photonics Products and services in the field of diffractive micro-optics Spatial Light Modulator (SLM) for the industrial research R&D in the field of diffractive optics Micro-display technologies
More informationLAGRANGIAN PARTICLE TRACKING IN ISOTROPIC TURBULENT FLOW VIA HOLOGRAPHIC AND INTENSITY BASED STEREOSCOPY. By Kamran Arjomand
LAGRANGIAN PARTICLE TRACKING IN ISOTROPIC TURBULENT FLOW VIA HOLOGRAPHIC AND INTENSITY BASED STEREOSCOPY By Kamran Arjomand I. Background A. Holographic Imaging 1. Acquire Hologram 3. Numerical Reconstruction
More informationSIMULATION AND VISUALIZATION IN THE EDUCATION OF COHERENT OPTICS
SIMULATION AND VISUALIZATION IN THE EDUCATION OF COHERENT OPTICS J. KORNIS, P. PACHER Department of Physics Technical University of Budapest H-1111 Budafoki út 8., Hungary e-mail: kornis@phy.bme.hu, pacher@phy.bme.hu
More informationLecture Notes on Wave Optics (04/23/14) 2.71/2.710 Introduction to Optics Nick Fang
.7/.70 Introduction to Optics Nic Fang Outline: Fresnel Diffraction The Depth of Focus and Depth of Field(DOF) Fresnel Zones and Zone Plates Holography A. Fresnel Diffraction For the general diffraction
More informationDiffractive optics: learning by computer experiments
Invited Paper Diffractive optics: learning by computer experiments Hartmut Bartelt*, Institute for Physical High Technology, Jena I Germany ABSTRACT Simulation and modeling software is especially useful
More informationZero Order Correction of Shift-multiplexed Computer Generated Fourier Holograms Recorded in Incoherent Projection Scheme
VII International Conference on Photonics and Information Optics Volume 2018 Conference Paper Zero Order Correction of Shift-multiplexed Computer Generated Fourier Holograms Recorded in Incoherent Projection
More informationFourier, 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 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 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 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 informationComputed tomography (Item No.: P )
Computed tomography (Item No.: P2550100) Curricular Relevance Area of Expertise: Biology Education Level: University Topic: Modern Imaging Methods Subtopic: X-ray Imaging Experiment: Computed tomography
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 informationarxiv: v1 [physics.optics] 1 Mar 2015
Random phase-free kinoform for large objects arxiv:1503.00365v1 [physics.optics] 1 Mar 2015 Tomoyoshi Shimobaba, 1 Takashi Kakue, 1 Yutaka Endo, 1 Ryuji Hirayama, 1 Daisuke Hiyama, 1 Satoki Hasegawa, 1
More informationLASer Cavity Analysis and Design
The unique combination of simulation tools for LASer Cavity Analysis and Design During the last 15 years LASCAD has become industry-leading so ware for LASer Cavity Analysis and Design. The feedback from
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 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 informationInvited Paper. Nukui-Kitamachi, Koganei, Tokyo, , Japan ABSTRACT 1. INTRODUCTION
Invited Paper Wavefront printing technique with overlapping approach toward high definition holographic image reconstruction K. Wakunami* a, R. Oi a, T. Senoh a, H. Sasaki a, Y. Ichihashi a, K. Yamamoto
More informationRedundant Data Elimination for Image Compression and Internet Transmission using MATLAB
Redundant Data Elimination for Image Compression and Internet Transmission using MATLAB R. Challoo, I.P. Thota, and L. Challoo Texas A&M University-Kingsville Kingsville, Texas 78363-8202, U.S.A. ABSTRACT
More informationAn Intuitive Explanation of Fourier Theory
An Intuitive Explanation of Fourier Theory Steven Lehar slehar@cns.bu.edu Fourier theory is pretty complicated mathematically. But there are some beautifully simple holistic concepts behind Fourier theory
More informationFundamentals of Digital Image Processing
\L\.6 Gw.i Fundamentals of Digital Image Processing A Practical Approach with Examples in Matlab Chris Solomon School of Physical Sciences, University of Kent, Canterbury, UK Toby Breckon School of Engineering,
More information3. Image formation, Fourier analysis and CTF theory. Paula da Fonseca
3. Image formation, Fourier analysis and CTF theory Paula da Fonseca EM course 2017 - Agenda - Overview of: Introduction to Fourier analysis o o o o Sine waves Fourier transform (simple examples of 1D
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 informationUNIT 102-9: INTERFERENCE AND DIFFRACTION
Name St.No. - Date(YY/MM/DD) / / Section Group # UNIT 102-9: INTERFERENCE AND DIFFRACTION Patterns created by interference of light in a thin film. OBJECTIVES 1. Understand the creation of double-slit
More information(Refer Slide Time 00:17) Welcome to the course on Digital Image Processing. (Refer Slide Time 00:22)
Digital Image Processing Prof. P. K. Biswas Department of Electronics and Electrical Communications Engineering Indian Institute of Technology, Kharagpur Module Number 01 Lecture Number 02 Application
More informationFRED Slit Diffraction Application Note
FRED Slit Diffraction Application Note The classic problem of diffraction through a slit finds one of its chief applications in spectrometers. The wave nature of these phenomena can be modeled quite accurately
More informationDietrich Paulus Joachim Hornegger. Pattern Recognition of Images and Speech in C++
Dietrich Paulus Joachim Hornegger Pattern Recognition of Images and Speech in C++ To Dorothea, Belinda, and Dominik In the text we use the following names which are protected, trademarks owned by a company
More informationPROCEEDINGS OF SPIE. Simulation of lateral color for a hybrid refractive-diffractive eyepiece by field tracing methods
PROCEEDINGS OF SPIE SPIEDigitalLibrary.org/conference-proceedings-of-spie Simulation of lateral color for a hybrid refractive-diffractive eyepiece by field tracing methods D. Batte, M. Kuhn, F. Wyrowski
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 informationVivekananda. Collegee of Engineering & Technology. Question and Answers on 10CS762 /10IS762 UNIT- 5 : IMAGE ENHANCEMENT.
Vivekananda Collegee of Engineering & Technology Question and Answers on 10CS762 /10IS762 UNIT- 5 : IMAGE ENHANCEMENT Dept. Prepared by Harivinod N Assistant Professor, of Computer Science and Engineering,
More informationChapter - 2 : IMAGE ENHANCEMENT
Chapter - : IMAGE ENHANCEMENT The principal objective of enhancement technique is to process a given image so that the result is more suitable than the original image for a specific application Image Enhancement
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 informationWhat is Frequency Domain Analysis?
R&D Technical Bulletin P. de Groot 9/3/93 What is Frequency Domain Analysis? Abstract: The Zygo NewView is a scanning white-light interferometer that uses frequency domain analysis (FDA) to generate quantitative
More informationDigital Image Processing
Digital Image Processing Third Edition Rafael C. Gonzalez University of Tennessee Richard E. Woods MedData Interactive PEARSON Prentice Hall Pearson Education International Contents Preface xv Acknowledgments
More informationHANDBOOK OF THE MOIRE FRINGE TECHNIQUE
k HANDBOOK OF THE MOIRE FRINGE TECHNIQUE K. PATORSKI Institute for Design of Precise and Optical Instruments Warsaw University of Technology Warsaw, Poland with a contribution by M. KUJAWINSKA Institute
More informationHyperspectral interferometry for single-shot absolute measurement of 3-D shape and displacement fields
EPJ Web of Conferences 6, 6 10007 (2010) DOI:10.1051/epjconf/20100610007 Owned by the authors, published by EDP Sciences, 2010 Hyperspectral interferometry for single-shot absolute measurement of 3-D shape
More informationWinter College on Optics in Environmental Science February Adaptive Optics: Introduction, and Wavefront Correction
2018-23 Winter College on Optics in Environmental Science 2-18 February 2009 Adaptive Optics: Introduction, and Wavefront Correction Love G. University of Durham U.K. Adaptive Optics: Gordon D. Love Durham
More informationShading of a computer-generated hologram by zone plate modulation
Shading of a computer-generated hologram by zone plate modulation Takayuki Kurihara * and Yasuhiro Takaki Institute of Engineering, Tokyo University of Agriculture and Technology, 2-24-16 Naka-cho, Koganei,Tokyo
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 informationUniversity of Huddersfield Repository
University of Huddersfield Repository Muhamedsalih, Hussam, Jiang, Xiang and Gao, F. Interferograms analysis for wavelength scanning interferometer using convolution and fourier transform. Original Citation
More informationComputed tomography of simple objects. Related topics. Principle. Equipment TEP Beam hardening, artefacts, and algorithms
Related topics Beam hardening, artefacts, and algorithms Principle The CT principle is demonstrated with the aid of simple objects. In the case of very simple targets, only a few images need to be taken
More informationannual report 2011 / 2012 INSTITUT FÜR TECHNISCHE OPTIK UNIVERSITÄT STUTTGART
annual report 2011 / 2012 INSTITUT FÜR TECHNISCHE OPTIK UNIVERSITÄT STUTTGART INSTITUT FÜR TECHNISCHE OPTIK UNIVERSITÄT STUTTGART Prof. Dr. W. Osten Pfaffenwaldring 9 D-70569 Stuttgart Tel.: +49(0)711
More informationLens Design I. Lecture 11: Imaging Herbert Gross. Summer term
Lens Design I Lecture 11: Imaging 2015-06-29 Herbert Gross Summer term 2015 www.iap.uni-jena.de 2 Preliminary Schedule 1 13.04. Basics 2 20.04. Properties of optical systrems I 3 27.05. 4 04.05. Properties
More informationspecular diffuse reflection.
Lesson 8 Light and Optics The Nature of Light Properties of Light: Reflection Refraction Interference Diffraction Polarization Dispersion and Prisms Total Internal Reflection Huygens s Principle The Nature
More informationarxiv: v2 [physics.comp-ph] 1 Aug 2018
arxiv:1807.11901v2 [physics.comp-ph] 1 Aug 2018 GDoeSII: Open source software for design of diffractive optical elements and conversion to GDSII lithography files Raghu Dharmavarapu a,b,, Shanti Bhattacharya
More informationRange Sensors (time of flight) (1)
Range Sensors (time of flight) (1) Large range distance measurement -> called range sensors Range information: key element for localization and environment modeling Ultrasonic sensors, infra-red sensors
More informationNormal vector method for the RCWA with automated vector field generation
P. Götz, T. Schuster, K. Frenner, S. Rafler, W. Osten Normal vector method for the RCWA with automated vector field generation Stuttgart, May 2011 Institute of Applied Optics, University of Stuttgart,
More informationUnit - I Computer vision Fundamentals
Unit - I Computer vision Fundamentals It is an area which concentrates on mimicking human vision systems. As a scientific discipline, computer vision is concerned with the theory behind artificial systems
More informationDownloaded From: on 05/11/2015 Terms of Use:
Texture-Based hologram generation using triangles Matthias Koenig a,oliver Deussen b, Thomas Strothotte a a University of Magdeburg, FIN-ISG, 39016 Magdeburg, Germany b Dresden University of Technology,
More informationTechniques of Noninvasive Optical Tomographic Imaging
Techniques of Noninvasive Optical Tomographic Imaging Joseph Rosen*, David Abookasis and Mark Gokhler Ben-Gurion University of the Negev Department of Electrical and Computer Engineering P. O. Box 653,
More informationSLM Developer System Liquid Crystal Spatial Light Modulator and dual head PC system for smart optics developers
SLM Developer System Liquid Crystal Spatial Light Modulator and dual head PC system for smart optics developers Modulator Features High spatial resolution - SXGA (1280 x 1024 pixels) Very high fill factor
More informationEECS 556 Image Processing W 09. Image enhancement. Smoothing and noise removal Sharpening filters
EECS 556 Image Processing W 09 Image enhancement Smoothing and noise removal Sharpening filters What is image processing? Image processing is the application of 2D signal processing methods to images Image
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 informationChapter 3 Geometric Optics
Chapter 3 Geometric Optics [Reading assignment: Goodman, Fourier Optics, Appendix B Ray Optics The full three dimensional wave equation is: (3.) One solution is E E o ûe i ωt± k r ( ). This is a plane
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 informationMetrology and Sensing
Metrology and Sensing Lecture 4: Fringe projection 2016-11-08 Herbert Gross Winter term 2016 www.iap.uni-jena.de 2 Preliminary Schedule No Date Subject Detailed Content 1 18.10. Introduction Introduction,
More informationMetrology and Sensing
Metrology and Sensing Lecture 4: Fringe projection 2017-11-09 Herbert Gross Winter term 2017 www.iap.uni-jena.de 2 Preliminary Schedule No Date Subject Detailed Content 1 19.10. Introduction Introduction,
More informationComputer Vision 2. SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung. Computer Vision 2 Dr. Benjamin Guthier
Computer Vision 2 SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung Computer Vision 2 Dr. Benjamin Guthier 1. IMAGE PROCESSING Computer Vision 2 Dr. Benjamin Guthier Content of this Chapter Non-linear
More informationSUPPLEMENTARY INFORMATION DOI: /NPHOTON
DOI:.38/NPHOTON.2.85 This supplement has two parts. In part A, we provide the rigorous details of the wavefront correction algorithm and show numerical simulations and experimental data for the cases of
More informationELEC Dr Reji Mathew Electrical Engineering UNSW
ELEC 4622 Dr Reji Mathew Electrical Engineering UNSW Review of Motion Modelling and Estimation Introduction to Motion Modelling & Estimation Forward Motion Backward Motion Block Motion Estimation Motion
More informationEnhancing the pictorial content of digital holograms at 100 frames per second
Enhancing the pictorial content of digital holograms at 100 frames per second P.W.M. Tsang, 1 T.-C Poon, 2 and K.W.K. Cheung 1 1 Department of Electronic Engineering, City University of Hong Kong, Hong
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 informationAliasing. Can t draw smooth lines on discrete raster device get staircased lines ( jaggies ):
(Anti)Aliasing and Image Manipulation for (y = 0; y < Size; y++) { for (x = 0; x < Size; x++) { Image[x][y] = 7 + 8 * sin((sqr(x Size) + SQR(y Size)) / 3.0); } } // Size = Size / ; Aliasing Can t draw
More informationExploiting scattering media for exploring 3D objects
Exploiting scattering media for exploring 3D objects Alok Kumar Singh 1, Dinesh N Naik 1,, Giancarlo Pedrini 1, Mitsuo Takeda 1, 3 and Wolfgang Osten 1 1 Institut für Technische Optik and Stuttgart Research
More informationMetrology and Sensing
Metrology and Sensing Lecture 4: Fringe projection 2018-11-09 Herbert Gross Winter term 2018 www.iap.uni-jena.de 2 Schedule Optical Metrology and Sensing 2018 No Date Subject Detailed Content 1 16.10.
More informationIntroduction to Computer-Based Holography
Mitglied der Helmholtz-Gemeinschaft Introduction to Computer-Based Holography October 29, 2013 Carsten Karbach, Jülich Supercomputing Centre (JSC) Why computergenerated holography? Applications Source:
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 informationHolographic Method for Extracting Three-Dimensional Information with a CCD Camera. Synopsis
Mem. Fac. Eng., Osaka City Univ., Vol. 36,pp. 1-11.(1995) Holographic Method for Extracting Three-Dimensional Information with a CCD Camera by Hideki OKAMOTO*, Hiroaki DEDA*, Hideya TAKAHASHI**, and Eiji
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 informationBabu Madhav Institute of Information Technology Years Integrated M.Sc.(IT)(Semester - 7)
5 Years Integrated M.Sc.(IT)(Semester - 7) 060010707 Digital Image Processing UNIT 1 Introduction to Image Processing Q: 1 Answer in short. 1. What is digital image? 1. Define pixel or picture element?
More information15c Fourier Optics Lab (Exercises to be Submitted are listed at the End) 1. Review of geometric optics NOT 2. Wave Optics Introduction
15c Fourier Optics Lab (Exercises to be Submitted are listed at the End) 1. Review of geometric optics a. Missions: i. Image the ceiling fluorescent lights onto the floor using a lens ii. You have a flashlight
More informationAUTOMATED 4 AXIS ADAYfIVE SCANNING WITH THE DIGIBOTICS LASER DIGITIZER
AUTOMATED 4 AXIS ADAYfIVE SCANNING WITH THE DIGIBOTICS LASER DIGITIZER INTRODUCTION The DIGIBOT 3D Laser Digitizer is a high performance 3D input device which combines laser ranging technology, personal
More informationIMAGE DE-NOISING IN WAVELET DOMAIN
IMAGE DE-NOISING IN WAVELET DOMAIN Aaditya Verma a, Shrey Agarwal a a Department of Civil Engineering, Indian Institute of Technology, Kanpur, India - (aaditya, ashrey)@iitk.ac.in KEY WORDS: Wavelets,
More informationDigital Image Fundamentals
Digital Image Fundamentals Image Quality Objective/ subjective Machine/human beings Mathematical and Probabilistic/ human intuition and perception 6 Structure of the Human Eye photoreceptor cells 75~50
More informationCompression of RADARSAT Data with Block Adaptive Wavelets Abstract: 1. Introduction
Compression of RADARSAT Data with Block Adaptive Wavelets Ian Cumming and Jing Wang Department of Electrical and Computer Engineering The University of British Columbia 2356 Main Mall, Vancouver, BC, Canada
More informationDigital Image Processing. Image Enhancement in the Frequency Domain
Digital Image Processing Image Enhancement in the Frequency Domain Topics Frequency Domain Enhancements Fourier Transform Convolution High Pass Filtering in Frequency Domain Low Pass Filtering in Frequency
More informationEffects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines
Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Velocity models used for wavefield-based seismic
More informationCHAPTER 2: THREE DIMENSIONAL TOPOGRAPHICAL MAPPING SYSTEM. Target Object
CHAPTER 2: THREE DIMENSIONAL TOPOGRAPHICAL MAPPING SYSTEM 2.1 Theory and Construction Target Object Laser Projector CCD Camera Host Computer / Image Processor Figure 2.1 Block Diagram of 3D Areal Mapper
More informationLecture 8 Object Descriptors
Lecture 8 Object Descriptors Azadeh Fakhrzadeh Centre for Image Analysis Swedish University of Agricultural Sciences Uppsala University 2 Reading instructions Chapter 11.1 11.4 in G-W Azadeh Fakhrzadeh
More informationTRUE 3D DIGITAL HOLOGRAPHIC TOMOGRAPHY FOR VIRTUAL REALITY APPLICATIONS. Thesis. Submitted to. The School of Engineering of the UNIVERSITY OF DAYTON
TRUE 3D DIGITAL HOLOGRAPHIC TOMOGRAPHY FOR VIRTUAL REALITY APPLICATIONS Thesis Submitted to The School of Engineering of the UNIVERSITY OF DAYTON In Partial Fulfillment of the Requirements for The Degree
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 informationAdvanced Image Reconstruction Methods for Photoacoustic Tomography
Advanced Image Reconstruction Methods for Photoacoustic Tomography Mark A. Anastasio, Kun Wang, and Robert Schoonover Department of Biomedical Engineering Washington University in St. Louis 1 Outline Photoacoustic/thermoacoustic
More informationImage Processing Lecture 10
Image Restoration Image restoration attempts to reconstruct or recover an image that has been degraded by a degradation phenomenon. Thus, restoration techniques are oriented toward modeling the degradation
More informationX-ray Crystallography
X-ray Crystallography Rhodes, Chapters 2, 5 Chapters 3 and 4 include some general considerations MacPherson (on reserve in library), Chapters 1, 3, and 4 What does anything look like? We see objects by
More informationImage Processing. Filtering. Slide 1
Image Processing Filtering Slide 1 Preliminary Image generation Original Noise Image restoration Result Slide 2 Preliminary Classic application: denoising However: Denoising is much more than a simple
More informationContext based educational Java Applets using consumer products
Center for Active Learning of Microelectronics and Photonics Context based educational Java Applets using consumer products A. N. Cartwright, Chu R. Wie, Pratibha Gopalam, N. Liu, Z. Yuan,T. Tang State
More informationComputer Graphics. Sampling Theory & Anti-Aliasing. Philipp Slusallek
Computer Graphics Sampling Theory & Anti-Aliasing Philipp Slusallek Dirac Comb (1) Constant & δ-function flash Comb/Shah function 2 Dirac Comb (2) Constant & δ-function Duality f(x) = K F(ω) = K (ω) And
More informationInvestigation on reconstruction methods applied to 3D terahertz computed Tomography
Investigation on reconstruction methods applied to 3D terahertz computed Tomography B. Recur, 3 A. Younus, 1, P. Mounaix 1, S. Salort, 2 B. Chassagne, 2 P. Desbarats, 3 J-P. Caumes, 2 and E. Abraham 1
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 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 informationReduction of reconstructed particle elongation using iterative min-max filtering in holographic particle image velocimetry
Reduction of reconstructed particle elongation using iterative min-max filtering in holographic particle image velocimetry Yohsuke Tanaka 1, *, Shigeru Murata 1 1: Department of Mechanical System Engineering,
More informationCoE4TN3 Medical Image Processing
CoE4TN3 Medical Image Processing Image Restoration Noise Image sensor might produce noise because of environmental conditions or quality of sensing elements. Interference in the image transmission channel.
More informationTHREE-DIMENSIONA L ELECTRON MICROSCOP Y OF MACROMOLECULAR ASSEMBLIE S. Visualization of Biological Molecules in Their Native Stat e.
THREE-DIMENSIONA L ELECTRON MICROSCOP Y OF MACROMOLECULAR ASSEMBLIE S Visualization of Biological Molecules in Their Native Stat e Joachim Frank CHAPTER 1 Introduction 1 1 The Electron Microscope and
More informationMODERN DESCRIPTIVE GEOMETRY SUPPORTED BY 3D COMPUTER MODELLING
International Conference on Mathematics Textbook Research and Development 2014 () 29-31 July 2014, University of Southampton, UK MODERN DESCRIPTIVE GEOMETRY SUPPORTED BY 3D COMPUTER MODELLING Petra Surynková
More informationMulti-frame blind deconvolution: Compact and multi-channel versions. Douglas A. Hope and Stuart M. Jefferies
Multi-frame blind deconvolution: Compact and multi-channel versions Douglas A. Hope and Stuart M. Jefferies Institute for Astronomy, University of Hawaii, 34 Ohia Ku Street, Pualani, HI 96768, USA ABSTRACT
More informationEvaluation of Two-Dimensional Phase Unwrapping Algorithms for Interferometric Adaptive Optics Utilizing Liquid-Crystal Spatial Light Modulators
48 The Open Optics Journal, 008,, 48-5 Open Access Evaluation of Two-Dimensional Phase Unwrapping Algorithms for Interferometric Adaptive Optics Utilizing Liquid-Crystal Spatial Light Modulators K.L. Baker
More informationOptical Design with Zemax for PhD
Optical Design with Zemax for PhD Lecture : Physical Optics 06-03-3 Herbert Gross Winter term 05 www.iap.uni-jena.de Preliminary Schedule No Date Subject Detailed content.. Introduction 0.. Basic Zemax
More informationALGORITHM OF PRELIMINARY SESSILE DROP PROCESSING FOR NOISE SUPPRESSION X
ALGORITHM OF PRELIMINARY SESSILE DROP PROCESSING FOR NOISE SUPPRESSION X Y.P. KULYNYCH, L.I. MURAVSKY, R.S. BATCHEVSKY, O.G. KUTS Physic-Mechanical Institute of the National Academy of Sciences of Ukraine,
More informationMichael Moody School of Pharmacy University of London 29/39 Brunswick Square London WC1N 1AX, U.K.
This material is provided for educational use only. The information in these slides including all data, images and related materials are the property of : Michael Moody School of Pharmacy University of
More information