Поиск звездообразных объектов на кадре

Сначала, замечание: Практически при каждом запуске различных команд DAOPHOT'а будут запрашиваться названия входных и выходных файлов, однако, после запуска команды attach при запросе имен файлов будут предлагаться варианты по умолчанию (отличающиеся от имени исходного файла только по расширению (.coo,.ap,.grp,.psf и т.д.)) Поэтому в дальнейшем все, что касается запроса командами имен входных и выходных файлов, упоминаться не будет. Также будем считать, что мы обрабатываем файл cluster.bdf .

Для поиска звездообразных объектов в DAOPHOT существует команда FIND. Это одна из основных программ в DAOPHOT'е. Важно то, что на ее работу существенным образом влияют значения параметров, описанных в пункте 2.1. После запуска команды FIND, пользователя запросят значения двух параметров $s$ и $a$(вводить их надо через запятую), если обрабатываемое изображение $Image=\frac{1}{a}*\sum\limits_{i=1}^{i=a}(\sum\limits_{j=1}^{j=s}Image_{ij})$ , где $Image_{ij}$ различные индивидуальные изображения. После задания этих параметров запускается процесс поиска объектов.

После прохода по всем строчкам (поиск осуществляется именно построчно) программа выводит список обнаруженных объектов c указанием номера, приблизительных X и Y координат центра объекта, звездной величины по отношению к пределу обнаружения, а также значение параметров SHARP и ROUND для каждого объекта. Параметр SHARP равен отношению интенсивности двумерной дельта-функции вписанной в объект и интенсивности двумерной вписанной гауссианы (допустимое значение по умолчанию $0.2 \le SHARP \le 1.00$Кстати, объекты со значением параметра SHARP, существенно большим нуля являются, либо галактиками, либо неразрешенными двойными. Объекты с SHARP, существенно меньшим нуля являются, либо космическими частицами, либо дефектными пикселями. Параметр ROUND равен разности центральных яркостей двух одномерных гауссиан, вписанных, соответственно, вдоль Х и У, поделенной на среднюю центральную яркость вышеупомянутых одномерных гауссиан (допустимое по умолчанию значение $-1.0 \le ROUND \le 1.0$.

После вывода списка обнаруженных объектов, пользователю будет задан вопрос о том, удовлетворен ли он результатами работы программы FIND (количеством объектов). Если пользователь ответит, что нет, то пользователю будет предложено изменить значение параметра THRESHOLD (см. п. 2.1.5)

Формат таблицы, создаваемой командой FIND (cluster.coo):

1. Собственный порядковый номер объекта (присваивается DAOPHOT'ом обьекту раз и навсегда, т.е. все запущенные впоследствии команды будут оперировать с созданной командой FIND нумерацией; нумерация будет изменена, только при следующем запуске команды FIND для этого файла).
2. X координата центра объекта.
3. Y координата центра объекта.
4. Очень грубое значение звездной величины, по отношению к нуль-пункту, определенному значением параметра THRESHOLD.
5. Значение параметра SHARP для данного объекта.
6. Значение параметра ROUND для данного объекта.

Кстати, замечание о структуре ASCII таблиц создаваемых DAOPHOT'ом. В заголовке каждой такой таблицы выписаны, в частности, значения основных параметров DAOPHOT'а. Заголовки разных файлов немного отличаются Например:

NL	NX	NY	LOWBAD	HIGHBAD	THRESH	AP1	PH/ADU	RNOISE	FRAD
1	1752	1751	657.2	65500.0	90.68	2.00	1.34	6.00	5.50

Итак структура заголовка:

• NL Это число, характеризующее ту команду, которая создала данную таблицу. Поэтому у таблиц с разными расширениями (als, coo , lst и проч.) должны быть разные значения NL.
• NX Это количество элементов в изображении по оси X.
• NY Это количество элементов в изображении по оси Y.
• LOWBAD Это значение параметра LOW GOOD DATUM в единицах ADU.
• HIGHBAD Это значение параметра HIGH GOOD DATUM в единицах ADU.
• THRESH Это значение параметра THRESHOLD в единицах ADU.
• AP1 Это радиус первой апертуры (A1).
• PH/ADU Это Gain.
• RNOISE Это величина шума считывания в единицах ADU.
• FRAD Это значение параметра FITTING RADIUS.