Optimum Notch Filtering - I

Similar documents
Image and Multidimensional Signal Processing

ECE472/572 - Lecture 11. Roadmap. Roadmap. Image Compression Fundamentals and Lossless Compression Techniques 11/03/11.

Image Compression. Fundamentals: Coding redundancy. The gray level histogram of an image can reveal a great deal of information about the image

Midterm Summary Fall 08. Yao Wang Polytechnic University, Brooklyn, NY 11201

Image Compression. Qiaoyong Zhong. November 19, CAS-MPG Partner Institute for Computational Biology (PICB)

Reduce the amount of data required to represent a given quantity of information Data vs information R = 1 1 C

Color Science Light & Spectra

Compression. What. Why. Reduce the amount of information (bits) needed to represent image Video: 720 x 480 res, 30 fps, color

Department of Electrical Engineering, Polytechnic University, Brooklyn Fall 05 EL DIGITAL IMAGE PROCESSING (I) Final Exam 1/5/06, 1PM-4PM

Number Representation and Waveform Quantization

IMAGE COMPRESSION-II. Week IX. 03/6/2003 Image Compression-II 1

Digital Image Processing Lectures 25 & 26

Lecture 2: Introduction to Audio, Video & Image Coding Techniques (I) -- Fundaments

Fresnel Equations cont.

Computer Vision & Digital Image Processing

Lecture 2: Introduction to Audio, Video & Image Coding Techniques (I) -- Fundaments. Tutorial 1. Acknowledgement and References for lectures 1 to 5

Multimedia Information Systems

Digital communication system. Shannon s separation principle

Digital Image Processing

Intensity Transformations and Spatial Filtering: WHICH ONE LOOKS BETTER? Intensity Transformations and Spatial Filtering: WHICH ONE LOOKS BETTER?

at Some sort of quantization is necessary to represent continuous signals in digital form

EE5356 Digital Image Processing

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Basic Principles of Video Coding

IMAGE COMPRESSION IMAGE COMPRESSION-II. Coding Redundancy (contd.) Data Redundancy. Predictive coding. General Model

Run-length & Entropy Coding. Redundancy Removal. Sampling. Quantization. Perform inverse operations at the receiver EEE

COMPRESSIVE (CS) [1] is an emerging framework,

L. Yaroslavsky. Fundamentals of Digital Image Processing. Course

Module 2 LOSSLESS IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

6.003: Signals and Systems. Sampling and Quantization

EECS490: Digital Image Processing. Lecture #11

Transform Coding. Transform Coding Principle

Multimedia communications

Use the slope-intercept form to graph the equation. 8) 6x + y = 0

Screen-space processing Further Graphics

Multimedia Communications. Scalar Quantization

SYDE 575: Introduction to Image Processing. Image Compression Part 2: Variable-rate compression

Vector Quantization and Subband Coding

Image Acquisition and Sampling Theory

Wavelet Scalable Video Codec Part 1: image compression by JPEG2000

EE5585 Data Compression April 18, Lecture 23

Visual Imaging and the Electronic Age Color Science

EE376A: Homework #2 Solutions Due by 11:59pm Thursday, February 1st, 2018

Scalable color image coding with Matching Pursuit

Proyecto final de carrera

Overview. Analog capturing device (camera, microphone) PCM encoded or raw signal ( wav, bmp, ) A/D CONVERTER. Compressed bit stream (mp3, jpg, )

Visual Imaging and the Electronic Age Color Science

LORD: LOw-complexity, Rate-controlled, Distributed video coding system

Image Enhancement: Methods. Digital Image Processing. No Explicit definition. Spatial Domain: Frequency Domain:

Histogram Processing

EE5356 Digital Image Processing. Final Exam. 5/11/06 Thursday 1 1 :00 AM-1 :00 PM

Source Coding for Compression

NEW PIXEL SORTING METHOD FOR PALETTE BASED STEGANOGRAPHY AND COLOR MODEL SELECTION

<Outline> JPEG 2000 Standard - Overview. Modes of current JPEG. JPEG Part I. JPEG 2000 Standard

Image Data Compression

Colour Part One. Energy Density CPSC 553 P Wavelength 700 nm

Pulse-Code Modulation (PCM) :

ECE Information theory Final (Fall 2008)

encoding without prediction) (Server) Quantization: Initial Data 0, 1, 2, Quantized Data 0, 1, 2, 3, 4, 8, 16, 32, 64, 128, 256

