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

[Ivan7enych]

Уф, дописал таки сложение sd mask.
Писал сегодня в 7 утра в самолете, подлетая к москве...

Тренировался складывать кадры выровненные по комете так, чтобы по максимуму исчезли все звезды.

2 результата для сравнения.
1-й сложен с сигмой 2.5
2-й сложен в 2 дополнительных прохода sd mask с той же сигмой.

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

Конечно, сложение кадров от кенона 450d, где каждый кадр = 140МБ (цветной float fits) занимает весьма продолжительное время...

Версия еще не сильно проверенная.
qhy.narod.ru/FTS/FITStackerV8.zip

Да, учел несколько пожеланий.

Список фитсов показывает имя фильтра, дает по нему сортировать, дает выделить кучу кадров и кнопкой del убрать из списка.

Добавил баланс белого (для просмотра)

[Ruslan 7]

Как складывать по ядру кометы? Выделяю зелёным квадратом фон, красным ядро кометы. Складывает по звёздам.

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

выровняйте кадры в максиме по ядру

[Ruslan 7]

Ясно. я думал прога сама выравнивает.

[lazyBSD]

Пока, надеюсь, не умеет.

Вань, смайл адресован тебе.

[Грин]

В фильме он сказал, что не будет этим заниматься. Алигните в Максиме, благо он делает это прекрасно. Да и по сложению у меня пока не получилось получить S\N в фитстакере лучше, чем в максиме, правда, на неплохих исходниках...

[martalex]

CCDstack (http://ccdware.com/products/ccdstack/) никто не смотрел?

Беглое ознакомление показало:
1) способность проги оперировать большими объемами данных
2) сильно наворочанные алгоритмы отброса данных - сигма клиппинг для маньяков, причем вроде бы с таким же анализом бэкграунда, как у Ивана.
3) ну и вообще настроечек там много, тянет на серьезный софт.

очень клевая программа... пользовался ей, пока 30 дней не истекли...
функций действительно куча... и очень классных...

собираюсь ее купить...

P.S. А Ивану спасибо за программку 

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

Иван, если позволишь (или поясни, в чём я заблуждаюсь или что делаю не так, плиз):
1. выставлять значения Local Sigma Scale для трёх цветовых каналов раздельно для каждого снимка вручную мышкой неудобно и трудоёмко.
Можно добавить некое средство, которое как минимум для выбранного кадра скопирует значение Local Sigma Scale из настраиваемого цветового канала в 2 других? См. эксиз внизу, возможно, впрочем, что кнопка - не лучшее решение.
2. Также было бы неплохо добавить возможность ввода значения Local Sigma Scale с клавиатуры, а не только мышкой, как сделано сейчас.
3. Не уверен, как вычисляется значение median, от которого по сигма-допуску фильтуются значения в точках в индивидуальных кадрах, но - не имеет ли смысл добавить флажок "Skip for Master Median" или что-то в этом роде, который первоначально бы исключал конкретный кадр из вычисления "опорного", но затем анализировал бы значения его точек и добавлял бы их в финальное значение с необходимыми коэффициентами? Может пригодиться, когда есть несколько "сильно бракованных кадров" с тянутыми звёздами, подвижками и т.д., для того, чтобы они не влияли на структуру изображения, но добавили бы свои "5 копеек" в сигнал.

Спасибо!

[Ivan7enych]

Иван, если позволишь (или поясни, в чём я заблуждаюсь или что делаю не так, плиз):
1. выставлять значения Local Sigma Scale для трёх цветовых каналов раздельно для каждого снимка вручную мышкой неудобно и трудоёмко.

Плодить кнопки не хочется, их и так много.
Я могу объединить Local Sigma Scale для трёх цветов намертво. По моему опыту, настраивать эту штуку приходится очень редко, а по разному для каждого цвета - еще реже.

2. Также было бы неплохо добавить возможность ввода значения Local Sigma Scale с клавиатуры, а не только мышкой, как сделано сейчас.

можно

3. Не уверен, как вычисляется значение median, от которого по сигма-допуску фильтуются значения в точках в индивидуальных кадрах

Это просто среднее всех кадров, никакой медианы там нет.

не имеет ли смысл добавить флажок "Skip for Master Median" или что-то в этом роде, который первоначально бы исключал конкретный кадр из вычисления "опорного", но затем анализировал бы значения его точек и добавлял бы их в финальное значение с необходимыми коэффициентами? Может пригодиться, когда есть несколько "сильно бракованных кадров" с тянутыми звёздами, подвижками и т.д., для того, чтобы они не влияли на структуру изображения, но добавили бы свои "5 копеек" в сигнал.

