Друк

[Гість]

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

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

Увійти

[Новини]

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

Наш канал

Наші сайти

Астро-сайти

Програми

Наш канал

Наші сайти

Астро-сайти

Програми

[SP]

. Предоставить ряд из 426 чисел S/N, с определением как считалось. Результирующая функция, кстати,  необязательно имеет один экстремум.

Это интересно. Результаты можно будет выложить потом? Кстати, а возможно ли создание утилиты, которая будет анализировать  серию кадров на предмет нужности их для складывания (с анализом SNR, с учётом экстремумов и алгоритмов сигма-клиппинга?)

[MaxMan]

Обследовал 10-минутные L кадры (319 штук). Заявленного бешеного изменения S/N (в 6 раз) не увидел. Что и как ты мерял?

Не помню что я вчера мерял. Мне вчера ночью Yvk прислал свои исходники и я поменял расположение участков под его кадр. С таким расположением получается максимум 31,67 а минимум 4,16. Разница в 7,6 раза!!!

[tlgleonid]

Пусть у нас есть N кадров. На каждом кадре есть некое значение яркости полезного сигнала Si и величина шума Дi. i=1,2..N Для конкретного кадра соотношение сигнал/шум равно  Si/Дi. Если мы складываем эти N кадров, то Соотношение сигнал/шум суммы равно  Sum(Si,i=1..N)/Sqrt(Sum(Дi^2,i=1..N)).
Но ведь каждый кадр можно масштабировать, то есть сигнал и шум умножать на определенную константу u, то есть фактически складывать кадры со своими весовыми коэффициентами. Как их найти? Нужно искать их такими, что бы итоговое соотношение сигнал/шум было максимальным.
Пусть у нас есть сумма k кадров. соотношение сигнал/шум равно Z/sqrt(M)   Z=Sum(uj*Sj,j=1..k)  M=Sum((uj*Дj)^2,j=1..k). При добавлении еще одного кадра (сигнал a, шум r) с новым весовым коэффициентом u получим новое соотношение сигнал/шум равное f(u)=(Z+u*a)/sqrt(M+r^2). Максимум f(u) найти просто. Найдем u при котом f'(u)=0, f''(u)<0. Первое равенство достигается, когда a*M=u*r^2*Z или u=a*M/(Z*r^2). Несложно убедится, что f''(u)<0 всегда будет истинно, поскольку все величины положительны.

Таким образом соотношение сигнал/шум правильно сложенных N кадров всегда будет не хуже, чем правильно сложенная сумма любого подмножества этих кадров. Складывать нужно все кадры, но правильно выбирать весовые коэффициенты.

Если же коэффициенты выбирать неправильно, например просто тупо складывать хорошие кадры с плохими, можно получить и ухудшение.
Пример:
Пусть на всех кадрах шум равен 1, а сигнал на 130 кадрах равен 8, а на 131-ом только 2
Тогда сумма 130 кадров даст соотношение сигнал/шум, равный 130*8/sqrt(130)=91.214, а если тупо добавить плохой 131 кадр, то (130*8+2)/sqrt(131)=91.04  окажется меньше. Но если плохой кадр взять с весом u=2*130/(130*8 )=1/4, то соотношение сигнал/шум станет равным (130*8+2/4)/sqrt(130+(1/4)^2)=91.236>91.214 То есть плохой кадр все равно добавляет полезную информацию.

Подякували

3 :serhio_67, MaxMan, Yvk

[MaxMan]

Спасибо Лёня!  Практика подтверждается теорией. Собственно именно по этой причине Ваня Ионов и написал свой фитстеккер. Именно там он и реализовал эти весовые коэффициенты.

Подякували

1 :tlgleonid

[M_M]

MaxMan, а ты уверен в том, что сравнивать S/N B и L каналов в данном контексте корректно? И еще, из скриншотя неясно, что за программа и как S/N считается. Расположение красных и зеленых квадратиков весьма подозрительно.

tlgleonid, в случае применения весовых коэффициентов всегда можно на плохой кадр выдать коэффициент, равный нулю, для того, чтобы доказать тезис "всегда можно сделать не хуже". Вопрос в том, применяется ли весовой коэффициент в данном случае и о его реальном значении. А также о самих понятиях сигнал и шум. Есть подозрение, что четкого определения у многих просто нет.

