03. Diffraction Integrals (Fresnel + Fraunhofer propagation, Point Spread Function, Fourier optics)

Lecture notes: drive.google.com/drive/folder...
ERRATA: at 7:29, the condition for Fresnel diffraction should read 1/2*(rho/lambda)^4/3 less than z/lambda. (Thanks to Mohammed ElKabbash for pointing this out)
0:54 Huygens' principle and the Rayleigh-Sommerfeld integral
3:00 Fresnel and Fraunhofer approximations
5:58 Overview of diffraction integrals
6:40 Validity of the Fresnel and Fraunhofer approximations, Fresnel number
8:59 Relation between ray optics and wave optics
9:46 Deriving the transmission function of a thin convex lens
10:50 Deriving the thin lens equation using wave optics
11:44 Deriving the magnification using wave optics
13:18 Introducing the pupil function, deriving the Point Spread Function (PSF)
16:16 Effect of the PSF on the imaging resolution
18:25 Coherent and incoherent imaging by convolving the object with the PSF
20:56 Coherent and incoherent imaging using transfer functions
22:34 The Fourier transforming property of a lens
24:59 Fourier filtering, dark field imaging, phase contrast imaging
Other sources:
Introduction to Fourier Optics, Joseph W. Goodman
Microscopy: Measuring The Point Spread Function (Nico Stuurman)
• Microscopy: Measuring ...

