Hyperspectral image compression approaches

1 downloads 0 Views 1MB Size Report
Received XX Month XXXX; revised XX Month, XXXX; accepted XX Month XXXX; posted XX ..... CALIC) [49] is another extension of CALIC, it modified interband.
Hyperspectral image compression approaches: opportunities, challenges and future directions discussion RUI DUSSELAAR1, MANORANJAN PAUL2 1School 2School

of Computing and Mathematics, Charles Sturt University, Bathurst, NSW, Australia.2795 of Computing and Mathematics, Charles Sturt University, Bathurst, NSW, Australia.2795 * Corresponding author:[email protected]; mpaul@ csu.edu.au

Received XX Month XXXX; revised XX Month, XXXX; accepted XX Month XXXX; posted XX Month XXXX (Doc. ID XXXXX); published XX Month XXXX

This paper establishes the review of the recent study in the field of hyperspectral (HS) image compression approaches. Lately, image compression techniques have achieved significant advances from diverse types of coding standards/approaches. HS image compression requires an unconventional coding technique because of its unique, multiple dimensional structures. The data redundancy exists in both inter-band and in intra-band. The survey summarises current literature in inter- and intra-band compression methods. The challenges, opportunities and future research possibilities regarding HS image compression are further discussed. The experimental results are also provided for validity and applicability of the existing HS image compression techniques. OCIS codes: (100.0100 ) Image processing, (100.6890) Three-dimensional image processing; (100.2000) Digital image processing; (280.0280) Remote sensing and sensors;

Introduction Multispectral remote captures image data in a few relatively broad wavelength bands, typically 3 to 15 of bands. By comparison, HS remote sensors collect image data simultaneously in dozens or hundreds of narrow, adjacent spectral bands that cover a wide range of wavelengths extending from 400 to 2500 nanometre (nm) [1]. A normal visible light spectrum for human-vision is about 400 nm to 700 nm in wavelength. HS images contain an extensive range of spectral information this can provide observation power beyond the capability of human vision. However, this benefit also comes at a price, which is a wealthy spectrum of information also contains a huge amount of data. The problem is how to process enormous HS images while preserving the required useful information. Of course this is a fundamental problem for compression techniques. As we know data compression can be performed because there is unwanted redundancy in the data. High data redundancy may obtain a high compression ratio. On the contrary, lower data redundancy may get a low compression ratio. In the case of an HS image, it typically contains large amounts of image data redundancy. There are two types of redundancy that generally exists in an HS image, (i) Spatial redundancy: It is a redundancy corresponding to statistical dependencies among pixels; (ii) Spectral redundancy: It is a redundancy corresponding to pixels with the same spatial location in adjacent spectral bands. Spectral redundancy is a particular feature of HS image. In the past two decades research efforts have been directed at the reduction of spatial redundancy. Several compression standards have matured and formed

worldwide image compression standards like JPEG[2] or JPEG2000[3]. Nowadays, compression techniques aimed at eliminating spectral redundancy are still being researched. A desired HS image compression technique must be based on eliminating spatial redundancy as well as taking full consideration of spectral redundancy. The purpose is to reduce the spatial and spectral correlation consecutively and therefore achieve an optimised compression result.

Methodology 1.

Lossless and lossy compression methods Lossless and lossy compression methods, depending on whether the original image can be precisely re-generated from the compression data [4]. A tree structure of different HS image compression strategies is summarised in Fig 1. Lossless compression has no signal loss. The decoder is able to recover the compressed data back to the original data. Lossy compression method selectively discards part of signals. Decoder usually cannot recover data exactly but it can lead to relatively higher compression ratios. Lossless compression is used for applications that require the reconstructed image to restore to the original signal with high precision. Because of the intrinsic entropy of the data, lossless compression algorithm generally achieves modest compression ratio and cost more storage space for no loss of image data. By contrast, a lossy compression algorithm of HS image technique is for preserving essential spectral information of target objects, which attempt to balance of compression

