Abstract: Generating a multispectral or hyperspectral image of an image source with an optical system having an adjustable, wavenumber-dependent point spread function, by collecting panchromatic images of the image source, each of which corresponds to a selected point spread function and includes a measured intensity data set corresponding to a range of wavelengths, transforming the panchromatic images into the spatial frequency domain by using a Fourier transform, solving a matrix equation at each spatial frequency, in which a vector of the transformed panchromatic images is equal to the product of a predetermined matrix of discrete weighting coefficients and a vector representing a wavenumber content of the image source at each spatial frequency, resulting in a determined wavenumber content of the image source in the spatial frequency domain, and inverse transforming the determined wavenumber content of the image source from the spatial frequency domain into the image domain, resulting in the multispectral or hypersp