A method for processing seismic data. The method includes performing a plurality of stochastic simulations for one or more rock model parameters to generate one or more anisotropic parameters for a subsurface area of the earth. The method then derives one or more joint multi-dimensional probability density functions for the anisotropic parameters. Using the joint multi-dimensional probability density functions and measured well log data, the method computes one or more posterior probability density functions. The method then includes deriving one or more anisotropic profiles from the posterior probability density functions and generating a seismic image from the anisotropic profiles.