This Article 
   
 Share 
   
 Bibliographic References 
   
 Add to: 
 
Digg
Furl
Spurl
Blink
Simpy
Google
Del.icio.us
Y!MyWeb
 
 Search 
   
High-Order Differentiation Filters that Work
July 1994 (vol. 16 no. 7)
pp. 734-739

Reliable derivatives of digital images have always been hard to obtain, especially (but not only) at high orders. We analyze the sources of errors in traditional filters, such as derivatives of the Gaussian, that are used for differentiation. We then study a class of filters which is much more suitable for our purpose, namely filters that preserve polynomials up to a given order. We show that the errors in differentiation can be corrected using these filters. We derive a condition for the validity domain of these filters, involving some characteristics of the filter and of the shape. Our experiments show a very good performance for smooth functions.

[1] P. J. Besl, J. B. Birch, and L. T. Watson, "Robust window operations," inProc. IEEE 2nd Int. Conf. Comp. vision, Dec. 1988, pp. 591-600.
[2] A. Blake and A. Zisserman,Visual Reconstruction. Cambridge, MA: MIT Press, 1987.
[3] P. Dierckx, "An algorithm for least square fitting of cubic spline surfaces to functions on a rectilinear mesh over a rectangle,"J. Comput. Appl. Math., vol. 3, p. 2, 1977.
[4] W. E. L. Grimson, "The implicit constraint of the primal sketch," MIT, A.I. Memo No. 663, 1981.
[5] R.M. Haralick, "Digital step edges from zero crossings of second directional derivatives",IEEE Trans. Pattern Anal. Machine Intell., vol. 6, pp. 58-68, 1984.
[6] M. Hashimoto and J. Sklansky, "Multiple-order derivatives for detecting local image characteristics,"Comput. Vision, Graphics, and Image Processing, vol. 39, pp. 28-55, 1987.
[7] B. K. P. Horn, "The curve of least energy,"ACM Trans. Math. Software, vol. 9, pp. 441-460, 1983.
[8] M. H. Hueckel, "A local visual operator which recognizes edges and lines,"J. Assoc. Comput. Machine, vol. 20, pp. 634-647, 1973.
[9] R. A. Hummel, "Feature detection using basis functions,"Comput. Graphics and Image Processing, vol. 9, pp. 40-55, 1979.
[10] M. Kass, A. Witkin, and D. Terzopoulos, "Snakes: Active contour models",Int. J. Computer Vision, vol. 1, pp. 321-331, 1988.
[11] P. Meer and I. Weiss, "Smoothed differentiation filters for images,"J. Visual Commun. and Image Representat., vol. 3, pp. 58-72, 1992.
[12] J. M. Ortega and W. G. Poole,An Introduction to Numerical Analysis for Differential Equations. Marshfield, MA: Pitman, 1981.
[13] T. Poggio,Proc. DARPA Image Understanding Workshop, Los Angeles, CA, 1987.
[14] T. Poggio, V. Torre, and C. Koch, "Computational vision and regularization theory,"Nature, vol. 317, pp. 314-319, 1985.
[15] T. Poggio, H. Voorhees, and A. Yuille, "A regularized solution to edge detection,"J. Complexity, vol. 4, pp. 106-123, 1988.
[16] G. Suria and M. Subbarao, "Depth from focus by changing camera aperture: a spatial domain approach," inProc. IEEE Comput. Vision Pattern Recognit., NY, June 1993, pp. 61-67.
[17] D. Terzopoulos, "The computation of visible surface representation",IEEE Trans. Pattern Anal. Machine. Intell., vol. PAMI-10, pp. 417-438, 1988.
[18] I. Weiss, "Projective invariants of shapes," inProc. DARPA Image Understanding Workshop, Cambridge, MA, 1988.
[19] I. Weiss, "Higher order differentiation filters that work," Univ. of Maryland CS-TR-2639, 1991.
[20] I. Weiss, "Geometric invariants and object recognition,"Int. J. Comput. Vision, vol. 10, no. 3, pp. 201-231, 1993.
[21] I. Weiss, "Noise resistant invariants of curves,"IEEE Pattern Anal. Machine Intell., vol. 15, pp. 943-948, Sept. 1993.

Index Terms:
image processing; filtering and prediction theory; polynomials; high-order differentiation filters; digital images; error source; polynomials; validity domain; smooth functions; image smoothing; noise suppression
Citation:
I. Weiss, "High-Order Differentiation Filters that Work," IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 16, no. 7, pp. 734-739, July 1994, doi:10.1109/34.297955
Usage of this product signifies your acceptance of the Terms of Use.