Fig.1. Block diagram of general categorisation of compression methods efficiency and loss of information. Lossy image compression mainly uses predictive coding methods, transform coding etc. to compress data reducing correlations of pixels. Lossy compression is used when the user can tolerate some signal loss. Such as medical images, materials analysis commonly requires lossless compression; general web applications often use lossy compression. We’ll have briefly looked at the major types of compression algorithms in this subsection. Further introduction and discussion about the technical characteristics of typical compression methods are provided in subsequent paragraphs. The lossless compression methods are generally grouped into dictionary-based schemes and statistical schemes [5]. However dictionary-based schemes can be used for lossy compression [6], by means of quantization or entropy coding. Dictionary-based schemes do not require a prior knowledge of the source symbol probability distribution; they encode variable-length strings of symbols as single tokens, each token retrieve for a phrase dictionary such as Lempel-ZivWelch (LZW) compression algorithm [7]. Statistical-based schemes require distribution knowledge where the compression takes place based on the frequency of input characters. The most well-known statistic-based algorithms are Huffman Coding [25] and Arithmetic Coding [26, 27]. Moreover, another type of compression technique used Lookup Tables (LUT) [28, 29]. The LUT searches the previous band for a pixel equal to the current band in the same position, called a predictor. The predictor is used as a key to search Lookup tables to speed up the search process.

several groups of bands (GOBs). The authors apply intra-band prediction to the first band in each GOB. The hypothesis is that two blocks (8x8 pixels) located in the same position of adjacent HS bands are highly similar. Likewise, the authors of [10] observe that in HS images exists strong similarities between adjacent bands. They use the first band for intra-band prediction encoded and the remaining bands are inter-band encoded using fractal encoded. The encoding process of fractal coding is to use the similar interpolation method. An original image is divided into sub-images and then in the fractal set to find a similar sub-image. Decoding is completed by an iteration procedure. The encoding time is long because of the search used to find the self-similarities. Conversely, decoding time is relatively shorter. The image becomes resolution independent after being converted to fractal code. In summary of BS algorithm, selecting a more efficient GOBs measure is still vital to performance. There is no universal metric of GOBs that is applicable to all HS images in different wavelengths. Moreover, in some experimental results [9, 10] show that the algorithms can only achieve effective compression at a low bit rate.

Channel

2.

Inter-band / Intra-band compression methods Currently, available literature presents various methods regarding HS image compression studies. As mentioned previously, HS image is a three-dimensional data structure, capturing both spatial and spectral information. In general, two key steps are included into HS image compression processes. The first step of HS image compression is normally spectral decorrelation. We call it inter-band compression. This step is working on dimensional reduction. It is a very important stage to eliminate spectral redundancy. The second stage, mainly focuses on various types of compression encoders. It can be 2D or 3D compression encoder. That is considered as the intra-band compression method. Both inter-band compression and intra-band compression actually should be considered. In some cases, those two procedures can be combined as one step. A general HS image compression procedure is demonstrated at Fig 2. Band-to-band correlation is usually very high in HS image, removing such redundant information can achieve a significant compression ratio [8]. Band selection-based (BS) method is to select a subset of bands from HS image. Zhao, et al.[9] introduce an algorithm based on intra-band prediction and inter-band fractal encoding. HS bands are partitioned into

Fig.2. Spectral/spatial HS image compression procedure 3.

Compression methods based on different coding algorithms We have given a preliminary categorization of HS compression methods. Next, based on the specific different compression methods, we’ll have a review. (A) PCA (Principle Component Analysis) transform based inter-band compression PCA compression is probably the most popular and commonly used in HS compression. PCA is developed based on the principle of the Karhunen-Loêve transform (KLT). KLT converts M dimensional data to low-order uncorrelated components. The largest variance is

concentrated in the first N components (M0, Ѱ(x,y,z)= 4P(x,y-1,z); an estimation of local differences 𝑑̂ z(t) can be the inner product of a weight vector 𝑑̂ z(t)=W (t)Uz(t) where the difference between the local sum and the corresponding scaled original pixels is tracked and stored in a vector Uz(t) named, local difference vector. Then in the encoding stage, to encode the prediction residuals, the sample-adaptive entropy coding approach used. A predictor followed by an encoder of the prediction residuals is introduced. For an end user to use the CCSDS-123 standard they need to select a suitable combination of parameters for a specific application scenario. We have listed a summary of HS compression literature based on the above discussion in Table 1.

Experiments To observe the validity and suitability of different image compression approaches, four completely different datasets were used for experiments. We used Peak Signal to Noise Ratio (PSNR) to compare restoration results. The multiple experiments were performed by the relevant Benchmark encoders such as JPEG [2, 12], PCA-DCT [14], HEVC [62, 63] and other wavelet-based encoders such as JPEG2000 [16], 3DSPIHT [40], ASWDR [74]and EZW[40, 41]. They all are very competitive encoders for HS compression. At present, many new developed methods are actually still based on these algorithms with partial improvements. In order to compare with the state-of-the-art schemes, we also included AT3DSPECK [76], AT-3DSPIHT [46], and F. Zhao's [60] algorithms. They all standouts in term of PSNR for HS compression. Dataset Publicly available HS datasets can be a good source for HS image research and result verification. Specifically, Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) reflectance data is a benchmark imaging system. AVIRIS data is mainly collected for identifying, measuring, and monitoring of the earth's surface and atmosphere based on molecular absorption and particle scattering signatures.

Table 1. Summary of major contributions and challenges of existing HS image compression approaches. Algorithms Authors Contribution Inter-band compression KLT/PCA Transform

Du and Fowler 2007, Wang, Wu et al. 2009,

Pairwise orthogonal transform (POT),

Nian, Liu et al. 2016

Band-selection-based inter-bands compression

Wang, Gong et al. 2015, Zhu, Zhao et al. 2015, Zhao, Zhu et al. 2016

Tensor

Zhang, Zhang et al. 2015, Veganzones, Cohen et al. 2016)

