More and more medical devices can capture different features of human body and form three dimensional (3D) images. In clinical applications, usually it is required to fuse multiple source images containing different and crucial information into one for the purpose of assisting medical treatment. However, traditional image fusion methods are normally designed for two dimensional (2D) images and will lead to loss of the third dimensional information if directly applied to 3D data. Therefore, a novel 3D magnetic image fusion method was proposed based on the combination of newly invented beyond wavelet transform, called 3D band limited shearlet transformand (BLST), and four groups of traditional fusion rules. The proposed method was then compared with the 2D and 3D wavelet and dual-tree complex wavelet transform fusion methods through 4 groups of human brain T2* and quantitative susceptibility mapping (QSM) images. The experiments indicated that the performance of the method based on 3D transform was generally superior to the existing methods based on 2D transform. Taking advantage of direction representation, shearlet transform could effectively improve the performance of conventional fusion method based on 3D transform. It is well concluded, therefore, that the proposed method is the best among the methods based on 2D and 3D transforms.