Method of image processing of magnetic resonance images for creating de-noised MR images, comprising the steps of - providing image data sets including multiple complex MR images, - subjecting the MR images to a wavelet decomposition for creating coefficient data sets of wavelet coefficients ( S n,m) representing the MR images in a wavelet frequency domain, - calculating normalized coefficient data sets of wavelet coefficients S n, m norm, wherein the coefficient data sets are normalized with a quantitative amount σ m noise of variation of noise contributions included in the coefficient data sets ( S n,m), - averaging the wavelet coefficients of each coefficient data set for providing averaged wavelet coefficients ( S n) of the coefficient data sets, - calculating phase difference maps (Δ φ n,m) for all coefficient data sets, wherein the phase difference maps provide phase differences between the phase of each wavelet coefficient and the phase of the averaged wavelet coefficients ( S n), - calculating scaled averaged coefficient data sets of wavelet coefficients by scaling the averaged wavelet coefficients ( S n) with scaling factors ( C n,m), which are obtained by comparing parts of the normalized wavelet coefficients of the normalized coefficient data sets S n, m norm that are in phase with the averaged wavelet coefficients ( S n) (S22, S24), - calculating rescaled coefficient data sets of wavelet coefficients S n, m new by applying a transfer function T | S n, m norm |, Δ φ n, m on the coefficient data sets ( S n,m) and on the scaled averaged coefficient data sets, wherein the transfer function includes combined amplitude and phase filters, each depending on the normalized coefficient data sets S n, m norm and the phase difference maps (Δ φ n,m), resp., and - subjecting the rescaled coefficient data sets to a wavelet reconstruction S n, m new for providing the denoised MR images. - providing image data sets including multiple complex MR images, - subjecting the MR images to a wavelet decomposition for creating coefficient data sets of wavelet coefficients ( S n,m) representing the MR images in a wavelet frequency domain, - calculating normalized coefficient data sets of wavelet coefficients S n, m norm, wherein the coefficient data sets are normalized with a quantitative amount σ m noise of variation of noise contributions included in the coefficient data sets ( S n,m), - averaging the wavelet coefficients of each coefficient data set for providing averaged wavelet coefficients ( S n) of the coefficient data sets, - calculating phase difference maps (Δ φ n,m) for all coefficient data sets, wherein the phase difference maps provide phase differences between the phase of each wavelet coefficient and the phase of the averaged wavelet coefficients ( S n), - calculating scaled averaged coefficient data sets of wavelet coefficients by scaling the averaged wavelet coefficients ( S n) with scaling factors ( C n,m), which are obtained by comparing parts of the normalized wavelet coefficients of the normalized coefficient data sets S n, m norm that are in phase with the averaged wavelet coefficients ( S n) (S22, S24), - calculating rescaled coefficient data sets of wavelet coefficients S n, m new by applying a transfer function T | S n, m norm |, Δ φ n, m on the coefficient data sets ( S n,m) and on the scaled averaged coefficient data sets, wherein the transfer function includes combined amplitude and phase filters, each depending on the normalized coefficient data sets S n, m norm and the phase difference maps (Δ φ n,m), resp., and - subjecting the rescaled coefficient data sets to a wavelet reconstruction S n, m new for providing the denoised MR images.