Vector quantization transform

Zhang, Li et al. 2012, Zhao and Jing 2013, Li, Ren et al. 2014

Discrete cosine transform/3DCT

Rasti, Sveinsson et al. 2012, Qiao, Ren et al. 2014)

3D Wavelets Transforms

Tang and Pearlman 2006, Shingate, Sontakke et al. 2010, Zala and Parmar 2013, George and Manimekalai 2014 (Lin and Hwang 2010, Castrodad, Xing et al. 2011, Lin and Hwang 2011, Ülkü and Töreyin 2015)

PCA-based compression techniques discard part of the image signal. This may not be suitable for specific applications. Also, one of the major difficulties to apply PCA is to determine the optimal number of PCs for implementation. It overcomes the problem of KLT, such as bit depth expansion, lack of scalability and reduces memory requirements Selecting efficient GOBs measure is still vital to choose. There is no universal metric of GOBs applicable to all HS images in different wavelengths. Moreover, some of the experimental results [15, 16] show that the algorithms can only achieve effective compression at low bitrate. Tensor can well keep the spatial structure of data. But a huge amount of spectral information contained in an HS image to generate a tensor by using a multiplicative update algorithm makes computation extremely complicated and sometimes even unfeasible. To reduce the computational cost and speed up the convergence is very important for using tensor.

Intra-band compression

Dictionary learning and prediction based compression

Video Coding H.264/HEVC

Gao, Ji et al. 2014, Dusselaar, Paul et al. 2015, Podder, Paul et al. 2015

VQ can achieve a compression result by reducing colour distribution. However, the design of a codebook is considered to be a challenging problem. One of the traditional compression techniques. It can well remove the intra-band correlation but inter-band correlation is left untreated. Therefore, 3D-DCT approach is introduced. It’s an enhancement of DCT to extend the multiple dimensional process level. However, the reconstructed image gets blocking effects which severely degrades the visual quality. 3D transform is not specifically designed for features of HS image. It doesn’t fully utilise the common background of inner bands and spectral redundancy. Dictionary-based learning methods take consideration of inherent characteristics of the input data. Using only a few numbers of dictionary elements to present the dataset. It’s a data-driving method. Using a dictionary generating process is a key step to determine the compression performance. HEVC standard has shown significant improvement over state of the art transformation based still-image coding standard with spectral prediction modelling.

HS library

Provider

Table 2 publicly available HS datasets Spectral Range

AVIRIS Airborne Visible Infrared Imaging spectrometer HYDICE The Hyperspectral Digital Imagery Collection Experiment Natural scenes dataset

NASA

0.4~2.5µm

Naval Research Laboratory

0.4~2.5µm

University of Manchester Foster, D.H. etc.

0.4~0.7µm

Purdue’s Indianan Indian Pine test site

The Purdue University Research Repository

0.4~2.5µm

Scene data

SCIEN, Stanford University

0.4~2.5µm

USGS Digital Spectral Library

USGS Spectroscopy Lab

0.4~2.5µm

Description

AVIRIS is the most widely used HS database. AVIRIS data is basically related to the global environment and climate change. http://aviris.jpl.nasa.gov/ HYDICE is a 206 channel imaging spectrometer. The data was collected on ERIM CV 580 aircraft and a C-141 aircraft at altitudes up to 40,000 feet [76]. The natural scenes dataset consist of a mixture of rural scenes from the Minho region of Portugal, containing rocks, trees, leaves and grass, and urban scenes from the cities of Porto and Braga. http://personalpages.manchester.ac.uk/staff/david.foster /Hyperspectral_images_of_natural_scenes_04.html Purdue’s Indianan Indian Pine test site is another most widely used HS image. It provides 12 band moderate dimension Images and 220 band HS images. https://engineering.purdue.edu/~biehl/MultiSpec/hyper spectral.html The dataset contains images of faces, landscapes and buildings. The images have been recorded with two HySpex line-scan imaging spectrometers covering the spectral ranges 0.4 to 1 micrometres and 1 to 2.5 micrometres. https://scien.stanford.edu/index.php/hyperspectralimage-data/ USGS is a very useful laboratory data library provided various remote sensing spectrometers. https://speclab.cr.usgs.gov/spectral-lib.html

