Biomedical Image Processing

Size: px
Start display at page:

Download "Biomedical Image Processing"

Transcription

1 Biomedical Image Processing Jason Thong Gabriel Grant 1

2 2

3 Motivation from the Medical Perspective MRI, CT and other biomedical imaging devices were designed to assist doctors in their diagnosis and treatment of patients As engineers, we must design tools with the end-user in mind (medical professionals, not engineers) Need to know criteria of a good medical image this depends on the application Can optimize image processing to take advantage of the specific technology (MRI is better for soft tissues and has higher contrast resolution whereas CT is better for bones and has higher spatial resolution) Recent increase in computational power facilitates exploration of new image processing algorithms 3

4 Segmentation Divide the image by identifying different physical objects Three broad types: Automatic Semi-automatic Atlas-based (will be discussed later) 4

5 Segmentation > Automatic > Edge Detection Gray Level First Derivative 5

6 Segmentation > Automatic > Thresholding By inspecting image histogram, a threshold point can be selected 6

7 Segmentation > Automatic > Thresholding At higher- and lower-than-optimal points image quality suffers. Original, 20%, 40%, 55%, 75% Thresholds 7

8 Segmentation > Automatic > Thresholding Threshold point is not always so well defined Original, Histogram, 25%, 50%, 75% Thresholds 8

9 Segmentation > Semi-automatic > Region Growing Start from a seed a small group of pixels Grow the seed by adding neighboring pixels, based on edge rules Problems: Selection of a good seed Creation of effective edge rules Stop Condition(s) 9

10 Segmentation > Semi-automatic > Active contours (snakes) Allows user to roughly outline an object Computer corrects the outline 10

11 Noise Generally, noise in electronic sensors can come from electrons deviating from their theoretical path, caused mainly by heat Noise in biomedical images can also come from photon scattering (Rayleigh, Compton) The raw image is blurry This noise can be filtered Manipulating the pixels directly is called filtering in the time or spatial domain A more powerful method is to filter the image in the frequency domain, specifically, to modify the Fourier Transform of an image 11

12 Fourier Series and Transform The Fourier Series is a way to convert any signal into a summation of sinusoids But a n cos (nt) + b n sin (nt) can be expressed as c n cos (nt + θ). The Fourier Transform is kind of like c n as a function of n. The Inverse Fourier Transform 12

13 Interpreting the Fourier Transform In this example, there is a large peak at ω = 4 in the Fourier Transform This means the signal is close to but not exactly a sinusoid with a frequency of 4 rad/s 13

14 DFT (Discrete Fourier Transform) A computer works only on discrete data Input x n and transform X k are discrete and finite in length No negative frequency, thus summation starts at n=0 Cannot analyze frequency higher than data sampling rate, thus upper limit of summation is finite Inverse DFT Multidimensional DFT (images have 2 dimensions) 14

15 FFT (Fast Fourier Transform) A simple implementation is O(n 2 ) Many different FFT algorithms, most common is Cooley-Turkey Divide and conquer method the idea is to calculate the Fourier Transform from only the even index inputs, then calculate it for the odd index inputs, then add the two results Done recursively, thus requires log 2 n stages Total run time is O(n log 2 n) 15

16 FFT (continued) = WTF? M = N/2, so you have: (even) + (complex coefficient)*(odd) (even) (complex coefficient)*(odd) k in first half k in second half Important: 1 multiplication generates both terms 16

17 FFT In Block Form x n inputs X f FFT output Even index transform Odd index transform Just the final stage in FFT 17

18 Basic unit of computation The Butterfly Involves only 1 complex multiplication and 2 complex additions Complex coefficients can sometimes be factorized, thus requiring the minimum amount of computation 18

19 FFT Implementation The order of the inputs have changed Entering the final stage, the even indexes are on the top 4 lines and the odd on the bottom. Repeating this idea back towards the initial stage, the initial index is the reverse in binary of the final index (ex. s 4 -> S 1 ) 19

20 Filtering in the Frequency Domain Input image Pre-processing Fourier Transform Filter function Inverse Fourier Transform Post-processing Filtered image Time domain Frequency domain Time domain 20

21 Pre-Processing and Post-Processing Done in the time or spatial domain acts directly on the pixels Typically an image is multiplied by (-1) x+y so that the Fourier Transform of the new image is centered in the frequency domain (low frequencies are now in the center, high frequencies outside) The post-processing must inverse the pre-processing FT without (-1) x+y FT with (-1) x+y 21

22 Frequency Filtering Sharp transitions (like the edges between different tissues) are the high frequencies in a Fourier Transform these are the details in an image The general appearance (like the general shape of an object) are the low frequencies these come from the smooth areas of an image A high pass filter (in the frequency domain) will enhance the boundaries between different tissues in a biomedical image, but it will also enhance noise A low pass filter will smooth out an image, thus making it blurry Filtering in the frequency domain is done by multiplying the FT by a filter function. The equivalent filtering in the time domain would require convolution. 22

