Author: Deborah Baines, Juan-Carlos Segovia, Raúl Gutierrez
This use case is partially based on the paper by Brown et. al. 2004, AJ 127, 2738 'RR Lyrae Stars in the Andromeda Halo from Deep Imaging with the Advanced Camera for Surveys'.
Workflow:
Additional:
# Import the relevant modules
from astroquery.utils.tap.core import TapPlus
import numpy as np
from astropy.stats import LombScargle
%matplotlib inline
import matplotlib.pyplot as plt
from matplotlib.colors import LogNorm
from astropy.io import fits
from astropy import visualization
# suppress warnings
import warnings
warnings.filterwarnings('ignore')
In this example, the RR Lyrae of interest is matchid = 106259089 (or V10 in Brown+ 2004).
# Call the HST TAP and print the available tables
hst = TapPlus(url="http://hst.esac.esa.int/tap-server/tap")
tables = hst.load_tables(only_names=True)
for table in (tables):
print(table.get_qualified_name())
Created TAP+ (v1.0.1) - Connection: Host: hst.esac.esa.int Use HTTPS: False Port: 80 SSL Port: 443 Retrieving tables... Parsing tables... Done. public.public.dual tap_schema.tap_schema.columns tap_schema.tap_schema.key_columns tap_schema.tap_schema.keys tap_schema.tap_schema.schemas tap_schema.tap_schema.tables hsc_v2.hsc_v2.hubble_sc2 hsc_v2.hsc_v2.hubble_sc2_ra_dec hsc.hsc.hubble_sc hsc.hsc.hubble_sc_ra_dec hcv.hcv.hcv hcv.hcv.hsc_aux hcv.hcv.mv_hcv_agg ehst.ehst.artifact ehst.ehst.energy ehst.ehst.instrument ehst.ehst.observation ehst.ehst.obs_simple_obs_composite ehst.ehst.plane ehst.ehst.position ehst.ehst.proposal ehst.ehst.publication ehst.ehst.publication_proposal ehst.ehst.target
# Inspect the columns of the HCV table
hst = TapPlus(url="http://hst.esac.esa.int/tap-server/tap")
hcv_table = hst.load_table('hcv.hcv')
for column in (hcv_table.columns):
print(column.name)
Created TAP+ (v1.0.1) - Connection: Host: hst.esac.esa.int Use HTTPS: False Port: 80 SSL Port: 443 Retrieving table 'hcv.hcv' Parsing table 'hcv.hcv'... Done. chi2 ci_d ci_v d_d dec d_v expert_class filter groupid lightcurve_cm lightcurve_d lightcurve_e lightcurve_i lightcurve_m lightcurve_r mad matchid pipeline_class ra
Query the HCV via an Astronomical Data Query Language (ADQL) query. See here for more information and examples of ADQL queries.
# Query the HCV for the RR Lyrae variable star with matchid = 106259089 and for the filter ACS_F814W:
job1 = hst.launch_job_async("SELECT * FROM hcv.hcv WHERE matchid = 106259089 AND filter LIKE '%ACS_F814W%'")
hcv_f814w = job1.get_results()
job1.get_data()
# Brings back a table of 33 rows.
Query finished.
chi2 | ci_d | ci_v | d_d | dec | d_v | expert_class | filter | groupid | lightcurve_cm | lightcurve_d | lightcurve_e | lightcurve_i | lightcurve_m | lightcurve_r | mad | matchid | pipeline_class | ra |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
days | days | deg | milliarcsec | mag | days | mag | mag | deg | ||||||||||
object | float64 | float64 | float64 | float64 | float64 | int32 | object | int32 | float64 | float64 | float64 | object | float64 | bool | float64 | int32 | int32 | float64 |
17.9332353549977981 | 52611.26638278016 | 1.0012962818145752 | 52611.26638278016 | 40.684165954589844 | 0.22348839044570923 | 1 | ACS_F814W | 1043756 | 25.4023624858617 | 52611.26638278016 | 0.052900001 | hst_9453_02_acs_wfc_f814w | 25.407801 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52611.57040482061 | 0.9920370578765869 | 52611.57040482061 | 40.684165954589844 | 0.029162408784031868 | 1 | ACS_F814W | 1043756 | 25.07910966437235 | 52611.57040482061 | 0.0403 | hst_9453_04_acs_wfc_f814w | 25.0867 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52612.09501705086 | 1.0814814567565918 | 52612.09501705086 | 40.684165954589844 | 2.103837013244629 | 1 | ACS_F814W | 1043756 | 25.119242325934767 | 52612.09501705086 | 0.041900001 | hst_9453_06_acs_wfc_f814w | 25.120701 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52614.096183951246 | 1.122592568397522 | 52614.096183951246 | 40.684165954589844 | 6.931029319763184 | 1 | ACS_F814W | 1043756 | 25.253824783127854 | 52614.096183951246 | 0.047699999 | hst_9453_08_acs_wfc_f814w | 25.253 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52614.7064349039 | 1.0025925636291504 | 52614.7064349039 | 40.684165954589844 | 3.0737504959106445 | 1 | ACS_F814W | 1043756 | 25.40671489103681 | 52614.7064349039 | 0.052299999 | hst_9453_10_acs_wfc_f814w | 25.4039 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52615.44050916098 | 1.1100925207138062 | 52615.44050916098 | 40.684165954589844 | 6.40604829788208 | 1 | ACS_F814W | 1043756 | 25.356878937073308 | 52615.44050916098 | 0.0526 | hst_9453_12_acs_wfc_f814w | 25.353399 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52616.09770231927 | 1.070462942123413 | 52616.09770231927 | 40.684165954589844 | 1.0694479942321777 | 1 | ACS_F814W | 1043756 | 25.264063289637335 | 52616.09770231927 | 0.047499999 | hst_9453_14_acs_wfc_f814w | 25.261 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52616.44162601838 | 1.0275925397872925 | 52616.44162601838 | 40.684165954589844 | 0.1264243721961975 | 1 | ACS_F814W | 1043756 | 25.2004714130637 | 52616.44162601838 | 0.044199999 | hst_9453_16_acs_wfc_f814w | 25.197001 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52618.30910287611 | 0.9604629278182983 | 52618.30910287611 | 40.684165954589844 | 1.1107144355773926 | 1 | ACS_F814W | 1043756 | 24.907989059070406 | 52618.30910287611 | 0.034299999 | hst_9453_20_acs_wfc_f814w | 24.905199 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
17.9332353549977981 | 52634.534542709356 | 1.0248147249221802 | 52634.534542709356 | 40.684165954589844 | 1.0793579816818237 | 1 | ACS_F814W | 1043756 | 25.29780057809714 | 52634.534542709356 | 0.048900001 | hst_9453_48_acs_wfc_f814w | 25.296 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52636.535665495554 | 0.9396296143531799 | 52636.535665495554 | 40.684165954589844 | 1.7726587057113647 | 1 | ACS_F814W | 1043756 | 25.34867465169749 | 52636.535665495554 | 0.051899999 | hst_9453_50_acs_wfc_f814w | 25.3444 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52639.632320828736 | 0.9510185122489929 | 52639.632320828736 | 40.684165954589844 | 1.160384178161621 | 1 | ACS_F814W | 1043756 | 24.907127633888535 | 52639.632320828736 | 0.034200002 | hst_9453_42_acs_wfc_f814w | 24.9067 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52644.20626751962 | 1.0006481409072876 | 52644.20626751962 | 40.684165954589844 | 1.3576633930206299 | 1 | ACS_F814W | 1043756 | 25.386741185391116 | 52644.20626751962 | 0.051899999 | hst_9453_54_acs_wfc_f814w | 25.3859 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52644.339669304434 | 1.0690741539001465 | 52644.339669304434 | 40.684165954589844 | 1.6359196901321411 | 1 | ACS_F814W | 1043756 | 25.06148505670546 | 52644.339669304434 | 0.041299999 | hst_9453_52_acs_wfc_f814w | 25.0674 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52646.140631039394 | 1.0164815187454224 | 52646.140631039394 | 40.684165954589844 | 1.0026894807815552 | 1 | ACS_F814W | 1043756 | 25.350860450582495 | 52646.140631039394 | 0.051100001 | hst_9453_56_acs_wfc_f814w | 25.354799 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52646.27403380303 | 0.9950925707817078 | 52646.27403380303 | 40.684165954589844 | 1.4151310920715332 | 1 | ACS_F814W | 1043756 | 25.377541952522026 | 52646.27403380303 | 0.052299999 | hst_9453_57_acs_wfc_f814w | 25.374701 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52646.407442505704 | 0.9368519186973572 | 52646.407442505704 | 40.684165954589844 | 3.7254488468170166 | 1 | ACS_F814W | 1043756 | 25.003620345310782 | 52646.407442505704 | 0.037700001 | hst_9453_58_acs_wfc_f814w | 24.999001 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52648.34181157872 | 0.9812037348747253 | 52648.34181157872 | 40.684165954589844 | 1.5127025842666626 | 1 | ACS_F814W | 1043756 | 25.434809714509754 | 52648.34181157872 | 0.054200001 | hst_9453_60_acs_wfc_f814w | 25.4345 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
17.9332353549977981 | 52650.40960101038 | 1.0606480836868286 | 52650.40960101038 | 40.684165954589844 | 3.3990793228149414 | 1 | ACS_F814W | 1043756 | 25.49112826041771 | 52650.40960101038 | 0.057 | hst_9453_59_acs_wfc_f814w | 25.488199 | False | 0.16130864386038724 | 106259089 | 2 | 11.502689361572266 |
# Query the HCV for the RR Lyrae variable star with matchid = 106259089 and for the filter ACS_F606W:
job2 = hst.launch_job_async("SELECT * FROM hcv.hcv WHERE matchid = 106259089 AND filter LIKE '%ACS_F606W%'")
hcv_f606w = job2.get_results()
job2.get_data()
# Brings back a table of 29 rows.
Query finished.
chi2 | ci_d | ci_v | d_d | dec | d_v | expert_class | filter | groupid | lightcurve_cm | lightcurve_d | lightcurve_e | lightcurve_i | lightcurve_m | lightcurve_r | mad | matchid | pipeline_class | ra |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
days | days | deg | milliarcsec | mag | days | mag | mag | deg | ||||||||||
object | float64 | float64 | float64 | float64 | float64 | int32 | object | int32 | float64 | float64 | float64 | object | float64 | bool | float64 | int32 | int32 | float64 |
60.0598271277244109 | 52610.161817079876 | 1.0532407760620117 | 52610.161817079876 | 40.684165954589844 | 1.876853346824646 | 1 | ACS_F606W | 1043756 | 25.194610885020445 | 52610.161817079876 | 0.028000001 | hst_9453_01_acs_wfc_f606w | 25.195101 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52611.43696170254 | 1.087592601776123 | 52611.43696170254 | 40.684165954589844 | 1.9586243629455566 | 1 | ACS_F606W | 1043756 | 24.949582840792328 | 52611.43696170254 | 0.023600001 | hst_9453_03_acs_wfc_f606w | 24.9569 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52611.703842436196 | 0.9923148155212402 | 52611.703842436196 | 40.684165954589844 | 0.8318588733673096 | 1 | ACS_F606W | 1043756 | 25.386888549877014 | 52611.703842436196 | 0.032900002 | hst_9453_05_acs_wfc_f606w | 25.395599 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52612.33606466558 | 0.9299073815345764 | 52612.33606466558 | 40.684165954589844 | 3.672166585922241 | 1 | ACS_F606W | 1043756 | 25.358675223412227 | 52612.33606466558 | 0.031500001 | hst_9453_07_acs_wfc_f606w | 25.356199 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52614.30502305669 | 0.9236111044883728 | 52614.30502305669 | 40.684165954589844 | 2.7837789058685303 | 1 | ACS_F606W | 1043756 | 25.122607184615944 | 52614.30502305669 | 0.0263 | hst_9453_09_acs_wfc_f606w | 25.119301 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52615.09689807985 | 1.0269443988800049 | 52615.09689807985 | 40.684165954589844 | 2.8672263622283936 | 1 | ACS_F606W | 1043756 | 25.336226176420528 | 52615.09689807985 | 0.031399999 | hst_9453_11_acs_wfc_f606w | 25.337601 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52615.707459393656 | 0.9747222661972046 | 52615.707459393656 | 40.684165954589844 | 3.2435302734375 | 1 | ACS_F606W | 1043756 | 25.23670991332121 | 52615.707459393656 | 0.028200001 | hst_9453_13_acs_wfc_f606w | 25.232201 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52616.306695474545 | 0.9244444370269775 | 52616.306695474545 | 40.684165954589844 | 2.2560760974884033 | 1 | ACS_F606W | 1043756 | 25.023412784196488 | 52616.306695474545 | 0.023800001 | hst_9453_15_acs_wfc_f606w | 25.0245 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52616.57511546649 | 1.010555624961853 | 52616.57511546649 | 40.684165954589844 | 1.6625279188156128 | 1 | ACS_F606W | 1043756 | 25.453671853616125 | 52616.57511546649 | 0.034400001 | hst_9453_17_acs_wfc_f606w | 25.449301 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
60.0598271277244109 | 52632.099195476854 | 1.0450925827026367 | 52632.099195476854 | 40.684165954589844 | 1.1953842639923096 | 1 | ACS_F606W | 1043756 | 25.03668623089168 | 52632.099195476854 | 0.0242 | hst_9453_39_acs_wfc_f606w | 25.038799 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52632.39987250231 | 0.9700926542282104 | 52632.39987250231 | 40.684165954589844 | 0.8609312772750854 | 1 | ACS_F606W | 1043756 | 25.439492511152125 | 52632.39987250231 | 0.033799998 | hst_9453_43_acs_wfc_f606w | 25.4391 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52632.53329266794 | 0.9382407069206238 | 52632.53329266794 | 40.684165954589844 | 3.5272231101989746 | 1 | ACS_F606W | 1043756 | 25.61165283094627 | 52632.53329266794 | 0.038899999 | hst_9453_45_acs_wfc_f606w | 25.609699 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52633.53394658561 | 0.9938889145851135 | 52633.53394658561 | 40.684165954589844 | 2.1130869388580322 | 1 | ACS_F606W | 1043756 | 25.12785408906986 | 52633.53394658561 | 0.0261 | hst_9453_47_acs_wfc_f606w | 25.129101 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52636.40226843394 | 0.9246296286582947 | 52636.40226843394 | 40.684165954589844 | 2.82668137550354 | 1 | ACS_F606W | 1043756 | 25.25209888468125 | 52636.40226843394 | 0.0287 | hst_9453_49_acs_wfc_f606w | 25.2523 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52636.7406423497 | 1.1949999332427979 | 52636.7406423497 | 40.684165954589844 | 4.103062629699707 | 1 | ACS_F606W | 1043756 | 25.445038014728738 | 52636.7406423497 | 0.037300002 | hst_9453_51_acs_wfc_f606w | 25.437799 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52640.40434069326 | 0.9737963080406189 | 52640.40434069326 | 40.684165954589844 | 2.4760165214538574 | 1 | ACS_F606W | 1043756 | 25.11754497263483 | 52640.40434069326 | 0.0263 | hst_9453_41_acs_wfc_f606w | 25.1166 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52644.473090491956 | 1.0048147439956665 | 52644.473090491956 | 40.684165954589844 | 0.4666833281517029 | 1 | ACS_F606W | 1043756 | 25.04464239639793 | 52644.473090491956 | 0.024499999 | hst_9453_53_acs_wfc_f606w | 25.040501 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52650.24182896991 | 1.008796215057373 | 52650.24182896991 | 40.684165954589844 | 3.2974252700805664 | 1 | ACS_F606W | 1043756 | 25.49795708409161 | 52650.24182896991 | 0.036200002 | hst_9453_55_acs_wfc_f606w | 25.496099 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 55204.29098817939 | 0.9488888382911682 | 55204.29098817939 | 40.684165954589844 | 4.4939727783203125 | 1 | ACS_F606W | 1043756 | 25.00825036271675 | 55204.29098817939 | 0.0337 | hst_11684_04_acs_wfc_f606w | 25.008699 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
# Perform exactly the same query as above, but this time using the astroquery.esa_hubble module
# Import ESAHubble from astroquery.esa_hubble
from astroquery.esa_hubble import ESAHubble
job = ESAHubble.query_hst_tap("SELECT * FROM hcv.hcv WHERE matchid = 106259089 AND filter LIKE '%ACS_F606W%'")
hcv_f606w = job.get_results()
job.get_data()
Created TAP+ (v1.0.1) - Connection: Host: hst.esac.esa.int Use HTTPS: False Port: 80 SSL Port: 443
chi2 | ci_d | ci_v | d_d | dec | d_v | expert_class | filter | groupid | lightcurve_cm | lightcurve_d | lightcurve_e | lightcurve_i | lightcurve_m | lightcurve_r | mad | matchid | pipeline_class | ra |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
days | days | deg | milliarcsec | mag | days | mag | mag | deg | ||||||||||
object | float64 | float64 | float64 | float64 | float64 | int32 | object | int32 | float64 | float64 | float64 | object | float64 | bool | float64 | int32 | int32 | float64 |
60.0598271277244109 | 52610.161817079876 | 1.0532407760620117 | 52610.161817079876 | 40.684165954589844 | 1.876853346824646 | 1 | ACS_F606W | 1043756 | 25.194610885020445 | 52610.161817079876 | 0.028000001 | hst_9453_01_acs_wfc_f606w | 25.195101 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52611.43696170254 | 1.087592601776123 | 52611.43696170254 | 40.684165954589844 | 1.9586243629455566 | 1 | ACS_F606W | 1043756 | 24.949582840792328 | 52611.43696170254 | 0.023600001 | hst_9453_03_acs_wfc_f606w | 24.9569 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52611.703842436196 | 0.9923148155212402 | 52611.703842436196 | 40.684165954589844 | 0.8318588733673096 | 1 | ACS_F606W | 1043756 | 25.386888549877014 | 52611.703842436196 | 0.032900002 | hst_9453_05_acs_wfc_f606w | 25.395599 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52612.33606466558 | 0.9299073815345764 | 52612.33606466558 | 40.684165954589844 | 3.672166585922241 | 1 | ACS_F606W | 1043756 | 25.358675223412227 | 52612.33606466558 | 0.031500001 | hst_9453_07_acs_wfc_f606w | 25.356199 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52614.30502305669 | 0.9236111044883728 | 52614.30502305669 | 40.684165954589844 | 2.7837789058685303 | 1 | ACS_F606W | 1043756 | 25.122607184615944 | 52614.30502305669 | 0.0263 | hst_9453_09_acs_wfc_f606w | 25.119301 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52615.09689807985 | 1.0269443988800049 | 52615.09689807985 | 40.684165954589844 | 2.8672263622283936 | 1 | ACS_F606W | 1043756 | 25.336226176420528 | 52615.09689807985 | 0.031399999 | hst_9453_11_acs_wfc_f606w | 25.337601 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52615.707459393656 | 0.9747222661972046 | 52615.707459393656 | 40.684165954589844 | 3.2435302734375 | 1 | ACS_F606W | 1043756 | 25.23670991332121 | 52615.707459393656 | 0.028200001 | hst_9453_13_acs_wfc_f606w | 25.232201 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52616.306695474545 | 0.9244444370269775 | 52616.306695474545 | 40.684165954589844 | 2.2560760974884033 | 1 | ACS_F606W | 1043756 | 25.023412784196488 | 52616.306695474545 | 0.023800001 | hst_9453_15_acs_wfc_f606w | 25.0245 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52616.57511546649 | 1.010555624961853 | 52616.57511546649 | 40.684165954589844 | 1.6625279188156128 | 1 | ACS_F606W | 1043756 | 25.453671853616125 | 52616.57511546649 | 0.034400001 | hst_9453_17_acs_wfc_f606w | 25.449301 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
60.0598271277244109 | 52632.099195476854 | 1.0450925827026367 | 52632.099195476854 | 40.684165954589844 | 1.1953842639923096 | 1 | ACS_F606W | 1043756 | 25.03668623089168 | 52632.099195476854 | 0.0242 | hst_9453_39_acs_wfc_f606w | 25.038799 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52632.39987250231 | 0.9700926542282104 | 52632.39987250231 | 40.684165954589844 | 0.8609312772750854 | 1 | ACS_F606W | 1043756 | 25.439492511152125 | 52632.39987250231 | 0.033799998 | hst_9453_43_acs_wfc_f606w | 25.4391 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52632.53329266794 | 0.9382407069206238 | 52632.53329266794 | 40.684165954589844 | 3.5272231101989746 | 1 | ACS_F606W | 1043756 | 25.61165283094627 | 52632.53329266794 | 0.038899999 | hst_9453_45_acs_wfc_f606w | 25.609699 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52633.53394658561 | 0.9938889145851135 | 52633.53394658561 | 40.684165954589844 | 2.1130869388580322 | 1 | ACS_F606W | 1043756 | 25.12785408906986 | 52633.53394658561 | 0.0261 | hst_9453_47_acs_wfc_f606w | 25.129101 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52636.40226843394 | 0.9246296286582947 | 52636.40226843394 | 40.684165954589844 | 2.82668137550354 | 1 | ACS_F606W | 1043756 | 25.25209888468125 | 52636.40226843394 | 0.0287 | hst_9453_49_acs_wfc_f606w | 25.2523 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52636.7406423497 | 1.1949999332427979 | 52636.7406423497 | 40.684165954589844 | 4.103062629699707 | 1 | ACS_F606W | 1043756 | 25.445038014728738 | 52636.7406423497 | 0.037300002 | hst_9453_51_acs_wfc_f606w | 25.437799 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52640.40434069326 | 0.9737963080406189 | 52640.40434069326 | 40.684165954589844 | 2.4760165214538574 | 1 | ACS_F606W | 1043756 | 25.11754497263483 | 52640.40434069326 | 0.0263 | hst_9453_41_acs_wfc_f606w | 25.1166 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52644.473090491956 | 1.0048147439956665 | 52644.473090491956 | 40.684165954589844 | 0.4666833281517029 | 1 | ACS_F606W | 1043756 | 25.04464239639793 | 52644.473090491956 | 0.024499999 | hst_9453_53_acs_wfc_f606w | 25.040501 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 52650.24182896991 | 1.008796215057373 | 52650.24182896991 | 40.684165954589844 | 3.2974252700805664 | 1 | ACS_F606W | 1043756 | 25.49795708409161 | 52650.24182896991 | 0.036200002 | hst_9453_55_acs_wfc_f606w | 25.496099 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
60.0598271277244109 | 55204.29098817939 | 0.9488888382911682 | 55204.29098817939 | 40.684165954589844 | 4.4939727783203125 | 1 | ACS_F606W | 1043756 | 25.00825036271675 | 55204.29098817939 | 0.0337 | hst_11684_04_acs_wfc_f606w | 25.008699 | False | 0.16173090767108178 | 106259089 | 2 | 11.502689361572266 |
#Plotting the light curves (note, ignoring the one F606W value at MJD = 55204.29)
x = hcv_f814w['lightcurve_d']
y = hcv_f814w['lightcurve_cm']
yerr = hcv_f814w['lightcurve_e']
x2 = hcv_f606w['lightcurve_d']
y2 = hcv_f606w['lightcurve_cm']
yerr2 = hcv_f606w['lightcurve_e']
plt.figure(figsize=(8, 8))
plt.scatter(x, y, color='g', alpha=0.5)
plt.errorbar(x, y, yerr, fmt='o', color='g', alpha=0.5)
plt.scatter(x2, y2, color='b', alpha=0.5)
plt.errorbar(x2, y2, yerr2, fmt='o', color='b', alpha=0.5)
plt.xlim(52600.0, 52660.0)
plt.ylim(25.7, 24.7) # flip the y axis
plt.xlabel('MJD (days)')
plt.ylabel('Corrected Magnitude')
plt.legend(["F814W", "F606W"])
plt.show()
Generate a periodgram using the Astropy 'LombScargle' function and obtain the period based on the best frequency. More information can be found here.
# For the ACS_F814W filter:
hcv_f814w['lightcurve_d'].unit = None
hcv_f814w['lightcurve_cm'].unit = None
hcv_f814w['lightcurve_e'].unit = None
#generate periodogram
frequency, power = LombScargle(hcv_f814w['lightcurve_d'], \
hcv_f814w['lightcurve_cm'], \
hcv_f814w['lightcurve_e']).autopower(minimum_frequency=1.0,\
maximum_frequency=5.0)
best_frequency = frequency[np.argmax(power)]
best_power = power[np.argmax(power)]
plt.figure(figsize=(8, 8))
plt.scatter([best_frequency], [best_power], facecolors='none', edgecolors='r', s=80)
plt.plot(frequency, power)
plt.xlabel('Frequency')
plt.ylabel('Relative Power')
#print(best_frequency)
period=1/best_frequency
print("Best found period: " + str(period) + " days")
Best found period: 0.6874079029387833 days
This result compares well with Brown+ 2004, who found a period of 0.687 days.
phase = (hcv_f814w['lightcurve_d'] / period) % 1
plt.figure(figsize=(8, 8))
plt.axes().invert_yaxis()
plt.scatter(phase, hcv_f814w['lightcurve_cm'], color='b', alpha=0.5)
plt.errorbar(phase, hcv_f814w['lightcurve_cm'], hcv_f814w['lightcurve_e'], fmt='o', color='b', alpha=0.5)
plt.scatter(phase+1, hcv_f814w['lightcurve_cm'], color='b', alpha=0.5)
plt.errorbar(phase+1, hcv_f814w['lightcurve_cm'], hcv_f814w['lightcurve_e'], fmt='o', color='b', alpha=0.5)
plt.xlim(0.0, 2.0)
#plt.ylim(25.7, 24.7) # flip the y axis
plt.xlabel('Phase')
plt.ylabel('Corrected Magnitude')
plt.show()
Check the results by increasing and decreasing the period:
from ipywidgets import widgets
from ipywidgets import interact
def fn(x):
print("period=" + str(x))
phase = (hcv_f814w['lightcurve_d'] / x) % 1
#print(phase)
plt.figure(figsize=(8, 8))
plt.axes().invert_yaxis()
plt.scatter(phase, hcv_f814w['lightcurve_cm'], color='b', alpha=0.3)
plt.errorbar(phase, hcv_f814w['lightcurve_cm'], hcv_f814w['lightcurve_e'], fmt='o', color='b', alpha=0.3)
plt.xlabel('Phase')
plt.ylabel('Corrected Magnitude')
plt.show()
#interact(fn, x=period)
#interact(fn, x=period)
max = (period + period*0.01)
min = (period - period*0.01)
interact(fn, x=widgets.FloatSlider(min=min,max=max,step=0.00005,value=period));
hst = TapPlus(url="http://hst.esac.esa.int/tap-server/tap")
# Performing a cone search of radius 3 arcminutes:
job3 = hst.launch_job_async("SELECT * FROM hcv.hcv \
WHERE 1=CONTAINS( \
POINT('ICRS', ra, dec), \
CIRCLE('ICRS', 11.50269, +40.68417, 0.05))")
hcv = job3.get_results()
job3.get_data()
# Returns 3335 rows. Note the number of rows is larger than the number of sources since the filter information
# and magnitudes (light curves) are all given in the one (flattened) table.
Created TAP+ (v1.0.1) - Connection: Host: hst.esac.esa.int Use HTTPS: False Port: 80 SSL Port: 443 Query finished.
chi2 | ci_d | ci_v | d_d | dec | d_v | expert_class | filter | groupid | lightcurve_cm | lightcurve_d | lightcurve_e | lightcurve_i | lightcurve_m | lightcurve_r | mad | matchid | pipeline_class | ra |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
days | days | deg | milliarcsec | mag | days | mag | mag | deg | ||||||||||
object | float64 | float64 | float64 | float64 | float64 | int32 | object | int32 | float64 | float64 | float64 | object | float64 | bool | float64 | int32 | int32 | float64 |
3.09628240816075495 | 52611.26638278016 | 1.059999942779541 | 52611.26638278016 | 40.70929718017578 | 4.754214763641357 | 2 | ACS_F814W | 1043756 | 24.553740054403757 | 52611.26638278016 | 0.0266 | hst_9453_02_acs_wfc_f814w | 24.5581 | False | 0.03035902804316848 | 98239091 | 1 | 11.494104385375977 |
3.09628240816075495 | 52614.7064349039 | 1.1368519067764282 | 52614.7064349039 | 40.70929718017578 | 3.1014959812164307 | 2 | ACS_F814W | 1043756 | 24.454888789590164 | 52614.7064349039 | 0.0243 | hst_9453_10_acs_wfc_f814w | 24.451599 | False | 0.03035902804316848 | 98239091 | 1 | 11.494104385375977 |
3.09628240816075495 | 52616.09770231927 | 1.1049073934555054 | 52616.09770231927 | 40.70929718017578 | 7.896210670471191 | 2 | ACS_F814W | 1043756 | 24.47698945440076 | 52616.09770231927 | 0.0272 | hst_9453_14_acs_wfc_f814w | 24.477301 | False | 0.03035902804316848 | 98239091 | 1 | 11.494104385375977 |
3.09628240816075495 | 52626.490075099515 | 1.1442593336105347 | 52626.490075099515 | 40.70929718017578 | 4.49559211730957 | 2 | ACS_F814W | 1043756 | 24.50593209383585 | 52626.490075099515 | 0.0363 | hst_9453_61_acs_wfc_f814w | 24.4995 | False | 0.03035902804316848 | 98239091 | 1 | 11.494104385375977 |
3.09628240816075495 | 52631.39916067012 | 1.074814796447754 | 52631.39916067012 | 40.70929718017578 | 7.790589809417725 | 2 | ACS_F814W | 1043756 | 24.507348482443927 | 52631.39916067012 | 0.026000001 | hst_9453_38_acs_wfc_f814w | 24.5056 | False | 0.03035902804316848 | 98239091 | 1 | 11.494104385375977 |
3.09628240816075495 | 52644.20626751962 | 0.9599074125289917 | 52644.20626751962 | 40.70929718017578 | 7.809937477111816 | 2 | ACS_F814W | 1043756 | 24.59013931092299 | 52644.20626751962 | 0.0265 | hst_9453_54_acs_wfc_f814w | 24.590799 | False | 0.03035902804316848 | 98239091 | 1 | 11.494104385375977 |
3.09628240816075495 | 52646.27403380303 | 1.060092568397522 | 52646.27403380303 | 40.70929718017578 | 8.634530067443848 | 2 | ACS_F814W | 1043756 | 24.526681052691977 | 52646.27403380303 | 0.026000001 | hst_9453_57_acs_wfc_f814w | 24.5263 | False | 0.03035902804316848 | 98239091 | 1 | 11.494104385375977 |
5.55475441578078222 | 52611.26638278016 | 0.9527778029441833 | 52611.26638278016 | 40.7080078125 | 0.7770352363586426 | 2 | ACS_F814W | 1043756 | 23.584190315361923 | 52611.26638278016 | 0.012 | hst_9453_02_acs_wfc_f814w | 23.589199 | False | 0.021911071428846185 | 85054530 | 1 | 11.492891311645508 |
5.55475441578078222 | 52611.57040482061 | 0.9789814352989197 | 52611.57040482061 | 40.7080078125 | 0.48111963272094727 | 2 | ACS_F814W | 1043756 | 23.577598627816684 | 52611.57040482061 | 0.012 | hst_9453_04_acs_wfc_f814w | 23.582899 | False | 0.021911071428846185 | 85054530 | 1 | 11.492891311645508 |
... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
10.5620890982210387 | 52632.099195476854 | 1.0135184526443481 | 52632.099195476854 | 40.72599792480469 | 3.0463757514953613 | 1 | ACS_F606W | 1043756 | 25.912892871453717 | 52632.099195476854 | 0.0528 | hst_9453_39_acs_wfc_f606w | 25.9156 | False | 0.09967318683652593 | 58227924 | 1 | 11.4900484085083 |
10.5620890982210387 | 52632.39987250231 | 0.9375925660133362 | 52632.39987250231 | 40.72599792480469 | 4.641005992889404 | 1 | ACS_F606W | 1043756 | 25.906379075083034 | 52632.39987250231 | 0.052000001 | hst_9453_43_acs_wfc_f606w | 25.907301 | False | 0.09967318683652593 | 58227924 | 1 | 11.4900484085083 |
10.5620890982210387 | 52632.53329266794 | 1.0228703022003174 | 52632.53329266794 | 40.72599792480469 | 2.021033763885498 | 1 | ACS_F606W | 1043756 | 25.55147467856299 | 52632.53329266794 | 0.039000001 | hst_9453_45_acs_wfc_f606w | 25.5527 | False | 0.09967318683652593 | 58227924 | 1 | 11.4900484085083 |
10.5620890982210387 | 52633.53394658561 | 0.9825000166893005 | 52633.53394658561 | 40.72599792480469 | 5.543918609619141 | 1 | ACS_F606W | 1043756 | 25.906709789001926 | 52633.53394658561 | 0.053399999 | hst_9453_47_acs_wfc_f606w | 25.908199 | False | 0.09967318683652593 | 58227924 | 1 | 11.4900484085083 |
10.5620890982210387 | 52636.40226843394 | 0.9477777481079102 | 52636.40226843394 | 40.72599792480469 | 3.5023529529571533 | 1 | ACS_F606W | 1043756 | 25.896994696105338 | 52636.40226843394 | 0.0524 | hst_9453_49_acs_wfc_f606w | 25.895599 | False | 0.09967318683652593 | 58227924 | 1 | 11.4900484085083 |
10.5620890982210387 | 52636.7406423497 | 0.9320370554924011 | 52636.7406423497 | 40.72599792480469 | 1.2833119630813599 | 1 | ACS_F606W | 1043756 | 25.714898499205404 | 52636.7406423497 | 0.044 | hst_9453_51_acs_wfc_f606w | 25.7113 | False | 0.09967318683652593 | 58227924 | 1 | 11.4900484085083 |
10.5620890982210387 | 52640.40434069326 | 1.1238888502120972 | 52640.40434069326 | 40.72599792480469 | 1.8247007131576538 | 1 | ACS_F606W | 1043756 | 25.905899646319572 | 52640.40434069326 | 0.0526 | hst_9453_41_acs_wfc_f606w | 25.9074 | False | 0.09967318683652593 | 58227924 | 1 | 11.4900484085083 |
10.5620890982210387 | 52644.473090491956 | 1.0012037754058838 | 52644.473090491956 | 40.72599792480469 | 1.8107696771621704 | 1 | ACS_F606W | 1043756 | 25.790836267475434 | 52644.473090491956 | 0.047400001 | hst_9453_53_acs_wfc_f606w | 25.7946 | False | 0.09967318683652593 | 58227924 | 1 | 11.4900484085083 |
10.5620890982210387 | 52650.24182896991 | 1.0362036228179932 | 52650.24182896991 | 40.72599792480469 | 4.1564555168151855 | 1 | ACS_F606W | 1043756 | 25.630031976721146 | 52650.24182896991 | 0.041999999 | hst_9453_55_acs_wfc_f606w | 25.628099 | False | 0.09967318683652593 | 58227924 | 1 | 11.4900484085083 |
10.5620890982210387 | 55204.29098817939 | 0.9985185265541077 | 55204.29098817939 | 40.72599792480469 | 5.1229472160339355 | 1 | ACS_F606W | 1043756 | 25.920661378733925 | 55204.29098817939 | 0.078000002 | hst_11684_04_acs_wfc_f606w | 25.920601 | False | 0.09967318683652593 | 58227924 | 1 | 11.4900484085083 |
# The associated HLA images are found in the HCV column lightcurve_i:
print (hcv_f814w['lightcurve_i'])
lightcurve_i ------------------------- hst_9453_02_acs_wfc_f814w hst_9453_04_acs_wfc_f814w hst_9453_06_acs_wfc_f814w hst_9453_08_acs_wfc_f814w hst_9453_10_acs_wfc_f814w hst_9453_12_acs_wfc_f814w hst_9453_14_acs_wfc_f814w hst_9453_16_acs_wfc_f814w hst_9453_20_acs_wfc_f814w hst_9453_22_acs_wfc_f814w ... hst_9453_46_acs_wfc_f814w hst_9453_48_acs_wfc_f814w hst_9453_50_acs_wfc_f814w hst_9453_42_acs_wfc_f814w hst_9453_54_acs_wfc_f814w hst_9453_52_acs_wfc_f814w hst_9453_56_acs_wfc_f814w hst_9453_57_acs_wfc_f814w hst_9453_58_acs_wfc_f814w hst_9453_60_acs_wfc_f814w hst_9453_59_acs_wfc_f814w Length = 33 rows
# Download a single image from the ESA Hubble Science Archive using the astroquery.esa_hubble module.
# A single observation has the prefix 'drz' (drizzle). Therefore the files to download are:
# hst_9453_**_acs_wfc_f814w_drz.fits, where ** is a number
from astroquery.esa_hubble import ESAHubble
ESAHubble.get_artifact("hst_9453_02_acs_wfc_f814w_drz.fits")
image_file = "hst_9453_02_acs_wfc_f814w_drz.fits"
http://archives.esac.esa.int/ehst-sl-server/servlet/data-action?ARTIFACT_ID=hst_9453_02_acs_wfc_f814w_drz.fits File hst_9453_02_acs_wfc_f814w_drz.fits downloaded to current directory
# Check the headers
hdu_list = fits.open(image_file)
hdu_list.info()
hdu_list[0].header
Filename: hst_9453_02_acs_wfc_f814w_drz.fits No. Name Ver Type Cards Dimensions Format 0 PRIMARY 1 PrimaryHDU 861 () 1 SCI 1 ImageHDU 86 (5185, 5166) float32 2 WHT 1 ImageHDU 46 (5185, 5166) float32 3 CTX 1 ImageHDU 39 (5185, 5166) int32 4 HDRTAB 1 BinTableHDU 631 8R x 311C [9A, 3A, J, D, D, D, D, D, D, D, D, D, D, D, D, D, J, 3A, 9A, 7A, 18A, 4A, D, D, D, D, 3A, D, D, D, D, D, D, D, D, D, D, D, D, J, 8A, 23A, D, D, D, D, J, J, J, 8A, J, 23A, 9A, 19A, J, 4A, J, J, J, J, J, J, 23A, D, D, D, D, J, J, 3A, 3A, 4A, 4A, J, D, D, D, 3A, 1A, J, D, D, D, 13A, 3A, 4A, 4A, 12A, 12A, 23A, 8A, 23A, 10A, 10A, D, D, 3A, 3A, 23A, 4A, 8A, 7A, 23A, D, J, D, 6A, 9A, 8A, D, D, J, 4A, 18A, 3A, J, 7A, 5A, 3A, D, 13A, 8A, 4A, 3A, J, J, J, J, J, J, J, D, D, D, D, D, D, 3A, 1A, D, 23A, D, D, D, 3A, 23A, J, 1A, 3A, 1A, D, 3A, 6A, J, D, D, D, D, D, D, D, D, D, D, 23A, D, D, D, D, 3A, D, D, D, 1A, J, J, J, J, J, J, 23A, J, 5A, 7A, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, 13A, D, 24A, 23A, D, 1A, 1A, D, J, D, D, 1A, 1A, D, 4A, J, D, J, 8A, D, J, D, J, J, 23A, 23A, D, 8A, D, 29A, D, 3A, D, J, D, D, 3A, 6A, 5A, 2A, D, 3A, J, 1A, 1A, 1A, 1A, D, D, D, D, D, D, 4A, D, 4A, D, 4A, J, 4A, 3A, 1A, J, J, J, 1A, D, D, D, D, J, 3A, J, J, 6A, J, D, D, 11A, 1A, 1A, 3A, 8A, 1A, D, J, D, J, J, 5A, 4A, J, D, D, D, D, D, D, D, D, D, D, D, D, D, D, 8A, 8A, 3A, 3A] 5 EXP 1 ImageHDU 8 (5185, 5166) float32
SIMPLE = T / Fits standard BITPIX = 16 / Bits per pixel NAXIS = 0 / Number of axes EXTEND = T / File may contain extensions ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator DATE = '2017-08-23 ' / Date FITS file was generated IRAF-TLM= '2017-08-23T16:02:22' / Time of last modification COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H TELESCOP= 'HST' / telescope used to acquire data INSTRUME= 'ACS ' / identifier for instrument used to acquire data / DATA DESCRIPTION KEYWORDS IMAGETYP= 'EXT ' / type of exposure identifier PRIMESI = 'ACS ' / instrument designated as prime / TARGET INFORMATION TARGNAME= 'NGC224-HALO ' / proposer's target name / PROPOSAL INFORMATION PROPOSID= 9453 / PEP proposal identifier / EXPOSURE INFORMATION SUNANGLE= 128.445786 / angle between sun and V1 axis MOONANGL= 137.613449 / angle between moon and V1 axis SUN_ALT = 24.114336 / altitude of the sun above Earth's limb FGSLOCK = 'FINE ' / commanded FGS lock (FINE,COARSE,GYROS,UNKNOWN) GYROMODE= '3' / number of gyros scheduled, T=3+OBAD DATE-OBS= '2002-12-03' / UT date of start of observation (yyyy-mm-dd) EXPSTART= 5.261119857568E+04 / exposure start time (Modified Julian Date) EXPEND = 52611.33419003 / exposure end time (Modified Julian Date) EXPTIME = 5060.0 / exposure duration (seconds)--calculated TEXPTIME= 5060.0 EXPFLAG = 'NORMAL ' / Exposure interruption indicator DARKTIME= 1233.247202 / fiducial pixel dark time (secs) / POINTING INFORMATION / TARGET OFFSETS (POSTARGS) POSTARG1= -0.227366 / POSTARG in axis 1 direction POSTARG2= 2.969150 / POSTARG in axis 2 direction / DIAGNOSTIC KEYWORDS CSYS_VER= 'hstdp-2016.2' / Calibration software system version id / SCIENCE INSTRUMENT CONFIGURATION OBSTYPE = 'IMAGING ' / observation type - imaging or spectroscopic OBSMODE = 'ACCUM ' / operating mode CTEIMAGE= 'NONE' / type of Charge Transfer Image, if applicable SCLAMP = 'NONE ' / lamp status, NONE or name of lamp which is on NRPTEXP = 1 / number of repeat exposures in set: default 1 SUBARRAY= F / data from a subarray (T) or full frame (F) DETECTOR= 'WFC' / detector in use: WFC, HRC, or SBC FILTER1 = 'CLEAR1L ' / element selected from filter wheel 1 FILTER2 = 'F814W ' / element selected from filter wheel 2 FW1OFFST= 0 / computed filter wheel offset FW1ERROR= F / filter wheel position error flag FW2OFFST= 0 / computed filter wheel offset FW2ERROR= F / filter wheel position error flag FWSOFFST= 0 / computed filter wheel offset FWSERROR= F / filter wheel position error flag LRFWAVE = 0.000000 / proposed linear ramp filter wavelength APERTURE= 'WFC ' / aperture name PROPAPER= 'WFC ' / proposed aperture name DIRIMAGE= 'NONE ' / direct image for grism or prism exposure CTEDIR = 'NONE ' / CTE measurement direction: serial or parallel CRSPLIT = 1 / number of cosmic ray split exposures / CALIBRATION SWITCHES: PERFORM, OMIT, COMPLETE WRTERR = T / write out error array extension DQICORR = 'COMPLETE' / data quality initialization ATODCORR= 'OMIT ' / correct for A to D conversion errors BLEVCORR= 'COMPLETE' / subtract bias level computed from overscan img BIASCORR= 'COMPLETE' / Subtract bias image FLSHCORR= 'OMIT ' / post flash correction CRCORR = 'OMIT ' / combine observations to reject cosmic rays EXPSCORR= 'COMPLETE' / process individual observations after cr-reject SHADCORR= 'OMIT ' / apply shutter shading correction DARKCORR= 'COMPLETE' / Subtract dark image FLATCORR= 'COMPLETE' / flat field data PHOTCORR= 'COMPLETE' / populate photometric header keywords DRIZCORR= 'COMPLETE' / drizzle processing / CALIBRATION REFERENCE FILES BPIXTAB = 'jref$n9p09145j_bpx.fits' / bad pixel table CCDTAB = 'jref$o151506dj_ccd.fits' / CCD calibration parameters ATODTAB = 'N/A ' / analog to digital correction file OSCNTAB = 'jref$lch1459bj_osc.fits' / CCD overscan table BIASFILE= 'jref$mc91206ej_bia.fits' / bias image file name FLSHFILE= 'N/A ' / post flash correction file name CRREJTAB= 'N/A ' / cosmic ray rejection parameters SHADFILE= 'N/A ' / shutter shading correction file DARKFILE= 'jref$mc91206cj_drk.fits' / dark image file name DFLTFILE= 'N/A ' / delta flat field file name LFLTFILE= 'N/A ' / low order flat PHOTTAB = 'N/A ' / Photometric throughput table GRAPHTAB= 'N/A ' / the HST graph table COMPTAB = 'N/A ' / the HST components table IDCTAB = 'jref$0461802ej_idc.fits' / image distortion correction table DGEOFILE= 'jref$qbu16429j_dxy.fits' / Distortion correction image MDRIZTAB= 'jref$1371334kj_mdz.fits' / MultiDrizzle parameter table CFLTFILE= 'N/A ' / Coronagraphic spot image SPOTTAB = 'N/A ' / Coronagraphic spot offset table IMPHTTAB= 'jref$08b18470j_imp.fits' / Image Photometry Table / COSMIC RAY REJECTION ALGORITHM PARAMETERS MEANEXP = 0.000000 / reference exposure time for parameters SCALENSE= 0.000000 / multiplicative scale factor applied to noise SKYSUB = ' ' / sky value subtracted (MODE or NONE) SKYSUM = 0.0 / sky level from the sum of all constituent image CRRADIUS= 0.000000 / rejection propagation radius (pixels) CRTHRESH= 0.000000 / rejection propagation threshold BADINPDQ= 0 / data quality flag bits to reject REJ_RATE= 0.0 / rate at which pixels are affected by cosmic ray CRMASK = F / flag CR-rejected pixels in input files (T/F) / OTFR KEYWORDS / PATTERN KEYWORDS PATTERN1= 'NONE ' / primary pattern type P1_SHAPE= ' ' / primary pattern shape P1_PURPS= ' ' / primary pattern purpose P1_NPTS = 0 / number of points in primary pattern P1_PSPAC= 0.000000 / point spacing for primary pattern (arc-sec) P1_LSPAC= 0.000000 / line spacing for primary pattern (arc-sec) P1_ANGLE= 0.000000 / angle between sides of parallelogram patt (deg) P1_FRAME= ' ' / coordinate frame of primary pattern P1_ORINT= 0.000000 / orientation of pattern to coordinate frame (deg P1_CENTR= ' ' / center pattern relative to pointing (yes/no) PATTSTEP= 0 / position number of this point in the pattern / POST FLASH PARAMETERS FLASHDUR= 0.0 / Exposure time in seconds: 0.1 to 409.5 FLASHCUR= 'OFF ' / Post flash current: OFF, LOW, MED, HIGH FLASHSTA= 'NOT PERFORMED ' / Status: SUCCESSFUL, ABORTED, NOT PERFORMED SHUTRPOS= 'MULTIPLE' / Shutter position: A or B / ENGINEERING PARAMETERS CCDGAIN = 1 / commanded gain of CCD CCDOFSTA= 3 / commanded CCD bias offset for amplifier A CCDOFSTB= 3 / commanded CCD bias offset for amplifier B CCDOFSTC= 3 / commanded CCD bias offset for amplifier C CCDOFSTD= 3 / commanded CCD bias offset for amplifier D / CALIBRATED ENGINEERING PARAMETERS ATODGNA = 9.9989998E-01 / calibrated gain for amplifier A ATODGNB = 9.7210002E-01 / calibrated gain for amplifier B ATODGNC = 1.0107000E+00 / calibrated gain for amplifier C ATODGND = 1.0180000E+00 / calibrated gain for amplifier D READNSEA= 4.9699998E+00 / calibrated read noise for amplifier A READNSEB= 4.8499999E+00 / calibrated read noise for amplifier B READNSEC= 5.2399998E+00 / calibrated read noise for amplifier C READNSED= 4.8499999E+00 / calibrated read noise for amplifier D / ASSOCIATION KEYWORDS CRDS_CTX= 'hst_0515.pmap' CRDS_VER= '7.0.1, opus_2016.1-universal, af27872' UPWCSVER= '1.2.3.dev' / Version of STWCS used to updated the WCS PYWCSVER= '1.2.1 ' / Version of PYWCS used to updated the WCS HISTORY CCD parameters table: HISTORY reference table jref$o151506dj_ccd.fits HISTORY inflight HISTORY June 2002 HISTORY DQICORR complete ... HISTORY values checked for saturation HISTORY DQ array initialized ... HISTORY reference table jref$n9p09145j_bpx.fits HISTORY BIASCORR complete ... HISTORY reference image jref$mc91206ej_bia.fits HISTORY INFLIGHT 09/11/2002 04/12/2002 HISTORY Superbias created by Doug Van Orsow from proposal 9647 HISTORY BLEVCORR complete; bias level from overscan was subtracted. HISTORY BLEVCORR does not include correction for drift along lines. HISTORY Overscan region table: HISTORY reference table jref$lch1459bj_osc.fits HISTORY Uncertainty array initialized. HISTORY CCD parameters table: HISTORY reference table jref$o151506dj_ccd.fits HISTORY inflight HISTORY June 2002 PCTETRSH= -1.000000000000E+01 / PCTE over subtraction threshold HISTORY PCTECORR complete ... HISTORY reference table jref$xa81724cj_cte.fits DISTNAME= 'j8f802lcq_0461802ej-02c1450rj-02c1450oj' SIPNAME = 'j8f802lcq_0461802ej' RULESVER= 1.1 / Version ID for header kw rules file BLENDVER= '1.2.0 ' / Version of blendheader software used RULEFILE= '/home/pipemgr/anaconda2/lib/python2.7/site-packages/fitsblender/acs&' CONTINUE '_header.rules&' CONTINUE '' / Name of header kw rules file NEXTEND = 4 FILENAME= 'hst_9453_02_acs_wfc_f814w_run_1_drz.fits' PROD_VER= 'DrizzlePac 2.1.6' ROOTNAME= 'hst_9453_02_acs_wfc_f814w_run_1' ASN_MTYP= 'PROD-DTH' NDRIZIM = 8 / Drizzle, No. images drizzled onto output D001OUDA= 'hst_9453_02_acs_wfc_f814w_run_1_drz.fits' / Drizzle, output data imag D001VER = 'Callable C-based DRIZZLE Version 0.8 (20th M' / Drizzle, task version D001SCAL= 0.05 / Drizzle, pixel size (arcsec) of output image D001COEF= 'SIP ' / Drizzle, source of coefficients D001OUWE= 'hst_9453_02_acs_wfc_f814w_run_1_drz_wht.fits' / Drizzle, output weigh D001OUCO= 'hst_9453_02_acs_wfc_f814w_run_1_drz_ctx.fits' / Drizzle, output conte D001WTSC= 1 / Drizzle, weighting factor for input image D001MASK= 'j8f802lcq_sci1_final_mask.fits' / Drizzle, input weighting image D001FVAL= 'INDEF ' / Drizzle, fill value for zero weight output pix D001WKEY= '' / Input image WCS Version used D001OUUN= 'cps ' / Drizzle, units of output image - counts or cps D001KERN= 'square ' / Drizzle, form of weight distribution kernel D001GEOM= 'wcs ' / Drizzle, source of geometric information D001ISCL= 0.05 / Drizzle, default IDCTAB pixel size(arcsec) D001PIXF= 1.0 / Drizzle, linear size of drop D001DATA= 'j8f802lcq_flc.fits[sci,1]' / Drizzle, input data image D001DEXP= 1230.0 / Drizzle, input image exposure time (s) D002OUDA= 'hst_9453_02_acs_wfc_f814w_run_1_drz.fits' / Drizzle, output data imag D002VER = 'Callable C-based DRIZZLE Version 0.8 (20th M' / Drizzle, task version D002SCAL= 0.05 / Drizzle, pixel size (arcsec) of output image D002COEF= 'SIP ' / Drizzle, source of coefficients D002OUWE= 'hst_9453_02_acs_wfc_f814w_run_1_drz_wht.fits' / Drizzle, output weigh D002OUCO= 'hst_9453_02_acs_wfc_f814w_run_1_drz_ctx.fits' / Drizzle, output conte D002WTSC= 1 / Drizzle, weighting factor for input image D002MASK= 'j8f802lcq_sci2_final_mask.fits' / Drizzle, input weighting image D002FVAL= 'INDEF ' / Drizzle, fill value for zero weight output pix D002WKEY= '' / Input image WCS Version used D002OUUN= 'cps ' / Drizzle, units of output image - counts or cps D002KERN= 'square ' / Drizzle, form of weight distribution kernel D002GEOM= 'wcs ' / Drizzle, source of geometric information D002ISCL= 0.05 / Drizzle, default IDCTAB pixel size(arcsec) D002PIXF= 1.0 / Drizzle, linear size of drop D002DATA= 'j8f802lcq_flc.fits[sci,2]' / Drizzle, input data image D002DEXP= 1230.0 / Drizzle, input image exposure time (s) D003OUDA= 'hst_9453_02_acs_wfc_f814w_run_1_drz.fits' / Drizzle, output data imag D003VER = 'Callable C-based DRIZZLE Version 0.8 (20th M' / Drizzle, task version D003SCAL= 0.05 / Drizzle, pixel size (arcsec) of output image D003COEF= 'SIP ' / Drizzle, source of coefficients D003OUWE= 'hst_9453_02_acs_wfc_f814w_run_1_drz_wht.fits' / Drizzle, output weigh D003OUCO= 'hst_9453_02_acs_wfc_f814w_run_1_drz_ctx.fits' / Drizzle, output conte D003WTSC= 1 / Drizzle, weighting factor for input image D003MASK= 'j8f802loq_sci1_final_mask.fits' / Drizzle, input weighting image D003FVAL= 'INDEF ' / Drizzle, fill value for zero weight output pix D003WKEY= '' / Input image WCS Version used D003OUUN= 'cps ' / Drizzle, units of output image - counts or cps D003KERN= 'square ' / Drizzle, form of weight distribution kernel D003GEOM= 'wcs ' / Drizzle, source of geometric information D003ISCL= 0.05 / Drizzle, default IDCTAB pixel size(arcsec) D003PIXF= 1.0 / Drizzle, linear size of drop D003DATA= 'j8f802loq_flc.fits[sci,1]' / Drizzle, input data image D003DEXP= 1230.0 / Drizzle, input image exposure time (s) D004OUDA= 'hst_9453_02_acs_wfc_f814w_run_1_drz.fits' / Drizzle, output data imag D004VER = 'Callable C-based DRIZZLE Version 0.8 (20th M' / Drizzle, task version D004SCAL= 0.05 / Drizzle, pixel size (arcsec) of output image D004COEF= 'SIP ' / Drizzle, source of coefficients D004OUWE= 'hst_9453_02_acs_wfc_f814w_run_1_drz_wht.fits' / Drizzle, output weigh D004OUCO= 'hst_9453_02_acs_wfc_f814w_run_1_drz_ctx.fits' / Drizzle, output conte D004WTSC= 1 / Drizzle, weighting factor for input image D004MASK= 'j8f802loq_sci2_final_mask.fits' / Drizzle, input weighting image D004FVAL= 'INDEF ' / Drizzle, fill value for zero weight output pix D004WKEY= '' / Input image WCS Version used D004OUUN= 'cps ' / Drizzle, units of output image - counts or cps D004KERN= 'square ' / Drizzle, form of weight distribution kernel D004GEOM= 'wcs ' / Drizzle, source of geometric information D004ISCL= 0.05 / Drizzle, default IDCTAB pixel size(arcsec) D004PIXF= 1.0 / Drizzle, linear size of drop D004DATA= 'j8f802loq_flc.fits[sci,2]' / Drizzle, input data image D004DEXP= 1230.0 / Drizzle, input image exposure time (s) D005OUDA= 'hst_9453_02_acs_wfc_f814w_run_1_drz.fits' / Drizzle, output data imag D005VER = 'Callable C-based DRIZZLE Version 0.8 (20th M' / Drizzle, task version D005SCAL= 0.05 / Drizzle, pixel size (arcsec) of output image D005COEF= 'SIP ' / Drizzle, source of coefficients D005OUWE= 'hst_9453_02_acs_wfc_f814w_run_1_drz_wht.fits' / Drizzle, output weigh D005OUCO= 'hst_9453_02_acs_wfc_f814w_run_1_drz_ctx.fits' / Drizzle, output conte D005WTSC= 1 / Drizzle, weighting factor for input image D005MASK= 'j8f802m1q_sci1_final_mask.fits' / Drizzle, input weighting image D005FVAL= 'INDEF ' / Drizzle, fill value for zero weight output pix D005WKEY= '' / Input image WCS Version used D005OUUN= 'cps ' / Drizzle, units of output image - counts or cps D005KERN= 'square ' / Drizzle, form of weight distribution kernel D005GEOM= 'wcs ' / Drizzle, source of geometric information D005ISCL= 0.05 / Drizzle, default IDCTAB pixel size(arcsec) D005PIXF= 1.0 / Drizzle, linear size of drop D005DATA= 'j8f802m1q_flc.fits[sci,1]' / Drizzle, input data image D005DEXP= 1300.0 / Drizzle, input image exposure time (s) D006OUDA= 'hst_9453_02_acs_wfc_f814w_run_1_drz.fits' / Drizzle, output data imag D006VER = 'Callable C-based DRIZZLE Version 0.8 (20th M' / Drizzle, task version D006SCAL= 0.05 / Drizzle, pixel size (arcsec) of output image D006COEF= 'SIP ' / Drizzle, source of coefficients D006OUWE= 'hst_9453_02_acs_wfc_f814w_run_1_drz_wht.fits' / Drizzle, output weigh D006OUCO= 'hst_9453_02_acs_wfc_f814w_run_1_drz_ctx.fits' / Drizzle, output conte D006WTSC= 1 / Drizzle, weighting factor for input image D006MASK= 'j8f802m1q_sci2_final_mask.fits' / Drizzle, input weighting image D006FVAL= 'INDEF ' / Drizzle, fill value for zero weight output pix D006WKEY= '' / Input image WCS Version used D006OUUN= 'cps ' / Drizzle, units of output image - counts or cps D006KERN= 'square ' / Drizzle, form of weight distribution kernel D006GEOM= 'wcs ' / Drizzle, source of geometric information D006ISCL= 0.05 / Drizzle, default IDCTAB pixel size(arcsec) D006PIXF= 1.0 / Drizzle, linear size of drop D006DATA= 'j8f802m1q_flc.fits[sci,2]' / Drizzle, input data image D006DEXP= 1300.0 / Drizzle, input image exposure time (s) D007OUDA= 'hst_9453_02_acs_wfc_f814w_run_1_drz.fits' / Drizzle, output data imag D007VER = 'Callable C-based DRIZZLE Version 0.8 (20th M' / Drizzle, task version D007SCAL= 0.05 / Drizzle, pixel size (arcsec) of output image D007COEF= 'SIP ' / Drizzle, source of coefficients D007OUWE= 'hst_9453_02_acs_wfc_f814w_run_1_drz_wht.fits' / Drizzle, output weigh D007OUCO= 'hst_9453_02_acs_wfc_f814w_run_1_drz_ctx.fits' / Drizzle, output conte D007WTSC= 1 / Drizzle, weighting factor for input image D007MASK= 'j8f802mfq_sci1_final_mask.fits' / Drizzle, input weighting image D007FVAL= 'INDEF ' / Drizzle, fill value for zero weight output pix D007WKEY= '' / Input image WCS Version used D007OUUN= 'cps ' / Drizzle, units of output image - counts or cps D007KERN= 'square ' / Drizzle, form of weight distribution kernel D007GEOM= 'wcs ' / Drizzle, source of geometric information D007ISCL= 0.05 / Drizzle, default IDCTAB pixel size(arcsec) D007PIXF= 1.0 / Drizzle, linear size of drop D007DATA= 'j8f802mfq_flc.fits[sci,1]' / Drizzle, input data image D007DEXP= 1300.0 / Drizzle, input image exposure time (s) D008OUDA= 'hst_9453_02_acs_wfc_f814w_run_1_drz.fits' / Drizzle, output data imag D008VER = 'Callable C-based DRIZZLE Version 0.8 (20th M' / Drizzle, task version D008SCAL= 0.05 / Drizzle, pixel size (arcsec) of output image D008COEF= 'SIP ' / Drizzle, source of coefficients D008OUWE= 'hst_9453_02_acs_wfc_f814w_run_1_drz_wht.fits' / Drizzle, output weigh D008OUCO= 'hst_9453_02_acs_wfc_f814w_run_1_drz_ctx.fits' / Drizzle, output conte D008WTSC= 1 / Drizzle, weighting factor for input image D008MASK= 'j8f802mfq_sci2_final_mask.fits' / Drizzle, input weighting image D008FVAL= 'INDEF ' / Drizzle, fill value for zero weight output pix D008WKEY= '' / Input image WCS Version used D008OUUN= 'cps ' / Drizzle, units of output image - counts or cps D008KERN= 'square ' / Drizzle, form of weight distribution kernel D008GEOM= 'wcs ' / Drizzle, source of geometric information D008ISCL= 0.05 / Drizzle, default IDCTAB pixel size(arcsec) D008PIXF= 1.0 / Drizzle, linear size of drop D008DATA= 'j8f802mfq_flc.fits[sci,2]' / Drizzle, input data image D008DEXP= 1300.0 / Drizzle, input image exposure time (s) PA_V3 = 255.070404 / position angle of V3-axis of HST (deg) HISTORY CCD parameters table: HISTORY reference table jref$o151506dj_ccd.fits HISTORY inflight HISTORY June 2002 HISTORY DARKCORR complete ... HISTORY reference image jref$z552106tj_dkc.fits HISTORY INFLIGHT 24/09/2002 07/10/2002 HISTORY CTE corrected dark for data taken after Oct 07 2002---------------- HISTORY FLATCORR complete ... HISTORY reference image jref$n6u12592j_pfl.fits HISTORY INFLIGHT 18/04/2002 - 09/05/2002 HISTORY Ground LP-flat corrected by L-flat from program 9018 HISTORY PHOTCORR complete ... HISTORY reference table jref$08b18470j_imp.fits HISTORY EXPSCORR complete ... HISTORY ============================================================ HISTORY Header Generation rules: HISTORY Rules used to combine headers of input files HISTORY Start of rules... HISTORY ------------------------------------------------------------ HISTORY !VERSION = 1.1 HISTORY !INSTRUMENT = ACS HISTORY ROOTNAME HISTORY EXTNAME HISTORY EXTVER HISTORY A_0_2 HISTORY A_0_3 HISTORY A_0_4 HISTORY A_1_1 HISTORY A_1_2 HISTORY A_1_3 HISTORY A_2_0 HISTORY A_2_1 HISTORY A_2_2 HISTORY A_3_0 HISTORY A_3_1 HISTORY A_4_0 HISTORY ACQNAME HISTORY A_ORDER HISTORY APERTURE HISTORY ASN_ID HISTORY ASN_MTYP HISTORY ASN_TAB HISTORY ATODCORR HISTORY ATODGNA HISTORY ATODGNB HISTORY ATODGNC HISTORY ATODGND HISTORY ATODTAB HISTORY B_0_2 HISTORY B_0_3 HISTORY B_0_4 HISTORY B_1_1 HISTORY B_1_2 HISTORY B_1_3 HISTORY B_2_0 HISTORY B_2_1 HISTORY B_2_2 HISTORY B_3_0 HISTORY B_3_1 HISTORY B_4_0 HISTORY BADINPDQ HISTORY BIASCORR HISTORY BIASFILE HISTORY BIASLEVA HISTORY BIASLEVB HISTORY BIASLEVC HISTORY BIASLEVD HISTORY BINAXIS1 HISTORY BINAXIS2 HISTORY BITPIX HISTORY BLEVCORR HISTORY B_ORDER HISTORY BPIXTAB HISTORY BUNIT HISTORY CAL_VER HISTORY CBLKSIZ HISTORY CCDAMP HISTORY CCDCHIP HISTORY CCDGAIN HISTORY CCDOFSTA HISTORY CCDOFSTB HISTORY CCDOFSTC HISTORY CCDOFSTD HISTORY CCDTAB HISTORY CD1_1 HISTORY CD1_2 HISTORY CD2_1 HISTORY CD2_2 HISTORY CENTERA1 HISTORY CENTERA2 HISTORY CFLTFILE HISTORY COMPTAB HISTORY COMPTYP HISTORY CRCORR HISTORY CRMASK HISTORY CRPIX1 HISTORY CRPIX2 HISTORY CRRADIUS HISTORY CRREJTAB HISTORY CRSIGMAS HISTORY CRSPLIT HISTORY CRTHRESH HISTORY CRVAL1 HISTORY CRVAL2 HISTORY CTE_NAME HISTORY CTE_VER HISTORY CTEDIR HISTORY CTEIMAGE HISTORY CTYPE1 HISTORY CTYPE2 HISTORY D2IMFILE HISTORY DARKCORR HISTORY DARKFILE HISTORY DATE HISTORY DATE-OBS HISTORY DEC_APER HISTORY DEC_TARG HISTORY DETECTOR HISTORY DFLTFILE HISTORY DGEOFILE HISTORY DIRIMAGE HISTORY DQICORR HISTORY DRIZCORR HISTORY DRKCFILE HISTORY EQUINOX HISTORY ERRCNT HISTORY EXPEND HISTORY EXPFLAG HISTORY EXPNAME HISTORY EXPSCORR HISTORY EXPSTART HISTORY EXPTIME HISTORY EXTEND HISTORY FGSLOCK HISTORY FILENAME HISTORY FILETYPE HISTORY FILLCNT HISTORY FILTER1 HISTORY FILTER2 HISTORY FLASHCUR HISTORY FLASHDUR HISTORY FLASHSTA HISTORY FLATCORR HISTORY FLSHCORR HISTORY FLSHFILE HISTORY FW1ERROR HISTORY FW1OFFST HISTORY FW2ERROR HISTORY FW2OFFST HISTORY FWSERROR HISTORY FWSOFFST HISTORY GCOUNT HISTORY GLINCORR HISTORY GLOBLIM HISTORY GLOBRATE HISTORY GOODMAX HISTORY GOODMEAN HISTORY GOODMIN HISTORY GRAPHTAB HISTORY GYROMODE HISTORY IDCSCALE HISTORY IDCTAB HISTORY IDCTHETA HISTORY IDCV2REF HISTORY IDCV3REF HISTORY IMAGETYP HISTORY IMPHTTAB HISTORY INHERIT HISTORY INITGUES HISTORY INSTRUME HISTORY JWROTYPE HISTORY LFLGCORR HISTORY LFLTFILE HISTORY LINENUM HISTORY LOSTPIX HISTORY LRC_FAIL HISTORY LRC_XSTS HISTORY LRFWAVE HISTORY LTM1_1 HISTORY LTM2_2 HISTORY LTV1 HISTORY LTV2 HISTORY MDECODT1 HISTORY MDECODT2 HISTORY MDRIZSKY HISTORY MDRIZTAB HISTORY MEANBLEV HISTORY MEANDARK HISTORY MEANEXP HISTORY MEANFLSH HISTORY MLINTAB HISTORY MOFFSET1 HISTORY MOFFSET2 HISTORY MOONANGL HISTORY MTFLAG HISTORY NAXIS HISTORY NAXIS1 HISTORY NAXIS2 HISTORY NCOMBINE HISTORY NEXTEND HISTORY NGOODPIX HISTORY NPOLFILE HISTORY NRPTEXP HISTORY OBSMODE HISTORY OBSTYPE HISTORY OCD1_1 HISTORY OCD1_2 HISTORY OCD2_1 HISTORY OCD2_2 HISTORY OCRPIX1 HISTORY OCRPIX2 HISTORY OCRVAL1 HISTORY OCRVAL2 HISTORY OCTYPE1 HISTORY OCTYPE2 HISTORY OCX10 HISTORY OCX11 HISTORY OCY10 HISTORY OCY11 HISTORY ONAXIS1 HISTORY ONAXIS2 HISTORY OORIENTA HISTORY OPUS_VER HISTORY ORIENTAT HISTORY ORIGIN HISTORY OSCNTAB HISTORY P1_ANGLE HISTORY P1_CENTR HISTORY P1_FRAME HISTORY P1_LSPAC HISTORY P1_NPTS HISTORY P1_ORINT HISTORY P1_PSPAC HISTORY P1_PURPS HISTORY P1_SHAPE HISTORY PA_APER HISTORY PATTERN1 HISTORY PATTSTEP HISTORY PA_V3 HISTORY PCOUNT HISTORY PCTECORR HISTORY PCTEFRAC HISTORY PCTENSMD HISTORY PCTERNCL HISTORY PCTESHFT HISTORY PCTESMIT HISTORY PCTETAB HISTORY PFLTFILE HISTORY PHOTBW HISTORY PHOTCORR HISTORY PHOTFLAM HISTORY PHOTMODE HISTORY PHOTPLAM HISTORY PHOTTAB HISTORY PHOTZPT HISTORY PODPSFF HISTORY POSTARG1 HISTORY POSTARG2 HISTORY PRIMESI HISTORY PR_INV_F HISTORY PR_INV_L HISTORY PR_INV_M HISTORY PROCTIME HISTORY PROPAPER HISTORY PROPOSID HISTORY QUALCOM1 HISTORY QUALCOM2 HISTORY QUALCOM3 HISTORY QUALITY HISTORY RA_APER HISTORY RA_TARG HISTORY READNSEA HISTORY READNSEB HISTORY READNSEC HISTORY READNSED HISTORY REFFRAME HISTORY REJ_RATE HISTORY RPTCORR HISTORY SCALENSE HISTORY SCLAMP HISTORY SDQFLAGS HISTORY SHADCORR HISTORY SHADFILE HISTORY SHUTRPOS HISTORY SIMPLE HISTORY SIZAXIS1 HISTORY SIZAXIS2 HISTORY SKYSUB HISTORY SKYSUM HISTORY SNRMAX HISTORY SNRMEAN HISTORY SNRMIN HISTORY SOFTERRS HISTORY SPOTTAB HISTORY STATFLAG HISTORY STDCFFF HISTORY STDCFFP HISTORY SUBARRAY HISTORY SUN_ALT HISTORY SUNANGLE HISTORY TARGNAME HISTORY TDDALPHA HISTORY TDDBETA HISTORY TELESCOP HISTORY TIME-OBS HISTORY T_SGSTAR HISTORY VAFACTOR HISTORY WCSAXES HISTORY WCSCDATE HISTORY WFCMPRSD HISTORY WRTERR HISTORY XTENSION HISTORY WCSNAMEO HISTORY WCSAXESO HISTORY LONPOLEO HISTORY LATPOLEO HISTORY RESTFRQO HISTORY RESTWAVO HISTORY CD1_1O HISTORY CD1_2O HISTORY CD2_1O HISTORY CD2_2O HISTORY CDELT1O HISTORY CDELT2O HISTORY CRPIX1O HISTORY CRPIX2O HISTORY CRVAL1O HISTORY CRVAL2O HISTORY CTYPE1O HISTORY CTYPE2O HISTORY CUNIT1O HISTORY CUNIT2O HISTORY APERTURE APERTURE multi HISTORY DETECTOR DETECTOR first HISTORY EXPEND EXPEND max HISTORY EXPSTART EXPSTART min HISTORY EXPTIME TEXPTIME sum HISTORY EXPTIME EXPTIME sum HISTORY FILTER1 FILTER1 multi HISTORY FILTER2 FILTER2 multi HISTORY GOODMAX GOODMAX max HISTORY GOODMEAN GOODMEAN mean HISTORY GOODMIN GOODMIN min HISTORY INHERIT INHERIT first # maintain IRAF compatibility HISTORY INSTRUME INSTRUME first HISTORY LRFWAVE LRFWAVE first HISTORY NCOMBINE NCOMBINE sum HISTORY MDRIZSKY MDRIZSKY mean HISTORY PHOTBW PHOTBW mean HISTORY PHOTFLAM PHOTFLAM mean HISTORY PHOTMODE PHOTMODE first HISTORY PHOTPLAM PHOTPLAM mean HISTORY PHOTZPT PHOTZPT mean HISTORY PROPOSID PROPOSID first HISTORY SNRMAX SNRMAX max HISTORY SNRMEAN SNRMEAN mean HISTORY SNRMIN SNRMIN min HISTORY TARGNAME TARGNAME first HISTORY TELESCOP TELESCOP first HISTORY ATODCORR ATODCORR multi HISTORY ATODGNA ATODGNA first HISTORY ATODGNB ATODGNB first HISTORY ATODGNC ATODGNC first HISTORY ATODGND ATODGND first HISTORY ATODTAB ATODTAB multi HISTORY BADINPDQ BADINPDQ sum HISTORY BIASCORR BIASCORR multi HISTORY BIASFILE BIASFILE multi HISTORY BLEVCORR BLEVCORR multi HISTORY BPIXTAB BPIXTAB multi HISTORY CCDCHIP CCDCHIP first HISTORY CCDGAIN CCDGAIN first HISTORY CCDOFSTA CCDOFSTA first HISTORY CCDOFSTB CCDOFSTB first HISTORY CCDOFSTC CCDOFSTC first HISTORY CCDOFSTD CCDOFSTD first HISTORY CCDTAB CCDTAB multi HISTORY CFLTFILE CFLTFILE multi HISTORY COMPTAB COMPTAB multi HISTORY CRCORR CRCORR multi HISTORY CRMASK CRMASK first HISTORY CRRADIUS CRRADIUS first HISTORY CRREJTAB CRREJTAB multi HISTORY CRSPLIT CRSPLIT first HISTORY CRTHRESH CRTHRESH first HISTORY CTEDIR CTEDIR multi HISTORY CTEIMAGE CTEIMAGE first HISTORY DARKCORR DARKCORR multi HISTORY DARKFILE DARKFILE multi HISTORY DATE-OBS DATE-OBS first HISTORY DEC_APER DEC_APER first HISTORY DFLTFILE DFLTFILE multi HISTORY DGEOFILE DGEOFILE multi HISTORY DIRIMAGE DIRIMAGE multi HISTORY DQICORR DQICORR multi HISTORY DRIZCORR DRIZCORR multi HISTORY EXPFLAG EXPFLAG multi HISTORY EXPSCORR EXPSCORR multi HISTORY FGSLOCK FGSLOCK multi HISTORY FLASHCUR FLASHCUR multi HISTORY FLASHDUR FLASHDUR first HISTORY FLASHSTA FLASHSTA first HISTORY FLATCORR FLATCORR multi HISTORY FLSHCORR FLSHCORR multi HISTORY FLSHFILE FLSHFILE multi HISTORY FW1ERROR FW1ERROR multi HISTORY FW1OFFST FW1OFFST first HISTORY FW2ERROR FW2ERROR multi HISTORY FW2OFFST FW2OFFST first HISTORY FWSERROR FWSERROR multi HISTORY FWSOFFST FWSOFFST first HISTORY GRAPHTAB GRAPHTAB multi HISTORY GYROMODE GYROMODE multi HISTORY IDCTAB IDCTAB multi HISTORY IMAGETYP IMAGETYP first HISTORY IMPHTTAB IMPHTTAB multi HISTORY LFLGCORR LFLGCORR multi HISTORY LFLTFILE LFLTFILE multi HISTORY LTM1_1 LTM1_1 float_one HISTORY LTM2_2 LTM2_2 float_one HISTORY MDRIZTAB MDRIZTAB multi HISTORY MEANEXP MEANEXP first HISTORY MOONANGL MOONANGL first HISTORY NRPTEXP NRPTEXP first HISTORY OBSMODE OBSMODE multi HISTORY OBSTYPE OBSTYPE first HISTORY OSCNTAB OSCNTAB multi HISTORY P1_ANGLE P1_ANGLE first HISTORY P1_CENTR P1_CENTR multi HISTORY P1_FRAME P1_FRAME multi HISTORY P1_LSPAC P1_LSPAC first HISTORY P1_NPTS P1_NPTS first HISTORY P1_ORINT P1_ORINT first HISTORY P1_PSPAC P1_PSPAC first HISTORY P1_PURPS P1_PURPS multi HISTORY P1_SHAPE P1_SHAPE multi HISTORY P2_ANGLE P2_ANGLE first HISTORY P2_CENTR P2_CENTR multi HISTORY P2_FRAME P2_FRAME multi HISTORY P2_LSPAC P2_LSPAC first HISTORY P2_NPTS P2_NPTS first HISTORY P2_ORINT P2_ORINT first HISTORY P2_PSPAC P2_PSPAC first HISTORY P2_PURPS P2_PURPS multi HISTORY P2_SHAPE P2_SHAPE multi HISTORY PATTERN1 PATTERN1 multi HISTORY PATTERN2 PATTERN2 multi HISTORY PATTSTEP PATTSTEP first HISTORY PHOTCORR PHOTCORR multi HISTORY PHOTTAB PHOTTAB multi HISTORY POSTARG1 POSTARG1 first HISTORY POSTARG2 POSTARG2 first HISTORY PRIMESI PRIMESI multi HISTORY PROPAPER PROPAPER multi HISTORY RA_APER RA_APER first HISTORY READNSEA READNSEA first HISTORY READNSEB READNSEB first HISTORY READNSEC READNSEC first HISTORY READNSED READNSED first HISTORY REJ_RATE REJ_RATE first HISTORY SCALENSE SCALENSE first HISTORY SCLAMP SCLAMP multi HISTORY SHADCORR SHADCORR multi HISTORY SHADFILE SHADFILE multi HISTORY SHUTRPOS SHUTRPOS multi HISTORY SKYSUB SKYSUB multi HISTORY SKYSUM SKYSUM sum HISTORY SPOTTAB SPOTTAB multi HISTORY SUBARRAY SUBARRAY first HISTORY SUNANGLE SUNANGLE first HISTORY SUN_ALT SUN_ALT first HISTORY WRTERR WRTERR multi HISTORY ------------------------------------------------------------ HISTORY End of rules... HISTORY ============================================================ HISTORY AstroDrizzle processing performed using: HISTORY AstroDrizzle Version 2.1.6 HISTORY Numpy Version 1.11.3 HISTORY PyFITS Version 1.3 HISTORY All refereed publications based on data obtained from the HLA must HISTORY carry the following footnote: HISTORY HISTORY "Based on observations made with the NASA/ESA Hubble Space Telescope, HISTORY and obtained from the Hubble Legacy Archive, which is a collaboration HISTORY between the Space Telescope Science Institute (STScI/NASA), the Space HISTORY Telescope European Coordinating Facility (ST-ECF/ESA) and the HISTORY Canadian Astronomy Data Centre (CADC/NRC/CSA)." HISTORY HISTORY One copy of each paper resulting from data obtained from the HLA HISTORY should be sent to the STScI. HISTORY HISTORY In addition, publications of research supported by an STScI grant HISTORY must carry the following acknowledgment: HISTORY HISTORY "Support for Program number HISTORY ____________ was provided by NASA through a grant from the Space HISTORY Telescope Science Institute, which is operated by the Association of HISTORY Universities for Reasearch in Astronomy, Incorporated, under NASA HISTORY contract NAS5-26555."
# Inspect the array
image_data = fits.getdata(image_file)
print(type(image_data))
print(image_data.shape)
<class 'numpy.ndarray'> (5166, 5185)
from astropy.wcs import WCS
wcs = WCS(hdu_list[1].header)
pix = wcs.wcs.cdelt[1]*3600.0 # pixel size in arcsec
from astropy.visualization import (MinMaxInterval, SqrtStretch, ImageNormalize, ManualInterval)
# Create an ImageNormalize object
norm = ImageNormalize(image_data, interval = ManualInterval(-0.1,0.5))
print(wcs)
# Display the image
fig = plt.figure(figsize=(10,10),dpi=100)
ax = fig.add_subplot(111,projection=wcs)
im = ax.imshow(image_data, cmap='gray', origin='lower', norm=norm)
p1 = ax.scatter(hcv['ra'],hcv['dec'],transform=ax.get_transform('world'), \
s=30, edgecolor='salmon', facecolor='none', label='HCV')
p2 = ax.scatter(hcv_f814w['ra'],hcv_f814w['dec'],transform=ax.get_transform('world'), \
s=30, edgecolor='green', facecolor='none', label='RR Lyrae')
fig.colorbar(im)
ax.set_xlabel("RA")
ax.set_ylabel("Dec")
ax.legend(["HCV Sources", "RR Lyrae star"])
WCS Keywords Number of WCS axes: 2 CTYPE : 'RA---TAN' 'DEC--TAN' CRVAL : 11.52852151701503 40.70915097493346 CRPIX : 2592.999999999834 2583.499999999835 CD1_1 CD1_2 : -1.3888888888888e-05 0.0 CD2_1 CD2_2 : 0.0 1.38888888888896e-05 NAXIS : 5185 5166
<matplotlib.legend.Legend at 0x8285bae10>
In this example we are calling the HLA cutout service from the STScI to plot the cutouts (within a notebook markdown cell). The RR Lyrae star is in the centre of the image.
The brightest corrected magnitude for F814W is : lightcurve_cm = 24.820026461811118; and corresponds to image lightcurve_i = hst_9453_61_acs_wfc_f814w :
And the faintest corrected magnitude for F814W is : lightcurve_cm = 25.49112826041771; and corresponds to image lightcurve_i = hst_9453_59_acs_wfc_f814w
The brightest corrected magnitude for F606W is : lightcurve_cm = 24.827273229936768; and corresponds to image lightcurve_i = hst_9453_27_acs_wfc_f606w :
The faintest corrected magnitude for F606W is : lightcurve_cm = 25.61165283094627; and corresponds to image lightcurve_i = hst_9453_45_acs_wfc_f606w :