esbo_ds/ESBO-ETC
esbo_ds/ESBO-ETC
11
0
Fork 1
Code Issues 1 Pull Requests Releases 1 Wiki Activity

Bugfix: calculate exposure time correctly

Browse Source
This commit is contained in:
Lukas Klass 2020-05-19 10:45:42 +02:00
parent a2fbbfc49d
commit 4ff5f63992
1 changed files with 2 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
esbo_etc/classes/sensor/Imager.py
View File

@ -144,17 +144,14 @@ class Imager(ASensor):
signal_current, background_current, read_noise, dark_current = self.__exposePixels() signal_current, background_current, read_noise, dark_current = self.__exposePixels()
# Calculate the electron currents for all pixels # Calculate the electron currents for all pixels
signal_current_tot = signal_current.sum() signal_current_tot = signal_current.sum()
background_current_tot = background_current.sum()
read_noise_tot = read_noise.sum()
dark_current_tot = dark_current.sum()
# Fix the physical units of the SNR # Fix the physical units of the SNR
snr = snr * u.electron ** 0.5 snr = snr * u.electron ** 0.5
# Calculate the ratio of the background- and dark-current to the signal current as auxiliary variable # Calculate the ratio of the background- and dark-current to the signal current as auxiliary variable
current_ratio = (background_current_tot + dark_current_tot) / signal_current_tot current_ratio = (background_current.sum() + dark_current.sum()) / signal_current_tot
# Calculate the necessary exposure time as inverse of the CCD-equation # Calculate the necessary exposure time as inverse of the CCD-equation
exp_time = snr ** 2 * ( exp_time = snr ** 2 * (
1 + current_ratio + np.sqrt((1 + current_ratio) ** 2 + 4 * read_noise_tot ** 2 / snr ** 2)) / ( 1 + current_ratio + np.sqrt((1 + current_ratio) ** 2 + 4 * (read_noise ** 2).sum() / snr ** 2)) / (
2 * signal_current_tot) 2 * signal_current_tot)
# Print information # Print information
if exp_time.size > 1: if exp_time.size > 1: