<< Chapter < Page Chapter >> Page >
This module extends the ideas of the Discrete Fourier Transform (DFT) into two-dimensions, which is necessary for any image processing.

2d dft

To perform image restoration (and many other useful image processing algorithms) in a computer, we need a FourierTransform (FT) that is discrete and two-dimensional.

F k l u 2 k N v 2 l N F u v
for k 0 N 1 and l 0 N 1 .
F u v m n f m n u m v m
F k l m N 1 0 n N 1 0 f m n 2 k m N 2 l n N
where the above equation ( ) has finite support for an N x N image.

Inverse 2d dft

As with our regular fourier transforms, the 2D DFT also has an inverse transform that allows us to reconstruct an imageas a weighted combination of complex sinusoidal basis functions.

f m n 1 N 2 k N 1 0 l N 1 0 F k l 2 k m N 2 l n N

Periodic extensions

Illustrate the periodic extension of images.
Got questions? Get instant answers now!

2d dft and convolution

The regular 2D convolution equation is

g m n k N 1 0 l N 1 0 f k l h m k n l

Below we go through the steps of convolving two two-dimensional arrays. You can think of f as representing an image and h represents a PSF, where h m n 0 for m n 1 and m n 0 . h h 0 0 h 0 1 h 1 0 h 1 1 f f 0 0 f 0 N 1 f N 1 0 f N 1 N 1 Step 1 (Flip h ):

h m n h 1 1 h 1 0 0 h 0 1 h 0 0 0 0 0 0
Step 2 (Convolve):
g 0 0 h 0 0 f 0 0
We use the standard 2D convolution equation ( ) to find the first element of our convolved image. In order to better understand what ishappening, we can think of this visually. The basic idea is to take h m n and place it "on top" of f k l , so that just the bottom-right element, h 0 0 of h m n overlaps with the top-left element, f 0 0 , of f k l . Then, to get the next element of our convolved image, we slide the flipped matrix, h m n , over one element to the right and get the following result: g 0 1 h 0 0 f 0 1 h 0 1 f 0 0 We continue in this fashion to find all of the elements ofour convolved image, g m n . Using the above method we define the general formula to find a particular element of g m n as:
g m n h 0 0 f m n h 0 1 f m n 1 h 1 0 f m 1 n h 1 1 f m 1 n 1

Got questions? Get instant answers now!

Circular convolution

What does H k l F k l produce?

2D Circular Convolution

g ~ m n IDFT H k l F k l circularconvolutionin2D

Got questions? Get instant answers now!

Due to periodic extension by DFT ( ):

Based on the above solution, we will let

g ~ m n IDFT H k l F k l
Using this equation, we can calculate the value for each position on our final image, g ~ m n . For example, due to the periodic extension of the images, when circular convolution is applied we willobserve a wrap-around effect.
g ~ 0 0 h 0 0 f 0 0 h 1 0 f N 1 0 h 0 1 f 0 N 1 h 1 1 f N 1 N 1
Where the last three terms in are a result of the wrap-around effect caused by the presence of the images copies located all around it.

Zero padding

If the support of h is M x M and f is N x N , then we zero pad f and h to M N 1 x M N 1 (see ).

Circular Convolution = Regular Convolution

Computing the 2d dft

F k l m N 1 0 n N 1 0 f m n 2 k m N 2 l n N
where in the above equation, n N 1 0 f m n 2 l n N is simply a 1D DFT over n . This means that we will take the 1D FFT of each row; if wehave N rows, then it will require N N operations per row. We can rewrite this as
F k l m N 1 0 f m l 2 k m N
where now we take the 1D FFT of each column, which means that if we have N columns, then it requires N N operations per column.
Therefore the overall complexity of a 2D FFT is O N 2 N where N 2 equals the number of pixels in the image.

Questions & Answers

