Open access peer-reviewed chapter

Compression of High-Resolution Satellite Images Using Optical Image Processing

Written By

Anirban Patra, Arijit Saha, Debasish Chakraborty and Kallol Bhattacharya

Submitted: 04 March 2020 Reviewed: 22 September 2020 Published: 05 February 2021

DOI: 10.5772/intechopen.94147

From the Edited Volume

Satellite Systems - Design, Modeling, Simulation and Analysis

Edited by Tien Nguyen

Chapter metrics overview

903 Chapter Downloads

View Full Metrics


This chapter presents a novel method for compressing satellite imagery using phase grating to facilitate the optimization of storage space and bandwidth in satellite communication. In this research work, each Satellite image is first modulated with high grating frequency in a fixed orientation. Due to this modulation, three spots (spectrum) have been generated. From these three spots, by applying Inverse Fourier Transform in any one band, we can recover the image. Out of these three spots, one is center spectrum spot and other spots represent two sidebands. Care should be taken during the spot selection is to avoid aliasing effect. At the receiving end, to recover image we use only one spectrum. We have proved that size of the extracted image is less than the original image. In this way, compression of satellite image has been performed. To measure quality of the output images, PSNR value has been calculated and compared this value with previous techniques. As high-resolution satellite image contains a lot of information, therefore to get detail information from extracted image, compression ratio should be as minimum as possible.


  • image compression
  • grating
  • Fourier transform
  • image retrieval
  • LISS- III sensor

1. Introduction

Compression of images is an important application in the field of satellite image processing as it is suitable for optimization of storage space and sharing over internet with optimum bandwidth utilization. For compression of satellite images, it is performed either directly from the image or from transformed part of the images. As discussed in [1] the compression of satellite images is based on Block Truncation Coding (BTC) technique. It first converts RGB satellite image into HSV planes. After that, each of the H and S planes are encoded using block truncation coding with quad clustering and V plane is encoded with BTC based bi-clustering or tri clustering depending on the edge information present in the plane. This method is better than previous BTC methods compared to visual quality of the output image. [2] discussed the image compression method based on evidence theory and k-Nearest Neighbor (KNN) algorithm. The main drawback is that the information loss is large [2]. To improve the quality of the output image, Fourier Transform and Huffman Coding is used for modification the previous technique. In both method, visual quality of satellite image is poor. [3] discusses the use of integer wavelet regression by increasing temporal correlation, which consequently improves the compression gain. [4] discusses a satellite image compression technique using discrete wavelet transform for noise removal to compress satellite images. [5] has discussed the use only hardware-based solutions in this lossless compression technique of X Sat images. [6] has analyzed the use of Discrete Wavelet Transform in their lossy image compression work and performance of different wavelets for satellite image compression. [7] has used the conventional Discrete Cosine Transform system for lossless image compression. [8] proposes an image compression technique using multiplexing and encryption by optical grating method [9, 10, 11, 12].

In this chapter, we proposed a scheme to compress multiple high-resolution satellite images by using phase grating. Each image is modulated by applying high value of spatial frequency and a fixed orientation angles. For each image, multiple bands have been generated due to modulation which are placed in the same spectrum plane (only three bands are clearly visible). The spectrum is encoded and filtered using Gauss filtering. To detect the location of maximum image information, an intensity graph has been plotted in the decrypted plane. All stored images can be securely and efficiently retrieved by applying inverse Regional Fourier Transform operation. This proposed technique is simple and suitable for optimization of storage space and bandwidth in satellite communication.

The chapter is organized as follows:

  • Section I describes the location and data used in this chapter

  • Section II discusses the proposed compression method

  • Section III presents proposed methodology of our research work

  • Section IV provides result of our research work and PSNR value of the extracted images

  • Section V concludes the paper and discuss why low compression ratio is desired for land cover analysis.


2. Description of location and data used

Images which are used in our research work, collected from Regional Remote Sensing Centre (East). The images are satellite picture of different areas in Kolkata Metropolitan Area.

All satellite images used in this paper are captured by LISS III 23 m sensor. LISS- III sensor is an optical sensor working in four spectral bands (Green, Red, Near Infrared and Short-Wave Infrared). It covers a 141 km- wide swath with a resolution of 23 meters in all spectral bands [13, 14].


3. Proposed methodology

3.1 Frequency and orientation angle selection for phase grating

According to rule of phase grating, value of the grating frequency (𝑢0) should be high. Low grating frequency is creating aliasing problem and therefore it would be very difficult to reconstruct the original image. In this chapter, we propose to select 𝑢0 = 1400, which is sufficient for filtering. In grating, value of the orientation angle (θ) varies from 0 to 3600. We have worked with 0 deg orientation angle.

The diffraction gratings used are illustrated in Figures 1 and 2.

Figure 1.

Sinusoidal phase grating along x axis.

Figure 2.

Spectrum of an image due to modulation.

3.2 Mathematical expression

Let us assume that image, denoted by f1xy, is modulated by two sinusoidal phase gratings along the reference abscissa and of the form expim2sin2πu1x where u1 is the spatial frequency and m is the phase contrast. In such a case, the first modulated object may be expressed as:


Invoking the well-known convolution theorem, the Fourier transform of this modulated object is given by,


whereJqistheqthorderBesselFunctionofthe 1st kind. The diffraction order is represented by the parameter ‘q’.

The diffraction pattern, as given by Eq. (2) is a series of diffraction spots each containing the object spectrum. Considering the zero order and the first two orders of the spectrum, Eq. (2) may be represented by,


3.3 Retrieval of images from spectral band

