A method of reconstructing an (n+l)-dimensional image functionf representing a region of investigation (1) comprises the step of determining the image function f from n-dimensional or less dimensional Radon data comprising a plurality of projection functions p?(t) measured corresponding to a plurality of predetermined projection directions (T), wherein the image function f is determined as a sum polynomials multiplied with values of the projection functions p?(t). An imaging method and an imaging device, like e.g. a computer tomography device, using this reconstruction method are described.