Three charges q_{1}=+3\mu C, q_{2}=+6\mu C and q_{3}=+8\mu C are located at (2,0)m (0,0)m and (0,3) coordinates respectively. Find the magnitude and direction acted upon q_{2} by the two other charges.Draw the correct graphical illustration of the problem above showing the direction of all forces.
Kate Reply
To solve this problem, we need to first find the net force acting on charge q_{2}. The magnitude of the force exerted by q_{1} on q_{2} is given by F=\frac{kq_{1}q_{2}}{r^{2}} where k is the Coulomb constant, q_{1} and q_{2} are the charges of the particles, and r is the distance between them.
Muhammed
What is the direction and net electric force on q_{1}= 5µC located at (0,4)r due to charges q_{2}=7mu located at (0,0)m and q_{3}=3\mu C located at (4,0)m?
Kate Reply
what is the change in momentum of a body?
Eunice Reply
what is a capacitor?
Raymond Reply
Capacitor is a separation of opposite charges using an insulator of very small dimension between them. Capacitor is used for allowing an AC (alternating current) to pass while a DC (direct current) is blocked.
Gautam
A motor travelling at 72km/m on sighting a stop sign applying the breaks such that under constant deaccelerate in the meters of 50 metres what is the magnitude of the accelerate
Maria Reply
please solve
Sharon
8m/s²
Aishat
What is Thermodynamics
Muordit
velocity can be 72 km/h in question. 72 km/h=20 m/s, v^2=2.a.x , 20^2=2.a.50, a=4 m/s^2.
Mehmet
A boat travels due east at a speed of 40meter per seconds across a river flowing due south at 30meter per seconds. what is the resultant speed of the boat
Saheed Reply
50 m/s due south east
Someone
which has a higher temperature, 1cup of boiling water or 1teapot of boiling water which can transfer more heat 1cup of boiling water or 1 teapot of boiling water explain your . answer
Ramon Reply
I believe temperature being an intensive property does not change for any amount of boiling water whereas heat being an extensive property changes with amount/size of the system.
Someone
Scratch that
Someone
temperature for any amount of water to boil at ntp is 100⁰C (it is a state function and and intensive property) and it depends both will give same amount of heat because the surface available for heat transfer is greater in case of the kettle as well as the heat stored in it but if you talk.....
Someone
about the amount of heat stored in the system then in that case since the mass of water in the kettle is greater so more energy is required to raise the temperature b/c more molecules of water are present in the kettle
Someone
definitely of physics
Haryormhidey Reply
how many start and codon
Esrael Reply
what is field
Felix Reply
physics, biology and chemistry this is my Field
ALIYU
field is a region of space under the influence of some physical properties
Collete
what is ogarnic chemistry
WISDOM Reply
determine the slope giving that 3y+ 2x-14=0
WISDOM
Another formula for Acceleration
Belty Reply
a=v/t. a=f/m a
IHUMA
innocent
Adah
pratica A on solution of hydro chloric acid,B is a solution containing 0.5000 mole ofsodium chlorid per dm³,put A in the burret and titrate 20.00 or 25.00cm³ portion of B using melting orange as the indicator. record the deside of your burret tabulate the burret reading and calculate the average volume of acid used?
Nassze Reply
how do lnternal energy measures
Esrael
Two bodies attract each other electrically. Do they both have to be charged? Answer the same question if the bodies repel one another.
JALLAH Reply
No. According to Isac Newtons law. this two bodies maybe you and the wall beside you. Attracting depends on the mass och each body and distance between them.
Dlovan
Are you really asking if two bodies have to be charged to be influenced by Coulombs Law?
Robert
like charges repel while unlike charges atttact
Raymond
What is specific heat capacity
Destiny Reply
Specific heat capacity is a measure of the amount of energy required to raise the temperature of a substance by one degree Celsius (or Kelvin). It is measured in Joules per kilogram per degree Celsius (J/kg°C).
AI-Robot
specific heat capacity is the amount of energy needed to raise the temperature of a substance by one degree Celsius or kelvin
ROKEEB
Got questions? Join the online conversation and get instant answers!
Jobilize.com Reply

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Intro to digital signal processing. OpenStax CNX. Jan 22, 2004 Download for free at http://cnx.org/content/col10203/1.4
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Intro to digital signal processing' conversation and receive update notifications?

Ask