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

[NKV]

Иван, есть ли возможность осуществить/добавить алгоритм сложения SD Mask ?

Дело в том что когда космики или медленно летящий объект засоряет один и тот же пиксель на нескольких кадрах Сигма клипинг несправляется. Вернее либо (при низком значении сигмы) режет полезный сигнал, либо (при высоком значении сигмы) пролазят эти бяки. А SD Mask легко, т.к. делает это в несколько подходов.

Пример: во вложении кино "Булыжник в Плеядах" короткий AVI из 7и кадров.
1й кадр SDMask (MaximDL)
2-5й кадры сырые после калибровки (всего таких 17штук х 300s сложено)
7й кадр Sigma (FITStackerV7)

Разница между 1м и 7м колосальна!!!

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

Сугубо ИМХО, но как по мне, след от будыжника - это ж даже прикольно, не надо его давить уж совсем!
Если что, штамп в фотошопе справлялся и не с такими булыжниками.
Тут всё-таки важен максимальный сигнал/шум после сложения серии. 

[NKV]

Сугубо ИМХО, но как по мне, след от будыжника - это ж даже прикольно

Прикольно пока Ч/Б... а в цвете это новогодняя герлянда или НЛО с мигалками
А на большой матрице таких прикольных с десяток одновременно фигачит в разные стороны.
Но основной нюанс не в булыжниках, а в том что подлые космики на длительных выдержках... вобщем пословица про бомбу и воронку им побарабану.

Тут всё-таки важен максимальный сигнал/шум после сложения серии. 

Одно другому разьве мешает?

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

Тут всё-таки важен максимальный сигнал/шум после сложения серии. 

Одно другому разьве мешает?

Ничего не могу сказать, не знаю, как устроен SD Mask алгоритм...
Возможно, Ivan7enych прольёт свет на этот вопрос
Думаю, что скомбинировать можно почти всё что угодно

[Ivan7enych]

Тут всё-таки важен максимальный сигнал/шум после сложения серии. 

Одно другому разьве мешает?

Ничего не могу сказать, не знаю, как устроен SD Mask алгоритм...
Возможно, Ivan7enych прольёт свет на этот вопрос
Думаю, что скомбинировать можно почти всё что угодно

Я могу сделать по простому, суммирование с клиппиногом посчитать средним файлом, и прогнать суммирование с клиппиногом на основе этого нового файла. Время работы возрастет в 2-3 раза (для 2-3 итераций). Работать будет, но не быстро.

Я пока собираю все предложения, через неделю уеду в командировку и там по вечерам займусь доработками.

[NKV]

Я могу сделать по простому, суммирование с клиппиногом посчитать средним файлом, и прогнать суммирование с клиппиногом на основе этого нового файла.

Немного не то... есть нюансы: http://www.rc-astro.com/resources/sigma_reject.html

The Doubly-Iterated Sigma-Reject Algorithm

This is an extension of the standard sigma-reject algorithm, and is what is implemented in SigmaReject. It produces both better rejection of artifacts, and better signal-to-noise ratio. Rejection properties are especially improved for relatively faint artifacts, such as satellite trails.

   1. Compute the mean and standard deviation as above.
   2. Reject spurious values and compute a new mean value as above.
   3. Compute a new standard deviation based on this new mean, and using only the accepted pixel values from step 2.
   4. Using this new, better estimate of the standard deviation, repeat the pixel rejection process, and generate a new mean value from the accepted pixels.
   5. Repeat for every pixel in the final image.

The idea behind iterating twice is that the original standard deviation computation is wrong in the presence of an artifact; it is too large. Re-computing the standard deviation after rejecting spurious pixel values results in a better estimate of the "true" standard deviation, which can then be applied to the second rejection process.

Количество итераций задаётся.

Вот фильмец про всякие сигмы. там вся теория на примерах (и не только в максиме) изюмительно разжевана:
http://files.mail.ru/Q8277X
57Mb

[Ivan7enych]