Пікірлер: 37

  • @joem8251
    @joem82513 жыл бұрын

    I am blown away by how quickly, clearly, and concisely you present this material!

  • @shane3379
    @shane33793 жыл бұрын

    I am only 5 minutes into this video and have already learned more than in my entire advanced optics module.... brilliant video, this is the future of learning

  • @mathematicality
    @mathematicality3 жыл бұрын

    I'm still learning a lot of fourier optics concepts...going over LARGE, demoralizing textbooks😑 and you just managed to compress so much into one video. Incredible!

  • @jacobvandijk6525
    @jacobvandijk65255 ай бұрын

    @ 0:55 Here, in H's Principle, the source of light is called a field too. One could call this starting field U(0,0). @ 1:40 If you also put in U(0,0), then you have an expression for "the propagator" from QFT. Here, U(0,0) = 1.

  • @ovieigherebuo9317
    @ovieigherebuo93174 жыл бұрын

    Thanks a whole semesters lecture in one coherent video

  • @Kolibril
    @Kolibril3 жыл бұрын

    Thank you, this is really condensed and purified content

  • @martinmillischer6083
    @martinmillischer60833 жыл бұрын

    Last year I took a class at photonics west on Fourier Optics. It was 8h long and $600. In 25 min I learn more with your video. Thank you so much for your work.

  • @Lovearbre
    @Lovearbre6 жыл бұрын

    Great Video! Thank you!

  • @iosianchad2180
    @iosianchad21802 жыл бұрын

    Good material for reviewing wave optics after dropping it for 20 years.

  • @rajatsaxena1703
    @rajatsaxena1703 Жыл бұрын

    Great video!

  • @esmirhodzic981
    @esmirhodzic9815 жыл бұрын

    very nice video

  • @mementomori6734
    @mementomori6734 Жыл бұрын

    18:06 from angular spectrum method we have that the resolution is limited to lambda/2. But if we have a lens with a NA>1 we can have a resolution

  • @SanderKonijnenberg

    @SanderKonijnenberg

    Жыл бұрын

    Generally, NA is defined NA=n*sin(theta), where n is the refractive index of the immersion medium. So NA>1 is typically achieved by immersing the sample in a medium with high refractive index, e.g. water or immersion oil. If the refractive index of a medium is n, then the wavelength of light in that medium is lambda/n (where lambda is the wavelength in vacuum). By decreasing the wavelength we increase the resolution. Therefore, the resolution is still larger than lambda/2 if we define lambda to be the wavelength in the immersion medium, but it can be smaller than lambda/2 if we define lambda to be the wavelength in vacuum.

  • @adurgh
    @adurgh3 жыл бұрын

    any reference textbooks on this that can be recommended?

  • @SanderKonijnenberg

    @SanderKonijnenberg

    3 жыл бұрын

    I'd go with J. Goodman's 'Introduction to Fourier Optics'

  • @daguaishouxd
    @daguaishouxd4 жыл бұрын

    Just a minor question, at 15:30, I thought that in F.T. shift theorem the shift in F.T. has the same sign as the exponent of the phase shift, i.e. F{f(x)e^(-2i pi ax)}(x',y') = F{f(x)}(x'-a, y'-a). I guess the later equations will still hold if we define M as |z_i|/|z_0| without the negative sign... Any thoughts?

  • @SanderKonijnenberg

    @SanderKonijnenberg

    4 жыл бұрын

    If the Fourier transform is defined as F{f(x)} (x') = int f(x)e^(-2 pi i x x') dx then F{f(x) e^(- 2 pi i ax)} (x') = int f(x)e^(-2 pi i x (x'+a)) dx = F{f(x)} (x'+a). Perhaps what you're thinking of is the inverse Fourier transform? If you multiply F{f(x)} (x') with e^(-2 pi i a x'), then inverse Fourier transforming get you f(x-a). Physically, we also know that a real image generated by a single lens is inverted, so one way or another, M must turn out to be negative.

  • @daguaishouxd

    @daguaishouxd

    4 жыл бұрын

    @@SanderKonijnenberg Yes you are definitely right. It's easier to think of (x'+a) as a new dummy variable then it's quite obvious. I did confused it with an e^(-2 pi i a x') shift in x' domain which will lead to f(x-a) in x domain. Thanks for clarifying!

  • @user-si1zn3ir7x
    @user-si1zn3ir7x2 жыл бұрын

    What does global phase factor in 10:36 mean? Thank you for this wonderful video!

  • @SanderKonijnenberg

    @SanderKonijnenberg

    2 жыл бұрын

    Thanks for your comment. 'Global phase factor' means in this case a phase factor that does not depend on (x,y), and therefore does not affect the way that the field propagates.

  • @arturocatalanogonzaga7874
    @arturocatalanogonzaga78747 ай бұрын

    Hi! Thanks a lot for the extremely helpful videos! Could you please clarify how and why, at minute 21:16, you went from having an Object in the (xo,yo) object plane to the (x/M,y/M) plane? why is it the lens plane that gets demagnified? shouldn't it be (xo/M,yo/M)?

  • @SanderKonijnenberg

    @SanderKonijnenberg

    7 ай бұрын

    My apologies for the confusion: in this slide, (x,y) denotes the image plane coordinates (not the lens plane). If we ignore the effects of the PSF (i.e. assume that it's a delta peak), then the image I(x,y) should be a magnified copy of |O(x_o,y_o)|^2. That is, I(x,y)=|O(x/M,y/M)|^2, so that the object field at object coordinates (x_o,y_o) is mapped to image coordinates (x,y)=(M*x_o,M*y_o).

  • @arturocatalanogonzaga7874

    @arturocatalanogonzaga7874

    6 ай бұрын

    @@SanderKonijnenberg thank you so much for the reply!

  • @nebulae_wanderer
    @nebulae_wanderer2 жыл бұрын

    My 5am savior

  • @thetatheta9926
    @thetatheta99263 жыл бұрын

    Whats the difference between Rayleigh-Sommerfeld integral and Huygen-Fresnel integral?

  • @SanderKonijnenberg

    @SanderKonijnenberg

    3 жыл бұрын

    That's a good question, I never thought about it before. Section 3.7 of Joseph Goodman's 'Introduction to Fourier Optics' (can be found online) says: "The Huygens-Fresnel principle, as predicted by the first Rayleigh-Sommerfeld solution (see Eq. (3-40)), can be expressed mathematically as follows: [...]". It appears that the Huygens-Fresnel integral is a heuristically derived principle, whereas the Rayleigh-Sommerfeld solution is more rigorously derived from the wave equation, and from which the Huygens-Fresnel principle follows. But ultimately, it appears they describe the same thing.

  • @thetatheta9926

    @thetatheta9926

    3 жыл бұрын

    @@SanderKonijnenberg Thanks for the reply. I'm actually trying to accurately determine how the field of a Gaussian beam is transformed after impinging on a reflective metasurface (i.e. is it still a Gaussian beam?). I applied The Huygen-Fresnel diffraction integral mainly because it appeared to be a simpler version of the RS integral. I guess I have to checkout Goodman.

  • @grandaurore
    @grandaurore8 ай бұрын

    3:58 why in approximation on complex exponential, we should change the sign?

  • @SanderKonijnenberg

    @SanderKonijnenberg

    8 ай бұрын

    Because exp(ik * lambda/2)=exp(i*pi)=-1 (recall k=2pi/lambda)

  • @grandaurore

    @grandaurore

    8 ай бұрын

    thanks, very clear, the minus come from exp(i*pi). I was misleading to expansion generating the minus sign.

  • @user-jq5ut4jg5d
    @user-jq5ut4jg5d5 ай бұрын

    Can you share the MATLAB code?

  • @SanderKonijnenberg

    @SanderKonijnenberg

    5 ай бұрын

    The Google drive link in the video description also contains a zip-file with the Matlab codes. Hopefully those include what you're looking for.

  • @mightbin
    @mightbin2 жыл бұрын

    at 2:23 , i think the complex form of the wave function i thinkl is : z/r2 * ei(kr−ωt), why do you say it is z/r2 * e ikr, and ommit the −ωt part ? Thanks a ot.

  • @SanderKonijnenberg

    @SanderKonijnenberg

    2 жыл бұрын

    Yes, technically the field has a time-dependence, but the advantage of complex notation is that it allows you to omit the time-dependent complex exponential e^{-i\omega t}, see my video kzread.info/dash/bejne/ZWVkmZSjhsmvhag.html at 7:39 - 9:33 . For a monochromatic field, we typically don't care about how each point oscillates in time, but only about the phase difference between fields at different points in space.

  • @mightbin

    @mightbin

    2 жыл бұрын

    really helpful. Thank you so much.@@SanderKonijnenberg

  • @superdogmeatmeat
    @superdogmeatmeat4 жыл бұрын

    Nou Sander.. optica is toch niet mijn ding :O

  • @ErichAmMeer
    @ErichAmMeer2 жыл бұрын

    why are that many books so long and boring without clear explanation,disappointing until I found this