The reconstruction of images is achieved by decrypting the encrypted plane followed by plotting an intensity graph of the spectrum horizontally. This intensity graph is generated using the intensity values with respect to its location (pixel position). In the intensity graph, peak values indicate the brightest spots, where information of images are maximum. To avoid any human intervention, system automatically finds the intensity level horizontally, vertically and diagonally. In addition, regional Inverse Fourier Transform is applied around the peaks by using a fixed cut-off frequency. As only one sideband is enough for image reconstruction, hence we select one side-band. This spectrum is filtering using a cut-off value of 120. Proper operating frequency is selected and there is no overlapping during the band selection, so the extracted images are free from aliasing problem. As center band contains the information of all images, therefore filtering is not applied across center band.

3.4 Quality checking of filtered image

PSNR value is calculated using Mean Square Error (MSE) approach:


Where fxyandbxy are original image and retrieved image respectively and m, n denotes size of images.


4. Results and discussions

Satellite Images chosen for testing the algorithm described in Section III are shown in Figure 3(a)-(c). The dimensions of the selected images are 512 x 512.

Figure 3.

Three high Resolution Satellite Images.

During zonal filtering operation, images 𝑓1(𝑥, 𝑦), 𝑓2(𝑥, 𝑦) and 𝑓3 (𝑥, 𝑦) have been extracted by Regional Inverse Fourier Transform taking upper spectrum from horizontal direction. Extracted images are shown in Figure 4(a)-(c), respectively. A summary of the PSNR calculation is presented in Table 1.

Figure 4.

Three extracted images.

ImageOriginal image size (kB)Compression ratio (compressed image/original image)PSNR
𝑓1(𝑥, 𝑦)2420.8432.5
𝑓2(𝑥, 𝑦)2400.8131.73
𝑓3 (𝑥, 𝑦)2540.8631.6

Table 1.

Summary of PSNR calculation.


5. Conclusion

In this chapter, phase grating technique has been proposed for compressing the high- resolution satellite images in frequency domain. The original image is retrieved by applying Inverse Fourier Transform from the respective spectrum of the image. As presented here, since we have taken only few coefficients from the spectrum, the size of the output image is less than the main original image. It should be mentioned that in high resolution satellite image, compression should be as minimum as possible. The main reason for the requirement of low compression ratio is mainly due to large geographical area representation (as these images contain a lot of information). High compression ratio is not suitable for accurate land cover analysis. To maintain the same dimension with original image and to avoid aliasing effect, spectral area is carefully selected. Compared with earlier methods, visual quality of the selected satellite images is very good as it is captured by optical LISS-III sensor. Our proposed technique is simple and suitable for optimization of storage space and bandwidth in satellite communication.


  1. 1. B. Balakrishnan, S. H. Darsana1, J Mathews, M S. Nair; “Satellite/Aerial Image Compression Using Adaptive Block Truncation Coding Technique”; Journal of Indian Society of Remote Sensing, July 2018
  2. 2. K Sahnoun, N Benabadji ; “Satellite Image Compression Technique Based On The Evidence Theory”; Advanced Computing: An International Journal (ACIJ), Vol.5, No.1, January 2014
  3. 3. K Sahnoun, N Benabadji ; “Satellite Image Compression Algorithm Based On The FFT” ; The International Journal of Multimedia & Its Applications (IJMA) Vol.6, No.1, February 2014
  4. 4. Md. Al Mamun and Md. Ali Hossain ; “Satellite Image Compression Using Integer Wavelet Regression” ; International Conference on Electrical, Computer and Communication Engineering (ECCE) ; 2017
  5. 5. K Sahnoun, N Benabadji ; “Satellite Image Compression Technique Using Noise Bit Removal and Discrete Wavelet Transform”; International Journal of Imaging and Robotics; Volume 15, Issue Number 3 , 2015
  6. 6. R.M. Susilo , T.R. Bretschneider; “On the real time satellite image compression of X-Sat” ; Fourth International Conference on Information, Communications and Signal Processing, 2003
  7. 7. T Memane and S D Ruikar ; “Selection Of Wavelet For Satellite Image Compression Using Picture Quality Measures”; International Conference on Communication and Signal Processing ; 2014
  8. 8. I. Hacihaliloglu and M. Kartal ; “DCT and Wavelet Based Image Compression In Satellite Images”; Recent Advances in Space Technologies, 2003. RAST’03
  9. 9. A Patra, A Saha, K Bhattacharya ; “Multiplexing and encryption of images using phase grating and random phase mask” ; Optical Engineering 59 (3), 033105-033111; 2020
  10. 10. A Patra, D Chakraborty, A Saha, K Bhattacharya ; “ Compression of Satellite Images using Sinusoidal Amplitude Grating”: ; IJEE 11(1); 2019
  11. 11. A Patra, S Bandyopadhyay, D Chakraborty, A Saha ; “ A Novel Approach to Compression of Satellite Images Using Butterworth Filtering”; Springer LNSS; 2020
  12. 12. A Patra, D Sutradhar, S Nath, S poddar, R Nandi ; “ Study on Video Compression and Reconstruction using Downsampling and Upsampling”; International Journal of Advance Research in Science and Engineering 9(5); pp-12-17; 2020
  13. 13. Madhubala M , S.K.Mohan Rao, G. Ravindra Babu; Classification of IRS LISS-III Images by using Artificial Neural Networks; IJCA Special Issue on “Recent Trends in Image Processing and Pattern Recognition” RTIPPR, 2010
  14. 14. R Naidu, M.V.S.S Giridhar ; Un-Supervised Classification of Rice Crop using IRS LISS III Satellite Images for Wazirabad Command Area; International Journal Of Engineering Development And Research

Written By

Anirban Patra, Arijit Saha, Debasish Chakraborty and Kallol Bhattacharya

Submitted: 04 March 2020 Reviewed: 22 September 2020 Published: 05 February 2021