Вообще-то для этого сделан многопроходный сигма-клиппинг (sd-mask), он полностью выполняет ту же задачу - отбрасывает всё сильно бракованное из суммы, и берет результат в качестве опорного для следующей итерации сложения.

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

1. выставлять значения Local Sigma Scale для трёх цветовых каналов раздельно для каждого снимка вручную мышкой неудобно и трудоёмко.

Я могу объединить Local Sigma Scale для трёх цветов намертво. По моему опыту, настраивать эту штуку приходится очень редко, а по разному для каждого цвета - еще реже.

Ок, объединение навсегда - то, что нужно:)

2. Также было бы неплохо добавить возможность ввода значения Local Sigma Scale с клавиатуры, а не только мышкой, как сделано сейчас.

можно

Спасибо!

Вообще-то для этого сделан многопроходный сигма-клиппинг (sd-mask), он полностью выполняет ту же задачу - отбрасывает всё сильно бракованное из суммы, и берет результат в качестве опорного для следующей итерации сложения.

Ок, не пробовал ещё последнюю версию, работал с v.7. Посмотрю, как этот самый sd-mask работает.

FITStacker падает, если по ошибке в него сбросить файлы в формате, отличном от FITS. Я его случайно убил с помощью PIC от IRIS.

[Pilgrim]

Иван, во-первых спасибо за прогу ещё раз, во всю использую! 
Особенно сейчас - накопил тонны разнородного лунно-городского материала.

Маленькая просьба - а можно в списке по кнопке "list..." показывать коэффициенты сложения для каждого файла?
Просто для любопытства - знать "по какому курсу пошла вчерашняя ночь". Или например принять решение выкинуть файл целиком (например сильный смаз заставляет задирать сигму для всей серии, а толку от кадра около нуля - удаляем)

[Ivan7enych]

Маленькая просьба - а можно в списке по кнопке "list..." показывать коэффициенты сложения для каждого файла?
Просто для любопытства - знать "по какому курсу пошла вчерашняя ночь". Или например принять решение выкинуть файл целиком (например сильный смаз заставляет задирать сигму для всей серии, а толку от кадра около нуля - удаляем)

Этот коэфициент по сути уже показывается. Выделив 2 прямоугольника и нажав List ты получаешь список с колонкой S/N (сигнал/шум), по ней можно отсортировать все файлы и выкинуть нижние.

Коэфициент при сложении = S/N в квадрате (деленный на сумму коэфициентов для всех файлов).

[Pilgrim]

Маленькая просьба - а можно в списке по кнопке "list..." показывать коэффициенты сложения для каждого файла?
Просто для любопытства - знать "по какому курсу пошла вчерашняя ночь". Или например принять решение выкинуть файл целиком (например сильный смаз заставляет задирать сигму для всей серии, а толку от кадра около нуля - удаляем)

Этот коэфициент по сути уже показывается. Выделив 2 прямоугольника и нажав List ты получаешь список с колонкой S/N (сигнал/шум), по ней можно отсортировать все файлы и выкинуть нижние.

Коэфициент при сложении = S/N в квадрате (деленный на сумму коэфициентов для всех файлов).

Сортировка, это качественное, а хочется количественно знать, сколько времени ушло коту под хвост при крайней фотосессии
Вот этот "S/N в квадрате (деленный на сумму коэфициентов для всех файлов)" отдельной колонкой и показать, было бы крайне удобно, а то утром после наблюдений думать умом трудно 

[Ivan7enych]

Сортировка, это качественное, а хочется количественно знать, сколько времени ушло коту под хвост при крайней фотосессии

Ладно, сделаю.

[Севастьянов Игорь (RotT_)]

...
Я правильно понимаю технологию:
Колибрую все снимки в IRISe, перегоняю в RGB, Регистрирую (совмещаю), подрезаю ... а уже потом из CFA (*.pic) преобразовываю в ФИТСы ? (вот как это сделать ?)
И только после этого запускаю в FITStacker ?

[Pilgrim]

...
Я правильно понимаю технологию:
Колибрую все снимки в IRISe, перегоняю в RGB, Регистрирую (совмещаю), подрезаю ... а уже потом из CFA (*.pic) преобразовываю в ФИТСы ? (вот как это сделать ?)
И только после этого запускаю в FITStacker ?