Спасибо, идею вспомнил.

[NKV]

Ещё предложение:
1) Отображать в % кол-во отброшенных сигмой пикселей.
2) Значение пикселя (который под курсором) отображать более детально: с плавающей точкой или ввести настройку для кол-ва отображаемых знаков после запятой. Надо для файлов в которых диапазон значений от 0 до 1. Например PixInsight любит от 0 до 1 работать.

[pavelL]

В качестве идеи:
что если отбрасываемые пиксели заменить средним по всем кадрам (или лучшим, или...), сумму кадров не делать, а перезаписать исходные только с замененными пикселями. После сложить методом drizzle (в 5ом максиме уже заработал).
Просто складывая кадры методом drizzle выходит очень шумная картинка.
PS
программу испытал в выходные - работает, немного не привычное управление, но через 2-3 попытки, все понятно. Спасибо!
Да, результат в среднем +20% S/N по сравнению с обычным сложением.

[Тот самый Кеплер]

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

[Ivan7enych]

Чтобы работать с дриззлингом надо реализовывать регистрацию по звездам, разные алгоритмы выравнивания. Я сейчас точно не возьмусь за это.

Просто складывая кадры методом drizzle выходит очень шумная картинка.

А что вы хотели? при 2х дриззлинге на каждый выходной пиксель попадает в 4 раза меньше входных пикселей.

Ещё предложение:
1) Отображать в % кол-во отброшенных сигмой пикселей.

можно

2) Значение пикселя (который под курсором) отображать более детально: с плавающей точкой или ввести настройку для кол-ва отображаемых знаков после запятой. Надо для файлов в которых диапазон значений от 0 до 1. Например PixInsight любит от 0 до 1 работать.

ок

[lazyBSD]

Чтобы работать с дриззлингом надо реализовывать регистрацию по звездам, разные алгоритмы выравнивания. Я сейчас точно не возьмусь за это.

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

[pavelL]

Просто складывая кадры методом drizzle выходит очень шумная картинка.

А что вы хотели? при 2х дриззлинге на каждый выходной пиксель попадает в 4 раза меньше входных пикселей.

Уточняю - возможностей дриззлинга в Максиме достаточно. Идея отдавать файлы в максим с замененными пикселями (те которые отбраковываються менять на усредненные). Тогда после дриззлинга картинка будет менее шумная.
И до кучи - не менять заголовок фитса если будет такая замена отбракованных пикселей.

[NKV]

Combination Tests (Sum, Average, Median, Sigma Clip and SD Mask)
http://www.dens-astropics.org.uk/page%2019.htm

Some of the combining algorithms can get rid of the normal interference in the form of satellite and aircraft trails, hot pixels and cosmic ray hits.

All were combined in different ways using Maxim DL

[onwlad]

Иван, а под Вистой твоя программа работает?
А то при попытке открыть фитсовский файл пишет - ошибка библиотеки FITSio и все закрывается  .

dll от visual studio 2005 скачана.

[Ivan7enych]

Иван, а под Вистой твоя программа работает?
А то при попытке открыть фитсовский файл пишет - ошибка библиотеки FITSio и все закрывается  .

dll от visual studio 2005 скачана.

Под Vista 64x пробовали, все работает.

Может FITs какой-то особенный?

[onwlad]

Да вроде нет.
Что из моего Сбига-2000ХМ выходит, то и пытаюсь сложить.
Да и Максим с ДСС их открывают без проблем.

[Ivan7enych]

Да вроде нет.
Что из моего Сбига-2000ХМ выходит, то и пытаюсь сложить.
Да и Максим с ДСС их открывают без проблем.

У меня дома нет ни вашего файла ни висты. Как я могу что-то проверить?

[Anton]

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

Вах!
Какой вкюсьный штючка!
Учитывая наличие дотнетовского враппера - это прямо праздник какой то.

[lazyBSD]

Посмотри на ютьюбе, что народ с ней вытворяет.