Hi all,

I am having a hard time understanding the meaning, of the PSF and the
convolution term in optic .
I would be very glad if anyone can help.

- what is the meaning of PSF in optical system
- does the PSF is a function that I can create s only when I try to
create an image with a lens, and a point source?
(is there a PSF for other source?)

and how you calculate it?

what is the meaning of convolution in Optic ??

"Marc Reinig" wrote

The PSF (Point Spread Function) of an optical system is what the system
does to a point source object in the process of imaging it.  For a round
pupil with no aberrations in the system it will be a sync squared function
( ( sin(x)^2 ) / ( x^2 ) whose extent is a function of the diameter of the
pupil and the focal length of the system.  This function has a central
bright core and a series of dark rings and decreasing intensity bright
rings around it.

Oops...

Pls have a look athttp://en.wikipedia.org/wiki/Airy_disk

(among many such references, links etc.)

no sine (or sinc) here but the Bessel function J1   (J sub 1) etc.

The above applies to a square or rectangular pupil (also with the
appropriate variables x and y).

As told before, see Hecht... although here it's more likely to be a
"slip of the keyboard" and not a mistake given the author's expertise
in the field.    

The PSF (Point Spread Function) of an optical system is what the system
does to a point source object in the process of imaging it. For a round
pupil with no aberrations in the system it will be a sync squared function
( ( sin(x)^2 ) / ( x^2 ) whose extent is a function of the diameter of the
pupil and the focal length of the system. This function has a central
bright core and a series of dark rings and decreasing intensity bright
rings around it.

odprt...@gmail.com> wrote in message

news:06c55b30-32f7-4b13-b928-02a63ccda6d1@w4g2000prd.googlegroups.com...

Hi all,

I am having a hard time understanding the meaning, of the PSF and the
convolution term in optic  .
I would be very glad if anyone can help.

- what is the meaning of PSF in optical system
- does the PSF is a function that I can create  s only when I try to
create an image with  a lens, and a point source?
(is there a PSF for other source?)

The PSF (Point Spread Function) of an optical system is what the system does
to a point source object in the process of imaging it.  For a round pupil
with no aberrations in the system it will be a sync squared function  ( (
sin(x)^2 ) / ( x^2 ) whose extent is a function of the diameter of the pupil
and the focal length of the system.  This function has a central bright core
and a series of dark rings and decreasing intensity bright rings around it..
Any aberrations in the system will make the image of a point source
different and you will have a less than ideal PSF, i.e., the image of a
point in a non ideal system will be spread out compared to an ideal system..

If you are trying to image two point sources close together, each will
produce its own PSF in the image plane.  As the points sources get closer
and closer the PSFs will get closer and closer.  At some point, you will not
be able to tell if you have two points or an extended source that is larger
than a point.  Generally, when   the peak of on PSF is at the same position
as the first null of another PSF, that is considered the point beyond which
you cannot tell you have two points.

Generally, four your purpose, since it is a "point" spread function, it is
only used to describe the effect on point sources.

and how you calculate it?

See below for calculations.  However, for complex instruments, you would
usually measure it rather than calculate it.

what is the meaning of convolution in Optic ??

Light passing through a focusing lens also is passing through a pupil,  area
at the lens through which light is allowed to pass.  Think of the pupil as a
function that is 1 inside and 0 everywhere outside.  This assumes the pupil
has no aberrations in either intensity or phase, i.e., light passes through
uniformly in intensity inside the aperture and is slowed down uniformly
inside the aperture.

Now think of the light coming to the aperture.  The light on one side is
L0(x,y).  The aperture function is A(x,y).  After passing through the
aperture the light is L1 (x,y) = L0(x,y) A(x,y).  i.e., at each point x,y
the light after the aperture will either be L0(x,y) because it is multiplied
by 1; or 0 because it was multiplied by 0.

The light at the image plane turns out to be the Fourier transform of the
light at the pupil plane.  An important relationship in math is that the
Fourier transform of a product is the Fourier transform of the convolution
of the elements and vice versa.  So,

    F( X Y ) = F( X * Y )  and F( X * Y ) = F( X Y )  where * indicates the
convolution.

So, the amplitude of the light at the image plane is the Fourier transform
of the light at the pupil plane:

    Image(x,y) = F( L1(x,y) ) = F( L0(x,y) A(x,y) ) = F( L0(x,y) ) * F(
A(x,y) )

It is the convolution of the Fourier transform of the light at the pupil
plane with the Fourier transform of the Pupil function.

If the source light is a point source, the PSF would be the intensity in the
image plane.  Intensity if the square of the light amplitude, we perceive
intensity, not whether the E field of the light is + or -.

For an example of how this can be used, let's imagine a camera with a non
perfect lens and/or pupil aperture.  The image taken will be less than
ideal, but it will be the convolution of the Fourier transform of the bad
pupil/lens with the Fourier transform of the light at the pupil.  If you
could take the camera in isolation you could characterize it with a point
source to measure the bad PSF and develop a mathematical expression for it..

You know the image is the convolution of the bad PSF with the Fourier
transform of the light at the pupil.  You could then take the image and
deconvolve it (using the expression of the bad PSF) to get an expression for
the light at the pupil when the image was taken.  Then you could take that
expression and process it with the expression for a perfect PSF and get the
image as it would have been recorded with a perfect camera.

The theory is better than the practice and it is very computationally
intense.  However, it has been used to sharpen the original Hubble images
taken with the original improperly designed lens, it is used to enhance
images by law enforcement to correct for all manner of defects in
surveillance video, etc.  It is however a slow process, does not produce
enhancement back to having a perfect lens in the original camera, and not
the way you see it in CSI, NCIS, .... ;=)

Marco
________________________
Marc Reinig
UCO/Lick Observatory
Laboratory for Adaptive Optics

"Marc Reinig" wrote

The PSF (Point Spread Function) of an optical system is what the system
does to a point source object in the process of imaging it. For a round
pupil with no aberrations in the system it will be a sync squared
function ( ( sin(x)^2 ) / ( x^2 ) whose extent is a function of the
diameter of the pupil and the focal length of the system. This function
has a central bright core and a series of dark rings and decreasing
intensity bright rings around it.

Oops...

Pls have a look at http://en.wikipedia.org/wiki/Airy_disk

(among many such references, links etc.)

no sine (or sinc) here but the Bessel function J1 (J sub 1) etc.

The above applies to a square or rectangular pupil (also with the
appropriate variables x and y).

As told before, see Hecht... although here it's more likely to be a
"slip of the keyboard" and not a mistake given the author's expertise
in the field. :-)

Dear No Name,

You should check out Hecht's book Optics, and refer to his discussion
of the point spread function. Its a bit heavier on description and
lighter on equations. A good resource to understand the equations
better is Jack Gaskill's book Linear Systems, Fourier Transforms, and
Optics.

The PSF is the resulting image of a point through an optical system.
The PSF can be theoretically calculated with parameters for a
theoretical optical system, and can be demonstrated with a real (hold
it in your hands) optical system. PSFs are often measured in
laboratories of optical systems for test.

How do you calculate it? Read Gaskill's book or Goodman... I prefer
Gaskill.

In fourier space you can convolve the optical transfer function of one
optical system or lens with another to derive a complete description
of some real or theoretical optical system to compute the PSF of.

I hope the book title recommendations will steer you to some better
explanations than what you've been reading.

-Philip

    F( X Y ) = F( X * Y )  and F( X * Y ) = F( X Y )  where * indicates the
convolution.

So, the amplitude of the light at the image plane is the Fourier transform
of the light at the pupil plane:

    Image(x,y) = F( L1(x,y) ) = F( L0(x,y) A(x,y) ) = F( L0(x,y) ) * F(
A(x,y) )

    F( X Y ) = F( X * Y )  and F( X * Y ) = F( X Y )  where * indicates the
convolution.

So, the amplitude of the light at the image plane is the Fourier transform
of the light at the pupil plane:

    Image(x,y) = F( L1(x,y) ) = F( L0(x,y) A(x,y) ) = F( L0(x,y) ) * F(
A(x,y) )

 something is not right in this equations, I think there is a mistake,
tell me if I am wrong here...
it should be :

  F( X)F( Y ) = F( X * Y )

what do u say?

also after reading few things, and reading this comment again.

in the end the PSF is nothing but a simple transform Fourier over the
aperture , and the light that hit it, no ?

also tell me if i understand what u say correctly
if you have the PSF of a system, then you can calculate how any signal
will look, by convulting this signal with the transform the PSF .
is that right? :)

thanks

On Apr 16, 9:39 pm, "odprt...@gmail.com" <odprt...@gmail.com> wrote:



    F( X Y ) = F( X * Y )  and F( X * Y ) = F( X Y )  where * indicates the
convolution.

So, the amplitude of the light at the image plane is the Fourier transform
of the light at the pupil plane:

    Image(x,y) = F( L1(x,y) ) = F( L0(x,y) A(x,y) ) = F( L0(x,y) ) * F(
A(x,y) )

 something is not right in this equations, I think there is a mistake,
tell me if I am wrong here...
it should be :

  F( X)F( Y ) = F( X * Y )

what do u say?

also after reading few things, and reading this comment again.

in the end the PSF is nothing but a simple transform Fourier over the
aperture , and the light that hit it, no ?

also tell me if i understand what u say correctly
if you have the PSF of a system, then you can calculate how any signal
will look, by convulting this signal with the transform the PSF .
is that right? :)

