Сравнительная фотометрия большого числа объектов

[1212Lupus]

Есть следующая задача: получаем в результате обработки (нормировка, сложение и т.д.) фотографию, содержащую большое колличество звёзд. Возникает задача -- нужно задать опорные звёзды (несколько штук), указать их звёздную величину, а затем программа должна определить звёздную для всех остальных звёзд (понятно, что кликанье мышкой как в TeleAuto не подходит -- слишком много звёзд), видимых на фотографии (естесственно, вывод значений лучше осуществлять не на экран, а в файл -- а то мешанина получится). Затем программа должна отождествить данный участок неба с каталогом звёзд (естессвенно, заранее нужно задать положение центра кадра), и сравнить полученные значения звёздных величин со значениями из каталога. Если отклонение превысит, скажем, 0,5^m, то пишется лог-файл, где указывается положение этой звезды в пикселах.

У кого есть какие идеи -- как и при помощи какого софта это можно организовать? Поделитесь, пожалуйста.

Каталог USNO-A2.0 (6 ГБ) есть скачаный -- в NASA'вском стандарте (отсюдова -- ftp://ftp.nofs.navy.mil/usnoa/).

[Garmisch]

Вкратце, это можно сделать в IRAF при помощи пакета DAOPHOT. Из коммерческих продуктов стоит, наверное, обратить внимание на Mira Pro http://www.mirametrics.com/mira_pro.htm

[1212Lupus]

А что скажете насчёт PIXY2?
http://www.aerith.net/misao/pixy/index.html

[Garmisch]

А что скажете насчёт PIXY2?

Не знаю, нужно ставить и смотреть

[1212Lupus]

Вчера поставил PIXY2. Она написана на java, так что по-определению является кросплатформенной. Проверил, как работает пол Linux OS и Windows. Для затравки загрузил наш снимок квазара 3С273 (сложение и обработка 5 фотографий), сделанный из Минска. Там проработались звёзды до 14^m. Предварительно выкачал несколько каталогов, в том числе Tycho2. Звёзды отождествляютсянормально. Только одна трабла есть -- не открывается FITS 3000х2000 пикселей -- пришлось кропить, картинки поменьше открывает шустро. Возможно, я просто не дождался момента открытия -- слишком мало памяти выделил для работы проги (память можно выделять как тебе захочется -- хоть 200 МБ, хоть 500 -- смотря, сколько имеется в распоряжении).

[Евгений Ромас(BrainBug)]

Вкратце, это можно сделать в IRAF при помощи пакета DAOPHOT.

Юра, а можно поподробнее? Где почитать? HELP по DAOPHOT не предлагать...  

[Stefan Stefens]

Есть следующая задача: получаем в результате обработки (нормировка, сложение и т.д.) фотографию, содержащую большое колличество звёзд. Возникает задача -- нужно задать опорные звёзды (несколько штук), указать их звёздную величину, а затем программа должна определить звёздную для всех остальных звёзд (понятно, что кликанье мышкой как в TeleAuto не подходит -- слишком много звёзд), видимых на фотографии (естесственно, вывод значений лучше осуществлять не на экран, а в файл -- а то мешанина получится). Затем программа должна отождествить данный участок неба с каталогом звёзд (естессвенно, заранее нужно задать положение центра кадра), и сравнить полученные значения звёздных величин со значениями из каталога. Если отклонение превысит, скажем, 0,5m, то пишется лог-файл, где указывается положение этой звезды в пикселах.

У кого есть какие идеи -- как и при помощи какого софта это можно организовать? Поделитесь, пожалуйста.

Натолкнулся на такой же вопрос. Посмотрел разные программы, в основном с кликаньем мышкой по всему полю.
Сам я программист и не хочется верить, что нельзя проше. Вот начал писать сам. В программе всего три кнопки:
Выбор бызы, коррекция в списке снимков ( например на положение опорных звезд  / случаи неточностей монтировки) и интерпретация.
есть конечно вопросы. было бы здорово, если бы кто-нибудь помог. вопрос первый: как из интенсивности в пикселях посчитать звездную величину не теоретически , а практически? Буду очень благодарен за помошь. Часто на форуме возникают ненужные дискуссии, прошу их избежать в этом случае. Программа рождается на Яве, каталоги пока локально не подсоединяю: в раздумиях - надо ли. Итак: как реально из RGB суммы пикселей получить звездную величину одиночно стояшей звезды без перенасышения и с перенасышением? Спасибо заранее 

Все параметры аппроксимируюшего звезное изображение Гаусиана посчитаны с довольно высокой точность (в рамках DLSR сьемки). Их хотелось бы дл я вычисления и применить. Конечно тоже сделано для звезд сравнения ( для база ).

[Евгений Ромас(BrainBug)]

http://www.physics.hmc.edu/Astronomy/Iphot.html

Photometry with daofind

Daofind is a useful task for examining star clusters, because IRAF will automatically find stars using an algorithm and do photometry on them. To access this task, open the digiphot and then apphot packages.

...и далее.

[Stefan Stefens]

Спасибо за ссылку. Есть пакеты и с использованием нейронных сетей -    SExtractor . Все это не отличается простотой. А главное не отвечает на вопрос, какой алгоритм/формула используется для перевода интенсивности в звездную величину? Надеюсь кто-нибудь поможет ...

