A method to select a representative subset of a plurality of horizon surfaces or surface patches from geophysical subsurface imaging data, including: defining a score function on one or more horizon surfaces or surface patches; calculating, by a computer, the score for each of the plurality of horizon surfaces or surface patches with regard to other horizon surfaces or surface patches and whether the other horizon surfaces or surface patches have been selected or not for inclusion or exclusion in the subset of the plurality of horizon surfaces; selecting, by a computer, one or more of the plurality of horizon surfaces or surface patches to be included in the subset of the plurality of horizon surfaces or surface patches or excluded from the subset of the plurality of horizon surfaces or surface patches based on their respective scores; iteratively repeating the selecting and calculating steps until a stopping condition is reached and the subset of the plurality of horizon surfaces or surface patches is determined; and performing interpretation on the subset of the plurality of horizon surfaces or surfaces patches.