Mathematics of Image Processing  (MA439-01)  -  S. Allen Broughton, Fall 2002

Course Materials Page

[Class Handouts] [Computer Resources] [Homework List] [Homework Policies] [ Course home page]

Contact Instructor: allen.broughton@rose-hulman.edu -  last update 2 Mar 02


MA439 Class Handouts

Notes

No longer here

Worksheets & Projects

Quizzes and Exams


MA439 Computer Resources

Updating Maple and Installing Matlab

Maple and Matlab Scripts

 
Script Name  Demo/Template/Utility/Review Purpose
MapleIntrRGL.mws review for MA221 Basic Maple Introduction 
calcrev.mws review for MA221 Review of Maple commands from Calculus
matrices.mws demo (review from MA221) matrix operations
matlabeg.m intro & review simple intro to Matlab
cosinesinesums.mws demo and template analysis and synthesis of discrete and continuous signals
datasamp.mws data file for worksheet #1
cossinmat.m
cossinmat2.m
cossinmat3.m
utilities & templates shows how to calculate the trigonometric wave form matrix
Ncossinmat.m
Ncossinmat2.m
Ncossinmat3.m
Ncossinmat4.m
utilities & templates

shows how to calculate the trigonometric wave form matrix
matrices are structurally different than the the examples above

 

transformpairs.m demo shows signal and transform pairs for cosine sine transform
zerocount.m
zerocountdemo.m
utility
demo
counts the number of interpolated zeros between samples of a function, zerocount = 2x frequency, demo shows graphs of functions
wave1d.m demo show real and imaginary parts of exponential wave forms
show continuous and discrete samples
wave1dsteps.m
wave1dmovie.m
demos - show real and imaginary parts of exponential wave forms one wave at a time, show "frequencies"
- show real and imaginary parts of exponential wave forms as a "movie" of varying frequency"
wave1dfreqdemo.m
wave1dsymfreqdemo.m
demos

- graphically show how wave forms depend on the frequencies in the range 0..N-1
- graphically show how wave forms depend on the symmetric frequencies in the range -N/2 .. N/2

wave2d.m
wave2dcoldemo.m
demos

- show a 2d exponential wave form as an image and a surface
- as above but vary the color scheme

RGBOpic.m demo show the colour components of an image.
colormapdisplay demo show the R,G,B and average curves of a colormap
wave2dtrigcc.m demo show a factorized 2D trigonometric wave
wave2d.m
wave2dcoldemo.m
demos - show the real and imaginary parts of a 2D exponential wave
- show the real and imaginary parts of a 2D exponential wave, play with colormap
wavefft2lband.m
wavefft2cband.m
demos
- show a "linear" band of waveforms
- show a "circular" band of waveforms
dotanalysis.m demo - show how Fourier coefficients can be calculated by dot products
dft1demo1.m
dft1demo2.m
dft1demo3.m
demos - show signal and DFT pairs in Cartesian and polar form
- show signal and shifted DFT pairs in Cartesian and polar form
- movie of a DFT computation
dft2demo1.m
dft2demo2.m
dft2demo3.m
demos - show image, DFT and log(1+|DFT|) triples
- show image, log(1+|DFT|) and shifted log(1+|DFT|) triples
- show before and after pictures of an image and DFT when a mask is applied to a DFT
acomp.m
icomp.m
demos - compress an audio signal using global DFT
- compress an image using global DFT
energy.m utility compute the energy of a signal or image
jpeg2grayscale.m utility convert a .jpg file and convert it to grayscale image in uint8 format
symdft1demo.m demo show symmetric extension of a function, its DFT and derived (pseudo)-DCT
dctdemo.m demo properties of small DCT matrices
analsynDCT.m demo movie of DCT analysis and partial reconstruction of a signal
blkdct2.m
blkidct2.m
utilities construct mxn block DCT's and their inverses
circconv.m
circconv2.m
circconvs.m
circconv2s.m
utilities

circular convolution scripts
the two scripts circconvs.m, circconv2s.m allow for circular
convolution that is centred around 0 or (0,0).

conv1demo1.m
conv1demo2.m
demos - demonstrate before and after pictures of a filtered signal and the DFT of the mask
- movie of filtering
convtheorem.m demo demonstrate filtering exponential waveforms
alohi.m demo show and hear averaging and differencing of an audio signal
edgedet.m demo show edge detection by filtering and non-linear superposition
lohi1d.m demo show low and high pass filtering on a 1D signal
lohi2d.m
lohi2d2.m
demos show lo-lo, lo-hi, hi-lo, and hi-hi pass filtering on a natural image
show lo-lo, lo-hi, hi-lo, and hi-hi pass filtering on a smooth image defined by a function graph
dwt1demo.m demo show 1-level and 2-level DWT on a 1D signal

