The three dimensional (3D) innervation zone (IZ) imaging approach (3DIZI) has been developed in our group to localize the IZ of a particular motor unit (MU) from its motor unit action potentials decomposed from high-density surface electromyography (EMG) recordings. In this study, the developed 3DIZI approach was combined with electrical stimulation to investigate global distributions of IZs in muscles from electrically elicited M-wave recordings. Electrical stimulations were applied to the musculocutaneous nerve to activate supramaximal muscle response of the biceps brachii in one healthy subject, and high-density (128 channels) surface EMG signals of the biceps brachii muscles were recorded. The 3DIZI approach was then employed to image the IZ distribution of IZs in the 3D space of the biceps brachii. The performance of the M-wave based 3DIZI approach was evaluated with different stimulation intensities. Results show that the reconstructed IZs under supramaximal stimulation are spatially distributed in the center region of muscle belly which is consistent with previous studies. With sub-maximal stimulation intensity, the imaged IZ centers became more proximally and deeply located. The proposed M-wave based 3DIZI approach demonstrated its capability of imaging global distribution of IZs in muscles, which provide valuable information for clinical applications such as guiding botulinum toxin injection in treating muscle spasticity.