Телескопы покупают здесь


A A A A Автор Тема: Сравнительная фотометрия большого числа объектов  (Прочитано 3541 раз)

0 Пользователей и 1 Гость просматривают эту тему.

Оффлайн 1212LupusАвтор темы

  • *****
  • Сообщений: 3 078
  • Рейтинг: +194/-4
  • Мне стал не очень нравиться этот форум...
    • Show only replies by 1212Lupus
    • http://belastro.net
Есть следующая задача: получаем в результате обработки (нормировка, сложение и т.д.) фотографию, содержащую большое колличество звёзд. Возникает задача -- нужно задать опорные звёзды (несколько штук), указать их звёздную величину, а затем программа должна определить звёздную для всех остальных звёзд (понятно, что кликанье мышкой как в TeleAuto не подходит -- слишком много звёзд), видимых на фотографии (естесственно, вывод значений лучше осуществлять не на экран, а в файл -- а то мешанина получится). Затем программа должна отождествить данный участок неба с каталогом звёзд (естессвенно, заранее нужно задать положение центра кадра), и сравнить полученные значения звёздных величин со значениями из каталога. Если отклонение превысит, скажем, 0,5m, то пишется лог-файл, где указывается положение этой звезды в пикселах.

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

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


Если утро наступает в три -
Через два часа уже зажгут фонари.
Уже кончился день, а я только встал,
А я только что встал и уже устал.
(с) НОЛЬ

Оффлайн Garmisch

  • ASTRONOMY.RU
  • *****
  • Сообщений: 8 016
  • Рейтинг: +293/-1
  • Награды Неоднократный победитель конкурса астрофото
    • Show only replies by Garmisch
    • http://500px.com/ybeletsky
Вкратце, это можно сделать в IRAF при помощи пакета DAOPHOT. Из коммерческих продуктов стоит, наверное, обратить внимание на Mira Pro http://www.mirametrics.com/mira_pro.htm

Экспедиция в Атакаму: https://astronomy.ru/forum/index.php/topic,117810.0.html
Моя галерея: http://500px.com/ybeletsky
Моя страница в Google+: http://google.com/+YuriBeletsky
Twitter: https://twitter.com/YBeletsky

Оффлайн 1212LupusАвтор темы

  • *****
  • Сообщений: 3 078
  • Рейтинг: +194/-4
  • Мне стал не очень нравиться этот форум...
    • Show only replies by 1212Lupus
    • http://belastro.net
А что скажете насчёт PIXY2?
http://www.aerith.net/misao/pixy/index.html
Радиоастрономы-любители -- объединяемся!


Если утро наступает в три -
Через два часа уже зажгут фонари.
Уже кончился день, а я только встал,
А я только что встал и уже устал.
(с) НОЛЬ

Оффлайн Garmisch

  • ASTRONOMY.RU
  • *****
  • Сообщений: 8 016
  • Рейтинг: +293/-1
  • Награды Неоднократный победитель конкурса астрофото
    • Show only replies by Garmisch
    • http://500px.com/ybeletsky
А что скажете насчёт PIXY2?

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

Экспедиция в Атакаму: https://astronomy.ru/forum/index.php/topic,117810.0.html
Моя галерея: http://500px.com/ybeletsky
Моя страница в Google+: http://google.com/+YuriBeletsky
Twitter: https://twitter.com/YBeletsky

Оффлайн 1212LupusАвтор темы

  • *****
  • Сообщений: 3 078
  • Рейтинг: +194/-4
  • Мне стал не очень нравиться этот форум...
    • Show only replies by 1212Lupus
    • http://belastro.net
Вчера поставил PIXY2. Она написана на java, так что по-определению является кросплатформенной. Проверил, как работает пол Linux OS и Windows. Для затравки загрузил наш снимок квазара 3С273 (сложение и обработка 5 фотографий), сделанный из Минска. Там проработались звёзды до 14m. Предварительно выкачал несколько каталогов, в том числе Tycho2. Звёзды отождествляютсянормально. Только одна трабла есть -- не открывается FITS 3000х2000 пикселей -- пришлось кропить, картинки поменьше открывает шустро. Возможно, я просто не дождался момента открытия -- слишком мало памяти выделил для работы проги (память можно выделять как тебе захочется -- хоть 200 МБ, хоть 500 -- смотря, сколько имеется в распоряжении).
Радиоастрономы-любители -- объединяемся!


Если утро наступает в три -
Через два часа уже зажгут фонари.
Уже кончился день, а я только встал,
А я только что встал и уже устал.
(с) НОЛЬ

