One possible approach is to perform a local plane-wave decomposition, integrating the component of the vector-wavefield in the local polarization direction () along different directions p in order to build a function A(p,x,z,t), which represents the phase amplitude for each point of the space-time as a function of the horizontal slowness, and a polarization vector function . Time averaging should in this case provide three expectations:, , and . Next we need to use the Christoffel equation with the background elastic parameters used in the migration to find the wave type associated with each polarization vector. Then, an anisotropic inversion scheme can estimate the elastic parameters that best fit the scattering relation involving for the different modes and , as a function of the horizontal slowness .