一维灰度值一阶导二阶导对应卷积掩码推导

一维边缘提取过程中，计算离散一维灰度值剖面的导数：

对应两个卷积掩码为：

1/2(1 0 -1) 和1/2(1 -2 1)

推导过程：

三点拟合抛物线，抛物线公式：y = ax^2+bx+c,表示为矩阵方式为：

其中e(x)为误差。

 

在1X3的领域内，中心点为(0)的3个点[-1 0 1]，,有：

即：f+e = Ap，用最小二乘法求最优解

，叫做广义逆矩阵，只与A相关，计算得：

对F(x) = ax^2+bx+c 求导 f’(x) = 2ax+b。

x = 0, 故：f^’(x)  =  b  =-0.5*f_{i-1 +} 0*f_i+0.5* f_i+1得：

F_i^’=1/2(f_i+1– f_i-1)

对F(x) = ax^2+bx+c 求二阶导得：f^’’(x) = 2a，故

F_i^’’=0.5*f_i+1 –1*f_i+ 0.5*f_i-1 得：

F_i^’’=1/2(f_i+1–2f_i + f_i-1)

证毕。

 

个人觉得比较难理解的是x = 0，求导主要是等间隔(x)条件下,f(x)的变换率，所以x的值都可以+/-偏移量化解为[-1 0 1]的形式，不影响结果。