Register to reply 
A basic question: jinc function in coherent and incoherent system 
Share this thread: 
#1
Feb913, 11:27 PM

P: 37

I am learning Fourier optics recently and I have a problem of the jinc function.
In optical systems, digital image is blurred with the kernel of jinc function [itex]h(x,y)=jinc(r)=\frac{J_1(2\pi r / \lambda)}{ 2\pi r / \lambda}[/itex] in the coherent system, the blurred image [itex] g(x,y) = h(x,y) \star f(x,y)^2 [/itex] where f(x,y) is the unblurred image and [itex]\star[/itex] indicates the convolution. I assume we should normalize h, so that we have [itex] \sum_{x,y} h(x,y) = 1[/itex] in a discrete form. And in the incoherent system, the blurred image [itex] g(x,y) = h^2 * f(x,y)^2 [/itex] Do we need to normalize h differently as? [itex] \sum_{x,y} h(x,y)^2 = 1[/itex] if doing so, it seems that the blurred image is darker in the coherent system. if not, the blurred image in the incoherent system is darker. which one is correct? thanks 


#2
Feb1113, 01:24 PM

Sci Advisor
PF Gold
P: 2,082

I'm not an optics guy, so take my comments for what they're worth.
Rather than imposing an ad hoc normalization on the jinc function, I would think that the normalization is imposed on the system level by conservation of energythat is, all of the light energy entering the aperture, minus losses in the optical system, must be present in the image. If you use a point source (uniform aperture illumination) and assume that losses are negligible, then you should get a normalization similar to your second expression but with scale factors that include the light gathering power (aperture area) of your system. After all, intensity should increase with aperture area. The expression should involve an integral, from which you can take the appropriate discrete limit. As for the intensity of a coherent system, remember that intensity is the squared magnitude of the amplitudewhich takes you back to the incoherent expression. It doesn't make sense to talk of them separately, then, since the intensity expressions are the same for both. 


#3
Feb1113, 08:21 PM

Sci Advisor
P: 5,539




#4
Feb1113, 10:16 PM

P: 37

A basic question: jinc function in coherent and incoherent system
[itex]h =\left(\begin{array}{ccc} 1/12 & 1/12 & 1/12 \\ 1/12 & 1/3 & 1/12\\ 1/12 & 1/12 & 1/12 \\ \end{array}\right) [/itex] which is normalized with [itex] \sum_{xy} h(x,y)=1 [/itex] and assume we have two points of light source, say, located in (0,0) and (0,1) and both with intensity value 1. Then the intensity of the aerial image at point (1,1), blurred due to the diffraction effect, is [itex] (1/12+1/12)^2 = 1/36 [/itex] for a coherent system and [itex] (1/12)^2 + (1/12)^2 = 1/72 [/itex] for an incoherent system so it seems that the aerial image for a incoherent system is always darker than the coherent system. And the energy conservation is not guarantted in the incoherent system. That is why I am thinking we should normalize the point spread function (jinc function) for the incoherent system in a different way. still puzzled....... 


#5
Feb1213, 06:54 PM

Sci Advisor
P: 5,539

That post was very helpful, thanks. It is often surprising that incoherent imaging yields a 'better' image than coherent imaging in your case, the amount of blur is less for incoherent imaging than coherent imaging. Similarly, the cutoff frequency for incoherent imaging is twice that for coherent imaging.
Even so, comparing incoherent and coherent imaging is not trivial. Rather than be unclear, I'll simply point to the 'gold standard' explanation, in section 6.5 http://books.google.com/books?id=ow5...herent&f=false Does this help? 


#6
Feb1213, 09:30 PM

P: 37

So when I do simulations for incoherent system, I need to normalize the point spread function as [itex] \sum_{x,y} h(x,y)^2 = 1 [/itex] which make sense because it conserve the total energy. Is that correct? thanks 


#7
Feb1313, 11:10 AM

Sci Advisor
P: 5,539

I think so part of my confusion results from applying the concept to a discretized system.
To summarize my understanding, if the coherent PSF h_c(x,y) = D/λ jinc(Dr/λ) then h_c is normalized. The incoherent PSF h_i(x,y) = (D/λ)^2 jinc^2(Dr/λ) is also normalized. The image intensity isI_i = h_c * U_o^2 for the coherent case (U_o is the object *field*) and I_i = h_i * I_o (* = convolution) for the incoherent case. I think that makes everything selfconsistent. h_c may be complex and negativevalued, but h_i must be real and positive valued everywhere. 


#8
Feb1413, 02:28 AM

P: 37




#9
Feb1413, 09:09 PM

Sci Advisor
P: 5,539

I'm confused as well I ran a quick test to compare the coherent and incoherent case, using a 3x3 blur kernel: [[0.09, 0.09, 0.09],[0.09, 0.3, 0.09],[0.09, 0.09, 0.09]] and what I assumed was the square: [[0.01, 0.01, 0.01],[0.01, 0.9, 0.01],[0.01, 0.01, 0.01]]. The convolution was done on a 10 x 10 pixel array with one central pixel set white against a black background. For the coherent case, the unblurred object white value was 255, while for the incoherent unblurred, it was 65025.
The unblurred and incoherent blurred images gave the same integrated density = 65025, while the coherent case gave an integrated density of 9700. The only thing I can think of is that the coherent blur is really a complexvalued convolution and the phase information is lost during the computation. I say this because the coherent blur h*U conserves brightness while (h*U)^2 does not. Curious.... 


#10
Feb1513, 08:23 AM

P: 37

But I dont think the phase can explain this phenomena because the problem is still there if you considre a "monophase" system...... 


#11
Feb1513, 12:16 PM

Sci Advisor
PF Gold
P: 2,082




#12
Feb1513, 01:49 PM

Sci Advisor
P: 5,539




#13
Feb1513, 02:12 PM

Sci Advisor
PF Gold
P: 2,082




Register to reply 
Related Discussions  
Question about basic concept of solving system of equations  General Math  5  
Coherent versus incoherent averaging of data  Engineering, Comp Sci, & Technology Homework  0  
Basic Two State System Question  Advanced Physics Homework  6  
Generate coherent photons from incoherent plasmons?  General Physics  1  
Number system: a basic question  General Math  7 