Photometry in generic data. Help.
27 views
Skip to first unread message
Astro Nomia
Mar 19, 2021, 3:47:00 PM3/19/21
to photometrypipeline_users
Greetings.
Very grateful for the availability of the beautiful code.
I have a good data set for the asteroid (5482) Korankei. There are 20 photos clicked with Appo TOA 150.
Based on the data, it was possible to calculate the orbit using find_orb as we can see in the prints.
I would like to photometry this set of photos to learn more about its pipeline.
I am still starting in Astronomy and I think there is a problem with the header of my .fits.
You can give me a light.
Thanks!
Cicero Lopes Neto
Obs. Const Sul Brazil.
Data:
## Info add Find_Orb:
COD XXX
COM Long. 142 48 22.7 E, Lat. 43 41 06.0 N, Alt. 1300m
ACK MPCReport file updated 2021.03.14 09:17:35
NET UCAC-4
123456789|123456789|123456789|123456789|123456789|123456789|123456789|123456789|
<-ObjDesig->*nnYYYY MM DD.DDDDD HH MM SS.SSSsdd mm ss.ss<blanks >MM.MMBz<ref>COD
05482 C2017 01 27.51714 03 52 43.04 +23 59 43.9 18.2 G XXX
05482 C2017 01 27.52307 03 52 43.07 +23 59 43.2 18.1 G XXX
05482 C2017 01 27.52901 03 52 43.13 +23 59 42.7 18.1 G XXX
05482 C2017 01 27.53495 03 52 43.16 +23 59 41.8 18.0 G XXX
05482 C2017 01 27.54089 03 52 43.21 +23 59 41.3 18.5 G XXX
05482 C2017 01 27.54682 03 52 43.31 +23 59 41.2 18.0 G XXX
05482 C2017 01 27.55514 03 52 43.35 +23 59 40.1 18.4 G XXX
05482 C2017 01 27.56108 03 52 43.41 +23 59 39.7 18.3 G XXX
05482 C2017 01 27.56702 03 52 43.44 +23 59 39.1 18.4 G XXX
05482 C2017 01 27.57297 03 52 43.50 +23 59 38.7 18.4 G XXX
05482 C2017 01 27.57891 03 52 43.52 +23 59 38.1 18.1 G XXX
05482 C2017 01 27.58484 03 52 43.60 +23 59 37.4 18.3 G XXX
05482 C2017 01 27.59078 03 52 43.65 +23 59 36.7 18.2 G XXX
05482 C2017 01 27.59672 03 52 43.70 +23 59 36.2 18.3 G XXX
05482 C2017 01 27.60267 03 52 43.74 +23 59 35.9 18.3 G XXX
05482 C2017 01 27.60895 03 52 43.78 +23 59 35.3 18.3 G XXX
05482 C2017 01 27.61489 03 52 43.85 +23 59 35.0 18.1 G XXX
05482 C2017 01 27.62083 03 52 43.89 +23 59 34.0 18.1 G XXX
05482 C2017 01 27.62677 03 52 43.95 +23 59 33.3 18.4 G XXX
05482 C2017 01 27.63271 03 52 43.98 +23 59 32.6 18.6 G XXX
#### Resume Find_Orb:
Orbital elements: (5482)
Perihelion 2015 Jul 6.612520914524794 +/- 309 TT = 14:42:01.8070149421 (JD 2457210.112520914500000)
Epoch 2017 Jan 28.0 TT = JDT 2457781.5 Earth MOID: 0.7158 Find_Orb
M 149.527193406700750 +/- 90 (J2000 ecliptic)
n 0.261691407109638900 +/- 0.311 Peri. 335.972319826813530 +/- 100
a 2.420711728275770200 +/- 1.92 Node 300.306503319154610 +/- 51
e 0.28484496431441025 +/- 0.243 Incl. 4.552527958607501 +/- 3.3
P 3.77 H 13.2 G 0.15 U 12.0 SR
q 1.731184182419584200 +/- 1.7 Q 3.110239274131955800 +/- 102
From 20 observations 2017 Jan. 27 (2.8 hr); mean residual 0".19
# State vector (heliocentric equatorial J2000):
# +0.610246165632 +2.674498735820 +1.335996038829 AU
# -8.015860844735 +2.649626362428 +0.666231659660 mAU/day
# MOIDs: Me 1.269501 Ve 1.003920 Ea 0.715808 Ma 0.274533
# MOIDs: Ju 2.027599 Sa 5.919215 Ur 15.528666 Ne 26.785658
# Elements written: 15 Mar 2021 16:11:44 (JD 2459289.174815)
# Full range of obs: 2017 Jan. 27 (2.8 hr) (20 observations)
# Find_Orb ver: Feb 14 2021 16:54:02
# Perturbers: 00000001 ; JPL DE-430
# Diameter 9436.7 meters (assuming 10% albedo)
# Score: -0.470735
# $Name=(5482) $Ty=2015 $Tm=07 $Td=06.612520 $MA=149.52719
# $ecc=0.2848450 $Eqnx=2000.
# $a=2.4207117 $Peri=335.97232 $Node=300.30650 $Incl=4.55253
# $EpJD=2457781.500 $q=1.731184 $T=2457210.112521 $H=13.2
# Sigmas avail: 3
Astrometry in Gaia DR2:
Field 1 did not solve (index index-5006-00.fits, field objects 1-10).
log-odds ratio 105.267 (5.21237e+45), 13 match, 0 conflict, 20 distractors, 17 index.
RA,Dec = (58.1764,23.9909), pixel scale 0.995211 arcsec/pix.
Hit/miss: +++++-+-+++-++------+-----------+(best)-------------------------------------------------------------------
Field 1: solved with index index-5006-00.fits.
Field 1 solved: writing to file C:/tycho-tracker/Data_Externo/01_x/ccdcieltmp.solved to indicate this.
Field: C:/tycho-tracker/Data_Externo/01_x/ccdcieltmp.fits
Field center: (RA,Dec) = (58.18, 23.99) deg.
Field center: (RA H:M:S, Dec D:M:S) = (03:52:42.314, +23:59:27.564).
Field size: 25.6581 x 17.4794 arcminutes
Below is an installation suggestion for Ubuntu 20.04 or Windows WSL:
Greetings.
Very grateful for the availability of the beautiful code.
I have a good data set for the asteroid (5482) Korankei. There are 20 photos clicked with Appo TOA 150.
Based on the data, it was possible to calculate the orbit using find_orb as we can see in the prints.
I would like to photometry this set of photos to learn more about its pipeline.
I am still starting in Astronomy and I think there is a problem with the header of my .fits.
You can give me a light.
Thanks!
Cicero Lopes Neto
Obs. Const Sul Brazil.
Data:
## Info add Find_Orb:
COD XXX
COM Long. 142 48 22.7 E, Lat. 43 41 06.0 N, Alt. 1300m
ACK MPCReport file updated 2021.03.14 09:17:35
NET UCAC-4
123456789|123456789|123456789|123456789|123456789|123456789|123456789|123456789|
<-ObjDesig->*nnYYYY MM DD.DDDDD HH MM SS.SSSsdd mm ss.ss<blanks >MM.MMBz<ref>COD
05482 C2017 01 27.51714 03 52 43.04 +23 59 43.9 18.2 G XXX
05482 C2017 01 27.52307 03 52 43.07 +23 59 43.2 18.1 G XXX
05482 C2017 01 27.52901 03 52 43.13 +23 59 42.7 18.1 G XXX
05482 C2017 01 27.53495 03 52 43.16 +23 59 41.8 18.0 G XXX
05482 C2017 01 27.54089 03 52 43.21 +23 59 41.3 18.5 G XXX
05482 C2017 01 27.54682 03 52 43.31 +23 59 41.2 18.0 G XXX
05482 C2017 01 27.55514 03 52 43.35 +23 59 40.1 18.4 G XXX
05482 C2017 01 27.56108 03 52 43.41 +23 59 39.7 18.3 G XXX
05482 C2017 01 27.56702 03 52 43.44 +23 59 39.1 18.4 G XXX
05482 C2017 01 27.57297 03 52 43.50 +23 59 38.7 18.4 G XXX
05482 C2017 01 27.57891 03 52 43.52 +23 59 38.1 18.1 G XXX
05482 C2017 01 27.58484 03 52 43.60 +23 59 37.4 18.3 G XXX
05482 C2017 01 27.59078 03 52 43.65 +23 59 36.7 18.2 G XXX
05482 C2017 01 27.59672 03 52 43.70 +23 59 36.2 18.3 G XXX
05482 C2017 01 27.60267 03 52 43.74 +23 59 35.9 18.3 G XXX
05482 C2017 01 27.60895 03 52 43.78 +23 59 35.3 18.3 G XXX
05482 C2017 01 27.61489 03 52 43.85 +23 59 35.0 18.1 G XXX
05482 C2017 01 27.62083 03 52 43.89 +23 59 34.0 18.1 G XXX
05482 C2017 01 27.62677 03 52 43.95 +23 59 33.3 18.4 G XXX
05482 C2017 01 27.63271 03 52 43.98 +23 59 32.6 18.6 G XXX
#### Resume Find_Orb:
Orbital elements: (5482)
Perihelion 2015 Jul 6.612520914524794 +/- 309 TT = 14:42:01.8070149421 (JD 2457210.112520914500000)
Epoch 2017 Jan 28.0 TT = JDT 2457781.5 Earth MOID: 0.7158 Find_Orb
M 149.527193406700750 +/- 90 (J2000 ecliptic)
n 0.261691407109638900 +/- 0.311 Peri. 335.972319826813530 +/- 100
a 2.420711728275770200 +/- 1.92 Node 300.306503319154610 +/- 51
e 0.28484496431441025 +/- 0.243 Incl. 4.552527958607501 +/- 3.3
P 3.77 H 13.2 G 0.15 U 12.0 SR
q 1.731184182419584200 +/- 1.7 Q 3.110239274131955800 +/- 102
From 20 observations 2017 Jan. 27 (2.8 hr); mean residual 0".19
# State vector (heliocentric equatorial J2000):
# +0.610246165632 +2.674498735820 +1.335996038829 AU
# -8.015860844735 +2.649626362428 +0.666231659660 mAU/day
# MOIDs: Me 1.269501 Ve 1.003920 Ea 0.715808 Ma 0.274533
# MOIDs: Ju 2.027599 Sa 5.919215 Ur 15.528666 Ne 26.785658
# Elements written: 15 Mar 2021 16:11:44 (JD 2459289.174815)
# Full range of obs: 2017 Jan. 27 (2.8 hr) (20 observations)
# Find_Orb ver: Feb 14 2021 16:54:02
# Perturbers: 00000001 ; JPL DE-430
# Diameter 9436.7 meters (assuming 10% albedo)
# Score: -0.470735
# $Name=(5482) $Ty=2015 $Tm=07 $Td=06.612520 $MA=149.52719
# $ecc=0.2848450 $Eqnx=2000.
# $a=2.4207117 $Peri=335.97232 $Node=300.30650 $Incl=4.55253
# $EpJD=2457781.500 $q=1.731184 $T=2457210.112521 $H=13.2
# Sigmas avail: 3
Astrometry in Gaia DR2:
Field 1 did not solve (index index-5006-00.fits, field objects 1-10).
log-odds ratio 105.267 (5.21237e+45), 13 match, 0 conflict, 20 distractors, 17 index.
RA,Dec = (58.1764,23.9909), pixel scale 0.995211 arcsec/pix.
Hit/miss: +++++-+-+++-++------+-----------+(best)-------------------------------------------------------------------
Field 1: solved with index index-5006-00.fits.
Field 1 solved: writing to file C:/tycho-tracker/Data_Externo/01_x/ccdcieltmp.solved to indicate this.
Field: C:/tycho-tracker/Data_Externo/01_x/ccdcieltmp.fits
Field center: (RA,Dec) = (58.18, 23.99) deg.
Field center: (RA H:M:S, Dec D:M:S) = (03:52:42.314, +23:59:27.564).
Field size: 25.6581 x 17.4794 arcminutes
Below is an installation suggestion for Ubuntu 20.04 or Windows WSL:
# Photometry Pipeline Michael Mommert
Support:
C:\linux\diversos_Docs\Astro\01_photometrypipeline_Michael_Mommert\00a_photometrypipeline-readthedocs-io-en-latest.pdf
# Windows Subsystem for Linux (WSL)
wsl -l -v
NAME STATE VERSION
* Debian Stopped 1
wsl.exe --set-version Ubuntu-20.04 1
wsl.exe --set-version Debian 1
wsl.exe --set-version Debian 2 # (2) old CONVERTION
# Installation and Setup
3.1.1 Installation
PP is available from github. You can get the source code by typing into your terminal:
# This will create a photometrypipeline/ directory in your current directory.
Software Requirements
PP only runs on Python 3. It furthermore requires requires git, a number of non-standard Python modules (available
from the Python Package Index through pip):
• numpy # https://numpy.org/
• scipy # https://www.scipy.org/
• astropy # https://www.astropy.org/
• astroquery (version >= 0.3.9) # https://github.com/astropy/astroquery
• matplotlib # https://matplotlib.org/
• future # http://python-future.org/
• skimage # https://scikit-image.org/
• pandas # https://pandas.pydata.org/
and some freely available software:
• imagemagick # https://imagemagick.org/index.php
• Source Extractor # http://www.astromatic.net/software/sextractor
• SCAMP (please download the latest development version) # http://www.astromatic.net/software/scamp
# Installation Instructions for Ubuntu 16.04+, Ubuntu 20.04, Debian 10 Buster
Clone the PP github repo:
cd ~/
# Install software requirements for SCAMP, Source Extractor and imagemagick:
sudo apt-get install -y build-essential libssl-dev libffi-dev git sextractor wget imagemagick curl libplplot-dev libshp-dev libcurl4-gnutls-dev liblapack3 liblapack-dev liblapacke liblapacke-dev libfftw3-3 libfftw3-dev libfftw3-single3 libatlas-base-dev scamp
# ok all
sudo apt install qpdfview # 0.4.18 para abrir as imagens.
sudo apt install -y python-numpy
sudo apt-get -y install python3-matplotlib
sudo apt-get -y install python3-scipy
sudo apt-get -y install python3-pandas
sudo apt-get -y install python3-future
sudo apt-get -y install python3-astropy
sudo apt update
sudo apt install -y python3-pip
python3 -m pip install -U pip
python3 -m pip install -U scikit-image
python3 -m pip install -U astroquery
# pip install --upgrade --user numpy scipy astropy astroquery matplotlib pandas future scikit-image
#Add these lines to the .bashrc file in your home directory and replace <path> with the actual path to the PP
## directory:
'''
# photometry pipeline setup
export PHOTPIPEDIR=<path>/photometrypipeline
export PATH=$PATH:~<path>/photometrypipeline/
'''
# photometry pipeline setup
export PHOTPIPEDIR=/home/teste/photometrypipeline
export PATH=$PATH:/home/teste/photometrypipeline
# export PATH=$PATH:/home/teste/.local/bin
# Kudos to towicode for figuring out the SCAMP requirements.
# Parabéns ao towicode por descobrir os requisitos do SCAMP.
## install in MAC...
# 3.1.2 Update your Version of PP
# In order to update your version of PP, simply change into photometrypipeline/ and type:
git pull
## You should do this regularly as PP is still under constant development.
# 3.1.3 Example Data
'''
The PP github clone comes with some sample data that can be used to test if the pipeline works properly. The data
were taken with the VATT4k camera on the VATT and can be found in example_data/vatt4k. In order to run
the pipeline on these images, copy them to a new directory, change there, and run pp_run mscience*fits. If
everything works out properly, the results (photometry_3552.dat) should resemble those in example_data/
vatt4k/LOG.
O clone PP do github vem com alguns dados de amostra que podem ser usados para testar se o pipeline funciona corretamente. Os dados
foram tiradas com a câmera VATT4k no VATT e podem ser encontradas em example_data / vatt4k. Para correr
o pipeline nessas imagens, copie-as para um novo diretório, mude lá e execute pp_run mscience * fits. Se
tudo funcionar corretamente, os resultados (photometry_3552.dat) devem ser semelhantes aos de example_data /
vatt4k / LOG.
'''
cp -r /home/teste/photometrypipeline/example_data/vatt4k /home/teste/science/photometry/vatt4k
cp -r /home/teste/photometrypipeline/example_data/vatt4k /home/teste/science/photometry/vatt4k_2
cp -r /home/teste/photometrypipeline/example_data/vatt4k /home/teste/science/photometry/vatt4k_3
cp -r /home/teste/photometrypipeline/example_data/vatt4k /home/teste/science/photometry/vatt4k_4
cp -r /home/teste/photometrypipeline/example_data/vatt4k /home/teste/science/photometry/vatt4k_5
cp -r /home/teste/photometrypipeline/example_data/vatt4k /home/teste/science/photometry/vatt4k_6
cd /home/teste/science/photometry/vatt4k_6
ls
LOG
mscience0217.fits
mscience0218.fits
cat LOG
VATT/VATT4k observations of asteroid 3552
the resulting photometry_3552.dat file should look like this:
# filename julian_dateast_mag ast_sig ast_ra ast_dec [1] [2] [3] [4] [5] ZP ZP_sig inst_mag in_sig [6] [7] [8] [9]
mscience0217.ldac 2457786.9566337 19.4248 0.0308 179.27990133 +0.42829624 -0.00 0.01 0.00 0.00 180.00 29.3772 0.0283 -9.9525 0.0122 SDSS-R9_transformed V 0 VATT4K
mscience0218.ldac 2457786.9590839 19.4644 0.0311 179.27962720 +0.42829370 -0.03 -0.01 0.00 0.00 180.00 29.5133 0.0290 -10.0489 0.0111 SDSS-R9_transformed V 0 VATT4K
#
# [1]: Horizons_RA - image_RA [arcsec]
# [2]: Horizons_DDec - image_Dec [arcsec]
# [3,4]: manual target offsets in RA and DEC [arcsec]
# [5]: exposure time (s)
# [6]: photometric catalog
# [7]: photometric band
# [8]: Source Extractor flag
# [9]: telescope/instrument
mv LOG LOG_OLD
pp_run mscience * fits
teste@teste:~/science/vatt4k_2$ pp_run mscience * fits
|==========================================| 1.3k/1.3k (100.00%) 0s
ERROR: cannot open file LOG
ERROR: cannot open file mscience
ERROR: cannot open file fits
run photometry pipeline on 2 VATT4K V frames
preparing mscience0217.fits
preparing mscience0218.fits
creating diagnostic output
----- run image registration
* extract sources from 2 frames
313 sources extracted from frame mscience0218.fits
282 sources extracted from frame mscience0217.fits
query Vizier for GAIA at 179.534/+0.332 in a 0.64 deg radius 1788 sources retrieved.
> WARNING: This executable has been compiled using a version of the ATLAS library without support for multithreading. Performance will be degraded.
----- SCAMP 2.0.4 started on 2021-03-19 at 06:57:49 with 4 threads
----- 2 inputs:
mscience0217.ldac: "3552 " no ext. header 1 set 218 detections
mscience0218.ldac: "3552 " no ext. header 1 set 233 detections
----- 451 detections loaded
----- 1 instrument found for astrometry:
Instrument A1 :
1 extensions
FILTER = 'TOP 2 BOT 1'
QRUNID =
----- 1 instrument found for photometry:
Instrument P1 :
FILTER = 'TOP 2 BOT 1'
----- 1 field group found:
Group 1: 2 fields at 11:58:07.99 +00:19:56.7 with radius 8.985'
instruments epoch center coordinates radius scale
mscience0217.ldac A1 P1 2017.1 11:58:07.99 +00:19:56.7 8.985' 0.3750"
mscience0218.ldac A1 P1 2017.1 11:58:07.99 +00:19:56.7 8.985' 0.3750"
----- Reference catalogs:
> Loading Catalog GAIA.cat...
> WARNING: FLAGS parameter not found in catalog GAIA.cat
654 astrometric references loaded from GAIA.cat
Group 1: 654 standards found in file (1 band)
----- Astrometric matching:
Group 1: 654 standards in file (band 1)
instruments pos.angle scale cont. shift cont.
mscience0217.ldac A1 P1 -1.66 deg 0.3749" 14.1 -7.9e+02"+3.3e+02" 9.67
mscience0218.ldac A1 P1 -1.66 deg 0.3750" 13.0 -8.0e+02"+3.3e+02" 9.36
----- Astrometric clipping:
Group 1: 36/414 detections removed
----- Astrometric stats (internal) :
All detections | High S/N
dAXIS1 dAXIS2 chi2 ndets | dAXIS1 dAXIS2 chi2 ndets
Group 1: 0.0178" 0.0177" 2.4e-07 268 0.0122" 0.0119" 1.1e-07 118
----- Astrometric stats (external):
All detections | High S/N
dAXIS1 dAXIS2 chi2 nstars | dAXIS1 dAXIS2 chi2 nstars
Group 1: 0.014" 0.0159" 2.1e-07 32 0.0142" 0.0158" 2e-07 31
----- Photometric stats (internal):
All detections | High S/N
Instru mag RMS chi2 ndets | mag RMS chi2 ndets
Group 1: P1 0.0713 1.3e+02 256 0.0403 2.8e+02 118
----- Photometric stats (external):
All detections | High S/N
Instru mag RMS chi2 nstars | mag RMS chi2 nstars
> All done (in 33 s: 0.06 field/s , 14 detections/s)
################################# REGISTRATION SUMMARY:
###
### 2/2 images have been registered successfully
###
######################################################
creating diagnostic output
----- derive optimum photometry aperture
* extract sources from 2 images using 20 apertures
601 sources extracted from frame mscience0217.fits
705 sources extracted from frame mscience0218.fits
processing curve-of-growth for frame mscience0217.fits
processing curve-of-growth for frame mscience0218.fits
creating diagnostic output
#################################### PHOTOMETRY SUMMARY:
###
### best-fit aperture radius 4.11 (px)
###
#####################################################
* extract sources from 2 images using aperture radius 4.11px
601 sources extracted from frame mscience0217.fits
705 sources extracted from frame mscience0218.fits
----- run photometric calibration
(432, 29) (sources, columns) read from mscience0217.fits
(499, 29) (sources, columns) read from mscience0218.fits
query Vizier for PANSTARRS at 179.312/+0.424 in a 0.14 deg radius 210 sources retrieved.
210 sources downloaded from PANSTARRS
PANSTARRS V
PANSTARRS transformed to V-band: 207 sources
zeropoint for mscience0217.ldac: 29.333+-0.052 (59/118 reference stars)
zeropoint for mscience0218.ldac: 29.472+-0.049 (59/119 reference stars)
write calibrated data into database files
----- distill photometry results
#------ Identify Targets
# check JPL Horizons for primary target...
3552 identified
# pick control star... done!
#-----------------------
2 potential target(s) per frame identified.
write photometry results for 3552_Don_Quixote_(1983_SA)
write photometry results for Control Star
extracting thumbnail images
Done!
teste@teste:~/science/vatt4k_2$
## OK!!!
# 3.3 Supported Observatories:
Telescope Instrument PP Keyword
Generic Telescope any GENERIC
3.2.4 Resultados
PP deriva a fotometria calibrada para o alvo que encontra na palavra-chave do cabeçalho OBJECT, bem como uma
'estrela de controle' brilhante que é usada para verificar a consistência da calibração fotométrica. Os resultados são gravados em arquivos
photometry_ <objectname> .dat no respectivo diretório de filtro.
Embora o PP seja projetado para ser executado quase automaticamente, algum bom senso é necessário para garantir que os resultados sejam
de confiança.
# The pipeline is currently set up to work on data from the following observatories/instruments:
Generic Telescope any GENERIC
(*): Esta câmera é um instrumento multi-detector; é recomendado dividir quadros FITS de multi-extensão a partir deste
instrumento em imagens FITS individuais de extensão única e para executar o pipeline nessas imagens FITS individuais.
(**): O cabeçalho da imagem SOAR Goodman tem várias palavras-chave não padrão que forçam o astropy a travar; em ordem
para usar PP, remova as palavras-chave do cabeçalho PARAM0, PARAM61, PARAM62 e PARAM63 antes de executar o PP, por exemplo,
usando o comando delhead fornecido pelo WCSTools.
(***): os wcs fornecidos com os dados do telescópio são geralmente excelentes; neste caso, sugiro pular o registro e usar
pp_run com a opção -keep_wcs.
Recomenda-se juntar imagens de câmeras com vários amplificadores usando as orientações corretas.
Além disso, é recomendável
pp_run mscience * fits # no exemplo
pp_run -keep_wcs mscience * fits
cd /home/teste/science/photometry/01x_fits
# Error:
pp_run mscience * fits
ERROR: cannot open file fits
Traceback (most recent call last):
File "/home/teste/photometrypipeline/pp_run", line 435, in <module>
rerun_registration, asteroids, keep_wcs)
File "/home/teste/photometrypipeline/pp_run", line 112, in run_the_pipeline
raise KeyError('cannot identify telescope/instrument; please update' +
KeyError: 'cannot identify telescope/instrument; please update_pp_conf.instrument_keys accordingly'
cd /home/teste/photometrypipeline/
cat -n pp_run |less
435 rerun_registration, asteroids, keep_wcs)
112 raise KeyError('cannot identify telescope/instrument; please update' +
113 '_pp_conf.instrument_keys accordingly')
pp_run -keep_wcs mscience * fits
# Ler cabeçalho fits:
mscience0217.fits
SIMPLE = T / conforms to FITS standard
BITPIX = -64 / array data type
NAXIS = 2 / number of array dimensions
NAXIS1 = 2032
NAXIS2 = 2032
EXTEND = T
NEXTEND = 2 / Number of extensions
DETSIZE = '[1:4064,1:4064]' / Detector size
CCDSUM = '2 2 ' / CCD pixel summing
CCDBIN1 = 2 / Binning factor along axis 1
CCDBIN2 = 2 / Binning factor along axis 2
FILENAME= 'science0217.fits' / base filename at acquisition
NCCDS = 1 / Number of CCDs
NAMPS = 2 / Number of amplifiers
COMMENT ==================================================================
COMMENT Controller
COMMENT ==================================================================
COMMENT Filter Top 5: B
DARKTIME= 180.135 / Dark time (seconds)
DATE = '2017-02-02' / UTC date and time file writtten
DATE-OBS= '2017-02-02' / UTC shutter opened
DETGAIN = 2 / Video gain setting
EXPREQ = 180.0 / Exposure time requested (seconds)
EXPTIME = 180.0 / PP: copied
IMAGETYP= 'object ' / Image type
INSTRUME= '= "Vatt4k"' / Instrument name
LOCTIME = '03:56:03' / Local time at start of exposure
OBJECT = '3552 '
OBSERVAT= '= "Vatican Observatory, Mt.Graham"' / Observatory
OBSERVER= '= "Mommert"' / Observer
PCIFILE = 'pci2.lod' / PCI board DSP code filename
TELESCOP= '= "VATT "' / Telescope
TIME-OBS= '10:56:03.149' / UTC at start of exposure
TIMESYS = 'UTC ' / Time system
TIMEZONE= 7 / Local time zone
TIMFILE = 'tim2.lod' / Timing board DSP code filename
UT = '10:56:03.149' / UTC at start of exposure
UTC-OBS = '10:56:03.149' / UTC at start of exposure
UTILFILE= 'util2.lod' / Utility board DSP code filename
VIDGAIN = 2 / Video gain setting
VIDSPEED= 2 / Video speed setting
COMMENT ==================================================================
COMMENT ITL Focal plane
COMMENT ==================================================================
ITL-HEAD= 'OK ' / ITL Header flag
N-AMPS-X= 2 / Number of amplifiers in X
N-AMPS-Y= 1 / Number of amplifiers in Y
N-DET-X = 1 / Number of detectors in X
N-DET-Y = 1 / Number of detectors in Y
NUM-AMPX= 2 / Number of amplifiers in X
NUM-AMPY= 1 / Number of amplifiers in Y
NUM-DETX= 1 / Number of detectors in X
NUM-DETY= 1 / Number of detectors in Y
REF-PIX1= 2032 / Reference pixel 1
REF-PIX2= 2032 / Reference pixel 2
SYS-NAME= 'Vatt4k ' / System name
COMMENT ==================================================================
COMMENT Instrument
COMMENT ==================================================================
COMMENT ==================================================================
COMMENT Telescope
COMMENT ==================================================================
AIRMASS = '1.19 ' / PP: copied
AZIMUTH = '193.7 ' / azimuth
DEC = '+00:19:56.5' / declination
ELEVAT = '57.1 ' / elevation
EPOCH = 2000 / equinox of RA and DEC
EQUINOX = '2000.0 ' / / Mean equinox
FILTER = 'TOP 2 BOT 1' / PP:copied
HA = '+00:29:36' / hour angle
JULIAN = '2457787.0' / julian date
LST-OBS = '12:27:44' / local siderial time
MOTION = '0 ' / motion flag
RA = '11:58:08.00' / right ascension
ROTANGLE= '0.0 ' / IIS rotation angle
ST = '12:27:44' / local siderial time
COMMENT ==================================================================
COMMENT Temperature
COMMENT ==================================================================
COMMENT (c) 2010-2018 Sorbon
COMMENT
CAMTEMP = -112.612 / Camera temperature in C
DEWTEMP = -195.81 / Dewar temperature in C
CCDSIZE = '[1:4064,1:4064]'
SECPIX = 0.375
RADECSYS= 'ICRS ' / copied from RADESYS
MIDTIMJD= 2457786.956633669 / PP: obs midtime
PHOTFLAG= 'F ' / PP: data is not photometric (SCAMP)
PHOT_K = 0.05 / PP: assumed extinction coefficient
SECPIXX = 0.375 / PP: x pixscale after binning
SECPIXY = 0.375 / PP: y pixscale after binning
TELINSTR= 'VATT/VATT4k' / PP: tel/instr name
TEL_KEYW= 'VATT4K ' / pipeline telescope keyword
RADESYS = 'ICRS ' / / Astrometric system
CTYPE1 = 'RA---TAN' / / WCS projection type for this axis
CTYPE2 = 'DEC--TAN' / / WCS projection type for this axis
CRVAL1 = '179.3125194277' / / World coordinate on this axis
CRVAL2 = '0.4249153774139' / / World coordinate on this axis
CRPIX1 = '1016.0 ' / / Reference pixel on this axis
CRPIX2 = '1016.0 ' / / Reference pixel on this axis
CD1_1 = '-0.0001041065570659' / / Linear projection matrix
CD1_2 = '-3.016194122478e-06' / / Linear projection matrix
CD2_1 = '-3.016194122477e-06' / / Linear projection matrix
CD2_2 = '0.0001041065570659' / / Linear projection matrix
CUNIT1 = 'deg ' / / Axis unit
CUNIT2 = 'deg ' / / Axis unit
PV1_0 = '1.173434865238e-05' / / Projection distortion parameter
PV1_1 = '0.9998827777594' / / Projection distortion parameter
PV1_2 = '0.0001094169028045' / / Projection distortion parameter
PV1_4 = '-0.002493883364172' / / Projection distortion parameter
PV1_5 = '0.000213515928354' / / Projection distortion parameter
PV1_6 = '-0.002015255140113' / / Projection distortion parameter
PV2_0 = '-1.888417253359e-05' / / Projection distortion parameter
PV2_1 = '0.9999629019103' / / Projection distortion parameter
PV2_2 = '-8.549018962087e-05' / / Projection distortion parameter
PV2_4 = '0.0008095629793874' / / Projection distortion parameter
PV2_5 = '0.0008532330639964' / / Projection distortion parameter
PV2_6 = '0.001388164148975' / / Projection distortion parameter
FGROUPNO= '1.0 ' / / SCAMP field group label
ASTIRMS1= '3.400238428119e-06' / / Astrom. dispersion RMS (intern., high S/N)
ASTIRMS2= '3.302501381191e-06' / / Astrom. dispersion RMS (intern., high S/N)
ASTRRMS1= '4.493482304298e-06' / / Astrom. dispersion RMS (ref., high S/N)
ASTRRMS2= '5.368961294597e-06' / / Astrom. dispersion RMS (ref., high S/N)
ASTINST = '1.0 ' / / SCAMP astrometric instrument label
FLXSCALE= '0.0 ' / / SCAMP relative flux scale
MAGZEROP= '0.0 ' / / SCAMP zero-point
PHOTIRMS= '0.04026301' / / mag dispersion RMS (internal, high S/N)
PHOTINST= '1.0 ' / / SCAMP photometric instrument label
PHOTLINK= 'F ' / / True if linked to a photometric field
REGCAT = 'GAIA ' / catalog used in WCS registration
PHOTMODE= 'APER ' / PP photometry mode
APRAD = 4.105263157894736 / aperture phot radius (px)
APIDX = 5 / optimum aprad index
HISTORY Astrometric solution
END
my fits: 01x_fits:
SIMPLE = T / file does conform to FITS standard
BITPIX = 16
NAXIS = 2 / number of data axes
NAXIS1 = 1546 / length of data axis 1
NAXIS2 = 1054 / length of data axis 2
EXTEND = T / FITS dataset may contain extensions
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H
BZERO = 3.3536024381605E+04
BSCALE = 9.6389964520067E-01
ISOSPEED= 800
EXPTIME = 487.649993896484 / Exposure time (in seconds)
EXPOSURE= 487.649993896484 / Exposure time (in seconds)
SOFTWARE= 'D '
DATE-OBS= '2017-01-27T20:28:30'
HISTORY CONVERTER IMAGEM RGB EM P&B
COMMENT Processado com Fitswork 05/03/2021
CTYPE1 = 'RA---TAN-SIP' / TAN (gnomic) projection + SIP distortions
CTYPE2 = 'DEC--TAN-SIP' / TAN (gnomic) projection + SIP distortions
EQUINOX = 2000.0 / Equatorial coordinates definition (yr)
LONPOLE = 180.0 / no comment
LATPOLE = 0.0 / no comment
CUNIT1 = 'deg ' / X pixel scale units
CUNIT2 = 'deg ' / Y pixel scale units
CRPIX1 = 774 / X reference pixel
CRPIX2 = 528 / Y reference pixel
CRVAL1 = 58.1761219178 / RA of reference point
CRVAL2 = 23.991196879 / DEC of reference point
CD1_1 = 2.17791775157E-07 / Transformation matrix
CD1_2 = -0.000276385613711 / no comment
CD2_1 = 0.000276576371102 / no comment
CD2_2 = 4.63278566891E-08 / no comment
A_ORDER = 3 / Polynomial order, axis 1
A_0_2 = -2.29110556876E-06 / no comment
A_1_1 = 2.27939343663E-06 / no comment
A_2_0 = -4.88687073143E-07 / no comment
B_ORDER = 3 / Polynomial order, axis 2
B_0_2 = -1.30117689593E-06 / no comment
B_1_1 = 5.41988070939E-07 / no comment
B_2_0 = 1.93947340231E-07 / no comment
AP_ORDER= 3 / Inv polynomial order, axis 1
AP_0_2 = 2.29842123072E-06 / no comment
AP_1_1 = -2.28171712278E-06 / no comment
AP_2_0 = 4.87872421697E-07 / no comment
BP_ORDER= 3 / Inv polynomial order, axis 2
BP_0_2 = 1.31098620931E-06 / no comment
BP_1_1 = -5.45445478928E-07 / no comment
BP_2_0 = -1.94783117368E-07 / no comment
A_1_2 = 5.55349883345E-09 / no comment
A_2_1 = -2.1237531711E-09 / no comment
B_1_2 = 3.73687497411E-09 / no comment
B_2_1 = -1.27254846514E-09 / no comment
AP_1_2 = -5.58076634657E-09 / no comment
AP_2_1 = 2.12956873917E-09 / no comment
BP_1_2 = -3.77446343593E-09 / no comment
BP_2_1 = 1.28204049743E-09 / no comment
A_0_3 = -4.79576487651E-09 / no comment
A_3_0 = 1.46865144669E-10 / no comment
B_0_3 = -8.60779049322E-09 / no comment
B_3_0 = 1.74376212235E-10 / no comment
AP_0_3 = 4.83650795331E-09 / no comment
AP_3_0 = -1.47368119505E-10 / no comment
BP_0_3 = 8.68605068521E-09 / no comment
BP_3_0 = -1.75541523953E-10 / no comment
END
O nome do catálogo entre colchetes é o identificador usado pelo PP; por exemplo, se você deseja usar URAT-1 para o registro de seu
imagens, use a opção -cat URAT-1 com o comando pp_register. O registro?" coluna identifica se este
O catálogo é apoiado pelo SCAMP a fim de usá-lo para registro de imagens. Todas as bandas de filtro marcadas com um asterisco (*)
são obtidos por meio de transformações de outras bandas; o respectivo sistema fotométrico é mostrado entre parênteses.
Se você estiver interessado em usar catálogos além dos listados, entre em contato.
Reply all
Reply to author
Forward
0 new messages