1) Калибруешь в Ирис.
2) Переводишь из CFA в RGB (дебаеризация, Digital Photo->Sequence CFA conversion)
3) Выравниваешь
4) Переключаешь Iris в режим FITS (File->Settings или Ctrl+R, затем чекбокс "FIT", Ок)
5) Конвертируешь последовательность PICов в последовательность фитсов (File->Select files..., чекбокс PIC, имя последовательности, перетаскиваешь из любого файл-менеджера ПИКи в окошко)

Затем FITs Stacker.

6) На выходе сумма в виде цветного FITsа с плавающей точкой. Ирис такие открывает только как ЧБ, поэтому я обычно открываю Максимом и пересохраняю в целочисленный ФИТС или сразу в 16-битный линейный тиф и на пост-обработку в ФШ.

З.Ы. В связи с этим просьба, если у кого-то получается более коротким путём - поделитесь!

З.З.Ы Иван, можно ли сделать опцию сохранения в целочисленный ФИТс, с нормированием с диапазона в отрезок [0...32767] по максимальному значению в файле, специально для заскоков Ириса? Или я проглядел и она уже есть?

[Ivan7enych]

5) Конвертируешь последовательность PICов в последовательность фитсов (File->Select files..., чекбокс PIC, имя последовательности, перетаскиваешь из любого файл-менеджера ПИКи в окошко)

А сразу работать в FITS неудобно? Или промежуточные слишком много места занимать начинают?

З.З.Ы Иван, можно ли сделать опцию сохранения в целочисленный ФИТс, с нормированием с диапазона в отрезок [0...32767] по максимальному значению в файле, специально для заскоков Ириса? Или я проглядел и она уже есть?

Так есть же Ctrl + Save (-32768...32767) - сохраннение для ириса и Shift+ Save (0...65535) - для максима.

[Pilgrim]

5) Конвертируешь последовательность PICов в последовательность фитсов (File->Select files..., чекбокс PIC, имя последовательности, перетаскиваешь из любого файл-менеджера ПИКи в окошко)

А сразу работать в FITS неудобно? Или промежуточные слишком много места занимать начинают?

Дело в том, что ирис работает с цифрозеркальными баеризованными картинками только в PIC. Упс, сейчас проверил в крайней версии - уже не только 
Возможно стоит попробовать работать исключительно в фитсах. Просто полноценная работа с цветными фитсами в Ирис появились относительно недавно, и есть глюки - например зеркальное отображение по вертикали, знаковый диапазон (-32768...32767) вместо Максимовского беззнакового 0..65535 без возможности конвертации при импорте.  Вобщем мне субьективно кажется, что могут вылезти ещё подводные камни, поэтому по-старинке, предпочитаю "родной" ирисовский pic.

З.З.Ы Иван, можно ли сделать опцию сохранения в целочисленный ФИТс, с нормированием с диапазона в отрезок [0...32767] по максимальному значению в файле, специально для заскоков Ириса? Или я проглядел и она уже есть?

Так есть же Ctrl + Save (-32768...32767) - сохраннение для ириса и Shift+ Save (0...65535) - для максима.

О! Спасибо! Помню, что "про что-то такое с Control-shift в ветке было", но тогда не оценил полезности и в памяти не отфиксировал. Уж очень мудрёная придумка с клавишами, в жизни не догадаешься.  ИМХУ стоит спрашивать при сохранении, или лучше радио-баттоны сделать в диалоге возле кнопки сэйв. Прости за мелкий интерфейсный придиразм, но многим будущим пользователям будет полезно. 

[Ivan7enych]

ИМХУ стоит спрашивать при сохранении, или лучше радио-баттоны сделать в диалоге возле кнопки сэйв. Прости за мелкий интерфейсный придиразм, но многим будущим пользователям будет полезно. 

Интерфейс - это да, по хорошему переделывать надо. Програмка из мелкой постепенно разрослась и облепилась фичами...

[Севастьянов Игорь (RotT_)]

1) Калибруешь в Ирис.
2) Переводишь из CFA в RGB (дебаеризация, Digital Photo->Sequence CFA conversion)
3) Выравниваешь
4) Переключаешь Iris в режим FITS (File->Settings или Ctrl+R, затем чекбокс "FIT", Ок)
5) Конвертируешь последовательность PICов в последовательность фитсов (File->Select files..., чекбокс PIC, имя последовательности, перетаскиваешь из любого файл-менеджера ПИКи в окошко)

Затем FITs Stacker.
......

Огромное спасибо за подробное описание!
однако ...
Выполнил п.5 получил последовательность файлов *n.fit
Делаю add в FITstacker-e получаю ошибку:
FITSIO Library Error
Какую библиотеку надо? И куда ?