In this paper, we develop a homogenization methodology for the two-dimensional low-order quasi-diffusion equations for full-core reactor calculations that is based on a family of spatially consistent coarse-mesh discretization methods. The coarse-mesh solution generated by these methods preserves a number of spatial moments of the fine-mesh transport solution over each assembly. The proposed method reproduces accurately the complicated large-scale behavior of the transport solution within assemblies. To demonstrate the performance of the developed methodology, we present the numerical results of several test problems that simulate mixed-oxide-uranium and assembly-reflector interfacial effects.