© 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)
обработка
больших
массивов
данных,
в
т.ч.
автоматизированный
пакетный
анализ
большого
количества
спектров,
например,
для
обнаружения
экзопланет
с
учетом
разлучных
инструментальных
и
астрономических поправок и др.