Projection Moiré Profilometry using Liquid Crystal Digital Gratings

Similar documents
Optimization of frequency quantization. VN Tibabishev. Keywords: optimization, sampling frequency, the substitution frequencies.

Introduction to Three-phase Circuits. Balanced 3-phase systems Unbalanced 3-phase systems

Methods for Determination of Mean Speckle Size in Simulated Speckle Pattern

Figure 1a. A planar mechanism.

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

SPH3U1 Lesson 06 Kinematics

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

Interference is when two (or more) sets of waves meet and combine to produce a new pattern.

"NEET / AIIMS " SOLUTION (6) Avail Video Lectures of Experienced Faculty.

Solution to HW14 Fall-2002

sin sin Reminder, repetition Image formation by simple curved surface (sphere with radius r): The power (refractive strength):

Physics 2010 Motion with Constant Acceleration Experiment 1

Unit code: H/ QCF level: 5 Credit value: 15 OUTCOME 3 - STATIC AND DYNAMIC FLUID SYSTEMS TUTORIAL 3 - VISCOSITY

Cambridge Assessment International Education Cambridge Ordinary Level. Published

A Primer on Dispersion in Waveguides

Complete Path Planning for a Planar 2-R Manipulator With Point Obstacles

Numerical Simulation of the Thermal Resposne Test Within the Comsol Multiphysics Environment

Chapter 3: Cluster Analysis

Hubble s Law PHYS 1301

ENSC Discrete Time Systems. Project Outline. Semester

Bootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) >

General Chemistry II, Unit I: Study Guide (part I)

CONSTRUCTING STATECHART DIAGRAMS

Sections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic.

ENGINEERING COUNCIL CERTIFICATE LEVEL THERMODYNAMIC, FLUID AND PROCESS ENGINEERING C106 TUTORIAL 5 THE VISCOUS NATURE OF FLUIDS

making triangle (ie same reference angle) ). This is a standard form that will allow us all to have the X= y=

Inertial Mass of Charged Elementary Particles

1. Calculating and displaying the electric field of a single charged particle. r p. q p. r o

Synchronous Motor V-Curves

Modulational instabilities and Fermi Pasta Ulam recurrence. in a coupled long wave short wave system, with a mismatch in.

Higher. Specimen NAB Assessment

Chapter Summary. Mathematical Induction Strong Induction Recursive Definitions Structural Induction Recursive Algorithms

Determining the Accuracy of Modal Parameter Estimation Methods

Signal reconstruction algorithm based on a single intensity in the Fresnel domain

Kinetics of Particles. Chapter 3

3D FE Modeling Simulation of Cold Rotary Forging with Double Symmetry Rolls X. H. Han 1, a, L. Hua 1, b, Y. M. Zhao 1, c

Chapter 3 Kinematics in Two Dimensions; Vectors

Study Group Report: Plate-fin Heat Exchangers: AEA Technology

Review of the Roll-Damping, Measurements in the T-38 Wind Tunnel

FIZIKA ANGOL NYELVEN JAVÍTÁSI-ÉRTÉKELÉSI ÚTMUTATÓ

Yeu-Sheng Paul Shiue, Ph.D 薛宇盛 Professor and Chair Mechanical Engineering Department Christian Brothers University 650 East Parkway South Memphis, TN

Determination of ionic product constant of water (K w ) Handout 2014

Quantum Harmonic Oscillator, a computational approach

Effect of Paramagnetic Ions in Aqueous Solution for Precision Measurement of Proton Gyromagnetic Ratio

MODULE 1. e x + c. [You can t separate a demominator, but you can divide a single denominator into each numerator term] a + b a(a + b)+1 = a + b

NATIONAL RADIO ASTRONOMY OBSERVATORY COMPUTER DIVISION INTERNAL REPORT

APPLICATION OF THE BRATSETH SCHEME FOR HIGH LATITUDE INTERMITTENT DATA ASSIMILATION USING THE PSU/NCAR MM5 MESOSCALE MODEL

LAB III. CONDUCTIVITY & THE HALL EFFECT

6. Frequency Response

Q1. A string of length L is fixed at both ends. Which one of the following is NOT a possible wavelength for standing waves on this string?

Non-Linear Dynamic Behavior of Thin Rectangular Plates Parametrically Excited Using the Asymptotic Method, Part 1: Computation of the Amplitude

MATHEMATICS Higher Grade - Paper I

POLARISATION VISUAL PHYSICS ONLINE. View video on polarisation of light

TOPPER SAMPLE PAPER 2 Class XII- Physics

ECE 5318/6352 Antenna Engineering. Spring 2006 Dr. Stuart Long. Chapter 6. Part 7 Schelkunoff s Polynomial