[Евгений Ромас(BrainBug)]

Спасибо за ссылку. Есть пакеты и с использованием нейронных сетей -    SExtractor . Все это не отличается простотой. А главное не отвечает на вопрос, какой алгоритм/формула используется для перевода интенсивности в звездную величину? Надеюсь кто-нибудь поможет ...

http://astrozone.ru/um/3kurs/photometry/bookpm.pdf
как-то так... а вообще гугаль рулед. 

[Stefan Stefens]

Спасибо. Попробую погрызть и это. Немецкие издания уже погрыз, но что интересно, теории сколько хочешь, а как реально считать - ни гу гу. А так в целом получается. Вот только карты сравнения через  Web сервисы негде пока взять. Как уже писал, нет желания цеплять каталоги напрямую. Чуть что изменилось, и нельзя пользоваться.

[krussh]

http://www.adass.org/adass/proceedings/adass98/mighellkj/

здесь все есть.

[Stefan Stefens]

Спасибо за ссылку. Хорошая ссылка. Алгоритм поиска центра я давно уже встроил, распознавание обьекта - звезда не звезда тоже. А вот интенсивность в звездную величину так и остается теорией. В чем проблемма? если я имею несколько звезд сравнения в примерно равных условиях, то пересчет нелинеен. Например для слабых звезд ошибка столь велика, что нет смысла мерить, для ярких - они перенасышены - можно выташить правильную интенсивность или нельзя - не знаю. Где лежат границы на конкретном снимке, какова нелинейность (тут я уже некоторые аппроксимации сам обнаружил)? Как-то же в других пакетах все это учтено?  Или нет? Я пробую все на JPG . Это мой цель формат. Все должно быть просто и доступно всем. И тем кто первый раз снимает Алголь и тем кто хочет сранивать серии снимков.

[krussh]

Пересчитать интенсивность (сумму отсчетов в апертуре) в инструментальную звездную величину не сложно. Сложно формулу вставить в сообщение форума). Скажите куда выслать три пдфа с подробными пояснениями.

Если есть опорные звезды, то не сложно посчитать и абсолютную зв. величину. Но тут все не так просто. Лучше использовать ансамбль близких звезд сравнения с весами пропорциональными их ошибкам. Перечет линеен, в этом и есть прелесть ПЗС. Отсчет пропорционален потоку.

Слабые звезды можно выкинуть - ваше право. Яркие содержат отсчеты вблизи верхней границы линейного диапазона ПЗС. Для каждой камеры это свое число. Яркие надо выкинуть однозначно.

Ну и жпг это будет не очень хорошо.

[krussh]

http://www.astro.louisville.edu/software/astroimagej/

[Stefan Stefens]

Этот ресурс мне хорошо известен. Все равно спасибо.Прежде чем задуматься писать, я конечно многое испробовал. И везде было что-нибудь не так. Наибольшие проблеммы были с оборудованием для сьемки и конфигурацией программ под это оборудование. Вот и вздумалось все упростить. JPG я взял также из-за простоты и компактности в больших сериях. По моим ( может быть неправильным) рассуждениям - берем зеленый канал и имеем максимум 255 значений по RGB. 50 - может быть фон. Остается 200. Предел яркости на 120мм рефракторе и Canon 450D получается 9 звездной величине  (the exposition = 1min  - пересышение видно) и минимум мыслимый 14 звездная величина ( to see 15.5m on the picture). И так 5 звезных величин на 200, это порядка 0.025 и если использовать уточнение профиля с интерполяцией, то 0.015-0.02 - реальность. Или что-то в моих рассуждениях выпушено?

[1212Lupus]

Все параметры аппроксимируюшего звезное изображение Гаусиана посчитаны с довольно высокой точность (в рамках DLSR сьемки). Их хотелось бы дл я вычисления и применить. Конечно тоже сделано для звезд сравнения ( для база ).

Скажите, а фотомтерию треков на основе фотометрии опорных точечных звёзд, как ТУТ, нет желания сделать?

[krussh]

И так 5 звезных величин на 200, это порядка 0.025 и если использовать уточнение профиля с интерполяцией, то 0.015-0.02 - реальность. Или что-то в моих рассуждениях выпушено?

5 зв. величин, это отличие по потоку в 100 раз, шкала зв. величин логарифмическая.

[Stefan Stefens]

Скажите, а фотомтерию треков на основе фотометрии опорных точечных звёзд, как ТУТ, нет желания сделать?

Нет. Так легко дойти до сейсмологии звезд. Время трека - соизмеримо со временами интенсивных искажений, которые в ссылке перечисленны. Поэтоми считаю задачу безнадежной, а вот звучашие звезды в ультразвуке имеют меньше естественных припятствий кроме АПЕРТУРЫ ( ФЭУ подходят)

[Stefan Stefens]

5 зв. величин, это отличие по потоку в 100 раз, шкала зв. величин логарифмическая.

Тогда очень важен выбор экспозиции для того чтобы попасть на мах чувствительность. Интересно. как автоматически подсказать - туда нельзя, а здесь можно?