这是IEEE Image Processing, vol.9, No.6, 2000.上一篇文章的代码