Super-Resolution for Hyperspectral Remote Sensing Images with Parallel Level Sets
|
|
- Curtis Johnston
- 6 years ago
- Views:
Transcription
1 Super-Resolution for Hyperspectral Remote Sensing Images with Parallel Level Sets Leon Bungert 1, Matthias J. Ehrhardt 2, Marc Aurele Gilles 3, Jennifer Rasch 4, Rafael Reisenhofer 5 (1) Applied Mathematics 1, Friedrich-Alexander Universität Erlangen-Nürnberg (2) Department for Applied Mathematics and Theoretical Physics, University of Cambridge (3) Center for Applied Mathematics, Cornell University (4) Fraunhofer Heinrich Hertz Institute, Berlin (5) AG Computational Data Analysis, Universität Bremen AIP 2017, Zhejiang University 31 May 2017
2 Outline Problem Model Algorithmic Realization Numerical Results
3 Two Imaging Modalities low-res hyperspectral image (126 channels, px) channel 20 channel 47 channel 63 channel 90 high-res RGB image ( px) Eventually, biologists aim to correctly segment different kinds of trees. 3/31
4 Basic Approach Goal: Increase the resolution of the hyperspectral image to the resolution of the RGB image. Assumption: Structures that are visible in the hyperspectral image are also visible in the RGB image. We express the structural information of the RGB image in terms of level sets or gradients and use this as additional a-priori knowledge in the super-resolution process. 4/31
5 Model We consider solutions (u, k ) of the minimization problem where 1 argmin (u,k) U K 2 S C k u f 2 + α dtv v (u) + β TV(k), N l, N h, N k N denote numbers of pixels, U = [0, 1] N h, K is the N k -dimensional unit simplex, C k is the convolution operator with respect to a kernel k, S is a subsampling operator, f R N l is a single channel of the hyperspectral image, dtv v denotes directional total variation penalty with respect to the vector-field v, TV denotes the total variation penalty and α, β > 0 are regularization parameters. 5/31
6 Model - Forward Operator argmin (u,k ) U K 1 S Ck u f 2 + α dtvv (u ) + β TV(k ) 2 Cyclic convolution operator with a kernel k K Ck : RNh RNh, F 2D (Ck u )n = F 2D (u )n F 2D (k )n Subsampling operator S : RNh RNl, (S u )n = uπ [n ] with a mapping π {1,..., Nh }Nl. high-res channel u Ck (u ) S Ck (u ) 6/31
7 Model - Directional Total Variation 1 argmin (u,k) U K 2 S C k u f 2 + αdtv v (u) + β TV(k) Let : R N h R N h 2 denote the discrete gradient operator based on forward differences. For a paremter η > 0 and a high-res image u R N h, we consider u the vectorfield v = u. 2 +η 2 image u vectorfield v Ehrhardt & Betcke, /31
8 Model - Directional Total Variation 1 argmin (u,k) U K 2 S C k u f 2 + αdtv v (u) + β TV(k) Let a vectorfield v be given as in the previous slide and denote P vn = I v n v t n where I is the 2 2 identity matrix, then the directional TV penalty is given by dtv v (u) = N h N h P vn u n = u n v n, u n v n. n=1 n=1 Ehrhardt & Betcke, /31
9 Model - Regularization of the Kernel 1 argmin (u,k) U K 2 S C k u f 2 + α dtv v (u) + βtv(k) We assume that the kernel k has a sparse gradient and include the total variation penalty TV(k) = N h k n. n=1 Additionally, we assume that the kernel satisfies N k k n = 1 and k n 0, n N k. n=1 Rudin, Osher & Fatemi, /31
10 Model 1 argmin (u,k) U K 2 S C k u f 2 + α dtv v (u) + β TV(k) The kernel k is mainly associated to the PSF of the imaging device but can also be used to model motion blur stemming from the movement of the aircraft and to resolve translational differences between the hyperspectral image and the RGB image. The data fidelty term, dtv v as well as TV are convex and the energy functional is convex in u and convex in k. However, the energy functional is not convex in (u, k). 10/31
11 Algorithmic Realization We apply inertial alternating linearized minimization (ipalm; Pock & Sabach, 2016) which was recently proposed for a general class of nonconvex and nonsmooth problems. Guarantees convergence to a critical point for problems of the form argmin H(x, y) + f 1 (x) + f 2 (y), (x,y) R N 1 R N 2 where H(x, y) is smooth and its gradient Lipschitz continuous in both x and y and f 1 and f 2 can be general nonsmooth and nonconvex functions albeit with efficiently computable proximal mappings. 1 argmin (u,k) U K 2 S C k u f 2 }{{} H(u,k) + α dtv v (u) }{{} f 1 (u) + β TV(k) }{{} f 2 (k) 11/31
12 Algorithmic Realization ipalm includes inertial term in proximal alternating linearized minimization (PALM; Bolte, Sabach & Teboulle, 2014) The update step j j + 1 for the variable u is given by ũ = u (j ) + α (j ) (u (j ) u (j 1) ) ŭ = u (j ) + β (j ) (u (j ) u (j 1) ) ( ) ũ 1 u S C k (j ) ŭ f 2 u (j +1) prox dtv v τ (j ) 2τ (j ) 1 where α (j ) and β (j ) are inertial parameters and τ (j ) controls the step size and should be chosen proportional to the Lipschitz constant of u S C k (j ) ŭ f 2. The update step for the variable k is given analogously. Pock & Sabach, /31
13 Algorithmic Realization In order to apply the update step ũ = u (j ) + α (j ) (u (j ) u (j 1) ) ŭ = u (j ) + β (j ) (u (j ) u (j 1) ) ( ) ũ 1 u S C k (j ) ŭ f 2, u (j +1) prox dtv v τ (j ) We need to compute 2τ (j ) 1 prox dtv v τ (j ) (and prox TV τ (j ) for k) u S C k (j ) ŭ f 2 (and k S C k (j ) ŭ f 2 ) and the Lipschitz constant of u S C k (j ) ŭ f 2. 13/31
14 Algorithmic Realization u (j +1) prox dtv v τ (j ) ( ũ 1 2τ (j ) u S Ck (j ) ŭ f ) 2 1 prox dtv v and prox TV can by computed via fast gradient τ (j ) τ (j ) projection (Beck & Teboulle, 2009) u S C k (j ) ŭ f 2 = 2 C k S (S C k u f )) Let L be the Lipschitz constant of u S C k (j ) ŭ f 2 then it holds that L max F 2D 2 (k) n n Alternatively, a backtracking scheme can be applied to obtain smaller values for τ (j ) and thereby larger step sizes. 14/31
15 Numerical Results - Test Case 1 low-res hyperspectral channel ( px) high-res RGB image ( px) 15/31
16 Numerical Results - Test Case 1 input f argmin (u,k ) U K vectorfield v (η = 0.001) 1 S Ck u f 2 + α dtvv (u ) + β TV(k ), 2 with α = 0.5, β = 10 16/31
17 Numerical Results - Test Case 1 (u, k ) after 5000 iterations of ipalm (with backtracking) 17/31
18 Numerical Results - Test Case 1 Comparison of the RGB image, f and the recovered image u : RGB image f u 18/31
19 Numerical Results - Test Case 1 Comparison of the RGB image, f and u (zoom): RGB image f u 19/31
20 Numerical Results - Test Case 1 How well does the solution u fit the given data f? f S C k u f S C k u 20/31
21 Numerical Results - Test Case 2 low-res hyperspectral channel ( px) high-res RGB image ( px) 21/31
22 Numerical Results - Test Case 2 (u, k ) after 5000 iterations of ipalm (with backtracking) 22/31
23 Numerical Results - Test Case 2 Comparison of the RGB image, f and the recovered image u : RGB image f u 23/31
24 Numerical Results - Test Case 2 Comparison of the RGB image, f and u (zoom): RGB image f u 24/31
25 Numerical Results - Test Case 2 How well does the solution u fit the given data f? f S C k u f S C k u 25/31
26 Numerical Results - Blind Deconvolution vs. Fixed Kernel with blind deconvolution with fixed Gaussian kernel 26/31
27 Numerical Results - TV vs. dtv dtv TV 27/31
28 Conclusion and Outlook dtv regularization can be used to combine two imaging modalities in a super resolution task. Including blind deconvolution in the optimization process can be crucial in correctly linking the two modalities. Reconstructing one pixel channel in Matlab currently takes about 2 hours. This could be improved by a C++ or GPU-based implementation. Our solutions still suffer from high-frequency point artifacts. 28/31
29 Thank you! 29/31
30 Shameless Advertising 30/31
Super-resolution on Text Image Sequences
November 4, 2004 Outline Outline Geometric Distortion Optical/Motion Blurring Down-Sampling Total Variation Basic Idea Outline Geometric Distortion Optical/Motion Blurring Down-Sampling No optical/image
More informationGuided Image Super-Resolution: A New Technique for Photogeometric Super-Resolution in Hybrid 3-D Range Imaging
Guided Image Super-Resolution: A New Technique for Photogeometric Super-Resolution in Hybrid 3-D Range Imaging Florin C. Ghesu 1, Thomas Köhler 1,2, Sven Haase 1, Joachim Hornegger 1,2 04.09.2014 1 Pattern
More informationACCELERATED DUAL GRADIENT-BASED METHODS FOR TOTAL VARIATION IMAGE DENOISING/DEBLURRING PROBLEMS. Donghwan Kim and Jeffrey A.
ACCELERATED DUAL GRADIENT-BASED METHODS FOR TOTAL VARIATION IMAGE DENOISING/DEBLURRING PROBLEMS Donghwan Kim and Jeffrey A. Fessler University of Michigan Dept. of Electrical Engineering and Computer Science
More informationEfficient MR Image Reconstruction for Compressed MR Imaging
Efficient MR Image Reconstruction for Compressed MR Imaging Junzhou Huang, Shaoting Zhang, and Dimitris Metaxas Division of Computer and Information Sciences, Rutgers University, NJ, USA 08854 Abstract.
More informationThe Trainable Alternating Gradient Shrinkage method
The Trainable Alternating Gradient Shrinkage method Carlos Malanche, supervised by Ha Nguyen April 24, 2018 LIB (prof. Michael Unser), EPFL Quick example 1 The reconstruction problem Challenges in image
More informationImage Restoration. Diffusion Denoising Deconvolution Super-resolution Tomographic Reconstruction
Image Restoration Image Restoration Diffusion Denoising Deconvolution Super-resolution Tomographic Reconstruction Diffusion Term Consider only the regularization term E-L equation: (Laplace equation) Steepest
More informationConvex optimization algorithms for sparse and low-rank representations
Convex optimization algorithms for sparse and low-rank representations Lieven Vandenberghe, Hsiao-Han Chao (UCLA) ECC 2013 Tutorial Session Sparse and low-rank representation methods in control, estimation,
More informationApplication of Proximal Algorithms to Three Dimensional Deconvolution Microscopy
Application of Proximal Algorithms to Three Dimensional Deconvolution Microscopy Paroma Varma Stanford University paroma@stanford.edu Abstract In microscopy, shot noise dominates image formation, which
More informationConvex Optimization MLSS 2015
Convex Optimization MLSS 2015 Constantine Caramanis The University of Texas at Austin The Optimization Problem minimize : f (x) subject to : x X. The Optimization Problem minimize : f (x) subject to :
More informationOptimization for Machine Learning
Optimization for Machine Learning (Problems; Algorithms - C) SUVRIT SRA Massachusetts Institute of Technology PKU Summer School on Data Science (July 2017) Course materials http://suvrit.de/teaching.html
More informationSome Blind Deconvolution Techniques in Image Processing
Some Blind Deconvolution Techniques in Image Processing Tony Chan Math Dept., UCLA Joint work with Frederick Park and Andy M. Yip IPAM Workshop on Mathematical Challenges in Astronomical Imaging July 26-30,
More informationTotal Variation Denoising with Overlapping Group Sparsity
1 Total Variation Denoising with Overlapping Group Sparsity Ivan W. Selesnick and Po-Yu Chen Polytechnic Institute of New York University Brooklyn, New York selesi@poly.edu 2 Abstract This paper describes
More informationNumerical Methods on the Image Processing Problems
Numerical Methods on the Image Processing Problems Department of Mathematics and Statistics Mississippi State University December 13, 2006 Objective Develop efficient PDE (partial differential equations)
More informationProximal operator and methods
Proximal operator and methods Master 2 Data Science, Univ. Paris Saclay Robert M. Gower Optimization Sum of Terms A Datum Function Finite Sum Training Problem The Training Problem Convergence GD I Theorem
More informationDeep Learning and Its Applications
Convolutional Neural Network and Its Application in Image Recognition Oct 28, 2016 Outline 1 A Motivating Example 2 The Convolutional Neural Network (CNN) Model 3 Training the CNN Model 4 Issues and Recent
More informationTotal variation tomographic inversion via the Alternating Direction Method of Multipliers
Total variation tomographic inversion via the Alternating Direction Method of Multipliers Landon Safron and Mauricio D. Sacchi Department of Physics, University of Alberta Summary Geophysical inverse problems
More informationSparse Optimization Lecture: Proximal Operator/Algorithm and Lagrange Dual
Sparse Optimization Lecture: Proximal Operator/Algorithm and Lagrange Dual Instructor: Wotao Yin July 2013 online discussions on piazza.com Those who complete this lecture will know learn the proximal
More informationICCV 2013 Supplementary Material Super-resolution via Transform-invariant Group-sparse Regularization
ICCV 2013 Supplementary Material Super-resolution via Transform-invariant Group-sparse Regularization October 4, 2013 1 Contents This is the supplementary material for the paper Super-resolution via Transform-invariant
More informationThe Benefit of Tree Sparsity in Accelerated MRI
The Benefit of Tree Sparsity in Accelerated MRI Chen Chen and Junzhou Huang Department of Computer Science and Engineering, The University of Texas at Arlington, TX, USA 76019 Abstract. The wavelet coefficients
More informationCompressive Sensing Algorithms for Fast and Accurate Imaging
Compressive Sensing Algorithms for Fast and Accurate Imaging Wotao Yin Department of Computational and Applied Mathematics, Rice University SCIMM 10 ASU, Tempe, AZ Acknowledgements: results come in part
More informationA Patch Prior for Dense 3D Reconstruction in Man-Made Environments
A Patch Prior for Dense 3D Reconstruction in Man-Made Environments Christian Häne 1, Christopher Zach 2, Bernhard Zeisl 1, Marc Pollefeys 1 1 ETH Zürich 2 MSR Cambridge October 14, 2012 A Patch Prior for
More informationSpicyMKL Efficient multiple kernel learning method using dual augmented Lagrangian
SpicyMKL Efficient multiple kernel learning method using dual augmented Lagrangian Taiji Suzuki Ryota Tomioka The University of Tokyo Graduate School of Information Science and Technology Department of
More informationRobust Single Image Super-resolution based on Gradient Enhancement
Robust Single Image Super-resolution based on Gradient Enhancement Licheng Yu, Hongteng Xu, Yi Xu and Xiaokang Yang Department of Electronic Engineering, Shanghai Jiaotong University, Shanghai 200240,
More informationGRID WARPING IN TOTAL VARIATION IMAGE ENHANCEMENT METHODS. Andrey Nasonov, and Andrey Krylov
GRID WARPING IN TOTAL VARIATION IMAGE ENHANCEMENT METHODS Andrey Nasonov, and Andrey Krylov Lomonosov Moscow State University, Moscow, Department of Computational Mathematics and Cybernetics, e-mail: nasonov@cs.msu.ru,
More informationBlind Image Deblurring Using Dark Channel Prior
Blind Image Deblurring Using Dark Channel Prior Jinshan Pan 1,2,3, Deqing Sun 2,4, Hanspeter Pfister 2, and Ming-Hsuan Yang 3 1 Dalian University of Technology 2 Harvard University 3 UC Merced 4 NVIDIA
More informationMathematics in Image Processing
Mathematics in Image Processing Michal Šorel Department of Image Processing Institute of Information Theory and Automation (ÚTIA) Academy of Sciences of the Czech Republic http://zoi.utia.cas.cz/ Mathematics
More informationVariational Networks: Connecting Variational Methods and Deep Learning
Variational Networks: Connecting Variational Methods and Deep Learning Erich Kobler 1, Teresa Klatzer 1, Kerstin Hammernik 1 and Thomas Pock 1,2 1 Institute of Computer Graphics and Vision, Graz University
More informationLearning Data Terms for Non-blind Deblurring Supplemental Material
Learning Data Terms for Non-blind Deblurring Supplemental Material Jiangxin Dong 1, Jinshan Pan 2, Deqing Sun 3, Zhixun Su 1,4, and Ming-Hsuan Yang 5 1 Dalian University of Technology dongjxjx@gmail.com,
More informationConvex Optimization / Homework 2, due Oct 3
Convex Optimization 0-725/36-725 Homework 2, due Oct 3 Instructions: You must complete Problems 3 and either Problem 4 or Problem 5 (your choice between the two) When you submit the homework, upload a
More informationComputer Vision 2. SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung. Computer Vision 2 Dr. Benjamin Guthier
Computer Vision 2 SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung Computer Vision 2 Dr. Benjamin Guthier 1. IMAGE PROCESSING Computer Vision 2 Dr. Benjamin Guthier Content of this Chapter Non-linear
More informationA primal-dual framework for mixtures of regularizers
A primal-dual framework for mixtures of regularizers Baran Gözcü baran.goezcue@epfl.ch Laboratory for Information and Inference Systems (LIONS) École Polytechnique Fédérale de Lausanne (EPFL) Switzerland
More informationOptical Flow. Adriana Bocoi and Elena Pelican. 1. Introduction
Proceedings of the Fifth Workshop on Mathematical Modelling of Environmental and Life Sciences Problems Constanţa, Romania, September, 200, pp. 5 5 Optical Flow Adriana Bocoi and Elena Pelican This paper
More informationImage denoising using TV-Stokes equation with an orientation-matching minimization
Image denoising using TV-Stokes equation with an orientation-matching minimization Xue-Cheng Tai 1,2, Sofia Borok 1, and Jooyoung Hahn 1 1 Division of Mathematical Sciences, School of Physical Mathematical
More informationImage Processing and related PDEs Lecture 1: Introduction to image processing
Image Processing and related PDEs Lecture 1: Introduction to image processing Yves van Gennip School of Mathematical Sciences, University of Nottingham Minicourse on Image Processing and related PDEs University
More informationBandwidth Selection for Kernel Density Estimation Using Total Variation with Fourier Domain Constraints
IEEE SIGNAL PROCESSING LETTERS 1 Bandwidth Selection for Kernel Density Estimation Using Total Variation with Fourier Domain Constraints Alexander Suhre, Orhan Arikan, Member, IEEE, and A. Enis Cetin,
More informationIterative regularization in intensity-modulated radiation therapy optimization. Carlsson, F. and Forsgren, A. Med. Phys. 33 (1), January 2006.
Iterative regularization in intensity-modulated radiation therapy optimization Carlsson, F. and Forsgren, A. Med. Phys. 33 (1), January 2006. 2 / 15 Plan 1 2 3 4 3 / 15 to paper The purpose of the paper
More informationImage deconvolution using a characterization of sharp images in wavelet domain
Image deconvolution using a characterization of sharp images in wavelet domain Hui Ji a,, Jia Li a, Zuowei Shen a, Kang Wang a a Department of Mathematics, National University of Singapore, Singapore,
More informationA Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation
, pp.162-167 http://dx.doi.org/10.14257/astl.2016.138.33 A Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation Liqiang Hu, Chaofeng He Shijiazhuang Tiedao University,
More informationOptimal Sampling Geometries for TV-Norm Reconstruction of fmri Data
Optimal Sampling Geometries for TV-Norm Reconstruction of fmri Data Oliver M. Jeromin, Student Member, IEEE, Vince D. Calhoun, Senior Member, IEEE, and Marios S. Pattichis, Senior Member, IEEE Abstract
More informationImage restoration. Lecture 14. Milan Gavrilovic Centre for Image Analysis Uppsala University
Image restoration Lecture 14 Milan Gavrilovic milan@cb.uu.se Centre for Image Analysis Uppsala University Computer Assisted Image Analysis 2009-05-08 M. Gavrilovic (Uppsala University) L14 Image restoration
More informationImage Processing. Filtering. Slide 1
Image Processing Filtering Slide 1 Preliminary Image generation Original Noise Image restoration Result Slide 2 Preliminary Classic application: denoising However: Denoising is much more than a simple
More informationSupplemental Material for Efficient MR Image Reconstruction for Compressed MR Imaging
Supplemental Material for Efficient MR Image Reconstruction for Compressed MR Imaging Paper ID: 1195 No Institute Given 1 More Experiment Results 1.1 Visual Comparisons We apply all methods on four 2D
More informationCS 565 Computer Vision. Nazar Khan PUCIT Lectures 15 and 16: Optic Flow
CS 565 Computer Vision Nazar Khan PUCIT Lectures 15 and 16: Optic Flow Introduction Basic Problem given: image sequence f(x, y, z), where (x, y) specifies the location and z denotes time wanted: displacement
More informationBilevel Sparse Coding
Adobe Research 345 Park Ave, San Jose, CA Mar 15, 2013 Outline 1 2 The learning model The learning algorithm 3 4 Sparse Modeling Many types of sensory data, e.g., images and audio, are in high-dimensional
More informationImage Restoration under Significant Additive Noise
1 Image Restoration under Significant Additive Noise Wenyi Zhao 1 and Art Pope Sarnoff Corporation 01 Washington Road, Princeton, NJ 08540, USA email: { wzhao, apope }@ieee.org Tel: 408-53-178 Abstract
More informationAlgorithms for convex optimization
Algorithms for convex optimization Michal Kočvara Institute of Information Theory and Automation Academy of Sciences of the Czech Republic and Czech Technical University kocvara@utia.cas.cz http://www.utia.cas.cz/kocvara
More informationA Proximal Algorithm for Total Variation-Based Phase Retrieval
A Proximal Algorithm for Total Variation-Based Phase Retrieval Anqi Fu Department of Electrical Engineering Stanford University 350 Serra Mall, Stanford, CA 94305 Email: anqif@stanford.edu Abstract In
More informationBlind Image Deconvolution Technique for Image Restoration using Ant Colony Optimization
Blind Image Deconvolution Technique for Image Restoration using Ant Colony Optimization Amandeep Kaur CEM Kapurthala, Punjab Vinay Chopra DAVIET Jalandhar Punjab ABSTRACT Image Restoration is a field of
More informationADVANCED IMAGE PROCESSING METHODS FOR ULTRASONIC NDE RESEARCH C. H. Chen, University of Massachusetts Dartmouth, N.
ADVANCED IMAGE PROCESSING METHODS FOR ULTRASONIC NDE RESEARCH C. H. Chen, University of Massachusetts Dartmouth, N. Dartmouth, MA USA Abstract: The significant progress in ultrasonic NDE systems has now
More informationIterative Shrinkage/Thresholding g Algorithms: Some History and Recent Development
Iterative Shrinkage/Thresholding g Algorithms: Some History and Recent Development Mário A. T. Figueiredo Instituto de Telecomunicações and Instituto Superior Técnico, Technical University of Lisbon PORTUGAL
More informationSuperresolution software manual
Superresolution software manual S. Villena, M. Vega, D. Babacan, J. Mateos R. Molina and A. K. Katsaggelos Version 1.0 Contact information: Rafael Molina Departamento de Ciencias de la Computación e I.
More informationArtificial Intelligence for Robotics: A Brief Summary
Artificial Intelligence for Robotics: A Brief Summary This document provides a summary of the course, Artificial Intelligence for Robotics, and highlights main concepts. Lesson 1: Localization (using Histogram
More informationIterative CT Reconstruction Using Curvelet-Based Regularization
Iterative CT Reconstruction Using Curvelet-Based Regularization Haibo Wu 1,2, Andreas Maier 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science, 2 Graduate School in
More informationOptic Flow and Basics Towards Horn-Schunck 1
Optic Flow and Basics Towards Horn-Schunck 1 Lecture 7 See Section 4.1 and Beginning of 4.2 in Reinhard Klette: Concise Computer Vision Springer-Verlag, London, 2014 1 See last slide for copyright information.
More informationImage Restoration using Accelerated Proximal Gradient method
Image Restoration using Accelerated Proximal Gradient method Alluri.Samuyl Department of computer science, KIET Engineering College. D.Srinivas Asso.prof, Department of computer science, KIET Engineering
More informationDynamic Imaging from Incomplete Data Applications in MRI and PET
Dynamic Imaging from Incomplete Data Applications in MRI and PET Martin Benning, and many many others Department of Applied Mathematics & Theoretical Physics Challenges in Dynamic Imaging Data @ Isaac
More informationx 1 SpaceNet: Multivariate brain decoding and segmentation Elvis DOHMATOB (Joint work with: M. EICKENBERG, B. THIRION, & G. VAROQUAUX) 75 x y=20
SpaceNet: Multivariate brain decoding and segmentation Elvis DOHMATOB (Joint work with: M. EICKENBERG, B. THIRION, & G. VAROQUAUX) L R 75 x 2 38 0 y=20-38 -75 x 1 1 Introducing the model 2 1 Brain decoding
More informationInvestigating Image Inpainting via the Alternating Direction Method of Multipliers
Investigating Image Inpainting via the Alternating Direction Method of Multipliers Jonathan Tuck Stanford University jonathantuck@stanford.edu Abstract In many imaging applications, there exists potential
More informationAn edge-preserving multilevel method for deblurring, denoising, and segmentation
An edge-preserving multilevel method for deblurring, denoising, and segmentation S. Morigi 1, L. Reichel 2, and F. Sgallari 1 1 Dept. of Mathematics-CIRAM, University of Bologna, Bologna, Italy {morigi,sgallari}@dm.unibo.it
More informationRegularized Energy Minimization Models in Image Processing. Tibor Lukić. Faculty of Technical Sciences, University of Novi Sad, Serbia
20150715 NOVI SAD Regularized Energy Minimization Models in Image Processing Tibor Lukić Faculty of Technical Sciences, University of Novi Sad, Serbia SSIP Szeged, 2015 NOVI SAD NOVI SAD UNIVERSITY OF
More informationAliasing. Can t draw smooth lines on discrete raster device get staircased lines ( jaggies ):
(Anti)Aliasing and Image Manipulation for (y = 0; y < Size; y++) { for (x = 0; x < Size; x++) { Image[x][y] = 7 + 8 * sin((sqr(x Size) + SQR(y Size)) / 3.0); } } // Size = Size / ; Aliasing Can t draw
More informationComposite Self-concordant Minimization
Composite Self-concordant Minimization Volkan Cevher Laboratory for Information and Inference Systems-LIONS Ecole Polytechnique Federale de Lausanne (EPFL) volkan.cevher@epfl.ch Paris 6 Dec 11, 2013 joint
More informationSEMI-BLIND IMAGE RESTORATION USING A LOCAL NEURAL APPROACH
SEMI-BLIND IMAGE RESTORATION USING A LOCAL NEURAL APPROACH Ignazio Gallo, Elisabetta Binaghi and Mario Raspanti Universitá degli Studi dell Insubria Varese, Italy email: ignazio.gallo@uninsubria.it ABSTRACT
More informationCombining Total Variation and Nonlocal Means Regularization for Edge Preserving Image Deconvolution
Electronic Letters on Computer Vision and Image Analysis 10(1):42-53, 2011 Combining Total Variation and Nonlocal Means Regularization for Edge Preserving Image Deconvolution Binbin Hao and Jianguang Zhu
More informationFeature Tracking and Optical Flow
Feature Tracking and Optical Flow Prof. D. Stricker Doz. G. Bleser Many slides adapted from James Hays, Derek Hoeim, Lana Lazebnik, Silvio Saverse, who 1 in turn adapted slides from Steve Seitz, Rick Szeliski,
More informationReconstructing Reflective and Transparent Surfaces from Epipolar Plane Images
Sven Wanner and Bastian Goldlücke Heidelberg Collaboratory for Image Processing 1. Light Fields and Epipolar Plane Images 1. Light Fields and Epipolar Plane Images Light Fields as a dense sampling of a
More informationScaling Distributed Machine Learning
Scaling Distributed Machine Learning with System and Algorithm Co-design Mu Li Thesis Defense CSD, CMU Feb 2nd, 2017 nx min w f i (w) Distributed systems i=1 Large scale optimization methods Large-scale
More informationImage Deconvolution.
Image Deconvolution. Mathematics of Imaging. HW3 Jihwan Kim Abstract This homework is to implement image deconvolution methods, especially focused on a ExpectationMaximization(EM) algorithm. Most of this
More informationRegularization and Markov Random Fields (MRF) CS 664 Spring 2008
Regularization and Markov Random Fields (MRF) CS 664 Spring 2008 Regularization in Low Level Vision Low level vision problems concerned with estimating some quantity at each pixel Visual motion (u(x,y),v(x,y))
More informationFeature Tracking and Optical Flow
Feature Tracking and Optical Flow Prof. D. Stricker Doz. G. Bleser Many slides adapted from James Hays, Derek Hoeim, Lana Lazebnik, Silvio Saverse, who in turn adapted slides from Steve Seitz, Rick Szeliski,
More informationSuper Resolution Using Graph-cut
Super Resolution Using Graph-cut Uma Mudenagudi, Ram Singla, Prem Kalra, and Subhashis Banerjee Department of Computer Science and Engineering Indian Institute of Technology Delhi Hauz Khas, New Delhi,
More informationEdges, interpolation, templates. Nuno Vasconcelos ECE Department, UCSD (with thanks to David Forsyth)
Edges, interpolation, templates Nuno Vasconcelos ECE Department, UCSD (with thanks to David Forsyth) Edge detection edge detection has many applications in image processing an edge detector implements
More informationLecture 19: November 5
0-725/36-725: Convex Optimization Fall 205 Lecturer: Ryan Tibshirani Lecture 9: November 5 Scribes: Hyun Ah Song Note: LaTeX template courtesy of UC Berkeley EECS dept. Disclaimer: These notes have not
More informationCOMPUTER VISION > OPTICAL FLOW UTRECHT UNIVERSITY RONALD POPPE
COMPUTER VISION 2017-2018 > OPTICAL FLOW UTRECHT UNIVERSITY RONALD POPPE OUTLINE Optical flow Lucas-Kanade Horn-Schunck Applications of optical flow Optical flow tracking Histograms of oriented flow Assignment
More informationRobust Principal Component Analysis (RPCA)
Robust Principal Component Analysis (RPCA) & Matrix decomposition: into low-rank and sparse components Zhenfang Hu 2010.4.1 reference [1] Chandrasekharan, V., Sanghavi, S., Parillo, P., Wilsky, A.: Ranksparsity
More informationSYSTEMS OF NONLINEAR EQUATIONS
SYSTEMS OF NONLINEAR EQUATIONS Widely used in the mathematical modeling of real world phenomena. We introduce some numerical methods for their solution. For better intuition, we examine systems of two
More informationBlur Space Iterative De-blurring
Blur Space Iterative De-blurring RADU CIPRIAN BILCU 1, MEJDI TRIMECHE 2, SAKARI ALENIUS 3, MARKKU VEHVILAINEN 4 1,2,3,4 Multimedia Technologies Laboratory, Nokia Research Center Visiokatu 1, FIN-33720,
More informationFast Monotonic Blind Deconvolution Algorithm for Constrained TV Based Image Restoration
Fast Monotonic Blind Deconvolution Algorithm for Constrained TV Based Image Restoration Haiying Liu, W.-S. Lu, Yanan Fu, Yu Cheng, Tingfang Yan, Teng Li and Max Q.-H. Meng School of Control Science and
More informationGaussian Processes for Robotics. McGill COMP 765 Oct 24 th, 2017
Gaussian Processes for Robotics McGill COMP 765 Oct 24 th, 2017 A robot must learn Modeling the environment is sometimes an end goal: Space exploration Disaster recovery Environmental monitoring Other
More informationOn the regularizing power of multigrid-type algorithms
On the regularizing power of multigrid-type algorithms Marco Donatelli Dipartimento di Fisica e Matematica Università dell Insubria www.mat.unimi.it/user/donatel Outline 1 Restoration of blurred and noisy
More informationIntroduction to Image Super-resolution. Presenter: Kevin Su
Introduction to Image Super-resolution Presenter: Kevin Su References 1. S.C. Park, M.K. Park, and M.G. KANG, Super-Resolution Image Reconstruction: A Technical Overview, IEEE Signal Processing Magazine,
More information3D Computer Vision. Dense 3D Reconstruction II. Prof. Didier Stricker. Christiano Gava
3D Computer Vision Dense 3D Reconstruction II Prof. Didier Stricker Christiano Gava Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche Intelligenz http://av.dfki.de
More informationOpen Access Noisy Radial-blurred Images Restoration Using the BM3D Frames
Send Orders for Reprints to reprints@benthamscience.ae 54 The Open Cybernetics & Systemics Journal, 2014, 8, 54-60 Open Access Noisy Radial-blurred Images Restoration Using the BM3D Frames Li Chonglun*,
More informationOptimization Plugin for RapidMiner. Venkatesh Umaashankar Sangkyun Lee. Technical Report 04/2012. technische universität dortmund
Optimization Plugin for RapidMiner Technical Report Venkatesh Umaashankar Sangkyun Lee 04/2012 technische universität dortmund Part of the work on this technical report has been supported by Deutsche Forschungsgemeinschaft
More informationImage Processing. Traitement d images. Yuliya Tarabalka Tel.
Traitement d images Yuliya Tarabalka yuliya.tarabalka@hyperinet.eu yuliya.tarabalka@gipsa-lab.grenoble-inp.fr Tel. 04 76 82 62 68 Noise reduction Image restoration Restoration attempts to reconstruct an
More informationExpectation Maximization and Total Variation Based Model for Computed Tomography Reconstruction from Undersampled Data
Expectation Maximization and Total Variation Based Model for Computed Tomography Reconstruction from Undersampled Data Ming Yan and Luminita A. Vese Department of Mathematics, University of California,
More informationA Multilevel Proximal Gradient Algorithm for a Class of Composite Optimization Problems
A Multilevel Proximal Gradient Algorithm for a Class of Composite Optimization Problems Panos Parpas May 9, 2017 Abstract Composite optimization models consist of the minimization of the sum of a smooth
More informationREGULARIZATION PARAMETER TRIMMING FOR ITERATIVE IMAGE RECONSTRUCTION
REGULARIZATION PARAMETER TRIMMING FOR ITERATIVE IMAGE RECONSTRUCTION Haoyi Liang and Daniel S. Weller University of Virginia, Department of ECE, Charlottesville, VA, 2294, USA ABSTRACT Conventional automatic
More informationFinding a Needle in a Haystack: An Image Processing Approach
Finding a Needle in a Haystack: An Image Processing Approach Emily Beylerian Advisor: Hayden Schaeffer University of California, Los Angeles Department of Mathematics Abstract Image segmentation (also
More informationAn Adaptive Directional Haar Framelet-based Reconstruction Algorithm for Parallel Magnetic Resonance Imaging
An Adaptive Directional Haar Framelet-based Reconstruction Algorithm for Parallel Magnetic Resonance Imaging Yan-Ran Li Raymond H. Chan Lixin Shen Yung-Chin Hsu Wen-Yih Isaac Tseng February, 206 Abstract
More informationLeaves Machine Learning and Optimization Library
Leaves Machine Learning and Optimization Library October 7, 07 This document describes how to use LEaves Machine Learning and Optimization Library (LEMO) for modeling. LEMO is an open source library for
More informationReducing Effects of Bad Data Using Variance Based Joint Sparsity Recovery
Journal of Scientific Computing manuscript No. (will be inserted by the editor) Reducing Effects of Bad Data Using Variance Based Joint Sparsity Recovery Anne Gelb Theresa Scarnati Received: date / Accepted:
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
More informationA Derivative-Free Approximate Gradient Sampling Algorithm for Finite Minimax Problems
1 / 33 A Derivative-Free Approximate Gradient Sampling Algorithm for Finite Minimax Problems Speaker: Julie Nutini Joint work with Warren Hare University of British Columbia (Okanagan) III Latin American
More information4D Computed Tomography Reconstruction from Few-Projection Data via Temporal Non-local Regularization
4D Computed Tomography Reconstruction from Few-Projection Data via Temporal Non-local Regularization Xun Jia 1, Yifei Lou 2, Bin Dong 3, Zhen Tian 1,4, and Steve Jiang 1 1 Department of Radiation Oncology
More informationDense Image-based Motion Estimation Algorithms & Optical Flow
Dense mage-based Motion Estimation Algorithms & Optical Flow Video A video is a sequence of frames captured at different times The video data is a function of v time (t) v space (x,y) ntroduction to motion
More informationADAPTIVE REGULARIZATION OF THE NL-MEANS: APPLICATION TO IMAGE AND VIDEO DENOISING 1
ADAPTIVE REGULARIZATION OF THE NL-MEANS: APPLICATION TO IMAGE AND VIDEO DENOISING 1 Adaptive regularization of the NL-means: Application to image and video denoising Camille Sutour, Jean-François Aujol,
More informationConvex Optimization. Lijun Zhang Modification of
Convex Optimization Lijun Zhang zlj@nju.edu.cn http://cs.nju.edu.cn/zlj Modification of http://stanford.edu/~boyd/cvxbook/bv_cvxslides.pdf Outline Introduction Convex Sets & Functions Convex Optimization
More informationThe p-laplacian on Graphs with Applications in Image Processing and Classification
The p-laplacian on Graphs with Applications in Image Processing and Classification Abderrahim Elmoataz 1,2, Matthieu Toutain 1, Daniel Tenbrinck 3 1 UMR6072 GREYC, Université de Caen Normandie 2 Université
More informationDigital Image Processing ERRATA. Wilhelm Burger Mark J. Burge. An algorithmic introduction using Java. Second Edition. Springer
Wilhelm Burger Mark J. Burge Digital Image Processing An algorithmic introduction using Java Second Edition ERRATA Springer Berlin Heidelberg NewYork Hong Kong London Milano Paris Tokyo 12.1 RGB Color
More information