Pixel images {circumflex over (f)} are created from projections (q1 . . . qp) by backprojecting selected projections to produce intermediate images (Il,m), and performing digital image coordinate transformations and/or resampling on selected intermediate images. The digital image coordinate transformations are chosen to account for view angles of the constituent projections of the intermediate images and for their Fourier characteristics, so that the intermediate images may be accurately represented by sparse samples. The resulting intermediate images are aggregated into subsets, and this process is repeated in a recursive manner until sufficient projections and intermediate images have been processed and aggregated to form the pixel image {circumflex over (f)}. Digital image coordinate transformation can include rotation, shearing, stretching, contractions. Resampling can include up-sampling, down-sampling, and the like. Projections can be created from a pixel image ({circumflex over (f)}), by performing digital image coordinate transformation and/or resampling and/or decimation re-projecting the final intermediate image.