23 2D Fourier Transforms (A MatLab Experiment) Fourier Transform Real data Inverse Fourier Transform Complex data (MatLab function only shows real component of data) 23

24 IFT IFT 24

25 IFT IFT 25

26 Ideal Filters in Frequency Domain F(ω) is the Fourier Transform For a low pass filter: H(ω) 1 H(ω) = 0 F(ω) < F 0 1 F(ω) > F 0 F 0 is the cutoff frequency A high pass filter is simply 1 minus the low pass filter: H(ω) = 1 F(ω) < F0 0 F(ω) > F0 H(ω) 1 F 0 F(ω) F 0 F(ω) 26

27 Ideal Filters in Frequency Domain An ideal filter in the frequency domain is non-ideal in the time domain. The filtered image has artifacts. It turns out that an ideal filter in the time domain requires a non-ideal filter in the frequency domain. 27

28 Gaussian Filters H(ω) The Gaussian distribution is: 1 The Gaussian low pass filter is H(ω) = exp(-f 2 (ω) / 2F 02 ) F 0 = cutoff frequency H(ω) F 0 F(ω) The high pass filter is 1 minus the low pass filter 1 H(ω) = 1 - exp(-f 2 (ω) / 2F 02 ) F 0 The Gaussian filter does not create artifacts. F(ω) 28

29 A Filtered Image 29

30 Registration Matches two images of the same object or scene 30

31 Registration > Uses To fuse different sensors or times Matching PET to MR fmri To detect changes over time patient returns for second round of imaging To form one large image from several smaller ones To recognize an object by registering it's image to a model (atlasbased segmentation, mentioned above) 31

32 Registration > Method Transform sensed (target) image to match reference (source) image: preprocess determine transformation function (by applying to sensed image, deformed result has geometry of the reference image) resample the sensed image in the new geometry 32

33 Registration > Transformations Rigid body Translation Rotation Affine (adds scaling) Feature-based Elastic 33

34 Registration > Transformations > Rigid Body Principal Axes Registration (PAR) for full volume images (sensed and reference image must have same field-of-view) must compute: centroid principal axes 34

35 Registration > Transformations > Rigid Body Iterative Principal Axes Registration (IPAR) Can be used with partial volumes Useful for PET where field of view is less than MR 35

36 Registration > Transformations Iterative Feature-Based Registration matches "landmark points" from source and target Elastic Deformation Based Registration considers the target image to be deformable attempts to minimize deformation 36

37 Research at McMaster Jie Wu - Segmentation and Recognition Analysis of Human Organs in Magnetic Resonance Image Weiguang Guan, Research Engineer - Visualization Support with RHPCS C. Lambacher; W.F.S. Poehlman; M. Kamath; M. Noseworthy - Using a Plugin Architecture for Medical Image Processing 37

38 References Lathi, B. P. Linear Systems and Signals. New York: Oxford University Press, Pham, Dzung L. et al. A Survey of Current Methods in Medical Image Segmentation. Annual Review of Biomedical Imaging. January Zill, Dennis G. and Cullen, Michael R. Advanced Engineering Mathematics, Second Edition. Sudbury: Jones and Baerlett,

MEDICAL IMAGE ANALYSIS

MEDICAL IMAGE ANALYSIS SECOND EDITION MEDICAL IMAGE ANALYSIS ATAM P. DHAWAN g, A B IEEE Engineering in Medicine and Biology Society, Sponsor IEEE Press Series in Biomedical Engineering Metin Akay, Series Editor +IEEE IEEE PRESS

More information

Practical Image and Video Processing Using MATLAB

Practical Image and Video Processing Using MATLAB Practical Image and Video Processing Using MATLAB Chapter 18 Feature extraction and representation What will we learn? What is feature extraction and why is it a critical step in most computer vision and

More information

Image Segmentation and Registration

Image Segmentation and Registration Image Segmentation and Registration Dr. Christine Tanner (tanner@vision.ee.ethz.ch) Computer Vision Laboratory, ETH Zürich Dr. Verena Kaynig, Machine Learning Laboratory, ETH Zürich Outline Segmentation

More information

EE795: Computer Vision and Intelligent Systems

EE795: Computer Vision and Intelligent Systems EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 WRI C225 Lecture 04 130131 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Histogram Equalization Image Filtering Linear

More information

Computer 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 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 information

Vivekananda. Collegee of Engineering & Technology. Question and Answers on 10CS762 /10IS762 UNIT- 5 : IMAGE ENHANCEMENT.

Vivekananda. Collegee of Engineering & Technology. Question and Answers on 10CS762 /10IS762 UNIT- 5 : IMAGE ENHANCEMENT. Vivekananda Collegee of Engineering & Technology Question and Answers on 10CS762 /10IS762 UNIT- 5 : IMAGE ENHANCEMENT Dept. Prepared by Harivinod N Assistant Professor, of Computer Science and Engineering,

More information

Segmentation of Images

