Какова связь между сигмой в лапласианском гауссиане и двумя сигмами в разнице гауссианов?

12

Я понимаю, что фильтр Лапласа Гаусса может быть аппроксимирован фильтром Разности Гаусса, и что отношение двух сигм для последнего должно быть 1: 1,6 для наилучшего приближения. Однако, я не уверен, как две сигмы в Разнице Гауссианов относятся к сигме для лапласиана Гаусса. Является ли меньшая сигма в первом равно сигме второго? Большая сигма? Или отношения это что-то еще?

зрительно-кинетический
источник
> Я понимаю, что фильтр Лапласа Гаусса может быть аппроксимирован фильтром Разности Гаусса, и что отношение двух сигм для последнего должно быть 1: 1,6 для наилучшего приближения. извините с какой ссылкой вы знали это?
Привет, я думаю, что этот вопрос подойдет здесь - area51.stackexchange.com/proposals/86832/… Это также поддержит сообщество. Спасибо.
Рой

Ответы:

10

Я понимаю, что фильтр Лапласа Гаусса может быть аппроксимирован фильтром Разности Гаусса, и что отношение двух сигм для последнего должно быть 1: 1,6 для наилучшего приближения

Теоретически, чем меньше соотношение между двумя сигмами, тем лучше приближение. На практике в какой-то момент вы получите числовые ошибки, но пока вы используете числа с плавающей запятой, меньшие значения, чем 1.6, обеспечат вам лучшее приближение.

Чтобы проиллюстрировать это, я построил сечение LoG и DoG для нескольких значений k в Mathematica:

введите описание изображения здесь

Как видите, k = 1.6 не является идеальным приближением. Например, k = 1,1 даст гораздо более близкое приближение.

Но обычно вы хотите рассчитать аппроксимации LoG для диапазона сигм. (В противном случае, зачем вообще прибегать к приближению DoG? Вычисление одного фильтрованного изображения LoG не дороже, чем вычисление одного фильтрованного изображения DoG.) Поэтому значение k обычно выбирается так, что вы можете рассчитать серию гауссовских фильтров изображения с сигмами s, s k, s k ^ 2, s * k ^ 3 ..., а затем вычислить различия между соседними гауссианами. Поэтому, если вы выберете меньшее k, вам придется рассчитать больше «слоев» гауссианов для того же сигма-диапазона. k = 1.6 - это компромисс между желанием приблизительного приближения и нежеланием вычислять слишком много разных гауссиан.

Однако, я не уверен, как две сигмы в Разнице Гауссианов относятся к сигме для лапласиана Гаусса. Является ли меньшая сигма в первом равно сигме второго?

t=σ2σ2+Δtσ2ΔtΔt0

σLaplace=σ1+k22

Ники Эстнер
источник
я прошу прощения , если я не прав, но это не то , что на самом деле вычисления LOG является более дорогим , чем DoG. поскольку гауссов можно разделить на 2 1D фильтра, то есть сложность будет линейной O (2n) вместо полинома O (n ^ 2)
user1916182
@ user1916182: Правда, фильтр LoG сам по себе не отделим. Но ни один из них не является фильтром DoG. Но оба они представляют собой суммы двух разделяемых фильтров (два гауссиана с разным масштабом для DoG, два гауссовых производных фильтра 2-го порядка для LoG). Вы делаете сэкономить время с собачкой , если вы можете использовать «больше» два гауссовые для следующего уровня шкалы, так что вы должны вычислить п + 1 Гаусс , для п чешуи, в отличии от 2 * N гауссовых производных фильтров для п LOG чешуи ,
Ники Эстнер
3

Может быть, формулы здесь могут помочь вам.

Поскольку представление масштабного пространства удовлетворяет уравнению диффузии, LoG может быть вычислено как разность между двумя срезами масштабного пространства.

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

ЛИБОР
источник
Спасибо, но я уже посмотрел на них. Похоже, они не говорят мне, является ли сигма или k * sigma значением, соответствующим параметру t (которое совпадает со значением сигмы для уравнения Лапласа Гаусса).
визуально-кинетическая
1
Это где-то посередине: s <t <k * s. Поскольку разность (y (a) - y (b)) / (ba) приближает (когда b - a -> 0) производную при (a + b) / 2. Однако, поскольку вы не берете предел k-> 1, это только приблизительное значение, и вы не сможете точно определить лучшую сигму (если вы не определите конкретный критерий оптимизации).
nimrodm