This paper presents a 2D transposition of the Hilbert-Huang Transform (HHT), an empirical data analysis method designed for studying instantaneous amplitudes and phases of non-stationary data. The principle is to adaptively decompose an image into oscillating parts called Intrinsic Mode Functions (IMFs) using an Empirical Mode Decomposition method (EMD), and then to perform Hilbert spectral analysis on the IMFs in order to recover local amplitudes and phases. For the decomposition step, we propose a new 2D mode decomposition method based on non-smooth convex optimization, while for the instantaneous spectral analysis, we use a 2D transposition of Hilbert spectral analysis called monogenic analysis, based on Riesz transform and allowing to extract instantaneous amplitudes, phases, and orientations. The resulting 2D-HHT is validated on simulated data.

This paper provides an extension of the 1D Hilbert Huang transform for the analysis of images using recent optimization techniques. The proposed method consists of: 1) adaptively decomposing an image into oscillating parts called intrinsic mode functions (IMFs) using a mode decomposition procedure and 2) providing a local spectral analysis of the obtained IMFs in order to get the local amplitudes, frequencies, and orientations. For the decomposition step, we propose two robust 2D mode decompositions based on nonsmooth convex optimization: 1) a genuine 2D approach, which constrains the local extrema of the IMFs and 2) a pseudo-2D approach, which separately constrains the extrema of lines, columns, and diagonals. The spectral analysis step is an optimization strategy based on Prony annihilation property and applied on small square patches of the IMFs. The resulting 2D Prony–Huang transform is validated on simulated and real data.