Application and Optimisation of the Spatial Phase Shifting ...
Application and Optimisation of the Spatial Phase Shifting ...
Application and Optimisation of the Spatial Phase Shifting ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
3.2 <strong>Phase</strong>-shifting ESPI 63<br />
∞<br />
~ ~<br />
S'( ) I ( )<br />
*<br />
0 = ν S ( ν ) dν<br />
∫<br />
−∞<br />
∞<br />
x x x<br />
~ ~<br />
C'( ) I ( )<br />
*<br />
0 = ν C ( ν ) dν ,<br />
∫<br />
−∞<br />
x x x<br />
(3.35)<br />
where tilde denotes <strong>the</strong> Fourier transforms <strong>and</strong> <strong>the</strong> sign convention [Bro87]<br />
∞<br />
~ f ( ν x ): = ∫ f ( x)exp( + 2 π i ν x x)<br />
dx<br />
(3.36)<br />
−∞<br />
is adopted, i.e. <strong>the</strong> phase runs forward in <strong>the</strong> Fourier transform.<br />
It is seen from (3.35) that <strong>the</strong> spectrum <strong>of</strong> I(x) is weighted, or filtered, by <strong>the</strong> spectra <strong>of</strong> S(x) <strong>and</strong> C(x),<br />
which is why we have called <strong>the</strong>m filter functions. We will <strong>the</strong>refore refer to S<br />
~ ( ν ) <strong>and</strong> C<br />
~ ( ν ) as filter<br />
spectra. Since I(x), S(x) <strong>and</strong> C(x) are real functions with Hermitian Fourier transforms, we can simplify<br />
<strong>the</strong> integrals (3.35) to [Fre90a]<br />
x<br />
x<br />
∞<br />
∫<br />
*<br />
x x x<br />
S' ( 0) = 2 Re I ~ ( ν ) S ~ ( ν ) dν<br />
0<br />
∞<br />
∫<br />
*<br />
x x x<br />
C' ( 0) = 2 Re I ~ ( ν ) C ~ ( ν ) dν .<br />
0<br />
(3.37)<br />
This is, <strong>the</strong> filter outputs are indeed composed <strong>of</strong> all input spatial frequencies that may be present in I(x),<br />
with weights determined by <strong>the</strong> moduli <strong>of</strong> <strong>the</strong> filter spectra, S<br />
~ ( ν ) <strong>and</strong> C<br />
~ ( ν )<br />
latter also as filter responses. With our initial choice <strong>of</strong> S(x) <strong>and</strong> C(x), we have<br />
∞<br />
∫<br />
S'( 0) = I ( x) ⋅ − sin 2πν<br />
xdx<br />
−∞<br />
∞<br />
∫<br />
0x<br />
C'( 0)<br />
= I ( x) cos2πν xdx ,<br />
−∞<br />
0x<br />
x<br />
x<br />
; we will refer to <strong>the</strong>se<br />
(3.38)<br />
being proportional to <strong>the</strong> Fourier sine <strong>and</strong> cosine transforms [Bra87, p. 17] <strong>of</strong> I(x) at <strong>the</strong> frequency ν 0x ,<br />
<strong>and</strong> <strong>the</strong> spectral descriptions read:<br />
∞<br />
∫<br />
S'( 0) = 2Re I ~ ( ν ) iδ ( ν −ν ) dν<br />
0<br />
∞<br />
∫<br />
x x 0x x<br />
C'( 0) = 2Re I ~ ( ν ) δ ( ν −ν ) dν .<br />
0<br />
x x 0x x<br />
(3.39)<br />
This models <strong>the</strong> ideal case that we can evaluate <strong>the</strong> signal over an infinite amount <strong>of</strong> space, which leads to<br />
unity filter responses at <strong>the</strong> nominal frequency ν 0x <strong>and</strong> perfect suppression <strong>of</strong> all o<strong>the</strong>r ν x .