Оптимальное сложение снимков + Программа FITStacker

[Ivan7enych]

Я тут проясняю для себя некоторые пробелы. Возможно вопрос покажется странным, но... как ФитСтэкер считает S/N? Предположу, что среднее значение в красном прямоугольнике делит на стандартное отклонение в зеленом, так?

Почти.
(Avr(red rect) - Avr(green rect)) / Stdev (green rect)

При этом есть нюансы, Stdev (green rect) считается с выбрасыванием значений вне 3х сигм (выкидываются горячие пиксели, случайно попавшие в прямоугольник слабые звезды). Иначе любой очень горячий пиксель попавший в прямоугольник мог "увеличить" шум фона раза в 2.

[Юрий Звёздный]

Почти.
(Avr(red rect) - Avr(green rect)) / Stdev (green rect)

А зачем из красного вычитается avr(green)? Что бы измерять отношение сигнала к шуму ,а не сигнал+шум к шуму?

При этом есть нюансы, Stdev (green rect) считается с выбрасыванием значений вне 3х сигм (выкидываются горячие пиксели, случайно попавшие в прямоугольник слабые звезды). Иначе любой очень горячий пиксель попавший в прямоугольник мог "увеличить" шум фона раза в 2.

Интересно, результаты замеров Stdev(green) весьма значительно различаются в зависимости от того где я провожу измерения в Пиксе или в ФитСтекере.

[Ivan7enych]

Почти.
(Avr(red rect) - Avr(green rect)) / Stdev (green rect)

А зачем из красного вычитается avr(green)? Что бы измерять отношение сигнала к шуму ,а не сигнал+шум к шуму?

Чтобы не зависеть от яркости фона.

К примеру, фон кадра находится примерно на 1000 (из-за оффсета камеры или яркости неба, неважно), то есть Avr(green) = 1000
Объект (центр галактики) ярче фона на 100, то есть Avr(red rect) = 1100
Яркость галактики и есть разность между этими 2мя средними.

[Ivan7enych]

При этом есть нюансы, Stdev (green rect) считается с выбрасыванием значений вне 3х сигм (выкидываются горячие пиксели, случайно попавшие в прямоугольник слабые звезды). Иначе любой очень горячий пиксель попавший в прямоугольник мог "увеличить" шум фона раза в 2.

Интересно, результаты замеров Stdev(green) весьма значительно различаются в зависимости от того где я провожу измерения в Пиксе или в ФитСтекере.

Лучше поиграйся прямоугольником в писке (или максиме) и фитстекере, и оцени на сколько в рамках одной программы сильно прыгает оценка шума в зависимости от выбора прямоугольника (попал ли в него горячий пиксель или слабая звезда). Понятно, чем меньше прямоугольник, тем хуже точность оценки шума, поэтому желательно сохранять размер одинаковым и только менять его расположение в кадре.

Вот попробовал оценку шума в максиме и фитстекере на одном и том же кадре. Слева в прямоугольник попал горячий пиксель, справа фон чистый. В максиме оценка шума взлетела от 40 до 175, в фитстекере от 36 до 39. Если бы фитстекер оценивал шум так же как в максиме, этот кадр был бы практически полностью выкинут из сложения.

[Alex1958]

При этом есть нюансы, Stdev (green rect) считается с выбрасыванием значений вне 3х сигм (выкидываются горячие пиксели, случайно попавшие в прямоугольник слабые звезды). Иначе любой очень горячий пиксель попавший в прямоугольник мог "увеличить" шум фона раза в 2.

Интересно, результаты замеров Stdev(green) весьма значительно различаются в зависимости от того где я провожу измерения в Пиксе или в ФитСтекере.

Лучше поиграйся прямоугольником в писке (или максиме) и фитстекере, и оцени на сколько в рамках одной программы сильно прыгает оценка шума в зависимости от выбора прямоугольника (попал ли в него горячий пиксель или слабая звезда). Понятно, чем меньше прямоугольник, тем хуже точность оценки шума, поэтому желательно сохранять размер одинаковым и только менять его расположение в кадре.

