# Towards Optimal Compression of Meteorological Data: a Case Study of Using Interval-Motivated Overestimators in Global Optimization

Total Page：16

File Type：pdf, Size：1020Kb

Towards Optimal Compression of Meteorological Data: A Case Study of Using Interval-Motivated Overestimators in Global Optimization Olga Kosheleva Department of Electrical and Computer Engineering and Department of Teacher Education, University of Texas, El Paso, TX 79968, USA, [email protected] 1 Introduction The existing image and data compression techniques try to minimize the mean square deviation between the original data f(x; y; z) and the compressed- decompressed data fe(x; y; z). In many practical situations, reconstruction that only guaranteed mean square error over the data set is unacceptable. For example, if we use the meteorological data to plan a best trajectory for a plane, then what we really want to know are the meteorological parameters such as wind, temperature, and pressure along the trajectory. If along this line, the values are not reconstructed accurately enough, the plane may crash; the fact that on average, we get a good reconstruction, does not help. In general, what we need is a compression that guarantees that for each (x; y), the di®erence jf(x; y; z)¡fe(x; y; z)j is bounded by a given value ¢, i.e., that the actual value f(x; y; z) belongs to the interval [fe(x; y; z) ¡ ¢; fe(x; y; z) + ¢]: In this chapter, we describe new e±cient techniques for data compression under such interval uncertainty. 2 Optimal Data Compression: A Problem 2.1 Data Compression Is Important At present, so much data is coming from measuring instruments that it is necessary to compress this data before storing and processing. We can gain some storage space by using lossless compression, but often, this gain is not su±cient, so we must use lossy compression as well. 2 Olga Kosheleva 2.2 Successes of 2-D Image Compression In the last decades, there has been a great progress in image and data com- pression. In particular, the JPEG2000 standard (see, e.g., [TM02]) uses the wavelet transform methods together with other e±cient compression tech- niques to provide a very e±cient compression of 2D images. Within this standard, we can select di®erent bitrates (i.e., number of bits per pixel that is required, on average, for the compressed image), and depend- ing on the bitrate, get di®erent degrees of compression. When we select the highest possible bitrate, we get the lossless compres- sions that enables us to reconstruct the original image precisely. When we decrease the bitrate, we get a lossy compression; the smaller the bitrate, the more the compressed/decompressed image will di®er from the original image. In principle, it is possible to use these compression techniques to compress 2D measurement data as well. 2.3 Compressing 3-D Data: Layer-by-Layer Approach It is also possible to compress 3D measurement data f(x; y; z), e.g., meteoro- logical measurements taken in di®erent places (x; y) at di®erent heights z. One possibility is simply to apply the 2D JPEG2000 compression to each horizontal layer f(x; y; z0). 2.4 Compressing 3-D Data: An Approach that Uses KLT Transform Another possibility, in accordance with Part 2 of JPEG2000 standard [TM02], is to ¯rst apply the Karhunen-Lo`eve (KLT) transform to each vertical line. Speci¯cally, we: ² compute the average value 1 X f¹(z) = ¢ f(x; y; z) N x;y of the analyzed quantity at a given height z, where N is the overall number of horizontal points (x; y); ² compute the covariances between di®erent heights: 1 X V (z ; z ) = ¢ (f(x; y; z ) ¡ f¹(z )) ¢ (f(x; y; z ) ¡ f¹(z )); 1 2 N 1 1 2 2 x;y ² ¯nd the eigenvalues ¸k and the eigenvectors ek(z) of the covariance matrix V (z1; z2); we sort these eigenvectors into a sequence e1(z); e2(z);::: so that j¸1j ¸ j¸2j ¸ ::: ; Towards Optimal Compression of Meteorological Data 3 ² ¯nally, we represent the original 3D data values f(x; y; z) as a linear com- bination X ¹ f(x; y; z) = f(z) + ak(x; y) ¢ ek(z) k of the eigenvectors ek(z). 2.5 An Approach that Uses KLT Transform: Details How can we implement this approach? ² The ¯rst two steps are straightforward. ² The third step (computing eigenvalues and eigenvectors) is a known com- putational problem for which many standard software packages provides an e±cient algorithmic solution. ² So, to specify how this method can be implemented, we must describe how the last step can be implemented, i.e., how we can represent the original 3D data as a linear combination of the eigenvectors. First, why is such a representation possible at all? It is known (see, e.g., [TM02]), that in general, the eigenvalues of a covariance matrix form a ortho- normal basis in the space of all Nz-dimensional vectors e = fe(z)g = (e(1); e(2); : : : ; e(Nz)): By de¯nition of a basis this means, in particular, for every (x; y), the corre- sponding di®erence vector def ¹ dx;y(z) = f(x; y; z) ¡ f(z) can be represented as a linear combination of these eigenvalues, i.e., that for every z, we have X ¹ dx;y(z) = f(x; y; z) ¡ f(z) = ak(x; y) ¢ ek(z); k where by ak(x; y), we denoted the coe±cient at ek(z) in the corresponding expansion. From this formula, we get the desired representation of f(x; y; z). How can we actually compute this representation? Since the vectors ek(z) form an orthonormal basis, we can conclude that for the expansion of the vector dx;y, each coe±cient ak(x; y) at ek(z) is a dot (scalar) product of the vector dx;y and the k-th vector ek(z) from the basis, i.e., that def X ak(x; y) = dx;y ¢ ek = dx;y(z) ¢ ek(z): z Substituting the expression for dx;y(z) into this formula, we conclude that X ¹ ak(x; y) = (f(x; y; z) ¡ f(z)) ¢ ek(z): z 4 Olga Kosheleva Comment. Actually, instead of using this explicit (thus fast-to-compute) for- mula, we can use even faster formulas of the so-called fast KLT transform (see, e.g., [TM02]). 2.6 KLT Transform: Result As a result of the KLT approach, we represent the original 3-D data as a sequence of the horizontal slices ak(x; y): ² the ¯rst slice a1(x; y) corresponds to the main (1-st) eigenvalue; ² the second slice a2(x; y) corresponds to the next (2-nd) eigenvalue; ² etc., with the overall intensity decreasing from one slice to the next one. Next, to each of these slices ak(x; y), we apply a 2D JPEG2000 compression with the appropriate bit rate bk depending on the slice k. 2.7 Decompressing 3-D Data: KLT-Based Approach To reconstruct the data, we so the following: ² First, we apply JPEG2000 decompression to each slice; as a result, we get [bk] the values eak (x; y). ² Second, based on these reconstructed slices, we now reconstruct the origi- nal 3-D data data as X e ¹ [bk] f(x; y; z) = f(z) + eak (x; y) ¢ ek(z): k 2.8 Data Compression as an Optimization Problem: an Informal Introduction Di®erent bit rate allocations bk lead to di®erent quality of (lossy) compression, i.e., to di®erent quality of the corresponding decompressed data. Among all possible compression schemes that guarantee the given recon- struction accuracy (in some reasonable sense), we want to ¯nd the scheme for which the average bitrate def 1 X b = ¢ bk Nz k is the smallest possible. In some cases, the bandwidth is limited, i.e., we know the largest possible average bitrate b0. In such cases, among all compression schemes with b · b0, we must ¯nd a one for which the compression/decompression error is the smallest possible. Towards Optimal Compression of Meteorological Data 5 2.9 JPEG2000 Approach Is Tailored Towards Image Processing and Towards Mean Square Error (MSE) In the JPEG2000-based approach, for compressing measurement data, we use image compression techniques. The main objective of image compression is to retain the quality of the image. From the viewpoint of visual image quality, the image distortion can be reasonably well described by the mean square di®erence MSE (a.k.a. L2- norm) between the original image I(x; y) and the compressed-decompressed image Ie(x; y). As a result, sometimes, under the L2-optimal compression, an image may be vastly distorted at some points (x; y), and this is OK as long as the overall mean square error is small. 2.10 For Data Compression, MSE May Be a Bad Criterion When we compress measurement results, our objective is to be able to repro- duce each individual measurement result with a certain guaranteed accuracy. In such a case, reconstruction that only guaranteed mean square error over the data set is unacceptable: for example, if we use the meteorological data to plan a best trajectory for a plane, what we really want to know are the meteorological parameters such as wind, temperature, and pressure along the trajectory. If along this line, the values are not reconstructed accurately enough, the plane may crash. The fact that on average, we get a good reconstruction, does not help. 2.11 An Appropriate Criterion for Data Compression What we need is a compression that guarantees the given accuracy for all pixels, i.e., that guarantees that the L1-norm max jf(x; y; z) ¡ fe(x; y; z)j x;y;z is small.

Recommended publications
• Bosch Intelligent Streaming
Intelligent streaming | Bitrate-optimized high-quality video 1 | 22 Intelligent streaming Bitrate-optimized high-quality video – White Paper Data subject to change without notice | May 31st, 2019 BT-SC Intelligent streaming | Bitrate-optimized high-quality video 2 | 22 Table of contents 1 Introduction 3 2 Basic structures and processes 4 2.1 Image capture process ..................................................................................................................................... 4 2.2 Image content .................................................................................................................................................... 4 2.3 Scene content .................................................................................................................................................... 5 2.4 Human vision ..................................................................................................................................................... 5 2.5 Video compression ........................................................................................................................................... 6 3 Main concepts and components 7 3.1 Image processing .............................................................................................................................................. 7 3.2 VCA ..................................................................................................................................................................... 7 3.3 Smart Encoder
[Show full text]
• Isize Bitsave:High-Quality Video Streaming at Lower Bitrates
Datasheet v20210518 iSIZE BitSave: High-Quality Video Streaming at Lower Bitrates iSIZE’s innovative BitSave technology comprises an AI-based perceptual preprocessing solution that allows conventional, third-party encoders to produce higher quality video at lower bitrate. BitSave achieves this by preprocessing the video content before it reaches the video encoder (Fig. 1), such that the output after compressing with any standard video codec is perceptually optimized with less motion artifacts or blurring, for the same or lower bitrate. BitSave neural networks are able to isolate areas of perceptual importance, such as those with high motion or detailed texture, and optimize their structure so that they are preserved better at lower bitrate by any subsequent encoder. As shown in our recent paper, when integrated as a preprocessor prior to a video encoder, BitSave is able to reduce bitrate requirements for a given quality level by up to 20% versus that encoder. When coupled with open-source AVC, HEVC and AV1 encoders, BitSave allows for up to 40% bitrate reduction at no quality compromise in comparison to leading third-party AVC, HEVC and AV1 encoding services. Lower bitrates equate to smaller filesizes, lower storage, transmission and distribution costs, reduced energy consumption and more satisfied customers. Source BitSave Encoder Delivery AI-based pre- Single pass processing processing prior One frame latency per content for an to encoding (AVC, entire ABR ladder HEVC, VP9, AV1) Improves encoding quality Integrated within Intel as measured by standard OpenVINO, ONNX and Dolby perceptual quality metrics Vision, easy to plug&play (VMAF, SSIM, VIF) within any existing workflow Figure 1.
[Show full text]
• Adaptive Bitrate Streaming Over Cellular Networks: Rate Adaptation and Data Savings Strategies
Adaptive Bitrate Streaming Over Cellular Networks: Rate Adaptation and Data Savings Strategies Yanyuan Qin, Ph.D. University of Connecticut, 2021 ABSTRACT Adaptive bitrate streaming (ABR) has become the de facto technique for video streaming over the Internet. Despite a ﬂurry of techniques, achieving high quality ABR streaming over cellular networks remains a tremendous challenge. First, the design of an ABR scheme needs to balance conﬂicting Quality of Experience (QoE) metrics such as video quality, quality changes, stalls and startup performance, which is even harder under highly dynamic bandwidth in cellular network. Second, streaming providers have been moving towards using Variable Bitrate (VBR) encodings for the video content, which introduces new challenges for ABR streaming, whose nature and implications are little understood. Third, mobile video streaming consumes a lot of data. Although many video and network providers currently oﬀer data saving options, the existing practices are suboptimal in QoE and resource usage. Last, when the audio and video arXiv:2104.01104v2 [cs.NI] 14 May 2021 tracks are stored separately, video and audio rate adaptation needs to be dynamically coordinated to achieve good overall streaming experience, which presents interesting challenges while, somewhat surprisingly, has received little attention by the research community. In this dissertation, we tackle each of the above four challenges. Firstly, we design a framework called PIA (PID-control based ABR streaming) that strategically leverages PID control concepts and novel approaches to account for the various requirements of ABR streaming. The evaluation results demonstrate that PIA outperforms state-of-the-art schemes in providing high average bitrate with signif- icantly lower bitrate changes and stalls, while incurring very small runtime overhead.
[Show full text]
Bitrate Requirements of Non-Panoramic VR Remote Rendering Viktor Kelkkanen Markus Fiedler David Lindero Blekinge Institute of Technology Blekinge Institute of Technology Ericsson Karlskrona, Sweden Karlshamn, Sweden Luleå, Sweden [email protected] [email protected] [email protected] ABSTRACT This paper shows the impact of bitrate settings on objective quality measures when streaming non-panoramic remote-rendered Virtual Reality (VR) images. Non-panoramic here refers to the images that are rendered and sent across the network, they only cover the viewport of each eye, respectively. To determine the required bitrate of remote rendering for VR, we use a server that renders a 3D-scene, encodes the resulting images using the NVENC H.264 codec and transmits them to the client across a network. The client decodes the images and displays them in the VR headset. Objective full-reference quality measures are taken by comparing the image before encoding on the server to the same image after it has been decoded on the client. By altering the average bitrate setting of the encoder, we obtain objective quality Figure 1: Example of the non-panoramic view and 3D scene scores as functions of bitrates. Furthermore, we study the impact that is rendered once per eye and streamed each frame. of headset rotation speeds, since this will also have a large effect on image quality. We determine an upper and lower bitrate limit based on headset 1 INTRODUCTION rotation speeds. The lower limit is based on a speed close to the Remote rendering for VR is commonly done by using panorama average human peak head-movement speed, 360°/s.
[Show full text]
• An Analysis of VP8, a New Video Codec for the Web Sean Cassidy
[Show full text]
• Encoding Parameters Prediction for Convex Hull Video Encoding
Encoding Parameters Prediction for Convex Hull Video Encoding Ping-Hao Wu∗, Volodymyr Kondratenko y, Gaurang Chaudhari z, and Ioannis Katsavounidis z Facebook Inc., 1 Hacker Way, Menlo Park, CA 94025 Email: [email protected]∗, [email protected], [email protected], [email protected] Abstract—Fast encoding parameter selection technique have BD-rate loss. It uses a faster encoder or the same encoder been proposed in the past. Leveraging the power of convex hull with faster speed settings to perform analysis on the video video encoding framework, an encoder with a faster speed setting, sequence ﬁrst, to determine the optimal encoding parameters or faster encoder such as hardware encoder, can be used to determine the optimal encoding parameters. It has been shown to be used with a slower encoder with better compression that one can speed up 3 to 5 times, while still achieve 20-40% efﬁciency. Although this fast selection technique can achieve BD-rate savings, compared to the traditional ﬁxed-QP encodings. similar savings with a fraction of the computations, limitations Such approach presents two problems. First, although the are still present. First of all, the few percent loss in BD-rate is speedup is impressive, there is still ∼3% loss in BD-rate. not negligible in the case of higher-quality contents, where a Secondly, the previous approach only works with encoders implementing standards that have the same quantization scheme much higher target quality and/or target bitrate is often used. and QP range, such as between VP9 and AV1, and not in the Secondly, the quantization scheme of the fast encoder, used scenario where one might want to use a much faster H.264 for analysis, and that of the slow encoder, used for producing encoder, to determine and predict the encoding parameters for the actual deliverable encodings, need to be the same in order VP9 or AV1 encoder.
[Show full text]
• H.264 Codec Comparison
MPEG-4 AVC/H.264 Video Codecs Comparison Short version of report Project head: Dr. Dmitriy Vatolin Measurements, analysis: Dmitriy Kulikov, Alexander Parshin Codecs: XviD (MPEG-4 ASP codec) MainConcept H.264 Intel H.264 x264 AMD H.264 Artemis H.264 December 2007 CS MSU Graphics&Media Lab Video Group http://www.compression.ru/video/codec_comparison/index_en.html [email protected] VIDEO MPEG-4 AVC/H.264 CODECS COMPARISON MOSCOW, DEC 2007 CS MSU GRAPHICS & MEDIA LAB VIDEO GROUP SHORT VERSION Contents 1 Acknowledgments ................................................................................................4 2 Overview ..............................................................................................................5 2.1 Difference between Short and Full Versions............................................................ 5 2.2 Sequences ............................................................................................................... 5 2.3 Codecs ..................................................................................................................... 6 3 Objectives and Testing Rules ..............................................................................7 3.1 H.264 Codec Testing Objectives.............................................................................. 7 3.2 Testing Rules ........................................................................................................... 7 4 Comparison Results.............................................................................................9
[Show full text]
• Technologies for 3D Mesh Compression: a Survey
J. Vis. Commun. Image R. 16 (2005) 688–733 www.elsevier.com/locate/jvci Technologies for 3D mesh compression: A survey Jingliang Peng a, Chang-Su Kim b,*, C.-C. Jay Kuo a a Integrated Media Systems Center, Department of Electrical Engineering, University of Southern California, Los Angeles, CA 90089-2564, USA b Department of Information Engineering, The Chinese University of Hong Kong, Shatin, Hong Kong Received 1 February 2004; accepted 5 March 2005 Available online 16 April 2005 Abstract Three-dimensional (3D) meshes have been widely used in graphic applications for the rep- resentation of 3D objects. They often require a huge amount of data for storage and/or trans- mission in the raw data format. Since most applications demand compact storage, fast transmission, and eﬃcient processing of 3D meshes, many algorithms have been proposed to compress 3D meshes eﬃciently since early 1990s. In this survey paper, we examine 3D mesh compression technologies developed over the last decade, with the main focus on triangular mesh compression technologies. In this eﬀort, we classify various algorithms into classes, describe main ideas behind each class, and compare the advantages and shortcomings of the algorithms in each class. Finally, we address some trends in the 3D mesh compression tech- nology development. Ó 2005 Elsevier Inc. All rights reserved. Keywords: 3D mesh compression; Single-rate mesh coding; Progressive mesh coding; MPEG-4 1. Introduction Graphics data are more and more widely used in various applications, including video gaming, engineering design, architectural walkthrough, virtual reality, e-com- * Corresponding author. E-mail addresses: [email protected] (J.
[Show full text]
• Delivering Live and On-Demand Smooth Streaming
Delivering Live and On-Demand Smooth Streaming Video delivered over the Internet has become an increasingly important element of the online presence of many companies. Organizations depend on Internet video for everything from corporate training and product launch webcasts, to more complex streamed content such as live televised sports events. No matter what the scale, to be successful an online video program requires a carefully-planned workflow that employs the right technology and tools, network architecture, and event management. This deployment guide was developed by capturing real world experiences from iStreamPlanet® and its use of the Microsoft® platform technologies, including Microsoft® SilverlightTM and Internet Information Services (IIS) Smooth Streaming, that have powered several successful online video events. This document aims to deliver an introductory view of the considerations for acquiring, creating, delivering, and managing live and on-demand SD and HD video. This guide will touch on a wide range of new tools and technologies which are emerging and essential components to deployed solutions, including applications and tools for live and on- demand encoding, delivery, content management (CMS) and Digital Asset Management (DAM) editing, data synchronization, and advertisement management. P1 | Delivering Live and On-Demand Streaming Media A Step-by-Step Development Process This paper was developed to be a hands-on, practical 1. Production document to provide guidance, best practices, and 2. Content Acquisition resources throughout every step of the video production 3. Encoding and delivery workflow. When delivering a live event over 4. Content Delivery the Internet, content owners need to address six key 5. Playback Experience areas, shown in Figure 1: 6.
[Show full text]
• Lossless Coding of Point Cloud Geometry Using a Deep Generative Model
1 Lossless Coding of Point Cloud Geometry using a Deep Generative Model Dat Thanh Nguyen, Maurice Quach, Student Member, IEEE, Giuseppe Valenzise, Senior Member, IEEE, Pierre Duhamel, Life Fellow, IEEE Abstract—This paper proposes a lossless point cloud (PC) been optimized over several decades. On the other hand, G- geometry compression method that uses neural networks to PCC targets static content, and the geometry and attribute estimate the probability distribution of voxel occupancy. First, information are independently encoded. Color attributes can to take into account the PC sparsity, our method adaptively partitions a point cloud into multiple voxel block sizes. This be encoded using methods based on the Region Adaptive partitioning is signalled via an octree. Second, we employ a Hierarchical Transform (RAHT) [4], Predicting Transform or deep auto-regressive generative model to estimate the occupancy Lifting Transform [3]. Coding the PC geometry is particularly probability of each voxel given the previously encoded ones. We important to convey the 3D structure of the PC, but is also then employ the estimated probabilities to code efﬁciently a block challenging, as the non-regular sampling of point clouds using a context-based arithmetic coder. Our context has variable size and can expand beyond the current block to learn more makes it difﬁcult to use conventional signal processing and accurate probabilities. We also consider using data augmentation compression tools. In this paper, we focus on lossless coding techniques to increase the generalization capability of the learned of point cloud geometry. probability models, in particular in the presence of noise and In particular, we consider the case of voxelized point clouds.
[Show full text]
• Integrated Lossy, Near-Lossless, and Lossless Compression of Medical Volumetric Data
Integrated Lossy, Near-lossless, and Lossless Compression of Medical Volumetric Data Sehoon Yeaa, Sungdae Chob and William A.Pearlmana aCenter for Image Processing Research Electrical, Computer & Systems Engineering Dept. Rensselaer Polytechnic Institute 110 8th St. Troy, NY 12180-3590 bSamsung Electronics, Suwon, South Korea ABSTRACT We propose an integrated, wavelet based, two-stage coding scheme for lossy, near-lossless and lossless compression of medical volumetric data. The method presented determines the bit-rate while encoding for the lossy layer and without any iteration. It is in the spirit of \lossy plus residual coding" and consists of a wavelet-based lossy layer followed by an arithmetic coding of the quantized residual to guarantee a given pixel-wise maximum error bound. We focus on the selection of the optimum bit rate for the lossy coder to achieve the minimum total (lossy plus residual) bit rate in the near-lossless and the lossless cases. We propose a simple and practical method to estimate online the optimal bit rate and provide a theoretical justiﬁcation for it. Experimental results show that the proposed scheme provides improved, embedded lossy, and lossless performance competitive with the best results published so far in the literature, with an added feature of near-lossless coding. Keywords: near-lossless coding, SPIHT, lossy-to-lossless, medical image 1. INTRODUCTION Most of today's medical imaging techniques produce a huge amount of 3D data. Examples include Magnetic Res- onance(MR),Computerized Tomography(CT), Position Emission Tomography(PET) and 3D ultrasound. Con- siderable research eﬀorts have been made for their eﬃcient compression to facilitate storage and transmission.
[Show full text]
• Video Coding Using Geometry Based Block Partitioning and Reordering Discrete Cosine Transform*
Zhang et al. / J Zhejiang Univ-Sci C (Comput & Electron) 2012 13(1):71-82 71 Journal of Zhejiang University-SCIENCE C (Computers & Electronics) ISSN 1869-1951 (Print); ISSN 1869-196X (Online) www.zju.edu.cn/jzus; www.springerlink.com E-mail: [email protected] Video coding using geometry based block partitioning and * reordering discrete cosine transform Yi-xiong ZHANG†1, Jiang-hong SHI1, Wei-dong WANG2 (1Department of Communication Engineering, Xiamen University, Xiamen 361005, China) (2Department of Information Science and Electronic Engineering, Zhejiang University, Hangzhou 310027, China) †E-mail: [email protected] Received July 27, 2011; Revision accepted Nov. 28, 2011; Crosschecked Dec. 8, 2011 Abstract: Geometry based block partitioning (GBP) has been shown to achieve better performance than the tree structure based block partitioning (TSBP) of H.264. However, the residual blocks of GBP mode after motion compensation still present some non-vertical/non-horizontal orientations, and the conventional discrete cosine transform (DCT) may generate many high-frequency coefficients. To solve this problem, in this paper we propose a video coding approach by using GBP and reordering DCT (RDCT) techniques. In the proposed approach, GBP is first applied to partition the macroblocks. Then, before performing DCT, a reordering operation is used to adjust the pixel positions of the residual macroblocks based on the partition information. In this way, the partition information of GBP can be used to represent the reordering information of RDCT, and the bitrate can be reduced. Experimental results show that, compared to H.264/AVC, the proposed method achieves on average 6.38% and 5.69% bitrate reductions at low and high bitrates, respectively.
[Show full text]