thanks

meant:
if you have the PSF of a system, then you can calculate how any signal
 will look, by convulting this signal with  the PSF .
A

On Apr 16, 9:39 pm, "odprt...@gmail.com" <odprt...@gmail.com> wrote:



    F( X Y ) = F( X * Y )  and F( X * Y ) = F( X Y )  where * indicates the
convolution.

So, the amplitude of the light at the image plane is the Fourier transform
of the light at the pupil plane:

    Image(x,y) = F( L1(x,y) ) = F( L0(x,y) A(x,y) ) = F( L0(x,y) ) * F(
A(x,y) )

 something is not right in this equations, I think there is a mistake,
tell me if I am wrong here...
it should be :

  F( X)F( Y ) = F( X * Y )

what do u say?

also after reading few things, and reading this comment again.

in the end the PSF is nothing but a simple transform Fourier over the
aperture , and the light that hit it, no ?

also tell me if i understand what u say correctly
if you have the PSF of a system, then you can calculate how any signal
will look, by convulting this signal with the transform the PSF .
is that right? :)

thanks

meant:
if you have the PSF of a system, then you can calculate how any signal
 will look, by convulting this signal with  the PSF .

On Apr 16, 9:39 pm, "odprt...@gmail.com" <odprt...@gmail.com> wrote:



    F( X Y ) = F( X * Y )  and F( X * Y ) = F( X Y )  where * indicates the