f(u)=(Z+u*a)/sqrt(M+r^2)

Кажется, тут пропущен коэффициент u перед r².

По поводу примера. Просто для оценки величины коэффифициента, у нас несколько другой случай. Сигнал постоянный, растет шум вместе с засветкой. Предположим, 130 кадров с S=8, N=1, 131-й кадр с S=8, N=4. какой коэффициент u для оптимального сложения 131-го кадра?

Подякували

1 :tlgleonid

[SP]

Про коэффициенты "весомости" известно давно. А вот где их взять для каждого кадра из серии, и так чтобы это было оптимально, не ясно.

[Jeff]

Про коэффициенты "весомости" известно давно. А вот где их взять для каждого кадра из серии, и так чтобы это было оптимально, не ясно.

Например: http://pixinsight.com/doc/tools/ImageIntegration/ImageIntegration.html#description_002

Подякували

1 :M_M

[M_M]

Интересно, что за уфон на одном из кадров?

[MaxMan]

http://www.astroclub.kiev.ua/forum/index.php?topic=36682.msg440348#msg440348 - я уже публиковал. 

[tlgleonid]

tlgleonid, в случае применения весовых коэффициентов всегда можно на плохой кадр выдать коэффициент, равный нулю, для того, чтобы доказать тезис "всегда можно сделать не хуже". Вопрос в том, применяется ли весовой коэффициент в данном случае и о его реальном значении. А также о самих понятиях сигнал и шум. Есть подозрение, что четкого определения у многих просто нет.

Весовой коэффициент будет равным нулю только в том случае, когда полезный сигнал просто отсутствует.  Еще одна проблема возникает тогда, когда хотя бы на одном из кадров есть только полезный сигнал и полностью отсутствует шум. Сигнал - это значения яркости измеряемого объекта. Модель не зависит от выбранного на кадре объекта.

f(u)=(Z+u*a)/sqrt(M+r^2)

Кажется, тут пропущен коэффициент u перед r².

Верно! Подправил.

По поводу примера. Просто для оценки величины коэффифициента, у нас несколько другой случай. Сигнал постоянный, растет шум вместе с засветкой. Предположим, 130 кадров с S=8, N=1, 131-й кадр с S=8, N=4. какой коэффициент u для оптимального сложения 131-го кадра?

Z=130*8, M=130, a=8, r=4, значит u=8*130/(130*8*4^2)=1/16
Вообще, для упрощения нам ничто не мешает выровнять кадры как по уровню сигнала (меняются шумы), так и по уровню шумов (меняется уровень сигнала).
Если кадры выровнены по уровню сигнала, то весовые коэффициенты для каждого кадра обратно пропорциональны квадрату уровня шумов.
Если кадры выровнены по уровню шумов, то весовые коэффициенты пропорциональны уровню сигнала.
Коэффициент пропорциональности может быть выбран произвольно.

[tlgleonid]

Спасибо Лёня!  Практика подтверждается теорией. Собственно именно по этой причине Ваня Ионов и написал свой фитстеккер. Именно там он и реализовал эти весовые коэффициенты.

Юра, а можно поподробнее про этот фитстеккер и его матмодель. Об программе я что-то краем уха слыхал, но думал, что это какая-то упрощенная реализация ириса.
Вообще странно, что идею сложения с весовыми коэффициентами не реализовали ни в Ирисе, ни в Ауделе, ни в Максиме. А ведь идея то лежит на поверхности и формула для коэффициентов вычисляются очень просто. Деконволюция куда сложнее. 

[MaxMan]

По этому вопросу лучше к разработчику. Он теорию знает лучше меня в разы. Или полистать нему на Звездочёте. Я пользуюсь 12-й версией.
Очень много пользователей уже оценили прелесть разработки которая улучшает сумму на 20-30% по соотношению с/ш в сравнении с ирисами, максимами и прочим софтом.

[M_M]

Если кадры выровнены по уровню сигнала, то весовые коэффициенты для каждого кадра обратно пропорциональны квадрату уровня шумов.

Вот и добрались до сути...

Уровень сигнала условно можно считать постоянным, редко кто снимает объекты от горизонта до горизонта через зенит. В основном меняется засветка. А шум от засветки пропорционален корню квадратному от ее амплитуды.

