A computer-implemented method of seismic imaging, comprising: a. receiving (10, 11), at one or more computer processors, a field seismic dataset representative of a subsurface volume of interest and a synthetic seismic dataset wherein the synthetic seismic dataset is generated with no or known attribute variation with angle or other spatial coordinates, wherein the synthetic seismic dataset includes synthetic seismic data generated by forward modeling through an earth model having layer boundaries; b. applying (12, 13), via the one or more computer processors, seismic processing to the field seismic dataset and the synthetic seismic dataset to generate a processed field seismic dataset and a processed synthetic seismic dataset wherein the processed field seismic dataset and the processed synthetic seismic dataset each comprise output images and contain attribute distortion due to the seismic processing, wherein the processed synthetic seismic dataset comprises synthetic image responses of individual layer boundaries of the earth model computed individually, free from interference by image artifacts of other layer boundaries; c. extracting, at the one or more computer processors, the position, phase and amplitude from the synthetic image response of each of the individual layer boundaries; and d. correcting, via the one or more computer processors, the attribute distortion in the field seismic dataset based on the distortion assessed from the processed synthetic seismic dataset using deviations in the extracted position, phase and amplitude in the synthetic image responses and the expected position, phase and amplitude from the earth model to generate a corrected processed field dataset, and wherein: the field seismic dataset includes a plurality of traces recorded at a plurality of seismic sensors for each of a plurality of seismic shots. a. receiving (10, 11), at one or more computer processors, a field seismic dataset representative of a subsurface volume of interest and a synthetic seismic dataset wherein the synthetic seismic dataset is generated with no or known attribute variation with angle or other spatial coordinates, wherein the synthetic seismic dataset includes synthetic seismic data generated by forward modeling through an earth model having layer boundaries; b. applying (12, 13), via the one or more computer processors, seismic processing to the field seismic dataset and the synthetic seismic dataset to generate a processed field seismic dataset and a processed synthetic seismic dataset wherein the processed field seismic dataset and the processed synthetic seismic dataset each comprise output images and contain attribute distortion due to the seismic processing, wherein the processed synthetic seismic dataset comprises synthetic image responses of individual layer boundaries of the earth model computed individually, free from interference by image artifacts of other layer boundaries; c. extracting, at the one or more computer processors, the position, phase and amplitude from the synthetic image response of each of the individual layer boundaries; and d. correcting, via the one or more computer processors, the attribute distortion in the field seismic dataset based on the distortion assessed from the processed synthetic seismic dataset using deviations in the extracted position, phase and amplitude in the synthetic image responses and the expected position, phase and amplitude from the earth model to generate a corrected processed field dataset,