Table 3. Rate distortion performance of four HS images using benchmark methods. Low Bit rates 0.5 bpppb Images (35 PSNR or more) F.zhao’s[60] Cuprite_1 Jasper Ridge Lunar Lake

52.80 52.76 52.22

Cuprite_1 Jasper Ridge Lunar Lake

55.49 56.34 54.71

AT3DSPIHT[46] 51.40 51.16 51.74

AT3DSPECK[77] 52.07 51.11 52.01

EZW

ASWDR

3D-SPIHT

39.39 36.47 33.56

39.51 38.65 29.87

39.57 39.85 31.24

47.24 45.56 44.74

48.44 48.57 49.69

High bit rate 1 bpppb 54.90 56.08 54.27

54.32 55.93 55.01

44.35 43.27 43.45

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Fig.4 Rate distortion performance of nine HS images using the JPEG2000, JPEG, PCA-DCT, and HEVC encoder techniques.

They deliver calibrated images with 224 contiguous spectral bands and approximately 10 nm spectral resolutions covering the 0.4 – 2.5-μm. For more details of HS datasets can refer to Table 2. HS images for our experiments were captured using a Brimrose HS camera at Griffith University, Australia. The resolution of the images is 1392×1204×61 i.e. 61 wavelength bands, starting from 400nm to 1000nm wavelength bands, each with approximately 10nm interval steps. On top of the dataset collected from the Brimrose HS camera, we also used natural scenes [75], AVIRIS and Indian Pine dataset. Performance evaluation used four HS datasets covering both airborne and groundbased types of HS images. The PSNR gain is drawn based on the performance of the coding standard. The PSNR is defined as: PSNR=10 × log10 (

𝑃𝑒𝑎𝑘𝑖2 √𝑀𝑆𝐸

)

(1)

Where, 𝑝𝑒𝑎𝑘𝑖 is the maximum possible pixel value of the HS image, MSE is the mean square error. In Fig.4, Flower, Vase, Scene5 and Scene6 represent HS laboratory data, which are the ground-based types HS images. Those images have a relatively shorter wavelength in comparison with the AVIRIS dataset. It is clearly seen from Fig. 4 that in Flower, Vase, Scene 5, Indian-pines HS images, JPEG2000 show less rate-distortion in the same level of bit per pixel per band (bpppb). HEVC encoder performs solidly throughout for Cuprites datasets. HEVC achieved the best result in Scene 6. To summarise, JPEG2000 encoder has more obvious advantages in the ground-based HS images. Its overall performance maintains at relatively high levels. HEVC also have a superior performance, particularly in the airborne type of HS images. PCA-DCT reconstructs HS image to a lower dimensional image, uses an orthogonal transformation. It normally uses only first few PCs and discards the rest of the components for compression. PCA-DCT-based approach selects fewer principal components and changed the physical structure of the HS image. As a consequence, PCA-DCT could not obtain higher PSNR. But PCA is still able to obtain reasonably good image quality by judging of human vision system. It has been widely used for utilities of HS image lossy compression for a number of application areas. The recently high performance methods: AT-3DSPECK [76], AT3DSPIHT [46], F. Zhao's [60] algorithm, as well as popular wavelet-based encoders 3D-SPIHT, ASWDR and EZW have been evaluated in Table 3. We focused on PSNR evaluations of recently introduced compression methods. The chosen encoders: adaptively scanned wavelet difference reduction (ASWDR), EZW and 3D-SPIHT, ASWDR is one of the most recent image compression algorithms [74]. The rate-distortion results are summarised in two categories: low bit rates at 0.5bpppb and high bit rates at 1 bpppb. F. Zhao's algorithm outperforms than other methods. Both AT-3DSPECK and AT-3DSPIHT can provide competitive results. EZW, ASWDR couldn’t get the superior performance mainly because the 2-D tree structure which couldn’t eliminate inter-band redundancy very well.

Conclusion and future research direction In this review paper, we extensively reviewed current HS image compression approaches including current standard. Because of unique features of HS images, both intra-band and inter-band compression need to be considered based on the application specific demands. The applicability of the encoder for the purpose of HS image compression is validated by our experimental results. Overall, both HEVC and JPEG2000 achieves high performance. JPEG2000 performs well in the ground-based HS images and HEVC provides an outstanding performance in the airborne type of HS images. It is worth noting HEVC is more appropriate for high-definition HS images, also it can better maintain the brightness scale of HS images due to the increase of the prediction directions. HEVC allows high-precision reconstruction along different directions. As JPEG2000 is based on the wavelet transform, the

