Abstract This article presents a two-step procedure for the computation of radiative heat transfer with anisotropic scattering and reflection. It is based on a concept that the coincident processes of absorption and scattering/reflection can be separated factitiously. All medium elements and wall surfaces are supposed to be pure-absorbing when receiving incident radiation. Afterwards they emit the scattered/reflected radiations. The absorption of both the initial and the secondary radiations can be assessed by the direct exchange area. It is needed to repeat the processes for a few times until the radiations are substantially absorbed. For anisotropic scattering/reflection, a vector summation obtains the directional distribution of emissive power. The method is validated by several benchmark computations in terms of emissive power and heat transfer coefficients. It is shown that the method gives more accurate solution than the isotropic scaling for the heat transfer in anisotropically scattering media.