Rate-Distortion Comparisons
zfp is a very fast compressor and decompressor for floating-point data that achieves up to 2 gigabytes/second throughput. zfp also has excellent rate-distortion performance, or quality per bit of compressed storage. The images below compare the speed and quality of several published compression schemes. In addition to the compressed size in bits per value, the colored bars and numbers represent the compression ratio (uncompressed to compressed size; higher is better), peak signal to noise ratio (PSNR) in decibels (higher is better), and compression/decompression throughput in uncompressed megabytes per second (higher is better).
The original data is represented in double-precision floating point, i.e. using 64 bits/value. Notice the lack of artifacts in the zfp compressed data set, in spite of zfp using fewer bits per value than the other compressors, while also compressing and decompressing the data between 6 and 170 times faster.
The images below show using color coding the maximum (in magnitude) signed error computed along the z (vertical) direction for the data set above. This max reduction collapses the 3D volume to a 2D image. Thanks in part to zfp's fixed-accuracy mode that enforces an absolute error tolerance, it achieves low maximum errors as well.
The plot below shows rate-distortion curves for the 3D data set. Every PSNR increase of 20 dB corresponds to an additional decimal digit of accuracy. Thus zfp's 268 dB improvement over ISABELA represents an increase in accuracy (or reduction in RMS error) of more than 10 trillion times. Comparisons are also made with conversion from double to single (float) precision; truncation (zeroing) of least significant mantissa bits; and uniform scalar quantization.
The three green curves correspond to different modes of using zfp. In fixed-rate mode, the 3D array is partitioned into blocks of 4x4x4 values, and each block is compressed to a fixed number of bits, which enables random access at block granularity. In fixed-precision mode, the number of uncompressed bits per value is fixed, resulting in bounded relative error. In fixed-accuracy mode, each reconstructed value is guaranteed to fall within a user-specified absolute error tolerance, which usually results in the smallest RMS error and therefore highest peak signal to noise ratio for the same rate.
References
- ISABELA 0.2.1: Lakshminarasimhan et al., "Compressing the incompressible with ISABELA: In-situ reduction of spatio-temporal data," Euro-Par Parallel Processing 2011.
- fpzip 1.1.0: Lindstrom and Isenburg, "Fast and efficient compression of floating-point data," IEEE Transactions on Visualization and Computer Graphics, 12(5), 2006.
- HVQ 1.31: Schneider and Westermann, "Compression domain volume rendering," IEEE Visualization 2003.
- VAPOR 2.4.2: Clyne et al., "Interactive desktop analysis of high resolution simulations: Application to turbulent plume dynamics and current sheet formation," New Journal of Physics, 9(8), 2007.
- SQ: Iverson et al., "Fast and effective lossy compression algorithms for scientific datasets," Euro-Par Parallel Processing 2012.
- SZ 1.3: Di and Cappello, "Fast error-bounded lossy HPC data compression with SZ," IEEE International Parallel and Distributed Processing Symposium 2016.
- JPEG 2000: Woodring et al., "Revisiting wavelet compression for large-scale climate data using JPEG 2000 and ensuring data precision," LDAV 2011.
- zfp 0.5.0: Lindstrom, "Fixed-rate compressed floating-point arrays," IEEE Transactions on Visualization and Computer Graphics, 20(12), 2014.