Вы правильно поняли, просто используйте фильтры Гаусса, когда вычисляете пространственные производные.
Наконец, получите области, заполненные градиентом, выбрав темные пиксели величины градиента ориентационного изображения.
EDIT
Статья в Википедии о гауссовых фильтрах .
В моем коде я использую следующие ядра свертки (горизонтальные и вертикальные производные):
{
{0.00623304, 0.012975, 0., -0.012975, -0.00623304},
{0.0259501, 0.054019, 0., -0.054019, -0.0259501},
{0.0581332, 0.121013, 0., -0.121013, -0.0581332},
{0.0259501, 0.054019, 0., -0.054019, -0.0259501},
{0.00623304, 0.012975, 0., -0.012975, -0.00623304}
}
{
{0.00623304, 0.0259501, 0.0581332, 0.0259501, 0.00623304},
{0.012975, 0.054019, 0.121013, 0.054019, 0.012975},
{0., 0., 0., 0., 0.},
{-0.012975, -0.054019, -0.121013, -0.054019, -0.012975},
{-0.00623304, -0.0259501, -0.0581332, -0.0259501, -0.00623304}
}
В качестве доказательства концепции, это может выглядеть в Mathematica, используя тестовое изображение, которое содержит наклонную область градиента и область вертикального градиента: