Compression of High-Dimensional Multispectral Image Time Series Using Tensor Decomposition Learning

Multispectral imaging is widely used in many fields, such as in medicine and earth observation, as it provides valuable spatial, spectral and temporal information about the scene. It is of paramount importance that the large amount of images collected over time, and organized in multidimensional arrays known as tensors, be efficiently compressed in order to be stored or transmitted. In this paper, we present a compression algorithm which involves a training process and employs a symbol encoding dictionary. During training, we derive specially structured tensors from a given image time sequence using the CANDECOMP/PARAFAC (CP) decomposition. During runtime, every new image time sequence is quantized and encoded into a vector of coefficients corresponding to the learned CP decomposition. Experimental results on sequences of real satellite images demonstrate that we can efficiently handle higher-order tensors and obtain the decompressed data by composing the learned tensors by means of the received vector of coefficients, thus achieving a high compression ratio.

[1]  Bo Du,et al.  Compression of hyperspectral remote sensing images by tensor approach , 2015, Neurocomputing.

[2]  Panagiotis Tsakalides,et al.  Matrix and Tensor Signal Modelling in Cyber Physical Systems , 2018 .

[3]  D. Tralli,et al.  Satellite remote sensing of earthquake, volcano, flood, landslide and coastal inundation hazards , 2005 .

[4]  Tian Siva Tian,et al.  Functional Data Analysis in Brain Imaging Studies , 2010, Front. Psychology.

[5]  Alan Wells,et al.  Time series modeling of live-cell shape dynamics for image-based phenotypic profiling. , 2015, Integrative biology : quantitative biosciences from nano to macro.

[6]  Jingwen Yan,et al.  Karhunen-Loève transform for compressive sampling hyperspectral images , 2015 .

[7]  Magdy A. Bayoumi,et al.  Three-dimensional discrete wavelet transform architectures , 2002, IEEE Trans. Signal Process..

[8]  J. Chang,et al.  Analysis of individual differences in multidimensional scaling via an n-way generalization of “Eckart-Young” decomposition , 1970 .

[9]  Qian Du,et al.  Hyperspectral Image Compression Using JPEG2000 and Principal Component Analysis , 2007, IEEE Geoscience and Remote Sensing Letters.

[10]  N. Sidiropoulos,et al.  On the uniqueness of multilinear decomposition of N‐way arrays , 2000 .

[11]  T. Kolda Multilinear operators for higher-order decompositions , 2006 .

[12]  Tamara G. Kolda,et al.  Tensor Decompositions and Applications , 2009, SIAM Rev..

[13]  Jocelyn Chanussot,et al.  Compression-based nonnegative tensor CP decomposition of hyperspectral big data , 2015 .

[14]  Nikos D. Sidiropoulos,et al.  Tensor Decomposition for Signal Processing and Machine Learning , 2016, IEEE Transactions on Signal Processing.

[15]  Panagiotis Tsakalides,et al.  Matrix and Tensor Completion on a Human Activity Recognition Framework , 2017, IEEE Journal of Biomedical and Health Informatics.

[16]  Thomas Guyet,et al.  Long term analysis of time series of satellite images , 2016, Pattern Recognit. Lett..

[17]  A. Aggoun Compression of 3D Integral Images Using 3D Wavelet Transform , 2011, Journal of Display Technology.

[18]  Bo Du,et al.  PLTD: Patch-Based Low-Rank Tensor Decomposition for Hyperspectral Images , 2017, IEEE Transactions on Multimedia.

[19]  Leyuan Fang,et al.  CP tensor-based compression of hyperspectral images. , 2017, Journal of the Optical Society of America. A, Optics, image science, and vision.

[20]  Johan Håstad,et al.  Tensor Rank is NP-Complete , 1989, ICALP.

[21]  Wei Wu,et al.  JPEG2000 Still Image Coding Quality , 2013, Journal of Digital Imaging.

[22]  Andrzej Cichocki,et al.  Tensor Decompositions for Signal Processing Applications: From two-way to multiway component analysis , 2014, IEEE Signal Processing Magazine.