In the present paper, a method was proposed to represent the twophase, homogeneous random field. In the first place, the framework of twophase random field was put forward and the homogeneous and isotropic assumption was given as the premise of the work. Then, the stochastic harmonic function (SHF) representation for the Gaussian random field was provided as the intermedia for the twophase random field. With the SHF based random field, the Nataf’s transformation was applied to convert Gaussian random field into the target random field. The Hermite polynomials were used as the numerical solution for the bivariate Gaussian pdf involved in the Nataf’s transformation. Finally, several numerical examples were given to clarify the validity of stochastic harmonic function based twophase random field.