Друк

[Гість]

Ласкаво просимо, Гість. Будь ласка, увійдіть або зареєструйтеся.
23 Листопада 2024, 00:42:13

1 година 1 день 1 тиждень 1 місяць Назавжди

Увійти

[Новини]

•    Новини
•      Астрономія
•      Астрозаходи
•    Життя клуба
•      Новини
•      Проекти
•      Заходи
•      Звіти про спостереження
•      Документи
•    Статті
•      Обладнання та телескопи
•      Спостереження
•      Астрофотографія
•      Астрософт
•      Астрономія в Україні
•      Астрономи шуткують
•      Інше
•     
•      Файловий архів

Наш канал

Наші сайти

Астро-сайти

Програми

Наш канал

Наші сайти

Астро-сайти

Програми

[Серега]

Трапилось мені посилання на методику визначення Gain та readnoise камери.
Стало цікаво, спробував. паралельно трапилось ще пару методик.
І все вони дають різні результати  деякі дають близькі, а деякі взагалі і близько не співпадають.

Приведу з них формули розрахунку.

Вони потребують двох (або більше) флетів з різним насиченням. Я намагався зробити з мінімальним та максимальним. Робив навіть декілька по всьому діапазону.

Який з методів більш правильний??

Ось з хелпа програми Iraf:
https://iraf.net/irafhelp.php?val=obsutil.findgain&help=Help+Page

flatdif = flat1 - flat2

    zerodif = zero1 - zero2


       gain = ((mean(flat1) + mean(flat2)) - (mean(zero1) + mean(zero2))) /
         ((sigma(flatdif))**2 - (sigma(zerodif))**2 )


   readnoise = gain * sigma(zerodif) / sqrt(2)

zero = bias
sigma = Std Div


І ще одна:
https://www.photometrics.com/learn/imaging-topics/gain

A simple method to calculate the system gain is shown below:

Collect a bias image (zero-integration dark image) and label it “bias”.
Collect two even-illumination images and label them “flat1” and “flat2”.
Calculate a difference image: diff = flat2 – flat1.
Calculate the standard deviation of the central 100 x 100 pixels in the difference image.
Calculate the variance by squaring the standard deviation and dividing by 2 (variance adds per image, so the variance of the difference image is the sum of the variance of flat1 and flat2).
Calculate a bias-corrected image by subtracting the bias from one of the flat images and label it corr: corr = flat1 – bias.
Obtain the mean illumination level by calculating the mean of the central 100 x 100 region of the corr image.
The mean divided by the variance equals the gain: gain = mean /variance.

Ось метода по Крейгу:
https://www.cloudynights.com/articles/cat/column/fishing-for-photons/signal-to-noise-part-3-measuring-your-camera-r1929

[mpyat2]

Второй метод, мне кажется, дает лучший результат. Я делал для своего кенона (когда разбирался и переводил руководство по фотометрии, вот, кстати, кусок об этом https://drive.google.com/open?id=1S_TxAHShy2DFs1oSnmExxnR3uFy9eoo8 , там я и свои данные вставил).
Потом то же делал для asi120mm.
Между прочим, SharpCap имеет функцию измерения параметров сенсора камеры, в том числе и гейна, работает очень классно (мне так показалось). Меряет гейн, шум считывания и некоторые другие параметры.

Подякували

1 :Серега

[Серега]

Мені не потрібен кращий результат. Потрібна правильна метода-формула 
Бо я вже перестаю розуміти взагалі принцип - заварив кашу в голові


в мене камера не Байеровську структуру має, тому ніякий софт з тим не може справитись
Але навчився в Iris відокремлювати 50% зеленого каналу, з ним і працював.

[Серега]

2. Виміряйте медіану всіх пікселів; це і буде значення системного зміщення.

саме медіану? чи mean?

[Серега]

Второй метод, мне кажется, дает лучший результат. Я делал для своего кенона (когда разбирался и переводил руководство по фотометрии, вот, кстати, кусок об этом

так, формула така сама і доресі, вона мені теж більше сподобалась. Принаймні тим що в статті є посилання на літературу.  Тільки у вашому руководстві флети однакові по часу експозиції та насиченності. А всюди  або не вказано, або вказано що має бути різниця.
Що ж, спробую з однаковими параметрами.

[Серега]

Що я роблю не так? Зробив все по написаному для ISO200, отримав результат біля 20e/ADU, а для ISO2000 біля 2,5e/ADU 
Не може такого бути, ну чесно

там точно для обох флетів потрібна однакова експозиція?
Виміряне системне зміщення якось застосовується до флетів та біасів?

[mpyat2]

2. Виміряйте медіану всіх пікселів; це і буде значення системного зміщення.

саме медіану? чи mean?

Медіана тут краще працює.

[mpyat2]

Що я роблю не так? Зробив все по написаному для ISO200, отримав результат біля 20e/ADU, а для ISO2000 біля 2,5e/ADU 
Не може такого бути, ну чесно

там точно для обох флетів потрібна однакова експозиція?
Виміряне системне зміщення якось застосовується до флетів та біасів?

Серьога, давно я це робив, спробую повторити, тоді щось скажу, бо підзабув, чесно.

Щодо ваших результатів,то виходить, що ISO 2000 має приблизно в 10 разів більший коефіцієнт підсилення, ніж 200, але абсолютні значення... Якісь завеликі, дійсно.

[mpyat2]

Тільки у вашому руководстві флети однакові по часу експозиції та насиченності. А всюди  або не вказано, або вказано що має бути різниця.
Що ж, спробую з однаковими параметрами.

Чекайте, в другому способі ми фіксуємо ISO і знімаємо пари флетів з різними експозиціями. Тож насиченість буде різна. Для кожної експозиції два однакових флети. Один з флетів з пари використовуємо для вимірювання середнього ADU. Вимірюємо в маленькій ділянці, щоб вплив нерівномірності освітлення був мінімальним.
Далі вимірюємо дисперсію. Тут ми можемо взяти різницю між однаковими флетами, що компенсує нерівномірності, що можуть вплинути на вимір дисперсії. Але при цьому квадрат дисперсії зростає в два рази відносно одиничного зображення, тому ділимо квадрат дисперсії на два.

А в методі IRAF (перший спосіб) те ж саме, але вимірюється по одній експозиції. Тут нам треба прибрати нульовий рівень (біас), а з дисперсії прибрати частку, яка вноситься шумом зчитування. Однак вимірювання по одній точці не дуже точне.

[Серега]

По IRAF треба флет калібрувати біасом, чи просто знайдене нульове значення відняти?

Що значить "по точці міряти"? Я міряю 100х100

Роблю все згідно формули, пари біасів, флетів, відокремлюю канал, заміряю все що треба, і отримую те що отримую..
Ну не може ж бути така матриця?

[Серега]

" другий спосіб" то який? Що в першому пості моєму другий? Чи в вашому пдф, той що по нахилу графіка? Якщо пдф, то зрозуміло, щвісно, що різні експозиції.

[Серега]

Розібрався Смішной вийшло. Математичні операції проводив не з кадрами, а з цифрами... Тому і отримував хто зна що.

Подякували

1 :mpyat2

[mpyat2]

Щоб трохи підсумувати.
Ідея в тому, що стандартне відхилення сігма кількості фотонів, зафіксованих сенсором, тобто електронів, що реєструються, дорівнює корню квадратному від середнього значення  цієї кількості.