© 2019 IEEE. Image classification has been at the core of remote sensing applications. Optical remote sensing imaging systems naturally acquire images with spectral features corresponding to pixels. Spectral classification ignores the spatial distribution of the data which is becoming more relevant with the development of spatial resolution sensors, and many works aim to incorporate spatial features based on neighborhood through for example, Mathematical Morphology (MM). Additionally, one could stack multiple morphological transformations of the image resulting in a highly complex block of data. since classification is a tool that requires a matrix of samples and features, and simply stacking the different sets of features can lead to the problem of high dimensionality, we propose a way to create a matrix of low dimensional feature space by modeling the data as tensors and thanks to Canonical Polyadic (CP) decomposition. Experiments on real image show the effectiveness of the proposed method.