Fully smoothed L1-TV energies : properties of the optimal solutions, fast minimization and applications. Mila Nikolova
|
|
- Jane Laurel Atkinson
- 5 years ago
- Views:
Transcription
1 Fully smoothed L1-TV energies : properties of the optimal solutions, fast minimization and applications Mila Nikolova CMLA, ENS Cachan, CNRS 61 Av. President Wilson, F Cachan, FRANCE nikolova@cmla.ens-cachan.fr PGMO group MAthematics and Optimization for Imaging (MAORI) ENSTA Paristech, October 4th 2013
2 A fully smoothed l 1 -TV model [Nikolova, Wen and Chan, JMIV 2012] J(u, f) := Ψ(u, f) + βφ(u), β > 0 Ψ(u, f) := i I n ψ(u[i] f[i]) and Φ(u) := i I r φ(g i u) {g i R 1 n : i I r } forward discretization, options: N4 Only vertical and horizontal differences; N8 Diagonal differences are added. G := [ g T 1,, g T r ] T R r n. ψ( ) := ψ(, α 1 ) : R R and φ( ) := φ(, α 2 ) : R R where (α 1, α 2 ) > 0. They belong to the family of functions θ(, α) : R R satisfying (we set θ (t, α) := d θ(t, α)): dt H 1 For any α > 0 fixed, θ(, α) is C s 2 -continuous, even and θ (t, α) > 0, t R. H 2 For any α > 0 fixed, θ (t, α) < 1 and for t > 0 fixed, it is strictly decreasing in α > 0 with α > 0 lim t θ (t, α) = 1, t R lim α 0 θ (t, α) = 1 and lim α θ (t, α) = 0. 2
3 f1 f2 f3 θ θ θ t t2 + α t2 + α ( α log cosh t α log ( )) t α ( 1 + t α ) ( ) t 1 tanh α α t α + t ( 1 α ( t2 + α ) 3 ( ( )) ) 2 t tanh α α (α + t ) 2 Table 1: Relevant choices for θ(, α) obeying H1 and H θ(t) = t 2 + α θ (t) = t t 2 +α ξ(y) = (θ ) 1 (y) = y α 1 y 2 Assumptions H1 and H2. Function f1: plots for α = 0.05 ( ) and for α = 0.5 ( ). 3
4 Rationale for the choice of J: û tends to decrease the quantization noise Real-valued original f quantized on {0,, 15} Restored f In general fully smoothed L1 - TV are slow to minimize with high precision... 4
5 3. Preliminary facts on J The properties below, derived in [Nikolova, Wen, Chan 12] and [Bauss, Nikolova, Steidl, JIMV 2013] using H1 and H2, underly the proposed algorithm. For any β > 0 and any f, J(, f) has a unique minimizer û. For any β > 0 and any f K n, where K n is a dense open subset of R n, we have û[i] û[j] i, j I n, i j and û[i] f[i] i I n. (1) For any f R n, û satisfies (1) with a very high probability. There is an inverse function ξ(, α 1, ) := (ψ ) 1 (, α 1 ) : ( 1, 1) R which is odd and C s 1. α 1 ξ(y, α 1 ) is strictly increasing on (0, + ), for any y (0, 1). Set η := G 1. If βη < 1 then û f ξ ( ) βη, α 1. This bound is independent of f. And, û f ξ ( ) βη, α 1 as α2 0. 5
6 Semi-Explicit formula for the minimizer J(û, f) = 0 for any i I n ψ (û[i] f[i]) = β φ (g j û)g j [i] û[i] = f[i] + ξ β φ (g j û)g j [i]. j I r j Ir f1 f2 f3 ξ ξ α α y 1 y 2 ( 1 y 2 ) 3 α 2 ln 1 + y α 1 y 1 y 2 αy 1 y α 1 y Table 2: The inverse function ξ(y, α) = (θ ) 1 (y, α) and its derivative ξ with respect to y for all functions in Table 1. 6
7 4. A fast fixed point (FP) algorithm [Nikolova 13] Initialization: u 0 = f for k = 1, 2, u k+1 = X (u k ), X (u) := f ξ ( β Φ(u) ), with ξ applied componentwise. Theorem 1 Let (α 1, α 2 ) > and β > 0 be chosen so that βη < 1 and Then the iteration X converges. β ξ (βη) φ (0, α 2 ) G T G < 1. (2) Corollary 1 (Nikolova, Steidl 13) Let (α 1, α 2 ) > 0 and β > 0 be such that βη < 1. Then u k f ξ(βη, α 1 ) All iterates live in the same ball 7
8 Some values: If G corresponds only to vertical and horizontal differences, then η = G 1 = 4 and G T G = 8. For ψ and φ given by f1 in Table 1 and α 1 = 0.05, α 2 = 0.3 and β = 0.1 X (u) and û f If G also includes diagonal differences, η = G 1 = 8 and G T G = 16. For ψ and φ given by f1 and α 1 = 0.02, α 2 = 0.4 and β = 0.07 we have X (u) and û f FP algorithms are not reputated for their speed... Here we take advantage of two facts: We know the right initialization; All expressions in the algorithm are explicitly known. 8
9 4. Application: Histogram Specification (HS) f input digital gray value M N image / stored as an n := MN vector f[i] {0,, L 1} i I n := {1,, n} 8-bit image L = 256 Histogram H f [k] := 1 n {f[i] = k : i I n} k I L Target histogram ζ = (ζ[1],, ζ[l]) Goal of HS: convert f into f so that H f = ζ order the pixels in f: i j if f[i] < f[j] i 1 i 2 i ζ[1] }{{} ζ[1] i n ζ[l]+1 i n }{{} ζ[l] Ill-posed problem for digital (quantized) images since n L An issue: obtain a meaningful total strict ordering of all pixels in f Histogram equalization (HE) is a particular case of HS where ζ[k] = n/l k I L 9
10 Histogram Equalization original image HE by histeq HE by sort HE ours
11 Modern sorting algorithms (and a smoothed L1-TV) For any pixel f[i], extract K auxiliary information, a k [i], k I K, from f. Set a 0 := f. Then i j if f[i] f[j] and a k [i] < a k [j] for some k {0,, K}. Local Mean Algorithm (LM) [Coltuc, Bolon, Chassery, IEEE TIP, 2006] Idea: use the local mean intensity of a pixel s neighborhood. If two pixels have the same intensity and their local mean is the same, take a larger neighborhood. Repeat until all pixels are ordered. Edges are smoothed. Authors: 6 iterations are enough. Wavelet Approach (WA) [Wan and Shi IEEE TIP, 2007] Idea: use wavelet coefficients from different subbands to order the pixels. If two pixels have the same intensity and the wavelet coefficients in the finest subband are the same, go to next coarser subband. If wavelet coefficients are different, order according to their absolute value. Repeat until all pixels are ordered. Unstable. Authors: 9 iterations are enough. Variational Approach (VA)- fully smoothed L1-TV [Nikolova, Chan, Steidl, Wen 12, 13] Choose α 1, α 2, β so that f û 0.1 and FP converges. Compute û the minimizer of J Sort the pixels of û. No special storage requirements! LM and WA can not cope with large images, e.g (Canon EOS 5D MkII) storage. 11
12 Comparison with the state-of-the-art sorting algorithms The numerical results in [Nikolova, Wen, Chan 12] have shown that SVA clearly outperforms its main competitors LM and WA in terms of quality and memory needs but not in speed. VA G N4, ψ(t)= t 2 +α 1, φ(t)= t 2 +α 2, α 1 =α 2 =0.05, β =0.1., then f û [Nikolova 13]: for FP algor theithm - stopping rule J(u k, f) [Nikolova, Wen, Chan 12]: Polak-Ribière (PR) CG minimization, stopping rule J(u k, f) 10 6 and maximum iterations 35; LM algorithm [Coltuc, Bolon, Chassery 06] for K = 6 WA algorithm [Wan, Shi, 07] for Haar wavelet and K = 9. A much faster FP algorithm for sorting - [Nikolova, Steidll 13 (report)] 12
13 airplane(512 2 ) baril (512 2 ) car (512 2 ) clock (512 2 ) couple (256 2 ) F-16 (512 2 ) house (512 2 ) Lena (512 2 ) man ( ) raffia (512 2 ) sand (512 2 ) stream (512 2 ) Digital images with gray values in {0,, 255} used to compare the algorithms. 13
14 Fail % CPU SVA LM WA SVA LM WA SVA SVA J 10 7 Image PR FP PR FP airplane baril car clock couple F houseb Lena man raffia sand stream means Fail denotes the percentage of pixels that could not be sorted in a strict way. CPU is in seconds. 14
15 Faithful total strict ordering: VA outperforms by far LM and WA. Memory requirements: VA needs only 1 image to sort the pixels; LM - needs 6 and WA - 9. Extensive tests (contrast compression restoration & HE inversion) [Nikolova, Wen, Chan 12] The proposed FP minimization much faster and precise than PR - CG. LM WA Ours Man - zoom of the residuals in HE Inversion. 15
16 Enhancement of color images [Nikolova 13] Strict ordering for the intensity channel HS for the intensity channel Hue and Range-preserving color assignment Fast algorithm in RGB 16
17 Input image ( ) HS by [Han, Yang, Lee 11] ACE (IPOL) HS - ours 17
18 Input image ( ) HS by [Han, Yang, Lee 11] ACE (IPOL) HS - ours 18
Exact Histogram Specification for Digital Images Using a Variational Approach
Exact Histogram Specification for Digital Images Using a Variational Approach Mila Nikolova, You-Wei Wen, Raymond Chan To cite this version: Mila Nikolova, You-Wei Wen, Raymond Chan. Exact Histogram Specification
More informationAuthor's personal copy
DOI 10.1007/s10851-012-0401-8 Exact Histogram Specification for Digital Images Using a Variational Approach Mila Nikolova You-Wei Wen Raymond Chan Springer Science+Business Media New York 2012 Abstract
More informationContre-examples for Bayesian MAP restoration. Mila Nikolova
Contre-examples for Bayesian MAP restoration Mila Nikolova CMLA ENS de Cachan, 61 av. du Président Wilson, 94235 Cachan cedex (nikolova@cmla.ens-cachan.fr) Obergurgl, September 26 Outline 1. MAP estimators
More informationWavelets and Multiresolution Processing
Wavelets and Multiresolution Processing Wavelets Fourier transform has it basis functions in sinusoids Wavelets based on small waves of varying frequency and limited duration In addition to frequency,
More informationFunctionals for signal and image reconstruction: properties of their minimizers and applications
Functionals for signal and image reconstruction: properties of their minimizers and applications Research report Mila NIKOLOVA Centre de Mathématiques et de Leurs Applications (CMLA) - CNRS-UMR 8536 -
More informationEnergy minimization methods
Handbook of Mathematical Methods in Imaging, Springer 211, 1 st edition Editor: Otmar SCHERZER Chapter 5, p. 138 186 Energy minimization methods Mila NIKOLOVA CMLA, ENS Cachan, CNRS, UniverSud, 61 Av.
More informationOver-enhancement Reduction in Local Histogram Equalization using its Degrees of Freedom. Alireza Avanaki
Over-enhancement Reduction in Local Histogram Equalization using its Degrees of Freedom Alireza Avanaki ABSTRACT A well-known issue of local (adaptive) histogram equalization (LHE) is over-enhancement
More informationTWO-PHASE APPROACH FOR DEBLURRING IMAGES CORRUPTED BY IMPULSE PLUS GAUSSIAN NOISE. Jian-Feng Cai. Raymond H. Chan. Mila Nikolova
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume 00, Number 0, Xxxx XXXX pp. 000 000 TWO-PHASE APPROACH FOR DEBLURRING IMAGES CORRUPTED BY IMPULSE PLUS GAUSSIAN NOISE Jian-Feng
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationarxiv: v1 [cs.oh] 3 Oct 2014
M. Prisheltsev (Voronezh) mikhail.prisheltsev@gmail.com ADAPTIVE TWO-DIMENSIONAL WAVELET TRANSFORMATION BASED ON THE HAAR SYSTEM 1 Introduction arxiv:1410.0705v1 [cs.oh] Oct 2014 The purpose is to study
More informationLeast squares regularized or constrained by L0: relationship between their global minimizers. Mila Nikolova
Least squares regularized or constrained by L0: relationship between their global minimizers Mila Nikolova CMLA, CNRS, ENS Cachan, Université Paris-Saclay, France nikolova@cmla.ens-cachan.fr SIAM Minisymposium
More informationA Randomized Algorithm for the Approximation of Matrices
A Randomized Algorithm for the Approximation of Matrices Per-Gunnar Martinsson, Vladimir Rokhlin, and Mark Tygert Technical Report YALEU/DCS/TR-36 June 29, 2006 Abstract Given an m n matrix A and a positive
More informationModule 4. Multi-Resolution Analysis. Version 2 ECE IIT, Kharagpur
Module 4 Multi-Resolution Analysis Lesson Multi-resolution Analysis: Discrete avelet Transforms Instructional Objectives At the end of this lesson, the students should be able to:. Define Discrete avelet
More informationWavelets and Image Compression. Bradley J. Lucier
Wavelets and Image Compression Bradley J. Lucier Abstract. In this paper we present certain results about the compression of images using wavelets. We concentrate on the simplest case of the Haar decomposition
More information1 The Continuous Wavelet Transform The continuous wavelet transform (CWT) Discretisation of the CWT... 2
Contents 1 The Continuous Wavelet Transform 1 1.1 The continuous wavelet transform (CWT)............. 1 1. Discretisation of the CWT...................... Stationary wavelet transform or redundant wavelet
More informationL. Yaroslavsky. Fundamentals of Digital Image Processing. Course
L. Yaroslavsky. Fundamentals of Digital Image Processing. Course 0555.330 Lec. 6. Principles of image coding The term image coding or image compression refers to processing image digital data aimed at
More informationLecture 2: Haar Multiresolution analysis
WAVELES AND MULIRAE DIGIAL SIGNAL PROCESSING Lecture 2: Haar Multiresolution analysis Prof.V. M. Gadre, EE, II Bombay 1 Introduction HAAR was a mathematician, who has given an idea that any continuous
More informationAssignment #09 - Solution Manual
Assignment #09 - Solution Manual 1. Choose the correct statements about representation of a continuous signal using Haar wavelets. 1.5 points The signal is approximated using sin and cos functions. The
More informationSatellite image deconvolution using complex wavelet packets
Satellite image deconvolution using complex wavelet packets André Jalobeanu, Laure Blanc-Féraud, Josiane Zerubia ARIANA research group INRIA Sophia Antipolis, France CNRS / INRIA / UNSA www.inria.fr/ariana
More informationFast Multilevel Algorithm for a Minimization Problem in Impulse Noise Removal
Fast Multilevel Algorithm for a Minimization Problem in Impulse Noise Removal Raymond H. Chan and Ke Chen September 25, 7 Abstract An effective 2-phase method for removing impulse noise was recently proposed.
More informationInverse Problems in Image Processing
H D Inverse Problems in Image Processing Ramesh Neelamani (Neelsh) Committee: Profs. R. Baraniuk, R. Nowak, M. Orchard, S. Cox June 2003 Inverse Problems Data estimation from inadequate/noisy observations
More informationModeling Multiscale Differential Pixel Statistics
Modeling Multiscale Differential Pixel Statistics David Odom a and Peyman Milanfar a a Electrical Engineering Department, University of California, Santa Cruz CA. 95064 USA ABSTRACT The statistics of natural
More informationInverse modeling in inverse problems using optimization
Summer School 214: Inverse Problem and Image Processing Inverse modeling in inverse problems using optimization Mila Nikolova CMLA CNRS, ENS Cachan, France nikolova@cmla.ens-cachan.fr http://mnikolova.perso.math.cnrs.fr
More informationWavelets Marialuce Graziadei
Wavelets Marialuce Graziadei 1. A brief summary 2. Vanishing moments 3. 2D-wavelets 4. Compression 5. De-noising 1 1. A brief summary φ(t): scaling function. For φ the 2-scale relation hold φ(t) = p k
More informationNegative Association, Ordering and Convergence of Resampling Methods
Negative Association, Ordering and Convergence of Resampling Methods Nicolas Chopin ENSAE, Paristech (Joint work with Mathieu Gerber and Nick Whiteley, University of Bristol) Resampling schemes: Informal
More informationSignal Analysis. Filter Banks and. One application for filter banks is to decompose the input signal into different bands or channels
Filter banks Multi dimensional Signal Analysis A common type of processing unit for discrete signals is a filter bank, where some input signal is filtered by n filters, producing n channels Channel 1 Lecture
More informationNon-negative Quadratic Programming Total Variation Regularization for Poisson Vector-Valued Image Restoration
University of New Mexico UNM Digital Repository Electrical & Computer Engineering Technical Reports Engineering Publications 5-10-2011 Non-negative Quadratic Programming Total Variation Regularization
More informationEBCOT coding passes explained on a detailed example
EBCOT coding passes explained on a detailed example Xavier Delaunay d.xav@free.fr Contents Introduction Example used Coding of the first bit-plane. Cleanup pass............................. Coding of the
More informationSYDE 575: Introduction to Image Processing. Image Compression Part 2: Variable-rate compression
SYDE 575: Introduction to Image Processing Image Compression Part 2: Variable-rate compression Variable-rate Compression: Transform-based compression As mentioned earlier, we wish to transform image data
More informationDefn 3.1: An n-manifold, M, is a topological space with the following properties:
Chapter 1 all sections 1.3 Defn: M is locally Euclidean of dimension n if for all p M, there exists an open set U p such that p U p and there exists a homeomorphism f p : U p V p where V p R n. (U p, f)
More informationWavelets and multiresolution representations. Time meets frequency
Wavelets and multiresolution representations Time meets frequency Time-Frequency resolution Depends on the time-frequency spread of the wavelet atoms Assuming that ψ is centred in t=0 Signal domain + t
More informationImage restoration: edge preserving
Image restoration: edge preserving Convex penalties Jean-François Giovannelli Groupe Signal Image Laboratoire de l Intégration du Matériau au Système Univ. Bordeaux CNRS BINP 1 / 32 Topics Image restoration,
More informationIntroduction - Motivation. Many phenomena (physical, chemical, biological, etc.) are model by differential equations. f f(x + h) f(x) (x) = lim
Introduction - Motivation Many phenomena (physical, chemical, biological, etc.) are model by differential equations. Recall the definition of the derivative of f(x) f f(x + h) f(x) (x) = lim. h 0 h Its
More informationInexact Alternating Direction Method of Multipliers for Separable Convex Optimization
Inexact Alternating Direction Method of Multipliers for Separable Convex Optimization Hongchao Zhang hozhang@math.lsu.edu Department of Mathematics Center for Computation and Technology Louisiana State
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More information6 The SVD Applied to Signal and Image Deblurring
6 The SVD Applied to Signal and Image Deblurring We will discuss the restoration of one-dimensional signals and two-dimensional gray-scale images that have been contaminated by blur and noise. After an
More informationOptimum Notch Filtering - I
Notch Filtering 1 Optimum Notch Filtering - I Methods discussed before: filter too much image information. Optimum: minimizes local variances of the restored image fˆ x First Step: Extract principal freq.
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationMaximum Entropy - Lecture 8
1 Introduction Maximum Entropy - Lecture 8 In the previous discussions we have not been particularly concerned with the selection of the initial prior probability to use in Bayes theorem, as the result
More informationON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX- ON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. HOANG AND A. G. RAMM (Communicated
More informationConvergence of Newton's Method for a Minimization Problem in Impulse Noise Removal Raymond H. Chan, Chung-Wa Ho, and Mila Nikolova y February 5, 2004
Convergence of Newton's Method for a Minimization Problem in Impulse Noise Removal Raymond H. Chan, Chung-Wa Ho, and Mila Nikolova y February 5, 004 Abstract Recently, two-phase schemes for removing salt-and-pepper
More information8 The SVD Applied to Signal and Image Deblurring
8 The SVD Applied to Signal and Image Deblurring We will discuss the restoration of one-dimensional signals and two-dimensional gray-scale images that have been contaminated by blur and noise. After an
More information8 The SVD Applied to Signal and Image Deblurring
8 The SVD Applied to Signal and Image Deblurring We will discuss the restoration of one-dimensional signals and two-dimensional gray-scale images that have been contaminated by blur and noise. After an
More informationHorizon fluff. A semi-classical approach to (BTZ) black hole microstates
Horizon fluff A semi-classical approach to (BTZ) black hole microstates 1705.06257 with D. Grumiller, M.M. Sheikh-Jabbari and H. Yavartanoo Hamid R. Afshar 11 May 2017 - Recent Trends in String Theory
More information1 Directional Derivatives and Differentiability
Wednesday, January 18, 2012 1 Directional Derivatives and Differentiability Let E R N, let f : E R and let x 0 E. Given a direction v R N, let L be the line through x 0 in the direction v, that is, L :=
More informationWavelet Packet Based Digital Image Watermarking
Wavelet Packet Based Digital Image ing A.Adhipathi Reddy, B.N.Chatterji Department of Electronics and Electrical Communication Engg. Indian Institute of Technology, Kharagpur 72 32 {aar, bnc}@ece.iitkgp.ernet.in
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationMagnetic wells in dimension three
Magnetic wells in dimension three Yuri A. Kordyukov joint with Bernard Helffer & Nicolas Raymond & San Vũ Ngọc Magnetic Fields and Semiclassical Analysis Rennes, May 21, 2015 Yuri A. Kordyukov (Ufa) Magnetic
More informationLarge-Scale L1-Related Minimization in Compressive Sensing and Beyond
Large-Scale L1-Related Minimization in Compressive Sensing and Beyond Yin Zhang Department of Computational and Applied Mathematics Rice University, Houston, Texas, U.S.A. Arizona State University March
More informationMidterm Summary Fall 08. Yao Wang Polytechnic University, Brooklyn, NY 11201
Midterm Summary Fall 8 Yao Wang Polytechnic University, Brooklyn, NY 2 Components in Digital Image Processing Output are images Input Image Color Color image image processing Image Image restoration Image
More informationDifferential Attribute Profiles in Remote Sensing
Differential Attribute Profiles in Remote Sensing Georgios K. Ouzounis DigitalGlobe, Inc. 1601 Dry Creek Drive, Suite 260, Longmont, CO 80503, USA. Image: Dubai July 5, 2011 50cm Contents Overview Connected
More informationMultiresolution schemes
Multiresolution schemes Fondamenti di elaborazione del segnale multi-dimensionale Multi-dimensional signal processing Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Elaborazione
More informationFuzzy quantization of Bandlet coefficients for image compression
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 2013, 4(2):140-146 Fuzzy quantization of Bandlet coefficients for image compression R. Rajeswari and R. Rajesh ISSN:
More informationMultiresolution schemes
Multiresolution schemes Fondamenti di elaborazione del segnale multi-dimensionale Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Elaborazione dei Segnali Multi-dimensionali e
More informationManifold Coarse Graining for Online Semi-supervised Learning
for Online Semi-supervised Learning Mehrdad Farajtabar, Amirreza Shaban, Hamid R. Rabiee, Mohammad H. Rohban Digital Media Lab, Department of Computer Engineering, Sharif University of Technology, Tehran,
More informationencoding without prediction) (Server) Quantization: Initial Data 0, 1, 2, Quantized Data 0, 1, 2, 3, 4, 8, 16, 32, 64, 128, 256
General Models for Compression / Decompression -they apply to symbols data, text, and to image but not video 1. Simplest model (Lossless ( encoding without prediction) (server) Signal Encode Transmit (client)
More informationAccelerated Dual Gradient-Based Methods for Total Variation Image Denoising/Deblurring Problems (and other Inverse Problems)
Accelerated Dual Gradient-Based Methods for Total Variation Image Denoising/Deblurring Problems (and other Inverse Problems) Donghwan Kim and Jeffrey A. Fessler EECS Department, University of Michigan
More informationImage and Multidimensional Signal Processing
Image and Multidimensional Signal Processing Professor William Hoff Dept of Electrical Engineering &Computer Science http://inside.mines.edu/~whoff/ Image Compression 2 Image Compression Goal: Reduce amount
More informationHistogram Processing
Histogram Processing The histogram of a digital image with gray levels in the range [0,L-] is a discrete function h ( r k ) = n k where r k n k = k th gray level = number of pixels in the image having
More informationBilevel Image Denoising Using Gaussianity Tests
Bilevel Image Denoising Using Gaussianity Tests Jérôme Fehrenbach, Mila Nikolova, Gabriele Steidl, Pierre Weiss To cite this version: Jérôme Fehrenbach, Mila Nikolova, Gabriele Steidl, Pierre Weiss. Bilevel
More information= F (b) F (a) F (x i ) F (x i+1 ). a x 0 x 1 x n b i
Real Analysis Problem 1. If F : R R is a monotone function, show that F T V ([a,b]) = F (b) F (a) for any interval [a, b], and that F has bounded variation on R if and only if it is bounded. Here F T V
More informationSingular Value Decompsition
Singular Value Decompsition Massoud Malek One of the most useful results from linear algebra, is a matrix decomposition known as the singular value decomposition It has many useful applications in almost
More informationUnconstrained optimization
Chapter 4 Unconstrained optimization An unconstrained optimization problem takes the form min x Rnf(x) (4.1) for a target functional (also called objective function) f : R n R. In this chapter and throughout
More informationMean-field SDE driven by a fractional BM. A related stochastic control problem
Mean-field SDE driven by a fractional BM. A related stochastic control problem Rainer Buckdahn, Université de Bretagne Occidentale, Brest Durham Symposium on Stochastic Analysis, July 1th to July 2th,
More informationProgramming, numerics and optimization
Programming, numerics and optimization Lecture C-3: Unconstrained optimization II Łukasz Jankowski ljank@ippt.pan.pl Institute of Fundamental Technological Research Room 4.32, Phone +22.8261281 ext. 428
More information2D Wavelets. Hints on advanced Concepts
2D Wavelets Hints on advanced Concepts 1 Advanced concepts Wavelet packets Laplacian pyramid Overcomplete bases Discrete wavelet frames (DWF) Algorithme à trous Discrete dyadic wavelet frames (DDWF) Overview
More informationJanuary 29, Non-linear conjugate gradient method(s): Fletcher Reeves Polak Ribière January 29, 2014 Hestenes Stiefel 1 / 13
Non-linear conjugate gradient method(s): Fletcher Reeves Polak Ribière Hestenes Stiefel January 29, 2014 Non-linear conjugate gradient method(s): Fletcher Reeves Polak Ribière January 29, 2014 Hestenes
More informationHalf-Pel Accurate Motion-Compensated Orthogonal Video Transforms
Flierl and Girod: Half-Pel Accurate Motion-Compensated Orthogonal Video Transforms, IEEE DCC, Mar. 007. Half-Pel Accurate Motion-Compensated Orthogonal Video Transforms Markus Flierl and Bernd Girod Max
More informationSet theory and topology
arxiv:1306.6926v1 [math.ho] 28 Jun 2013 Set theory and topology An introduction to the foundations of analysis 1 Part II: Topology Fundamental Felix Nagel Abstract We provide a formal introduction into
More informationAdaptive Corrected Procedure for TVL1 Image Deblurring under Impulsive Noise
Adaptive Corrected Procedure for TVL1 Image Deblurring under Impulsive Noise Minru Bai(x T) College of Mathematics and Econometrics Hunan University Joint work with Xiongjun Zhang, Qianqian Shao June 30,
More informationNumber Representation and Waveform Quantization
1 Number Representation and Waveform Quantization 1 Introduction This lab presents two important concepts for working with digital signals. The first section discusses how numbers are stored in memory.
More informationThe Ordering Code for various standard model Recorders with an AC supply and without any additional options are as follows:
The basic PC software is supplied free with the recorder. There is an additional charge for the extensive Data Acquisition Software supplied with communication of RS-232/422/485 or Ethernet. The Ordering
More information6. Iterative Methods for Linear Systems. The stepwise approach to the solution...
6 Iterative Methods for Linear Systems The stepwise approach to the solution Miriam Mehl: 6 Iterative Methods for Linear Systems The stepwise approach to the solution, January 18, 2013 1 61 Large Sparse
More informationANALYTICAL BOUNDS ON THE MINIMIZERS OF (NONCONVEX) REGULARIZED LEAST-SQUARES. Mila Nikolova. (Communicated by Jackie (Jianhong) Shen)
Inverse Problems and Imaging Volume 1, No. 4, 2007, 661 677 Web site: http://www.aimsciences.org ANALYTICAL BOUNDS ON THE MINIMIZERS OF NONCONVEX) REGULARIZED LEAST-SQUARES Mila Nikolova CMLA, ENS Cachan,
More informationThe Application of Legendre Multiwavelet Functions in Image Compression
Journal of Modern Applied Statistical Methods Volume 5 Issue 2 Article 3 --206 The Application of Legendre Multiwavelet Functions in Image Compression Elham Hashemizadeh Department of Mathematics, Karaj
More informationScalable color image coding with Matching Pursuit
SCHOOL OF ENGINEERING - STI SIGNAL PROCESSING INSTITUTE Rosa M. Figueras i Ventura CH-115 LAUSANNE Telephone: +4121 6935646 Telefax: +4121 69376 e-mail: rosa.figueras@epfl.ch ÉCOLE POLYTECHNIQUE FÉDÉRALE
More informationSolutions to Problem Set 5 for , Fall 2007
Solutions to Problem Set 5 for 18.101, Fall 2007 1 Exercise 1 Solution For the counterexample, let us consider M = (0, + ) and let us take V = on M. x Let W be the vector field on M that is identically
More informationarxiv: v1 [cs.cv] 7 May 2016
FAS BILAERAL FILERIG OF VECOR-VALUED IMAGES Sanay Ghosh and Kunal. Chaudhury Department of Electrical Engineering, Indian Institute of Science. arxiv:1605.02164v1 [cs.cv] 7 May 2016 ABSRAC In this paper,
More informationImage Compression Using the Haar Wavelet Transform
College of the Redwoods http://online.redwoods.cc.ca.us/instruct/darnold/laproj/fall2002/ames/ 1/33 Image Compression Using the Haar Wavelet Transform Greg Ames College of the Redwoods Math 45 Linear Algebra
More informationRegularization Theory
Regularization Theory Solving the inverse problem of Super resolution with CNN Aditya Ganeshan Under the guidance of Dr. Ankik Kumar Giri December 13, 2016 Table of Content 1 Introduction Material coverage
More informationDigital Affine Shear Filter Banks with 2-Layer Structure
Digital Affine Shear Filter Banks with -Layer Structure Zhihua Che and Xiaosheng Zhuang Department of Mathematics, City University of Hong Kong, Tat Chee Avenue, Kowloon Tong, Hong Kong Email: zhihuache-c@my.cityu.edu.hk,
More informationFrom Fourier Series to Analysis of Non-stationary Signals - II
From Fourier Series to Analysis of Non-stationary Signals - II prof. Miroslav Vlcek October 10, 2017 Contents Signals 1 Signals 2 3 4 Contents Signals 1 Signals 2 3 4 Contents Signals 1 Signals 2 3 4 Contents
More informationComputational Linear Algebra
Computational Linear Algebra PD Dr. rer. nat. habil. Ralf-Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2018/19 Part 4: Iterative Methods PD
More informationComputation Theory Finite Automata
Computation Theory Dept. of Computing ITT Dublin October 14, 2010 Computation Theory I 1 We would like a model that captures the general nature of computation Consider two simple problems: 2 Design a program
More informationMultimedia Databases. Previous Lecture. 4.1 Multiresolution Analysis. 4 Shape-based Features. 4.1 Multiresolution Analysis
Previous Lecture Multimedia Databases Texture-Based Image Retrieval Low Level Features Tamura Measure, Random Field Model High-Level Features Fourier-Transform, Wavelets Wolf-Tilo Balke Silviu Homoceanu
More informationTikhonov Regularization of Large Symmetric Problems
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2000; 00:1 11 [Version: 2000/03/22 v1.0] Tihonov Regularization of Large Symmetric Problems D. Calvetti 1, L. Reichel 2 and A. Shuibi
More informationAN INTRODUCTION TO COMPRESSIVE SENSING
AN INTRODUCTION TO COMPRESSIVE SENSING Rodrigo B. Platte School of Mathematical and Statistical Sciences APM/EEE598 Reverse Engineering of Complex Dynamical Networks OUTLINE 1 INTRODUCTION 2 INCOHERENCE
More informationDetecting LSB Matching by Applying Calibration Technique for Difference Image
Detecting LSB atching by Applying Calibration Technique for Difference Image iaolong Li Institute of Computer Science & Technology Peking University 87, Beijing, P. R. China lixiaolong@icst.pku.edu.cn
More informationINTRODUCTION TO. Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR)
INTRODUCTION TO WAVELETS Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR) CRITICISM OF FOURIER SPECTRUM It gives us the spectrum of the
More information<Outline> JPEG 2000 Standard - Overview. Modes of current JPEG. JPEG Part I. JPEG 2000 Standard
JPEG 000 tandard - Overview Ping-ing Tsai, Ph.D. JPEG000 Background & Overview Part I JPEG000 oding ulti-omponent Transform Bit Plane oding (BP) Binary Arithmetic oding (BA) Bit-Rate ontrol odes
More informationSignal Denoising with Wavelets
Signal Denoising with Wavelets Selin Aviyente Department of Electrical and Computer Engineering Michigan State University March 30, 2010 Introduction Assume an additive noise model: x[n] = f [n] + w[n]
More informationMultimedia Databases. Wolf-Tilo Balke Philipp Wille Institut für Informationssysteme Technische Universität Braunschweig
Multimedia Databases Wolf-Tilo Balke Philipp Wille Institut für Informationssysteme Technische Universität Braunschweig http://www.ifis.cs.tu-bs.de 4 Previous Lecture Texture-Based Image Retrieval Low
More informationAPPLICATIONS OF DIFFERENTIABILITY IN R n.
APPLICATIONS OF DIFFERENTIABILITY IN R n. MATANIA BEN-ARTZI April 2015 Functions here are defined on a subset T R n and take values in R m, where m can be smaller, equal or greater than n. The (open) ball
More informationAn Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems
An Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems Justin Romberg Georgia Tech, School of ECE ENS Winter School January 9, 2012 Lyon, France Applied and Computational
More informationImage Processing by the Curvelet Transform
Image Processing by the Curvelet Transform Jean Luc Starck Dapnia/SEDI SAP, CEA Saclay, France. jstarck@cea.fr http://jstarck.free.fr Collaborators: D.L. Donoho, Department of Statistics, Stanford E. Candès,
More informationLinear Diffusion and Image Processing. Outline
Outline Linear Diffusion and Image Processing Fourier Transform Convolution Image Restoration: Linear Filtering Diffusion Processes for Noise Filtering linear scale space theory Gauss-Laplace pyramid for
More informationDesign and Implementation of Multistage Vector Quantization Algorithm of Image compression assistant by Multiwavelet Transform
Design and Implementation of Multistage Vector Quantization Algorithm of Image compression assistant by Multiwavelet Transform Assist Instructor/ BASHAR TALIB HAMEED DIYALA UNIVERSITY, COLLEGE OF SCIENCE
More informationarxiv: v1 [cs.it] 21 Feb 2013
q-ary Compressive Sensing arxiv:30.568v [cs.it] Feb 03 Youssef Mroueh,, Lorenzo Rosasco, CBCL, CSAIL, Massachusetts Institute of Technology LCSL, Istituto Italiano di Tecnologia and IIT@MIT lab, Istituto
More informationMultimedia Databases. 4 Shape-based Features. 4.1 Multiresolution Analysis. 4.1 Multiresolution Analysis. 4.1 Multiresolution Analysis
4 Shape-based Features Multimedia Databases Wolf-Tilo Balke Silviu Homoceanu Institut für Informationssysteme Technische Universität Braunschweig http://www.ifis.cs.tu-bs.de 4 Multiresolution Analysis
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More information