MULTI-REGION SEGMENTATION
|
|
- Beverly Kelley
- 5 years ago
- Views:
Transcription
1 MULTI-REGION SEGMENTATION USING GRAPH-CUTS Johannes Ulén Abstract This project deals with multi-region segmenation using graph-cuts and is mainly based on a paper by Delong and Boykov [1]. The difference between a multi-label model and multi-region model is that certain geometrical constraints can be added with the multi-region models.
2 1 Introduction A segmentation is a partitioning of an image into multiple segments. This can be the partitioning of an image into one segment belonging to a object in the image and one segment belonging to the background in the image. Each pixel is given a label corresponding to which segment it is part of. In multi-label segmentation more than one object is searched for in the image and in multi-region segmentation each pixel can be part of more than one label. In applications information about an image can be known a priori. When segmenting an image of the heart you know before you start that there are two chambers (in a healthy patient) and that the chambers are surrounded by the heart muscle. This information can be used to improve the segmentation by imposing restriction in between the different regions. More on this later. We first need some background on graphs. A graph G (V, E) is a collection of data-points called vertices (V ) and connections between the vertices called edges (E). Each edge has a capacity c, where we let e(i, j) = c denote an added edge between vertice i and j with capacity c. In our treatment e (i, j) and e (j, i) may differ, e (i, j) 0 i, j and e(i, i) = 0, i.e we have no loops in the graph. s min-cut Figure 1: The source (s) and sink (t) is indicated by squares, the vertices by circles and edges by directed arrows. The value over and edge indicate the capacity of the edge. The Minimum s t cut is indicated by the thick line. Removing the three edges covered by the line will stop all flow from s to t. Graphs can be used to formulate discrete optimization problems such as the max-flow problem. In the problem flow of any kind, say oil, should be transported from a vertice called source (s) to a vertice called sink (t). The goal is too find a route which maximizes the flow, this formulation gives rise to a so called s t graph. The minimum s t cut, or min-cut, on the graph is the set of edges with minimal capacity such that there are no edges connecting the source with the sink, an example is given in Figure 1. The term graph-cuts comes from the search for the minimal cut in a s t graph. A well known result in graph theory is the max-flow/mincut theorem [2], which states the maximum flow through a s t graph is the same as the minimum cut in the graph. There are very effective algorithms for finding the maximum 2 5 t flow, and thus also the minimum cut. The one we will use [3] is specifically developed to solve the kind of graphs which arises in computer vision problems. It works by pushing flow from both the source and sink, instead of just pushing flow from the source like standard algorithms. 2 Segmentation 2.1 Binary segmentation In binary segmentation the image should be partitioned into two different segments, one belonging to the object and one belonging to the background. We introduce P as the set of pixel indices and a binary variable for each pixel such that x p B P for every p P. When x p = 0 pixel p is a background pixel and when x p = 1 the pixel is part of the object. The intensity of any pixel p is given by I (p). The segmentation problem can formulated in the form of minimizing the energy function: E (x) = data terms }} }} D p (x p ) + V p,q (x p, x q ), (1) p P p,q N i j regularization terms where there are many different chooses of data and regularization terms. Data terms adds a cost dependent on how different I (p) is compared to the expected value of either the background E b or the object E o. One possible data term is: (I (p) E b ) 2, if x p = 0 D p (x p ) = (I (p) E o ) 2 (2), if x p = 1. Regularization terms adds a cost if neighboring pixels are similar but have different labels. N can be any pixel connectivity and is usually chosen as the -neighborhood. A possible regularization term is 0, if xp = x q V p,q (x p, x q ) = ( λ 1 (I (p) I (q)) 2) if x p x q, (3) where λ can be used to regulate the influence from the regularization term on the energy function Graph construction For a min-cut on any graph we set x p = 0 for all vertices connected to the source and x p = 1 for all vertices con- 1
3 nected to the sink. The energy function E (x) can be reformulated into a min-cut problem where the minimal cut leads to the same labeling as the global minimum of E (x). In [] a general method for converting the energy minimization problem into a min-cut problem is presented, which works as long as E (x) is a submodular function of up to three variables. Definition 1. A function f of two variables with the domain 0, 1} is called submodular if f (0, 0) + f (1, 1) f (1, 0) + f (0, 1). () A function of n variables is called submodular if from fixing n 2 variables to any values it follows that the two free variables fulfill (). Definition 2. A supermodular function is the same a submodular function with the inequality reversed. It is possible to minimize E (x) via graphs even when the function is not submodular [5], however global optimum is not guaranteed. These methods lay outside the scope of this project. We now shortly describe how the graph is constructed for the binary segmentation problem. We begin with the data term, for each pixel p, add the edge e (s, p) = D p (1) D p (0), if D p (1) > D p (0) e(p, t) = D p (0) D p (1), otherwise. For the regularization term we begin by splitting V p,q into the four different combinations of labeling for neighboring pixels. V p,q = V p,q (0, 0) V p,q (0, 1) V p,q (1, 0) V p,q (1, 1) = A B C D We continue by rewriting the expression as A B C D = A C A C A + 0 D C 0 D C + 0 B + C A D 0 0 The first term is a positive constant so we don t need any edges for it as it will increase E (x) with the same amount no matter what labeling is chosen. The second term can be implemented by adding one edge e (s, p) = C A, if C > A e (p, t) = A C, otherwise. The third term can be handled just as easy by adding another edge Figure 2: Examples on how E (x) is turned into a graph. LEFT: Edge for D p (x p ) when D p (1) > D p (0) RIGHT: Edges for V pq (x p, x q ) when D C > 0 and C A < 0. e (s, q) = D C, e (q, t) = C D, if D > C otherwise. The last term could be very tricky if it wasn t for the submodularity, we know that it must be positive and it suffices to add the edge e(p, q) = B + C A D. An example on how the edges are added is given in figure 2. To better understand why this work it is recommended to work through the example given in the figure and validate that any cut indeed give the same cost as V p,q (x p, x q ) User seeds If we a priori know the label of some pixels we can use this information to calculate an estimation of E b and E o. The information can also be added to the graph by adding the edges: e (s, p) = e (p, t) = p known to be background, p known to be object. By adding these edges we prohibit cuts which will contradict our a priori knowledge. 2.2 Multi-region segmentation We now extend the binary segmentation to multiple regions. The notation from [1] is adopted. Define L to be the set of region indices. We define a new binary variables x B P L which we index as x i p with i L and p P. If x i p = 1 pixel p is interior in region i. Let x p denote the vector of all region variables that correspond to pixel p. If x p = 0, then pixel p belongs to the background. The basic idea in [1] is to let each pixel be represented by L vertices in the graph. For L = 1 this will lead to the same graph as in binary segmentation. For L = 2 a new 2
4 i contains j x i p x j q Wpq ij i excludes j x i p x j q Wpq ij i attracts j x i p x j q Wpq ij α Figure 3: Graph construction for a pixel image. Each layer i and j have 16 vertices corresponding to each pixel. The geometric interaction terms are constructed by adding edges in the inter-region neighborhood N ij which can be seen as an extension from the normal neighborhood in an image. layer of vertices is added along with edges which interacts between the different layers, figure 3. Define the multi-region energy function as data }} }} E (x) = D p (x p ) + V i ( x i) }} + W ij ( x i, x j), p P regularization i L geometric interaction i,j L i j below we explain each part of the energy function. Data terms defines the cost for every combination of regions, i.e each pixel can either be inside or outside each region i L leading up to 2 L different costs. For L = 2 the data terms can be reformulated into a graph exactly like the regularization term for binary segmentation. Regularization terms is independent for each region i and edges are only being added between vertices belonging to the same layer i. For each region i the regularization terms are equal to those of binary segmentation. Geometric interaction terms encodes all geometric interactions between regions i and j. The inter-region neighborhood N ij is the set of all pixel pairs (p, q) at which region i is assigned some geometric interaction with region j, see figure 3 for an example. Formally we define the terms as: W ij = pq N ij W ij pq ( x i p, x i ) q. In table 1, W ij pq is given for three different geometric interactions. These terms can be reformulated into the graph just like the regularization term for binary segmentation. Table 1: Energy terms for different geometric interactions. For the i contains j interaction an edge with infinite capacity is added between x i p to x j p prohibiting any cut which separates them. Equal statements can be made for the exclusion and attraction terms. The exclusion term is everywhere supermodular, we can for simple interactions just flip the meaning of x i p and turn the problem into a everywhere submodular problem Three regions We will now explain how to solve a three-region segmentation problem for a special segmentation problem which is useful in medical image analysis. Assume we have three regions 1 3 and that we know a priori that region 2 and 3 are interior of region 1 and that region 2 and 3 should not overlap. We enforce the interior geometric constraint by adding edges with infinite capacity going from the vertices used to represent region 1 to those used to represent region 2 and 3. The data term will now be a function of three variables. We introduce the function ( D = D p x 1 p, x 2 p, x 3 ) p as the cost for pixel p to be interior in any combination for the three regions. We can represent this function as in table 2. Table 2: D i,j,k (0, 0, 0) D i,j,k (0, 0, 1) D i,j,k (0, 1, 0) D i,j,k (0, 1, 1) D i,j,k (1, 0, 0) D i,j,k (1, 0, 1) D i,j,k (1, 1, 0) D i,j,k (1, 1, 1) We introduce = A C E G B D F H P = (A + D + F + G) (B + C + E + H), from [] we know that D can be turned into a graph-cut problem if for every pixel p the following holds. Case 1: For every pixel p where P 0 then the following must also hold: 3
5 A + D A + F A + G B + C B + E. (5) C + E Figure Case 2: For every pixel p where P < 0 then the following must also hold: E + H C + H B + H F + G D + G. (6) D + F The energies B, C, D can be set to any finite value as the geometric interaction terms will make sure that these combinations never happened. Unfortunately we cannot have both containment and exclusion interactions for the same pixel without turning the energy function supermodular [1], the trick with flipping the meaning of x i p does not work. This means that the cost of H does matter. Since we don t want region 2 and 3 to overlap our goal is to find a maximum H such that D is still submodular. We begin with case 1. To maximize H we directly see that A + D + F + G = B + C + E + H. (7) From (7) and (5) it is evident that we must maximize D in order to maximize H which gives Inserting (8) into (7) yields: A + D = B + C. (8) F + G = E + H H = F + G E. The value of B and C does not matter for the value of H, it must however fulfill (5) so we set: B = A + F + G E, C = A + F + G E, which hold as long as G, F 0. In case 2 we see that we cannot achieve a better value for H because of (6). For the details on how to construct the graph from this the reader is referred to [], and note that P = 0 for the maximum chose of H. 3 Experiments In all experiments the data term is chosen as (2) and the regularization terms like (3) with λ = 1 if nothing else is stated. (a) Motivational example image. (b) LEFT: Region 1 in blue overlay. RIGHT: Region 2 in red overlay. 3.1 Two regions Motivation We begin with motivational example for a two-region segmentation. We have the image in Figure 5a, and assume we are only interested in segmenting out gray regions when they are contained inside a black region. This requirement can be enforced with geometric interaction terms leading up to the resulting segmentation found in Figure 5b. Here no user seeds are used, only expected intensities of both regions MR-segmentation We continue with the more difficult problem of segmenting out both the myocardium (heart muscle) and the left and right ventricles (heart chambers) of a short-axis images (a direction in which the MR-images are captured). We know a priori that both chambers will be surrounded by myocardium. We set region 1 to be the myocardium and region 2 to be both the left and right ventricle, we then force region 2 to be contained inside region 1. We give the algorithm some user seeds, figure 6a. The resulting segmentation can be seen in Figure 6b, where λ = 1/5 was used as parameter for the regularization term and connected components with less than 50 pixels was filtered away in a post processing step.
6 Figure 6 Figure 5 (a) A simple image containing three regions. User seeds is shown in colors. (a) LEFT: A short axis image of the human heart. The lighter area around the letters LV and RV are the left and right ventricles. The grayish area around both ventricles is the myocardium. The myocardium is much thicker around left ventricle compared to the right ventricle. RIGHT: Seeds used to segment the heart. (b) For both images the user seeds are shown in a darker overlay. LEFT: Region 1, myocardium and the left and right ventricles shown in blue overlay. RIGHT: Region 2, the left and right ventricles shown in red overlay. (b) LEFT: Region 1 in blue overlay. MIDDLE: Region 2 in red overlay. RIGHT: Region 3 in green overlay. 3.2 Three regions Basic result We begin the experiments with three region segmentation with an example showing that it is possible to get excellent results with a single cut, to do this we use the synthetic image found in Figure 7a along with the user seeds given in the same image. In Figure 7b the resulting segmentation is shown Submodular limitations The restriction on the energy functions to be submodular can lead to troubles even for some very simple image such as the one found in Figure 7a. As was derived in Section the maximum value of H = F + G E. A normal choice of data terms for the regions we are interested in may be: A = D p (0, 0, 0) = κ + (I (p) E (Background)) 2 E = D p (1, 0, 0) = κ + (I (p) E (Region 1)) 2 F = D p (1, 0, 1) = κ + (I (p) E (Region 2)) 2 G = D p (1, 1, 0) = κ + (I (p) E (Region 3)) 2 H = D p (1, 1, 1) = F + G E, where E ( ) is the expected value of a region and κ any positive constant. And in here lies the problem if F + G > 5
7 2E and H < A; H will be the smallest energy for that pixel. In Figure 7(b)-(f), the energy associate with labeling each pixel to A and E H is given. As can be seen the lowest energy is attained if all white pixels are labeled to H. Increasing κ decreases the relative difference but destroys the contrast in energies MR segmentation We redo the MR segmentation in Section but we make the left and right ventricle into separate regions. Region 1 is the set to be the myocardium, region 2 the left ventricle and region 3 the right ventricle. We use geometric interaction terms to force region 2 and 3 to be inside of region 1. The user seeds are shown in Figure 8a. The resulting segmentation when λ = 1/5 is used as parameter for the regularization term can be seen in Figure 8b. Region 2 and Region 3 overlap due to the problem discussed in Section (a) Example image where submodularity causes problem. User seeds are given in yellow, blue, red and green. The scale is used in all images in this Figure. Discussion It would be desirable to model the short axis image as three regions where the left and right ventricle are to different regions; but the experiment in Section gave poor result. This is because we are limited to submodular energies. The results can be improved if supermodular energies are used which can be solved by Quadratic Pseudo-Boolean Optimization [5]. (b) A = D p (0, 0, 0) (c) E = D p (1, 0, 0) References [1] A. Delong and Y. Boykov, Globally Optimal Segmentation of Multi-Region Objects, ICCV, [2] L. Ford and D. Fulkerson, Maximal Flow Through a Network, Canadian Journal of Mathematics, (d) F = D p (1, 0, 1) (e) G = D p (1, 1, 0) [3] Y. Boykov and V. Kolmogorov, An experimental comparison of min-cut/max- flow algorithms for energy minimization in vision, Pattern Analysis and Machine Intelligence, IEEE Transactions on, 200. [] V. Kolmogorov and R. Zabih, What Energy Functions Can Be Minimized via Graph Cuts?, Pattern Analysis and Machine Intelligence, IEEE Transactions on, 200. [5] V. Kolmogorov and C. Rother, Minimizing Nonsubmodular Functions with Graph Cuts-a Review, IEEE transactions on pattern analysis and machine intelligence, (f) H = D p (1, 1, 1) Figure 7: The problem with a submodular constraint on H, H = F + G E. 6
8 (a) A short axis images of the human heart with user seeds given by the user in yellow, blue, red and green. (b) LEFT: Region 1 shown in blue overlay. MIDDLE: Region 2 shown in red overlay. RIGHT: Region 3 shown in green overlay. Figure 8: Segmentation of a short-axis image with three a region model. 7
Supplementary Material: Adaptive and Move Making Auxiliary Cuts for Binary Pairwise Energies
Supplementary Material: Adaptive and Move Making Auxiliary Cuts for Binary Pairwise Energies Lena Gorelick Yuri Boykov Olga Veksler Computer Science Department University of Western Ontario Abstract Below
More informationCS6670: Computer Vision
CS6670: Computer Vision Noah Snavely Lecture 19: Graph Cuts source S sink T Readings Szeliski, Chapter 11.2 11.5 Stereo results with window search problems in areas of uniform texture Window-based matching
More informationMEDICAL IMAGE COMPUTING (CAP 5937) LECTURE 10: Medical Image Segmentation as an Energy Minimization Problem
SPRING 06 MEDICAL IMAGE COMPUTING (CAP 97) LECTURE 0: Medical Image Segmentation as an Energy Minimization Problem Dr. Ulas Bagci HEC, Center for Research in Computer Vision (CRCV), University of Central
More informationMEDICAL IMAGE COMPUTING (CAP 5937) LECTURE 10: Medical Image Segmentation as an Energy Minimization Problem
SPRING 07 MEDICAL IMAGE COMPUTING (CAP 97) LECTURE 0: Medical Image Segmentation as an Energy Minimization Problem Dr. Ulas Bagci HEC, Center for Research in Computer Vision (CRCV), University of Central
More informationMaximum flows and minimal cuts. Filip Malmberg
Maximum flows and minimal cuts Filip Malmberg MSF cuts I A Minimum Spanning Forest with respect to a set of seed-points divides a graph into a number of connected components. Thus, it defines a cut on
More informationCS 5540 Spring 2013 Assignment 3, v1.0 Due: Apr. 24th 11:59PM
1 Introduction In this programming project, we are going to do a simple image segmentation task. Given a grayscale image with a bright object against a dark background and we are going to do a binary decision
More informationintro, applications MRF, labeling... how it can be computed at all? Applications in segmentation: GraphCut, GrabCut, demos
Image as Markov Random Field and Applications 1 Tomáš Svoboda, svoboda@cmp.felk.cvut.cz Czech Technical University in Prague, Center for Machine Perception http://cmp.felk.cvut.cz Talk Outline Last update:
More informationCAP5415-Computer Vision Lecture 13-Image/Video Segmentation Part II. Dr. Ulas Bagci
CAP-Computer Vision Lecture -Image/Video Segmentation Part II Dr. Ulas Bagci bagci@ucf.edu Labeling & Segmentation Labeling is a common way for modeling various computer vision problems (e.g. optical flow,
More informationSupervised texture detection in images
Supervised texture detection in images Branislav Mičušík and Allan Hanbury Pattern Recognition and Image Processing Group, Institute of Computer Aided Automation, Vienna University of Technology Favoritenstraße
More informationADAPTIVE GRAPH CUTS WITH TISSUE PRIORS FOR BRAIN MRI SEGMENTATION
ADAPTIVE GRAPH CUTS WITH TISSUE PRIORS FOR BRAIN MRI SEGMENTATION Abstract: MIP Project Report Spring 2013 Gaurav Mittal 201232644 This is a detailed report about the course project, which was to implement
More informationDiscrete Optimization. Lecture Notes 2
Discrete Optimization. Lecture Notes 2 Disjunctive Constraints Defining variables and formulating linear constraints can be straightforward or more sophisticated, depending on the problem structure. The
More informationCS4670 / 5670: Computer Vision Noah Snavely
{ { 11/26/2013 CS4670 / 5670: Computer Vision Noah Snavely Graph-Based Image Segmentation Stereo as a minimization problem match cost Want each pixel to find a good match in the other image smoothness
More informationGlobally Optimal Segmentation of Multi-Region Objects
In International Conference on Computer Vision (ICCV), Kyoto, October 2009 Globally Optimal Segmentation of Multi-Region Objects Andrew Delong University of Western Ontario Yuri oykov University of Western
More informationConnected components - 1
Connected Components Basic definitions Connectivity, Adjacency, Connected Components Background/Foreground, Boundaries Run-length encoding Component Labeling Recursive algorithm Two-scan algorithm Chain
More informationCS443: Digital Imaging and Multimedia Binary Image Analysis. Spring 2008 Ahmed Elgammal Dept. of Computer Science Rutgers University
CS443: Digital Imaging and Multimedia Binary Image Analysis Spring 2008 Ahmed Elgammal Dept. of Computer Science Rutgers University Outlines A Simple Machine Vision System Image segmentation by thresholding
More informationBrain tumor segmentation with minimal user assistance
Western University Scholarship@Western Electronic Thesis and Dissertation Repository January 2014 Brain tumor segmentation with minimal user assistance Liqun Liu The University of Western Ontario Supervisor
More informationMarkov/Conditional Random Fields, Graph Cut, and applications in Computer Vision
Markov/Conditional Random Fields, Graph Cut, and applications in Computer Vision Fuxin Li Slides and materials from Le Song, Tucker Hermans, Pawan Kumar, Carsten Rother, Peter Orchard, and others Recap:
More informationPRIMAL-DUAL INTERIOR POINT METHOD FOR LINEAR PROGRAMMING. 1. Introduction
PRIMAL-DUAL INTERIOR POINT METHOD FOR LINEAR PROGRAMMING KELLER VANDEBOGERT AND CHARLES LANNING 1. Introduction Interior point methods are, put simply, a technique of optimization where, given a problem
More informationGraph Based Image Segmentation
AUTOMATYKA 2011 Tom 15 Zeszyt 3 Anna Fabijañska* Graph Based Image Segmentation 1. Introduction Image segmentation is one of the fundamental problems in machine vision. In general it aims at extracting
More informationMinimum Cost Edge Disjoint Paths
Minimum Cost Edge Disjoint Paths Theodor Mader 15.4.2008 1 Introduction Finding paths in networks and graphs constitutes an area of theoretical computer science which has been highly researched during
More informationDigital Image Processing. Prof. P.K. Biswas. Department of Electronics & Electrical Communication Engineering
Digital Image Processing Prof. P.K. Biswas Department of Electronics & Electrical Communication Engineering Indian Institute of Technology, Kharagpur Image Segmentation - III Lecture - 31 Hello, welcome
More informationInteractive Differential Segmentation of the Prostate using Graph-Cuts with a Feature Detector-based Boundary Term
MOSCHIDIS, GRAHAM: GRAPH-CUTS WITH FEATURE DETECTORS 1 Interactive Differential Segmentation of the Prostate using Graph-Cuts with a Feature Detector-based Boundary Term Emmanouil Moschidis emmanouil.moschidis@postgrad.manchester.ac.uk
More informationFast Approximate Energy Minimization via Graph Cuts
IEEE Transactions on PAMI, vol. 23, no. 11, pp. 1222-1239 p.1 Fast Approximate Energy Minimization via Graph Cuts Yuri Boykov, Olga Veksler and Ramin Zabih Abstract Many tasks in computer vision involve
More informationSegmentation with non-linear constraints on appearance, complexity, and geometry
IPAM February 2013 Western Univesity Segmentation with non-linear constraints on appearance, complexity, and geometry Yuri Boykov Andrew Delong Lena Gorelick Hossam Isack Anton Osokin Frank Schmidt Olga
More informationModels for grids. Computer vision: models, learning and inference. Multi label Denoising. Binary Denoising. Denoising Goal.
Models for grids Computer vision: models, learning and inference Chapter 9 Graphical Models Consider models where one unknown world state at each pixel in the image takes the form of a grid. Loops in the
More informationSTIC AmSud Project. Graph cut based segmentation of cardiac ventricles in MRI: a shape-prior based approach
STIC AmSud Project Graph cut based segmentation of cardiac ventricles in MRI: a shape-prior based approach Caroline Petitjean A joint work with Damien Grosgeorge, Pr Su Ruan, Pr JN Dacher, MD October 22,
More informationDr. Ulas Bagci
CAP- Computer Vision Lecture - Image Segmenta;on as an Op;miza;on Problem Dr. Ulas Bagci bagci@ucf.edu Reminders Oct Guest Lecture: SVM by Dr. Gong Oct 8 Guest Lecture: Camera Models by Dr. Shah PA# October
More informationLecture 3: Art Gallery Problems and Polygon Triangulation
EECS 396/496: Computational Geometry Fall 2017 Lecture 3: Art Gallery Problems and Polygon Triangulation Lecturer: Huck Bennett In this lecture, we study the problem of guarding an art gallery (specified
More informationLecture 20 : Trees DRAFT
CS/Math 240: Introduction to Discrete Mathematics 4/12/2011 Lecture 20 : Trees Instructor: Dieter van Melkebeek Scribe: Dalibor Zelený DRAFT Last time we discussed graphs. Today we continue this discussion,
More informationCSE 417 Network Flows (pt 3) Modeling with Min Cuts
CSE 417 Network Flows (pt 3) Modeling with Min Cuts Reminders > HW6 is due on Friday start early bug fixed on line 33 of OptimalLineup.java: > change true to false Review of last two lectures > Defined
More information3 No-Wait Job Shops with Variable Processing Times
3 No-Wait Job Shops with Variable Processing Times In this chapter we assume that, on top of the classical no-wait job shop setting, we are given a set of processing times for each operation. We may select
More informationRandomized algorithms have several advantages over deterministic ones. We discuss them here:
CS787: Advanced Algorithms Lecture 6: Randomized Algorithms In this lecture we introduce randomized algorithms. We will begin by motivating the use of randomized algorithms through a few examples. Then
More informationContrast adjustment via Bayesian sequential partitioning
Contrast adjustment via Bayesian sequential partitioning Zhiyu Wang, Shuo Xie, Bai Jiang Abstract Photographs taken in dim light have low color contrast. However, traditional methods for adjusting contrast
More informationInteractive Brain Tumor Segmentation with Inclusion Constraints
Western University Scholarship@Western Electronic Thesis and Dissertation Repository April 2016 Interactive Brain Tumor Segmentation with Inclusion Constraints Danfeng Chen The University of Western Ontario
More informationFigure : Example Precedence Graph
CS787: Advanced Algorithms Topic: Scheduling with Precedence Constraints Presenter(s): James Jolly, Pratima Kolan 17.5.1 Motivation 17.5.1.1 Objective Consider the problem of scheduling a collection of
More informationBasic relations between pixels (Chapter 2)
Basic relations between pixels (Chapter 2) Lecture 3 Basic Relationships Between Pixels Definitions: f(x,y): digital image Pixels: q, p (p,q f) A subset of pixels of f(x,y): S A typology of relations:
More informationModern Medical Image Analysis 8DC00 Exam
Parts of answers are inside square brackets [... ]. These parts are optional. Answers can be written in Dutch or in English, as you prefer. You can use drawings and diagrams to support your textual answers.
More informationLecture 11: Maximum flow and minimum cut
Optimisation Part IB - Easter 2018 Lecture 11: Maximum flow and minimum cut Lecturer: Quentin Berthet 4.4. The maximum flow problem. We consider in this lecture a particular kind of flow problem, with
More informationFigurative Tours and Braids
Figurative Tours and Braids Robert Bosch Dept. of Mathematics Oberlin College rbosch@oberlin.edu Tom Wexler Dept. of Computer Science Oberlin College wexler@cs.oberlin.edu Abstract We start with a rectangular
More informationBiomedical Image Analysis. Mathematical Morphology
Biomedical Image Analysis Mathematical Morphology Contents: Foundation of Mathematical Morphology Structuring Elements Applications BMIA 15 V. Roth & P. Cattin 265 Foundations of Mathematical Morphology
More informationLOGISTIC REGRESSION FOR MULTIPLE CLASSES
Peter Orbanz Applied Data Mining Not examinable. 111 LOGISTIC REGRESSION FOR MULTIPLE CLASSES Bernoulli and multinomial distributions The mulitnomial distribution of N draws from K categories with parameter
More informationMethods and Models for Combinatorial Optimization Exact methods for the Traveling Salesman Problem
Methods and Models for Combinatorial Optimization Exact methods for the Traveling Salesman Problem L. De Giovanni M. Di Summa The Traveling Salesman Problem (TSP) is an optimization problem on a directed
More informationCHAPTER 4 VORONOI DIAGRAM BASED CLUSTERING ALGORITHMS
CHAPTER 4 VORONOI DIAGRAM BASED CLUSTERING ALGORITHMS 4.1 Introduction Although MST-based clustering methods are effective for complex data, they require quadratic computational time which is high for
More informationGraphs, graph algorithms (for image segmentation),... in progress
Graphs, graph algorithms (for image segmentation),... in progress Václav Hlaváč Czech Technical University in Prague Czech Institute of Informatics, Robotics and Cybernetics 66 36 Prague 6, Jugoslávských
More informationSuperpixels via Pseudo-Boolean Optimization
Superpixels via Pseudo-Boolean Optimization Yuhang Zhang, Richard Hartley The Australian National University {yuhang.zhang, richard.hartley}@anu.edu.au John Mashford, Stewart Burn CSIRO {john.mashford,
More informationGEOMETRIC SEGMENTATION USING GAUSSIAN WAVE PACKETS
Proceedings of the Project Review, Geo-Mathematical Imaging Group (Purdue University, West Lafayette IN), Vol. 1 (2011) pp. 57-64. GEOMETRIC SEGMENTATION USING GAUSSIAN WAVE PACKETS FREDRIK ANDERSSON AND
More informationIntroduction à la vision artificielle X
Introduction à la vision artificielle X Jean Ponce Email: ponce@di.ens.fr Web: http://www.di.ens.fr/~ponce Planches après les cours sur : http://www.di.ens.fr/~ponce/introvis/lect10.pptx http://www.di.ens.fr/~ponce/introvis/lect10.pdf
More informationLecture 16 Segmentation and Scene understanding
Lecture 16 Segmentation and Scene understanding Introduction! Mean-shift! Graph-based segmentation! Top-down segmentation! Silvio Savarese Lecture 15 -! 3-Mar-14 Segmentation Silvio Savarese Lecture 15
More informations-t Graph Cuts for Binary Energy Minimization
s-t Graph Cuts for Binary Energy Minimization E data term ( L) = D ( ) + ( ) P Lp λ Ι Lp Lq p prior term Now that we have an energy function, the big question is how do we minimize it? pq N n Exhaustive
More informationComputer Vision I - Basics of Image Processing Part 2
Computer Vision I - Basics of Image Processing Part 2 Carsten Rother 07/11/2014 Computer Vision I: Basics of Image Processing Roadmap: Basics of Digital Image Processing Computer Vision I: Basics of Image
More informationUsing Game Theory for Image Segmentation
Using Game Theory for Image Segmentation Elizabeth Cassell Sumanth Kolar Alex Yakushev 1 Introduction 21st March 2007 The goal of image segmentation, is to distinguish objects from background. Robust segmentation
More informationApproximate decoding: ICM, block methods, alpha-beta swap & alpha-expansion
Approximate decoding: ICM, block methods, alpha-beta swap & alpha-expansion Julieta Martinez University of British Columbia August 25, 2015 https://www.youtube.com/watch?v=kp3ik5f3-2c&t=18m36s Outline
More informationCOMP 558 lecture 22 Dec. 1, 2010
Binocular correspondence problem Last class we discussed how to remap the pixels of two images so that corresponding points are in the same row. This is done by computing the fundamental matrix, defining
More informationInteger Programming Theory
Integer Programming Theory Laura Galli October 24, 2016 In the following we assume all functions are linear, hence we often drop the term linear. In discrete optimization, we seek to find a solution x
More informationQuadratic Pseudo-Boolean Optimization(QPBO): Theory and Applications At-a-Glance
Quadratic Pseudo-Boolean Optimization(QPBO): Theory and Applications At-a-Glance Presented By: Ahmad Al-Kabbany Under the Supervision of: Prof.Eric Dubois 12 June 2012 Outline Introduction The limitations
More informationDiscrete Optimization Methods in Computer Vision CSE 6389 Slides by: Boykov Modified and Presented by: Mostafa Parchami Basic overview of graph cuts
Discrete Optimization Methods in Computer Vision CSE 6389 Slides by: Boykov Modified and Presented by: Mostafa Parchami Basic overview of graph cuts [Yuri Boykov, Olga Veksler, Ramin Zabih, Fast Approximation
More informationMaximal Monochromatic Geodesics in an Antipodal Coloring of Hypercube
Maximal Monochromatic Geodesics in an Antipodal Coloring of Hypercube Kavish Gandhi April 4, 2015 Abstract A geodesic in the hypercube is the shortest possible path between two vertices. Leader and Long
More informationby conservation of flow, hence the cancelation. Similarly, we have
Chapter 13: Network Flows and Applications Network: directed graph with source S and target T. Non-negative edge weights represent capacities. Assume no edges into S or out of T. (If necessary, we can
More informationImage Segmentation Using Iterated Graph Cuts Based on Multi-scale Smoothing
Image Segmentation Using Iterated Graph Cuts Based on Multi-scale Smoothing Tomoyuki Nagahashi 1, Hironobu Fujiyoshi 1, and Takeo Kanade 2 1 Dept. of Computer Science, Chubu University. Matsumoto 1200,
More informationDavid Van Brackle Chief Judge, Southeast USA Region, ICPC
David Van Brackle Chief Judge, Southeast USA Region, ICPC Division 1 (2 correct/30 submissions) Given a list of scheduled flights, airports, turn times, etc., determine the minimum number of planes needed.
More informationData Term. Michael Bleyer LVA Stereo Vision
Data Term Michael Bleyer LVA Stereo Vision What happened last time? We have looked at our energy function: E ( D) = m( p, dp) + p I < p, q > N s( p, q) We have learned about an optimization algorithm that
More informationExploring Domains of Approximation in R 2 : Expository Essay
Exploring Domains of Approximation in R 2 : Expository Essay Nicolay Postarnakevich August 12, 2013 1 Introduction In this paper I explore the concept of the Domains of Best Approximations. These structures
More informationAn Eternal Domination Problem in Grids
Theory and Applications of Graphs Volume Issue 1 Article 2 2017 An Eternal Domination Problem in Grids William Klostermeyer University of North Florida, klostermeyer@hotmail.com Margaret-Ellen Messinger
More informationDistributed minimum spanning tree problem
Distributed minimum spanning tree problem Juho-Kustaa Kangas 24th November 2012 Abstract Given a connected weighted undirected graph, the minimum spanning tree problem asks for a spanning subtree with
More informationComputer Vision : Exercise 4 Labelling Problems
Computer Vision Exercise 4 Labelling Problems 13/01/2014 Computer Vision : Exercise 4 Labelling Problems Outline 1. Energy Minimization (example segmentation) 2. Iterated Conditional Modes 3. Dynamic Programming
More informationComputer Vision I - Filtering and Feature detection
Computer Vision I - Filtering and Feature detection Carsten Rother 30/10/2015 Computer Vision I: Basics of Image Processing Roadmap: Basics of Digital Image Processing Computer Vision I: Basics of Image
More informationCharacterizing and Controlling the. Spectral Output of an HDR Display
Characterizing and Controlling the Spectral Output of an HDR Display Ana Radonjić, Christopher G. Broussard, and David H. Brainard Department of Psychology, University of Pennsylvania, Philadelphia, PA
More informationSupport Vector Machines.
Support Vector Machines srihari@buffalo.edu SVM Discussion Overview 1. Overview of SVMs 2. Margin Geometry 3. SVM Optimization 4. Overlapping Distributions 5. Relationship to Logistic Regression 6. Dealing
More informationA Novel Image Classification Model Based on Contourlet Transform and Dynamic Fuzzy Graph Cuts
Appl. Math. Inf. Sci. 6 No. 1S pp. 93S-97S (2012) Applied Mathematics & Information Sciences An International Journal @ 2012 NSP Natural Sciences Publishing Cor. A Novel Image Classification Model Based
More informationUsing the Kolmogorov-Smirnov Test for Image Segmentation
Using the Kolmogorov-Smirnov Test for Image Segmentation Yong Jae Lee CS395T Computational Statistics Final Project Report May 6th, 2009 I. INTRODUCTION Image segmentation is a fundamental task in computer
More information2 Review of Set Theory
2 Review of Set Theory Example 2.1. Let Ω = {1, 2, 3, 4, 5, 6} 2.2. Venn diagram is very useful in set theory. It is often used to portray relationships between sets. Many identities can be read out simply
More informationStructural and Syntactic Pattern Recognition
Structural and Syntactic Pattern Recognition Selim Aksoy Department of Computer Engineering Bilkent University saksoy@cs.bilkent.edu.tr CS 551, Fall 2017 CS 551, Fall 2017 c 2017, Selim Aksoy (Bilkent
More informationDiversity Coloring for Distributed Storage in Mobile Networks
Diversity Coloring for Distributed Storage in Mobile Networks Anxiao (Andrew) Jiang and Jehoshua Bruck California Institute of Technology Abstract: Storing multiple copies of files is crucial for ensuring
More informationOpen and Closed Sets
Open and Closed Sets Definition: A subset S of a metric space (X, d) is open if it contains an open ball about each of its points i.e., if x S : ɛ > 0 : B(x, ɛ) S. (1) Theorem: (O1) and X are open sets.
More informationLecture 3: Some Strange Properties of Fractal Curves
Lecture 3: Some Strange Properties of Fractal Curves I have been a stranger in a strange land. Exodus 2:22 1. Fractal Strangeness Fractals have a look and feel that is very different from ordinary curves.
More informationDiscrete mathematics , Fall Instructor: prof. János Pach
Discrete mathematics 2016-2017, Fall Instructor: prof. János Pach - covered material - Lecture 1. Counting problems To read: [Lov]: 1.2. Sets, 1.3. Number of subsets, 1.5. Sequences, 1.6. Permutations,
More informationFormations in flow fields
Formations in flow fields Rick van Meer June 11, 2015 1 Introduction Pathfinding and formations (i.e. an arrangement of agents) are both a part of artificial intelligence and can be used in games. However,
More informationComputation of the Constrained Infinite Time Linear Quadratic Optimal Control Problem
Computation of the Constrained Infinite Time Linear Quadratic Optimal Control Problem July 5, Introduction Abstract Problem Statement and Properties In this paper we will consider discrete-time linear
More informationCS521 \ Notes for the Final Exam
CS521 \ Notes for final exam 1 Ariel Stolerman Asymptotic Notations: CS521 \ Notes for the Final Exam Notation Definition Limit Big-O ( ) Small-o ( ) Big- ( ) Small- ( ) Big- ( ) Notes: ( ) ( ) ( ) ( )
More informationGraph Cuts in Vision and Graphics: Theories and Applications
In Handbook of Mathematical Models in Computer Vision, Springer, 2006 p.1 Graph Cuts in Vision and Graphics: Theories and Applications Yuri Boykov and Olga Veksler Computer Science, The University of Western
More informationEE 584 MACHINE VISION
EE 584 MACHINE VISION Binary Images Analysis Geometrical & Topological Properties Connectedness Binary Algorithms Morphology Binary Images Binary (two-valued; black/white) images gives better efficiency
More information6.2 Classification of Closed Surfaces
Table 6.1: A polygon diagram 6.1.2 Second Proof: Compactifying Teichmuller Space 6.2 Classification of Closed Surfaces We saw that each surface has a triangulation. Compact surfaces have finite triangulations.
More informationLecture 1: Introduction
Lecture 1 1 Linear and Combinatorial Optimization Anders Heyden Centre for Mathematical Sciences Lecture 1: Introduction The course and its goals Basic concepts Optimization Combinatorial optimization
More informationGRMA: Generalized Range Move Algorithms for the Efficient Optimization of MRFs
International Journal of Computer Vision manuscript No. (will be inserted by the editor) GRMA: Generalized Range Move Algorithms for the Efficient Optimization of MRFs Kangwei Liu Junge Zhang Peipei Yang
More informationSet Cover with Almost Consecutive Ones Property
Set Cover with Almost Consecutive Ones Property 2004; Mecke, Wagner Entry author: Michael Dom INDEX TERMS: Covering Set problem, data reduction rules, enumerative algorithm. SYNONYMS: Hitting Set PROBLEM
More information/ Approximation Algorithms Lecturer: Michael Dinitz Topic: Linear Programming Date: 2/24/15 Scribe: Runze Tang
600.469 / 600.669 Approximation Algorithms Lecturer: Michael Dinitz Topic: Linear Programming Date: 2/24/15 Scribe: Runze Tang 9.1 Linear Programming Suppose we are trying to approximate a minimization
More informationOBJECT detection in general has many applications
1 Implementing Rectangle Detection using Windowed Hough Transform Akhil Singh, Music Engineering, University of Miami Abstract This paper implements Jung and Schramm s method to use Hough Transform for
More informationSimultaneous Segmentation and Filtering via Reduced Graph Cuts
Simultaneous Segmentation and Filtering via Reduced Graph Cuts N. Lermé F. Malgouyres LAGA UMR CNRS 7539 IMT UMR CNRS 5219 LIPN UMR CNRS 7030 Université Paul Sabatier Université Paris 13 fmalgouy@math.univ-toulouse.fr
More informationCS599: Convex and Combinatorial Optimization Fall 2013 Lecture 1: Introduction to Optimization. Instructor: Shaddin Dughmi
CS599: Convex and Combinatorial Optimization Fall 013 Lecture 1: Introduction to Optimization Instructor: Shaddin Dughmi Outline 1 Course Overview Administrivia 3 Linear Programming Outline 1 Course Overview
More informationLecture 9. Semidefinite programming is linear programming where variables are entries in a positive semidefinite matrix.
CSE525: Randomized Algorithms and Probabilistic Analysis Lecture 9 Lecturer: Anna Karlin Scribe: Sonya Alexandrova and Keith Jia 1 Introduction to semidefinite programming Semidefinite programming is linear
More informationInteractive Treatment Planning in Cancer Radiotherapy
Interactive Treatment Planning in Cancer Radiotherapy Mohammad Shakourifar Giulio Trigila Pooyan Shirvani Ghomi Abraham Abebe Sarah Couzens Laura Noreña Wenling Shang June 29, 212 1 Introduction Intensity
More informationOn ɛ-unit distance graphs
On ɛ-unit distance graphs Geoffrey Exoo Department of Mathematics and Computer Science Indiana State University Terre Haute, IN 47809 g-exoo@indstate.edu April 9, 003 Abstract We consider a variation on
More informationNP-Hardness. We start by defining types of problem, and then move on to defining the polynomial-time reductions.
CS 787: Advanced Algorithms NP-Hardness Instructor: Dieter van Melkebeek We review the concept of polynomial-time reductions, define various classes of problems including NP-complete, and show that 3-SAT
More informationDelaunay Triangulations. Presented by Glenn Eguchi Computational Geometry October 11, 2001
Delaunay Triangulations Presented by Glenn Eguchi 6.838 Computational Geometry October 11, 2001 Motivation: Terrains Set of data points A R 2 Height ƒ(p) defined at each point p in A How can we most naturally
More informationPacific Journal of Mathematics
Pacific Journal of Mathematics SIMPLIFYING TRIANGULATIONS OF S 3 Aleksandar Mijatović Volume 208 No. 2 February 2003 PACIFIC JOURNAL OF MATHEMATICS Vol. 208, No. 2, 2003 SIMPLIFYING TRIANGULATIONS OF S
More informationThe Encoding Complexity of Network Coding
The Encoding Complexity of Network Coding Michael Langberg Alexander Sprintson Jehoshua Bruck California Institute of Technology Email: mikel,spalex,bruck @caltech.edu Abstract In the multicast network
More informationDerivation of an Algorithm for Calculation of the Intersection Area of a Circle with a Grid with Finite Fill Factor
Derivation of an Algorithm for Calculation of the Intersection Area of a Circle with a Grid with Finite Fill Factor Abstract: The problem deals with an exact calculation of the intersection area of a circle
More informationChapter 1 - Introduction
Chapter 1-lntroduction Chapter 1 - Introduction The aim of this chapter is to provide a background to topics which are relevant to the subject of this thesis. The motivation for writing a thesis regarding
More informationConstruction C : an inter-level coded version of Construction C
Construction C : an inter-level coded version of Construction C arxiv:1709.06640v2 [cs.it] 27 Dec 2017 Abstract Besides all the attention given to lattice constructions, it is common to find some very
More informationCombinatorial optimization and its applications in image Processing. Filip Malmberg
Combinatorial optimization and its applications in image Processing Filip Malmberg Part 1: Optimization in image processing Optimization in image processing Many image processing problems can be formulated
More information