A multi-dimensional empirical mode decomposition method is provided. The method can be applied in image texture analysis, such as medical image analysis. The method can adaptively decompose a three-dimensional image into a number of characteristic image layers and extract characteristic images showing more noticeable textures from the layers. The method uses the physical concept of field to perform the data mode decomposition to obtain envelope and tendency estimation of multi-dimensional data. The method can also be applied in time and frequency analysis of two-dimensional data or multi-channel data.