A high order method for determining the edges in the gradient of a function

Rishu Saxena, Anne Gelb, Hans Mittelmann

Research output: Contribution to journalArticlepeer-review

3 Scopus citations


Detection of edges in piecewise smooth functions is important in many applications. Higher order reconstruction algorithms in image processing and post processing of numerical solutions to partial differential equations require the identification of smooth domains, creating the need for algorithms that will accurately identify discontinuities in a given function as well as those in its gradient. This work expands the use of the polynomial annihilation edge detector, (Archibald, Gelb and Yoon, 2005), to locate discontinuities in the gradient given irregularly sampled point values of a continuous function. The idea is to preprocess the given data by calculating the derivative, and then to use the polynomial annihilation edge detector to locate the jumps in the derivative. We compare our results to other recently developed methods.

Original languageEnglish (US)
Pages (from-to)694-711
Number of pages18
JournalCommunications in Computational Physics
Issue number2-4
StatePublished - Feb 2009


  • Derivative discontinuities
  • Multivariate edge detection
  • Non-uniform grids
  • Piecewise smooth functions
  • Polynomial annihilation

ASJC Scopus subject areas

  • Physics and Astronomy (miscellaneous)


Dive into the research topics of 'A high order method for determining the edges in the gradient of a function'. Together they form a unique fingerprint.

Cite this