An automated method for recovering piecewise smooth functions on spheres free from Gibbs oscillations

Chris Blakely, Anne Gelb, Antonio Navarra

Research output: Contribution to journalArticle

2 Citations (Scopus)

Abstract

Spectral methods using spherical harmonic coefficients are widely used in applications for geophysics and atmospheric sciences. The major draw-back in using spherical harmonic spectral methods occurs when the under-lying function is piecewise smooth. In this case, the well-known Gibbs phenomenon reduces the order of accuracy to first order and produces spurious oscillations, particularly in regions near the discontinuities. The frequently studied Gegenbauer reconstruction method has been shown to alleviate the effects of the Gibbs phenomenon while restoring the exponential accuracy of the spectral approximation. Since each reconstruction must be implemented only within smooth regions, the jump discontinuities of the piece-wise smooth function must first be located by an edge detection method. This study combines the recent developments in both edge detection and Gegenbauer reconstruction methods to construct an automated procedure for recovering horizontal or two-dimensional geophysical global fields free from Gibbs oscillations and without compromising the high order convergence properties of spectral methods. Numerical efliciency and robustness are discussed.

Original languageEnglish (US)
Pages (from-to)323-346
Number of pages24
JournalSampling Theory in Signal and Image Processing
Volume6
Issue number3
StatePublished - Sep 2007

Fingerprint

Piecewise Smooth Functions
Edge detection
Spectral Methods
Gibbs Phenomenon
Spherical Harmonics
Edge Detection
Oscillation
Geophysics
Discontinuity
Spectral Approximation
Smooth function
Convergence Properties
Jump
Horizontal
Higher Order
Robustness
First-order
Coefficient

Keywords

  • Edge detection
  • Gegenbauer polynomials
  • Gibbs phenomenon
  • Spherical harmonics

ASJC Scopus subject areas

  • Algebra and Number Theory
  • Analysis
  • Computational Mathematics
  • Radiology Nuclear Medicine and imaging

Cite this

An automated method for recovering piecewise smooth functions on spheres free from Gibbs oscillations. / Blakely, Chris; Gelb, Anne; Navarra, Antonio.

In: Sampling Theory in Signal and Image Processing, Vol. 6, No. 3, 09.2007, p. 323-346.

Research output: Contribution to journalArticle

@article{7477461acd8c46fdb0429598df7edb11,
title = "An automated method for recovering piecewise smooth functions on spheres free from Gibbs oscillations",
abstract = "Spectral methods using spherical harmonic coefficients are widely used in applications for geophysics and atmospheric sciences. The major draw-back in using spherical harmonic spectral methods occurs when the under-lying function is piecewise smooth. In this case, the well-known Gibbs phenomenon reduces the order of accuracy to first order and produces spurious oscillations, particularly in regions near the discontinuities. The frequently studied Gegenbauer reconstruction method has been shown to alleviate the effects of the Gibbs phenomenon while restoring the exponential accuracy of the spectral approximation. Since each reconstruction must be implemented only within smooth regions, the jump discontinuities of the piece-wise smooth function must first be located by an edge detection method. This study combines the recent developments in both edge detection and Gegenbauer reconstruction methods to construct an automated procedure for recovering horizontal or two-dimensional geophysical global fields free from Gibbs oscillations and without compromising the high order convergence properties of spectral methods. Numerical efliciency and robustness are discussed.",
keywords = "Edge detection, Gegenbauer polynomials, Gibbs phenomenon, Spherical harmonics",
author = "Chris Blakely and Anne Gelb and Antonio Navarra",
year = "2007",
month = "9",
language = "English (US)",
volume = "6",
pages = "323--346",
journal = "Sampling Theory in Signal and Image Processing",
issn = "1530-6429",
publisher = "Sampling Publishing",
number = "3",

}

TY - JOUR

T1 - An automated method for recovering piecewise smooth functions on spheres free from Gibbs oscillations

AU - Blakely, Chris

AU - Gelb, Anne

AU - Navarra, Antonio

PY - 2007/9

Y1 - 2007/9

N2 - Spectral methods using spherical harmonic coefficients are widely used in applications for geophysics and atmospheric sciences. The major draw-back in using spherical harmonic spectral methods occurs when the under-lying function is piecewise smooth. In this case, the well-known Gibbs phenomenon reduces the order of accuracy to first order and produces spurious oscillations, particularly in regions near the discontinuities. The frequently studied Gegenbauer reconstruction method has been shown to alleviate the effects of the Gibbs phenomenon while restoring the exponential accuracy of the spectral approximation. Since each reconstruction must be implemented only within smooth regions, the jump discontinuities of the piece-wise smooth function must first be located by an edge detection method. This study combines the recent developments in both edge detection and Gegenbauer reconstruction methods to construct an automated procedure for recovering horizontal or two-dimensional geophysical global fields free from Gibbs oscillations and without compromising the high order convergence properties of spectral methods. Numerical efliciency and robustness are discussed.

AB - Spectral methods using spherical harmonic coefficients are widely used in applications for geophysics and atmospheric sciences. The major draw-back in using spherical harmonic spectral methods occurs when the under-lying function is piecewise smooth. In this case, the well-known Gibbs phenomenon reduces the order of accuracy to first order and produces spurious oscillations, particularly in regions near the discontinuities. The frequently studied Gegenbauer reconstruction method has been shown to alleviate the effects of the Gibbs phenomenon while restoring the exponential accuracy of the spectral approximation. Since each reconstruction must be implemented only within smooth regions, the jump discontinuities of the piece-wise smooth function must first be located by an edge detection method. This study combines the recent developments in both edge detection and Gegenbauer reconstruction methods to construct an automated procedure for recovering horizontal or two-dimensional geophysical global fields free from Gibbs oscillations and without compromising the high order convergence properties of spectral methods. Numerical efliciency and robustness are discussed.

KW - Edge detection

KW - Gegenbauer polynomials

KW - Gibbs phenomenon

KW - Spherical harmonics

UR - http://www.scopus.com/inward/record.url?scp=35348885564&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=35348885564&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:35348885564

VL - 6

SP - 323

EP - 346

JO - Sampling Theory in Signal and Image Processing

JF - Sampling Theory in Signal and Image Processing

SN - 1530-6429

IS - 3

ER -