A method is described for seismic imaging including receiving a digital seismic image representative of a subsurface volume of interest and an earth model; decomposing the seismic image into a set of image subsets based on an image-domain decomposition function; forward modeling each image subset to generate a set of synthetic data subsets; decomposing each data subset into a set of data sub-subsets based on a data-domain decomposition function; imaging data sub-subsets using the earth model to generate a set of image sub-subsets; and combining the set of image subsets based on a criterion to create a high resolution seismic image. The method may be executed by a computer system.