Друк

[Гість]

Ласкаво просимо, Гість. Будь ласка, увійдіть або зареєструйтеся.
27 Листопада 2024, 07:02:18

1 година 1 день 1 тиждень 1 місяць Назавжди

Увійти

[Новини]

•    Новини
•      Астрономія
•      Астрозаходи
•    Життя клуба
•      Новини
•      Проекти
•      Заходи
•      Звіти про спостереження
•      Документи
•    Статті
•      Обладнання та телескопи
•      Спостереження
•      Астрофотографія
•      Астрософт
•      Астрономія в Україні
•      Астрономи шуткують
•      Інше
•     
•      Файловий архів

Наш канал

Наші сайти

Астро-сайти

Програми

Наш канал

Наші сайти

Астро-сайти

Програми

[Олег Малий]

Для скачивания исходников, повторяем процедуру из Ответ #11 и т.д. выбирая в форме ... or select or enter one or more extensions:
формат RAW для скачивания.

[Олег Малий]

Для примера, помещаем 24 "сырых" фитса будущей мозаики в папку с адресом /media/astronom/astro/M51_raw/ (в вашем случае, адрес папки может быть другим).
Калибровочные файлы помещаются в папку с названием jref, которую мы можем поместить в папку с исходниками:

[Олег Малий]

Перед запуском обработки, указываем место нахождения папки jref с калибровочными файлами.
В моём случае, вот так я это делаю:

Код: [Select]

export jref=/media/astronom/astro/M51_raw/jref/
А теперь запускаем калибровку и смотрим отчёт о неудаче:

Код: [Select]

astronom@ubuntu:~$ export jref=/media/astronom/astro/M51_raw/jref/
astronom@ubuntu:~$ source activate iraf27

(iraf27) astronom@ubuntu:~$ pyraf
setting terminal type to xgterm...


   NOAO/IRAF PC-IRAF Revision 2.16 EXPORT Thu May 24 15:41:17 MST 2012
      This is the EXPORT version of IRAF V2.16 supporting PC systems.


  Welcome to IRAF.  To list the available commands, type ? or ??.  To get
  detailed information about a command, type `help <command>'.  To run  a
  command  or  load  a  package,  type  its name.   Type  `bye' to exit a
  package, or `logout' to get out  of the CL.    Type `news' to find  out
  what is new in the version of the system you are using. 

  Visit http://iraf.net if you have questions or to report problems.

  The following commands or packages are currently defined:

  (Updated on 2013-12-13)

clpackage/:
 adccdrom/      deitab/         language/       optic/          system/
 apropos        esowfi/         lists/          plot/           tables/
 cfh12k/        finder/         mem0/           proto/          ucsclris/
 cirred/        fitsutil/       mscdb/          rvsao/          upsqiid/
 clpackage/     gemini/         mscred/         softools/       user/
 ctio/          gmisc/          mtools/         song/           utilities/
 cutoutpkg/     guiapps/        nfextern/       sqiid/          vo/
 dataio/        images/         noao/           stecf/          xdimsum/
 dbms/          kepler/         obsolete/       stsdas/         xray/
PyRAF 2.2.dev Copyright (c) 2002 AURA
Python 2.7.12 Copyright (c) 2001-2016 Python Software Foundation.
Python/CL command line wrapper
  .help describes executive commands
--> cd /media/astronom/astro/M51_raw/
--> import acstools
/home/astronom/anaconda2/envs/iraf27/lib/python2.7/site-packages/matplotlib/font_manager.py:273: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.
  warnings.warn('Matplotlib is building the font cache using fc-list. This may take a moment.')
The following tasks in the acstools package can be run with TEAL:
          acs2d                  acs_destripe           acs_destripe_plus     
          acsccd                    acscte                    acsrej         
          acssum                    calacs
PixCteCorr is no longer supported. Please use acscte.
--> epar calacs

Task calacs is running...

CALACSBEG*** CALACS -- Version 8.3.3 (27-Jul-2015) ***
Begin    06-Jan-2017 21:19:07 WET
2017006211907-I--------------- CALACS started ----------------------------------
Initialized Association data ...

Input    j97c11kdq_raw.fits
GetAsnTable: incase = 1
Processing SINGLE image j97c11kdq_raw.fits...
SetInput: Determined what ASN table to process and how.
LoadAsn:  Processing SINGLE exposure
GetGlobalInfo: Ready to open primary header...
GetGlobalInfo: asn_table is j97c11kdq_raw.fits
Trying to open j97c11kdq_raw.fits...
Read in Primary header from j97c11kdq_raw.fits...
GetGlobalInfo: Read in header from Image
GetGlobalInfo: Detector and Instrument determined
LoadAsn:  Read in ASN table j97c11kdq_raw.fits

