A fast wavelet algorithm for multidimensional signals using polyharmonic splines