Outline: Sampling
Outline: Sampling
Outline
FT of comb function Sampling Nyquist Condition sinc interpolation Truncation Aliasing
Sampling
As we saw with the CCD camera and the pinhole imager, the detector plane is not a continuous mapping, but a discrete set of sampled points. This of course limits the resolution that can be observed.
Limits of Sampling Finite # of data points Finite field of view High spatial frequency features can be missed or recorded incorrectly
Formalism of Sampling
Comb(x) =
( x n)
n=
dx
1; x 1 TopHat (x ) = 0; x >1
22.058 - lecture 5, Sampling and the Nyquist Condition
Formalism of Sampling
n= fov
2
+ fov 2
f (x ) (x nx)
fn =
+ fov 2 fov 2
f (x )(x n x)dx
Formalism of Sampling
f(x)
Comb(x/x)
TopHat(2x/fov)
Formalism of Sampling
f(x)
Comb(x/x)
TopHat(2x/fov)
Formalism of Sampling
{fn}
Frequency of Sampling
data = Table @ If @ Mod@ x , 7 D < 2 , 1 , - 1 D * Exp @ - Hx - 128 L^ 2 1000 D, 8x , 1 , 256 <D; -> All , Axes -> False , AspectRatio -> 1 4 , ListPlot @ data , 8PlotJoined -> True , PlotRange PlotStyle -> Thickness @ 0.005 D<D
FT k
22.058 - lecture 5, Sampling and the Nyquist Condition
Sampling
1
0.5
2 -0.5
-1
In[19]:=
p2 = Plot @ Cos @ 2 Pi 0.2 Hx - 1 LD, 8x , 1 , 6 <, 8PlotPoints -> 512 , PlotStyle -> 8Thickness @ 0.01 D, RGBColor @ 1 , 0 , 0 D<<D
In[11]:=
Sampling
1 1
0.5
0.5
2 -0.5
6 -0.5
-1
-1
In[19]:=
p2 = Plot @ Cos @ 2 Pi 0.2 Hx - 1 LD, 8x , 1 , 6 <, 8PlotPoints -> 512 , PlotStyle -> 8Thickness @ 0.01 D, RGBColor @ 1 , 0 , 0 D<<D
In[11]:=
In[18]:=
p1 = Plot @ Cos @ 2 Pi 1.2 Hx - 1 LD, 8x , 1 , 6 <, 8PlotPoints -> 512 , PlotStyle -> 8Thickness @ 0.01 D, RGBColor @ 0 , 0 , 1 D<<D
Sampling
1 1
0.5
0.5
2 -0.5
6 -0.5
-1
-1
In[19]:=
p2 = Plot @ Cos @ 2 Pi 0.2 Hx - 1 LD, 8x , 1 , 6 <, 8PlotPoints -> 512 , PlotStyle -> 8Thickness @ 0.01 D, RGBColor @ 1 , 0 , 0 D<<D
In[11]:=
In[20]:=
p3 = Plot @ Cos @ 2 Pi 2.2 Hx - 1 LD, 8x , 1 , 6 <, 8PlotPoints -> 512 , PlotStyle -> 8Thickness @ 0.01 D, RGBColor @ 0 , 1 , 0 D<<D
Sampling
1 1
0.5
0.5
2 -0.5
6 -0.5
-1
-1
Nyquist theorem: to correctly identify a frequency you must sample twice a period. So, if x is the sampling, then /x is the maximum spatial frequency.
s1 (t ) = cos(2 f o t )
We know the Fourier Transform of this
cos(2 f o t ) ( ( f f o ) + ( f + f o ))
1 What happens when we sample this at a rate of t
1 where has the units of Hz t and t = the dwell of the sampled signal .
Note, a short cut has been taken, and left off the integrals that are needed to sample with the delta function.
22.058 - lecture 5, Sampling and the Nyquist Condition
fo <
fs
= f n = Nyquist frequency
2f0
-fs
fs
In this case, there is no overlap and regardless of the complexity of this spectrum (think of having a number or continuum of cosine functions), the frequency spectrum correctly portrays the time evolution of the signal.
22.058 - lecture 5, Sampling and the Nyquist Condition
Case 2
fo >
fs
-fs
0 2f0
Now the frequency spectrum overlaps and look what happens to our picture,
fs
f '0 = fs f 0
-fs
fs
Nyquist Theorem
Consider what happens when there is a complex spectrum. Then the entire spectrum overlaps. Either way, the frequency spectrum does not correspond to a correct picture of the dynamics of the original time domain signal. Nyquist Theorem: In order to correctly determine the frequency spectrum of a signal, the signal must be measured at least twice per period.
cos(2 f o nt ) ; t = 1
So
fs
2 f o n cos fs
Now let
fo > fs fo = fs ( fs fo ) 1 4 2 4 3 f
(cont)
2 n cos [ f s f ] fs 2 n f cos2 n fs f cos2 n fs
2 nf cos fs
Therefore we can see that it is not the Fourier Transform that fails to correctly portray the signal, but by our own sampling process we misrepresented the signal.
22.058 - lecture 5, Sampling and the Nyquist Condition
Limitations of Sampling
{fn}
Limitations of Sampling
perfect sampling
Limitations of Sampling
perfect sampling
Limitations of Sampling
perfect sampling
Reciprocal Space
real space
reciprocal space
Sampling
60
50
40
30
20
10
0 0 10 20 30 40 50 60
real space
under sampled
Sampling
120
100
80
60
40
20
0 0 20 40 60 80 100 120
real space
zero filled
Sinc interpolation
Filtering We can change the information content in the image by manipulating the information in reciprocal space.
Fourier Convolution
Deconvolution
I (x) = O(x) PSF (x) + N (x) {
noise
1 PSF (k )
Deconvolution (cont)
This appears to be a well-balanced function, but look what happens in a Fourier space however:
i(k ) = O(k ) PSF ( k ) 1 n( k ) + PSF (k ) PSF ( k )
Recall that the Fourier Transform is linear where PSF(k) << n(k), then noise is blown up. The inverse filter is ill-conditioned and greatly increases the noise particularly the high-frequency noise since
PSF (k ) 0 at high k 1444 4 24444 3 normally
( N (k ))
( N)
S
Electron beam
Deflection magnets
Electron beam is scanned over the screen by deflection magnets and the intensity of the beam is modulated to give the intensity at each point. The refresh rate is 30 frames/second. By the sampling theorem, if the frequency is f o , independent information 1 is available once every 2 f seconds. Need information
frames lines pixels 500 650 = 9.75 10 6 impulse sec sec frame line f 4.875 MHz
30
where
(x nX ) =
1, x = nX 0, x nX
One trick to this is to express the comb function as a Fourier series expansion, not the transform.
Comb(x) = where
i 2 nx 1 2 X dx a n = Comb(x)e X x 2 x
a e
n n =
i 2 nx
i 2 nx 1 2 X a n = (x )e dx X x 2 14 4 244 3
select the x = 0 po int, note this is only true of a Dirac delta function .
an =
1 X
1 Comb(x) = e X
n =
i 2 nx
2n 2 (x nX ) k X X
2 X
2 X