A simple analytical calculation scheme to determine near field radiation through decomposing an emission domain into lots of thin thermal current sheets is presented. Through finding the orthogonal modes of thermal current of each thin layer, the thin current sheets can be treated as radiation sources of electromagnetic waves with determined analytical solutions. The outgoing electromagnetic waves from each thin current sheets can be either in transverse electric (TE) or transverse magnetic (TM) modes depending on the orientations of the current in the thin current sheets with respect to the directions of amplitude modulations of the orthogonal modes. Electromagnetic waves arriving to a collection domain are related to the electromagnetic waves leaving from each thin current thermal sheet with transfer coefficients. Transfer coefficient for TE and TM waves can be determined analytically with transfer matrix method or scattering matrix methods. Compared with existing dyadic Green's function method, the new calculation scheme allows material and temperature variations along one direction of the emission domain based on determined analytical expressions of TE and TM waves leaving from each thin current sheets. The simple calculation scheme is especially useful in near field radiation of layered structures with different material such as hyperbolic material with negative refractive indices. With this new approach, we recovered analytical solutions of near field radiation between two semi-infinite domains with uniform temperature and derived closed form solution of near field radiation between two semi-infinite domains with temperature profiles with/without laminated structures.