Segmentation of Images Segmentation of Images SEGMENTATION If an image has been preprocessed appropriately to remove noise and artifacts, segmentation is often the key step in interpreting the image. Image segmentation is a

More information

Biomedical Image Analysis. Spatial Filtering

Biomedical Image Analysis. Spatial Filtering Biomedical Image Analysis Contents: Spatial Filtering The mechanics of Spatial Filtering Smoothing and sharpening filters BMIA 15 V. Roth & P. Cattin 1 The Mechanics of Spatial Filtering Spatial filter:

More information

Image Registration. Prof. Dr. Lucas Ferrari de Oliveira UFPR Informatics Department

Image Registration. Prof. Dr. Lucas Ferrari de Oliveira UFPR Informatics Department Image Registration Prof. Dr. Lucas Ferrari de Oliveira UFPR Informatics Department Introduction Visualize objects inside the human body Advances in CS methods to diagnosis, treatment planning and medical

More information

Interpolation of 3D magnetic resonance data

Interpolation of 3D magnetic resonance data Interpolation of 3D magnetic resonance data J. Mikulka 1, E. Gescheidtova 2 and K. Bartusek 3 1, 2 Department of Theoretical and Experimental Electrical Engineering, Brno University of Technology, Kolejni

More information

Computer Vision. Fourier Transform. 20 January Copyright by NHL Hogeschool and Van de Loosdrecht Machine Vision BV All rights reserved

Computer Vision. Fourier Transform. 20 January Copyright by NHL Hogeschool and Van de Loosdrecht Machine Vision BV All rights reserved Van de Loosdrecht Machine Vision Computer Vision Fourier Transform 20 January 2017 Copyright 2001 2017 by NHL Hogeschool and Van de Loosdrecht Machine Vision BV All rights reserved j.van.de.loosdrecht@nhl.nl,

More information

ECG782: Multidimensional Digital Signal Processing

ECG782: Multidimensional Digital Signal Processing Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spatial Domain Filtering http://www.ee.unlv.edu/~b1morris/ecg782/ 2 Outline Background Intensity

More information

Digital Image Processing. Lecture 6

Digital Image Processing. Lecture 6 Digital Image Processing Lecture 6 (Enhancement in the Frequency domain) Bu-Ali Sina University Computer Engineering Dep. Fall 2016 Image Enhancement In The Frequency Domain Outline Jean Baptiste Joseph

More information

Digital Image Processing. Prof. P. K. Biswas. Department of Electronic & Electrical Communication Engineering

Digital Image Processing. Prof. P. K. Biswas. Department of Electronic & Electrical Communication Engineering Digital Image Processing Prof. P. K. Biswas Department of Electronic & Electrical Communication Engineering Indian Institute of Technology, Kharagpur Lecture - 21 Image Enhancement Frequency Domain Processing

More information

ENT 315 Medical Signal Processing CHAPTER 3 FAST FOURIER TRANSFORM. Dr. Lim Chee Chin

ENT 315 Medical Signal Processing CHAPTER 3 FAST FOURIER TRANSFORM. Dr. Lim Chee Chin ENT 315 Medical Signal Processing CHAPTER 3 FAST FOURIER TRANSFORM Dr. Lim Chee Chin Outline Definition and Introduction FFT Properties of FFT Algorithm of FFT Decimate in Time (DIT) FFT Steps for radix

More information

ECG782: Multidimensional Digital Signal Processing

ECG782: Multidimensional Digital Signal Processing Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 03 Image Processing Basics 13/01/28 http://www.ee.unlv.edu/~b1morris/ecg782/

More information

Outlines. Medical Image Processing Using Transforms. 4. Transform in image space

Outlines. Medical Image Processing Using Transforms. 4. Transform in image space Medical Image Processing Using Transforms Hongmei Zhu, Ph.D Department of Mathematics & Statistics York University hmzhu@yorku.ca Outlines Image Quality Gray value transforms Histogram processing Transforms

More information

Combinatorial optimization and its applications in image Processing. Filip Malmberg

Combinatorial 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

BME I5000: Biomedical Imaging

BME I5000: Biomedical Imaging 1 Lucas Parra, CCNY BME I5000: Biomedical Imaging Lecture 11 Point Spread Function, Inverse Filtering, Wiener Filtering, Sharpening,... Lucas C. Parra, parra@ccny.cuny.edu Blackboard: http://cityonline.ccny.cuny.edu/

More information

Computational Aspects of MRI

Computational Aspects of MRI David Atkinson Philip Batchelor David Larkman Programme 09:30 11:00 Fourier, sampling, gridding, interpolation. Matrices and Linear Algebra 11:30 13:00 MRI Lunch (not provided) 14:00 15:30 SVD, eigenvalues.

More information

Lecture 6: Edge Detection

Lecture 6: Edge Detection #1 Lecture 6: Edge Detection Saad J Bedros sbedros@umn.edu Review From Last Lecture Options for Image Representation Introduced the concept of different representation or transformation Fourier Transform

More information

PROSTATE CANCER DETECTION USING LABEL IMAGE CONSTRAINED MULTIATLAS SELECTION