downc.m
downr.m

 

utilities

downsample columns
downsample rows
similar to dyadown in wavelet package

hiloanlc.m
hiloanlr.m
hiloanl2.m
hilosync.m
hilosynr.m
hilosyn2.m
utilities low pass high pass analysis of the columns of a matrix
low pass high pass analysis of the rows of a matrix
low pass high pass analysis of the columns and rows of a matrix
low pass high pass synthesis of the columns of a matrix
low pass high pass synthesis of the rows of a matrix
low pass high pass synthesis of the columns and rows of a matrix
uses utilities in previous row, plus circular convolution utilities
similar to periodic DWT, DWT2, IDWT and IDWT2

fbanlcol.m
fbanlrow.m
fbanl2.m
fbsyncol.m
fbsynrow.m
fbsyn2.m

utilities s-level periodic filter bank analysis of the columns of a matrix (s-level periodic DWT)
s-level periodic filter bank analysis of the rows of a matrix (s-level periodic DWT)
s-level periodic filter bank analysis of a matrix (s-level periodic 2D DWT)
s-level periodic filter bank synthesis of the columns of a matrix (s-level periodic DWT)
s-level periodic filter bank synthesis of the rows of a matrix (s-level periodic DWT)
s-level periodic filter bank synthesis of a matrix (s-level periodic 2D DWT)
uses utilities in previous row
similar to periodic wavedec and waverec
analsynDWT.m demo show anlyais and synthesis waveforms for DWT
dwt2demo.m
dwt2demo2.m
demos two slightly different wavelet transform of an image and various types of details
uses wavelet toolbox
flatness.m demo demonstrate relation between flatness of DFT of filter and killing polynomial signals
tap4show1d.m
tap4show2d.m
demo show the effect of flatness of the high pass filter on the 1D DWT for the family of 4 tap orthogonal filters
show the effect of flatness of the high pass filter on the 2D DWT for the family of 4 tap orthogonal filters
findfilters.mws demo Show how to find the synthesis filters for the JPEG2000 (5,3) filter.

Matlab Tutorials


MA439 Homework, Solutions, and Schedule

Policies

Homework List & Schedule (in left column, stuff in red is to be handed in, schedule items are in blue)

Assignment or Schedule item Date or Completion Date Source Solutions
Assignment 1 Monday, 9 Sep 02 Class Handouts  
Worksheet 1 Thursday, 12 Sep 02 worksheet#1, datasamp.mws  
Worksheet 1, problem 7 Friday, 13 Sep 02 worksheet#1, datasamp.mws problem7.mws

Matlab session 1 Bring computer

Friday, 13 Sep 02 class handout
mtlbsess02-1.pdf
 
Reading Monday, 16 Sep 02 notes Chapters 1 and 2 except the 2D stuff  
Matlab Assignment 1 Tuesday, 17 Sep 02 mtlbsess02-1.pdf mtlbsess02-1-sol.txt
Assignment 2 Friday, 20 Sep 02 Notes, p. 40+, #2.1, 2.4, 2.5, 2.6  
Assignment 3 Monday, 23 Sep 02 Notes, p. 41+, #2.10, 2.11, 2.15, 2.18  
Quiz #1 Tuesday, 24 Sep 02 Exponential wave forms, Hermitian inner products. ansquiz1.pdf
Matlab session 2 Bring computer Thursday, 26 Sep 02 mtlbsess02-2.pdf  
Read Notes Thursday, 26 Sep 02 Read notes on the 2D topics.  
Matlab Assignment 2 Monday, 30 Sep 02 mtlbsess02-2.pdf mtlbsess02-2-sol.txt
Read Chapter 3 Thursday, 3 Oct 02    
Assignment 4 Thursday, 10 Oct 02 3.3, 3.4, 3.10 ansassn02-04.pdf
Assignment 5 Tuesday, 15 Oct 02 assn02-05.pdf  
Assignment 6 Tuesday, 22 Oct 02 4.3, 4.4, 4.6, What is the frequency of C_N,k  
Midterm Project Friday, 25 Oct 02 midproj02.pdf  
Worksheet 4 Thursday, 31 Oct 02 worksheet#4  
Worksheet 5 Friday, 1 Nov 02 worksheet#5 (look at Problem 5.6 for constructing convolution matrices) answs04-02.pdf
Worksheet 6 Thursday, 7 Nov 02 worksheet#6 worksheet6.m
Bring Computer Friday, 8 Nov 02 worksheet#7 worksheet7.m
Bring Computer Monday, 11 Nov 02 worksheet#8 worksheet8.m
Bring Computer Tuesday, 12 Nov 02    
Final Exam Monday, 18 Nov 02 takehomefinal.pdf