mi anteriorpregunta (Función integral de intensidad en python)

Puedes ver el modelo de difracción en la imagen a continuación:

Quiero calcular la integral de intensidad en cada píxel (cuadrado), por lo que no puedo usar R y Theta como variable. ¿Cómo puedo hacer esto en coordenadas X-Y?

Nuestra función:

en lugar de pecado (theta) podemos usar:

sintheta= (np.sqrt((x)**2 + (y)**2)/(np.sqrt((x)**2 + (y)**2 + d**2)))

Otras constantes:

lamb=550*10**(-9)
k=2.0*np.pi/lamb
a=5.5*2.54*10**(-2)
d=2.8

cuando traza la función, el resultado es algo como esto: (La imagen de arriba es una vista desde arriba)

El método en el tema anterior: calcular la integración de la función en (0.0, dist) y después de eso * (2 * np.pix) que x = ka * np.sin (theta), Pero ahora quiero integrar en cada píxel. que el método anterior no funciona, porque esta es la coordenada X-Y, no polar.