Wavelets and Multiresolution Processing

Statistical Analysis and Distortion Modeling of MPEG-4 FGS

Compression methods: the 1 st generation

JPEG Standard Uniform Quantization Error Modeling with Applications to Sequential and Progressive Operation Modes

Module 4. Multi-Resolution Analysis. Version 2 ECE IIT, Kharagpur

BASICS OF COMPRESSION THEORY

Multimedia Networking ECE 599

Basics on 2-D 2 D Random Signal

Multidimensional Signal Processing

Huffman Coding. C.M. Liu Perceptual Lab, College of Computer Science National Chiao-Tung University

Version 087 EX4 ditmire (58335) 1

Objective: Reduction of data redundancy. Coding redundancy Interpixel redundancy Psychovisual redundancy Fall LIST 2

EE67I Multimedia Communication Systems

Subband Coding and Wavelets. National Chiao Tung University Chun-Jen Tsai 12/04/2014

18/10/2017. Image Enhancement in the Spatial Domain: Gray-level transforms. Image Enhancement in the Spatial Domain: Gray-level transforms

SIGNAL COMPRESSION. 8. Lossy image compression: Principle of embedding

Communication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi

C.M. Liu Perceptual Signal Processing Lab College of Computer Science National Chiao-Tung University

Image Compression using DPCM with LMS Algorithm

Digital Trimulus Color Image Enhancing and Quantitative Information Measuring

Shannon's Theory of Communication

IASI PC compression Searching for signal in the residuals. Thomas August, Nigel Atkinson, Fiona Smith

Digital Image Processing

Jim Hagerman 4/12/99

Lossless Compression of Dynamic PET Data

Analysis of Rate-distortion Functions and Congestion Control in Scalable Internet Video Streaming

CSE 408 Multimedia Information System Yezhou Yang

Image Degradation Model (Linear/Additive)

Lec 04 Variable Length Coding (VLC) in JPEG

EC 624 Digital Image Processing ( ) Class I Introduction. Instructor: PK. Bora

Fully smoothed L1-TV energies : properties of the optimal solutions, fast minimization and applications. Mila Nikolova

Visual Imaging and the Electronic Age Color Science Metamers & Chromaticity Diagrams. Lecture #6 September 7, 2017 Prof. Donald P.

The information loss in quantization

CS Color. Aditi Majumder, CS 112 Slide 1

Introduction p. 1 Compression Techniques p. 3 Lossless Compression p. 4 Lossy Compression p. 5 Measures of Performance p. 5 Modeling and Coding p.

CS681 Computational Colorimetry

홀로그램저장재료. National Creative Research Center for Active Plasmonics Applications Systems

Multimedia Systems Giorgio Leonardi A.A Lecture 4 -> 6 : Quantization

Multimedia & Computer Visualization. Exercise #5. JPEG compression

Navneet Agrawal Deptt. Of Electronics & Communication Engineering CTAE,MPUAT,Udaipur,India

4 An Introduction to Channel Coding and Decoding over BSC

Transcription:

Notch Filtering 1

Optimum Notch Filtering - I Methods discussed before: filter too much image information. Optimum: minimizes local variances of the restored image fˆ x First Step: Extract principal freq. components of the interference pattern. Place a notch pass filter at the location of each spike. N u v H u v G u v NP Second Step: Find corresponding pattern in spatial domain. n x 1 H NP u v G u v Eq.1 2

Optimum Notch Filtering - II Third Step: Subtract weighted noise estimation from nois image. fˆ x g x w x n x Eq.2 How to get wx? Select wx so that the variance of over a neighborhood. fˆ x is minimized Consider: neighborhood size: 2a+1 b 2b+1 Average: a b 1 f ˆ x fˆ x s t 2a 12b 1 satb Local variance: 2 x 1 2a 12b 1 a b satb fˆ x s t fˆ x 2 3

4 Optimum Notch Filtering - III a a s b b t x n w x x g t s x n w x t s x g b a x 2 2 1 12 2 1 0 2 w x x 2 2 x n x n x n x g x n x g x w To minimize We find Eq.3 Get noise reduced image from Eq. 1 2 3.

Image Reconstruction : projection Computed Tomograph CT X-Ras from different angles. Back projection Angle varied: 90 Sum of two back projections 5