PROSTATE CANCER DETECTION USING LABEL IMAGE CONSTRAINED MULTIATLAS SELECTION PROSTATE CANCER DETECTION USING LABEL IMAGE CONSTRAINED MULTIATLAS SELECTION Ms. Vaibhavi Nandkumar Jagtap 1, Mr. Santosh D. Kale 2 1 PG Scholar, 2 Assistant Professor, Department of Electronics and Telecommunication,

More information

Segmenting 2D Ultrasound Images using Seeded Region Growing

Segmenting 2D Ultrasound Images using Seeded Region Growing University of British Columbia Department of Electrical and Computer Engineering EECE 544 Medical Imaging Term Project Segmenting 2D Ultrasound Images using Seeded Region Growing David Boen #43589050 April

More information

Motivation. Gray Levels

Motivation. Gray Levels Motivation Image Intensity and Point Operations Dr. Edmund Lam Department of Electrical and Electronic Engineering The University of Hong ong A digital image is a matrix of numbers, each corresponding

More information

MR IMAGE SEGMENTATION

MR IMAGE SEGMENTATION MR IMAGE SEGMENTATION Prepared by : Monil Shah What is Segmentation? Partitioning a region or regions of interest in images such that each region corresponds to one or more anatomic structures Classification

More information

Edge and local feature detection - 2. Importance of edge detection in computer vision

Edge and local feature detection - 2. Importance of edge detection in computer vision Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature

More information

CHAPTER 2. Morphometry on rodent brains. A.E.H. Scheenstra J. Dijkstra L. van der Weerd

CHAPTER 2. Morphometry on rodent brains. A.E.H. Scheenstra J. Dijkstra L. van der Weerd CHAPTER 2 Morphometry on rodent brains A.E.H. Scheenstra J. Dijkstra L. van der Weerd This chapter was adapted from: Volumetry and other quantitative measurements to assess the rodent brain, In vivo NMR

More information

Digital Image Processing. Image Enhancement in the Frequency Domain

Digital Image Processing. Image Enhancement in the Frequency Domain Digital Image Processing Image Enhancement in the Frequency Domain Topics Frequency Domain Enhancements Fourier Transform Convolution High Pass Filtering in Frequency Domain Low Pass Filtering in Frequency

More information

IMAGE ENHANCEMENT in SPATIAL DOMAIN by Intensity Transformations

IMAGE ENHANCEMENT in SPATIAL DOMAIN by Intensity Transformations It makes all the difference whether one sees darkness through the light or brightness through the shadows David Lindsay IMAGE ENHANCEMENT in SPATIAL DOMAIN by Intensity Transformations Kalyan Kumar Barik

More information

Computational Medical Imaging Analysis Chapter 5: Processing and Analysis

Computational Medical Imaging Analysis Chapter 5: Processing and Analysis Computational Medical Imaging Analysis Chapter 5: Processing and Analysis Jun Zhang Laboratory for Computational Medical Imaging & Data Analysis Department of Computer Science University of Kentucky Lexington,

More information

Image processing in frequency Domain

Image processing in frequency Domain Image processing in frequency Domain Introduction to Frequency Domain Deal with images in: -Spatial domain -Frequency domain Frequency Domain In the frequency or Fourier domain, the value and location

More information

INTENSITY TRANSFORMATION AND SPATIAL FILTERING

INTENSITY TRANSFORMATION AND SPATIAL FILTERING 1 INTENSITY TRANSFORMATION AND SPATIAL FILTERING Lecture 3 Image Domains 2 Spatial domain Refers to the image plane itself Image processing methods are based and directly applied to image pixels Transform

More information

Outline 7/2/201011/6/

Outline 7/2/201011/6/ Outline Pattern recognition in computer vision Background on the development of SIFT SIFT algorithm and some of its variations Computational considerations (SURF) Potential improvement Summary 01 2 Pattern

More information

EXAM SOLUTIONS. Image Processing and Computer Vision Course 2D1421 Monday, 13 th of March 2006,

EXAM SOLUTIONS. Image Processing and Computer Vision Course 2D1421 Monday, 13 th of March 2006, School of Computer Science and Communication, KTH Danica Kragic EXAM SOLUTIONS Image Processing and Computer Vision Course 2D1421 Monday, 13 th of March 2006, 14.00 19.00 Grade table 0-25 U 26-35 3 36-45

More information

Lecture 4 Image Enhancement in Spatial Domain

Lecture 4 Image Enhancement in Spatial Domain Digital Image Processing Lecture 4 Image Enhancement in Spatial Domain Fall 2010 2 domains Spatial Domain : (image plane) Techniques are based on direct manipulation of pixels in an image Frequency Domain

More information

1 Introduction Motivation and Aims Functional Imaging Computational Neuroanatomy... 12

