мой предыдущийвопрос (Функция интеграла интенсивности в питоне)

Вы можете увидеть дифракционную модель на изображении ниже:

Я хочу вычислить интеграл интенсивности в каждом пикселе (квадрате), поэтому я не могу использовать R и Theta в качестве переменной. Как я могу сделать это в координате X-Y.

Наша функция:

вместо греха (тэта) мы можем использовать:

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

Другие константы:

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

когда вы строите график, результат выглядит примерно так: (изображение выше - вид сверху)

метод в предыдущем разделе: вычислить интегрирование функции в (0.0, dist) и после этого * (2 * np.pi).х) который х = кa * np.sin (theta), но теперь я хочу интегрировать в каждый пиксель. который предыдущий метод не работает, потому что это координата X-Y, а не полярная.