АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ, 2016, том 71, № 4, с. 532-542
УДК 52-17
АВТОМАТИЗИРОВАННАЯ СИСТЕМА ОБРАБОТКИ НАБЛЮДАТЕЛЬНЫХ ДАННЫХ НА РАДИОТЕЛЕСКОПЕ РАТАН-600
©2016 Р. Ю. Удовицкий1*, Ю. В. Сотникова1**, М. Г. Мингалиев1-2, П. Г. Цыбулев1, Г. В. Жеканис1, Н. А. Нижельский1
1Специальная астрофизическая обсерватория РАН, Нижний Архыз, 369167 Россия 2Казанский (Приволжский) федеральный университет, Казань, 420008 Россия Поступила в редакцию 23 марта 2016 года; принята в печать 13 сентября 2016 года
В работе представлена разработанная на РАТАН-600 и введенная в штатную эксплуатацию автоматизированная система оценки параметров радиоисточников, наблюдаемых на всех имеющихся радиометрах континуума (два приемных комплекса вторичных зеркал №1 и №2 с общим числом радиометров — 30). Система также применяется для мониторинга параметров антенных и приемных систем радиотелескопа РАТАН-600, который проводится с использованием текущих измерений калибровочных радиоисточников.
Ключевые слова: радиоконтинуум: общее — методы: наблюдательные — техники: радиоастрономические
1. ВВЕДЕНИЕ
Радиотелескоп РАТАН-600 — крупнейший в мире радиотелескоп с антенной переменного профиля [1, 2]. Особенности его геометрии позволяют за счет поперечных аберраций при малых потерях сигнала получать радиоспектр источника в течение одного наблюдения (1—2 минуты), что важно при исследованиях переменных источников. Современный уровень приемной аппаратуры обеспечивается сверхмалошумящими усилителями на транзисторах с высокой подвижностью электронов (HEMT) и цифровыми сигнальными процессорами в системе регистрации данных. Наблюдения проводятся одновременно в двух- и трехзеркальной конфигурациях радиотелескопа (на двух независимых секторах антенны, Северном и Южном), перекрывая широкий диапазон склонений исследуемых объектов — от +80° до —40°. При этом возможно получать наблюдательный материал для больших списков объектов (до 200 в сутки), оценивать и отбирать наиболее перспективные для исследований. На РАТАН-600 регулярно проводится мониторинг активных ядер галактик (АЯГ). В результате многолетних наблюдений накоплен многочастотный материал, например, [3—10], охватывающий большой список представителей АЯГ на длительных временных масштабах.
E-mail: nemiroff@sao.ru
E-mail: sjv@sao.ru
На протяжении 40 лет работы радиотелескопа уже создавались программы обработки наблюдательных данных континуума для упрощения процедуры оценок плотностей потоков. Во всех из них используются общепринятые методы обработки цифрового сигнала с разной степенью автоматизации и проработки пользовательского интерфейса. Среди них можно выделить:
• алгоритм потоковой многочастотной обработки, разработанный Бурсовым Н.Н. [11], в котором используются модули штатного пакета FADPS [12] в среде Linux; применяется для данных широкополосных радиометров континуума вторичного зеркала № 1, частично автоматизирован (до этапа получения антенных температур);
• пакет YURZUF для автоматизированного режима обработки мгновенных спектров радиоизлучения дискретных объектов [13], включающий в себя весь набор основных программ и командных файлов, необходимых для обработки сигнала: от расчета используемых параметров телескопа до формирования таблицы с окончательными результатами в виде спектров и плотностей потоков;
• полуавтоматический алгоритм потоковой обработки, созданный группой сотрудников ГАИШ МГУ для обработки обзоров радиоисточников на РАТАН-600 (Горшков А. Г., Конникова В. К.).
Актуальность решаемой нами задачи определяется прежде всего большим объемом данных,
Web browser -based GUI
Apache server
Data reduction server
System core
(PHP)
MySQL server
stages of reduction
Database
Рис. 1. Структура программного обеспечения системы наблюдательных данных РАТАН-600.
который требуется обработать для получения спектров исследуемых объектов. Так, при исследованиях АЯГ непрерывно используются два рабочих приемных комплекса вторичных зеркал №1 и №2 (девять радиометров). При среднем количестве наблюдений, равном 150 в сутки, суммарное число записей в месяц составляет 40.5 тысяч. Ручная обработка такого объема данных не представляется возможной. Кроме того, необходимо оперативно (ежедневно) получать информацию о спектральном поведении большого числа радиоисточников для быстрой оценки их параметров и текущего состояния систем радиотелескопа.
В работе описана общая характеристика системы и структура программного обеспечения, алгоритмы получения антенных температур и спектральных плотностей потоков, приведены соотношения для ошибок измерений, представлена реализация контроля параметров антенных и приемных систем.
2. ОБЩАЯ ХАРАКТЕРИСТИКА СИСТЕМЫ ОБРАБОТКИ
Система представляет собой комплекс программного обеспечения (ПО), объединяющий программные компоненты штатного пакета обработки наблюдений РАТАН-600 в среде Linux — FADPS [12, 14], а также программное обеспечение, разработанное авторами. Алгоритмы обработки данных реализованы на языках С и С++. В системе использована СУБД MySQL для универсальности доступа ко всем данным, как исходным, так и появляющимся в процессе обработки наблюдений, а также для хранения конфигурации системы, информации о калибровочных объектах, характеристик радиотелескопа и радиометров, хранения полученных результатов. Связь системы обработки с СУБД организована с использованием языка PHP. Для взаимодействия с системой пользователю предоставлен разработанный нами веб-интерфейс.
В процессе обработки наблюдатель задает входной массив данных, а также контролирует все
последующие этапы обработки: оценку качества наблюдательного материала, оценки параметров калибровочных радиоисточников, получение калибровочных характеристик радиотелескопа и оценки параметров исследуемых объектов. Реализована возможность активно вмешиваться в ход обработки на любом этапе при возникновении проблем, связанных, например, с качеством данных или изменившимися параметрами радиотелескопа. Конечный результат в виде спектров, кривых блеска и значений плотностей потоков может сохраняться в нескольких текстовых и графических форматах, а также экспортироваться и автоматически загружаться в интерактивный каталог BLcat [15]. Анализ и обработка наблюдательного сета длительностью в один месяц до получения окончательных результатов занимает около 20 минут.
3. СТРУКТУРА ПРОГРАММНОГО ОБЕСПЕЧЕНИЯ СИСТЕМЫ
Структура ПО системы обработки схематично показана на рис. 1. Для реализации была выбрана архитектура «клиент—сервер». При этом вся обработка инициируется пользователем на клиентской стороне через веб-интерфейс и выполняется на стороне сервера. Таким образом, специализированное ПО для клиентского компьютера полностью отсутствует (необходим только веб-браузер). Такая архитектура позволяет реализовать централизованное сопровождение ПО на сервере. В случае необходимости в роли и клиента, и сервера может выступать один и тот же компьютер.
Основным ядром сервера системы («System core» на рис. 1) является набор PHP-модулей, которые активируются пользователем через веб-интерфейс и выполняют все задачи, связанные с процедурой обработки: запуск внешних программ обработки сигналов, взаимодействие с пользователем через Apache-сервер (визуализация процесса обработки и представления результатов), взаимодействие с базой данных системы через MySQL-сервер.
Рис. 2. Структурная схема системы обработки.
Встроенная в систему база данных содержит следующие таблицы:
• «Входной набор данных» — список радиоисточников, подлежащих обработке в текущем сеансе; хранит всю информацию, полученную из FITS-заголовков входных файлов; ключевым полем таблицы являются уникальные имена источников, определяемые координатами и датой наблюдений;
• «Известные радиоисточники» — физические параметры радиоастрономических объектов базы данных системы: координаты, альтернативное имя, красное смещение, тип объекта и другая информация из литературы; система идентифицирует объекты «входного набора данных» с объектами этой таблицы по координатам с окном 30 х 30 секунд по прямому восхождению и склонению;
• «Параметры калибровочных объектов»;
• «Параметры антенны и радиометров»;
• «Служебные константы,» используемые при реализации алгоритмов обработки: значения рабочих окон записей, окон вычитания линейного и нелинейного фона и др.
Собственно обработка наблюдательных данных производится ядром системы в блоке «Stages of reduction» (рис. 1). Все этапы обработки более подробно описаны в последующих разделах.
4. ПОСЛЕДОВАТЕЛЬНОСТЬ ЭТАПОВ ОБРАБОТКИ ДАННЫХ
4.1. Определение входного набора данных
Структурная схема этапов обработки наблюдений показана на рис. 2. Номера (1)—(4) на рисунке соответствуют последовательности стандартной процедуры автоматической обработки. Пер-
вым этапом подготовки к обработке является создание входного набора данных — списка объектов на заданном интервале дат и времени (этап «Define the data set» на рис. 2) с учетом конфигурации антенны (заданный сектор радиотелескопа и номер вторичного зеркала). Кроме того, при создании этого массива, как правило, учитывается имя автора наблюдательной программы, а также требуемый частотный диапазон наблюдений.
При автоматизации процесса обработки наблюдений первый этап является самым трудоемким, поскольку одна из задач системы — обработка любых данных цифрового архива РАТАН-600, накопленного за более чем 30 лет наблюдений. На таком длительном интервале времени архив не является однородным, и требуется учитывать так называемые «исторические данные»:
• развитие регистрационного (архивного) формата данных: добавление новых необходимых параметров в информационных заголовках архивных файлов;
• несколько этапов модернизации систем регистрации данных, которые отражают изменения как в идеологии регистрации, так и в архитектуре самих радиометров континуума;
• различия в координатно-временной привязке наблюдений.
Поэтому первый этап включает в себя не только формирование заданного списка объектов, но и предварительную подготовку наблюдательного материала к обработке, в результате которой формируется однородный по параметрам и структуре массив. Остановимся подробнее на перечисленных выше факторах.
4.1.1. Форматы! данных радиометров континуума
В начале 90-х годов основным выходным форматом систем регистрации данных комплексов радио-
метров континуума стал RFLEX-формат, по своей структуре близкий к FITS-формату. Аббревиатура RFLEX означает RATAN-FLEX, где FLEX (FLexible EXchange) — формат гибкого обмена данными [14]. RFLEX является «многочастотным» (или многоканальным) форматом, когда в одном файле хранятся данные, полученные на разных радиометрах (всех или некоторых) конкретного радиометрического комплекса. Наблюдательные записи предыдущих лет (80-х и начала 90-х годов), хранимые в т.н. R-формате, были преобразованы в RFLEX-формат В. К. Кононовым. RFLEX-формат разрабатывался для хранения и передачи данных, но не предназначался для их обработки.
В процессе обработки сигналов, полученных от всех радиометров, требуется последовательная работа с одномерными временными рядами, которые соответствуют сигналам с выходов детекторов каждого радиометра. Для обработки данных был введен так называемый F-формат [16], который является подмножеством FITS-формата данных и представляет собой одномерный массив (временной ряд), соответствующий одному наблюдению на одной длине волны одним радиометром (в отличие от «многочастотного» RFLEX-формата данных). В системе FADPS существует программа перекодировки из RFLEX-формата в F-формат данных (программа fl2f). Параметрически текущая версия F-формата полностью согласована с соответствующей версией RFLEX-формата.
4.1.2. Изменения в регистрации данных и в радиометрах
Изменения в схемах построения радиометров привели к изменению систем регистрации, в результате чего архивный RFLEX-формат с 2010 г. хранит сигналы радиометров в измененной форме, отличной от более ранних способов представления. Если ранее одному выходу радиометра (одному детектору) соответствовал один одномерный временной ряд, то сейчас одному детектору соответствует два одномерных временных ряда (сигнала). Каждый из этих двух сигналов соответствует одному из полупериодов модуляции радиометра, условно названных сигналами «a» (antenna) и «e» (equivalent), как описано в работе [17]. Для построения окончательного сигнала радиометра требуется выполнить его синтез, который заключается в получении либо разностного сигнала a — e (если радиометр модуляционный), либо сигнала (a + e)/2 (если это «радиометр полной мощности»). Информация об архитектуре радиометра хранится в RFLEX-подзаголовке соответствующего сигнала.
4.1.3. Время и координаты
Эфемериды наблюдаемых радиоисточников исторически рассчитывались разными авторами наблюдательных программ на различные астрономические эпохи. Для работы в единой системе координат нами принято решение автоматически привести все данные к одной эпохе (нами выбрана эпоха Л2000) путем пересчета координат источника штатной программой пакета FADPS.
4.1.4. Подготовка исходных данных к обработке
Исходя из перечисленных выше факторов, связанных с первым этапом обработки, автоматизированные нами выбор и подготовка данных к обработке выглядят следующим образом:
• выборка данных из цифрового архива радиометров континуума РАТАН-600 согласно списку критериев: интервал дат—времени, автор наблюдательной программы, конфигурация антенны и приемных комплексов;
• приведение всех данных к последней версии RFLEX-формата с заполнением всех (насколько это возможно) новых полей в информационных заголовках RFLEX-файлов;
• распаковка многочастотных архивных данных: получение из каждого RFLEX-файла, содержащего N одномерных сигналов, N файлов F-формата, каждый из которых содержит один одномерный сигнал;
• синтез выходных сигналов радиометров согласно архитектуре последних;
• пересчет координат и привязка временной шкалы на эпоху Л2000 для всех объектов.
В результате создания исходного набора данных в базе данных системы автоматически заполняются таблицы, где каждая строка содержит все параметры из заголовков F-файлов (описывает полностью наблюдение на данной длине волны данным радиометром), а также имена директорий на диске компьютера, где физически хранятся исходные данные и будут храниться промежуточные и окончательные результаты обработки. Также на этом этапе подготовки к обработке в базе данных системы заполняются таблицы, относящиеся к наблюдаемым радиоисточникам, и автоматически производится идентификация радиоисточников по координатам с таблицей «Известные радиоисточники» (см. выше). Кроме того, формируется список калибровочных радиоисточников, которые наблюдались в программе. На этом этапе есть возможность ручной коррекции списка исследуемых объектов. Массив готов к следующему этапу автоматизированной обработки.
4.2. Измерение параметров радиоисточников
Исходный набор данных содержит два логических подмножества: калибровочные радиоисточники и исследуемые объекты. Задача дальнейшей обработки сводится к получению одних и тех же параметров для этих двух подмножеств объектов: амплитуды и полуширины гауссова отклика неподвижного радиотелескопа на прохождение радиоисточника через диаграмму направленности (ДН). Массивы этих двух подмножеств поочередно направляются на одну и ту же автоматическую процедуру оценки указанных параметров, показанную на рис. 2 в виде блока «Data reduction sequence.» На этом этапе производится ряд стандартных для РАТАН-600 процедур обработки наблюдений: выделение участка записи, содержащего радиоисточник (далее работа производится только с ним), вычитание на выбранном участке линейного тренда, а также фоновой кривой более высокого порядка (больше 2).
Как упоминалось выше, отклик радиотелескопа на прохождение компактного радиоисточника через неподвижную ДН является гауссовым, поэтому следующий шаг автоматической обработки — это аппроксимация зашумленной кривой прохождения радиоисточника гауссианой методом наименьших квадратов. Однако исходные записи содержат белую шумовую компоненту в диапазоне частот 0 — Fmax. Для систем регистрации радиометров континуума РАТАН-600 [17] частота Fmax примерно равна 9 Гц. В то же время Фурье-спектр гауссова отклика на прохождение источника через неподвижную ДН локализован в диапазоне 0—1.5 Гц (ширина этого спектра обратно пропорциональна длине волны). При обработке необходимо устранить избыточный шум для максимизации отношения сигнал/шум, что достигается путем применения фильтра нижних частот (ФНЧ). В нашей системе обработки в качестве такого ФНЧ применен разработанный нами цифровой фильтр с бесконечной импульсной характеристикой на основе фильтра Баттерворта с нулевой фазовой характеристикой [17]. Удобство применения данного фильтра состоит в том, что его входными параметрами являются частота среза и порядок. При таком подходе требуемые параметры фильтра синтезируются «на лету,» в начале его работы. Это позволяет менять частоту среза фильтра, поскольку она должна быть согласована с максимальной частотой в спектре радиоисточника при наблюдении на заданной высоте и на данной длине волны. Порядок фильтра выбран нами равным 8, что обеспечивает достаточную крутизну характеристики. Свойство нулевой фазовой характеристики ФНЧ обеспечивается путем двухкратной фильтрации сигнала фильтром вдвое меньшего порядка, причем перед
вторым шагом фильтрации производится «обращение времени» массива данных. Нулевая фазовая характеристика ФНЧ приводит к тому, что отклик от радиоисточника в результате фильтрации не искажается и не смещается по оси времени.
После НЧ-фильтрации производится аппроксимация отклика на источник гауссианой. В результате для всех объектов получены искомые оценки параметров: амплитуды (антенной температуры на частоте v) и полуширины гауссова отклика на данной длине волны фх,и (здесь «i» означает investigated source, а индекс x говорит о том, что сечение источника диаграмой направленности производится вдоль горизонтальной координаты x). Все полученные на данном этапе обработки оценки параметров источников заносятся в базу данных системы для использования их на следующем этапе обработки — калибровке измерений по спектральной плотности потока.
Программное обеспечение, разработанное авторами для данного этапа обработки, автоматизирует процесс вызова необходимых методов и подготовки входных данных для них, а также хранение всех промежуточных и выходных результатов. Кроме того, если пользователю требуется визуально контролировать процесс обработки, графический интерфейс отображает всю необходимую информацию: таблицы с оценками параметров радиоисточников, исходные и обработанные сигналы радиоисточников.
Все методы в блоке «Data reduction sequence,» за исключением НЧ-фильтрации, реализованы в системе FADPS и протестированы наблюдателями РАТАН-600 в течение более 15 лет. Разработанный нами НЧ-фильтр уже более пяти лет применяется также в составе ПО-систем регистрации данных РАТАН-600 [17].
4.3. Оценки спектральных плотностей потоков.
Ошибки измерений
Дальнейший этап — калибровка амплитуды исследуемых объектов по плотности потока излучения на основании полученных оценок для калибровочных радиоисточников. Стандартная процедура оценки спектральной плотности потока описана подробно в работах [3, 9, 18, 19]. Список калибровочных объектов, значения плотностей потоков, поправки на линейную поляризацию и угловой размер хранятся в базе данных системы (см. таблицу 1). Здесь: первая колонка — имя источника, (2) — имя по координатам RA и DEC (эпоха 2000.0), колонки (3)—(10) — принятые значения плотностей потоков на соответствующих частотах. В работе используется абсолютная шкала потоков [20—22], ее ошибки составляют
Характеристики калибровочных объектов
Источник Имя Частота, ГГц
(RA, DEC) 21.7 11.2 8.2 7.7 4.8 2.3 1.3 1.1
ЗС 48 013741 + 330940 1.21 2.43 3.32 3.57 5.62 10.77 17.13 22.33
J 0237- 23 024008 - 230916 1.00 1.18 1.55 1.64 2.66 4.74 5.78 6.41
4С+76.03 041045 + 765645 1.28 1.79 2.17 2.26 2.97 4.32 5.83 6.36
ЗС 138 052109+ 163822 1.07 1.93 2.45 2.62 3.75 6.15 8.59 10.23
ЗС 147 054236 + 495106 1.81 3.60 4.91 5.21 8.06 15.08 22.85 28.96
ЗС 161 062709 - 055311 1.25 2.79 3.94 4.21 6.79 12.97 19.88 24.10
J 1154-35 115421 - 350529 0.85 1.49 1.91 2.01 2.93 4.59 6.32 7.76
ЗС 286 133108 + 303033 2.61 4.23 5.32 5.52 7.52 11.50 15.56 17.20
ЗС 286 new 133108 + 303033 2.56 4.15 5.18 5.41 7.51 11.58 15.37 16.70
4С+12.50 134733+ 121739 1.12 1.76 2.04 2.12 2.95 4.13 5.26 6.02
ЗС 295 new 141120 + 521209 0.99 2.39 3.51 3.78 6.51 14.29 23.90 27.18
ЗС 309.1 145907+714020 1.01 1.78 2.29 2.41 3.47 5.85 8.49 9.41
NGC 7027 210701 + 421410 5.71 6.10 5.95 5.93 5.05 3.00 1.15 0.81
10% на 21.7 ГГц и 3—5% на остальных частотах РАТАН-600. Данные о линейной поляризации взяты из работы [23].
Для РАТАН-600 эффективная площадь Аед-) „(Ь) (на частоте V) зависит от высоты наблюдаемого источника [24]. Наблюдаемая спектральная плотность потока определяется из известного соотношения [25]:
2Ш
Sv,vis -
^eff,v (h)
где
gv(h) -
2k
^eff,v (h)'
Здесь Ta, v — антенная температура, k — постоянная Больцмана, gv(h) — коэффициент, отражающий изменение эффективной площади антенны с высотой h. Функция gv(h) (мы называем ее «калибровочной кривой») строится для набора высот, определяемых высотами используемых калибровочных источников и аппроксимируется полиномом второй степени (назовем эту функцию gv(h)) для всех возможных высот наблюдения (см. рис. 3).
При построении калибровочных кривых видимый поток Sv, vis от калибровочного радиоисточника пересчитывается к истинному радиопотоку Sv согласно формуле:
Sv = gi g2 Sv, vis,
где коэффициент д1 учитывает угловые размеры, а коэффициент д2 — поляризацию калибровочных источников.
Значения коэффициента д1 взяты нами из работы [26]. Исключение составляет планетарная туманность NGC7027, для которой считается [18], что ее изображение есть равномерный диск с диаметром Б = 8". Коэффициент для учета размеров источника определяется по формуле:
gi -
D
1 — exp(-x2)
1.2фРх
где фрх — горизонтальный угловой размер ДН.
Учет поляризации калибровочного источника производится согласно известной формуле [25]:
g2 -
1 — Р . 2 (Ф
- р cos -
Здесь ф — угол между плоскостью поляризации источника и антенны на сфере Пуанкаре. Величина д2 является коэффициентом связи волны от источника с антенной. Для источников, поляризованных только линейно, ф = 2(Аа — А;), где Аа и А; — позиционные углы поляризации источника и антенны.
В настоящее время на вторичном зеркале №1 на длинах волн 1.38, 2.7 и 3.7 см анализаторы линейной поляризации стоят на прием вертикальной поляризации (вектор ЕЕ вертикален). Аналогично для трехчастотного комплекса «Эридан»
2
x
Рис. 3. Калибровочные зависимости д„(Н) (точки) и (Н) (линии) для Северного сектора РАТАН-600, полученные для частот 1, 2.3, 4.8, 7.7, 11.2 и 21.7 ГГц в одном из наблюдательных циклов. Зависимости отражают учет изменения эффективной площади антенны Леяи атмосферного поглощения на коротких волнах от высоты радиоисточника при фактических абберациях из-за поперечного смещения первичных облучателей от электрической оси телескопа.
(вторичное зеркало №2). Для учета поляризации калибровочных источников угол ф рассчитывается исходя из того, что позиционные углы анализаторов радиометров Аа для Северного сектора и Южного сектора с Перископом составляют 90°.
В случае когда сами исследуемые источники не являются точечными, вводится поправка на «конечность» размеров источника, которая принимает значения больше 1 (угловые размеры источника больше размеров ДН телескопа ф\ » фр). Для источника с гауссовым распределением радиояркости и при аппроксимации гауссианой такая поправка К определяется соотношением:
K =
\
1 +
Фи
фрх
1 +
ф,
iy
фРу
где фх и фу — горизонтальный и вертикальный угловые размеры источника, фрх и фру — горизонтальный и вертикальный угловые размеры диаграммы, которые можно определить по известным точечным источникам. Значение коэффициента К
для измерений на РАТАН-600 зависит только от соотношения фх/фрх, так как вторым множителем под знаком корня можно пренебречь (фу/фРу ^ 1) из-за большого вертикального углового размера ДН РАТАН-600 в режиме наблюдений с одним сектором главного зеркала радиотелескопа
(фр » фру ).
Полученные калибровочные кривые для выбранной конфигурации радиотелескопа и для всех длин волн применяются при пересчете антенных температур в плотности потоков для всех исследуемых объектов данного наблюдательного сета (на этапе «Flux density estimation, final results,» рис. 2) согласно формуле:
Svi = Kgv (h) Tav,i,
где индекс «i» относится к параметрам исследуемого объекта.
Суммарная погрешность измерений учитывает погрешность определения антенной температуры в результате аппроксимации кривой прохождения гауссианой, ошибку оценки калибровочной кривой,
2
2
ошибку определения антенной температуры по эталонному источнику шумовой мощности в радиометре. Общая относительная среднеквадратичная ошибка спектральной плотности потока определяется из соотношения:
Sv
Oc
g (h)
+
Om
Ta
1.5
где ot — значение абсолютной погрешности спектральной плотности потока Sv, oc — стандартная ошибка калибровочной кривой, составляющая 2—15% на разных углах в зависимости от длины волны, om — стандартная ошибка измерения антенной температуры Ta,v, gv(h) — функция зависимости калибровочного коэффициента от высоты установки антенны на частоте v. Величина om является оценкой среднеквадратичного отклонения случайных флуктуаций антенной температуры радиоисточника в течение наблюдательного сета. Итоговая относительная среднеквадратичная ошибка измерения спектральной плотности потока исследуемых объектов 5—20% на 11.2, 7.7 и 4.8 ГГц и 10-25% на 1.1, 2.3 и 21.7 ГГц.
В случае использования алгоритмов векторного осреднения многих записей радиоисточника (модули aver или ravr [12]) для повышения отношения сигнал/шум ошибка уменьшается в корень квадратный из числа осредняемых записей.
4.4. Оценка параметров антенной и приемных систем
Система обработки предоставляет информацию о параметрах антенной системы, таких как полуширина главного лепестка ДН (ф^Х, см. формулу для коэффициента K в предыдущем разделе), изменение эффективной площади Aeff, v со временем и параметры радиометрического комплекса (блок «Telescope parameters estimation» на рис. 2). Для измерений параметров антенной системы используются наблюдения двенадцати стандартных точечных калибровочных объектов. Ежедневные измерения калибровочных объектов позволяют своевременно обнаруживать отклонение экспериментальных значений указанных величин от расчетных на каждой высоте (рис. 4 для величины ф^Х).
Теоретическая кривая полуширин получена с использованием набора расчетных диаграмм направленности методом, описанным в работе [27].
Также система обработки графически представляет нормированные значения сигнала калибровки и среднего уровня сигнала для каждого из радиометров полной мощности (рис. 5), что позволяет оперативно реагировать на возможные неисправности в радиометрах.
10 20 30 40 50 60 70 Elevation, deg
80 90
Рис. 4. Зависимость полуширины точечных калибровочных источников на частоте 11.2 ГГц (апрель 2015 г.) от высоты установки антенны (черная линия) в сравнении с теоретической кривой (серая линия), полученной с помощью набора расчетных диаграмм направленности [27] РАТАН-600.
Каждая точка на графиках рис. 5 соответствует одному из двух измерений значений У8у8 и Каиь (рис. 6).
Рассмотрим нижнюю кривую, представленную на рис. 5. При постоянном калибровочном сигнале и переменном усилении в радиометре напряжение на выходе радиометра полной мощности можно записать в форме:
^sys (t) — k T
sys
BjG(t),
где полное усиление С(£) меняется со временем. Здесь к — постоянная Больцмана, Т8у8 — температура системы «радиотелескоп + радиометр,» В — СВЧ-полоса радиометра, 7 — вольт-ваттная характеристика квадратичного детектора. При включении постоянного калибровочного сигнала Канъ напряжение на выходе радиометра получает добавку, равную:
КаИЪ (¿) = кТСаИЬ В7С(Ь).
Эта величина напряжения измеряется во время процедуры калибровки радиометра. Если зафиксировать в произвольный момент времени £ = ¿0 значение КаиЪ(£) (назовем его УсацЪ, 0) — это напряжение будет соответствовать какому-то конкретному усилению, назовем его О0. Тогда для любого другого момента времени можно записать:
Каиъ(£) = кТСаиъ В 7 Со д(£) = Усаиъ,о д(£),
где д(£) — переменная часть усиления. Нижняя кривая на рис. 5 представляет величину
frcalibffl Kcalib, 0
— g(t),
то есть демонстрирует относительные изменения усиления со временем (относительно усиления С0). Кривые на рис. 5 приведены для одного из
2
2
2
1.5 I I I I I I
57423
57424
57425 MJD
57426
57427
Рис. 5. Графики относительного уровня сигнала радиометра полной мощности (верхняя кривая) и относительного уровня калибровочного сигнала (нижняя кривая), построенные для выбранного наблюдательного сета. Пики на верхней кривой соответствуют большому отклонению величины ДТзуа^), от величины ДТауа, о(£), вызванному наблюдениями на разных высотах над горизонтом.
500
а S
a
m
C400
а
A
-1-
А - 6.25 cm Source NGC 7027
Kcaiib (CALIBR1,2)
FITS header data: CALIBR1—138.73 BASE1—355.87
Vsys (BASE1,2)
_I_
100
FITS header data: CALIBR2=138.80 BASE2=355.63
200
300
Time,
Рис. 6. Запись, соответствующая единичному наблюдению на одном из радиометров. В начале и в конце наблюдения включается калибровочный сигнал. Длительность записи соответствует времени, требуемому для того, чтобы радиоисточник зарегистрировали все радиометры, последовательно расположенные на фокальной линии вторичного зеркала.
неохлаждаемых и нетермостатированных радиометров полной мощности. Усиление в радиометре является функцией температуры окружающей среды (синусоидальные суточные колебания на графике относительного калибровочного сигнала, нижняя кривая на рис. 5). Гладкое поведение этой кривой говорит об отсутствии неисправностей в системе калибровки радиометра.
Суточный ход виден и на графике относительного уровня сигнала радиометра полной мощно-
сти (верхняя кривая на рис. 5). Здесь ситуация сложнее. Запишем выражение для измеряемого напряжения на выходе радиометра полной мощности с учетом изменений температуры системы «радиотелескоп + радиометр» (ДТ8у8) и усиления со временем:
Узув(^) = кТу3(Ь) БС(1)
Так же, как и для построения кривой относительного усиления, зафиксируем момент времени ¿0, которому будут соответствовать усиление С0
0
0
s
и температура системы Твув, 0 и выходное напряжение радиометра 0 = к Твув, 0 БС0 7. Тогда в любой другой момент времени усиление С(£) будет равно С0 д^), а температура системы будет равна Тзув, 0 + ДТвув(^). Подставив в формулу для Узув (см. первую формулу этого раздела) эти зависимости и разделив полученное выражение на Увув, 0, окончательно получим:
^sys(i) Л + A Tsys(i)
Vs
sys, 0
T
sys, 0
g(t).
Из этого выражения видно, что относительные изменения уровня сигнала радиометра полной мощности отражают произведение ДТзуз(1) д(£). Вклад каждого из сомножителей в это произведение представляет отдельный интерес для радиометристов и является предметом текущих исследований.
Измерения хода усиления радиометра, представленные на рис. 5, вместе с измерениями физической температуры в приемной кабине вторичного зеркала позволили численно оценить температурную зависимость усиления радиометров. Она примерно равна 1.1% на 1°С в рабочем температурном диапазоне радиометров (10—30°С).
5. РЕЗУЛЬТАТЫ
Представленная система обработки радиоастрономических наблюдений прошла успешную апробацию, результаты опубликованы, например, в [8—10, 15]. Впервые автоматизирована полная процедура обработки наблюдений компактных объектов, включающая:
► расчет калибровочных зависимостей на всех длинах волн;
измерение спектральных плотностей потоков радиоисточников на нескольких частотах одновременно, построение кривых блеска;
построение радиоспектров источников;
измерение основных параметров радиотелескопа: изменения эффективной площади Ае^; и и ширины главного лепестка ДН антенны ф^х;
► оценку основных параметров качества радиометров континуума.
Автоматизация процесса обработки наблюдательных данных радиотелескопа РАТАН-600 становится все более актуальной в связи с запланированным увеличением числа радиометров континуума.
БЛАГОДАРНОСТИ
Авторы выражают благодарность Ю. А. Ковалеву за полезные замечания и обсуждение результатов. Работа выполнена при финансовой поддержке Министерства образования и науки Российской Федерации (государственный контракт №14.518.11.7054) и Российского фонда фундаментальных исследований (проект №12-02-31649).
СПИСОК ЛИТЕРАТУРЫ
1. D. V. Korolkov and Yu. N. Pariiskii, Sky Telesc. 57, 324(1979).
2. Y. N. Parijskij, IEEE Antennas Propagation Mag. 35, 7(1993).
3. M. G. Mingaliev, Y. V. Sotnikova, N. N. Bursov, et al., Astronomy Reports 51, 343 (2007).
4. M. G. Mingaliev, Y. V. Sotnikova, N. S. Kardashev, and M. G. Larionov, Astronomy Reports 53, 487 (2009).
5. M. G. Mingaliev, Y. V. Sotnikova, M. G. Larionov, and A. K. Erkenov, Astronomy Reports 55, 187 (2011).
6. P. A. R. Ade et al. (Planck Collab.), Astron. and Astrophys. 536, A15 (2011).
7. P. A. R. Ade et al. (Planck Collab.), Astron. and Astrophys. 536, A18 (2011).
8. P. Giommi, G. Polenta, A. Lahteenmaki, et al., Astron. and Astrophys. 541, A160 (2012).
9. M. G. Mingaliev, Y. V. Sotnikova, I. Torniainen, et al., Astron. and Astrophys. 544, A25 (2012).
10. M. G. Mingaliev, Y. V. Sotnikova, T. V. Mufakharov, et al., Astrophysical Bulletin 68, 262 (2013).
11. N. N. Bursov, PhD Thesis (Spec. Astrophys. Obs., Nizhnij Arkhyz, 2002).
12. O. V. Verkhodanov, ASP Conf. Ser., 125,46(1997).
13. Y. A. Kovalev, Doctoral Dissertation in Mathematics and Physics (Physical Inst. RAS, Moscow, 2003).
14. O. V. Verkhodanov, V. K. Kononov, E. K. Majorova, and P. G. Tsybulev, Report No. 233, SAO RAS (Spec. Astrophys. Observ., Nizhnij Arkhyz, 1994).
15. M. G. Mingaliev, Y. V. Sotnikova, R. Y. Udovitskiy, et al., Astron. and Astrophys. 572, A59 (2014).
16. O. V. Verkhodanov, B. L. Erukhimov, O. P. Zhelenkova, et al., Report No. 200, SAO RAS (Spec. Astrophys. Observ., Nizhnij Arkhyz, 1991).
17. P. G. Tsybulev, Astrophysical Bulletin 66, 109 (2011).
18. K. D. Aliakberov, M. G. Mingaliev, M. N. Naugolnaya, et al., Астрофиз. исслед. (Известия Спец. астрофиз. обс.) 19, 60(1985).
19. M. G. Mingaliev, V. A. Stolyarov, R. D. Davies, et al., Astron. and Astrophys. 370, 78 (2001).
20. M. Ott, A. Witzel, A. Quirrenbach, et al., Astron. and Astrophys. 284, 331 (1994).
21. J. W. M. Baars, R. Genzel, I. I. K. Pauliny-Toth, and A. Witzel, Astron. and Astrophys. 61, 99 (1977).
22. R. A. Perley and B. J. Butler, Astrophys. J. Suppl. 204, 19(2013).
23. H. Tabara and M. Inoue, Astron. and Astrophys. Suppl. 39,379(1980).
24. Y. N. Pariiskii, Izv. Glavnoj Astron. Obs. 188, 149
Astron. and Astrophys. Suppl. 139,545(1999).
(1972).
25. J. D. Kraus, Radio Astronomy, 2nd ed., (Cygnus—
27. A. N. Korzhavin, Астрофиз. исслед. (Известия Спец. астрофиз. обс.) 9,71 (1977).
Quasar Books, Powell, 1986).
26. Y. Y. Kovalev, N. A. Nizhelsky, Y. A. Kovalev, et al.,
Automated System for Reduction of Observational Data on RATAN-600 Radio Telescope
R. Yu. Udovitskii, Yu. V. Sotnikova, M. G. Mingaliev, P. G. Tsybulev, G. V. Zhekanis, and
N. A. Nizhelskii
We present the automated system for estimating the parameters of radio sources observed on all available continuum radiometers (two receiving complexes of secondary mirrors No. 1 and No. 2 with a total of 30 radiometers) developed at RATAN-600 radio telescope and put into normal operation. The system is also used for the monitoring of the parameters of the antenna and receiving systems of RATAN-600 radio telescope, which is carried out using current measurements of calibration radio sources.
Keywords: radio continuum: general—methods: observational—techniques: radar astronomy