Archive for August, 2012

How To Measure “Photon Transfer Curve” (1) ?

Tuesday, August 28th, 2012

Over the last few years a lot of information about the Photon Transfer Curve (PTC) is being published in this blog.  The influence of many different sensor parameters on the PTC was studied and described.  For those of you that are still hungry to get more information, please remember that Jim Janesick wrote a complete book about the PTC.  It was Jim who developed this wonderful technique in the ‘80s for CCDs, but the PTC is perfectly applicable to CMOS imagers as well.  The great advantage of this measuring method is the fact that no absolute measurements of any light input is needed. 

In this blog, the application of the PTC on the data/images collected earlier will be described.  So the same images (with light on the sensor) as used in the previous blogs are being re-used.  What was done to generate these images : the sensor under test was uniformly exposed to light, while varying the exposure time.  For every exposure time applied, several images were grabbed, the average value of the obtained images as well as the temporal noise on pixel level was calculated. 

The measured signal of the sensor can be written as :

Stot = kNo + Soff    (1)

While the measured temporal noise on pixel level can be written as :

stot2 = k2sR2 + k2so2  (2)

with :

       Stot : measured output signal [DN],

       k : conversion gain [DN/e],

       No : number of optically generated electrons,

       Soff : offset signal [DN],

       stot : total temporal noise measured [DN],

       sR : temporal noise associated with the readout channel [e], also the noise floor in dark at 0 s exposure time,

       so : photon shot noise [e].

And with :

so2 = No = (Stot – Soff)/k       (3)

the total noise can be written as :

stot2 = k2sR2 + k(Stot – Soff).    (4)

The signal-to-noise ratio of the system under consideration can be written as :

                                    Stot/stot = (Stot-Soff)/( k2sR2 + k(Stot – Soff))0.5.     (5)

The obtained measurement data/images can be used to create :

1)    Standard deviation versus average effective signal, both on a log-scale (“Mean-Standard Deviation”).

 

The standard deviation is representing the temporal noise of the pixels, the effective signal is equal to the average output signal minus the DC-offset of the output signal.  The obtained curve is shown in Figure 1.

120827_blog_1

 Figure 1. Photon Transfer Curve : standard deviation versus effective output signal.

From left to right, one can recognize three parts in this PTC graph :

       A flat part in which the temporal noise of the readout channel or the noise floor at 0 s exposure time, is dominant, in the example shown this noise is equal to 100.605 DN = 4.03 DN.  This can also be derived from equation (4) and scrapping the photon shot noise :

           stot2 = k2sR2    or    stot = ksR.

       A part with increasing noise for increasing effective signal, in this region the photon shot noise is dominant, and the slope of the curve needs to be equal to 0.5 (or at least very close to this value).  The regression line drawn through this part of the curve is characterized by the relation shown in the graph.  Using equation (4) again, but this time focusing on the photon shot noise, it reduces to :

 

    stot2 = k(Stot – Soff)    or    k = 1/(Stot – Soff) at stot = 1 DN.

 

In this way the conversion gain k can be found by extrapolation of the regression line towards the crossing with the horizontal axis.  The value found for k is : 10-(0.2285/0.4781) = 0.33 DN/e.  Be aware that this method actually requires the curve to get a slope = 0.5 before any extrapolation can be done.  In this example the slope is only 0.4781.  Although this value is pretty close to 0.5, a substantial error can be introduced into the calculation of the conversion factor.  For that reason a better method to find the conversion factor is in the “Mean-Variance” technique, explained further in this blog.

 

       The part at which the noise reaches a maximum value and for any further increase in effective signal, the curve collapses.  This is the point at which the pixels reach their saturation level.  In the example shown, this happens for an effective output signal equal to 103.515 DN = 3273.0 DN.

 

2)    Variance versus average effective signal, both on a linear scale (“Mean-Variance”).

The curve obtained is illustrated in Figure 2.

120827_blog_2

Figure 2. Photon Transfer Curve : noise variance versus effective output signal.

Also this figure shows some interesting features from left to right :

       Crossing with the vertical axes (= 0 DN effective signal) for a noise variance equal to 15.61 DN.  This value corresponds to the temporal noise at very low signal values or noise floor in dark for 0 s exposure time, and actually should represent the variance of the noise associated with the readout channel.  The value found here is 15.610.5 DN = 3.95 DN.

 

       A linear part in the curve for which the relation is shown in the figure.  In this region again formula (4) is valid :

     stot2 = k2sR2 + k(Stot – Soff).

 

Taking into account the relationship of the regression line, it should be concluded that the conversion gain k (being the slope of the line) is equal to 0.23 DN/e.  This value is lower than the one previously found with the Mean-Standard Deviation method.  This is due to the fact that there is a relative large error introduced in the calculation because the slope of the noise versus signal curve is not perfectly equal to 0.5 !

 

       A saturation part with a collapsing curve for a maximum signal level of 3275 DN, being the saturation level of the pixels.

 

3)    Signal-to-noise ratio versus average effective signal, both on a log-scale (“Mean-Signal-to-Noise”).

The relation between signal-to-noise ration and effective signal is shown in Figure 3.

