Как реализовать анизотропное размытие из геометрических формул? ⇐ Python

[Anonymous]

Мне нужно применить размытие/фильтр к 2D-сетке, где параметры размытия варьируются для каждого пикселя на основе геометрических расчетов. У меня есть:

• Изображение: 2D-сетка значений
• Углы обзора: 2D-сетка, в которой каждый пиксель имеет угол θ (в диапазоне 22–25°)

Для параметров размытия мне нужно:

• Вертикальное направление (y): константа σ = 2 пикселя для всех местоположения
• Горизонтальное направление (x): варьируется в зависимости от пикселя, рассчитывается по формулам, показанным на прикрепленном изображении:

Формулы учитывают кривизну Земли и рассчитывают, как удлиняется зона действия спутникового датчика при разном просмотре углы. Дано:

Код: Выделить всё

R_E = 6371 km (Earth's radius)
h_VIIRS = 833 km (satellite altitude)
θ₁ = viewing angle from nadir
σ_nadir = 742 m (footprint at nadir)

Вот ссылка для иллюстрации математических расчетов и задачи.
Как правильно преобразовать эти геометрические формулы в код, чтобы применить пространственно изменяющееся анизотропное размытие по Гауссу?
Решения на Python, R, MATLAB или Julia приветствуются.>

Подробнее здесь: https://stackoverflow.com/questions/798 ... c-formulas