We propose a new  generic and flexible methodology for nonparametric function estimation  in which we first estimate the number and locations of any features that may be present in the function  and then estimate the function parametrically between each pair of neighbouring detected features