Transformada fraccional de Fourier en el caso de un plano imagen inclinado

Ingenierı́a y Ciencia, ISSN 1794–9165 Volumen 6, número 12, julio-diciembre de 2010, páginas 51–64 Transformada fraccional de Fourier en el caso de un plano imagen inclinado Fourier fracional transformar no caso de um plano imagem inclinado Fraccional Fourier transform in the case of an inclined image plane C. O. Torres1, L. Mattos2, C. Jiménez3, J. Castillo4 y Y. Torres5 Recepción:13-oct-2009 /Modificación:27-abr-2010 /Aceptación:27-abr-2010 Se aceptan comentarios y/o discusiones al artı́culo Resumen La conocida fórmula de difracción de Fresnel relaciona la distribución de amplitud compleja de una onda en el plano objeto (campo ondulatorio de entrada) con la distribución de amplitud compleja de la onda en el plano imagen (campo ondulatorio de salida) cuando se trata de propagación en el espacio libre; esto significa que si los planos objeto e imagen son paralelos entre sı́, el sistema imagen correspondiente se dice que es un sistema lineal invariante a desplazamiento (LSI). Esta propiedad ventajosa es esencial para el desarrollo de técnicas de imagen sensitivas a fase; sin embargo, si el plano imagen 1 Doctor en Ciencias fı́sicas, cesartorres@unicesar.edu.co, docente, grupo LOI, Universidad Popular del Cesar, Valledupar–Colombia. 2 Magı́ster en Ingenierı́a electrónica, lorenzo.mattos@gmx.fr, docente, grupo LOI, Universidad Popular del Cesar, Valledupar–Colombia. 3 Magı́ster en Matemática aplicada, carlosj114@gmail.com, docente, grupo GIMA, Universidad de La Guajira, Riohacha–Colombia. 4 Magı́ster en Matemática aplicada, jacas68@yahoo.es, docente, grupo GIMA, Universidad de La Guajira, Riohacha–Colombia. 5 Docteur en Optique et Traitement du Signal, ytorres@uis.edu.co, docente investigador, grupo GOTS, Universidad Industrial de Santander, Bucaramanga–Colombia. Universidad EAFIT 51| Transformada fraccional de Fourier en el caso de un plano imagen inclinado está inclinado con respecto al haz incidente, la distancia efectiva de propagación cambiará sobre el plano imagen, consecuentemente el sistema imagen será no invariante a desplazamiento. En este artı́culo es propuesta una extensión del formalismo de la difracción de Fresnel al caso de un plano imagen inclinado utilizando la transformada de Fourier de orden fraccional. Palabras claves: difracción de Fresnel, plano imagen inclinado, transformada fraccional de Fourier. Resumo A bem conhecida equação de difração de Fresnel estabelece uma relação entre a distribuição de amplitude complexo de uma onda no plano do objeto (campo de onda de entrada) com a distribuição de amplitude complexa da onda no plano da imagen (campo de onda de saı́da) no caso da propagação em espaço libre. Isto significa que se os planos objeto e imagen são paralelos uns aos outros, o sistema imagem correspondente é dito ser um sistema linear invariante ao deslocamento (LSI). Esta propriedade vantajosa é essencial para o desenvolvimento de técnicas de imagen sesı́vel á fase, No entanto, se o plano da imagem está inclinada respeito do feixe incidente, a distancia efetiva da propagação mudara sobre o plano imagen, consequentemente o sistema imagem não é invariante ao deslocamento. Neste trabalho uma extensão do formalismo da difração de Fresnel para o caso de um plano imagem inclinado é proposto usando a transformada de Fourier de ordem fracionária. Palavras chaves: difração de Fresnel, plano de imagem inclinada, Fourier fracionária transformar. Abstract The well-known Fresnel integral relates a known complex wave defined in the object plane (the input wave field) to the observable complex wave (the output wave field) defined in the image plane after free-space propagation; this means that if the object and image plane are parallel to each other, corresponding imaging system is said to be linear-shift-invariant (LSI). This advantageous property was essential for the development of phase sensitive imaging techniques; however, if the image plane is inclined with respect to the incident beam, the effective propagation distance will vary over the image plane, consequently, the imaging system is not shiftinvariant. In this paper an extension of the theoretical formalism of Fresnel diffraction to the case of an inclined image plane is proposed using the fractional Fourier transform. Key words: Fresnel diffraction, inclined image plane, fractional Fourier transform. |52 Ingenierı́a y Ciencia, ISSN 1794–9165 C. O. Torres, L. Mattos, C. Jiménez, J. Castillo y Y. Torres 1 Introducción Las mediciones de superficies rugosas de componentes metálicos son consideradas de gran importancia en el campo de maquinado mecánico. Desde 1972, cuando Sprague [1] aplicó técnicas de speckle blanco, muchos investigadores han estudiado este aspecto para realizar mediciones de superficies rugosas, [2, 3, 4, 5, 6, 7, 8]. Los patrones de speckle contienen información sobre las estructuras microscópicas de una superficie rugosa. Para estudiar la relación entre las propiedades estadı́sticas de los patrones de speckle, formados por superficies rugosas, se debe obtener la distribución de amplitud compleja del campo de speckle. Para lograr este propósito, las ecuaciones correspondientes a la difracción, que describen la propagación de las ondas difractadas de la superficie rugosa al plano de observación, deben ser construidas; en estudios preliminares, las ecuaciones de difracción de Fresnel y Fraunhofer fueron utilizadas para calcular el campo del patrón de speckle [6, 7], con la consideración de que el plano de observación es paralelo al plano del objeto en todos los montajes experimentales, a excepción de la referencia [8]. El propósito de este artı́culo es derivar las ecuaciones de difracción correspondientes para describir la difracción de una superficie inclinada rugosa reflectiva, bajo iluminación de una onda plana coherente en términos de la transformada fraccional de Fourier. Se considera la diferencia entre el modelo de difracción por superficie reflectiva (RSDM) y el modelo de difracción por transmisión de la apertura (TADM); encontrándose una ecuación general de difracción en términos de la transformada de Fourier de orden fraccional. Aplicando las aproximaciones de campo cercano y de campo lejano a una ecuación general de difracción con consideración de rotación del sistema coordenado, se obtienen las ecuaciones de difracción, tanto para el campo cercano como para el lejano. El significado fı́sico de estas ecuaciones y la obtención como casos particulares de las ecuaciones de Fresnel y Fraunhofer son también presentadas. 2 2.1 Integral básica de difracción para una superficie inclinada Integral de difracción de Fresnel–Kirchhoff La figura 1 representa el diagrama esquemático del modelo TADM. La figura 1(a) muestra que la luz de una fuente puntual P2 es difractada por transmisión Volumen 6, número 12 53| Transformada fraccional de Fourier en el caso de un plano imagen inclinado de la apertura, P0 es un punto de observación inicial en la región de difracción. P2 n n r 21 P1 P1 Σ Σ δ r 01 r 01 k21 a) P0 b) P0 Figura 1: modelo de difracción por transmisión de la apertura: a) Fuente puntual incidente. b) Luz paralela incidente Luego la integral de difracción de Fresnel–Kirchhoff se puede escribir como Z Z exp (ikr01 ) ds , (1) U ′ (P1 ) U (P0 ) = r01 Σ donde Σ es la superficie de la integral y se escoge como apertura ficticia  una  − → − → − →− → cos n , r ( 01 )−cos( n , r 21 ) 1 A exp (ikr21 ) ′ curvada. Matemáticamente, U (P1 ) = iλ r21 2 puede ser vista como una fuente puntual secundaria virtual en el punto P1 → → sobre Σ, r21 y r01 son los módulos de los radio vectores − r 21 y − r 01 , respecti→ − vamente; n es la normal a Σ en el punto P1 . La expresión (1) es únicamente aplicable para el caso de incidencia de fuente puntual (ver figura 1(a)) e inaplicable cuando una onda plana incidente ilumina una apertura de transmisión (ikr21 ) (ver figura 1(b)). Pero como el término A expr21 en U ′ (P1 ) es exactamente la distribución del campo luminoso de la superficie de la apertura Σ, y se representa por U (P1 ), entonces (1) toma la forma U (P0 ) = 1 iλ Z Z U (P1 ) Σ   → → → → n,− r 01 ) − cos (− n,− r 21 ) exp (ikr01 ) cos (− ds . r01 2 (2) La ecuación (2) puede ser usada en el caso de incidencia de luz paralela, es decir, la fuente de luz puntual P2 está ahora infinitamente alejada de la |54 Ingenierı́a y Ciencia, ISSN 1794–9165 C. O. Torres, L. Mattos, C. Jiménez, J. Castillo y Y. Torres apertura de transmisión. Visto desde la figura 1(b), el vector de onda incidente → − → k 21 tiene la misma dirección de − r 21 , luego, (2) se puede reescribir como U (P0 ) = 1 iλ Z Z  − →  − → → − − → cos ( n , r ) − cos n , k 21 01 exp (ikr01 )   ds . U (P1 ) r01 2 Σ  (3) Dando la expresión de U (P1 ), el campo ondulatorio U (P0 ) en el punto de observación P0 puede ser # obtenido. El ángulo δ en la figura 1(b) y el término " → − → − → − → − cos( n , r 01 )−cos n , k 21 en (3) representa el ángulo de difracción y el factor 2 de inclinación respectivamente para el modelo de difracción por transmisión de la apertura. 2.2 Integral básica de difracción de RSDM La figura 2 es la ilustración esquemática del modelo RSDM. Una onda plana está incidiendo sobre una superficie que produce scattering débil con un vector → − de onda k 21 y es reflejada y difundida. P0 r 01 k’21 k21 γ δ γ P1 Plano tangente Superficie real n Plano medio Figura 2: modelo de difraccón (RSDM) por superficie reflectiva → − El vector de onda reflejado en el punto P1 es k′ 21 , comenzando en un punto arbitrario sobre la superficie actual. La atención debe ser puesta sobre el ángulo δ en la figura 2, el cual es definido como el ángulo de difracción para el modelo de difracción por superficie reflectiva, pero es diferente Volumen 6, número 12 55| Transformada fraccional de Fourier en el caso de un plano imagen inclinado del definido anteriormente en el modelo de difracción por transmisión de la apertura; especı́ficamente, el ángulo δ en la figura 1(b) está definido entre el → − vector de onda incidente k 21 y la lı́nea P1 P0 , mientras que δ, en la figura → − 2, es el ángulo entre el vector de onda reflejado k′ 21 y la lı́nea P1 P0 . Como resultado de esto, el factor de inclinación utilizado para el modelo de difracción por transmisión de la apertura no puede ser directamente usado para el modelo de difracción por superficie reflectiva. Para obtener el factor de inclinación de una difracción por superficie reflectiva, el factor de inclinación # "  − →  → → → cos(− n ,− r 01 )−cos − n , k 21 para el modelo de difracción por transmisión de la 2 "  − → # → → → n , k ′ 21 cos(− n ,− r 01 )−cos − apertura debe ser reemplazado por . De la figura 2 se 2  −  − →  →  → → → n es la normal a la supern , k 21 , donde − nota que cos − n , k′ 21 = − cos − ficie en el punto P1 . Combinando con (3), se obtiene la integral de difracción básica para una superficie reflectiva 1 U (P0 ) = iλ Z Z   −  → → − − → → − exp (ikr01 )  cos ( n , r 01 ) + cos n , k 21  ds , U (P1 ) r01 2 Σ (4) → − → → → → → donde (− n , r01 ) es el ángulo entre − n y− r 01 , (− n , k 21 ) es el ángulo entre − n → − → − y k 21 , r01 es el módulo del radio vector r 01 y P0 es el punto de observación de la región de difracción. 3 3.1 Ecuación general de la difracción para un plano inclinado Distribución de amplitud compleja de una superficie clásica La figura 3 ilustra el caso cuando una onda plana coherente está incidiendo sobre una superficie rugosa reflectiva, con un ángulo θ con respecto al eje z. El plano medio del objeto es representado por el plano coordenado xOy, donde el origen O y Pm (x1 , y1 , 0) es la proyección de un punto genérico P1 (x1 , y1 , z1 ) situado sobre la superficie convencional, z1 es una medida de la altura de la superficie en P1 . El punto P0 (x0 , y0 , z0 ) es un punto de observación en el régimen de difracción, sin pérdida de generalidad, asumiendo que la luz incidente tiene amplitud constante A0 y paralela al plano xOz, se nota que es |56 Ingenierı́a y Ciencia, ISSN 1794–9165 C. O. Torres, L. Mattos, C. Jiménez, J. Castillo y Y. Torres P0(x 0, y0, z0) Luz incidente θ z y P1(x 1, y1, z1) Luz incidente y1 θ Pm(x 1, y1, 0) x1 0 Figura 3: onda plana incidiendo sobre una superficie reflectiva un objeto de reflexión de pura fase, ignorando la absorción de la luz y otros factores, y haciendo que la fase en el origen O sea cero, se tiene la fase en P1 ϕ (P1 ) = k (x1 sin θ − z1 cos θ) , donde k = 2π λ es el número de onda y λ es la longitud de la onda incidente. Aquı́ se obtiene la distribución de amplitud compleja de la luz de onda incidente sobre la superficie   2π (5) U (P1 ) = A0 exp i (x1 sin θ − z1 cos θ) . λ 3.2 Cálculo del factor de inclinación Suponiendo que la superficie es z = h (x, y), donde h (x, y) es el origen de la altura de la superficie del plano medio, se pueden obtener las componentes → del vector normal unitario descendente − n en el punto P1 (x1 , y1 , z1 ) − → n = {nx , ny , nz } = donde hx = ∂h(x1 ,y1 ) , ∂x      hy hz hx q , ,q ,q 2 2 2 2 2 hx + hy + 1 hx + hy + 1 hx + h2y + 1  hy = ∂h(x1 ,y1 ) ∂y notando que → − − → r 01 = {x1 − x0 , y1 − y0 , z1 − z0 } y k 21 = k {sin θ, 0, − cos θ} , Volumen 6, número 12 57| Transformada fraccional de Fourier en el caso de un plano imagen inclinado se tiene: → → cos (− n,− r 01 ) =         hx (x1 − x0 ) + hy (y1 − y0 ) − (z1 − z0 ) q , q 2 2 2 h2x + h2y + 1 (x1 − x0 ) + (y1 − y0 ) + (z1 − z0 )      h sin θ + cos θ  − → x → cos − n , k 21 = q .  h2 + h2 + 1  x y  Sustituyendo (6) y (7) en " #  − →  → → → cos(− n ,− r 01 )−cos − n , k 21 2 (6) (7) , se obtiene el factor de inclinación del modelo de difracción por superficie reflectiva. 3.3 Ecuación general de la difracción para una superficie inclinada Se tiene r01 = q (x1 − x0 )2 + (y1 − y0 )2 + (z1 − z0 )2 y el elemento de área en (3), en el punto P1 , q ds = h2x + h2y + 1dx1 dy1 . (8) (9) Introduciendo (5) y (6)-(9) en (4), se obtiene la ecuación general de la difracción para una superficie inclinada   Z Z 2π 1 (x1 sin θ − z1 cos θ) × A0 exp i i2λ λ  q D  2 2 2 exp ik (x1 − x0 ) + (y1 − y0 ) + (z1 − z0 ) q × 2 2 2 (x1 − x0 ) + (y1 − y0 ) + (z1 − z0 )   h (x − x ) + h (y − y ) − (z − z ) 0 y 1 0 1 0  qx 1 + hx sin θ + cos θ dx1 dy1 , 2 2 2 (x1 − x0 ) + (y1 − y0 ) + (z1 − z0 ) U (P0 ) = (10) donde D define la región de integración en el plano medio. |58 Ingenierı́a y Ciencia, ISSN 1794–9165 C. O. Torres, L. Mattos, C. Jiménez, J. Castillo y Y. Torres 4 4.1 Ecuación de difracción de campo cercano Arreglo geométrico del experimento La figura 4 muestra el arreglo geométrico del experimento para el análisis de la difracción de una superficie reflectiva de scattering débil. n P0 z Luz incidente y θ r 01 0’ L θ P1(x 1, y1, z1) r 0m ξ Pm(x 1, y1, 0) x Figura 4: configuración geométrica 4.2 Derivación de la ecuación de difracción de campo cercano A partir de (10) es difı́cil calcular la difracción de la luz que debe ser aproximadamente considerada para utilidades prácticas. Todas las aproximaciones son hechas sobre la base de la aproximación paraxial: (1) la distancia L es muy grande comparada con las dimensiones lineales máximas de la región de iluminación, y (2) el punto de observación P0 está cerrado sobre el eje z. Considerando la relación entre el plano de observación y el plano objeto, el factor r01 en el término exponencial exp ikr01 en la ecuación (10) puede ser expandido en serie como z12 − 2z1 z0 , 2r0m donde el término cuadrático retenido como k es un número grande, mientras los términos de mayor orden son despreciados, q r0m = (x1 − x0 )2 + (y1 − y0 )2 + z02 , r01 ≈ r0m + Volumen 6, número 12 59| Transformada fraccional de Fourier en el caso de un plano imagen inclinado es la distancia entre los puntos P0 y Pm . De acuerdo a la aproximación paraxial (1) y (2), r01 puede ser además refinado por las siguientes tres razonaz12 0m ; (b) rz0m ≈ zr00 ; (c) bles aproximaciones: (a) ignorando el término pequeño 2r0m p (x2 +y2 ) 1 y0 ) x20 + y02 + z02 es la distancia de donde r = r0m ≈ r0 + 12r0 1 − (x1 x0r+y 0 0 P0 del origen O. Luego r01 en el término exponencial se reduce a
x21 + y12 (x1 x0 + y1 y0 ) z0 z1 r01 ≈ r0 + − − . (11) 2r0 r0 r0 Si se asume que el área de iluminación sobre la superficie del objeto es pequeña, la siguiente desigualdad puede ser utilizada
x21 + y12 (x1 x0 + y1 y0 ) z0 z1 − − , (12) r0 ≫ 2r0 r0 r0 r01 en el denominador de (10) puede ser escrita como q r01 ≈ r0 = x20 + y02 + z02 . (13) Asumiendo que la superficie de scattering débil satisface la aproximación de pendiente pequeña, luego las inecuaciones | nx | ≪ | nz |, | ny | ≪ | nz | o equivalentemente | hx | ≪ 1, | hy | ≪ 1 son válidas. Los términos de la derivación parcial en (10) pueden ser simplificados como z0 hx (x1 − x0 ) + hy (y1 − y0 ) − (z1 − z0 ) q + hx sin θ + cos θ ≈ p 2 + cos θ . 2 2 2 x0 + y02 + z02 (x1 − x0 ) + (y1 − y0 ) + (z1 − z0 ) (14) Sustituyendo (11)-(14) en (10), se obtiene la expresión aproximada de la ecuación de difracción para una superficie inclinada reflectiva |60 # " z0 A0 p + cos θ × U (P0 ) = U (x0 , y0 , z0 ) = i2λ x20 + y02 + z02   p   exp ik x20 + y02 + z02 Z Z 2π p (x1 sin θ − z1 cos θ) × exp i λ x20 + y02 + z02 D ) ( ) ( 2π x1 x0 + y1 y0 x21 + y12 π exp −i p × exp i p λ x20 + y02 + z02 λ x20 + y02 + z02 Ingenierı́a y Ciencia, ISSN 1794–9165 C. O. Torres, L. Mattos, C. Jiménez, J. Castillo y Y. Torres ( z0 z1 2π exp −i p 2 λ x0 + y02 + z02 ) dx1 dy1 . (15) Aplicando la conversión de coordenadas a (15), x0 = L sin θ + ξ cos θ, y0 = η y z0 = L cos θ −ξ sin θ e introduciendo la función apertura P (x, y), (15) puede ser transformada en   p # " exp ik L2 + ξ 2 + η 2 A0 L cos θ − ξ sin θ p p + cos θ × U (ξ, η) = i2λ L2 + ξ 2 + η 2 L2 + ξ 2 + η 2   Z Z ∞ 2π P (x1 , y1 ) exp i (x1 sin θ) × λ −∞ " ( # ) ) ( L cos θ − ξ sin θ x21 + y12 π 2π cos θ + p z1 exp i p × exp −i λ λ L2 + ξ 2 + η 2 L2 + ξ 2 + η 2 ) ( 2π x1 (L sin θ + ξ cos θ) + y1 η p dx1 dy1 exp −i λ L2 + ξ 2 + η 2 (16) Esta es la ecuación de difracción de campo cercano para una superficie rugosa inclinada.     L √cos θ−ξ sin θ z1 , (16) pueSi U0 (x1 , y1 ) = P (x1 , y1 ) exp −i 2π λ cos θ + 2 2 2 L +ξ +η de ser reescrita como   p 2 + ξ 2 + η2 L exp ik A0 L cos θ − ξ sin θ p p + cos θ × U (ξ, η) = i2λ L2 + ξ 2 + η 2 L2 + ξ 2 + η 2 ) ( ) ( 2 2 π (L sin θ + ξ cos θ) + η 2 π (L sin θ + ξ cos θ) + η 2 p p exp i × exp −i λ λ L2 + ξ 2 + η 2 L2 + ξ 2 + η 2 ) (   Z Z ∞ 2π x21 + y12 π U0 (x1 , y1 ) exp i × (x1 sin θ) exp i p λ λ L2 + ξ 2 + η 2 −∞ ) ( 2π x1 (L sin θ + ξ cos θ) + y1 η p dx1 dy1 . exp −i λ L2 + ξ 2 + η 2 " # (17) La ecuación (17) corresponde a una transformación fraccional de Fourier de orden α dada por la expresión   p 2 + ξ 2 + η2 exp ik L A0 L cos θ − ξ sin θ p p U (ξ, η) = + cos θ × i2λ L2 + ξ 2 + η 2 L2 + ξ 2 + η 2 " Volumen 6, número 12 # 61| Transformada fraccional de Fourier en el caso de un plano imagen inclinado ( 2 π (L sin θ + ξ cos θ) + η 2 p exp −i λ L2 + ξ 2 + η 2 )    2π (x1 sin θ) , U0 (x1 , y1 ) exp i F λ α (18) es decir, la funciónde apertura o pupila se encuentra montada sobre una fase de la forma exp i 2π λ (x1 sin θ) correspondiendo a una onda esférica; y la distribución de amplitudd compleja en el plano de observación corresponde a una transformada fraccional de Fourier de orden fraccional de la distribución de amplitud compleja del plano objeto. Donde "
  π
# 2 + η2 − α exp i −i ξ 2 × exp F α [U0 (x1 , y1 )] = 2π sin α 2 tan α "
#   Z ∞Z ∞ −i x1 2 + y1 2 i (x1 ξ + y1 η) exp exp U0 (x1 , y1 ) dx1 dy1 2 tan α sin α −∞ −∞ Considerando en (18) las aproximaciones necesarias que permitan describir la difracción de campo cercano, se puede obtener la forma estándar de la ecuación de difracción de Fresnel para una superficie inclinada y reflectiva, es decir, # " A0 L cos θ − ξ sin θ p + cos θ × U (ξ, η) = i2λ L2 + ξ 2 + η 2   p ) ( 2 exp ik L2 + ξ 2 + η 2 π (L sin θ + ξ cos θ) + η 2 p p × exp −i λ L2 + ξ 2 + η 2 L2 + ξ 2 + η 2        2π L cos θ−ξ sin θ √ exp i z , cos θ + F α P (x1 , y1 ) exp −i 2π (x sin θ) 1 1 λ λ 2 2 2 L +ξ +η (19)
de la misma forma, considerando difracción de campo lejano x21 + y12 Máx ≪ λL π , se llega a la forma estándar de la ecuación difracción de Fraunhofer para una superficie inclinada y reflectiva, es decir,  p  2 + ξ 2 + η2 exp ik L A0 L cos θ − ξ sin θ p p + cos θ × U (ξ, η) = i2λ L2 + ξ 2 + η 2 L2 + ξ 2 + η 2        2π L cos θ−ξ sin θ √ F P (x1 , y1 ) exp −i 2π exp i z cos θ + (x sin θ) 1 1 λ λ 2 2 2 # " L +ξ +η |62 (20) Ingenierı́a y Ciencia, ISSN 1794–9165 C. O. Torres, L. Mattos, C. Jiménez, J. Castillo y Y. Torres 4.3 Comparaciones con las ecuaciones de difracción de Fresnel y Fraunhofer En el caso especial de incidencia normal, es decir, θ = 00 , (19) puede además ser redefinida utilizando las aproximaciones: o n p p 2 +η 2 2 + ξ 2 + η2 en L2 + ξ 2 + η 2 ≈ L + ξ 2L L , por el factor exp i 2π (a) λ (19). p (b) L2 + ξ 2 + η 2 ≈ L, para otros factores en (19). Luego, (19) puede ser transformada en      4π π ξ 2 + η2 A0 exp (ikL) α . (21) P (x1 , y1 ) exp −i z1 F exp −i U (ξ, η) = iλ L λ L λ De (20) se obtiene A0 exp (ikL) U (ξ, η) = F iλ L    4π P (x1 , y1 ) exp −i z1 . λ (22) Las ecuaciones (21) y (22) son justamente las expresiones clásicas correspondientes a la difracción de Fresnel y de Fraunhofer para plano de difracción y plano de observación no inclinados. 5 Conclusiones Partiendo de la integral de Difracción de Fresnel–Kirchhoff, se ha obtenido una integral básica de difracción para una superficie inclinada reflectiva donde la diferencia en el cálculo del factor de inclinación entre el modelo RSDM y el TADM ha sido tenido en cuenta. Como resultado fundamental se ha encontrado una ecuación general de difracción para una superficie inclinada y reflectiva, bajo iluminación de una onda plana coherente, considerando los efectos de rotación del sistema coordenado sobre la aproximación de campo cercano e introduciendo la consideración de inclinaciones pequeñas; la difracción de campo cercano y del campo lejano son obtenidas y presentadas como casos particulares de (8). Por tanto, la expresión derivada es más general que Volumen 6, número 12 63| Transformada fraccional de Fourier en el caso de un plano imagen inclinado las ecuaciones de difracción clásicas de Fresnel y Fraunhofer. Se debe resaltar que los campos de luz microscópicos emergentes están contenidos en estas ecuaciones de difracción, y hacen referencia a la distribución de altura de la superficie considerada. Referencias [1] Robert A. Sprague. Surface Roughness Measurement Using White Light Speckle. Applied Optics, ISSN 0003–6935, 11(12), 2811–2816 (1972). Referenciado en 53 [2] P. J. Chandley. Surface roughness measurements from coherent light scattering. Optical and Quantum Electronics, ISSN 0306–8919, 8(4), 323–327 (1976). Referenciado en 53 [3] H. M. Pedersen. Object–roughness dependence of partially developed speckle patterns in coherent light . Optics Communications, ISSN 0030–4018, 16(1), 63–67 (1976). Referenciado en 53 [4] H. Fuji, T. Asakura and Y. Shindo. Measurement of surface roughness properties by means of laser speckle techniques. Optics Communications, ISSN 0030–4018, 16(1), 68–72 (1976). Referenciado en 53 [5] Hitoshi Fujii and Toshimitsu Asakura. Roughness measurements of metal surfaces using laser speckle. Journal of the Optical Society of America B, ISSN 0740– 3224, 67(9), 1171–1176 (1977). Referenciado en 53 [6] Ulf Persson. Real time measurement of surface roughness on ground surfaces using speckle–contrast technique. Optics and Lasers in Engineering, ISSN 0143– 8166, 17(2), 61–67 (1992). Referenciado en 53 [7] Lisa C. Leonard and Vincent Toal. Roughness measurement of metallic surfaces based on the laser speckle contrast method . Optics and Lasers in Engineering, ISSN 0143–8166, 30(5), 433–440 (1998). Referenciado en 53 [8] S. L. Toh, H. M. Shang and C. J. Tay. Surface–roughness study using laser speckle method . Optics and Lasers in Engineering, ISSN 0143–8166, 29(2– 3), 217–225 (1998). Referenciado en 53 |64 Ingenierı́a y Ciencia, ISSN 1794–9165