Оффлайн Евгений Ромас(BrainBug)

  • Первооткрыватель астероидов
  • *****
  • Сообщений: 5 843
  • Рейтинг: +145/-5
  • КУ!
    • Skype - jeka_brainbug
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Show only replies by Евгений Ромас(BrainBug)
Вкратце, это можно сделать в IRAF при помощи пакета DAOPHOT.
Юра, а можно поподробнее? Где почитать? HELP по DAOPHOT не предлагать... ;)  ;D
Celestron SCT 6" XLT, EQ5 + EQDream - в продаже.
SkyWatcher BKP 200 DS, HEQ5 Pro(временно), SBIG ST-402ME, гид Datyson T7m(клон ASI120MM) + 60/240.
malum edo ipsum

Оффлайн Stefan Stefens

  • **
  • Сообщений: 98
  • Рейтинг: +4/-0
    • Show only replies by Stefan Stefens
Есть следующая задача: получаем в результате обработки (нормировка, сложение и т.д.) фотографию, содержащую большое колличество звёзд. Возникает задача -- нужно задать опорные звёзды (несколько штук), указать их звёздную величину, а затем программа должна определить звёздную для всех остальных звёзд (понятно, что кликанье мышкой как в TeleAuto не подходит -- слишком много звёзд), видимых на фотографии (естесственно, вывод значений лучше осуществлять не на экран, а в файл -- а то мешанина получится). Затем программа должна отождествить данный участок неба с каталогом звёзд (естессвенно, заранее нужно задать положение центра кадра), и сравнить полученные значения звёздных величин со значениями из каталога. Если отклонение превысит, скажем, 0,5m, то пишется лог-файл, где указывается положение этой звезды в пикселах.

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

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

Все параметры аппроксимируюшего звезное изображение Гаусиана посчитаны с довольно высокой точность (в рамках DLSR сьемки). Их хотелось бы дл я вычисления и применить. Конечно тоже сделано для звезд сравнения ( для база ).
« Последнее редактирование: 25 Апр 2013 [19:01:00] от Stefan Stefens »
TLAPO804, ZEQ25GT, ZWO Finder/Guider, Orion Flattener, MGEN-2

Оффлайн Евгений Ромас(BrainBug)

  • Первооткрыватель астероидов
  • *****
  • Сообщений: 5 843
  • Рейтинг: +145/-5
  • КУ!
    • Skype - jeka_brainbug
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Show only replies by Евгений Ромас(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.
...и далее.
Celestron SCT 6" XLT, EQ5 + EQDream - в продаже.
SkyWatcher BKP 200 DS, HEQ5 Pro(временно), SBIG ST-402ME, гид Datyson T7m(клон ASI120MM) + 60/240.
malum edo ipsum

Оффлайн Stefan Stefens

  • **
  • Сообщений: 98
  • Рейтинг: +4/-0
    • Show only replies by Stefan Stefens
Спасибо за ссылку. Есть пакеты и с использованием нейронных сетей -    SExtractor . Все это не отличается простотой. А главное не отвечает на вопрос, какой алгоритм/формула используется для перевода интенсивности в звездную величину? Надеюсь кто-нибудь поможет ...
TLAPO804, ZEQ25GT, ZWO Finder/Guider, Orion Flattener, MGEN-2

Оффлайн Евгений Ромас(BrainBug)

  • Первооткрыватель астероидов
  • *****
  • Сообщений: 5 843
  • Рейтинг: +145/-5
  • КУ!
    • Skype - jeka_brainbug
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Show only replies by Евгений Ромас(BrainBug)
Спасибо за ссылку. Есть пакеты и с использованием нейронных сетей -    SExtractor . Все это не отличается простотой. А главное не отвечает на вопрос, какой алгоритм/формула используется для перевода интенсивности в звездную величину? Надеюсь кто-нибудь поможет ...
http://astrozone.ru/um/3kurs/photometry/bookpm.pdf
как-то так... а вообще гугаль рулед.  ;D
Celestron SCT 6" XLT, EQ5 + EQDream - в продаже.
SkyWatcher BKP 200 DS, HEQ5 Pro(временно), SBIG ST-402ME, гид Datyson T7m(клон ASI120MM) + 60/240.
malum edo ipsum

Оффлайн Stefan Stefens

  • **
  • Сообщений: 98
  • Рейтинг: +4/-0
    • Show only replies by Stefan Stefens
Спасибо. Попробую погрызть и это. Немецкие издания уже погрыз, но что интересно, теории сколько хочешь, а как реально считать - ни гу гу. А так в целом получается. Вот только карты сравнения через  Web сервисы негде пока взять. Как уже писал, нет желания цеплять каталоги напрямую. Чуть что изменилось, и нельзя пользоваться.
« Последнее редактирование: 26 Апр 2013 [14:05:09] от Stefan Stefens »
TLAPO804, ZEQ25GT, ZWO Finder/Guider, Orion Flattener, MGEN-2

Оффлайн krussh

  • *****
  • Сообщений: 2 272
  • Рейтинг: +102/-5
  • http://www.soulfulbits.com/
    • Show only replies by krussh

Оффлайн Stefan Stefens

  • **
  • Сообщений: 98
  • Рейтинг: +4/-0
    • Show only replies by Stefan Stefens
Спасибо за ссылку. Хорошая ссылка. Алгоритм поиска центра я давно уже встроил, распознавание обьекта - звезда не звезда тоже. А вот интенсивность в звездную величину так и остается теорией. В чем проблемма? если я имею несколько звезд сравнения в примерно равных условиях, то пересчет нелинеен. Например для слабых звезд ошибка столь велика, что нет смысла мерить, для ярких - они перенасышены - можно выташить правильную интенсивность или нельзя - не знаю. Где лежат границы на конкретном снимке, какова нелинейность (тут я уже некоторые аппроксимации сам обнаружил)? Как-то же в других пакетах все это учтено?  Или нет? Я пробую все на JPG . Это мой цель формат. Все должно быть просто и доступно всем. И тем кто первый раз снимает Алголь и тем кто хочет сранивать серии снимков.
TLAPO804, ZEQ25GT, ZWO Finder/Guider, Orion Flattener, MGEN-2

Оффлайн krussh

  • *****
  • Сообщений: 2 272
  • Рейтинг: +102/-5
  • http://www.soulfulbits.com/
    • Show only replies by krussh
Пересчитать интенсивность (сумму отсчетов в апертуре) в инструментальную звездную величину не сложно. Сложно формулу вставить в сообщение форума). Скажите куда выслать три пдфа с подробными пояснениями.

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

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

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