1 Introduction Motivation and Aims Functional Imaging Computational Neuroanatomy... 12 Contents 1 Introduction 10 1.1 Motivation and Aims....... 10 1.1.1 Functional Imaging.... 10 1.1.2 Computational Neuroanatomy... 12 1.2 Overview of Chapters... 14 2 Rigid Body Registration 18 2.1 Introduction.....

More information

Image representation. 1. Introduction

Image representation. 1. Introduction Image representation Introduction Representation schemes Chain codes Polygonal approximations The skeleton of a region Boundary descriptors Some simple descriptors Shape numbers Fourier descriptors Moments

More information

Image Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi

Image Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi Image Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi 1. Introduction The choice of a particular transform in a given application depends on the amount of

More information

EE795: Computer Vision and Intelligent Systems

EE795: Computer Vision and Intelligent Systems EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 WRI C225 Lecture 02 130124 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Basics Image Formation Image Processing 3 Intelligent

More information

6 credits. BMSC-GA Practical Magnetic Resonance Imaging II

6 credits. BMSC-GA Practical Magnetic Resonance Imaging II BMSC-GA 4428 - Practical Magnetic Resonance Imaging II 6 credits Course director: Ricardo Otazo, PhD Course description: This course is a practical introduction to image reconstruction, image analysis

More information

Digital Image Processing

Digital Image Processing Digital Image Processing Third Edition Rafael C. Gonzalez University of Tennessee Richard E. Woods MedData Interactive PEARSON Prentice Hall Pearson Education International Contents Preface xv Acknowledgments

More information

Image features. Image Features

Image features. Image Features Image features Image features, such as edges and interest points, provide rich information on the image content. They correspond to local regions in the image and are fundamental in many applications in

More information

Review on Image Segmentation Techniques and its Types

Review on Image Segmentation Techniques and its Types 1 Review on Image Segmentation Techniques and its Types Ritu Sharma 1, Rajesh Sharma 2 Research Scholar 1 Assistant Professor 2 CT Group of Institutions, Jalandhar. 1 rits_243@yahoo.in, 2 rajeshsharma1234@gmail.com

More information

Image Processing Fundamentals. Nicolas Vazquez Principal Software Engineer National Instruments

Image Processing Fundamentals. Nicolas Vazquez Principal Software Engineer National Instruments Image Processing Fundamentals Nicolas Vazquez Principal Software Engineer National Instruments Agenda Objectives and Motivations Enhancing Images Checking for Presence Locating Parts Measuring Features

More information

Digital Image Processing, 2nd ed. Digital Image Processing, 2nd ed. The principal objective of enhancement

Digital Image Processing, 2nd ed. Digital Image Processing, 2nd ed. The principal objective of enhancement Chapter 3 Image Enhancement in the Spatial Domain The principal objective of enhancement to process an image so that the result is more suitable than the original image for a specific application. Enhancement

More information

Non-rigid Image Registration

Non-rigid Image Registration Overview Non-rigid Image Registration Introduction to image registration - he goal of image registration - Motivation for medical image registration - Classification of image registration - Nonrigid registration

More information

Methods for data preprocessing

Methods for data preprocessing Methods for data preprocessing John Ashburner Wellcome Trust Centre for Neuroimaging, 12 Queen Square, London, UK. Overview Voxel-Based Morphometry Morphometry in general Volumetrics VBM preprocessing

More information

Image Enhancement in Spatial Domain. By Dr. Rajeev Srivastava

Image Enhancement in Spatial Domain. By Dr. Rajeev Srivastava Image Enhancement in Spatial Domain By Dr. Rajeev Srivastava CONTENTS Image Enhancement in Spatial Domain Spatial Domain Methods 1. Point Processing Functions A. Gray Level Transformation functions for

More information

DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING DS7201 ADVANCED DIGITAL IMAGE PROCESSING II M.E (C.S) QUESTION BANK UNIT I 1. Write the differences between photopic and scotopic vision? 2. What

More information

Digital Image Processing COSC 6380/4393

Digital Image Processing COSC 6380/4393 Digital Image Processing COSC 6380/4393 Lecture 21 Nov 16 th, 2017 Pranav Mantini Ack: Shah. M Image Processing Geometric Transformation Point Operations Filtering (spatial, Frequency) Input Restoration/

More information

Noise Model. Important Noise Probability Density Functions (Cont.) Important Noise Probability Density Functions

Noise Model. Important Noise Probability Density Functions (Cont.) Important Noise Probability Density Functions Others -- Noise Removal Techniques -- Edge Detection Techniques -- Geometric Operations -- Color Image Processing -- Color Spaces Xiaojun Qi Noise Model The principal sources of noise in digital images

More information

Applying Catastrophe Theory to Image Segmentation

Applying Catastrophe Theory to Image Segmentation Applying Catastrophe Theory to Image Segmentation Mohamad Raad, Majd Ghareeb, Ali Bazzi Department of computer and communications engineering Lebanese International University Beirut, Lebanon Abstract

More information

Feature Extraction and Image Processing, 2 nd Edition. Contents. Preface

