AbstractIn this paper we address the problem of approximating functions with discontinuities via kernel-based methods. The main result is the construction of discontinuous kernel-based basis functions. The linear spaces spanned by these discontinuous kernels lead to a very flexible tool which sensibly or completely reduces the well-known Gibbs phenomenon in reconstructing functions with jumps. For the new basis we provide error bounds and numerical results that support our claims. The method is also effectively tested for approximating satellite images.