Астигматизм

Немного отличается от того, что у вас, но суть +/- та же - изображение смазывается.

Есть два вида:
1 - тангенциальный (вдоль радиуса)
2 - радиальный (он же саггитальный) - вдоль окружности

Принцип:
Как мы уже знаем, идеально все работает только вдоль оптической оси, вне ее лучи фокусируются как правило вне соответствующей оси (если провести от исходной точки через центр условной линзы) из-за чего объекты так или иначе искажаются (вспоминаем геометрические искажения) и становятся смазанными.

Можно было бы копнуть глубже и рассказать про S и T фокусы, но я тут про общие принципы и реализацию рассказываю, а не школьный курс физики)))

Алгоритм

Cчитаем параметры пространства, и когда считаем точки для свертки - применяем его

struct AstigCtx {
 float2 e_r;
 float2 e_t;
 float s_r;
 float s_t;
};

struct Parameters {
 float2 invSize;
 float2 opticalCenter;
 float astigK;
};

inline AstigCtx astigSetup(float2 uv, Parameters params) {
 AstigCtx C;
 float ar = params.invSize.y / params.invSize.x;
 float2 df = uv - params.opticalCenter;
 float2 dfIso= float2(df.x * ar, df.y);
 float rhoN = clamp(length(dfIso) * 2.0f, 0.0f, 1.0f);

 float2 rIso = (length(dfIso) > 1e-6f) ? normalize(dfIso) : float2(1.0f, 0.0f);
 C.e_r = normalize(float2(rIso.x / ar, rIso.y));
 C.e_t = float2(-C.e_r.y, C.e_r.x);

 float mag = fabs(params.astigK) * (rhoN * rhoN);
 float a = 1.0f + mag;
 C.s_r  = (params.astigK >= 0.0f) ? (1.0f / a) : a;
 C.s_t  = (params.astigK >= 0.0f) ? a : (1.0f / a);
 return C;
}

inline float2 astigApply(float2 d, AstigCtx C) {
 float dr = dot(d, C.e_r);
 float dt = dot(d, C.e_t);
 return dr * C.s_r * C.e_r + dt * C.s_t * C.e_t;
}

Разбираюсь с build123d, чтобы процедурно генерить ракушки

У меня на второй работе происходит жесткий церебральный коитус с экспортом и женитьбой разных кодеков на разных андроидофонах, поэтому пока посты про видео на паузе.

Вот зато ракушку сделал в перерыве на build123d

Кома

Смещение пучка из центра пятна боке - получаются свеого рода кометы.

Принцип

Из-за неидеальности оптики вне центральной оптической оси лучи фокусируются неравномерно (см на картинке).

Соотвественно, чем дальше от оси, тем сильнее эффект.

Алгоритм

1. Как и с астигматизмом, сначала рассчитываем параметры радиально-тангенциального поля в точке (можно расширить функцию для астигаматизма
2. Для каждой точки свертки пересчитываем, куда она смещается

struct AstigCtx {
    float2 e_r;  // radial unit in uv metric
    float2 e_t;  // tangential unit (perp)
    float  s_r;  // radial scale
    float  s_t;  // tangential scale
    float  rhoN; // normalized field height (0 center .. 1 edge)
};

inline AstigCtx astigSetup(float2 uv, constant Parameters& P) {
    AstigCtx C;

    float ar = P.invSize.y / P.invSize.x; // width/height
    float2 oc = P.opticalCenter;
    float2 df = uv - oc;
    float2 dfIso = float2(df.x * ar, df.y);
    float  rhoN = clamp(length(dfIso) * 2.0f, 0.0f, 1.0f);

    float2 rIso = (length(dfIso) > 1e-6f) ? normalize(dfIso) : float2(1.0f, 0.0f);
    float2 e_r  = normalize(float2(rIso.x / ar, rIso.y));
    float2 e_t  = float2(-e_r.y, e_r.x);

    float mag   = fabs(P.astigK) * (rhoN * rhoN);
    float a     = 1.0f + mag;
    float s_r   = (P.astigK >= 0.0f) ? (1.0f / a) : a;
    float s_t   = (P.astigK >= 0.0f) ? a : (1.0f / a);

    C.e_r = e_r;
    C.e_t = e_t;
    C.s_r = s_r;
    C.s_t = s_t;
    C.rhoN = rhoN;

    return C;
}

inline float2 comaApply(float2 d, AstigCtx C, float comaK) {
    if (fabs(comaK) < 1e-6f) return d;
    // Project onto radial/tangential axes
    float dr = dot(d, C.e_r);
    float dt = dot(d, C.e_t);
    // Magnitude grows ~ rho^3 across the field (third-order behavior)
    float mag = comaK * (C.rhoN * C.rhoN * C.rhoN);
    // Quadratic skew in pupil coordinates → "comet" shape in image PSF
    float2 delta = ( (dr*dr - 0.5f*dt*dt) * C.e_r
                   + (dr*dt) * C.e_t );
    return d + mag * delta;
}

Вообще, как видите, сверточные эффекты тоже опираются на линейную оптику, но, если для идеального фокуса мы можем использовать исключительно пересчет координат, то для объектов вне фокуса необходимо применять свертку и тут все крутится вокруг ядра этой свертки (ее формы, центра и тд).

Тут в соседнем чате скинули статью про настроечные таблицы и я вспомнил, что давно когда-то пробовал их использовать для реверсинжинирига цветокоррекций и эффектов, но они справлялись не очень (не для того их делали), поэтому пришлось сделать свою.

Из минусов - тут нет разного сатурейшена по цветам, поэтому для некоторых эффектов она не подходила, но подходил слайс (y)UV.

PS: серое поле между шахматной доской и угловым градиентом - тоже шахматная доска 1х1. Зачем-то она тоже была нужна, как и все остальные поля - не помню уже точно (лет 10 назад это было)

надо будет на досуге сделать новую версию - исправленную и дополненную

Полевая кривизна
При смещении от центра начинает усиливаться расфокус (на видео утрированно)

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

Алгоритм
Считаем радиус ядра свертки в зависимости от расстояния до центра и применяем его

Продольная цветовая абберация
Изменение оттенка у объектов до и после плоскости фокуса

Пример не совсем корректный, но для корректного примера нам нужна картинка/видео с картой глубин - тут мы просто делаем имитацию

Принцип
Из-за разного преломления лучей разной длины волны либо красные (недолет фокуса), либо синие (перелет) лучи будут вылетать из точки фокуса, что даст определнный оттенок.

Алгоритм
Считаем отдельные смещения для сэмплинга каждого канала. Потом собираем RGB

Яблоки релизнули иос26. Теперь телефон будет греться и садиться просто по факту своего существования 🤠

Сегодня в обед таких вот зверей наблюдал на отливе.

Я сейчас с работами разберусь и буду дальше про видео вещать

Сидел наблюдал за физарумом и нашел что он может объединяться вот в такие пузыри и ползать ими