wavelet coefficients of subbands describe the spatial frequency characteristics of the horizontal and vertical directions of the spectral bands. The wavelet coefficients of different subbands reflect the characteristics of different spatial resolution of HS spectral bands. Through multi-level wavelet decomposition, wavelet coefficients can represent both high-frequency information and low-frequency information well. In this way, JPEG2000 can keep more HS image details. That probably led to the superior performance of JPEG2000. 3D-SPIHT technique generates a relatively higher compression ratio compared with ASWDR and EZW. It’s probably because 3D-SPIHT can diminish intra-band redundancy to a smaller degree. JPEG2000 is suitable for intra-band compression method, both HEVC and 3D-SPIHT has intraband and inter-band prediction. Because SPIHT sorting from the root coefficients to their descendants. Sorting coefficients and transmission process is repeated recursively. For 3D-SPIHT algorithm, this process will require even higher memory. Besides, to adopt a floating point type of wavelet coefficient may cause a certain degree of image distortion which is not ideal for lossless compression. The future research about HS images is quite exciting. VTT Technical Research Centre of Finland has created the world's first hyperspectral mobile device by converting an iPhone camera into a new kind of optical sensor. This will bring the new possibilities of low-cost spectral imaging to consumer applications. Consumers will be able to use their mobile phones, for example, to sense food quality or monitor health [67]. Due to rapid HS image sensors development, HS videos are available in recent years. Imec (world-leading research in Nanoelectronics) presents a prototype HS video acquisition [78]. HS video compression extends to process 4 dimensional data, this is another interesting area to explore. Recently, a breakthrough in precision farming is to use multispectral crop monitoring. This technology combined a compact mapping drone with a five-band sensor and is a reliable platform to capture high-quality multi-spectral data for agricultural applications [79]. Extensively used HS videos will require further improvements in compression techniques. Acknowledgments. Special thanks to Anthony Dusselaar, Rosemarry Dusselaar, Jetze Dusselaar, and Xiao Xixun for their support.

References 1. R. B. Smith, "Introduction to hyperspectral imaging," Microimages.

Retrieved on June 30, 2008 (2006). 2. W. F. Good, G. S. Maitz, and D. Gur, "Joint photographic experts group (JPEG) compatible data compression of mammograms," Journal of digital imaging 7, 123 (1994). 3. A. Skodras, C. Christopoulos, and T. Ebrahimi, "The JPEG 2000 still image compression standard," IEEE Signal processing magazine 18, 36-58 (2001). 4. G. Motta, F. Rizzo, and J. A. Storer, Hyperspectral data compression (Springer Science & Business Media, 2006). 5. A. C. Bovik, Handbook of image and video processing (Academic press, 2010). 6. G. Dudek, P. Borys, and Z. J. Grzywna, "Lossy dictionary-based image compression method," Image and Vision Computing 25, 883-889 (2007). 7. G. G. King, C. C. Seldev, and N. A. Singh, "A Novel Compression Technique for Compound Images Using Parallel Lempel-Ziv-Welch Algorithm," Applied Mechanics and Materials 626, 44 (2014). 8. C.-I. Chang, Hyperspectral data processing: algorithm design and analysis (John Wiley & Sons, 2013). 9. D. Zhao, S. Zhu, and F. Wang, "Lossy hyperspectral image compression based on intra-band prediction and inter-band fractal encoding," Computers & Electrical Engineering 54, 494-505 (2016). 10. S. Zhu, D. Zhao, and F. Wang, "Hybrid prediction and fractal hyperspectral image compression," Mathematical Problems in Engineering 2015(2015).

11. X. Pan, R. Liu, and X. Lv, "Low-complexity compression method for hyperspectral images based on distributed source coding," IEEE Geoscience and Remote Sensing Letters 9, 224-227 (2012). 12. T. Qiao, J. Ren, M. Sun, J. Zheng, and S. Marshall, "Effective compression of hyperspectral imagery using an improved 3D DCT approach for landcover analysis in remote-sensing applications," International Journal of Remote Sensing 35, 7316-7337 (2014). 13. I. Blanes and J. Serra-Sagristà, "Cost and scalability improvements to the

