By stratifying a random scattering media into thin multi-layers in z direction,the first-order scattering solution of each thin layer with Fourier transform is employed to derive high-order scattering solution of the whole random media.An iterative approach to solve high-order scattering solution of a vector radiative transfer (VRT) equation is newly developed.Numerical results are well compared with the Mueller matrix solutions of the first-order for a layer medium,second-order for a half-space,and also the results of the discrete ordinates and eigen-analysis method.It demonstrates that our approach is as feasible,effective and especially applicable to high-order solutions of VRT for both bistatic scattering and thermal emission of inhomogeneous non-spherical scattering media.