A method for stratigraphic analysis of seismic data, known as GeoLeaf, comprises: selecting a seismic data volume comprising a plurality of samples; selecting a number of horizons (TH, BH1 TPH, BPH, SMH) within the seismic data volume; selecting a scale of interest within the seismic data volume, which ranges from 500 m up to the size of the entire seismic data volume; applying an attribute function to at least one of the selected number of horizons (TH, BH, TPH, BPH, SMH) at the selected scale of interest; and displaying a stratigraphic attribute, which represents a geologic feature, such as a seismic fades attribute (FM), an unconformity attribute (Zmax/Z) and/or a thinning attribute (Z1/Z2).