Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://comet.sai.msu.ru/~dmbiz/prac/next/zad17/node4.html
Дата изменения: Fri Sep 7 17:25:57 2001 Дата индексирования: Tue Oct 2 04:17:09 2012 Кодировка: koi8-r Поисковые слова: п п п п п п п п п п п п р п р п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п |
Для выполнения задачи необходимо иметь 28 файлов: 4 файла изображения галактики GRB980703 в фильтрах B, V, R, I (galb.bdf, galv.bdf, galr.bdf, gali.bdf), 8 файла площадок звездных стандартов RU149 (ru149b.bdf, ru149v.bdf, ru149r.bdf, ru149i.bdf), 12 файлов "плоского поля" матрицы (flat field) в фильтрах B, V, R, I (ffb1.bdf - ffb3.bdf, ffv1.bdf - ffv3.bdf, ffr1.bdf - ffr3.bdf, ffi1.bdf - ffi3.bdf - по три в каждом фильтре), 3 файла измерения темнового тока матрицы (dark1.bdf - dark3.bdf) и 5 файлов измерения уровня сигнала смещения матрицы (bias1.bdf - bias5.bdf).
Перед выполнением задачи обязательно ознакомьтесь
с описанием простейших команд Linux и MIDAS, приведенных в
Приложении.
Создать рабочую директорию, чтобы затем можно было продолжать в ней обработку [например, mkdir ivan_23022001].
Файлы, необходимые для работы, содержаться в директории initial_data/24/
Cкопировать 28 файлов в директорию, в которой будет выполняться работа, командой cp /initial_data/24/* /ваша_рабочая_директория
В X-Windows открыть окно и запустить MIDAS командой inmidas -p 0x , где x - цифра от 0 до 9. Данная команда позволяет запускать одновременно несколько сессий MIDAS на одной и той же машине.
С помощью команды create/icat создать файл- каталог байесов. Затем, подключив дополнительную библиотеку команд ccdred (с помощью команды set/context), построить супербайес, используя команду combine/ccd.
Исправить изображения галактики, звезд сравнения, плоских полей и темновых шумов (т.е. оставшиеся 23 файла) за байес, вычитая его из указанных изображений командой compute/pixel.
Итак, следующий шаг - с помощью команд create/icat и combine/ccd построить т.н. супердарк (объединенный файл темновых экспозиций).
Полученное изображение необходимо вычесть из изображений объекта, звезд сравнения и плоских полей (всего 20 файлов) с учетом времени их экспозиции. Для этого рекомендуется вычитать супердарк, нормированный на соответствующие времена экспозиций галактики, звезд сравнения и плоских полей с помощью команды compute/pixel. Время экспозиции изображений можно посмотреть в дескрипторах (заголовках) файлов (ключевое слово o_time) при помощи команды read/descr. Седьмая (последняя) цифра в массиве есть время экспозиции кадра в секундах.
Построить с помощью команд create/icat и combine/ccd изображения с "поскими полями", или т.н. суперфлеты. Они будут разными для разлиных фильтров (B, V, R, I).
Определить, используя команду statistics/image, медианное среднее каждого кадра плоского поля. Скомбинировать плоские поля, получая суперфлеты для каждого фильтра. Суммировать следует плоские поля, нормированные на их медианное среднее. (команда compute/pixel).
Разделить изображения галактики и звезд сравнения на нормированные суперфлеты для каждого из четырех фильтров (compute/pixel). Оценить градиент изменения фона вдоль и поперек кадров с изображением стандартов и объектов. Оценить (по порядку величины), с какой точностью можно получить оценку звездной величины объекта (эти оценки пригодятся в дальнейшем, когда мы будем оценивать ошибки, с которыми получены величины объектов).
В результате мы получим 8 изображенй: 4 файла с объектом и 4 файла с изображением звезд фотометрического стандарта, исправленные за сигнал смещения матрицы, темновой ток и плоское поле - по одному изображению в каждом фильтре.
Провести калибровку по звездам фотометрического стандарта. Для этого с помощью команды magnitude/circle определить интегральные интенсивности звезд RU149, RU149A, RU149B, RU149C, RU149D, RU149E, RU149F и RU149G на каждом изображении. Обратите внимание на то, что размер апертуры, в которой производится аппроксимация гауссианой, задается в команде в явном виде и не зависит от размера кружка, изображенного на экране. Размер апертур выбрать таким образом, чтобы звезда полностью попадала в нее, но не слишком большим, чтобы флуктуации фона не оказывали заметного влияния на значения интенсивности. Обратите внимание, что результаты можно автоматически заносить в таблицу.
Примечание: При выполнения этих операции удобно пользоваться редактором таблиц в MIDAS (команды create/table, create/column, edit/table, compute/table, name/column).
В результате мы получим калибровочные уравнения, которые позволят перевести суммарные потоки (отсчеты) в некоторых апертурах в звездные величины. При расчетах не забывайте приводить потоки от объектов к экспозиции в 1 сек.
Найти сумарные потоки от галактики в фильтрах B, V, R, I, используя команду magnitude/circle. Размеры диафрагмы должны быть одинаковыми в каждом из четырех изображений. Центры - также по возможности близкими.
Определить ошибки измерений звездных величин галактики. Учтите ошибки фотометрической привязки по соотвествующему цветовому уравнению, ошибки определения потоков от галактики. Оцените ошибки определения показателей цвета.
тип | U-B | B-V | V-R | R-I |
E | -0.52 | 0.83 | 1.49 | 1.57 |
S0 | -0.53 | 0.78 | 1.53 | 1.44 |
Sab | -0.50 | 0.38 | 1.60 | 1.46 |
Sbc | -0.62 | 0.72 | 0.81 | 1.22 |
Scd | -0.61 | 0.48 | 0.69 | 1.19 |
Im | -0.77 | 0.28 | 0.29 | 0.91 |
Примечание. К-поправки вычитаются!
Примечание. Полученные показатели цвета не должны совпадать со средними показателями цвета галактик различных морфологических типов из Табл. 1, поскольку исследуемая в задаче GRB980703 является галактикой со вспышкой звездообразования. Для определения морфологического типа галактики нужно использовать в первую очередь показатель цвета V-R.
Морфологический тип | К(U-B) | К(B-V) | К(V-R) | К(R-I) | К(V) |
E | -1.16 | -0.13 | 0.88 | 0.87 | 3.18 |
S0 | -0.95 | -0.07 | 0.99 | 0.83 | 2.82 |
Sab | -0.83 | -0.50 | 1.04 | 0.81 | 2.95 |
Sbc | -0.61 | -0.17 | 0.29 | 0.60 | 1.67 |
Scd | -0.53 | -0.02 | 0.19 | 0.62 | 1.12 |
Im | -0.42 | 0.01 | -0.02 | 0.58 | 0.29 |