The present disclosure provides a system for and a method of obtaining a magnetic resonance image by performing magnetic resonance imaging (MRI) at multiple slices simultaneously. The method comprises generating a multiband pulse sequence for spin-echo imaging, the pulse sequence comprising a multiband excitation pulse and at least one multiband refocusing pulse, wherein the multiband excitation pulse simultaneously excites multiple bands, wherein the at least one multiband refocusing pulse simultaneously refocuses the multiple bands, and wherein the phases of the bands excited by the multiband excitation pulse and the phases of the bands refocused by the at least one multiband refocusing pulse are set according to a single row of an orthogonal encoding matrix. The multiband excitation pulse and the at least one multiband refocusing pulse collectively form a multiband pulse pair.