Projection Projection using man angles: more true construction of the original image. Sum of 32 back projections. 6

Principles of CT G1: Pencil X-Ra beam; one detector; angle [0 ~ 180]. G2: Fan X-Ra beam; multiple detectors. G3: Wider X-Ra beam; a bank of detectors 1000. G4: Circular positioned detectors 5000. G5: No mechanical motion uses electron beams controlled electromagneticall. G6 G7.. 7

Color Fundamentals White color: composed of 6 visible colors. 8

Primar & Secondar Colors Green Red Blue Yellow Magenta Can 9

Characteristics of Colors 1. Brightness: how bright intensit the color is. [Perception.] 2. Hue: dominant wavelength in a mixture of light waves. 3. Saturation: the amount of white light mixed with a hue. Inversel proportional Hue + Saturation = Chromaticit 10

RGB Color Model Black 000 White 111 24-bit color cube. All R G B values are normalized to [0 1] range. 11

Safe Color Man sstems in use toda are limited to 256 colors. Fort 40 of these are processed differentl b various operating sstems. The rest 216 are common: de facto standard for safe colors. 6 3 = 216 12

13 CMY and CMYK Color Models B G R Y M C 1 1 1 Used as primar colors in printers. Equal amounts of C M Y should produce Black. But in practice it produces mudd-looking black. Four-color printing: add a fourth color: Black to CMY producing CMYK model.

HSI Color Model - I Hue Saturation Intensit. 14

HSI Color Model - II Triangular Plane Circular Plane 15

Fundamentals of Image Compression Relative data redundanc R 1 1 C Compression ratio: C b b Bits required in the method Bits required in the base method Average number of bits required to represent each pixel: L 1 LAvg l rk Pr rk k 0 nk Pr rk k 01 2... L 1 MN 16

Variable Length Coding Table 1. For code 1: l 1 r k = 8 bits for al r k. Average = 8. For code 2: L Avg 0.252 0.471 0.253 0.033 1.81 bits 256 2568 1 C 4.42 R 1 0. 774 256 2561.81 4.42 77.4% of the data in the original 8-bit intensit arra is redundant. 17

Irrelevant Information Spatial redundanc Figure 8.1 Histogram Histogram equalized image 18

Entrop Entrop: H L 1 k0 p r r k log 2 pr rk Entrop in Table 1.: 1.6614 bits / pixel Entrop of Figure 1 b = 8 bits / pixel Entrop of Figure 1 c = 1.566 bits / pixel Figure 1 a Higher entrop than Fig. 1 a Little or no information But comparable entrop with Fig. 1a! Entrop of an image is far from intuitive. 19

Fidelit Criteria - I Quantifing the nature of the loss after removal of noise. Objective: mathematical expression such as rms error SNR etc. e rms 1 MN M 1N 1 x0 0 fˆ x f x 2 1/ 2 Subjective: Human evaluation. 20

Fidelit Criteria - II Misleading image rms error: 5.17 15.67 14.17 Objective criteria fails. 21

Image Compression Models Quantizer: irreversible. 22

Some Compression Standards 23

Huffman Coding Encoded: 0101001111100 Decoded: a a a a 3 1 2 2a6 Average length of this code: 0.41 0.32 0.13 0.065 0.045 L Avg 2.2 bits/pixel 24

Golomb Coding G m n Step 1: Step 2: Step 3: Form the unar code of quotient n/m. Let k = log 2 m c = 2 k m r = n mod m compute r r truncatedto k 1 bits r r c truncat edtok bits 0 r c otherwis e Concatenated the results of steps 1 and 2. Example: Compute G 4 9. 9/4 = 2.25 = 2 Unar code: 110 k = 2 c = 2 2 4 = 0 r = 1 0001. After concatenating: G 4 9 = 11001 r = 01 truncated to 2 bits. 25

Exponential Golomb Coding G k exp n Step 1: Find an integer i >= 0 such that And form the unar code of i. i1 j0 2 jk n i j0 2 jk If k = 0 i = log 2 n+1 Step 2: Step 3: Truncate the binar representation of n i 1 j0 2 jk for k i Concatenated the results of steps 1 and 2. least significan t bits. Example: Compute G 0 exp8. i = 3 because k = 0. 1110 Check for the equation in Step 1. 8 3 1 j0 2 j0 Truncate 8710001 After concatenating: G 0 exp8 = 1110001 001 26