Mesure du gain

Calcul du gain a partir des images de labo et fait par procedure inter

 
Calcul du gain:
 
Pour niveau (median) = 996
gain        = 4.03 [ADU/e-]
conv factor = 0.25 [e-/ADU]

Pour niveau (median) = 1986
gain        = 4.10 [ADU/e-]
conv factor = 0.24 [e-/ADU]

Pour niveau (median) = 3018
gain        = 9.49 [ADU/e-]
conv factor = 0.11 [e-/ADU]

local image(6)=""
image(1)="/gls/data/raw/TestBigEye/Linearity/testLinearity_016280.fits" ! median 990
image(2)="/gls/data/raw/TestBigEye/Linearity/testLinearity_016540.fits" ! median 1003

image(3)="/gls/data/raw/TestBigEye/Linearity/testLinearity_034296.fits" ! median 1957
image(4)="/gls/data/raw/TestBigEye/Linearity/testLinearity_035656.fits" ! median 2026

image(5)="/gls/data/raw/TestBigEye/Linearity/testLinearity_053872.fits" ! median 3070
image(6)="/gls/data/raw/TestBigEye/Linearity/testLinearity_056096.fits" ! median 3193

local s(3)
local moy(2)
local b ! (rms)
local smoy ! (signal moyen)
local g ! (gain)
local c ! (conversion factor)

nx(1:2)=10 
ny(1:2)=10 
xstart(1:2)=1280/2 
ystart(1:2)=1024/2 
xstep(1:2)=1 
ystep(1:2)=1

local inda indb
write "Calcul du gain:"
write ""
do i=1,3
	inda=(i-1)*2+1	
	indb=inda+1
	
	fimage /read image(inda) [3]
	[1]=[3]    /world
	moy(1)=sum([1])/nx/ny
	
	fimage /read image(indb) [3]
	[2]=[3]    /world
	moy(2)=sum([2])/nx/ny
	
	b=stddev([1]-[2])/sqrt(2)
	smoy=(moy(1)+moy(2))/2
	c=(b*b)/smoy
	g=1/c
	
	write "Pour niveau (median) = " itoa(median(vec([1],[2])))
	write "gain        = " format(g,"f4.2") " [ADU/e-]"
	write "conv factor = " format(c,"f4.2") " [e-/ADU]"
	write ""
	
enddo