205MPa and a modulus of elasticity E 207 GPa. The critical load 75kN. Gravity is vertically downward and the weight of link 3 is W3

Practical results of GPS/IMU/camera system calibration. Eija Honkavaara*, Risto Ilves**, Juha Jaakkola*

Principal Components

convection coefficient. The different property values of water at 20 C are given by: u W/m K, h=8062 W/m K

Dispersion Ref Feynman Vol-I, Ch-31

LEARNING : At the end of the lesson, students should be able to: OUTCOMES a) state trigonometric ratios of sin,cos, tan, cosec, sec and cot

Internal vs. external validity. External validity. This section is based on Stock and Watson s Chapter 9.

Applying Kirchoff s law on the primary circuit. V = - e1 V+ e1 = 0 V.D. e.m.f. From the secondary circuit e2 = v2. K e. Equivalent circuit :

EEO 401 Digital Signal Processing Prof. Mark Fowler

, which yields. where z1. and z2

BASD HIGH SCHOOL FORMAL LAB REPORT

1996 Engineering Systems Design and Analysis Conference, Montpellier, France, July 1-4, 1996, Vol. 7, pp

5 th grade Common Core Standards

Image Processing Adam Finkelstein & Tim Weyrich Princeton University

Performance Bounds for Detect and Avoid Signal Sensing

CESAR Science Case The differential rotation of the Sun and its Chromosphere. Introduction. Material that is necessary during the laboratory

Equilibrium of Stress

Department of Electrical Engineering, University of Waterloo. Introduction

ABSOLUTE TECHNIQUE FOR NEUTRON SOURCE CALIBRATION BY RADIATION INDUCED ACTIVITY

Current/voltage-mode third order quadrature oscillator employing two multiple outputs CCIIs and grounded capacitors

Mingqing Xing 1 School of Economics and Management, Weifang University, Weifang ,

CHAPTER 6 WORK AND ENERGY

ANALYTICAL SOLUTIONS TO THE PROBLEM OF EDDY CURRENT PROBES

( ) kt. Solution. From kinetic theory (visualized in Figure 1Q9-1), 1 2 rms = 2. = 1368 m/s

Cairo University Faculty of Engineering. Outline

TABLE OF CONTENTS. SUMMARY (ii) ILLUSTRATIONS (iv) 1.0 INI RODUCTION STANDARD TEST METHODS 1

Modeling the Nonlinear Rheological Behavior of Materials with a Hyper-Exponential Type Function

Biplots in Practice MICHAEL GREENACRE. Professor of Statistics at the Pompeu Fabra University. Chapter 13 Offprint

initially lcated away frm the data set never win the cmpetitin, resulting in a nnptimal nal cdebk, [2] [3] [4] and [5]. Khnen's Self Organizing Featur

Free Vibrations of Catenary Risers with Internal Fluid

Proceedings of the 2006 IASME/WSEAS International Conference on Continuum Mechanics, Chalkida, Greece, May 11-13, 2006 (pp )

A - LEVEL MATHEMATICS 2018/2019

the results to larger systems due to prop'erties of the projection algorithm. First, the number of hidden nodes must

ANSWER KEY FOR MATH 10 SAMPLE EXAMINATION. Instructions: If asked to label the axes please use real world (contextual) labels

Transactions on Engineering Sciences vol 19, 1998 WIT Press, ISSN

COASTAL ENGINEERING Chapter 2

THE PARTITION OF ENERGY INTO WAVES AND CURRENTS

CHAPTER 2 Algebraic Expressions and Fundamental Operations

February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA

Key words Shock waves. Dusty gas. Solid particles. Shock jump relations. Mach number

Drought damaged area

Rigid Body Dynamics (continued)

The Sputtering Problem James A Glackin, James V. Matheson

CS 477/677 Analysis of Algorithms Fall 2007 Dr. George Bebis Course Project Due Date: 11/29/2007

Transcription:

0th IMEKO TC4 ymsium n Braunschweig, GERMAY, 20, etember 2-4 Prjectin Miré Prfilmetry using iquid Crystal Digital Gratings Fumi Kbayashi, Center fr Otical Research and Educatin, Utsunmiya University; Yukitshi Otani, chl f Bimedical Engineering; Tru Yshizawa, aitama Medical University Abstract: A rjectin miré rfilmetry using a air f liquid crystal digital gratings (CDG) is alied t determine a ste height using tw liquid crystal digital gratings (CDGs). Miré cnturs are filtered ver t remve the images f the riginal grating frm the miré cunter when a airs f grating images are traveled t the same directin in-hase electrically. The new idea f this rert is rsed t vercme the 2π ambiguity f miré cntur by changing the fringe interval f miré cnturs which are adjusted the erid f CDGs. The measured accuracy has been achieved t ±0.7mm in the range f 96 mm f stee height by rjecting tw different f grating erids as 6.7 and 3.3 time f riginal grating.. Intrductin In this rert, we rse the rjectin miré rfilmetry by means a hase-shifting methd using a air f liquid crystal digital gratings (CDG) t measure a stee height r searated areas. Miré cnturs are filtered ver t remve the images f the riginal grating frm the miré cunter when a airs f grating images are traveled t the same directin inhase electrically. The new idea f this rert is rsed that we vercme the 2π ambiguity f miré cntur by dual grating erid methd by changing the fringe interval f miré cnturs which are adjusted the erid f CDGs. Recently, there are requirements fr three-dimensinal rfile with measuring seed and high satial reslutin. One f the useful methds is sterescic analysis. A light-sectining methd and a grating-rjectin methd are ne f the ular three-dimensinal rfilmetry. A shadw miré tgrahy is smart t remve the image f the riginal grating frm the miré cunters by mving the grating[]. We have succeeded t measure a ste height shae by means the frequency mdulatin methd fr shadw miré tye f tgrahy by rtating the grating [2]. A rjectin tye f miré tgrahy requires nly a small system which is made f a rjectr and a camera with gratings. Hwever there are rblems f measurement time and accuracy because f the effect f riginal gratings. We have rsed a structured C digital grating, esecially fr the grating rjectin methd[3].

0th IMEKO TC4 ymsium n Braunschweig, GERMAY, 20, etember 2-4 We rsed the three-dimensinal shae measurement by a hase-shifting grating rjectin methd by the structured C digital grating. In this aer, we rse a rjectin miré rfilmetry by hase-shifting using a air f liquid crystal digital gratings (CDG) t measure stee height areas and searated areas. 2. Thery and exerimental setu Prjectin Miré Prfilmetry Figure shws an exerimental setu a rjectin miré rfilmetry by means the hase-shifting using CDGs. A air f liquid-crystal digital gratings (CDGs) and imaging lenses are set arallel t the tical axis. One CDG wrks as a rjectin grating and the ther is reference grating. A defrmed grating image that is illuminated by light surce t s samle is verlaed t the reference grating and is rduced miré cnturs. The intensity Ix f the miré cntur at the bject with height h frm the reference lane n the CCD camera is written as ), Fig. : Otical cnfiguratin f rjectin miré rfilmetry I x I 0 2π 2π ah = + γ cs x+ γ cs x+ 4 m 2 2π + γ cs 2x+ ah 2 2π + γ cs ( + h ) ( ) + ( + ) m h m h ah, () where, the maximum f intensity I 0,the visibility f miré cntur, the distance between imaging lenses a, the itch f CDG, the distance between the imaging lens and reference lane, the distance between the imaging lens and CDG b, and the magnitude f the imaging lens at the reference lane m (= /b ), resectively. The height h is exressed by the gemetrical relatin by the fringe rder as, h = m. (2) a m In Eq. (), the first term is the bias, and the secnd t furth terms are the fundamental and secnd harmnic frequencies f the riginal gratings, the s-called nise terms, and the last

0th IMEKO TC4 ymsium n Braunschweig, GERMAY, 20, etember 2-4 term is the miré cnturs. We remve the secnd t furth terms in Eq. () and extract the miré cnturs. When the images f the grating atterns f the CDGs travel in the x- directin, synchrnized electrically with velcity V and integrating the exsure time t f the CCD camera, we can remve the nise grating images f the cmnents frm the secnd t furth terms in Eq. (). 3. te height and searated area measurement by the dual grating-erid methd The fringe rder is exressed using the hase f miré cnturs by the hase-shifting methd. It is easy t aly hase shifting by changing the hase f the CDG. Hwever, there is ambiguity, such as the integer f the fringe rder. It can be determined the fringe rder using the dual grating-erid methd and the recise hase by the shrt- and lnggrating erids because the sensitivity f miré cnturs deends un the grating erid f miré cnturs. Figure 2 shws the dual grating-erid methd, which cnsists f shrt- and lng-erids f grating. The lng erid f miré cntur cvers the dynamic range f the measurement in ne fringe, r a few erids in the case f an easy unwra area. First, the fringe rder f the shrt erid f miré cntur is calculated, and the hase φ f the shrt erid f miré cnturs can be easily unwraed using the fringe number as fllws: φ = 2 π + φ, (3) where indicates the hase f the shrt erid f miré cnturs, which is analyzed by the hase-shifting methd. The fringe rder f the shrt erid f miré cnturs can be determined as.5 One erid 40 8 0.5 0 0 0 30 50 70 90 0 3-0.5 - Reference 5 sitin Phase (r Fig. 2: Prsed technique using shrt- and lng-grating erids = 2 φ π + 4 ϕ ϕ, (4)

0th IMEKO TC4 ymsium n Braunschweig, GERMAY, 20, etember 2-4 where the hase f the lng erid f miré cnturs φ, the erid itch f grating, the hase data f the shrt-grating erid ϕ,, and Gaussian bracket [ ]. T cmensate fr the errr ε f the wraed hase within ± π 2 resulting frm the hase errr between the hase f the lng-grating erid and the hase t change the hase f the shrt-grating erid ϕ t ( φ ) between the hases f shrt- and lng-grating erids as fllws. φ =, we cmare the hase difference s ε = (5) φ s φ The determinatin arameter k is shwn as, k π = 2.. (6) If ε is larger than k, we need t crrect the wraed hase. We assume the next equatin with crrectin as j = and withut crrectin as j = 0. Finally, we btain the crrected hase as fllws. φ ε ϕ ' φ = j π + π + π ε 2 2 ϕ ϕ (7) Here, ε means the discriminant term, which gives the sign f ε. ε Figure 3 shws a simulated result in a case where randm errr at the lng-grating erid f the miré cnturs is added. The simulatin arameters used were: itch 6.08 mm er line fr the lng-grating erid, and 0.456 mm er line fr the shrt-grating erid. The simulatin is used as a maximum hase errr with 7% fr the lng-grating erid and a maximum hase errr with 5% fr the shrt -grating erid. As the errr is much larger than the hase value f the shrt-erid grating erid, the unwraed errr is bserved. Hwever, there are n unwraed errrs using determinatin arameter k. Our rsed methd is a rbust enugh measurement t determine the fringe rder as a crrectin f the unwra errr.

0th IMEKO TC4 ymsium n Braunschweig, GERMAY, 20, etember 2-4 hase( = 0.456 m m ) h ase ( = 6.08 mm) hase unw raed( = 0.456 m m ) crrected hase unw raed( = 0.456 m m ) 40 30 = 0.456 mm unwraed ith Prsed technique (crrected hase unwraed) 位相 (rad) 20 0 = 6.08 mm with errr 0-0 -0 0 0 20 30 40 50-20 Distance 高さ (mm) (mm) Fig. 3: imulatin f errr reductin by rsed technique 4. Exerimental results f ste height measurement Figure 4 shws the measured results f a stair mdel with three miré cnturs fr different heights. We calculated wraed miré cnturs with a high-grating erid by using the lwgrating erid f the miré data. The lw-grating erid has a itch f 3.04 mm. The highgrating erid has a itch f 0.456 mm and a itch f 0.228 mm. In the figure, 3.04 shws the individual result f a grating itch f 3.04 mm, and 0.456 Grating itch (mm) shws the result using a lwgrating itch f 3.04 mm and 3.04 te 0.456 a high-grating itch f 0.456 0.228 mm, 0.228 shws the result using a lw-grating itch f te 2 3.04 mm and a high-grating itch f 0.228 mm. We te 3 cmare each height in Table. Frm measurement 2 in this table, we estimate that the maximum height errr is less than 0.6 mm. 90 540 590 640 690 740 790 te 4 Pixels Fig. 4 Exerimental result f ste height

0th IMEKO TC4 ymsium n Braunschweig, GERMAY, 20, etember 2-4 5. Cnclusin We have rsed hase-shifting rjectin miré tgrahy using a air f liquid crystal digital gratings (CDGs) t measure ste height r searated areas. Our nvel idea is that we vercme the 2π ambiguity f miré cnturs and determined the fringe number by changing the fringe interval f miré cnturs, which are adjusted t the CDG erid in this case by 6.7 and 3.3 times the itc.h he rsed methd is t determine the fringe number f miré cnturs by dual grating erid methd and the recise hase by lng and shrt grating erids because the sensitivity f miré cnturs deends n the grating erid f miré cnturs. Finally, we have checked that the measured accuracy has been achieved t ±0.7mm in the range f 96 mm f stee height by rjecting tw different f grating erids as 6.7 and 3.3 time f riginal grating. [] H.Takasaki ; Al. Ot. 9, 4 (970).467-472. [2].H.Jin, T.Yshizawa, Y.Otani ; Otical Engineering, 40, 7 (200). 383-386. [3] K.Yamatani, Y.Otani, T.Yshizawa, H.Fujita, et.al : Prc.PIE Vl.3782 (999). 29-296.