Karhunen–Loêve transform for remote-sensing image coding," IEEE Transactions on Geoscience and Remote Sensing 48, 2854-2863 (2010). 14. Y. Nian, Y. Liu, and Z. Ye, "Pairwise KLT-Based Compression for Multispectral Images," Sensing and Imaging 17, 1-15 (2016). 15. S. Shahriyar, M. Paul, M. Murshed, and M. Ali, "Lossless Hyperspectral Image Compression Using Binary Tree Based Decomposition," in Digital Image Computing: Techniques and Applications (DICTA), 2016 International Conference on, (IEEE, 2016), 1-8. 16. Q. Du and J. E. Fowler, "Hyperspectral image compression using JPEG2000 and principal component analysis," IEEE Geoscience and Remote Sensing Letters 4, 201-205 (2007). 17. L. Wang, J. Wu, L. Jiao, and G. Shi, "Lossy-to-lossless hyperspectral image compression based on multiplierless reversible integer TDLT/KLT," IEEE Geoscience and remote sensing letters 6, 587-591 (2009). 18. Q. Du, N. Ly, and J. E. Fowler, "An operational approach for hyperspectral image compression," in Geoscience and Remote Sensing Symposium (IGARSS), 2012 IEEE International, (IEEE, 2012), 1357-1360. 19. B. Penna, T. Tillo, E. Magli, and G. Olmo, "Transform coding techniques for lossy hyperspectral data compression," IEEE Transactions on Geoscience and Remote Sensing 45, 1408-1421 (2007). 20. F. L. Hitchcock, "The expression of a tensor or a polyadic as a sum of products," Studies in Applied Mathematics 6, 164-189 (1927). 21. L. R. Tucker, "Some mathematical notes on three-mode factor analysis," Psychometrika 31, 279-311 (1966). 22. A. Karami, R. Heylen, and P. Scheunders, "Hyperspectral Image Compression Optimized for Spectral Unmixing," IEEE Transactions on Geoscience and Remote Sensing 54, 5884-5894 (2016). 23. R. A. Harshman, "Foundations of the PARAFAC procedure: Models and conditions for an" explanatory" multi-modal factor analysis," (1970). 24. J. D. Carroll and J.-J. Chang, "Analysis of individual differences in multidimensional scaling via an N-way generalization of “Eckart-Young” decomposition," Psychometrika 35, 283-319 (1970). 25. L. Zhang, L. Zhang, D. Tao, X. Huang, and B. Du, "Compression of hyperspectral remote sensing images by tensor approach," Neurocomputing 147, 358-363 (2015). 26. M. A. Veganzones, J. E. Cohen, R. C. Farias, J. Chanussot, and P. Comon, "Nonnegative tensor CP decomposition of hyperspectral data," IEEE Transactions on Geoscience and Remote Sensing 54, 2577-2588 (2016). 27. J. Mielikainen and B. Huang, "Lossless compression of hyperspectral images using clustered linear prediction with adaptive prediction length," IEEE geoscience and remote sensing letters 9, 1118-1121 (2012). 28. N. R. Mat Noor and T. Vladimirova, "Investigation into lossless hyperspectral image compression for satellite remote sensing," International journal of remote sensing 34, 5072-5104 (2013). 29. S. Shukla and M. Prasad, Lossy image compression: domain decomposition-based algorithms (Springer Science & Business Media, 2011). 30. K. Sayood, Introduction to data compression (Newnes, 2012). 31. M. F. M. Salleh and J. Soraghan, "A new multistage lattice vector quantization with adaptive subband thresholding for image compression," EURASIP Journal on Advances in Signal Processing 2007, 092928 (2007). 32. J. Zhang, Y. Li, K. Wang, and H. Liu, "The vector quantization for AVIRIS hyperspectral imagery compression with fixed low bitrate," in Proc. of SPIE Vol, 2012), 85140W-85141. 33. X. J. Zhao and Y. F. Jing, "The Application of Vector Quantization Algorithm in Hyperspectral Image Compression," in Advanced Materials Research, (Trans Tech Publ, 2013), 1479-1483.

