Я пытаюсь запрограммировать этот интеграл:
E (s, t, z) - излучение фотона от сцены изображения, в В моем случае это массив 3D NumPy с различными значениями фотонов, размещенными по всему объему (сцена изображения в градациях серого). Идея заключается в том, что я должен быть в состоянии интегрировать по объему изображения и генерировать 4D светлое поле, которое будет наблюдаться от камеры светового поля на некотором расстоянии. Камера сфокусирована на плоскости z_0, например, средней плоскости массива NumPy. U, V представляет направление света angular, т. Е. Также может рассматриваться как выборка по основному объективу камеры. s ', t' - пространственные отсчеты на датчике камеры, M - увеличение, которое больше нуля.
Очевидно, что дельта-функция исключает интеграл, и вы просто суммируете по E с различными сдвигами в s и t для каждого z. Поскольку не все сдвиги в конечном итоге будут целыми числами, я не совсем уверен в наилучшем подходе к программированию этого интеграла в python. Спасибо!!
Источник статьи для интеграла https://www.osapublishing.org/DirectPDFAccess/54129FF9-D674-C25B-9FA4479D45E57A0E_396313/oe-26-18-22574.pdf?da=1&id=396313&seq=0&mobile=no