© Special Astrophysical Observatory of the Russian Academy of Sciences

Обзор

Для обработки данных FFOREST мы используем собственный программный пакет DECH. В отличие от многих других астрономических программных пакетов и систем обработки, DECH разработан для операционной системы Microsoft Windows. Однако это ограничение можно обойти с помощью эмуляторов Windows, что позволяет DECH беспрепятственно работать на платформах Linux и Mac OS. Пакет DECH был разработан по принципу "астроном создает программное обеспечение для своих коллег-астрономов", что привело к созданию дружественного и простого интерфейса. Использование командной строки сведено к минимуму и, в основном, ограничивается этапом предварительной обработки спектральных изображений, который происходит перед извлечением (щкстракцией) спектров. Извлечение и последующий анализ спектров выполняются в концепции “рабочий стол”, когда все необходимые инструменты собраны в одной/двух программах, что минимизирует необходимость взаимодействия с командной строкой.
В отличие от полностью автоматизированных программ обработки, часто называемых "pipelines", мы отказались от подхода "черного ящика", который обычно скрывает все промежуточные данные от пользователя. DECH упрощает каждый этап обработки и анализа спектральных данных, начиная от предварительной обработки изображений и извлечения спектров (включая спектры с наклонной щелью), до калибровки по длине волны с использованием двумерного полинома, нормализации континуума (как вручную, так и автоматически), измерения эквивалентных ширин и лучевых скоростей различными методами, а также кросс-корреляционного анализа и других методов. Программное обеспечение DECH предоставляет возможность выполнения высокоточных измерений лучевой скорости, необходимых для различных астрофизических задач, включая обнаружение и исследование экзопланет.

Наблюдения

Высокое качество данных наблюдений может быть утеряно из-за неадекватных или некачественных калибровочных данных. Для обеспечения оптимальных результатов даже при нестабильных погодных условиях примите во внимание несколько простых рекомендаций: Необходимое количество изображений bias зависит от желаемого соотношения сигнал/шум (S/N). Для большинства ПЗС-детекторов при стабильной работе и стандартном S/N~100-200, как правило, достаточно 10-20 изображений bias для последующего получения среднего изображения ( superbias) . Если требуется очень высокое соотношение сигнал/шум (800 - 1000+) необходимо получить 100- 150 изображений bias , чтобы свести к минимуму шум в усредненном superbias . Усредненный спектр плоского поля ( superflat ) также должен иметь соотношение сигнал/шум, по меньшей мере равное таковому в спектрах наблюдаемых объектов во всем диапазоне длин волн. Особое внимание следует уделять синей части спектра, где эффективность лабораторных источников света не очень высока. Как правило, для получения superflat с адекватным уровнем сигнала на всем диапазоне длин волн требуется минимум 10 изображений flat- field . Однако, если требуется отношение сигнал/шум ~800-1000 или выше, количество необходимых экспозиций плоского поля может превысить 100. Разделив исходные спектры (изображения) на усредненное плоское поле, можно эффективно устранить влияние неоднородности чувствительности пикселей и влияние интерференционных полос (fringes), которое особенно выражено в красной области спектра. Оптимальное время экспозиции. Как правило, для определения оптимального времени экспозиции астрономических объектов используются специальные онлайн-калькуляторы. Однако такая оценка времени экспозиции часто оказывается далекой от оптимального значения, как из-за неточных исходных данных, так и по причинам, связанным с погодными условиями. Мы рекомендуем следующий подход: в первую очередь следует сделать короткую экспозицию (одна минута или одна секунда, в зависимости от яркости объекта). Затем оценить уровень накопления сигнала, достигнутый при этой короткой экспозиции, а затем рассчитать оптимальное время экспозиции так, чтобы достичь примерно 70% полной емкости используемой ПЗС (обычно ~65000). Таким образом, максимальный накопленный сигнал не должен превышать 40000–45000 в пиковом значении. Не забудьте вычесть уровень подложки ( bias ) при оценке величины накопленного сигнала. Космические частицы. Чувствительность ПЗС-детекторов к космическим частицам не позволяет получать изображения с очень большой экспозицией Для эффективной очистки изображений от следов космических частиц рекомендуется наблюдать каждый объект не менее двух раз с одинаковой экспозицией. Это позволит более эффективно очищать спектры от следов космических частиц на этапе усреднения полученных спектров. Продолжительность каждой экспозиции не должна превышать ~60 минут. Такой подход обеспечивает надежную калибровку даже при сильном шуме детектора или изменениях погодных условий.

Предварительная обработка

