Тематические рабочие группы > ПО и методики
VAST - Программа для поиска переменных на ПЗС-снимках
BrainBug:
Лучше поздно, чем никогда! ;D
Проблема была в следующем: так как управление ПСЗ-камерой QHY6 велась в CCDcap-е, то заголовок FITS-файлов записывался в НЕСТАНДАРТНОМ виде! Из-за этого MuniWin не хотел "кушать" такие FITS-ы...
Файлы, полученные в MaximDL-е MuniWin хавает на УРА просто! :)
Вывод: пользуйте ПРАВИЛЬНЫЕ программы для управления CCD! 8)
BrainBug:
--- Код: ---cerebrum vast-1.0rc33 # ./vast
--==VaST 1.0rc33==--
Detection of variable objects.
Usage: ./vast [options] image1.fit image2.fit ... imageN.fit
Options:
-h or --help print this message
-f or --nofind run ./find_candidates manually
-d or --debug run in debug (verbose) mode
-t or --time use ANOVA & BLS period search algorithms
-s or --small use small comparison window (1 pix, useful for most CCD images)
-m or --medium use medium comparison window (6 pix)
-w or --wide use wide comparison window (10 pix, usful for scanned photographic plates)
-p or --poly do not use polynomal magnitude calibration (useful for good quality CCD images)
-P use default.sex.PHOTO parameter file instead of default.sex (for photographic plates)
-r or --norotation assume no rotation greater then 3 deg.
-l or --filter use sigma filter (slower computation)
Examples:
./vast ../data/ccd_image-001.fit ../data/ccd_image-*.fit # Typical CCD image reduction.
./vast -s ../data/ccd_image-001.fit ../data/ccd_image-*.fit # CCD image reduction using fixed size of comparison window.
./vast -Pt ../data/photoplate-001.fit ../data/photoplate-*.fit # Typical scanned photographic plate reduction,
# configuration file default.sex.PHOTO and period search algorithms will be used.
cerebrum vast-1.0rc33 #
--- Конец кода ---
--- Код: ---Total number of measurements 0 (0 measurements stored in RAM)
writing lightcurve (outNNNN.dat) files...
MATCH_SUCESS=0
Writing summary file: vast_summary.log ...Ok
Images processed 100
Images used for photometry 1
Ref. image: 2454839.264456 07.01.2009 18:20:19 /mnt/loop1/AstroPhoto/qhy6/HEQ5ProSyntrek/Search_VS/Orion/RA_03h53m24s__DE_040629/Search_VS-001_R_.fit
First image: 2454839.264456 07.01.2009 18:20:19 /mnt/loop1/AstroPhoto/qhy6/HEQ5ProSyntrek/Search_VS/Orion/RA_03h53m24s__DE_040629/Search_VS-001_R_.fit
Last image: 2454839.341690 07.01.2009 20:11:32
--- Конец кода ---
Почему он не видит звёзд на кадре? :'( ???
kirx:
--- Цитата: BrainBug от 10 января 2009, 15:20:43 ---
--- Код: ---Total number of measurements 0 (0 measurements stored in RAM)
writing lightcurve (outNNNN.dat) files...
MATCH_SUCESS=0
Writing summary file: vast_summary.log ...Ok
Images processed 100
Images used for photometry 1
Ref. image: 2454839.264456 07.01.2009 18:20:19 /mnt/loop1/AstroPhoto/qhy6/HEQ5ProSyntrek/Search_VS/Orion/RA_03h53m24s__DE_040629/Search_VS-001_R_.fit
First image: 2454839.264456 07.01.2009 18:20:19 /mnt/loop1/AstroPhoto/qhy6/HEQ5ProSyntrek/Search_VS/Orion/RA_03h53m24s__DE_040629/Search_VS-001_R_.fit
Last image: 2454839.341690 07.01.2009 20:11:32
--- Конец кода ---
Почему он не видит звёзд на кадре? :'( ???
--- Конец цитаты ---
Да странно, почему-то не отождествился ни один кадр. Мне приходят в голову вот какие возможности:
* А сколько звёзд видно на кадре (нужно порядка сотни или больше)?
* Насколько согласуется число звёзд найденных VaST'ом (пишет по ходу обработки) с числом звёзд которые видно на кадре глазами? Не хватает ли васт много шума/горячих пикселей вместо звёзд? Если так, нужно попробовать поизменять значения
--- Код: ---DETECT_MINAREA 2.0 # minimum number of pixels above threshold
DETECT_THRESH 2.0 # <sigmas> or <threshold>,<ZP> in mag.arcsec-2
ANALYSIS_THRESH 2.0 # <sigmas> or <threshold>,<ZP> in mag.arcsec-2
--- Конец кода ---
в файле настроек SExtractor'а - default.sex в папке васта.
* Нет ли в Makefile (в папке васта) незакоментированной строки
--- Код: ---CFLAGS = -O3
--- Конец кода ---
? Благодаря Андрею Самохвалову выяснилось, что собранный с этим флагом васт начинает работать зверски нестабильно. Если флаг там - нужно убрать его и снова сделать make
* Ну и если программа так и не заработает - мне самому хотелось бы взглянуть на снимки! Плиз, загрузите их по фтп сюда ftp://scan.sai.msu.ru/uploads/
BrainBug:
Меня вот это беспокоит:
--- Код: ---Processing /mnt/loop1/AstroPhoto/qhy6/HEQ5ProSyntrek/VarSrars/XY_Leo/XY_Leo-098_R_.fit
Exposure 7 sec, 07.01.2009 00:22:31 UT = JD 2454838.515671 (mid. exp.)
sex -c autodetect_aper.sex /mnt/loop1/AstroPhoto/qhy6/HEQ5ProSyntrek/VarSrars/XY_Leo/XY_Leo-098_R_.fit
sex -PHOT_APERTURES 12.4 /mnt/loop1/AstroPhoto/qhy6/HEQ5ProSyntrek/VarSrars/XY_Leo/XY_Leo-098_R_.fit
----- SExtractor 2.5.0 started on 2009-01-10 at 22:42:30 with 1 thread
Measuring from: " " / 752 x 572 / 16 bits INTEGER data
(M+D) Background: 10137.4 RMS: 100.992 / Threshold: 201.983
Objects: detected 169 / sextracted 54
> All done (in 1 s)
No similar triangles
0 stars detected, 0 stars matched
--- Конец кода ---
Почему "0 stars detected, 0 stars matched" и кто такой "No similar triangles" ???
Вываливается после прогона всех файлов вот с такой кракозябриной:
--- Код: ---Total number of measurements 0 (0 measurements stored in RAM)
writing lightcurve (outNNNN.dat) files...
MATCH_SUCESS=0
Writing summary file: vast_summary.log ...Ok
Images processed 99
Images used for photometry 1
Ref. image: 2454838.472569 06.01.2009 23:20:27 /mnt/loop1/AstroPhoto/qhy6/HEQ5ProSyntrek/VarSrars/XY_Leo/XY_Leo-001_R_.fit
First image: 2454838.472569 06.01.2009 23:20:27 /mnt/loop1/AstroPhoto/qhy6/HEQ5ProSyntrek/VarSrars/XY_Leo/XY_Leo-001_R_.fit
Last image: 2454838.516123 07.01.2009 00:23:10
./find_candidates
Wiping out bad (nan) measurements... done! =)
ERROR: There are no light curve files!!!
ERROR: There are no light curve files!!!
check by util/data_parser -
util/nopgplot.sh: line 12: data.m_sigma: Нет такого файла или каталога
mv: невозможно выполнить stat для `tmp': Нет такого файла или каталога
ERROR: Can't open data.m_sigma!
done =)
ERROR: Can't open file "data.m_sigma"!
--- Конец кода ---
kirx:
Ситуация проясняется: действительно, проблема в том, что мало звёзд - поэтому васту не удаётся отождествить кадры. А SExtractor видимо не оптимально настроен на обработку данной комбинации ПЗС+телескоп. Он пишет, что найдено 169 звёзд, а обработано только 54:
--- Код: ---Objects: detected 169 / sextracted 54
--- Конец кода ---
Т.е. 2/3 обнаруженных звёзд секстрактор посчитал шумом и отбросил не передав дальше васту.
Да, прости, забыл написать в прошлом посте - параметры DETECT_THRESH и ANALYSIS_THRESH в default.sex нужно менять вместе. Если они отличаются - это может вызвать разницу между detected и sextracted. А может быть большая часть того, что хватает трактор и правда не звёзды а шум... Если на кадре правда только 54 звезды, то васт может вообще не справится с такими данными, ему для нормальной работы нужна хотя бы сотня. :(
No similar triangles означает, что не удалось отождествить кадр, отождествление осуществляется поиском подобных треугольников.
0 stars detected, 0 stars matched - не знаю, возможно это означает что секстрактор всем найденным звёздам проставляет плохие флаги (последняя колонка в файле test.cat, если не равна нулю - значит есть проблемы, васт такие звёзды не берёт), а может нули выдаются просто потому что отождествление не удалось.
Остальная ругань васта нормальна для ситуации, когда не удалось отождествить и обработать достаточное количество кадров.
Короче говоря, проблема в том, что секстрактор плохо находит звёзды. Это происходит либо потому, что их там просто нет (в нужном количестве), либо - потому, что настройки в default.sex плохо подходят для данного инструмента (телескоп+ПЗС).
Такие дела... ::)
Навигация
Перейти к полной версии