34. X. Li, J. Ren, C. Zhao, T. Qiao, and S. Marshall, "Novel multivariate vector quantization for effective compression of hyperspectral imagery," Optics Communications 332, 192-200 (2014). 35. B. Rasti, J. R. Sveinsson, M. O. Ulfarsson, and J. A. Benediktsson, "Hyperspectral image denoising using 3D wavelets," in Geoscience and Remote Sensing Symposium (IGARSS), 2012 IEEE International, (IEEE, 2012), 1349-1352. 36. J. Gonzalez-Conejero, J. Bartrina-Rapesta, and J. Serra-Sagrista, "JPEG2000 encoding of remote sensing multispectral images with no-data regions," IEEE Geoscience and Remote Sensing Letters 7, 251-255 (2010). 37. H. Shen, W. D. Pan, and D. Wu, "Predictive Lossless Compression of Regions of Interest in Hyperspectral Images With No-Data Regions," IEEE Transactions on Geoscience and Remote Sensing 55, 173-182 (2017). 38. V. Shingate, T. Sontakke, and S. Talbar, "Still image compression using embedded zerotree wavelet encoding," International Journal of Computer Science & Communication 1, 21-24 (2010). 39. R. George and M. Manimekalai, "A novel approach for image compression using zero tree coding," in Electronics and Communication Systems (ICECS), 2014 International Conference on, (IEEE, 2014), 1-5. 40. M. Zala and S. Parmar, "3D Wavelet transform with SPIHT algorithm for image compression," International Journal of Application Or Innovation in Engineering & Management (IJAIEM) 2(2013). 41. X. Tang and W. A. Pearlman, "Three-dimensional wavelet-based compression of hyperspectral images," in Hyperspectral Data Compression (Springer, 2006), pp. 273-308. 42. N. Amrani, J. Serra-Sagristà, V. Laparra, M. W. Marcellin, and J. Malo, "Regression wavelet analysis for lossless coding of remote-sensing data," IEEE Transactions on Geoscience and Remote Sensing 54, 5616-5627 (2016). 43. A. Karami, M. Yazdi, and G. Mercier, "Compression of hyperspectral images using discerete wavelet transform and tucker decomposition," IEEE journal of selected topics in applied earth observations and remote sensing 5, 444450 (2012). 44. X. Tang, W. A. Pearlman, and J. W. Modestino, "Hyperspectral image compression using three-dimensional wavelet coding," in Proceedings of SPIE, 2003), 1037-1047. 45. H. H. Zayed, S. E. Kishk, and H. M. Ahmed, "3D wavelets with SPIHT coding for integral imaging compression," Int J Comput Sci Netw Secur 12, 126 (2012). 46. X. Tang, S. Cho, and W. A. Pearlman, "3D set partitioning coding methods in hyperspectral image compression," in Image Processing, 2003. ICIP 2003. Proceedings. 2003 International Conference on, (IEEE, 2003), II-239. 47. X. Wu and N. Memon, "Context-based lossless interband compressionextending CALIC," IEEE Transactions on Image Processing 9, 994-1001 (2000). 48. G. Ulacha and R. Stasiński, "New context-based adaptive linear prediction algorithm for lossless image coding," in Signals and Electronic Systems (ICSES), 2014 International Conference on, (IEEE, 2014), 1-4. 49. E. Magli, G. Olmo, and E. Quacchio, "Optimized onboard lossless and nearlossless compression of hyperspectral data using CALIC," IEEE Geoscience and remote sensing letters 1, 21-25 (2004). 50. Y. Qian, M. Ye, and J. Zhou, "Hyperspectral image classification based on structured sparse logistic regression and three-dimensional wavelet texture features," IEEE Transactions on Geoscience and Remote Sensing 51, 22762291 (2013). 51. İ. Ülkü and B. U. Töreyin, "Lossy compression of hyperspectral images using online learning based sparse coding," in Computational Intelligence for Multimedia Understanding (IWCIM), 2014 International Workshop on, (IEEE, 2014), 1-5. 52. Y. Qian and M. Ye, "Hyperspectral imagery restoration using nonlocal spectral-spatial structured sparse representation with noise estimation," IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing 6, 499-515 (2013).

