A computer-implemented method comprising the steps of: (a) receiving seismic data acquired by a seismic acquisition array comprising a plurality of receivers in a survey of a region; (b) performing a preliminary analysis on selected traces in the seismic data over the surveyed region to identify surface wave properties, the surface wave properties including wave modes, wave types, frequency and velocity ranges of the surface wave for each different wave mode and/or wave type and lateral variations within the selected traces; (c) for each of a plurality of locations associated with the region: (i) based on the preliminary assessment, selecting wave modes to be analyzed and for the selected wave modes, designing data processing geometries, wherein the data processing geometries include an azimuthal range and a source-receiver offset range that correspond to a sub-array of the overall seismic acquisition array at that location, wherein the azimuthal range defines the boundaries of a region surveyed by the sub-array at that location and wherein the offset range includes a frequency dependent minimum and maximum offset that takes near field and attenuation effects into account; (ii) selecting traces from the sub-array that conform to the defined azimuthal and offset ranges for the sub-array at the location and applying a wavefield transform to map data of the selected traces from the time-source offset domain into a frequency wavenumber domain; and (iii) determining a dispersion curve for the selected traces by computing the phase velocity at each frequency and determining an attenuation curve for the location by computing the attenuation coefficient at each frequency; and (d) interpolating the dispersion and attenuation curves to build pseudo- sections of surface wave propagation parameters for each location and integrating the pseudo sections for the plurality of locations to form a spatially continuous representation of the surface wave propagation parameter across the surveyed region. (a) receiving seismic data acquired by a seismic acquisition array comprising a plurality of receivers in a survey of a region; (b) performing a preliminary analysis on selected traces in the seismic data over the surveyed region to identify surface wave properties, the surface wave properties including wave modes, wave types, frequency and velocity ranges of the surface wave for each different wave mode and/or wave type and lateral variations within the selected traces; (c) for each of a plurality of locations associated with the region: (i) based on the preliminary assessment, selecting wave modes to be analyzed and for the selected wave modes, designing data processing geometries, wherein the data processing geometries include an azimuthal range and a source-receiver offset range that correspond to a sub-array of the overall seismic acquisition array at that location, wherein the azimuthal range defines the boundaries of a region surveyed by the sub-array at that location and wherein the offset range includes a frequency dependent minimum and maximum offset that takes near field and attenuation effects into account; (ii) selecting traces from the sub-array that conform to the defined azimuthal and offset ranges for the sub-array at the location and applying a wavefield transform to map data of the selected traces from the time-source offset domain into a frequency wavenumber domain; and (iii) determining a dispersion curve for the selected traces by computing the phase velocity at each frequency and determining an attenuation curve for the location by computing the attenuation coefficient at each frequency; and (i) based on the preliminary assessment, selecting wave modes to be analyzed and for the selected wave modes, designing data processing geometries, wherein the data processing geometries include an azimuthal range and a source-receiver offset range that correspond to a sub-array of the overall seismic acquisition array at that location, wherein the azimuthal range defines the boundaries of a region surveyed by the sub-array at that location and wherein the offset range includes a frequency dependent minimum and maximum offset that takes near field and attenuation effects into account; (ii) selecting traces from the sub-array that conform to the defined azimuthal and offset ranges for the sub-array at the location and applying a wavefield transform to map data of the selected traces from the time-source offset domain into a frequency wavenumber domain; and (iii) determining a dispersion curve for the selected traces by computing the phase velocity at each frequency and determining an attenuation curve for the location by computing the attenuation coefficient at each frequency; and (d) interpolating the dispersion and attenuation curves to build pseudo- sections of surface wave propagation parameters for each location and integrating the pseudo sections for the plurality of locations to form a spatially continuous representation of the surface wave propagation parameter across the surveyed region.