Abstract
We present a method that uses Fourier spectral data to locate jump discontinuities in the first derivatives of functions that are continuous with piecewise smooth derivatives. Since Fourier spectral methods yield strong oscillations near jump discontinuities, it is often difficult to distinguish true discontinuities from artificial oscillations. In this paper we show that by incorporating a local difference method into the global derivative jump function approximation, we can reduce oscillations near the derivative jump discontinuities without losing the ability to locate them. We also present an algorithm that successfully locates both simple and derivative jump discontinuities.
Similar content being viewed by others
References
Archibald, R., Chen, K., Gelb, A., Renaut, R.: Improving tissue segmentation of human brain MRI through pre-processing by the Gegenbauer reconstruction method. NeuroImage 20(1), 489–502 (2003)
Archibald, R., Gelb, A.: A method to reduce the gibbs ringing artifact in MRI scans while keeping tissue boundary integrity. IEEE Med. Imag. 21(4), 305–319 (2002)
Archibald, R., Gelb, A., Yoon, J.: Polynomial fitting for edge detection in irregularly sampled signals and images. SINUM 43, 259–279 (2005)
Banergee, N., Geer, J.: Exponentially accurate approximations to piecewise smooth periodic Lipschitz functions based on Fourier series partial sums. J. Sci. Comput. 13, 419–460 (1998)
Bary, N.: Treatise of Trigonometric Series. Macmillan, New York (1964)
Bauer, R.B.: Band filters for determining shock locations. Ph.D. thesis, Applied Mathematics, Brown University, Providence, RI (1995)
Boyd, J.P.: Chebyshev and Fourier spectral methods (2nd edn.). Dover Publications, Mineola, New York (2001)
Brezinski, C.: Extrapolation algorithms for filtering series of functions and treating the Gibbs phenomenon. Numer. Algorithms 36, 309–329 (2004)
Canny, J.: A computational approach to edge detection. IEEE Trans. PAMI, 8(6), 679–698 (1986)
Canuto, C., Hussaini, M.Y., Quarteroni, A., Zang, T.A.: Spectral Methods in Fluid Dynamics. Springer, Berlin (1988)
Curlander, J., McDonough, R.: Synthetic Aperture Radar Systems and Signal Processing. Wiley, New York (1991)
Driscoll, T.A., Fornberg, B.: A Padé-based algorithm for overcoming the gibbs phenomenon. Numer. Algorithms 26, 77–92 (2001)
Eckhoff, K.S.: On a high order numerical method for functions with singularities. Math. Comput. 67, 1063–1087 (1998)
Fornberg, B.: A Practical Guide to Pseudospectral Methods. Cambridge University Press, New York (1996)
Gacougnolle, J., Guilpin, C., Simon, Y.: The ε-algorithm allows to detect dirac delta function. Appl. Numer. Math. 48, 27–40 (2004)
Gelb, A., Tadmor, E.: Detection of edges in spectral data. Appl. Comput. Harmon. Anal. 7, 101–135 (1999)
Gelb, A., Tadmor, E.: Detection of edges in spectral data II – nonlinear enhancement. SINUM 38:4, 1389–1408 (2000)
Gelb, A., Tadmor, E.: Enhanced spectral viscosity approximation for conservation laws. Appl. Numer. Math. 33, 1–21 (2000)
Gelb, A., Tadmor, E.: Adaptive edge detectors for piecewise smooth data based on the minmod limiter. J. Sci. Comput. 28(2–3), 279–306 (2006)
Gelb, A., Tanner, J.: Robust reprojection methods for the resolution of the Gibbs phenomenon. ACHA 20(1), 3–25 (2006)
Hesthaven, J., Gottlieb, S., Gottlieb, D.: Spectral Methods for Time Dependent Problems. Cambridge University Press (2007)
Gottlieb, D., Orszag, S.: Numerical Analysis of Spectral Methods: Theory and Applications. SIAM, Philadelphia (1977)
Gottlieb, D., Shu, C.W.: On the gibbs phenomenon and its resolution. SIAM Rev. 30, 644–668 (1997)
Gottlieb, D., Tadmor, E.: Recovering pointwise values of discontinuous data within spectral accuracy, progress and supercomputing in computational fluid dynamics. In: Murman, E.M., Abarbanel, S.S. (eds.) Proceedings of a 1984 U.S.–Israel Workshop, Progress in Scientific Computing, vol. 6 (1985)
Hildreth, E.C., Marr, D.: Theory of edge detection. Roy. Soc. Lond. Ser. B. 207, 187–217 (1980)
Hwang, W., Mallat, S.: Singularity detection and processing with wavelets. IEEE Trans. Inf. Theory 38, 617–643 (1992)
Kvernadze, G.: Determination of the jump of a bounded function by its Fourier series. Approx, J. Theory 92, 167–190 (1998)
Lanczos, C.: Applied Analysis. Prentice Hall, Englewood Cliffs, New Jersey (1961)
Liang, Z., Lauterbur, P.: Principles of Magnetic Resonance Imaging: A Signal Processing Perspective. The Institute of Electrical and Electronics Engineers Press, New York (2000)
Oliver, C., Quegan, S.: Understanding Synthetic Aperture Radar Images. Artech House, Boston (1998)
Shizgal, B., Jung, J.H.: Towards the resolution of the gibbs phenomena. Comp. J. Appl. Math. 161, 41–65 (2003)
Tadmor, E., Tanner, J.: Adaptive mollifiers – high resolution recovery of piecewise smooth data from its spectral information. J. Foundations of Comp. Math. 2, 155–189 (2002)
Tanner, J.: Optimal filter and mollifier for piecewise smooth spectral data. Math. Comput. 75, 767–790 (2005)
Author information
Authors and Affiliations
Corresponding author
Additional information
This work was partially supported by NSF grants CNS 0324957 and DMS 0510813, and NIH grant EB 02553301 (AG).
Rights and permissions
About this article
Cite this article
Cates, D., Gelb, A. Detecting derivative discontinuity locations in piecewise continuous functions from Fourier spectral data. Numer Algor 46, 59–84 (2007). https://doi.org/10.1007/s11075-007-9127-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-007-9127-x