Предварительная обработка спектральных изображений включает следующие продедуры: (i) Получение superbias - арифметического среднего группы изображений без засветки ( bias) . Для очистки отдельных изображений от следов космических частиц перед усреднением изображений предварительно используется медианный фильтр; (ii) Получение superflat - арифметического среднего группы изображений спектрального плоского поля ( flat-field ). Для очистки отдельных изображений от следов космических частиц перед усреднением изображений также предварительно используется медианный фильтр; (iii) Вычитание superbias из всех изображений, включая superflat. Перед выполнением процедур усреднения рекомендуется предварительного контролировать однородность исходных данных - у усредняемых изображений должен быть примерно одинаковый уровень накопления. Для этого удобно использовать поперечные (перепендикулярно направлению основной дисперсии) разрезы изображений. Важное замечание касательно учета плоского поля: дело в том, что коррекция спектров астрономических объектов с помощью спектра плоского поля может выполняться двумя способами: (1) на этапе обработки изображений, путем деления спектральных изображений объектов на изображение нормализованного плоского поля или, (2) выполнять коррекцию плоским полем уже после экстракции спектров, путем деления экстрагированных спектров объектов на экстрагированный спектр superflat . Для щелевых спектрографов следует использовать первый способ. Это связано с переменным от объекта к объекту распределением сигнала по высоте щели. Характерный пример - спектрографы UVES (oбсерватория Paranal) или MIKE (обсерватория Las Campanas), где высота щели значительна и, как результат, ширина и/или положение спектральных порядков по высоте щели меняется в зависимости от яркости объекта, качества изображений на небе, качества гидирования и др. факторов. Для оптоволоконных спектрографов (включая FFOREST ) предпочтительнее второй способ. Для учета плоского поля первым способом, то есть до экстракции спектров, необходимо предварительно вычесть рассеянный свет из всех изображений объектов и спектра superflat. Очевидно, для этого сначала определяется положение и ширина спектральных порядков во всех спектрах. После вычитания рассеянного света, изображение superflat должно быть нормализовано. Дaлее изображения объектов с вычтенным рассеянным спектром делятся на нормализованное изображение superflat . Предварительная обработка данных, полученных с помощью оптоволоконных спектрографов, в том числе FFOREST, проще и выполняется в два этапа. Сначала получают superbias , затем его вычитают из всех остальных изображений. Следующий этап - экстракция спектров.

Экстракция спектров

Процедура экстракции (извлечения) предварительно обработанных спектральных изображений FFOREST включает следующие операции: создание маски, то есть определение положения (траектории вдоль дисперсии) и границ спектральных порядков (их ширины поперек оси дисперсии); вычитание рассеянного света; собственно экстракции, которая заключается в интегрировании накопленного сигнала в направлении, перпендикулярном направлению основной дисперсии или по наклонной, если наклон щели имеет место (см. ниже), то есть по ширине спектрального порядка, в границах, определяемых маской, и по всей длине порядка. Извлечение выполняется отдельно для каждого спектрального порядка, результат сохраняется в файле формата FITS. В некоторых эшелле-спектрографах спектральные линии могут отклоняться от перпендикулярности к основной дисперсии, причем степень отклонения может изменяться как вдоль порядков, так и вдоль направления основной дисперсии. Обычная экстракция спектра из таких изображений приводит к потере спектрального разрешения и искажению профиля спектральных линий в извлеченных спектрах. В программном пакете DECH предоставлено решение для корректного извлечения спектров с наклонными спектральными линиями, в т.ч. с переменным наклоном. Более подробную информацию см. в «Руководстве по обработке данных FFOREST».

Измерения и анализ

Оптоволоконный спектрограф FFOREST оптимизирован для высокоточной спектроскопии. Обработка данных FFOREST включает стандартные для подобных инструментов процедуры. После экстракции спектра все спектры объектов делятся на экстрагированный спектр superflat . Это исправляет попиксельные вариации чувствительности и интерференционный эффект fringes . Затем строится шкала длин волн для одного или, при необходимости, нескольких (всех) спектров сравнения. Спектр ThAr-лампы обеспечивает опорные линии излучения. Перед анализом спектры нормализуются на континуум с использованием интерактивной или автоматической аппроксимации кубическим сплайном. Программный пакет DECH предоставляет широкий выбор методов спектрального анализа и способов измерения спектральных линий. Например: (i) измерение эквивалентных ширин различными методами, в т.ч. путем деления бленд; (ii) определение инструментальной или абсолютной лучевой скорости как методом интерактивной подгонки прямого и зеркального профилей так и методом кросс-корреляции; (iii) определение лучевой концентрации для оптически тонких или умеренно насыщенных линий; (iv) обработка больших массивов данных, в т.ч. автоматизированный пакетный анализ большого количества спектров, например, для обнаружения экзопланет с учетом разлучных инструментальных и астрономических поправок и др.
Руководство по обработке данных FFOREST Руководство по обработке данных FFOREST
© Special Astrophysical Observatory of the Russian Academy of Sciences

