New rate adaptation method for JPEG2000-based SNR Scalable Video Coding with Integer Linear Programming models

Size: px
Start display at page:

Download "New rate adaptation method for JPEG2000-based SNR Scalable Video Coding with Integer Linear Programming models"

Transcription

1 New rate adaptaton method for JPEG2000-based SNR Scalable Vdeo Codng wth Integer Lnear Programmng models Lvo Lma #1, Renata Mansn, Rccardo Leonard #2 # Sgnal Processng Group Department of Electroncs for Automaton, Unversty of Bresca Va Branze, Bresca, Italy 1 lvo.lma@ng.unbs.t 2 rccardo.leonard@ng.unbs.t Operatons Research Group Department of Electroncs for Automaton, Unversty of Bresca Va Branze, Bresca, Italy renata.mansn@ng.unbs.t Abstract In the last few years scalable vdeo codng emerged as a promsng technology for effcent dstrbuton of vdeos through heterogeneous networks. In a heterogeneous envronment, the vdeo content needs to be adapted n order to meet dfferent end termnal capablty requrements (user adaptaton) or fluctuatons of the avalable bandwdth (network adaptaton). Consequently, the adaptaton problem s a crtcal ssue n scalable vdeo codng desgn. In ths paper we ntroduce a new adaptaton method for a proposed JPEG2000-based SNR scalable codec, that formulates and solves the adaptaton problem as an Integer Lnear Programmng problem. Index Terms Scalable vdeo codng, rate adaptaton, JPEG2000, Integer Lnear Programmng. I. INTRODUCTION In the last years, scalable vdeo codng emerged as a promsng technology for effcent dstrbuton of vdeos through heterogeneous networks, and t has been recently standardzed as scalable extenson of the H.264/AVC standard [1], hereafter ndcated as SVC. An useful overvew of the SVC extenson can be found n [2]. The man advantage of SVC s that t offers codng flexblty to decode dfferent workng ponts n terms of spatal, temporal and qualty resoluton from a unque coded representaton. In a heterogeneous envronment typcally the vdeo content needs to be adapted n order to meet dfferent end termnal capablty requrements (user adaptaton) or fluctuatons of the avalable bandwdth (network or rate adaptaton). In ths work we address only the rate adaptaton problem. The scalablty features gven by the scalable vdeo codng offer a very flexble way to perform the adaptaton, for example reducng the spatal resoluton or the vdeo qualty. In partcular, wth the SVC scalablty features, rate adaptaton can be effcently managed usng qualty scalablty,.e. coarse gran (CGS) or fne gran (FGS) scalablty. Although the ntrnsc support for adaptaton provded wth the scalablty, the problem of the adaptaton for scalable vdeo content s an open research ssue that s stll under nvestgaton. An exhaustve survey on the proposed approaches for solvng the adaptaton problem can be found n [3]. In [3] dfferent classfcaton crtera and propertes of the adaptaton methods are presented. For the purpose of the present work, three aspects have to be consdered when comparng dfferent adaptaton methods: the performance, evaluated n terms of decoded qualty (dependng on a partcular metrcs) of the extracted data. the satsfacton of addtonal constrants on the decoded vdeo sequences. the complexty of the adaptaton process. Among the dfferent approaches presented n [3] the frst and most attractve approach proposed for rate adaptaton wth SVC s that proposed n [4]. Recently, another nterestng approach has been proposed n [5]. In [5] the authors compare ther approach to that n [4] and they show that the two approaches have smlar extracton performance. Snce the mplementaton of the approach n [4] s ncluded n the SVC reference software (JSVM 9, verson 9.14) [6] we decde to use only ths approach as reference for evaluatng the performance of our method, descrbed n the followng sectons. Although SVC s the reference standard for scalable vdeo codng, other scalable vdeo codng approaches has been proposed n lterature. Most of them are based on the wavelet technology, that has natve spatal scalablty features. Furthermore, JPEG2000 s the state-of-the-art n stll mage compresson and t offers very effcent spatal and SNR scalablty. Inspred by the potentaltes of JPEG2000 and the prevous wavelet approaches we have proposed a very smple but effcent new scalable vdeo codng soluton based on JPEG2000. In Secton II we defne the archtecture of the proposed codec and the man features of JPEG2000 that are fundamental to understand