NUMBER of MEMBERS in TABLE: 1  PRODUCTS: 1  SUB-PRODUCTS: 1
Exposure-- Processing SINGLE Exposure j97c11kdq_raw.fits

CALACS: Detector ACS, type 1
CALACS: processing a CCD product
CALACS: processing CCD product 0
CALACS: processing posid = 1
GetSingle: Rootname: j97c11kdq, Output rootname: j97c11kdq_raw.fits

*** ERROR in HST I/O functions:
Filename jref$o151506fj_ccd.fits EXTNAME  EXTVER 0 CFITSIO status 104
failed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/o151506fj_ccd.fits[0]
Error opening image array.
Warning    CCDTAB `jref$o151506fj_ccd.fits' not found or can't open.

*** ERROR in HST I/O functions:
Filename jref$q860440tj_bpx.fits EXTNAME  EXTVER 0 CFITSIO status 104
failed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/o151506fj_ccd.fitsfailed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/q860440tj_bpx.fits[0]
Error opening image array.

*** ERROR in HST I/O functions:
Filename jref$p2k20392j_bia.fits EXTNAME  EXTVER 0 CFITSIO status 104
failed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/q860440tj_bpx.fitsfailed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/p2k20392j_bia.fits[0]
Error opening image array.
Warning    BIASFILE `jref$p2k20392j_bia.fits' not found or can't open.

*** ERROR in HST I/O functions:
Filename jref$xa81724cj_cte.fits EXTNAME  EXTVER 0 CFITSIO status 104
failed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/p2k20392j_bia.fitsfailed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/xa81724cj_cte.fits[0]
Error opening image array.
Warning    PCTETAB `jref$xa81724cj_cte.fits' not found or can't open.

*** ERROR in HST I/O functions:
Filename jref$p2k2039bj_drk.fits EXTNAME  EXTVER 0 CFITSIO status 104
failed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/xa81724cj_cte.fitsfailed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/p2k2039bj_drk.fits[0]
Error opening image array.
Warning    DARKFILE `jref$p2k2039bj_drk.fits' not found or can't open.

*** ERROR in HST I/O functions:
Filename jref$z5l1718qj_dkc.fits EXTNAME  EXTVER 0 CFITSIO status 104
failed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/p2k2039bj_drk.fitsfailed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/z5l1718qj_dkc.fits[0]
Error opening image array.
Warning    DRKCFILE `jref$z5l1718qj_dkc.fits' not found or can't open.

*** ERROR in HST I/O functions:
Filename jref$n6u12592j_pfl.fits EXTNAME  EXTVER 0 CFITSIO status 104
failed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/z5l1718qj_dkc.fitsfailed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/n6u12592j_pfl.fits[0]
Error opening image array.

*** ERROR in HST I/O functions:
Filename jref$lch1459bj_osc.fits EXTNAME  EXTVER 0 CFITSIO status 104
failed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/n6u12592j_pfl.fitsfailed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/lch1459bj_osc.fits[0]
Error opening image array.
Warning    OSCNTAB `jref$lch1459bj_osc.fits' not found or can't open.

*** ERROR in HST I/O functions:
Filename jref$08b18470j_imp.fits EXTNAME  EXTVER 0 CFITSIO status 104
failed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/lch1459bj_osc.fitsfailed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/08b18470j_imp.fits[0]
Error opening image array.
Warning    IMPHTTAB `jref$08b18470j_imp.fits' not found or can't open.
Warning    PFLTFILE, DFLTFILE, LFLTFILE are all blank
ERROR:    8 reference files were missing.
ERROR:    Couldn't process CCD data
ERROR:    CALACS processing NOT completed for j97c11kdq_raw.fits
ERROR:    CALACS processing NOT completed for j97c11kdq_raw.fits
ERROR:    Calibration file(s) missing.
ERROR:    Can't open file
-->
Причина неудачи в том, что программа не нашла калибровочных файлов в папке jref, где они должны были быть.
Большой плюс в том, что в отчёте программа нам указывает типы и названия нужных калибровочных файлов.

[Олег Малий]

Открываем в Интернете страницу HST Calibration Reference Data System (CRDS), на которой есть ссылки на нужные калибровочные файлы.

Программа сообщила, что не найден BIAS файл:

Код: [Select]

*** ERROR in HST I/O functions:
Filename jref$p2k20392j_bia.fits EXTNAME  EXTVER 0 CFITSIO status 104
failed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/q860440tj_bpx.fitsfailed to find or open the following file: (ffopen)/media/astronom/astro/M51_raw/jref/p2k20392j_bia.fits[0]
Error opening image array.
Warning    BIASFILE `jref$p2k20392j_bia.fits' not found or can't open.
Для его поиска и скачивания, открываем раздел biasfile ---- Bias Frame для камеры asc.
Подобным образом, мы находим скачиваем остальные калибровочные файлы, требуемые программой.

