Коллективная съёмка галактики М94

[Якут]

   Олег, у тебя 5555.

[mo]

Олег, у тебя 5555.

Этот косяк мы щас исправим

Сложил RGB без своего 10ч фита.
Вышло так. Так и оставлю, пожалуй.

Оба фита в архиве best-and-other-RGB а суммах от вчера http://download.milantiev.com/m94/_sum/2016-05-29/

... минут через пять дольётся.

[mo]

Приступил к сборке L-канала

Собрал other-L, то есть L-канал из худшей по FWHM половины кадров.

Собрал best-L, то есть L-канал из половины лучших кадров (по FWHM).
Т.к. там много леопарда, то собрал три варианта с разными локальными и глобальными сигмами.
SNR потерял незначительно и хотел было уже оставить самый агрессивный вариант с наименьшим лео, но на всяк не стал стирать все три.

best-L: как есть, сигма 3.8 без локальных уменьшений на лео;
best-L2: сигма 3.8 с локальным уменьшением на 20%;
best-L3: сигма 2.8 с локальным уменьшением на 40%.

SNR, соответственно, 368, 365 и ~363 попугая, то есть падение в рамках.

[mo]

Собрал other-short-L (короткие, половина с бОльшим FWHM). Там оказался заметный градиент из-за Володиных кадров (от Луны?). Что не большая проблема, т.к. в основном в этом срезе интересен центр.

Собрал best-short-L (короткие, лучшая половина с меньшим FWHM). Градиент есть и тут, но ожидаемо - меньше.

Все сборки, связанные с L-каналом в архиве L.zip (заливается, скоро будет доступен на download...).
Ещё раз, туда включены:
- other-short-L: худшая (по FWHM) половина кадров с короткой выдержкой;
- other-L: худшая (по FWHM) половина кадров с длинной выдержкой;
- best-short-L: лучшая половина кадров с короткой выдержкой;
- best-L: лучшая половина кадров с длинной выдержкой;
- best-L2: то же, что и best-L, но менее леопардово;
- best-L3: то же, что и best-L, но ещё менее леопардово, чем best-L2.

Напомню, что во всех сборках файлы other-* НЕ включают файлы best-*.
Однако, я дополнительно планирую сделать сборки с префиксом all-* с арифметической суммой other+best.

[mo]

Подготовил архив all.zip (поставил в очередь на загрузку в папку). Содержит три фита:

- all-L: сумма всех файлов L-канала с длинной выдержкой;
- all-short-L: сумма всех файлов L-канала с короткой выдержкой;
- all-RGB: сумма всех RGB (как чб камер с фильтрами, так и цветных камер), как с длинной, так и с короткой выдержками;

Суммы all-* собраны, как средне-арифметическое между соответствующими best-* и other-*.

Остался водород (H_a).

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

Суммы all-* собраны, как средне-арифметическое между соответствующими best-* и other-*.

А просто собрать сумму all в ImageIntegration не получается?
Ну или собрать сумму из всего, но без "леопардов"?

Почему собирается половина лучших по FWHM, а не 1/3, например?

Я так понимаю (по словам "локальное уменьшение сигма на лео"), что суммы складываются в FitsStacker'e?
...
Придётся мне продолжить свою работу по начальной де-леопардизации исходных снимков Никиты и подготовке сумм...

[mo]

Придётся мне продолжить свою работу по начальной де-леопардизации исходных снимков Никиты и подготовке сумм...

Ни малейших возражений. Только "за" без тени сомнений.
Буду только рад сравнить сделанным мною суммы с твоими.... когда они будут.

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

Олег, прокомментируй по техническим вопросам из предыдущего сообщения, если хочешь - в личке.
Работа трудозатратная, не хочется наступать на уже известные грабли.

[mo]

А просто собрать сумму all в ImageIntegration не получается?
Ну или собрать сумму из всего, но без "леопардов"?

Поставил складывать всё подряд из отдельных файлов, а не из сумм. Файлы all* будут заменены.
Единственное - ты верно заметил, я не пользуюсь сложением в пиксе, а складываю в фитстакере. Это не вопрос каких-либо убеждений, просто мне так удобней. Я пробовал в пиксе, но для меня не столь визуально что именно удалит сигма. Да, пробовал по кропу с рассматриванием итоговых rejection_high и _low, но мне так не удобно.

Так же актуален поднятый тобой вопрос неверного поиска SNR для разномасштабных работ. Мои кадры пикс оценивает в 20...40 раз лучше некоторых кадров Никиты. С фистакером если такая проблема и есть, то не столь явно выражена.

Почему собирается половина лучших по FWHM, а не 1/3, например?

Я хотел было поделить L на три части. И, возможно, так и стоит сделать. А вот насчёт остальных каналов и вовсе не уверен, что их стоило делить. Но пополам разделить не столь трудозатратно. Критерий "мы снимали, но половина вышла хуже, а половина лучше", мне показался вменяемым, пусть и не лучшим. Предложишь новый критерий - велкам.

В конце концов ... не мне (не только мне) обрабатывать эти фиты.

Я так понимаю (по словам "локальное уменьшение сигма на лео"), что суммы складываются в FitsStacker'e?

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

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

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

***
Мои планы пока выглядят так.

(кликните для показа/скрытия)

1. Убирание пятен с "леопардов" Никиты - в PixelMath, написал выраженьице, убирает с минимумом артефактов. Трудоёмко.

