Abstract: A method for generating a geophysical image of a subsurface region includes defining a computational sub-volume for the geophysical image including a predetermined number of seismic traces of a plurality of seismic traces and a predetermined number of samples per each one of the plurality of seismic traces, generating a data matrix corresponding to a first sub-volume of the subsurface region based on the defined computational sub-volume, the data matrix comprising the predetermined number of samples for the predetermined number of traces of a portion of a seismic dataset corresponding to the first sub-volume. The method also includes performing a singular value decomposition and estimating a coherence between the predetermined number of traces of the data matrix by performing a weighted sum of a variance of a plurality of right singular vectors of the data matrix, and assigning the estimated coherence to a location in the geophysical image.