07.11.2014 Views

Application and Optimisation of the Spatial Phase Shifting ...

Application and Optimisation of the Spatial Phase Shifting ...

Application and Optimisation of the Spatial Phase Shifting ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

3.4 <strong>Spatial</strong> phase shifting 91<br />

frequencies. The difference is solely that those regions <strong>of</strong> <strong>the</strong> frequency plane where <strong>the</strong>re is no signal,<br />

<strong>and</strong> hence relatively little spectral power, are quite a bit more noisy.<br />

ν N<br />

1.57<br />

0.785<br />

ν y<br />

0<br />

0<br />

ν<br />

0 1 2 3 x /ν 0x 4<br />

-0.785<br />

–ν N<br />

2 3 4|0 1 2 ν x /ν 0x -1.57<br />

Fig. 3.33: Left: bsc(ν x ,ν y ) for (3.18) as calculated from (3.73), with 0 black <strong>and</strong> 2π white; note that ν x =ν y = 0 is<br />

in <strong>the</strong> centre <strong>of</strong> <strong>the</strong> image. Right: bsc(ν x ) for (3.18) (black) <strong>and</strong> (3.19) (white); average <strong>of</strong> 50 rows from<br />

<strong>the</strong> small black frame on <strong>the</strong> left.<br />

To <strong>the</strong> right, <strong>the</strong> plot <strong>of</strong> bsc(ν x ) as output by (3.18) (black) <strong>and</strong> (3.19) (white) does indeed show that both<br />

compare quite well with <strong>the</strong> corresponding graph in Fig. 3.13. The high noise around ν x = 0 <strong>and</strong> ν x = 2ν N<br />

reflects <strong>the</strong> suppression <strong>of</strong> I b , i.e. <strong>the</strong> fact that ~ ~ S ( 0 ) = C(<br />

0)<br />

= 0.<br />

It is interesting to note that, in agreement<br />

with <strong>the</strong> larger absolute filter output <strong>of</strong> (3.18), <strong>the</strong> susceptibility to noise is indeed somewhat lower than<br />

for (3.19); but this affects a frequency range that produces large errors anyway, so that <strong>the</strong> difference in<br />

performance will be very small. For this case <strong>of</strong> α=90°, we <strong>the</strong>refore conclude that <strong>the</strong> utilisation <strong>of</strong><br />

different sets <strong>of</strong> pixels for different representations <strong>of</strong> <strong>the</strong> 3-sample 90° formula does not invalidate <strong>the</strong><br />

<strong>the</strong>oretical considerations in 3.2.2.3.<br />

For α=120°, we use an interferogram with a power spectrum as in Fig. 3.29 on <strong>the</strong> left; this time, <strong>the</strong><br />

signal sideb<strong>and</strong>s cover a smaller part <strong>of</strong> <strong>the</strong> frequency plane. When this interferogram is processed with<br />

(3.17), we can expect a qualitative behaviour resembling that in Fig. 3.31 because <strong>the</strong> sampling formulae<br />

are both derived from (3.15). As to be seen from Fig. 3.34, this is indeed <strong>the</strong> case; again <strong>the</strong> highfrequency<br />

preference <strong>of</strong> C<br />

~ ( ν ) is clearly visible. The distribution <strong>of</strong> <strong>the</strong> phase lag has a mean <strong>of</strong> 89.6° <strong>and</strong><br />

a st<strong>and</strong>ard deviation <strong>of</strong> 22.0°.<br />

x<br />

Fig. 3.34: From left to right: ~ 2 ~ ~<br />

2<br />

~ ~ 2<br />

( ν , ν ) ; I ( ν , ν ) ⋅ S ( ν ) <strong>of</strong> (3.17); I ( ν , ν ) ⋅ C(<br />

ν ) <strong>of</strong> (3.17); pixel<br />

I x y<br />

x y x<br />

histogram <strong>of</strong> phase lag between I(x,y)¡S x (n) <strong>and</strong> I(x,y)¡C x (n) <strong>of</strong> (3.17).<br />

x y x

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!