2. Исследование изменения SNR при масштабировании (уменьшении/увеличении) кадра - беру кадры Никиты, уменьшаю их, увеличиваю обратно, пытаюсь найти "гладкую функцию", описывающую, как изменяется при этом SNR при разном коэффициенте масштабирования.
Потом функцию использую при получении финальной суммы для "уменьшения" искусственно высокого значения SNR кадров с низких разрешением, таких, как твои, например, по отношению к кадрам с изначально высоким разрешением, на которых есть высокочастотные шумы, которых просто нет на снимках с более низким разрешением. - ЭТО НОВЫЙ ШАГ И НОВЫЕ ЗНАНИЯ, для меня сейчас он кажется критичным.

3. Выравниваю все кадры по опорнику (один из кадров Никиты). Собираю частичные суммы и подсуммы по авторам, на автомате, с Сигма клиппингом.
3.1. Выравниваю линейные градиенты фона.

4. Собираю финальные суммы из подсумм из п.3.1, используя коэффициенты SNR, скорректированные на значение из п.2, зависящие от величины масштабирования исходников подсуммы.
4.1. Собираю подсуммы из кадров с лучшим разрешением.

5. Делаю HDR из полученных сумм для длинных/коротких экспозиций. Всё выкладываю.
...
6. PROFIT!

Пока это то, что кажется мне оптимальным математически обоснованным алгоритмом для получения лучшей (без кавычек) суммы из разнородных исходников.
Шаг 2 - это то, чего не хватало при работе над М51. Сейчас, мне кажется, у меня есть хорошая идея, как решить проблему разнородных исходников.

Есть ещё какие-то идеи?

[mo]

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

Второй вариант суммы всегда полезен. Даже если он окажется хуже первого. Для развития навыка в "колхозной теме".
А уж если твой вариант окажется лучше, так вообще отлично! Однако, шляпу я снял входя в проект...

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

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

Исследование изменения SNR при масштабировании (уменьшении/увеличении) кадра - беру кадры Никиты, уменьшаю их, увеличиваю обратно, пытаюсь найти "гладкую функцию", описывающую, как изменяется при этом SNR при разном коэффициенте масштабирования.

Взял один из 20 мин L кадров Никиты, отмасштабировал (уменьшил) и записал его копии в масштабах 100%, 95%, 90% ... 50% от исходного размера.
Очевидно, при уменьшении кадра (если только на исходнике не было завидного oversampling'a) информация на нём теряется.

Выровнял уменьшенные копии по исходному 100% кадру, во время чего они интерполировались на лУчшее разрешение.

Посмотрел, какие относительные веса присваиваются при суммировании этим кадрам в
- PixInsight ImageIntegration;
- PixelMath SubframeSelector (скрипт);
- FitStacker (измеряемое значение SNR возведено в квадрат для получение "веса").

См. таблицу и графики.

В описанном эксперименте сильнее уменьшенным кадрам должен приписываться меньший вес при сложении, - они несут меньше информации. Кадр, уменьшенный в 2 раза от исходного (50%) должен иметь мЕньший "вес" при сложении.

Как видим (и это было неожиданно для меня) с задачей относительно корректно справляется только стандартный процесс ImageIntegration в Пиксе.

Веса кадров с худшим разрешением (бОльшим масштабом в угл.сек./пиксел), вычисляемые SuframeSelector'ом Пикса и FitStacker'ом, оказываются завышенными против ожидаемого. Эти инструменты не следует использовать при получении сумм в ситуации сложения кадров с разнородного оборудования, поскольку присваиваемые ими "веса" ведут к сумме, лучшей по SNR (из-за отсутствия высокочастотных шумов на кадрах низкого разрешения и/или с плохим FWHM), но заметно худшей по разрешению (деталям) и контрасту...

PS: немного путано, поскольку в случае "колхоза" вместо понятного критерия "оптимизации по SNR суммы" возникает новый критерий "сохранения и оптимизации "количества информации" с исходных снимков, вес кадра должен зависеть сразу от разрешения кадра, его SNR и FWHM...

[Световод]

Перекалибровал H-alpha Opossum исходники за 11 апреля 2016
Концентрические круги откалибровались.
Залил на ftp Олега

[mo]

Перекалибровал H-alpha Opossum исходники за 11 апреля 2016
Концентрические круги откалибровались.
Залил на ftp Олега

Перекинул в папку Никиты. Обновил архив opossum.zip

[mo]

Себя не похвалишь... Вырезал свою часть выступления на Астрофесте по этой теме



... правда, качество так себе. Что звука, что видео. Но... в следующий раз будет лучше, я уверен

[CoFF]

Попробовал обработать собранный материал:

[mo]

Попробовал обработать собранный материал:

Отлично!

Я не совсем про обработку, она чуть спорная на мой вкус, однако мне приятно, что пусть и не самая лучшая, но первая обработка полной суммы нашей коллективной съёмки, в виде готового жпега, наконец-то появилась. Однозначно, шаг вперёд. Спасибо.

С чистой совестью влепил "+" .

[Sanni]

Попробовал обработать собранный материал:

Нормально так вроде, но со звёздами не порядок на мой взгляд! и тоже плюс

[Opossum]

Попробовал обработать собранный материал:

кмк, очень ушло в синеву.
Про звёзды уже отметили

[Opossum]

призываю соблюдать правила ориентации - север сверху