Обзор

Для обработки данных FFOREST мы используем собственный программный пакет DECH. В отличие от многих других астрономических программных пакетов и систем обработки, DECH разработан для операционной системы Microsoft Windows. Однако это ограничение можно обойти с помощью эмуляторов Windows, что позволяет DECH беспрепятственно работать на платформах Linux и Mac OS. Пакет DECH был разработан по принципу "астроном создает программное обеспечение для своих коллег-астрономов", что привело к созданию дружественного и простого интерфейса. Использование командной строки сведено к минимуму и, в основном, ограничивается этапом предварительной обработки спектральных изображений, который происходит перед извлечением (щкстракцией) спектров. Извлечение и последующий анализ спектров выполняются в концепции “рабочий стол”, когда все необходимые инструменты собраны в одной/двух программах, что минимизирует необходимость взаимодействия с командной строкой.
В отличие от полностью автоматизированных программ обработки, часто называемых "pipelines", мы отказались от подхода "черного ящика", который обычно скрывает все промежуточные данные от пользователя. DECH упрощает каждый этап обработки и анализа спектральных данных, начиная от предварительной обработки изображений и извлечения спектров (включая спектры с наклонной щелью), до калибровки по длине волны с использованием двумерного полинома, нормализации континуума (как вручную, так и автоматически), измерения эквивалентных ширин и лучевых скоростей различными методами, а также кросс-корреляционного анализа и других методов. Программное обеспечение DECH предоставляет возможность выполнения высокоточных измерений лучевой скорости, необходимых для различных астрофизических задач, включая обнаружение и исследование экзопланет.
Руководство по обработке данных FFOREST Руководство по обработке данных FFOREST

Наблюдения

Высокое качество данных наблюдений может быть утеряно из-за неадекватных или некачественных калибровочных данных. Для обеспечения оптимальных результатов даже при нестабильных погодных условиях примите во внимание несколько простых рекомендаций: Необходимое количество изображений bias зависит от желаемого соотношения сигнал/шум (S/N). Для большинства ПЗС-детекторов при стабильной работе и стандартном S/N~100-200, как правило, достаточно 10-20 изображений bias для последующего получения среднего изображения ( superbias) . Если требуется очень высокое соотношение сигнал/шум (800 - 1000+) необходимо получить 100-150 изображений bias , чтобы свести к минимуму шум в усредненном superbias . Усредненный спектр плоского поля ( superflat ) также должен иметь соотношение сигнал/шум, по меньшей мере равное таковому в спектрах наблюдаемых объектов во всем диапазоне длин волн. Особое внимание следует уделять синей части спектра, где эффективность лабораторных источников света не очень высока. Как правило, для получения superflat с адекватным уровнем сигнала на всем диапазоне длин волн требуется минимум 10 изображений flat- field . Однако, если требуется отношение сигнал/шум ~800-1000 или выше, количество необходимых экспозиций плоского поля может превысить 100. Разделив исходные спектры (изображения) на усредненное плоское поле, можно эффективно устранить влияние неоднородности чувствительности пикселей и влияние интерференционных полос (fringes), которое особенно выражено в красной области спектра. Оптимальное время экспозиции. Как правило, для определения оптимального времени экспозиции астрономических объектов используются специальные онлайн-калькуляторы. Однако такая оценка времени экспозиции часто оказывается далекой от оптимального значения, как из-за неточных исходных данных, так и по причинам, связанным с погодными условиями. Мы рекомендуем следующий подход: в первую очередь следует сделать короткую экспозицию (одна минута или одна секунда, в зависимости от яркости объекта). Затем оценить уровень накопления сигнала, достигнутый при этой короткой экспозиции, а затем рассчитать оптимальное время экспозиции так, чтобы достичь примерно 70% полной емкости используемой ПЗС (обычно ~65000). Таким образом, максимальный накопленный сигнал не должен превышать 40000–45000 в пиковом значении. Не забудьте вычесть уровень подложки ( bias ) при оценке величины накопленного сигнала. Космические частицы. Чувствительность ПЗС- детекторов к космическим частицам не позволяет получать изображения с очень большой экспозицией Для эффективной очистки изображений от следов космических частиц рекомендуется наблюдать каждый объект не менее двух раз с одинаковой экспозицией. Это позволит более эффективно очищать спектры от следов космических частиц на этапе усреднения полученных спектров. Продолжительность каждой экспозиции не должна превышать ~60 минут. Такой подход обеспечивает надежную калибровку даже при сильном шуме детектора или изменениях погодных условий.

Предварительная обработка

