A novel reconstruction method was proposed for dual-band EPI in an MRI system equipped with four-channel receiver coils. This method was based on a conventional kernel method utilizing an iterative calculation with regional constraints in the image domain. The method significantly improves the quality of the reconstructed images, even in the regions with less coil sensitivity. The results showed higher signal-to-noise ratio, less signal leakage, and better long-term stability in repetitions in comparison to the conventional method. The proposed method can be applied to clinical systems that have relatively few receiver coils, as well as animal systems.