2 Fg. 1. Herarchcal B-pcture decomposton moton nfo ntra-frame MC-resdual moton nfo codng JPEG2000 codng Proposed JPEG2000-based codec archtecture Multplex the proposed adaptaton soluton. At present, the proposed archtecture s stll under nvestgaton relatvely to the spatal scalablty. Thus, n ths work, we only focus on the SNR scalablty, also avodng to consder the temporal scalablty for adaptaton purposes. The man dea behnd the proposed approach s that the rate adaptaton problem for SNR scalable codng can be seen as an optmal resources allocaton problem and thus formulated usng Integer Lnear Programmng (ILP) as follows: mnmze { subject to c T x Ax b x 0 nteger where c T x s the objectve functon and descrbes the target of the adaptaton problem, Ax b represents the constrants gven by the codng system and addtonal constrants on the decoded vdeo sequence and x s a vector of bnary unknowns descrbng how resources are allocated. In Secton III we descrbe the proposed nteger programmng problem for rate adaptaton and underlne some of ts propertes. Although NP-hard, the problem (1) can be effcently solved by means of a commercal software for mxed nteger lnear programmng such as CPLEX and explotng partcular propertes of the model. Fnally, Secton IV s devoted to compare our approach performance wth that of the method n [4], whereas n Secton V conclusons and possble future developments are drawn. II. PROPOSED CODEC SOLUTION The proposed scalable vdeo codng system, that at present only offers temporal and SNR scalablty features, s based on a very smple archtecture shown n Fgure 1. The nput vdeo sequence s temporally decomposed usng a smlar approach to Herarchcal B-pcture decomposton proposed for SVC (see [2] and [7]), that enables closed-loop moton estmaton and natve temporal scalablty. Herarchcal B-pcture decomposton s based on the defnton of key pcture and group of pcture (GOP). The frst pcture of a vdeo sequence s an ntra-coded pcture, whle the other key-pctures could be ntra-coded or nter coded usng the prevous key pctures as reference. Usually key-pctures are coded at regular ntervals, and the dstance between two key-pctures defnes the GOP length. In fact a key pcture and all pctures that are temporally located between the key pcture and the prevous key pcture are consdered to buld a group of pctures. Fgure 2 shows a typcal decomposton structure based on herarchcal B- pcture for a GOP wth length equal to 8. Ths partcular GOP structure enables 4 levels of temporal scalablty. In the frst (1) one, whch we call level 0 and represents the sequence at ts lowest avalable frame rate, we consder only the sequence made of key pctures, whle n a generc level ( 1) we consder the pctures used at lower temporal levels plus the pctures ndcated n Fgure 2 as B level pctures. For each GOP there s only 1 B level 1 pcture and t s predcted usng the key pcture of the prevous GOP for forward predcton and the key pcture of the same GOP for backward predcton. For > 1 the B level pctures are predcted usng the pctures belongng to the lower temporal resoluton. It should be notced that the herarchcal B-pcture structure enables closed-loop codng. In fact, the encodng order s dfferent from the dsplay order. The frst encoded pcture s the frst frame of the sequence, then for every GOP we encode the key pcture before moton estmaton and compensaton for the B level 1 pcture. Then every B level pcture s encoded before moton estmaton and compensaton of the B level + 1 pctures. Ths encodng order ensures that at each temporal level the moton estmaton and compensaton process uses the already encoded reference pctures. In the proposed scalable vdeo codec the herarchcal B-pcture decomposton s adopted wth the only constrant that the keypcture of each GOP can be only ntra-coded. After the temporal decomposton, for each GOP, the keypcture and the moto-compensaton resdual for all the B- pctures wthn a GOP are encoded wth a JPEG2000 framework as a sngle pcture. SNR scalablty s obtaned generatng JPEG2000 codestreams wth multple qualty layers. Snce the proposed work addresses only the SNR scalablty, n the followng we overvew the man features of JPEG2000 that are fundamental to understand the model proposed n Secton III. JPEG2000 s the state-of-the-art n mage compresson and s based on the Dscrete Wavelet Transform (DWT), together wth Embedded Block Codng wth Optmzed Truncaton (EBCOT) [8]. D stages of DWT analyss decompose the mage nto 3D+1 subbands, labeled LH d, HL d, HH d and LL D, for d = 1,..., D. An useful overvew of the man features of JPEG2000 can be found n [9], whle for a complete techncal descrpton the reader s referred to [10]. Each subband s parttoned nto rectangular blocks called code-blocks, each of whch s ndependently coded. Resoluton scalablty s obtaned by dscardng the code-blocks of detal subbands and omttng the fnal DWT synthess stage. Qualty scalablty s obtaned through a qualty layers abstracton. Each layer represents an ncremental contrbuton (possbly empty) from the embedded bt-stream assocated wth each code-block n the mage. Dscardng one or more layers (startng from the hghest one) produces a representaton of the code-block wth lower qualty. JPEG2000 also defnes collectons of spatally adjacent code-blocks as precncts. Each precnct of resoluton level LL d consst of the code-blocks correspondng to the same spatal regon wthn the subbands LH d+1, HL d+1 and HH d+1 f d < D, or wthn the subband LL D f d = D. The data-stream assocated wth each precnct s organzed as a

3 key-pc B-lev-3 B-lev-2 B-lev-3 B-lev-1 B-lev-3 B-lev-2 B-lev-3 key-pc f P f 1 f 2 f 3 f 4 f 5 f 6 f 7 f 8 GOP border GOP border Fg. 2. Herarchcal B-pctures temporal decomposton wth GOP sze equal to 8 collecton of packets, one for each qualty layer. III. RATE EXTRACTION OPTIMIZATION As prevously descrbed, the rate adaptaton problem can be formulated as a problem of optmal allocaton of resources. Snce n the proposed scalable vdeo codec the frames (ntraframes and moto-compensated resdual) are encoded usng JPEG2000, the resources that have to be allocated are the JPEG2000 codng elements. In the proposed nteger model the mnmum addressable codng element s the JPEG2000 packet, whereas we avod to consder for the adaptaton the moton nformaton. In the proposed approach the adaptaton s performed ndependently on each GOP (wth length equal to F frames and F s a power of 2), and consequently the optmal allocaton of the JPEG2000 packets has to be performed n order to mnmze the overall dstorton over the GOP, gven by: F D T OT = D t (2) t=1 where D t s the dstorton for the frame t evaluated as MSE between the orgnal and the decoded frame. The total dstorton (2) represents the objectve functon c T x of the model (1). The budget avalable for the allocaton s gven by the bandwdth or equvalently by the data-rate. Gven the GOP length and the frame-rate of the sequence, t s always possble to determne the amount of bts L that can be used to encode each GOP, where constant bt-rate transmsson s consdered. In secton III-A the expressons for the dstorton of the keypcture and the B-pctures wthn each GOP wll be provded, avodng the mathematcal detals that can be found n [11], whle n secton III-B the ILP model wll be presented. A. Dstorton computaton Frst, we ntroduce the dstorton contrbuton gven by a sngle JPEG2000 packet. Let I be the number of precncts n each frame (supposed to be constant over the frames), K the number of qualty layers ncluded n the data-stream assocated wth each precnct and F the GOP length (n frames). We defne the followng quanttes: P t a generc precnct, = 1,..., I, belongng to frame t, t = 1,..., F, P t,k the decoded verson of P t usng k, k = 1,..., K, qualty layers and Lt,k the sze (n bts) of the frst k packets related to the precnct P t. The dstorton ntroduced approxmatng Pt wth Pt,k s gven by D t,k = P t Pt,k 2. Introducng the rate-dstorton slope S t,n, defned as the rato D/ L related to each qualty layer of each code-block, the dstorton D t,k can be expressed as: D t,k = D t,k + S t,n L t,n (3) n=k+1 where D t,k s the dstorton experenced (possbly equal to 0) f all the qualty layers are consdered, L t,n = L t,n L t,n 1 s the sze (n bts) of the qualty layer n and S t,n L t,n s the dstorton contrbuton gven by the layer n. Unfortunately, the exact calculaton of the dstorton ntroduced for each precnct P t requre the knowledge of the rate-dstorton slopes S t,n. Typcally, n order to reduce the overhead requred to mantan ths slope nformaton, only the rate-dstorton slope threshold values used for the layer generaton are ncluded nto the JPEG2000 codestream header (see [10] for more detals). Ths means that the rate dstorton slope for a partcular layer n s consdered to be constant over the precncts and equal to the threshold Tt n. Usng the rate-dstorton slope thresholds the dstorton (3) can be approxmated as: ˆD t,k = D t,k + n=k+1 T n t L t,n (4) In order to estmate the overall dstorton D T OT, the man assumpton that s consdered hereafter s that the JPEG2000 packets are ndependent, that s approxmately true for the 9/7 tap borthogonal flters typcally used n JPEG2000. Ths leads to an addtve model for the dstorton of a frame so that the sum of all the precncts dstorton contrbuton I =1 Dt, where the sngle contrbuton Dt s estmated usng the equaton (4), can be consdered as a good approxmaton of D t. In order to estmate the dstorton of a generc frame t nsde

4 2 220 R-D Layer R-D Layer 3 18 R-D Layer proposed dstorton model measured dstorton proposed dstorton model measured dstorton proposed dstorton model measured dstorton rate [byte] rate [byte] rate [byte] (a) R-D performance of layer 1 (b) R-D performance of layer 3 (c) R-D performance of layer 7 Fg. 3. Typcal Rate-Dstorton performance of a key-pcture encoded wth 8 qualty layers Drft effect for temporal level 1 Drft effect for temporal level Drft effect temporal level proposed dstorton model measured dstorton 20 proposed dstorton model measured dstorton 60 proposed dstorton model measured dstorton dscarded qualty layers dscarded qualty layers dscarded qualty layers (a) Drft effect on temporal level 1 (b) Drft effect on temporal level 2 (c) Drft effect on temporal level 3 Fg. 4. Analyss of the drft dstorton contrbuton on B-frames at dfferent temporal levels each GOP, dfferent expressons has to be consdered for ntra-frames (key-pctures) and nter-frames (B-pctures). To correctly defne these quanttes n vew of ther use nsde the nteger problem, we need to ntroduce the problem varables as follows. Let x = (x 1,..., x F ) be a vector of bnary varables where x t s gven by: ( ) x t = x t,1 1,..., xt,k 1, x t,1 2,..., xt,k 2,..., x t,1,..., xt,k Each bnary varable x t,k assumes value 1 f the JPEG2000 packet made by the layer k of the precnct P t s consdered n the decodng process and 0 otherwse. Intra-frames are managed as stll pctures, and ths leads to a straghtforward expresson for the dstorton estmaton: ˆD F = D F,K + I =1 k=1 I I (1 x F,k )TF k L F,k (5) where we recall that wth the hrarchcal B-pcture decomposton the key-pcture s always the last pcture of each GOP and consequently t s n poston F, and that D F,K s the key pcture dstorton f all the qualty layers are decoded. The correctness of the expresson (5) has been expermentally valdated, and the results are reported n Fgure 3 whch s referred to a key-pcture encoded wth 8 qualty layers. Fgure 3 compares the dstorton model gven by the equaton (5) to the real dstorton expermentally measured. Fgures 3(a), 3(b) and 3(c) show respectvely the Rate-Dstorton performance for the qualty layer 1,3 and 7. For a partcular layer, each pont n the R-D curve s the contrbuton gven by a precnct,.e. the contrbuton of a packet, to the dstorton reducton of the whole frame. As t can be notced n Fgure 3 the dstorton model gven by equaton (5) s qute accurate, except for the very frst qualty layer. For the purpose of the rate adaptaton, ths dstorton model nconsstency n the lowest qualty layer s neglgble, snce n typcal adaptaton scenaro the lowest qualty layer s always fully consdered n the decodng process, and, as t can be notced from Fgure 3(a), the overall dstorton contrbuton gven by layer 1 s correctly estmated wth the model (5). The expresson of the dstorton for an nter-frame s not easy to derve, snce the herarchcal B-pcture decomposton wth closed-loop moton estmaton ntroduces complex relatonshps between the dstorton on the nter-frames at a partcular level of temporal decomposton and the dstorton on the reference frames of the hgher temporal levels. Brefly, the dstorton on the nter-frames depends on two factors: the dstorton on the moton-compensated resdual related to the nter-frame and the drft effect that depends on the fact that n SNR scalablty dfferent qualty verson of the reference frames could be used for moton compensaton. If the qualty verson that has been used at the encoder for the moton estmaton and compensaton process s not avalable at the decoder, for example because parts of the vdeo btstream are dscarded for rate adaptaton purpose, the drft s ntroduced. The drft problem can be reduced usng a low qualty verson of the reference frames for moton

5 ˆD τ = ntal dstorton y=1 moton-compensated resdual component {}}{{}}{ 2 τ 1 2 τ 1 I D λ,k + (1 x λ,k )Tλ k L λ,k + 2 τ 3 I =1 k=1 y=1 =1 k=1 (1 x F,k )T k F L F,k } {{ } drft component from GOP key-pcture + 2 τ 3 I τ 1 2 z 1 + z=1 w=1 =1 k=1 drft component from lower temporal layers {}}{ [ I 1 ] H 2 τ 2 z (1 x λz,k )Tλ k z L λz,k + (1 x P,k )TP k L P,k }{{} drft component from prevous GOP key-pcture =0 k=1 (6) estmaton and compensaton, but t s well know that ths choce decreases the predcton effcency and consequently the codng performance. Therefore, the common adopted soluton s to use the hghest avalable qualty verson,.e. consderng all the qualty layers, for the reference frames n moton estmaton and compensaton process, and to accept the drft effect. Furthermore, wth the herarchcal B-pcture temporal decomposton, the drft effect also depends on the temporal level. Intutvely, referrng to Fgure 2, t s clear as the drft component of the dstorton ntroduced on the B-level-1 pcture f 4 depends only on the key-pctures f P and f 8, the drft effect ntroduced on the B-level-2 pcture f 2 depends on the the key-pcture f P and the B-level-1 pcture f 4 (that s also affected by f 8 ), and so on. Although we ntentonally avod the mathematcal detals, t can be shown that a good approxmaton of the overall dstorton on the B-level-τ frames, τ = 1,..., log 2 F, s gven by equaton (6), where λ = s + (y 1)δ δ = F 2 τ 1 s = F 2 τ λ z = s z + (w 1)δ z δ z = F 2 z 1 s z = F 2 z It s mportant to note as n expresson (6) t has been consdered that the hghest avalable qualty verson of the reference frames has been used for moton estmaton and compensaton, as prevously descrbed. Ths means that the three drft components are affected by all the mssng qualty layers up to K. Although the proof of the expresson (6) can be found n [11], the expresson has been expermentally valdated, and the results are shown n Fgure 4. Fgure 4 shows the dstorton on the frames that belong to dfferent temporal levels, assumng that the moton-compensated resdual f fully consdered n the decodng process. Consequently, the dstorton s gven by the effect of the three drft components and the ntal dstorton (not equal to 0 snce lossy codng s consdered). In order to show the effect of the drft and ts approxmaton gven by equaton (6), we dscard whole qualty layers of the reference frames n the decodng process, from 0 dscarded layers (no drft) to 5 qualty layers, where the encodng process has been performed generatng 8 qualty layers. As t can be notced n Fgure 4 the dstorton model gven by equaton (6) s accurate assumng to dscard up to 3 qualty layers, a model error of approxmately 10% s evdenced for the fourth dscarded layer whle a greater error s ntroduced from the ffth dscarded layer. Nevertheless, some consderaton has to be done n order to justfy the error ntroduced. Frst, the am of the proposed work s not to gve a quanttatve expresson of the GOP dstorton, but to provde a consstent approxmaton n order to use t for the rate adaptaton model descrbed n the followng secton. Furthermore, the approxmaton error s ntroduced only f we dscard many qualty layers, that means that most of the orgnal vdeo btstream has to be dscarded for rate adaptaton purposes. Nevertheless, the typcal rate adaptaton scenaro that has been consdered n the smulatons performed assumes to dscard from 30% to 50% of the compressed vdeo btstream. It has been verfed that n ths range the proposed dstorton models are accurate. B. The ILP model We now descrbe how modelng the optmal extracton as an ILP problem. The vector x of bnary varables has already been defned. We use bnary varables snce, n order to optmze the decodng performance, we assume to decode only full packets and the stuaton n whch part of packets are forwarded or dscarded s not taken nto account. The objectve functon of the ILP problem s represented by the overall dstorton (2) that, usng the dstorton models ntroduced n the prevous secton, can be approxmated as ˆD T OT = F t=1 log ˆD t = ˆD 2 F F + τ=1 ˆD τ (7) where ˆD F and ˆD τ are respectvely gven by equatons (5) and (6). In the basc ILP model two types of constrants can be dentfed: the qualty layer constrants and the budget constrant. By the progressve layers generaton process of the EBCOT algorthm, the extracton has to verfy the followng condtons: x t,k x t,k+1 0 (8) wth = 1,..., I, k = 1,..., K 1 and t = 1,..., F. As prevously descrbed, the amount of packets that have to be dscarded to perform the adaptaton depends on the avalable bandwdth, that could be converted n terms of bts (L) avalable for each GOP. Consequently, the followng budget constrant has to be consdered: F I t=1 =1 k=1 x t,k L t,k L (9)

6 The exact soluton of the descrbed ILP problem provdes the Rate-Dstorton optmal way to adapt the btstream for each GOP. Addtonally, further constrants could be ntroduced n order to satsfy partcular decodng requrements. For example, n the performed tests the vdeo sequences have been encoded n order to enable near constant decoded qualty f the full vdeo btstream s decoded. Typcally, after the rate adaptaton process ths feature s not mantaned. In order to lmt the decoded qualty fluctuaton potentally ntroduced by the adaptaton, the followng constrants could be ntroduced n the ILP model: β ˆD T OT F ˆD t γ ˆD T OT t = 1,..., F (10) F where β 1 and γ 1. The constrants (10) controls the fluctuaton of the dstorton on the sngle frame wth respect to the mean dstorton over the GOP. Adjustng the values β and γ t s possble to control the level of the dstorton fluctuatons. Nevertheless, t s mportant to note that the ntroducton of the constrant (10) does not guarantee that the ILP problem have a soluton for all the values of β and γ. Furthermore, t s expected that the decodng performance wll be reduced n order to satsfy the constrants. The effect of the constrants (10) on the decodng performance wll be analyzed n secton IV. In contrast to Lnear Programmng (LP) problems, whch can be effcently solved, ILP problems are typcally NPhard thus requrng a computatonal tme whch ncreases exponentally wth the problem sze. The analyzed problem s NP-hard. Nevertheless, t has a partcular structure. The matrx assocated to the qualty layer constrants (8) can be shown to be Totally Unmodular (TUM) (a proof s provded n [11]). ILP problems wth TUM constrant matrx and nteger rghthand-sdes can be solved very effcently, snce the optmal soluton of the related LP problem (obtaned relaxng the constrant that varables x are nteger) corresponds to the optmal nteger soluton. However, t has to be noted that the constrants matrx of our problem s not TUM, snce we have the extra budget constrant (9) and eventually the constrants (10) for the dstorton control. Nevertheless, the constrant (9) s a classcal knapsack constrant, that can be effcently managed by common solvers for mxed nteger lnear programmng problems as CPLEX. Dfferent consderatons has to be done for the constrants (10), that, depends on the values of β and γ, could decrease the model resoluton effcency. IV. EXPERIMENTAL VALIDATION To valdate our method we compare t to the approach proposed n [4]. It s worth notcng that the two approaches are appled on dfferent scalable vdeo codec,.e. SVC and the proposed JPEG2000-based vdeo codec. In order to compare the dfferent codng performance of the two systems Fgures 5(a) and 5(b) show the comparson of the two codecs n sngle layer mode,.e. generatng a non-scalable btstream. As t can be notced, SVC has better compresson performance. A deeper nvestgaton shows that ths dfference depends on the features of the vdeo sequence. Two aspects are evaluated related to the SNR extracton performance: the mean PSNR over the frames obtaned extractng sub-btstream at dfferent data-rate from the full SNR scalable btstream, and the PSNR fluctuatons (evaluated as the PSNR standard devaton) between the frames. The PSNR fluctuatons are evaluated snce the method proposed n [4] operatonally starts by dscardng parts of the full scalable btstream from the lower temporal layers. Ths approach enables to maxmze the mean dstorton over the frames but could ntroduce a fluctuaton of the vdeo qualty that generates annoyng vsual artfacts. As prevously descrbed, due to the dscrepancy between the real dstorton and the proposed dstorton models, the rate adaptaton s performed n order to dscard up to the 50% of the full vdeo btstream. As shown n Fgures 5(c) and 5(d), the proposed extracton model enables comparable performance for Soccer sequence, and better performance for Harbour sequence. However, for both the sequences the codng performance n SNR extracton scenaro ncreases compared to the sngle layer scenaro. Furthermore, Fgures 5(e) and 5(f) show as the proposed approach mantans lmted the PSNR fluctuatons compared to [4], especally at hgher btrates. It has to be noted that at lower btrates the extracton performance of the proposed method decreases compared to SVC. Ths not only depends on the lower codng effcency of our codec compared to SVC, but also by the consderaton that the proposed dstorton models became naccurate when many qualty layers have to be dscarded, as descrbed n secton III-A, generatng a non-optmal extracton. Relatvely to the addtonal constrants (10), n Fgure 6 s shown the effect of the constrants on the extracton performance and on the PSNR fluctuatons, where the dstorton fluctuatons have been controlled settng β = 0 and γ = 1.1. Fgure 6(a) shows how the addtonal constrants decrease the extracton performance of approxmately 0, 5dB, but, as shown n Fgure 6(b), wth a consderable reducton of the PSNR fluctuatons. As a fnal remark, t s mportant to note as the computatonal tme of the proposed adaptaton approach s neglgble. For example, n the test reported n Fgure 5 we used a vdeo sequence at 4CIF resoluton (704x576 pxel), a GOP length equal to 8 frames, and we confgured the JPEG2000 encoder n order to have 5 levels of resoluton, precnct sze equal to pxel and 8 qualty layers. Ths leads to approxmately 3000 JPEG2000 packets (equal to the problem sze) for each GOP. For each GOP the optmal allocaton of the rate, obtaned solvng the ILP problem wth the CPLEX verson 8.1 [12], s performed n fractons of seconds, approxmately 2 or 3 tenths of a second. Ths s manly due the specal structure of the descrbed problem. Smlar computatonal tme can be obtaned also ncreasng the problem sze, for example for vdeo at Hgh Defnton (HD) resoluton. Furthermore, even n the test reported n Fgure 6 when the addtonal constrants (10) have been added to the model, the optmal allocaton s performed n smlar computatonal tme. Ths s manly due to the fact that the dstorton threshold γ = 1.1 s not very

7 σ-y-psnr SNR scalable codng wth adaptaton basc ILP model ILP model wth constrants (10) rate [kb/s] (a) Harbour SNR extracton basc ILP model ILP model wth constrants (10) Y-PSNR standard devaton Verson 4: Sept. 2005, Verson 5 and Verson 6: June 2006, Verson 7: Apr. 2007, Verson 8 (ncludng SVC extenson): Consented n July [2] H. Schwarz, D. Marpe, and T. Wegand, Overvew of the Scalable Vdeo Codng extenson of the H.264/AVC standard, IEEE Transacton on Crcut and System for Vdeo Technology, vol. 17, no. 9, pp , [3] T. Thang, J.-G. Km, J. W. Kang, and J.-J. Yoo, SVC adaptaton: Standard tools and supportng methods, Sgnal Processng: Image Communcaton, [4] I. Amonou, N. Cammas, S. Kervadec, and S. Pateux, Optmzed rate-dstorton extracton wth qualty layers n the Scalable extenson of H.264/AVC, IEEE Transacton on Crcut and System for Vdeo Technology, vol. 17, no. 9, pp , [5] J. Sun, G. Wen, D. Zhao, and W. L, On Rate-Dstorton Modellng and Extracton of H.264/SVC Fne-Granular Scalable Vdeo, IEEE Transacton on Crcut and System for Vdeo Technology, vol. 19, no. 3, pp , [6] ITU-T, JSVM 10 software, jvt-w203, Aprl [7] H. Schwarz, D. Marpe, and T. Wegand, Comparson of MCTF and closed-loop herarchcal B pcs, jvt nput document P059, July [8] D. Taubman, Hgh performance scalable mage compresson wth EBCOT, IEEE Transacton on Image Processng, vol. 9, no. 7, pp , July [9] C. Chrstopoulos, A. Skodras, and T. Ebrahm, The JPEG2000 stll mage codng system: An overvew, IEEE Transacton on Consumer Electroncs, vol. 46, no. 1, pp , November [10] D. Taubman and M. Marcelln, JPEG2000: Image Compresson Fundamentals, Standards and Practce. Kluwer Academc, Boston, MA, USA, [11] L. Lma, Scalablty of vsual nformaton for mproved communcaton, PhD thess, Shortly avalable at lvo.lma/. [12] I. INC., ILOG CPLEX 8.1 reference manual, Inclne Vllage: ILOG Inc., CPLEX Dv., (b) Harbour Y-PSNR standard devaton Fg. 6. Analyss of the effect of the constrants (10) restrctve. Decreasng ts value ncreases the computatonal tme requred to solve the model. V. CONCLUSIONS In ths work we propose an effcent method for SNR scalable vdeo adaptaton based on the formulaton of the adaptaton problem as an Integer Lnear Programmng problem and successfully appled to a JPEG2000-based scalable vdeo codec. The proposed approach shows two very nterestng features. Frst, t provdes a comparable performance wth respect to the adaptaton method used for SVC. Secondly, the TUM property of part of the constrants matrx of the proposed ILP problem can be exploted to fnd effcent approxmated solutons (for nstance, by means of Lagrangan Relaxaton) to more complex adaptaton problems where addtonal constrants such as a further control on the dstorton fluctuaton are ntroduced. REFERENCES [1] Advanced vdeo codng for generc audovsual servces, ITU-T Rec. H.264 and ISO/IEC (MPEG-4 AVC), ITU-T and ISO/IEC JTC 1, Verson 1: May 2003, Verson 2: May 2004, Verson 3: Mar. 2005,

8 42 Sngle layer codng 42 Sngle layer codng SVC proposed codec 36 SVC proposed codec rate [Kb/s] rate [Kb/s] (a) Soccer sngle layer (b) Harbour sngle layer SNR scalable codng wth adaptaton SNR scalable codng wth adaptaton proposed extracton SVC MGS wth QL 39 proposed extracton SVC MGS wth QL rate [kb/s] (c) Soccer SNR extracton rate [kb/s] (d) Harbour SNR extracton 4 Y-PSNR standard devaton 4 Y-PSNR standard devaton 3 3 σ-y-psnr 2 σ-y-psnr 2 1 proposed extracton SVC MGS wth QL 1 proposed extracton SVC MGS wth QL (e) Soccer Y-PSNR standard devaton (f) Harbour Y-PSNR standard devaton Fg. 5. Analyss of the extracton performance

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Novel Pre-Compression Rate-Distortion Optimization Algorithm for JPEG 2000

Novel Pre-Compression Rate-Distortion Optimization Algorithm for JPEG 2000 Novel Pre-Compresson Rate-Dstorton Optmzaton Algorthm for JPEG 2000 Yu-We Chang, Hung-Ch Fang, Chung-Jr Lan, and Lang-Gee Chen DSP/IC Desgn Laboratory, Graduate Insttute of Electroncs Engneerng Natonal

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

Smooth Extraction of SVC Fine-Granular SNR Scalable Videos with a Virtual-GOP-Based Rate-Distortion Modeling

Smooth Extraction of SVC Fine-Granular SNR Scalable Videos with a Virtual-GOP-Based Rate-Distortion Modeling F Smooth Extracton of SVC Fne-Granular SNR Scalable Vdeos wth a Vrtual-GOP-Based Rate-Dstorton Modelng * Jun a SunF, Wen Gao a, Debn Zhao b a Insttute of Dgtal Meda, EECS, Pekng Unversty, Bejng, 1871,

More information

Lossy Compression. Compromise accuracy of reconstruction for increased compression.

Lossy Compression. Compromise accuracy of reconstruction for increased compression. Lossy Compresson Compromse accuracy of reconstructon for ncreased compresson. The reconstructon s usually vsbly ndstngushable from the orgnal mage. Typcally, one can get up to 0:1 compresson wth almost

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

Microwave Diversity Imaging Compression Using Bioinspired

Microwave Diversity Imaging Compression Using Bioinspired Mcrowave Dversty Imagng Compresson Usng Bonspred Neural Networks Youwe Yuan 1, Yong L 1, Wele Xu 1, Janghong Yu * 1 School of Computer Scence and Technology, Hangzhou Danz Unversty, Hangzhou, Zhejang,

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

VQ widely used in coding speech, image, and video

VQ widely used in coding speech, image, and video at Scalar quantzers are specal cases of vector quantzers (VQ): they are constraned to look at one sample at a tme (memoryless) VQ does not have such constrant better RD perfomance expected Source codng

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

ECE559VV Project Report

ECE559VV Project Report ECE559VV Project Report (Supplementary Notes Loc Xuan Bu I. MAX SUM-RATE SCHEDULING: THE UPLINK CASE We have seen (n the presentaton that, for downlnk (broadcast channels, the strategy maxmzng the sum-rate

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

Calculation of time complexity (3%)

Calculation of time complexity (3%) Problem 1. (30%) Calculaton of tme complexty (3%) Gven n ctes, usng exhaust search to see every result takes O(n!). Calculaton of tme needed to solve the problem (2%) 40 ctes:40! dfferent tours 40 add

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING INTRODUCTION

CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING INTRODUCTION CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING N. Phanthuna 1,2, F. Cheevasuvt 2 and S. Chtwong 2 1 Department of Electrcal Engneerng, Faculty of Engneerng Rajamangala

More information

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing Advanced Scence and Technology Letters, pp.164-168 http://dx.do.org/10.14257/astl.2013 Pop-Clc Nose Detecton Usng Inter-Frame Correlaton for Improved Portable Audtory Sensng Dong Yun Lee, Kwang Myung Jeon,

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud Resource Allocaton wth a Budget Constrant for Computng Independent Tasks n the Cloud Wemng Sh and Bo Hong School of Electrcal and Computer Engneerng Georga Insttute of Technology, USA 2nd IEEE Internatonal

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals Smultaneous Optmzaton of Berth Allocaton, Quay Crane Assgnment and Quay Crane Schedulng Problems n Contaner Termnals Necat Aras, Yavuz Türkoğulları, Z. Caner Taşkın, Kuban Altınel Abstract In ths work,

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased

More information

Low Complexity Soft-Input Soft-Output Hamming Decoder

Low Complexity Soft-Input Soft-Output Hamming Decoder Low Complexty Soft-Input Soft-Output Hammng Der Benjamn Müller, Martn Holters, Udo Zölzer Helmut Schmdt Unversty Unversty of the Federal Armed Forces Department of Sgnal Processng and Communcatons Holstenhofweg

More information

Linear Regression Analysis: Terminology and Notation

Linear Regression Analysis: Terminology and Notation ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented

More information

Comparison of Regression Lines

Comparison of Regression Lines STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence

More information

On balancing multiple video streams with distributed QoS control in mobile communications

On balancing multiple video streams with distributed QoS control in mobile communications On balancng multple vdeo streams wth dstrbuted QoS control n moble communcatons Arjen van der Schaaf, José Angel Lso Arellano, and R. (Inald) L. Lagendjk TU Delft, Mekelweg 4, 68 CD Delft, The Netherlands

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Some modelling aspects for the Matlab implementation of MMA

Some modelling aspects for the Matlab implementation of MMA Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton

More information

Pulse Coded Modulation

Pulse Coded Modulation Pulse Coded Modulaton PCM (Pulse Coded Modulaton) s a voce codng technque defned by the ITU-T G.711 standard and t s used n dgtal telephony to encode the voce sgnal. The frst step n the analog to dgtal

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Negative Binomial Regression

Negative Binomial Regression STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...

More information

Homework Assignment 3 Due in class, Thursday October 15

Homework Assignment 3 Due in class, Thursday October 15 Homework Assgnment 3 Due n class, Thursday October 15 SDS 383C Statstcal Modelng I 1 Rdge regresson and Lasso 1. Get the Prostrate cancer data from http://statweb.stanford.edu/~tbs/elemstatlearn/ datasets/prostate.data.

More information

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors Multple Lnear and Polynomal Regresson wth Statstcal Analyss Gven a set of data of measured (or observed) values of a dependent varable: y versus n ndependent varables x 1, x, x n, multple lnear regresson

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

On a direct solver for linear least squares problems

On a direct solver for linear least squares problems ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

More information

Tornado and Luby Transform Codes. Ashish Khisti Presentation October 22, 2003

Tornado and Luby Transform Codes. Ashish Khisti Presentation October 22, 2003 Tornado and Luby Transform Codes Ashsh Khst 6.454 Presentaton October 22, 2003 Background: Erasure Channel Elas[956] studed the Erasure Channel β x x β β x 2 m x 2 k? Capacty of Noseless Erasure Channel

More information

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

More information

General theory of fuzzy connectedness segmentations: reconciliation of two tracks of FC theory

General theory of fuzzy connectedness segmentations: reconciliation of two tracks of FC theory General theory of fuzzy connectedness segmentatons: reconclaton of two tracks of FC theory Krzysztof Chrs Ceselsk Department of Mathematcs, West Vrgna Unversty and MIPG, Department of Radology, Unversty

More information

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem Appled Mathematcal Scences Vol 5 0 no 65 3 33 Interactve B-Level Mult-Objectve Integer Non-lnear Programmng Problem O E Emam Department of Informaton Systems aculty of Computer Scence and nformaton Helwan

More information

Tutorial 2. COMP4134 Biometrics Authentication. February 9, Jun Xu, Teaching Asistant

Tutorial 2. COMP4134 Biometrics Authentication. February 9, Jun Xu, Teaching Asistant Tutoral 2 COMP434 ometrcs uthentcaton Jun Xu, Teachng sstant csjunxu@comp.polyu.edu.hk February 9, 207 Table of Contents Problems Problem : nswer the questons Problem 2: Power law functon Problem 3: Convoluton

More information

x = , so that calculated

x = , so that calculated Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to

More information

Lecture 4: Universal Hash Functions/Streaming Cont d

Lecture 4: Universal Hash Functions/Streaming Cont d CSE 5: Desgn and Analyss of Algorthms I Sprng 06 Lecture 4: Unversal Hash Functons/Streamng Cont d Lecturer: Shayan Oves Gharan Aprl 6th Scrbe: Jacob Schreber Dsclamer: These notes have not been subjected

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information

Entropy Coding. A complete entropy codec, which is an encoder/decoder. pair, consists of the process of encoding or

Entropy Coding. A complete entropy codec, which is an encoder/decoder. pair, consists of the process of encoding or Sgnal Compresson Sgnal Compresson Entropy Codng Entropy codng s also known as zero-error codng, data compresson or lossless compresson. Entropy codng s wdely used n vrtually all popular nternatonal multmeda

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

An Admission Control Algorithm in Cloud Computing Systems

An Admission Control Algorithm in Cloud Computing Systems An Admsson Control Algorthm n Cloud Computng Systems Authors: Frank Yeong-Sung Ln Department of Informaton Management Natonal Tawan Unversty Tape, Tawan, R.O.C. ysln@m.ntu.edu.tw Yngje Lan Management Scence

More information

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification E395 - Pattern Recognton Solutons to Introducton to Pattern Recognton, Chapter : Bayesan pattern classfcaton Preface Ths document s a soluton manual for selected exercses from Introducton to Pattern Recognton

More information

Number of cases Number of factors Number of covariates Number of levels of factor i. Value of the dependent variable for case k

Number of cases Number of factors Number of covariates Number of levels of factor i. Value of the dependent variable for case k ANOVA Model and Matrx Computatons Notaton The followng notaton s used throughout ths chapter unless otherwse stated: N F CN Y Z j w W Number of cases Number of factors Number of covarates Number of levels

More information

TOPICS MULTIPLIERLESS FILTER DESIGN ELEMENTARY SCHOOL ALGORITHM MULTIPLICATION

TOPICS MULTIPLIERLESS FILTER DESIGN ELEMENTARY SCHOOL ALGORITHM MULTIPLICATION 1 2 MULTIPLIERLESS FILTER DESIGN Realzaton of flters wthout full-fledged multplers Some sldes based on support materal by W. Wolf for hs book Modern VLSI Desgn, 3 rd edton. Partly based on followng papers:

More information

On the Multicriteria Integer Network Flow Problem

On the Multicriteria Integer Network Flow Problem BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 5, No 2 Sofa 2005 On the Multcrtera Integer Network Flow Problem Vassl Vasslev, Marana Nkolova, Maryana Vassleva Insttute of

More information

Numerical Solution of Ordinary Differential Equations

Numerical Solution of Ordinary Differential Equations Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples

More information

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE- 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com

More information

Grover s Algorithm + Quantum Zeno Effect + Vaidman

Grover s Algorithm + Quantum Zeno Effect + Vaidman Grover s Algorthm + Quantum Zeno Effect + Vadman CS 294-2 Bomb 10/12/04 Fall 2004 Lecture 11 Grover s algorthm Recall that Grover s algorthm for searchng over a space of sze wors as follows: consder the

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA 4 Analyss of Varance (ANOVA) 5 ANOVA 51 Introducton ANOVA ANOVA s a way to estmate and test the means of multple populatons We wll start wth one-way ANOVA If the populatons ncluded n the study are selected

More information

An Interactive Optimisation Tool for Allocation Problems

An Interactive Optimisation Tool for Allocation Problems An Interactve Optmsaton ool for Allocaton Problems Fredr Bonäs, Joam Westerlund and apo Westerlund Process Desgn Laboratory, Faculty of echnology, Åbo Aadem Unversty, uru 20500, Fnland hs paper presents

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Lec 12 Rate-Distortion Optimization (RDO) in Video Coding-II

Lec 12 Rate-Distortion Optimization (RDO) in Video Coding-II Sprng 07: Multmeda Communcaton Lec ate-dstorton Optmzaton (DO) n Vdeo Codng-II Zhu L Course Web: http://l.web.umkc.edu/lzhu/ Z. L Multmeda Communcaton, Sprng 07 p. Outlne Lec ecap Lagrangan Method HW-3

More information

Annexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances

Annexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances ec Annexes Ths Annex frst llustrates a cycle-based move n the dynamc-block generaton tabu search. It then dsplays the characterstcs of the nstance sets, followed by detaled results of the parametercalbraton

More information

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M CIS56: achne Learnng Lecture 3 (Sept 6, 003) Preparaton help: Xaoyng Huang Lnear Regresson Lnear regresson can be represented by a functonal form: f(; θ) = θ 0 0 +θ + + θ = θ = 0 ote: 0 s a dummy attrbute

More information

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

A FAST HEURISTIC FOR TASKS ASSIGNMENT IN MANYCORE SYSTEMS WITH VOLTAGE-FREQUENCY ISLANDS

A FAST HEURISTIC FOR TASKS ASSIGNMENT IN MANYCORE SYSTEMS WITH VOLTAGE-FREQUENCY ISLANDS Shervn Haamn A FAST HEURISTIC FOR TASKS ASSIGNMENT IN MANYCORE SYSTEMS WITH VOLTAGE-FREQUENCY ISLANDS INTRODUCTION Increasng computatons n applcatons has led to faster processng. o Use more cores n a chp

More information

Linear Feature Engineering 11

Linear Feature Engineering 11 Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19

More information

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan Wnter 2008 CS567 Stochastc Lnear/Integer Programmng Guest Lecturer: Xu, Huan Class 2: More Modelng Examples 1 Capacty Expanson Capacty expanson models optmal choces of the tmng and levels of nvestments

More information

Supporting Information

Supporting Information Supportng Informaton The neural network f n Eq. 1 s gven by: f x l = ReLU W atom x l + b atom, 2 where ReLU s the element-wse rectfed lnear unt, 21.e., ReLUx = max0, x, W atom R d d s the weght matrx to

More information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011 Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected

More information

Chapter - 2. Distribution System Power Flow Analysis

Chapter - 2. Distribution System Power Flow Analysis Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load

More information

Direct Inter-Mode Selection for H.264 Video Coding using Phase Correlation

Direct Inter-Mode Selection for H.264 Video Coding using Phase Correlation Drect Inter-Mode Selecton for H.64 Vdeo Codng usng Phase Correlaton Manoranjan Paul, Member, IEEE, Wes Ln, Senor member, IEEE, Chew Tong Lau, Member, IEEE, and Bu-sung Lee, Member, IEEE Abstract The H.64

More information

Tracking with Kalman Filter

Tracking with Kalman Filter Trackng wth Kalman Flter Scott T. Acton Vrgna Image and Vdeo Analyss (VIVA), Charles L. Brown Department of Electrcal and Computer Engneerng Department of Bomedcal Engneerng Unversty of Vrgna, Charlottesvlle,

More information

Single-Facility Scheduling over Long Time Horizons by Logic-based Benders Decomposition

Single-Facility Scheduling over Long Time Horizons by Logic-based Benders Decomposition Sngle-Faclty Schedulng over Long Tme Horzons by Logc-based Benders Decomposton Elvn Coban and J. N. Hooker Tepper School of Busness, Carnege Mellon Unversty ecoban@andrew.cmu.edu, john@hooker.tepper.cmu.edu

More information

Conversion Scheme for DCT-domain Transcoding of MPEG-2 to H.264/AVC

Conversion Scheme for DCT-domain Transcoding of MPEG-2 to H.264/AVC Converson Scheme for DC-doman ranscodng of MPEG-2 to H.26/AVC Joo-kyong Lee 1, K-dong Chung 1 1, Dept. of Computer Engneerng, Pusan Natonal Unv., Jangjeon-dong, Geumjeong-gu, Busan, 609-735, Korea jklee@melon.cs.pusan.ac.kr,

More information

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS HCMC Unversty of Pedagogy Thong Nguyen Huu et al. A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS Thong Nguyen Huu and Hao Tran Van Department of mathematcs-nformaton,

More information

ADVANCED MACHINE LEARNING ADVANCED MACHINE LEARNING

ADVANCED MACHINE LEARNING ADVANCED MACHINE LEARNING 1 ADVANCED ACHINE LEARNING ADVANCED ACHINE LEARNING Non-lnear regresson technques 2 ADVANCED ACHINE LEARNING Regresson: Prncple N ap N-dm. nput x to a contnuous output y. Learn a functon of the type: N

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

Week 5: Neural Networks

Week 5: Neural Networks Week 5: Neural Networks Instructor: Sergey Levne Neural Networks Summary In the prevous lecture, we saw how we can construct neural networks by extendng logstc regresson. Neural networks consst of multple

More information

The Expectation-Maximization Algorithm

The Expectation-Maximization Algorithm The Expectaton-Maxmaton Algorthm Charles Elan elan@cs.ucsd.edu November 16, 2007 Ths chapter explans the EM algorthm at multple levels of generalty. Secton 1 gves the standard hgh-level verson of the algorthm.

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

Lecture 5 Decoding Binary BCH Codes

Lecture 5 Decoding Binary BCH Codes Lecture 5 Decodng Bnary BCH Codes In ths class, we wll ntroduce dfferent methods for decodng BCH codes 51 Decodng the [15, 7, 5] 2 -BCH Code Consder the [15, 7, 5] 2 -code C we ntroduced n the last lecture

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

S Advanced Digital Communication (4 cr) Targets today

S Advanced Digital Communication (4 cr) Targets today S.72-3320 Advanced Dtal Communcaton (4 cr) Convolutonal Codes Tarets today Why to apply convolutonal codn? Defnn convolutonal codes Practcal encodn crcuts Defnn qualty of convolutonal codes Decodn prncples

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6 Department of Quanttatve Methods & Informaton Systems Tme Seres and Ther Components QMIS 30 Chapter 6 Fall 00 Dr. Mohammad Zanal These sldes were modfed from ther orgnal source for educatonal purpose only.

More information

4DVAR, according to the name, is a four-dimensional variational method.

4DVAR, according to the name, is a four-dimensional variational method. 4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The

More information

Lecture 2: Gram-Schmidt Vectors and the LLL Algorithm

Lecture 2: Gram-Schmidt Vectors and the LLL Algorithm NYU, Fall 2016 Lattces Mn Course Lecture 2: Gram-Schmdt Vectors and the LLL Algorthm Lecturer: Noah Stephens-Davdowtz 2.1 The Shortest Vector Problem In our last lecture, we consdered short solutons to

More information

A LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS. Dr. Derald E. Wentzien, Wesley College, (302) ,

A LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS. Dr. Derald E. Wentzien, Wesley College, (302) , A LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS Dr. Derald E. Wentzen, Wesley College, (302) 736-2574, wentzde@wesley.edu ABSTRACT A lnear programmng model s developed and used to compare

More information

x i1 =1 for all i (the constant ).

x i1 =1 for all i (the constant ). Chapter 5 The Multple Regresson Model Consder an economc model where the dependent varable s a functon of K explanatory varables. The economc model has the form: y = f ( x,x,..., ) xk Approxmate ths by

More information

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm Desgn and Optmzaton of Fuzzy Controller for Inverse Pendulum System Usng Genetc Algorthm H. Mehraban A. Ashoor Unversty of Tehran Unversty of Tehran h.mehraban@ece.ut.ac.r a.ashoor@ece.ut.ac.r Abstract:

More information