The Simplified Spherical Harmonic (SPN) approximation was first introduced as a three-dimensional (3-D) extension of the plane-geometry Spherical Harmonic (PN) equations. A third order SPN (SP3) solver, recently implemented in the Nodal Expansion Method (NEM), has shown promising performance in the reactor core neutronics simulations. This work is focused on the development and implementation of the transport corrected interface and boundary conditions in NEM SP3 solver, following recent published work on the rigorous SPN theory for piecewise homogeneous regions. A streamlined procedure has been developed to generate the flux zero and second order/moment discontinuity factors (DFs) of the generalized equivalence theory to eliminate the error introduced by pin-wise homogenization. Moreover, several color set models with varying size and configuration are later explored for their capability of generating DFs that can produce results equivalent to that using the whole-core homogenization model for more practical implementations. The new developments are tested and demonstrated on the C5G7 benchmark. The results show that the transport corrected SP3 solver shows general improvements to power distribution prediction compared to the basic SP3 solver with no DFs or only zero order/moment DFs. The complete equivalent calculations using the DFs can almost reproduce transport solutions with high accuracy. The use of equivalent parameters from larger size color set models show better prediction in the whole-core calculations. By coupling different color set models DFs can offer the best accuracy at both eigenvalues and power distributions.