convolution.

So, the amplitude of the light at the image plane is the Fourier transform
of the light at the pupil plane:

    Image(x,y) = F( L1(x,y) ) = F( L0(x,y) A(x,y) ) = F( L0(x,y) ) * F(
A(x,y) )

 something is not right in this equations, I think there is a mistake,
tell me if I am wrong here...
it should be :

  F( X)F( Y ) = F( X * Y )

what do u say?

also after reading few things, and reading this comment again.

in the end the PSF is nothing but a simple transform Fourier over the
aperture , and the light that hit it, no ?

also tell me if i understand what u say correctly
if you have the PSF of a system, then you can calculate how any signal
will look, by convulting this signal with the transform the PSF .
is that right? :)

thanks

meant:
if you have the PSF of a system, then you can calculate how any signal
 will look, by convulting this signal with  the PSF .

odprt...@gmail.com wrote:
On Apr 16, 9:46 pm, "odprt...@gmail.com" <odprt...@gmail.com> wrote:
On Apr 16, 9:39 pm, "odprt...@gmail.com" <odprt...@gmail.com> wrote:

    F( X Y ) = F( X * Y )  and F( X * Y ) = F( X Y )  where * indicates the
convolution.
So, the amplitude of the light at the image plane is the Fourier transform
of the light at the pupil plane:
    Image(x,y) = F( L1(x,y) ) = F( L0(x,y) A(x,y) ) = F( L0(x,y) ) * F(
A(x,y) )
 something is not right in this equations, I think there is a mistake,
tell me if I am wrong here...
it should be :
  F( X)F( Y ) = F( X * Y )
what do u say?
also after reading few things, and reading this comment again.
in the end the PSF is nothing but a simple transform Fourier over the
aperture , and the light that hit it, no ?
also tell me if i understand what u say correctly
if you have the PSF of a system, then you can calculate how any signal
will look, by convulting this signal with the transform the PSF .
is that right?
thanks
meant:
if you have the PSF of a system, then you can calculate how any signal
 will look, by convulting this signal with  the PSF .

anyone ?

Convolution in the space domain is multiplication in the spatial
frequency domain, and vice versa.  So your first set of equations should
read

F(X*Y) = F(X)F(Y)    and F(XY) = F(X)*F(Y).

You have to be very careful to distinguish between the coherent transfer
function (CTF), which operates on field amplitudes, and the less well
defined and less well behaved optical transfer function (OTF), which
operates on intensity.  See lots of earlier threads on this newsgroup.

Cheers,

Phil Hobbs

On Apr 16, 9:46 pm, "odprt...@gmail.com" <odprt...@gmail.com> wrote:
On Apr 16, 9:39 pm, "odprt...@gmail.com" <odprt...@gmail.com> wrote:



F( X Y ) = F( X * Y ) and F( X * Y ) = F( X Y ) where * indicates the
convolution.
So, the amplitude of the light at the image plane is the Fourier transform
of the light at the pupil plane:
Image(x,y) = F( L1(x,y) ) = F( L0(x,y) A(x,y) ) = F( L0(x,y) ) * F(
A(x,y) )
something is not right in this equations, I think there is a mistake,
tell me if I am wrong here...
it should be :
F( X)F( Y ) = F( X * Y )
what do u say?
also after reading few things, and reading this comment again.
in the end the PSF is nothing but a simple transform Fourier over the
aperture , and the light that hit it, no ?
also tell me if i understand what u say correctly
if you have the PSF of a system, then you can calculate how any signal
will look, by convulting this signal with the transform the PSF .
is that right?
thanks
meant:
if you have the PSF of a system, then you can calculate how any signal
will look, by convulting this signal with the PSF .

anyone ?