53. A. S. Charles, B. A. Olshausen, and C. J. Rozell, "Learning sparse codes for hyperspectral imagery," IEEE Journal of Selected Topics in Signal Processing 5, 963-978 (2011). 54. İ. Ülkü and B. U. Töreyin, "Sparse coding of hyperspectral imagery using online learning," Signal, Image and Video Processing 9, 959-966 (2015). 55. A. Castrodad, Z. Xing, J. B. Greer, E. Bosch, L. Carin, and G. Sapiro, "Learning discriminative sparse representations for modeling, source separation, and mapping of hyperspectral imagery," IEEE Transactions on Geoscience and Remote Sensing 49, 4263-4281 (2011). 56. W. Jifara, F. Jiang, B. Zhang, H. Wang, J. Li, A. Grigorev, and S. Liu, "Hyperspectral image compression based on online learning spectral features dictionary," Multimedia Tools and Applications, 1-12 (2017). 57. C.-C. Lin and Y.-T. Hwang, "Lossless Compression of Hyperspectral Images Using Adaptive Prediction and Backward Search Schemes," J. Inf. Sci. Eng. 27, 419-435 (2011). 58. C.-C. Lin and Y.-T. Hwang, "An efficient lossless compression scheme for hyperspectral images using two-stage prediction," IEEE Geoscience and Remote Sensing Letters 7, 558-562 (2010). 59. L. Santos, S. Lopez, G. M. Callico, J. F. Lopez, and R. Sarmiento, "Performance evaluation of the H. 264/AVC video coding standard for lossy hyperspectral image compression," IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing 5, 451-461 (2012). 60. F. Zhao, G. Liu, and X. Wang, "An efficient macroblock-based diverse and flexible prediction modes selection for hyperspectral images coding," Signal Processing: Image Communication 25, 697-708 (2010). 61. F. Gao, X. Ji, C. Yan, and Q. Dai, "Compression of multispectral image using HEVC," in SPIE/COS Photonics Asia, (International Society for Optics and Photonics, 2014), 92732X-92732X-92710. 62. R. Dusselaar, M. Paul, and T. Bossomaier, "Hyperspectral image coding using Spectral Prediction Modelling in HEVC coding framework," in Image and Vision Computing New Zealand (IVCNZ), 2015 International Conference on, (IEEE, 2015), 1-6. 63. M. Paul, R. Xiao, J. Gao, and T. Bossomaier, "Reflectance Prediction Modelling for Residual-Based Hyperspectral Image Coding," PloS one 11, e0161212 (2016). 64. D. Marmanis, M. Datcu, T. Esch, and U. Stilla, "Deep learning earth observation classification using ImageNet pretrained networks," IEEE Geoscience and Remote Sensing Letters 13, 105-109 (2016). 65. W. Li, H. Fu, L. Yu, P. Gong, D. Feng, C. Li, and N. Clinton, "Stacked Autoencoder-based deep learning for remote-sensing image classification: a case study of African land-cover mapping," International Journal of Remote Sensing 37, 5632-5646 (2016). 66. N. Kussul, M. Lavreniuk, S. Skakun, and A. Shelestov, "Deep Learning Classification of Land Cover and Crop Types Using Remote Sensing Data," IEEE Geoscience and Remote Sensing Letters 14, 778-782 (2017). 67. A. Romero, C. Gatta, and G. Camps-Valls, "Unsupervised deep feature extraction for remote sensing image classification," IEEE Transactions on Geoscience and Remote Sensing 54, 1349-1362 (2016). 68. W. Zhao and S. Du, "Spectral–spatial feature extraction for hyperspectral image classification: A dimension reduction and deep learning approach," IEEE Transactions on Geoscience and Remote Sensing 54, 4544-4554 (2016). 69. Y. Chen, H. Jiang, C. Li, X. Jia, and P. Ghamisi, "Deep feature extraction and classification of hyperspectral images based on convolutional neural networks," IEEE Transactions on Geoscience and Remote Sensing 54, 62326251 (2016). 70. G. Toderici, D. Vincent, N. Johnston, S. J. Hwang, D. Minnen, J. Shor, and M. Covell, "Full resolution image compression with recurrent neural networks," arXiv preprint arXiv:1608.05148 (2016). 71. L. D. Compression, "CCSDS-121.0-B-1 Blue Book," (May, 1997). 72. B. Book, "Consultative Committee for Space Data Systems (CCSDS), Image Data Compression CCSDS 122.0-B-1, ser. Blue Book, Nov. 2005," (2005). 73. L. Multispectral and H. I. C. Standard, "CCSDS 123.0-B-1," Blue Book (2012).

74. S. Raja and A. Suruliandi, "Image compression using WDR & ASWDR techniques with different wavelet codecs," ACEEE Int. J. Inform. Technol 1, 23-26 (2011). 75. D. H. Foster, K. Amano, S. M. Nascimento, and M. J. Foster, "Frequency of metamerism in natural scenes," Josa a 23, 2359-2372 (2006). 76. P. J. Pabich, "Hyperspectral Imagery: Warfighting Through a Different Set of Eyes," (DTIC Document, 2002). 77. J. Wu, Z. Wu, and C. Wu, "Lossy to lossless compressions of hyperspectral images using three-dimensional set partitioning algorithm," Optical Engineering 45, 027005-027005-027008 (2006). 78. "Imec Hyperspectral Camera " (2013), retrieved http://www2.imec.be/be_en/press/imecnews/imechyperspectralcamera2 013.html. 79. "A Breakthrough in Precision Farming" (2017), retrieved 3017, https://www.questuav.com/media/case-study/multispectral-imagingquestuav-micasense-pix4dmapper-questuav-news/.