Вот попробовал оценку шума в максиме и фитстекере на одном и том же кадре. Слева в прямоугольник попал горячий пиксель, справа фон чистый. В максиме оценка шума взлетела от 40 до 175, в фитстекере от 36 до 39. Если бы фитстекер оценивал шум так же как в максиме, этот кадр был бы практически полностью выкинут из сложения.

Любопытно. Поэтому иногда суммы в Максиме приводят в недоумение .
А с DSS  такого сравнения нет? Хотя, как там можно оценить SNR...
Алгоритм Фитса понятен, а DSS - как черный ящик, складывает в основном хорошо, даже иногда не хуже Фитса, но сомнения всегда остаются...

[Ivan7enych]

Вылезла проблема при сложении кадров с разных телескопов, которую пришлось решать в программе.

При сложении черные поля выкидываются, но выравнивание кадров интерполирует черные поля и несколько размазывает их на соседние пиксели, в результате на сумме вылезают черные окантовки.

Проблему решил добавив новый параметр отступа от черных полей (задается в пикселях от 0 до... ).

1я картинка - сумма всех кадров

2я - для примера сумма 2х кадров с отступом = 0
3я - сумма тех же 2х кадров с отступом = 5

[Ivan7enych]

На 3й картинке выше осталась видна 2я проблема, черные поля ушли, но градиент на одном кадре не дал выровнять фон по краям - где-то край получился темнее, где-то светлее.

Я добавил галку - выравнивать градиент на одиночных кадрах перед сложением, стало заметно лучше. Для корректного убирания градиентов надо расставить 2-3 прямоугольника по фону в районе краев самого маленького кадра в серии. На картинке хорошо видно как их стоит поставить.

Сумма всех кадров (с 3х телескопов) теперь выглядит намного приличнее.

Архив с новой версией программы -
http://ivan.starlake.ru/FITStacker13.zip
Всех с наступающим!

[4D]

Обрастает программка. Растёт

[mo]

Супер! Спасибо.

А то приходилось складывать по частям, потом кропить, потом доскладывать. Ну и за убиралку градиентов спасибо. А то через Remove All приходилось в файлы *_gr, а потом их складывать.

С наступающим!
Пойду качать прогу.

[Dmitrii]

Архив с новой версией программы

Большое спасибо и с наступающим !!!

[Юрий Торопин]

Архив с новой версией программы -
http://ivan.starlake.ru/FITStacker13.zip
Всех с наступающим!

Спасибо, с наступающим!

[buklikovx]

За галку "Remove gradient" большое спасибо, Иван! Для исходников, полученных на плохом небе с сильным, неравномерным и подвижным градиентом, как у меня, очень полезная фича. С Новым Годом всех форумчан! 

[PavelGhost]

А что она делает эта галка ?
Всех с НасТупающим.

[Vladimir3621]

Докладываю: галка "Убрать градиент" убирает градиент. Ваш Кэп 
С Новым годом!

[PavelGhost]

Докладываю: галка "Убрать градиент" убирает градиент.

А для чего тогда кнопы Flatt и Grad ?

[Ihtamnet II]

Я что-то тоже не совсем понял из описания. Раньше убрать градиент на каждом снимке можно было в меню кнопки List - Remove gradient All. После этого все кадры сохранялись и их нужно было сложить отдельно. А теперь что делает кнопка Grad а что галка Remove Gradient?

[Ivan7enych]

Я что-то тоже не совсем понял из описания. Раньше убрать градиент на каждом снимке можно было в меню кнопки List - Remove gradient All. После этого все кадры сохранялись и их нужно было сложить отдельно. А теперь что делает кнопка Grad а что галка Remove Gradient?

Remove Gradient - выравнивает градиент на одиночных кадрах в процессе сложения, без необходимости 2х стадий обработки (сперва List - Remove gradient All, потом сложение)

Grad - как и раньше, работает с одним открытым кадром.

[PavelGhost]

Очень долго работает.

[mike]

Пережженые звезды. В ИРИСе понятно, как выйти из этого.
А как здесь можно "дожечь" центра звезд. Или только в фотошопе?

[PavelGhost]

Ну в том же ИРИСе можно.