For understanding physical mechanisms occurring in ferromagnetic nano-layers, micromagnetic simulations are commonly used . In order to perform an advanced analysis of complex magnetization dynamics induced by phenomena such as Spin-Transfer-Torque (STT) [1,2] or voltage-controlled anisotropy [3,4] changes, a spatial distribution of power spectral density is often needed [5] . We use a specially designed open-source tool for creating spectral density maps directly from a set of subsequent magnetization files [6], an output typically available in micromagnetic simulations [7] . By investigating multiple frequencies and areas, oscillation modes and nodes of spin waves can be precisely localized . In this work we discuss a ferromagnetic resonance in Magnetic Tunnel Junctions (MTJs) . The Fast Fourier Transform (FFT) can be computed once for the whole oscillator, which produces a single distribution consistent with experimental results . However, in such case the information about the localization of magnetization oscillations modes of is lost . Our approach allows for efficient creating of spectral density maps with frequency domain analysis conducted separately for each cell . Because of the significant size of the data involved, we implemented parallel computing methods to greatly reduce computations time . We examine the distribution of FFT amplitude for the respective frequencies of identified peaks in both free and reference layer of the MTJ-based nanooscillator . An example result calculated for an STT-induced resonance is presented in Fig . 1 . The differences in localizations of two modes are clearly visible, allowing for more detailed analysis of the junction behavior and peaks origin . Using such an approach we compare localization of frequency modes excited by: single pulses of magnetic field, alternating magnetic field, constant spin-polarized current flow, alternating current (via spin-diode effect) and alternating voltage which induces anisotropy changes . Obtained results will be presented and used in attempt to optimize excitation types as well as geometrical and material parameters of experimentally investigated MTJs . We also discuss opportunities for further analysis based on our approach combined with wavelet transform in case of non-stationary dynamic processes in MTJs . 7] M . J . Donahue, D . G . Porter, NIST report (1999) . Figure 1. Spectral analysis of a 250×150 nm spintorque oscillator based on MTJ under constant magnetic field. a) overall FFT spectrum b) FFT spectrum density for free layer, f=4.7 GHz. c) FFT spectrum density for free layer, f=6.2 GHz. d) FFT spectrum density for reference layer, f=4.7 GHz. e) FFT spectrum density for reference layer, f=6.2 GHz.