Оффлайн krussh

  • *****
  • Сообщений: 2 272
  • Рейтинг: +102/-5
  • http://www.soulfulbits.com/
    • Show only replies by krussh

Оффлайн Stefan Stefens

  • **
  • Сообщений: 98
  • Рейтинг: +4/-0
    • Show only replies by 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 - реальность. Или что-то в моих рассуждениях выпушено?
« Последнее редактирование: 26 Апр 2013 [19:42:41] от Stefan Stefens »
TLAPO804, ZEQ25GT, ZWO Finder/Guider, Orion Flattener, MGEN-2

Оффлайн 1212LupusАвтор темы

  • *****
  • Сообщений: 3 078
  • Рейтинг: +194/-4
  • Мне стал не очень нравиться этот форум...
    • Show only replies by 1212Lupus
    • http://belastro.net
Все параметры аппроксимируюшего звезное изображение Гаусиана посчитаны с довольно высокой точность (в рамках DLSR сьемки). Их хотелось бы дл я вычисления и применить. Конечно тоже сделано для звезд сравнения ( для база ).
Скажите, а фотомтерию треков на основе фотометрии опорных точечных звёзд, как ТУТ, нет желания сделать?
Радиоастрономы-любители -- объединяемся!


Если утро наступает в три -
Через два часа уже зажгут фонари.
Уже кончился день, а я только встал,
А я только что встал и уже устал.
(с) НОЛЬ

Оффлайн krussh

  • *****
  • Сообщений: 2 272
  • Рейтинг: +102/-5
  • http://www.soulfulbits.com/
    • Show only replies by krussh

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

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

Оффлайн Stefan Stefens

  • **
  • Сообщений: 98
  • Рейтинг: +4/-0
    • Show only replies by Stefan Stefens
Скажите, а фотомтерию треков на основе фотометрии опорных точечных звёзд, как ТУТ, нет желания сделать?
Нет. Так легко дойти до сейсмологии звезд. Время трека - соизмеримо со временами интенсивных искажений, которые в ссылке перечисленны. Поэтоми считаю задачу безнадежной, а вот звучашие звезды в ультразвуке имеют меньше естественных припятствий кроме АПЕРТУРЫ ( ФЭУ подходят) :)
TLAPO804, ZEQ25GT, ZWO Finder/Guider, Orion Flattener, MGEN-2

Оффлайн Stefan Stefens

  • **
  • Сообщений: 98
  • Рейтинг: +4/-0
    • Show only replies by Stefan Stefens
5 зв. величин, это отличие по потоку в 100 раз, шкала зв. величин логарифмическая.
Тогда очень важен выбор экспозиции для того чтобы попасть на мах чувствительность. Интересно. как автоматически подсказать - туда нельзя, а здесь можно?
TLAPO804, ZEQ25GT, ZWO Finder/Guider, Orion Flattener, MGEN-2