Feature Extraction and Image Processing, 2 nd Edition. Contents. Preface , 2 nd Edition Preface ix 1 Introduction 1 1.1 Overview 1 1.2 Human and Computer Vision 1 1.3 The Human Vision System 3 1.3.1 The Eye 4 1.3.2 The Neural System 7 1.3.3 Processing 7 1.4 Computer Vision

More information

Fourier Transform in Image Processing. CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012)

Fourier Transform in Image Processing. CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012) Fourier Transform in Image Processing CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012) 1D: Common Transform Pairs Summary source FT Properties: Convolution See book DIP 4.2.5:

More information

Registration-Based Segmentation of Medical Images

Registration-Based Segmentation of Medical Images School of Computing National University of Singapore Graduate Research Paper Registration-Based Segmentation of Medical Images by Li Hao under guidance of A/Prof. Leow Wee Kheng July, 2006 Abstract Medical

More information

Volumetry of hypothalamic substructures by multimodal morphological image registration

Volumetry of hypothalamic substructures by multimodal morphological image registration Volumetry of hypothalamic substructures by multimodal morphological image registration Dominik Löchel 14.09.2011 Institute for Applied and Numerical Mathematics KIT University of the State of Baden-Württemberg

More information

Biomagnetic inverse problems:

Biomagnetic inverse problems: Biomagnetic inverse problems: Magnetic resonance electrical property tomography (MREPT) and magnetoencephalography (MEG) 2018 Aug. 16 The University of Tokyo Takaaki Nara 1 Contents Measurement What is

More information

Computer Vision and Graphics (ee2031) Digital Image Processing I

Computer Vision and Graphics (ee2031) Digital Image Processing I Computer Vision and Graphics (ee203) Digital Image Processing I Dr John Collomosse J.Collomosse@surrey.ac.uk Centre for Vision, Speech and Signal Processing University of Surrey Learning Outcomes After

More information

2D Rigid Registration of MR Scans using the 1d Binary Projections

2D Rigid Registration of MR Scans using the 1d Binary Projections 2D Rigid Registration of MR Scans using the 1d Binary Projections Panos D. Kotsas Abstract This paper presents the application of a signal intensity independent registration criterion for 2D rigid body

More information

Segmentation and Grouping

Segmentation and Grouping Segmentation and Grouping How and what do we see? Fundamental Problems ' Focus of attention, or grouping ' What subsets of pixels do we consider as possible objects? ' All connected subsets? ' Representation

More information

Schedule for Rest of Semester

Schedule for Rest of Semester Schedule for Rest of Semester Date Lecture Topic 11/20 24 Texture 11/27 25 Review of Statistics & Linear Algebra, Eigenvectors 11/29 26 Eigenvector expansions, Pattern Recognition 12/4 27 Cameras & calibration

More information

Advanced Video Content Analysis and Video Compression (5LSH0), Module 4

