Como se menciona en la pregunta, una lente delgada producirá en su plano focal la transformada de Fourier del campo óptico en su pupila, posiblemente multiplicada por un término de fase cuadrático. Sin embargo, para entender cómo se relaciona esto con la formación de imágenes en la imagen de la óptica ondulatoria, tenemos que dar un paso atrás, y observar la situación de forma más general. Bajo la aproximación paraxial, la propagación de un campo óptico puede modelarse con la integral de difracción de Fresnel: \begin{equation} U^\prime(x,y) = \frac{e^{i k z}}{i \lambda z} \mathrm{exp}\left[ \frac{i \pi (x^2 + y^2)}{\lambda z}\right] \ldots\\ \times \iint_{-\infty}^{\infty} U(\xi, \eta) \mathrm{exp}\left[ \frac{i \pi (\xi^2 + \eta^2)}{\lambda z}\right] \mathrm{exp}\left[ \frac{-i 2 \pi (x \xi + y \eta)}{\lambda z}\right] d\xi d\eta \end{equation} donde $U(\xi,\eta)$ es un campo óptico, $U^\prime(x,y)$ es el campo después de la propagación por una distancia $z$ y $\lambda$ y $k$ son la longitud de onda y el número de onda, respectivamente.
En el caso de una lente fina, una transparencia en contacto con la lente, y una distancia de propagación igual a la distancia focal $f$ podemos representar el campo de entrada como $$ U(\xi, \eta) = t_A(\xi, \eta) \mathrm{exp}\left[ \frac{-i \pi (\xi^2 + \eta^2)}{\lambda f}\right] $$ donde $t_A$ es la transmisión de amplitud de la transparencia, y el término de fase cuadrática es la curvatura del frente de onda introducida por una lente delgada de longitud focal $f$ . Si se introduce esto en la integral de difracción anterior, se ve que, cuando $z = f$ la integral se reduce a una transformada de Fourier y tenemos \begin{aligned} U^\prime(x,y) &= \frac{e^{i k z}}{i \lambda z} \mathrm{exp}\left[ \frac{i \pi (x^2 + y^2)}{\lambda z}\right] \iint_{-\infty}^{\infty} t_A(\xi, \eta) \mathrm{exp}\left[ \frac{-i 2 \pi (x \xi + y \eta)}{\lambda z}\right] d\xi d\eta \\ {} &= \frac{e^{i k z}}{i \lambda z} \mathrm{exp}\left[ \frac{i \pi (x^2 + y^2)}{\lambda z}\right] \mathcal{F}[t_A](x,y) \end{aligned} donde $\mathcal{F}$ es la transformada de Fourier. No lo digo explícitamente, pero puedes suponer que las transformadas de Fourier que escribo están siempre debidamente escaladas. En este caso, si la FT se define para tomar funciones de $(\xi, \eta)$ y funciones de retorno de la frecuencia espacial $(\alpha, \beta)$ , debe asumir el escalado implícito $(\alpha, \beta) \rightarrow (\frac{x}{\lambda z},\frac{y}{\lambda z})$ .
Ahora, no lo derivaré aquí porque las integrales son enormes, pero si usas la primera ecuación que escribí para propagar algún campo de objetos por una distancia $f$ y luego aplicar la modificación del frente de onda mediante una lente delgada de longitud focal $f$ y se propagan a otra distancia $f$ Si no se puede hacer nada, se verá que los términos de fase cuadrática se cancelan entre sí, y el campo resultante es exactamente la transformada de Fourier del campo del objeto, sin ni siquiera el término de fase cuadrática que se obtiene si el objeto está directamente contra la lente. Si tienes problemas con esto, ten en cuenta la identidad de la transformada de Fourier para la doble FT de una función; esto hace que la derivación sea sencilla.
De forma más general, esta derivación puede aplicarse a una serie arbitraria de elementos ópticos y distancias de propagación. Con suficiente esfuerzo, se puede demostrar que, para un sistema óptico paraxial descrito por una matriz ABCD, un campo óptico se propaga a través del sistema por: \begin{equation} U^\prime(x,y) = \frac{e^{i k L_0}}{i \lambda B} \mathrm{exp}\left[ \frac{i \pi D(x^2 + y^2)}{\lambda B}\right] \ldots\\ \times \iint_{-\infty}^{\infty} U(\xi, \eta) \mathrm{exp}\left[ \frac{i \pi A (\xi^2 + \eta^2)}{\lambda B}\right] \mathrm{exp}\left[ \frac{-i 2 \pi (x \xi + y \eta)}{\lambda B}\right] d\xi d\eta \end{equation} donde $L_0$ es la longitud efectiva del camino óptico a través del eje óptico del sistema.
Esto, por supuesto, sigue siendo válido sólo para un sistema óptico coherente. Una forma de pensar en esto en el contexto de un sistema de imágenes (como un ojo o una cámara) es que sólo se aplica al campo debido a un único punto de la escena que se está fotografiando. La imagen final puede obtenerse propagando coherentemente el campo de cada punto del objeto, tomando la magnitud al cuadrado del campo resultante para obtener su intensidad, y sumando después las intensidades de cada punto del objeto.
Así, supongo que se podría afirmar que vemos una superposición de transformadas de Fourier de cada punto del objeto, en lugar de ver directamente una transformada de Fourier. De hecho, la imagen en su retina no mira como la imagen que obtienes si tomas una escena cotidiana y la transformas de Fourier en tu ordenador. Sin embargo, las lentes hacer realizar transformadas de Fourier en campos ópticos. Sin embargo, cuando se considera un sistema de imágenes, hay que tener en cuenta dónde está el campo que se transforma, en relación con el objetivo . En general, este campo es pas el campo en el objeto que está mirando; es el campo a cierta distancia delante de su pupila, y en una situación del mundo real, no es simplemente un campo coherente de un punto de origen, sino una superposición incoherente de campos de cada punto en su campo de visión.
En la práctica, esto significa que la imagen incoherente rara vez se simula con la integral ABCD anterior. Este tipo de cálculo es útil para los sistemas de imagen coherentes (un telescopio es un buen ejemplo, si sólo se habla de estrellas y no de objetos extendidos), pero en el caso incoherente es mucho más sencillo simular la imagen aplicando puramente la MTF/OTF como una convolución o un filtro lineal. Sin embargo, incluso en este caso, el cálculo sigue basándose en una transformada de Fourier.