Предварительная обработка спектральных изображений включает следующие продедуры: (i) Получение superbias - арифметического среднего группы изображений без засветки ( bias) . Для очистки отдельных изображений от следов космических частиц перед усреднением изображений предварительно используется медианный фильтр; (ii) Получение superflat - арифметического среднего группы изображений спектрального плоского поля ( flat- field ). Для очистки отдельных изображений от следов космических частиц перед усреднением изображений также предварительно используется медианный фильтр; (iii) Вычитание superbias из всех изображений, включая superflat. Перед выполнением процедур усреднения рекомендуется предварительного контролировать однородность исходных данных - у усредняемых изображений должен быть примерно одинаковый уровень накопления. Для этого удобно использовать поперечные (перепендикулярно направлению основной дисперсии) разрезы изображений. Важное замечание касательно учета плоского поля: дело в том, что коррекция спектров астрономических объектов с помощью спектра плоского поля может выполняться двумя способами: (1) на этапе обработки изображений, путем деления спектральных изображений объектов на изображение нормализованного плоского поля или, (2) выполнять коррекцию плоским полем уже после экстракции спектров, путем деления экстрагированных спектров объектов на экстрагированный спектр superflat . Для щелевых спектрографов следует использовать первый способ. Это связано с переменным от объекта к объекту распределением сигнала по высоте щели. Характерный пример - спектрографы UVES (oбсерватория Paranal) или MIKE (обсерватория Las Campanas), где высота щели значительна и, как результат, ширина и/или положение спектральных порядков по высоте щели меняется в зависимости от яркости объекта, качества изображений на небе, качества гидирования и др. факторов. Для оптоволоконных спектрографов (включая FFOREST ) предпочтительнее второй способ. Для учета плоского поля первым способом, то есть до экстракции спектров, необходимо предварительно вычесть рассеянный свет из всех изображений объектов и спектра superflat. Очевидно, для этого сначала определяется положение и ширина спектральных порядков во всех спектрах. После вычитания рассеянного света, изображение superflat должно быть нормализовано. Дaлее изображения объектов с вычтенным рассеянным спектром делятся на нормализованное изображение superflat . Предварительная обработка данных, полученных с помощью оптоволоконных спектрографов, в том числе FFOREST, проще и выполняется в два этапа. Сначала получают superbias , затем его вычитают из всех остальных изображений. Следующий этап - экстракция спектров.

Экстракция спектров

Процедура экстракции (извлечения) предварительно обработанных спектральных изображений FFOREST включает следующие операции: создание маски, то есть определение положения (траектории вдоль дисперсии) и границ спектральных порядков (их ширины поперек оси дисперсии); вычитание рассеянного света; собственно экстракции, которая заключается в интегрировании накопленного сигнала в направлении, перпендикулярном направлению основной дисперсии или по наклонной, если наклон щели имеет место (см. ниже), то есть по ширине спектрального порядка, в границах, определяемых маской, и по всей длине порядка. Извлечение выполняется отдельно для каждого спектрального порядка, результат сохраняется в файле формата FITS. В некоторых эшелле-спектрографах спектральные линии могут отклоняться от перпендикулярности к основной дисперсии, причем степень отклонения может изменяться как вдоль порядков, так и вдоль направления основной дисперсии. Обычная экстракция спектра из таких изображений приводит к потере спектрального разрешения и искажению профиля спектральных линий в извлеченных спектрах. В программном пакете DECH предоставлено решение для корректного извлечения спектров с наклонными спектральными линиями, в т.ч. с переменным наклоном. Более подробную информацию см. в «Руководстве по обработке данных FFOREST».

Измерения и анализ

Оптоволоконный спектрограф FFOREST оптимизирован для высокоточной спектроскопии. Обработка данных FFOREST включает стандартные для подобных инструментов процедуры. После экстракции спектра все спектры объектов делятся на экстрагированный спектр superflat . Это исправляет попиксельные вариации чувствительности и интерференционный эффект fringes . Затем строится шкала длин волн для одного или, при необходимости, нескольких (всех) спектров сравнения. Спектр ThAr-лампы обеспечивает опорные линии излучения. Перед анализом спектры нормализуются на континуум с использованием интерактивной или автоматической аппроксимации кубическим сплайном. Программный пакет DECH предоставляет широкий выбор методов спектрального анализа и способов измерения спектральных линий. Например: (i) измерение эквивалентных ширин различными методами, в т.ч. путем деления бленд; (ii) определение инструментальной или абсолютной лучевой скорости как методом интерактивной подгонки прямого и зеркального профилей так и методом кросс- корреляции; (iii) определение лучевой концентрации для оптически тонких или умеренно насыщенных линий; (iv) обработка больших массивов данных, в т.ч. автоматизированный пакетный анализ большого количества спектров, например, для обнаружения экзопланет с учетом разлучных инструментальных и астрономических поправок и др.