Соответственно, весовой коэффициент обратно пропорционален уровню засветки.

Отсюда следует, что для "вытаскивания" фоток из засветки нужно набирать количество "експочасов" пропорционально уровню засветки для того же результата.

Например, в Киеве и на даче у меня уровень засветки отличается в ~250 раз. Итого, час экспозиции при хорошей астропогоде на даче эквивалентен паре-тройке лет из города, учитывая продолжительность астроночи, наличие конкретного объекта в доступном секторе, погоду итд.

Хорошо бы вычленить это обсуждение в отдельную тему

[tlgleonid]

Например, в Киеве и на даче у меня уровень засветки отличается в ~250 раз. Итого, час экспозиции при хорошей астропогоде на даче эквивалентен паре-тройке лет из города, учитывая продолжительность астроночи, наличие конкретного объекта в доступном секторе, погоду итд.

Тут все не так однозначно. Если я луплю по переменкам с выдержками в 30 сек, то разница между Киевом и дачей для меня малосущественна. Если кто-то снимает в Аш-Альфа или О3, то там разница от засветки также будет не так велика

[MaxMan]

Галактика М33 в созвездии Треугольник (водородная версия)

CCD: AtiK 16HR
Exposure time:    Ha: 99x1200sec, 1x1 bin
                           
Exposure date: 06.10.2014 - 15.10.2014
Location: Pesochin Observatory "Astrobudka", Kharkov, Ukraine
Filter: Astronomik 1.25" Ha  CCD filter
Instrument: Newton 190mm (f-760mm) + Parracor Type 1 = f-875mm
Mount: WS-180, control block "EQDrive Standard"
Guide: Astro-Tech 72mm, f-400mm + QHY5-II

Большой размер и описание тут - http://astrobudka.blogspot.com/2014/10/33.html

Подякували

5 :OlegBr, trojn(Юрий), vakula, ARTY, Nao

[MaxMan]

И вот такая вольность автора - наложение цвета из прошлогоднего снимка на водородный кадр. 

Подякували

1 :Garuda

[malex]

ЧБ с узкополосниками  Фото уже претендует на некую научную, а не только художественную ценность. Очень интересная работа.

[Alone_Ghost]

Юрчик, а ты все разрываешь моск  астробудка-это сила

[MaxMan]

Галактика М33 (новый взгляд на старый объект)

CCD: AtiK 16HR
Exposure time:    L: 14x600sec, 1x1 bin
                            Ha: 99x1200sec, 1x1 bin
                            R: 14x280sec, 1x1 bin
                           G: 14x320sec, 1x1 bin
                           B: 14x360sec, 1x1 bin
Exposure date: 24.10.2014 - 25.10.2014
Location: Pesochin Observatory "Astrobudka", Kharkov, Ukraine
Filter: Astronomik 1.25" L R G B Ha  CCD filter
Instrument: Newton 190mm (f-760mm) + Parracor Type 1 = f-875mm
Mount: WS-180, control block "EQDrive Standard"
Guide: Astro-Tech 72mm, f-400mm + QHY5-II

Большой размер и разные вариации тут - http://astrobudka.blogspot.com/2014/10/33_25.html

Подякували

11 :Garuda, serhio_67, bigol, Diskus, Vladimirovich, OlegBr, Вадим, ds40a, novik, vakula, Tristan

[MaxMan]

Снимаем на полурабочем сетапе сквозь дым от костров под песню "Дым сигарет с ментолом" 

Галактика в Пегасе NGC 7741

CCD: AtiK 16HR
Exposure time:    L: 8x600sec, 1x1 bin
                         
Exposure date: 30.10.2014
Location: Pesochin Observatory "Astrobudka", Kharkov, Ukraine
Filter: Astronomik 1.25" L CCD filter
Instrument: Newton 190mm (f-760mm) + Parracor Type 1 = f-875mm
Mount: WS-180, control block "EQDrive Standard"
Guide: Astro-Tech 72mm, f-400mm + QHY5-II

Большой размер и описание тут - http://astrobudka.blogspot.com/2014/10/ngc-7741.html

Подякували

3 :OlegBr, Vladimir_sx, Nao