[Олег Малий]

Вот эти файлы для калибровки 24-х "сырых" фитсов  объекта M51 от камеры ACS/WFC, снятых через фильтр F814W.

[Олег Малий]

Вот как выглядит "сырой" фитс j97c11kdq_raw.fits из двух частей:

[Олег Малий]

Вид калиброванного j97c11kdq_flt.fits.
Также прикрепляю отчёт программы по калибровке этого исходника.

[Олег Малий]

После калибровки, необходимо внести изменения в хидеры, а также добавить специальные фитсы в каждый из 24-х нами обрабатываемых фитсов нашей будущей мозаики, согласно инструкции Hubble Space Telescope
ACS Distortion.

В нашем случем, для откалиброванного фитса j97c11kdq_flt.fits, мы сделаем вот так:

Код: [Select]

source activate iraf27
pyraf
cd /media/astronom/astro/M51_raw/flt
from astropy.io import fits
fits.setval('j97c11kdq_flt.fits', 'IDCTAB', value='/home/astronom/0461802ej_idc.fits')
fits.setval('j97c11kdq_flt.fits', 'NPOLFILE', value='/home/astronom/02c1450rj_npl.fits')
fits.setval('j97c11kdq_flt.fits', 'D2IMFILE', value='/home/astronom/02c1450rj_npl.fits')
Итак, изменения внесены, а теперь собственно обновляем хидер:

Код: [Select]

from stwcs import updatewcs
updatewcs.updatewcs('j97c11kdq_flt.fits')
Подобные действия проводим для всех нами откалиброванных фитсов.

В итоге, наш  j97c11kdq_flt.fits превращается в "монстрика" на 167,9Мб из 10 вложенных фитсов, необходимых для дальнейшего сложения мозаики:

[Олег Малий]

Теперь мы можем сложить эти 24 фитса в мозаику и получить такую красоту:

Подякували

1 :Konstantin

[OlegBr]

Бесконечное спасибо Олегу за найденные и подготовленные исходники 
Было очень интересно обрабатывать такой необычный материал.
За несколько дней сделал в Пиксе несколько вариантов комбинирования каналов.
К великому сожалению, чудом сохранился только один промежуточный вариант 
Во время очередного сохранения проекта в Пиксе я случайно обесточил компьютер 
Начинать все сначала уже нет ни времени, ни желания: я все свободное сейчас посвящаю обработке  моего собственного материала, отснятого в 2016 году.
Файл в полном разрешении здесь: http://olegbr.astroclub.kiev.ua/files/astrofoto/M51_HST_source/M51_HST_RHaGB_v3.jpg

Подякували

3 :Олег Малий, Jacques, serhio_67

[Олег Малий]

Бесконечное спасибо Олегу за найденные и подготовленные исходники 
Было очень интересно обрабатывать такой необычный материал.
За несколько дней сделал в Пиксе несколько вариантов комбинирования каналов.
К великому сожалению, чудом сохранился только один промежуточный вариант 
Во время очередного сохранения проекта в Пиксе я случайно обесточил компьютер 
Начинать все сначала уже нет ни времени, ни желания: я все свободное сейчас посвящаю обработке  моего собственного материала, отснятого в 2016 году.
Файл в полном разрешении здесь: http://olegbr.astroclub.kiev.ua/files/astrofoto/M51_HST_source/M51_HST_RHaGB_v3.jpg

Уникальность именно этого изображения в том, что оно есть полным с косыми краями, а не привычным прямоугольным кропом.

См похожую картинку от Роберта Гендлера http://www.robgendlerastropics.com/M51HST-GendlerSSS.html

[Oleg Astro]

Відновлю цю тему.
Отже, зараз всі ці мною описані процеси можна зібрати в один "блокнот", которий запускається в ізольованому середовищі.
Інакше кажучи, автоматична налагоджена лінія, котра прописана у "блокноті", сама все робить.

Кому із вас це цікаво?