Topics. View Aliasing. CT Sampling Requirements. Sampling Requirements in CT Sampling Theory Aliasing
|
|
- Morgan Paul
- 5 years ago
- Views:
Transcription
1 Topics Bioengineering 280A Principles of Biomedical Imaging Sampling Requirements in CT Sampling Theory Aliasing Fall Quarter 2010 CT/Fourier Lecture 5 CT Sampling Requirements View Aliasing What should the size of the detectors be? How many detectors do we need? How many views do we need? Kak and Slaney 1
2 Aliasing Artifacts Object Effect of Noise Aliasing due to insufficient number of detectors Aliasing due to insufficient number of views Kak and Slaney Analog vs. Digital The Process of Sampling The Analog World: Continuous time/space, continuous valued signals or images, e.g. vinyl records, photographs, -ray films. g() The Digital World: Discrete time/space, discrete-valued signals or images, e.g. CD-Roms, DVDs, digital photos, digital -rays, CT, MRI, ultrasound. sample g[n]=g(n ) 2
3 Questions Sampling in the Time Domain How finely do we need to sample? What happens if we don t sample finely enough? Can we reconstruct the original signal or image from its samples? Sampling in Image Space Sampling in k-space 3
4 Sampling in k-space Comb Function comb() = δ( n) n= Other names: Impulse train, bed of nails, shah function. Scaled Comb Function comb = n= δ( n) = δ( n ) n= = δ( n) n= 1D spatial sampling g S () = g() 1 comb Recall the sifting property = g() δ( n) n= = g(n)δ( n) n= g()δ( a) = g(a) But we can also write g(a)δ( a) = g(a) δ( a) = g(a) So, g()δ( a) = g(a)δ( a) 4
5 1D spatial sampling g() comb(/)/ Fourier Transform of comb() F[ comb() ] = comb(k ) = δ(k n) n= g S () F 1 comb( ) = 1 comb(k ) = δ(k n) n= = 1 n= δ(k n ) Fourier Transform of comb(/ ) comb(/ )/ F 1/ comb(k ) 1/ k Fourier Transform of g S () [ ] = F g() 1 comb F g S () = G(k ) F 1 comb = G(k ) 1 = 1 = 1 n= n= n= δ k n G(k ) δ k n G k n 5
6 Fourier Transform of g S () G(k ) Nyquist Condition G(k ) k -B B k G S (k ) 1/ G S (k ) k k =1/ 1/ To avoid overlap, we require that 1/>2B or > 2B where =1/ is the sampling frequency Eample Reconstruction from Samples Assume that the highest spatial frequency in an object is B = 2 cm -1. Thus, smallest spatial period is 0.5 cm.. Nyquist theorem says we need to sample with < 1/2B = 0.25 cm =1/ G S (k ) Multiply by (1/ ) rect(k / ) This corresponds to 2 samples per spatial period. (1/ ) G S (k )rect(k / ) =G(k ) 6
7 Eample Cosine Reconstruction cos(2π ) >2 Reconstruction from Samples If the Nyquist condition G ˆ S (k ) = 1 is met, then G S (k )rect(k / ) = G(k ) And the signal can be reconstructed by convolving the sample with a sinc function =2 g ˆ S () = g S () sinc(k s ) = g(nδx)δ( nδx) sinc(k s ) n= = g(nδx) sinc(k s ( n)) n= g() Reconstruction from Samples g S () Cosine Eample with =2 Sample at ˆ g S () sinc(k s ) = sinc( / ) 7
8 Eample with K s =4 Eample with K s =8 Aliasing Aliasing Eample -B B G(k ) k Aliasing occurs when the Nyquist condition is not satisfied. This occurs for 2B 8
9 Aliasing Eample Aliasing Eample cos(2π ) cos(2π ) = 2 > > Eample Detector Sampling Requirements 1. Consider the function g() = cos 2 ( 2π ). Sketch this function. You sample this signal in the spatial domain with a sampling rate =1/ (e.g. samples spaced at intervals of ). What is the minimum sampling rate that you can use without aliasing? Give an intuitive eplanation for your answer. Sampling interval Δr Beamwidth Δs 9
10 Smoothing of Projection Smoothing of Projection Projection g s (l,θ) = rect(l /Δs) g( l,θ) G s (k,θ) = Δssinc(k Δs)G(k,θ) Beam Width 2/(Δs) Smoothed Projection Sampling Requirements View Aliasing Smoothed Projection Detectors Δr Δs/2 Sampled Smooth Projection Kak and Slaney 10
11 View Sampling Requirements View Sampling -- how many views? Basic idea is that to make the maimum angular sampling the same as the projection sampling. πfov N views = Δr N views,360 = πfov Δr N views,180 = πn proj 2 = πn proj (for 360 degrees) (for 180 degrees) beamwidth Δs = 1 mm Eample Field of View (FOV) = 50 cm Δr = Δs/2 = 0.5 mm 500 mm/ 0.5 mm = N = 1000 detector samples π *N = 3146 views per 360 degrees 1500 views per 180 degrees CT "Rule of Thumb" N view = N det ectors = N piels 11
Topics. Projections. Review Filtered Backprojection Fan Beam Spiral CT Applications. Bioengineering 280A Principles of Biomedical Imaging
Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 24 X-Rays/CT Lecture 2 Topics Review Filtered Backprojection Fan Beam Spiral CT Applications s I θ (r) = I exp µ(x, y)ds Lr,θ = I exp µ(rcosθ
More informationTopics. Review Filtered Backprojection Fan Beam Spiral CT Applications. Bioengineering 280A Principles of Biomedical Imaging
Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 24 X-Rays/CT Lecture 2 Topics Review Filtered Backprojection Fan Beam Spiral CT Applications Projections I θ (r) = I exp Lr,θ µ(x,y)ds =
More informationAdvanced Computer Graphics. Aliasing. Matthias Teschner. Computer Science Department University of Freiburg
Advanced Computer Graphics Aliasing Matthias Teschner Computer Science Department University of Freiburg Outline motivation Fourier analysis filtering sampling reconstruction / aliasing antialiasing University
More informationFourier analysis and sampling theory
Reading Required: Shirley, Ch. 9 Recommended: Fourier analysis and sampling theory Ron Bracewell, The Fourier Transform and Its Applications, McGraw-Hill. Don P. Mitchell and Arun N. Netravali, Reconstruction
More informationSampling and Reconstruction
Page 1 Sampling and Reconstruction The sampling and reconstruction process Real world: continuous Digital world: discrete Basic signal processing Fourier transforms The convolution theorem The sampling
More informationCS354 Computer Graphics Sampling and Aliasing
Slide Credit: http://www.siggraph.org/education/materials/hypergraph/aliasing/alias0.htm CS354 Computer Graphics Sampling and Aliasing Qixing Huang February 12th 2018 Sampling and Aliasing Display is discrete
More informationTheoretically Perfect Sensor
Sampling 1/60 Sampling The ray tracer samples the geometry, only gathering information from the parts of the world that interact with a finite number of rays In contrast, a scanline renderer can push all
More informationReading. 2. Fourier analysis and sampling theory. Required: Watt, Section 14.1 Recommended:
Reading Required: Watt, Section 14.1 Recommended: 2. Fourier analysis and sampling theory Ron Bracewell, The Fourier Transform and Its Applications, McGraw-Hill. Don P. Mitchell and Arun N. Netravali,
More informationAliasing and Antialiasing. ITCS 4120/ Aliasing and Antialiasing
Aliasing and Antialiasing ITCS 4120/5120 1 Aliasing and Antialiasing What is Aliasing? Errors and Artifacts arising during rendering, due to the conversion from a continuously defined illumination field
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 informationImage Acquisition Systems
Image Acquisition Systems Goals and Terminology Conventional Radiography Axial Tomography Computer Axial Tomography (CAT) Magnetic Resonance Imaging (MRI) PET, SPECT Ultrasound Microscopy Imaging ITCS
More informationTheoretically Perfect Sensor
Sampling 1/67 Sampling The ray tracer samples the geometry, only gathering information from the parts of the world that interact with a finite number of rays In contrast, a scanline renderer can push all
More informationIndex. aliasing artifacts and noise in CT images, 200 measurement of projection data, nondiffracting
Index Algebraic equations solution by Kaczmarz method, 278 Algebraic reconstruction techniques, 283-84 sequential, 289, 293 simultaneous, 285-92 Algebraic techniques reconstruction algorithms, 275-96 Algorithms
More informationTo Do. Advanced Computer Graphics. Sampling and Reconstruction. Outline. Sign up for Piazza
Advanced Computer Graphics CSE 63 [Spring 207], Lecture 3 Ravi Ramamoorthi http://www.cs.ucsd.edu/~ravir Sign up for Piazza To Do Assignment, Due Apr 28. Anyone need help finding partners? Any issues with
More informationComputed Tomography (Part 2)
EL-GY 6813 / BE-GY 603 / G16.446 Medical Imaging Computed Tomography (Part ) Yao Wang Polytechnic School of Engineering New York University, Brooklyn, NY 1101 Based on Prince and Links, Medical Imaging
More informationLecture 6 Basic Signal Processing
Lecture 6 Basic Signal Processing Copyright c1996, 1997 by Pat Hanrahan Motivation Many aspects of computer graphics and computer imagery differ from aspects of conventional graphics and imagery because
More informationKey properties of local features
Key properties of local features Locality, robust against occlusions Must be highly distinctive, a good feature should allow for correct object identification with low probability of mismatch Easy to etract
More informationOutline. Sampling and Reconstruction. Sampling and Reconstruction. Foundations of Computer Graphics (Fall 2012)
Foundations of Computer Graphics (Fall 2012) CS 184, Lectures 19: Sampling and Reconstruction http://inst.eecs.berkeley.edu/~cs184 Outline Basic ideas of sampling, reconstruction, aliasing Signal processing
More informationSampling, Aliasing, & Mipmaps
Sampling, Aliasing, & Mipmaps Last Time? Monte-Carlo Integration Importance Sampling Ray Tracing vs. Path Tracing source hemisphere Sampling sensitive to choice of samples less sensitive to choice of samples
More informationComputer Graphics. Sampling Theory & Anti-Aliasing. Philipp Slusallek
Computer Graphics Sampling Theory & Anti-Aliasing Philipp Slusallek Dirac Comb (1) Constant & δ-function flash Comb/Shah function 2 Dirac Comb (2) Constant & δ-function Duality f(x) = K F(ω) = K (ω) And
More informationScaled representations
Scaled representations Big bars (resp. spots, hands, etc.) and little bars are both interesting Stripes and hairs, say Inefficient to detect big bars with big filters And there is superfluous detail in
More informationEECS 556 Image Processing W 09. Image enhancement. Smoothing and noise removal Sharpening filters
EECS 556 Image Processing W 09 Image enhancement Smoothing and noise removal Sharpening filters What is image processing? Image processing is the application of 2D signal processing methods to images Image
More informationThe Fundamental Theorem of Calculus Using the Rule of Three
The Fundamental Theorem of Calculus Using the Rule of Three A. Approimations with Riemann sums. The area under a curve can be approimated through the use of Riemann (or rectangular) sums: n Area f ( k
More informationOutline. Foundations of Computer Graphics (Spring 2012)
Foundations of Computer Graphics (Spring 2012) CS 184, Lectures 19: Sampling and Reconstruction http://inst.eecs.berkeley.edu/~cs184 Basic ideas of sampling, reconstruction, aliasing Signal processing
More informationAssignment 2. Due Feb 3, 2012
EE225E/BIOE265 Spring 2012 Principles of MRI Miki Lustig Assignment 2 Due Feb 3, 2012 1. Read Nishimura Ch. 3 2. Non-Uniform Sampling. A student has an assignment to monitor the level of Hetch-Hetchi reservoir
More informationComputer Graphics. Texture Filtering & Sampling Theory. Hendrik Lensch. Computer Graphics WS07/08 Texturing
Computer Graphics Texture Filtering & Sampling Theory Hendrik Lensch Overview Last time Texture Parameterization Procedural Shading Today Texturing Filtering 2D Texture Mapping Forward mapping Object surface
More informationImage Filtering, Warping and Sampling
Image Filtering, Warping and Sampling Connelly Barnes CS 4810 University of Virginia Acknowledgement: slides by Jason Lawrence, Misha Kazhdan, Allison Klein, Tom Funkhouser, Adam Finkelstein and David
More informationBasics. Sampling and Reconstruction. Sampling and Reconstruction. Outline. (Spatial) Aliasing. Advanced Computer Graphics (Fall 2010)
Advanced Computer Graphics (Fall 2010) CS 283, Lecture 3: Sampling and Reconstruction Ravi Ramamoorthi http://inst.eecs.berkeley.edu/~cs283/fa10 Some slides courtesy Thomas Funkhouser and Pat Hanrahan
More informationTopics. 2D Signal. a b. c d. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2004 Lecture 4 2D Fourier Transforms
Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2004 Lecture 4 2D Fourier Transforms Topics. 2D Signal Representations 2. 2D Fourier Transform 3. Transform Pairs 4. FT Properties 2D Signal
More informationLecture 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 informationEnhancement Image Quality of CT Using Single Slice Spiral Technique
Enhancement Image Quality of CT Using Single Slice Spiral Technique Doaa. N. Al Sheack 1 and Dr.Mohammed H. Ali Al Hayani 2 1 2 Electronic and Communications Engineering Department College of Engineering,
More informationImage Reconstruction from Projection
Image Reconstruction from Projection Reconstruct an image from a series of projections X-ray computed tomography (CT) Computed tomography is a medical imaging method employing tomography where digital
More informationSampling, Aliasing, & Mipmaps
Sampling, Aliasing, & Mipmaps Last Time? Monte-Carlo Integration Importance Sampling Ray Tracing vs. Path Tracing source hemisphere What is a Pixel? Sampling & Reconstruction Filters in Computer Graphics
More informationThe main goal of Computer Graphics is to generate 2D images 2D images are continuous 2D functions (or signals)
Motivation The main goal of Computer Graphics is to generate 2D images 2D images are continuous 2D functions (or signals) monochrome f(x,y) or color r(x,y), g(x,y), b(x,y) These functions are represented
More informationFilterbanks and transforms
Filterbanks and transforms Sources: Zölzer, Digital audio signal processing, Wiley & Sons. Saramäki, Multirate signal processing, TUT course. Filterbanks! Introduction! Critical sampling, half-band filter!
More informationMulti-slice CT Image Reconstruction Jiang Hsieh, Ph.D.
Multi-slice CT Image Reconstruction Jiang Hsieh, Ph.D. Applied Science Laboratory, GE Healthcare Technologies 1 Image Generation Reconstruction of images from projections. textbook reconstruction advanced
More informationMotivation. The main goal of Computer Graphics is to generate 2D images. 2D images are continuous 2D functions (or signals)
Motivation The main goal of Computer Graphics is to generate 2D images 2D images are continuous 2D functions (or signals) monochrome f(x,y) or color r(x,y), g(x,y), b(x,y) These functions are represented
More informationComputer 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 informationMRI Physics II: Gradients, Imaging
MRI Physics II: Gradients, Imaging Douglas C., Ph.D. Dept. of Biomedical Engineering University of Michigan, Ann Arbor Magnetic Fields in MRI B 0 The main magnetic field. Always on (0.5-7 T) Magnetizes
More informationBME 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 informationPhysics 2c Lecture 25. Chapter 37 Interference & Diffraction
Physics 2c Lecture 25 Chapter 37 Interference & Diffraction Outlook for rest of quarter Today: finish chapter 37 Tomorrow & Friday: E&M waves (Chapter 34) Next Monday, June 4 th : Quiz 8 on Chapter 37
More informationCh. 7.4, 7.6, 7.7: Complex Numbers, Polar Coordinates, ParametricFall equations / 17
Ch. 7.4, 7.6, 7.7: Complex Numbers, Polar Coordinates, Parametric equations Johns Hopkins University Fall 2014 Ch. 7.4, 7.6, 7.7: Complex Numbers, Polar Coordinates, ParametricFall equations 2014 1 / 17
More informationEE123 Digital Signal Processing
Multi-Dimensional Signals EE23 Digital Signal Processing Lecture Our world is more complex than D Images: f(x,y) Videos: f(x,y,t) Dynamic 3F scenes: f(x,y,z,t) Medical Imaging 3D Video Computer Graphics
More informationOutline. Visualization Discretization Sampling Quantization Representation Continuous Discrete. Noise
Fundamentals Data Outline Visualization Discretization Sampling Quantization Representation Continuous Discrete Noise 2 Data Data : Function dependent on one or more variables. Example Audio (1D) - depends
More informationconvolution shift invariant linear system Fourier Transform Aliasing and sampling scale representation edge detection corner detection
COS 429: COMPUTER VISON Linear Filters and Edge Detection convolution shift invariant linear system Fourier Transform Aliasing and sampling scale representation edge detection corner detection Reading:
More informationIntroduction to Medical Imaging. Lecture 6: X-Ray Computed Tomography. CT number (in HU) = Overview. Klaus Mueller
Overview Introduction to Medical Imaging Lecture 6: X-Ray Computed Tomography Scanning: rotate source-detector pair around the patient Klaus Mueller data Computer Science Department Stony Brook University
More informationModern CT system generations Measurement of attenuation
CT reconstruction repetition & hints Reconstruction in CT and hints to the assignments Jørgen Arendt Jensen October 4, 16 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering
More informationReconstruction in CT and hints to the assignments
Reconstruction in CT and hints to the assignments Jørgen Arendt Jensen October 24, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging
More informationImage Reconstruction 3 Fully 3D Reconstruction
Image Reconstruction 3 Fully 3D Reconstruction Thomas Bortfeld Massachusetts General Hospital, Radiation Oncology, HMS HST.S14, February 25, 2013 Thomas Bortfeld (MGH, HMS, Rad. Onc.) Image Reconstruction
More informationXI Signal-to-Noise (SNR)
XI Signal-to-Noise (SNR) Lecture notes by Assaf Tal n(t) t. Noise. Characterizing Noise Noise is a random signal that gets added to all of our measurements. In D it looks like this: while in D
More informationSynthesis Imaging. Claire Chandler, Sanjay Bhatnagar NRAO/Socorro
Synthesis Imaging Claire Chandler, Sanjay Bhatnagar NRAO/Socorro Michelson Summer Workshop Caltech, July 24-28, 2006 Synthesis Imaging 2 Based on the van Cittert-Zernike theorem: The complex visibility
More information2D Fan Beam Reconstruction 3D Cone Beam Reconstruction
2D Fan Beam Reconstruction 3D Cone Beam Reconstruction Mario Koerner March 17, 2006 1 2D Fan Beam Reconstruction Two-dimensional objects can be reconstructed from projections that were acquired using parallel
More informationFrequency analysis, pyramids, texture analysis, applications (face detection, category recognition)
Frequency analysis, pyramids, texture analysis, applications (face detection, category recognition) Outline Measuring frequencies in images: Definitions, properties Sampling issues Relation with Gaussian
More information2D Fan Beam Reconstruction 3D Cone Beam Reconstruction. Mario Koerner
2D Fan Beam Reconstruction 3D Cone Beam Reconstruction Mario Koerner Moscow-Bavarian Joint Advanced Student School 2006 March 19 2006 to March 29 2006 Overview 2D Fan Beam Reconstruction Shortscan Reconstruction
More informationTomographic reconstruction for multiguidestar. Don Gavel, Jose Milovich LLNL
Tomographic reconstruction for multiguidestar adaptive optics Don Gavel, Jose Milovich LLNL Tomographic Reconstruction Fundamental questions for the architecture of AO for etremely large telescopes Number
More informationIntroduction to Topics in Machine Learning
Introduction to Topics in Machine Learning Namrata Vaswani Department of Electrical and Computer Engineering Iowa State University Namrata Vaswani 1/ 27 Compressed Sensing / Sparse Recovery: Given y :=
More informationBME I5000: Biomedical Imaging
BME I5000: Biomedical Imaging Lecture 1 Introduction Lucas C. Parra, parra@ccny.cuny.edu 1 Content Topics: Physics of medial imaging modalities (blue) Digital Image Processing (black) Schedule: 1. Introduction,
More informationImage Sampling & Quantisation
Image Sampling & Quantisation Biomedical Image Analysis Prof. Dr. Philippe Cattin MIAC, University of Basel Contents 1 Motivation 2 Sampling Introduction and Motivation Sampling Example Quantisation Example
More informationChapter 9. Coherence
Chapter 9. Coherence Last Lecture Michelson Interferometer Variations of the Michelson Interferometer Fabry-Perot interferometer This Lecture Fourier analysis Temporal coherence and line width Partial
More informationBiomedical 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 informationCoding and Modulation in Cameras
Mitsubishi Electric Research Laboratories Raskar 2007 Coding and Modulation in Cameras Ramesh Raskar with Ashok Veeraraghavan, Amit Agrawal, Jack Tumblin, Ankit Mohan Mitsubishi Electric Research Labs
More informationSparse Reconstruction / Compressive Sensing
Sparse Reconstruction / Compressive Sensing Namrata Vaswani Department of Electrical and Computer Engineering Iowa State University Namrata Vaswani Sparse Reconstruction / Compressive Sensing 1/ 20 The
More informationCentral Slice Theorem
Central Slice Theorem Incident X-rays y f(x,y) R x r x Detected p(, x ) The thick line is described by xcos +ysin =R Properties of Fourier Transform F [ f ( x a)] F [ f ( x)] e j 2 a Spatial Domain Spatial
More informationDigital Image Processing
Digital Image Processing SPECIAL TOPICS CT IMAGES Hamid R. Rabiee Fall 2015 What is an image? 2 Are images only about visual concepts? We ve already seen that there are other kinds of image. In this lecture
More informationImage preprocessing in spatial domain
Image preprocessing in spatial domain Sampling theorem, aliasing, interpolation, geometrical transformations Revision:.4, dated: May 25, 26 Tomáš Svoboda Czech Technical University, Faculty of Electrical
More informationPhase problem and the Radon transform
Phase problem and the Radon transform Andrei V. Bronnikov Bronnikov Algorithms The Netherlands The Radon transform and applications Inverse problem of phase-contrast CT Fundamental theorem Image reconstruction
More informationImage preprocessing in spatial domain
Image preprocessing in spatial domain Sampling theorem, aliasing, interpolation, geometrical transformations Revision:.3, dated: December 7, 25 Tomáš Svoboda Czech Technical University, Faculty of Electrical
More informationA/D Converter. Sampling. Figure 1.1: Block Diagram of a DSP System
CHAPTER 1 INTRODUCTION Digital signal processing (DSP) technology has expanded at a rapid rate to include such diverse applications as CDs, DVDs, MP3 players, ipods, digital cameras, digital light processing
More informationPyramid Coding and Subband Coding
Pyramid Coding and Subband Coding Predictive pyramids Transform pyramids Subband coding Perfect reconstruction filter banks Quadrature mirror filter banks Octave band splitting Transform coding as a special
More informationDigital Image Processing. Image Enhancement - Filtering
Digital Image Processing Image Enhancement - Filtering Derivative Derivative is defined as a rate of change. Discrete Derivative Finite Distance Example Derivatives in 2-dimension Derivatives of Images
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 informationMEDICAL 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 informationComputer Vision I. Announcements. Fourier Tansform. Efficient Implementation. Edge and Corner Detection. CSE252A Lecture 13.
Announcements Edge and Corner Detection HW3 assigned CSE252A Lecture 13 Efficient Implementation Both, the Box filter and the Gaussian filter are separable: First convolve each row of input image I with
More informationApplications of Image Filters
02/04/0 Applications of Image Filters Computer Vision CS 543 / ECE 549 University of Illinois Derek Hoiem Review: Image filtering g[, ] f [.,.] h[.,.] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 90
More informationMéthodes d imagerie pour les écoulements et le CND
Méthodes d imagerie pour les écoulements et le CND Journée scientifique FED3G CEA LIST/Lab Imagerie Tomographie et Traitement Samuel Legoupil 15 juin 2012 2D/3D imaging tomography Example Petrochemical
More informationImage Coding and Data Compression
Image Coding and Data Compression Biomedical Images are of high spatial resolution and fine gray-scale quantisiation Digital mammograms: 4,096x4,096 pixels with 12bit/pixel 32MB per image Volume data (CT
More informationNon-Stationary CT Image Noise Spectrum Analysis
Non-Stationary CT Image Noise Spectrum Analysis Michael Balda, Björn J. Heismann,, Joachim Hornegger Pattern Recognition Lab, Friedrich-Alexander-Universität Erlangen Siemens Healthcare, Erlangen michael.balda@informatik.uni-erlangen.de
More informationDrawing a Triangle (and an introduction to sampling)
Lecture 4: Drawing a Triangle (and an introduction to sampling) Computer Graphics CMU 15-462/15-662, Spring 2017 Assignment 1 is out! https://15462-s17.github.io/asst1_drawsvg/ Let s draw some triangles
More informationFinal Review. Image Processing CSE 166 Lecture 18
Final Review Image Processing CSE 166 Lecture 18 Topics covered Basis vectors Matrix based transforms Wavelet transform Image compression Image watermarking Morphological image processing Segmentation
More informationCPSC 425: Computer Vision
CPSC 425: Computer Vision Image Credit: https://docs.adaptive-vision.com/4.7/studio/machine_vision_guide/templatematching.html Lecture 9: Template Matching (cont.) and Scaled Representations ( unless otherwise
More informationDrawing a Triangle (and an Intro to Sampling)
Lecture 4: Drawing a Triangle (and an Intro to Sampling) Computer Graphics CMU 15-462/15-662, Spring 2018 HW 0.5 Due, HW 1 Out Today! GOAL: Implement a basic rasterizer - (Topic of today s lecture) - We
More informationImage Sampling and Quantisation
Image Sampling and Quantisation Introduction to Signal and Image Processing Prof. Dr. Philippe Cattin MIAC, University of Basel 1 of 46 22.02.2016 09:17 Contents Contents 1 Motivation 2 Sampling Introduction
More informationAliasing And Anti-Aliasing Sampling and Reconstruction
Aliasing And Anti-Aliasing Sampling and Reconstruction An Introduction Computer Overview Intro - Aliasing Problem definition, Examples Ad-hoc Solutions Sampling theory Fourier transform Convolution Reconstruction
More informationRadon Transform and Filtered Backprojection
Radon Transform and Filtered Backprojection Jørgen Arendt Jensen October 13, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging Department
More informationCSE 417T: Introduction to Machine Learning. Lecture 6: Bias-Variance Trade-off. Henry Chai 09/13/18
CSE 417T: Introduction to Machine Learning Lecture 6: Bias-Variance Trade-off Henry Chai 09/13/18 Let! ", $ = the maximum number of dichotomies on " points s.t. no subset of $ points is shattered Recall
More informationRevisit of the Ramp Filter
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 62, NO. 1, FEBRUARY 2015 131 Revisit of the Ramp Filter Gengsheng L. Zeng, Fellow, IEEE Abstract An important part of the filtered backprojection (FBP) algorithm
More informationFrom multiple images to catalogs
Lecture 14 From multiple images to catalogs Image reconstruction Optimal co-addition Sampling-reconstruction-resampling Resolving faint galaxies Automated object detection Photometric catalogs Deep CCD
More informationGRAPHICAL USER INTERFACE (GUI) TO STUDY DIFFERENT RECONSTRUCTION ALGORITHMS IN COMPUTED TOMOGRAPHY
GRAPHICAL USER INTERFACE (GUI) TO STUDY DIFFERENT RECONSTRUCTION ALGORITHMS IN COMPUTED TOMOGRAPHY A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Engineering
More informationOptimised corrections for finite-difference modelling in two dimensions
Optimized corrections for 2D FD modelling Optimised corrections for finite-difference modelling in two dimensions Peter M. Manning and Gary F. Margrave ABSTRACT Finite-difference two-dimensional correction
More informationLectures 12 and 13 Dynamic programming: weighted interval scheduling
Lectures 12 and 13 Dynamic programming: weighted interval scheduling COMP 523: Advanced Algorithmic Techniques Lecturer: Dariusz Kowalski Lectures 12-13: Dynamic Programming 1 Overview Last week: Graph
More informationDEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE
Rajesh et al. : Proceedings of the National Seminar & Exhibition on Non-Destructive Evaluation DEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE Rajesh V Acharya, Umesh Kumar, Gursharan
More informationLessons in Spatial Sampling Or... Does Anybody Know the Shape of a Wavefield?
Lessons in Spatial Sampling Or... Does Anybody Know the Shape of a Wavefield? Norm Cooper - Mustagh Resources Ltd., Calgary, Canada In a homogeneous, isotropic, infinite half space (remember those regimes
More informationDigital 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 informationImaging and Deconvolution
Imaging and Deconvolution Urvashi Rau National Radio Astronomy Observatory, Socorro, NM, USA The van-cittert Zernike theorem Ei E V ij u, v = I l, m e sky j 2 i ul vm dldm 2D Fourier transform : Image
More informationHow and what do we see? Segmentation and Grouping. Fundamental Problems. Polyhedral objects. Reducing the combinatorics of pose estimation
Segmentation and Grouping Fundamental Problems ' Focus of attention, or grouping ' What subsets of piels do we consider as possible objects? ' All connected subsets? ' Representation ' How do we model
More information18th World Conference on Nondestructive Testing, April 2012, Durban, South Africa
18th World Conference on Nondestructive Testing, 16-20 April 2012, Durban, South Africa TOWARDS ESTABLISHMENT OF STANDARDIZED PRACTICE FOR ASSESSMENT OF SPATIAL RESOLUTION AND CONTRAST OF INTERNATIONAL
More informationBME I5000: Biomedical Imaging
1 Lucas Parra, CCNY BME I5000: Biomedical Imaging Lecture 4 Computed Tomography Lucas C. Parra, parra@ccny.cuny.edu some slides inspired by lecture notes of Andreas H. Hilscher at Columbia University.
More informationG Practical Magnetic Resonance Imaging II Sackler Institute of Biomedical Sciences New York University School of Medicine. Compressed Sensing
G16.4428 Practical Magnetic Resonance Imaging II Sackler Institute of Biomedical Sciences New York University School of Medicine Compressed Sensing Ricardo Otazo, PhD ricardo.otazo@nyumc.org Compressed
More informationEE538 - Final Report Design of Antenna Arrays using Windows
EE538 - Final Report Design of Antenna Arrays using Windows Mahadevan Srinivasan April 29, 2010 Abstract The Design of Uniformly-spaced Antenna Arrays has a significant similarity to the Design of FIR
More informationAdvanced Imaging Trajectories
Advanced Imaging Trajectories Cartesian EPI Spiral Radial Projection 1 Radial and Projection Imaging Sample spokes Radial out : from k=0 to kmax Projection: from -kmax to kmax Trajectory design considerations
More information