There is provided a computer-implemented method for solving inverse imaging problems to compensate for distortions in an image. The method comprises: minimizing a cost objective function containing a data fitting term and one or more image prior terms to each of the plurality of channels, the one or more image prior terms comprising cross-channel information for a plurality of channels derived from the image.