分数阶微分,Grumwald-Letnikov定义在图像的数值实现中更为准确:
Gamma函数:
若s(t)的持续期t [a,t],将函数持续期间[a,t]按单位间隔h=1进行等分,得到:
推到一元函数是s(t)的 v 阶分数阶微分差分表达式为:
这n个非零系数只有常数“1”,其他都是n-1个都是分数阶微分阶次的函数。n个系数分别为:
微分近似掩模算子
分数阶微分( 0.1阶 ~ 0.9阶), 则可以检测到或部分检测到那些模糊边界和细微的弱边界, 而且还可以减少图像噪声的增加。
下左图:原图,下右图:0.5阶分数阶微分处理结果(v=0.5)
void CShowPicView::OnFractionalCalculus( ) // 分数阶微积分 灰度图 单通道
{
CShowPicDoc* pDoc = GetDocument();
CDC* pDC=GetDC();
CString str;
int i,j;
double temp ;
int **g;