Advanced Video Content Analysis and Video Compression (5LSH0), Module 4 Advanced Video Content Analysis and Video Compression (5LSH0), Module 4 Visual feature extraction Part I: Color and texture analysis Sveta Zinger Video Coding and Architectures Research group, TU/e ( s.zinger@tue.nl

More information

CoE4TN4 Image Processing

CoE4TN4 Image Processing CoE4TN4 Image Processing Chapter 11 Image Representation & Description Image Representation & Description After an image is segmented into regions, the regions are represented and described in a form suitable

More information

Steen Moeller Center for Magnetic Resonance research University of Minnesota

Steen Moeller Center for Magnetic Resonance research University of Minnesota Steen Moeller Center for Magnetic Resonance research University of Minnesota moeller@cmrr.umn.edu Lot of material is from a talk by Douglas C. Noll Department of Biomedical Engineering Functional MRI Laboratory

More information

Image Processing. Application area chosen because it has very good parallelism and interesting output.

Image Processing. Application area chosen because it has very good parallelism and interesting output. Chapter 11 Slide 517 Image Processing Application area chosen because it has very good parallelism and interesting output. Low-level Image Processing Operates directly on stored image to improve/enhance

More information

Topic 6 Representation and Description

Topic 6 Representation and Description Topic 6 Representation and Description Background Segmentation divides the image into regions Each region should be represented and described in a form suitable for further processing/decision-making Representation

More information

Registration Techniques

Registration Techniques EMBO Practical Course on Light Sheet Microscopy Junior-Prof. Dr. Olaf Ronneberger Computer Science Department and BIOSS Centre for Biological Signalling Studies University of Freiburg Germany O. Ronneberger,

More information

An Intuitive Explanation of Fourier Theory

An Intuitive Explanation of Fourier Theory An Intuitive Explanation of Fourier Theory Steven Lehar slehar@cns.bu.edu Fourier theory is pretty complicated mathematically. But there are some beautifully simple holistic concepts behind Fourier theory

More information

Biophysical Techniques (BPHS 4090/PHYS 5800)

Biophysical Techniques (BPHS 4090/PHYS 5800) Biophysical Techniques (BPHS 4090/PHYS 5800) Instructors: Prof. Christopher Bergevin (cberge@yorku.ca) Schedule: MWF 1:30-2:30 (CB 122) Website: http://www.yorku.ca/cberge/4090w2017.html York University

More information

Medical Image Fusion using Rayleigh Contrast Limited Adaptive Histogram Equalization and Ant Colony Edge Method

Medical Image Fusion using Rayleigh Contrast Limited Adaptive Histogram Equalization and Ant Colony Edge Method Medical Image Fusion using Rayleigh Contrast Limited Adaptive Histogram Equalization and Ant Colony Edge Method Ramandeep 1, Rajiv Kamboj 2 1 Student, M. Tech (ECE), Doon Valley Institute of Engineering

More information

82 REGISTRATION OF RETINOGRAPHIES

82 REGISTRATION OF RETINOGRAPHIES 82 REGISTRATION OF RETINOGRAPHIES 3.3 Our method Our method resembles the human approach to image matching in the sense that we also employ as guidelines features common to both images. It seems natural

More information

HST.582J / 6.555J / J Biomedical Signal and Image Processing Spring 2007

HST.582J / 6.555J / J Biomedical Signal and Image Processing Spring 2007 MIT OpenCourseWare http://ocw.mit.edu HST.582J / 6.555J / 16.456J Biomedical Signal and Image Processing Spring 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.

More information

Chapter 11 Image Processing

Chapter 11 Image Processing Chapter Image Processing Low-level Image Processing Operates directly on a stored image to improve or enhance it. Stored image consists of a two-dimensional array of pixels (picture elements): Origin (0,

More information

Snakes, level sets and graphcuts. (Deformable models)

Snakes, level sets and graphcuts. (Deformable models) INSTITUTE OF INFORMATION AND COMMUNICATION TECHNOLOGIES BULGARIAN ACADEMY OF SCIENCE Snakes, level sets and graphcuts (Deformable models) Centro de Visión por Computador, Departament de Matemàtica Aplicada

More information

Examination in Image Processing

Examination in Image Processing Umeå University, TFE Ulrik Söderström 203-03-27 Examination in Image Processing Time for examination: 4.00 20.00 Please try to extend the answers as much as possible. Do not answer in a single sentence.

More information

Lecture 2: 2D Fourier transforms and applications

Lecture 2: 2D Fourier transforms and applications Lecture 2: 2D Fourier transforms and applications B14 Image Analysis Michaelmas 2017 Dr. M. Fallon Fourier transforms and spatial frequencies in 2D Definition and meaning The Convolution Theorem Applications

More information

Motivation. Intensity Levels

Motivation. Intensity Levels Motivation Image Intensity and Point Operations Dr. Edmund Lam Department of Electrical and Electronic Engineering The University of Hong ong A digital image is a matrix of numbers, each corresponding

More information

Mutual Information Based Methods to Localize Image Registration

Mutual Information Based Methods to Localize Image Registration Mutual Information Based Methods to Localize Image Registration by Kathleen P. Wilkie A thesis presented to the University of Waterloo in fulfilment of the thesis requirement for the degree of Master of

More information

Introduction to Medical Image Processing

Introduction to Medical Image Processing Introduction to Medical Image Processing Δ Essential environments of a medical imaging system Subject Image Analysis Energy Imaging System Images Image Processing Feature Images Image processing may be

More information

Analysis of Binary Images

Analysis of Binary Images Analysis of Binary Images Introduction to Computer Vision CSE 52 Lecture 7 CSE52, Spr 07 The appearance of colors Color appearance is strongly affected by (at least): Spectrum of lighting striking the

More information

Boundary descriptors. Representation REPRESENTATION & DESCRIPTION. Descriptors. Moore boundary tracking

Boundary descriptors. Representation REPRESENTATION & DESCRIPTION. Descriptors. Moore boundary tracking Representation REPRESENTATION & DESCRIPTION After image segmentation the resulting collection of regions is usually represented and described in a form suitable for higher level processing. Most important

More information

identified and grouped together.

identified and grouped together. Segmentation ti of Images SEGMENTATION If an image has been preprocessed appropriately to remove noise and artifacts, segmentation is often the key step in interpreting the image. Image segmentation is

More information

Texture Segmentation and Classification in Biomedical Image Processing

Texture Segmentation and Classification in Biomedical Image Processing Texture Segmentation and Classification in Biomedical Image Processing Aleš Procházka and Andrea Gavlasová Department of Computing and Control Engineering Institute of Chemical Technology in Prague Technická

More information

Biomedical Image Analysis. Point, Edge and Line Detection

Biomedical Image Analysis. Point, Edge and Line Detection Biomedical Image Analysis Point, Edge and Line Detection Contents: Point and line detection Advanced edge detection: Canny Local/regional edge processing Global processing: Hough transform BMIA 15 V. Roth

More information

Methodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion

Methodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion Methodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion Mattias P. Heinrich Julia A. Schnabel, Mark Jenkinson, Sir Michael Brady 2 Clinical

More information

Sampling and Reconstruction

Sampling and Reconstruction Sampling and Reconstruction Sampling and Reconstruction Sampling and Spatial Resolution Spatial Aliasing Problem: Spatial aliasing is insufficient sampling of data along the space axis, which occurs because

More information

Binary Image Processing. Introduction to Computer Vision CSE 152 Lecture 5

Binary Image Processing. Introduction to Computer Vision CSE 152 Lecture 5 Binary Image Processing CSE 152 Lecture 5 Announcements Homework 2 is due Apr 25, 11:59 PM Reading: Szeliski, Chapter 3 Image processing, Section 3.3 More neighborhood operators Binary System Summary 1.

More information

An Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy

An Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy An Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy Chenyang Xu 1, Siemens Corporate Research, Inc., Princeton, NJ, USA Xiaolei Huang,

More information

Short Survey on Static Hand Gesture Recognition

Short Survey on Static Hand Gesture Recognition Short Survey on Static Hand Gesture Recognition Huu-Hung Huynh University of Science and Technology The University of Danang, Vietnam Duc-Hoang Vo University of Science and Technology The University of

More information

2D Image Processing INFORMATIK. Kaiserlautern University. DFKI Deutsches Forschungszentrum für Künstliche Intelligenz

2D Image Processing INFORMATIK. Kaiserlautern University.   DFKI Deutsches Forschungszentrum für Künstliche Intelligenz 2D Image Processing - Filtering Prof. Didier Stricker Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche Intelligenz http://av.dfki.de 1 What is image filtering?

More information

Role of Parallel Imaging in High Field Functional MRI

Role of Parallel Imaging in High Field Functional MRI Role of Parallel Imaging in High Field Functional MRI Douglas C. Noll & Bradley P. Sutton Department of Biomedical Engineering, University of Michigan Supported by NIH Grant DA15410 & The Whitaker Foundation

More information

Reconstruction in CT and relation to other imaging modalities

Reconstruction in CT and relation to other imaging modalities Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 16, 2015 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound

More information

Ulrik Söderström 16 Feb Image Processing. Segmentation

Ulrik Söderström 16 Feb Image Processing. Segmentation Ulrik Söderström ulrik.soderstrom@tfe.umu.se 16 Feb 2011 Image Processing Segmentation What is Image Segmentation? To be able to extract information from an image it is common to subdivide it into background

More information

EXAM SOLUTIONS. Computer Vision Course 2D1420 Thursday, 11 th of march 2003,

EXAM SOLUTIONS. Computer Vision Course 2D1420 Thursday, 11 th of march 2003, Numerical Analysis and Computer Science, KTH Danica Kragic EXAM SOLUTIONS Computer Vision Course 2D1420 Thursday, 11 th of march 2003, 8.00 13.00 Exercise 1 (5*2=10 credits) Answer at most 5 of the following

More information

Volume 2, Issue 9, September 2014 ISSN

Volume 2, Issue 9, September 2014 ISSN Fingerprint Verification of the Digital Images by Using the Discrete Cosine Transformation, Run length Encoding, Fourier transformation and Correlation. Palvee Sharma 1, Dr. Rajeev Mahajan 2 1M.Tech Student

More information

Multimodality Imaging for Tumor Volume Definition in Radiation Oncology

Multimodality Imaging for Tumor Volume Definition in Radiation Oncology 81 There are several commercial and academic software tools that support different segmentation algorithms. In general, commercial software packages have better implementation (with a user-friendly interface

More information

Texture Segmentation by using Haar Wavelets and K-means Algorithm

Texture Segmentation by using Haar Wavelets and K-means Algorithm Texture Segmentation by using Haar Wavelets and K-means Algorithm P. Ashok Babu Associate Professor, Narsimha Reddy Engineering College, Hyderabad, A.P., INDIA, ashokbabup2@gmail.com Dr. K. V. S. V. R.

More information

Deformable Segmentation using Sparse Shape Representation. Shaoting Zhang

Deformable Segmentation using Sparse Shape Representation. Shaoting Zhang Deformable Segmentation using Sparse Shape Representation Shaoting Zhang Introduction Outline Our methods Segmentation framework Sparse shape representation Applications 2D lung localization in X-ray 3D

More information

Image Processing. Bilkent University. CS554 Computer Vision Pinar Duygulu

Image Processing. Bilkent University. CS554 Computer Vision Pinar Duygulu Image Processing CS 554 Computer Vision Pinar Duygulu Bilkent University Today Image Formation Point and Blob Processing Binary Image Processing Readings: Gonzalez & Woods, Ch. 3 Slides are adapted from

More information

Building a Panorama. Matching features. Matching with Features. How do we build a panorama? Computational Photography, 6.882

Building a Panorama. Matching features. Matching with Features. How do we build a panorama? Computational Photography, 6.882 Matching features Building a Panorama Computational Photography, 6.88 Prof. Bill Freeman April 11, 006 Image and shape descriptors: Harris corner detectors and SIFT features. Suggested readings: Mikolajczyk

More information