В википедии же прекрасно расписано: оператор Собеля — обычный градиент с весовыми коэффициентами. Получается сверткой двух матриц: градиента в нужном направлении (скажем, [-1/2 0 1/2] — Dx_{x,y} = (M_{x+1,y} - M_{x-1,y})/2) и весовых коэффициентов (скажем, [1; 2; 1]/4 = [1/4;1/2;1/4]) → [-1/8 0 1/8; -1/4 0 1/4; -1/8 0 1/8]. Естественно, строгого градиента здесь нет — вы вольны выбирать какие угодно весовые коэффициенты.
Можно и расширить матрицу градиента. Если мы возьмем ядро 5x1, получим: Dx_{x,y} = [(M_{x+1,y} - M_{x-1,y})/2) + (M_{x+2,y} - M_{x-2,y})/4)]/2 → [-1/8 -1/4 0 1/4 1/8]. Весами можем взять, например, матрицу [1;2;4;2;1].
Аналогично для матриц 7x7, 9x9 и больше.