120827_blog_3 Figure 3. Photon Transfer Curve : signal-to-noise ratio versus effective output signal.

A similar story as before can be repeated about three different regions in the curve :

       At the very left end a first part that can be linearized is indicated by the relation between signal-to-noise and average effective signal.  In this region the noise of the readout channel is dominant, and based on equation (5), the signal-to-noise ratio can be written as :

         Stot/stot = (Stot-Soff)/(ksR)    or    ksR=Stot-Soff at Stot/stot = 1.

In this example, the readout noise is equal to 1011.388/18.941 = 3.99 DN.

       Just before saturation, a second linear part can be recognized in the region where the photon shot noise is the dominant noise source, and based on equation (5), the signal-to-noise ratio can be simplified to :

          Stot/stot = ((Stot-Soff)/k)0.5    or    k = 1/(Stot-Soff) at Stot/stot = 1.

In this example, the conversion gain is equal to 104.7219/10.452 = 0.35 DN/e.  Also this value deviates from the previously calculated values, because the extrapolation of the data can introduce some nasty errors.

       Finally the signal-to-noise ratio is becoming extremely large, due to the fact that the temporal noise of the sensor is reduced to almost zero at saturation of the pixels for a value of 103.515 DN = 3273 DN.

In conclusion, parameters such as the readout noise or noise floor in dark, the conversion gain and the saturation level can be found by means of the Photon Transfer Curve.  But the most accurate and reliable data for (although 3 different curves are used, the input data for these curves is the same) :

       the conversion gain can be obtained via the Mean-Variance curve,

       the readout noise or noise floor in dark at 0 s exposure time, can be obtained via the Mean-Standard Deviation curve,

       the saturation level can be obtained by any form of PTC.  

The PTC curve can be constructed in various ways.  In this blog uniform illuminated sensor areas at various exposure times are used, in the next blog a few alternatives will be highlighted.

Albert, 27-08-2012.

How To Measure Temporal Noise with Light ?

Monday, August 13th, 2012

 

It is again time to switch on the light and to measure the temporal noise of a sensor with light input.  Actually, do not expect that much new information with this measurement as far as noise is concerned.  Once there is light shining on the sensor and once the sensor is delivering enough signal, the photon shot noise will be the largest noise component and the measurement results will be dominated by this photon shot noise .  And it is known that the photon shot noise (expressed in electrons) has a square-root relation with the amount of electrons generated in the pixels.  There is no real need to proof the Poisson statistics again, others have done this already millions of time. 

When measuring the temporal noise, it must be admitted that most interesting information is already found in the sessions without light input and with an exposure time equal to 0 s.  Nevertheless, to complete the story about temporal noise measurements, also the evaluation of temporal noise with light will be discussed in this blog. 

Figure 1 repeats the average output signal of the sensor, with light input, as a function of the exposure time.  This information is already shown in an earlier blog.

next_blog_1

 

Figure 1. Average output signal, with light input, as a function of the exposure time.

As can be expected, the various colours in Figure 1 represent the four colour channels of the device (red, green in red line, green in blue line and blue).  The following remarks can be made :

       Offset at 0 s exposure time, being equal for all colour channels, being 820 DN,

       Saturation level being equal for all colour channels, the saturation levels are equal to 4095 DN.  So the saturation level of the signals is determined by the maximum level of the ADC (being 212 -1 DN = 4095 DN).  With the correction of the offset, the sensor delivers maximum 3275 DN,

       Sensitivity of the blue channel is lowest, next is red, and the most sensitive channels are the two green ones.  This observation is the case for the light input parameters indicated, 5600K colour temperature, and 5 lux on the sensor.

Figure 2 illustrates the measured temporal noise under the same light and sensor conditions as the one used in Figure 1.

next_blog_2

 

Figure 2. Temporal noise of the sensor with light input, as a function of the exposure time.

Also from this graph a few interesting remarks can be made :

       For an exposure time 0 s, the minimum noise is equal to 3.92 DN.  This is the noise floor of the sensor in dark at 0 s exposure time,

       For increasing exposure times, the noise increases as well, but in a non-linear way.  This is due to the fact that for larger exposure times, the noise is dominated by the photon shot noise, and the latter has a square-root relation with the average signal (in the charge domain),

       Because the sensitivity of the green channels is the largest of all, and the sensitivity of the blue channel is the lowest, the same is true for the noise levels at a particular exposure time,

       As soon as the signals saturate, the temporal noise is quickly dropping to a very low level.  This is due to the fact that the ADC is saturating and fixing the output value of the signal to 4095 DN and any further noise content in the analog signal does not show up in the digital representation of the output signal.

Figure 3 shows the temporal noise histogram, obtained at an exposure time of 0.032 s, at which the green channels reach 25 % of saturation.

next_blog_3

 

Figure 3. Noise histogram for the sensor with light, only one of the green channels is shown.

 At the left side in Figure 3, the Poission distribution of the photon shot noise can be recognized.  The tail at the right side in the histogram is representing the pixels of which the temporal noise is dominated by 1/f and RTS noise.

 

Albert, 13-07-2012.