2028.sagews

import numpy as np
import pylab as lab

Esta función toma la dirección de una imagen que represente una placa con un agujero. La imagen debe estar en formato de escala de grises. La función muestra la imagen, y muestra el patron de difracción correspondiente.

def difracta(dir_placa):
    placa=lab.imread(dir_placa)#Se importa imagen placa como una matriz
    a=matrix_plot(placa)#Se grafica la placa
    a.show(figsize=10)#Se muestra la placa
    F=np.fft.fft2(placa)#Se calcula la transformada de Fourier
    Fs=np.fft.fftshift(F)#Se mueve la zona de frecuencia cero al centro de la imagen
    c=matrix_plot(abs(Fs))#Se grafica el modulo de la imagen
    c.show(figsize=10)#Se muestra

difracta(DATA+'cuad.png')

difracta(DATA+'circ.png')

difracta(DATA+'tria.png')

difracta(DATA+'equi.png')

difracta(DATA+'hexa.png')

difracta(DATA+'reji.png')

difracta(DATA+'nump.png')