MJD date.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from astropy.io import fits"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "spectro= '../quicklook/spectro/data/20190101/'\n",
    "redux = '../quicklook/spectro/redux/'\n",
    "f1 = fits.open('%s%s'%(spectro,'desi-00000003.fits.fz'))\n",
    "f2 = fits.open('%s%s'%(spectro,'fibermap-00000003.fits'))\n",
    "f3 = fits.open('%s%s'%(redux,'/calib2d/20190101/fiberflat-b0-00000001.fits'))\n",
    "f4 = fits.open('%s%s'%(redux,'/exposures/20190101/00000003/sframe-b0-00000003.fits'))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "SIMPLE  =                    T / conforms to FITS standard                      BITPIX  =                    8 / array data type                                NAXIS   =                    0 / number of array dimensions                     EXTEND  =                    T                                                  TILEID  =                    6 / DESI tile ID                                   EXPID   =                    3 / DESI exposure ID                               AIRMASS =                  1.0 / Airmass at middle of exposure                  DATE-OBS= '2019-01-01T22:00:00' / Start of exposure                             TELDEC  =                14.84 / Telescope pointing dec [degrees]               PROGRAM = 'dark    '                                                            NIGHT   = '20190101'           / Night of observation YEARMMDD                  TELRA   =               333.22 / Telescope pointing RA [degrees]                EXPTIME =                 1000 / Exposure time [sec]                            FLAVOR  = 'science '           / Flavor [arc, flat, science, ...]               DOSVER  = 'SIM     '                                                            FEEVER  = 'SIM     '                                                            DETECTOR= 'SIM     '                                                            DEPNAM00= 'python  '                                                            DEPVER00= '3.5.3   '                                                            DEPNAM01= 'numpy   '                                                            DEPVER01= '1.11.3  '                                                            DEPNAM02= 'scipy   '                                                            DEPVER02= '0.18.1  '                                                            DEPNAM03= 'astropy '                                                            DEPVER03= '1.3.3   '                                                            DEPNAM04= 'yaml    '                                                            DEPVER04= '3.12    '                                                            DEPNAM05= 'desiutil'                                                            DEPVER05= '1.9.6.dev526'                                                        DEPNAM06= 'desispec'                                                            DEPVER06= '0.15.2.dev1569'                                                      DEPNAM07= 'desitarget'                                                          DEPVER07= '0.14.0.dev1118'                                                      DEPNAM08= 'desimodel'                                                           DEPVER08= '0.7.0.dev204'                                                        DEPNAM09= 'desisim '                                                            DEPVER09= '0.20.0.dev942'                                                       DEPNAM10= 'specter '                                                            DEPVER10= '0.7.0dev1'                                                           DEPNAM11= 'speclite'                                                            DEPVER11= '0.5     '                                                            END                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              \n",
      "\n",
      "------------------------------------------------------------------------------------------\n",
      "SIMPLE  =                    T / conforms to FITS standard                      BITPIX  =                    8 / array data type                                NAXIS   =                    0 / number of array dimensions                     EXTEND  =                    T                                                  CHECKSUM= 'C4fFD4Z9C4fCC4Z9'   / HDU checksum updated 2017-07-31T12:21:56       DATASUM = '0       '           / data unit checksum updated 2017-07-31T12:21:56 END                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              \n",
      "\n",
      "------------------------------------------------------------------------------------------\n",
      "SIMPLE  =                    T / conforms to FITS standard                      BITPIX  =                  -32 / array data type                                NAXIS   =                    2 / number of array dimensions                     NAXIS1  =                 2701                                                  NAXIS2  =                  500                                                  EXTEND  =                    T                                                  EXTNAME = 'FIBERFLAT'          / name of this binary table extension            FLAVOR  = 'flat    '           / Flavor [arc, flat, science, ...]               PROGRAM = 'calib   '                                                            NIGHT   = '20190101'           / Night of observation YEARMMDD                  AIRMASS =                  1.0 / Airmass at middle of exposure                  EXPID   =                    1 / DESI exposure ID                               EXPTIME =                   10 / Exposure time [sec]                            TELRA   =               335.03 / Telescope pointing RA [degrees]                TILEID  =                    4 / DESI tile ID                                   TELDEC  =                19.88 / Telescope pointing dec [degrees]               DATE-OBS= '2019-01-01T22:00:00' / Start of exposure                             DOSVER  = 'SIM     '                                                            CAMERA  = 'b0      '                                                            FEEVER  = 'SIM     '                                                            DETECTOR= 'SIM     '                                                            GAIN1   =                  1.0                                                  GAIN2   =                  1.0                                                  GAIN3   =                  1.0                                                  GAIN4   =                  1.0                                                  RDNOISE1=                  3.0                                                  RDNOISE2=                  3.0                                                  RDNOISE3=                  3.0                                                  RDNOISE4=                  3.0                                                  PRESEC1 = '[1:4,1:2048]'                                                        DATASEC1= '[5:2052,1:2048]'                                                     BIASSEC1= '[2053:2102,1:2048]'                                                  CCDSEC1 = '[1:2048,1:2048]'                                                     PRESEC2 = '[4201:4204,1:2048]'                                                  DATASEC2= '[2153:4200,1:2048]'                                                  BIASSEC2= '[2103:2152,1:2048]'                                                  CCDSEC2 = '[2049:4096,1:2048]'                                                  PRESEC3 = '[1:4,2049:4096]'                                                     DATASEC3= '[5:2052,2049:4096]'                                                  BIASSEC3= '[2053:2102,2049:4096]'                                               CCDSEC3 = '[1:2048,2049:4096]'                                                  PRESEC4 = '[4201:4204,2049:4096]'                                               DATASEC4= '[2153:4200,2049:4096]'                                               BIASSEC4= '[2103:2152,2049:4096]'                                               CCDSEC4 = '[2049:4096,2049:4096]'                                               INHERIT =                    T                                                  DEPNAM00= 'python  '                                                            DEPVER00= '3.5.3   '                                                            DEPNAM01= 'numpy   '                                                            DEPVER01= '1.11.3  '                                                            DEPNAM02= 'scipy   '                                                            DEPVER02= '0.19.0  '                                                            DEPNAM03= 'astropy '                                                            DEPVER03= '1.3.3   '                                                            DEPNAM04= 'yaml    '                                                            DEPVER04= '3.12    '                                                            DEPNAM05= 'desiutil'                                                            DEPVER05= '1.9.6.dev526'                                                        DEPNAM06= 'desispec'                                                            DEPVER06= '0.15.1.dev1525'                                                      DEPNAM07= 'desitarget'                                                          DEPVER07= '0.14.0.dev1118'                                                      DEPNAM08= 'desimodel'                                                           DEPVER08= '0.7.0.dev204'                                                        DEPNAM09= 'desisim '                                                            DEPVER09= '0.20.0.dev942'                                                       DEPNAM10= 'specter '                                                            DEPVER10= '0.7.0dev1'                                                           DEPNAM11= 'speclite'                                                            DEPVER11= '0.5     '                                                            OVERSCN1=       111.4876171875                                                  OBSRDN1 =    3.015872958689156                                                  OVERSCN2=     117.178037109375                                                  OBSRDN2 =    3.028112366245134                                                  OVERSCN3=         133.89359375                                                  OBSRDN3 =    3.012665994933165                                                  OVERSCN4=      166.15365234375                                                  OBSRDN4 =    3.011381794220364                                                  CHECKSUM= '9UI1DR909RG0AR90'   / HDU checksum updated 2017-08-01T12:56:39       DATASUM = '3897283379'         / data unit checksum updated 2017-08-01T12:56:39 NSPEC   =                  500 / Number of spectra                              WAVEMIN =               3570.0 / First wavelength [Angstroms]                   WAVEMAX =               5730.0 / Last wavelength [Angstroms]                    WAVESTEP=                  0.8 / Wavelength step size [Angstroms]               SPECTER = '0.7.0dev1'          / https://github.com/desihub/specter             IN_PSF  = '...edison/desimodel/data/specpsf/psf-b.fits' / Input spectral PSF    IN_IMG  = '' / Input image                                                      FIBERMIN=                    0                                                  CHI2PDF =    1.296712143514523                                                  DEPNAM12= 'matplotlib'                                                          DEPVER12= '2.0.2   '                                                            END                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              \n",
      "\n",
      "------------------------------------------------------------------------------------------\n",
      "SIMPLE  =                    T / conforms to FITS standard                      BITPIX  =                  -32 / array data type                                NAXIS   =                    2 / number of array dimensions                     NAXIS1  =                 2701                                                  NAXIS2  =                  500                                                  EXTEND  =                    T                                                  EXTNAME = 'FLUX    '           / name of this binary table extension            TILEID  =                    6 / DESI tile ID                                   EXPID   =                    3 / DESI exposure ID                               AIRMASS =                  1.0 / Airmass at middle of exposure                  DATE-OBS= '2019-01-01T22:00:00' / Start of exposure                             TELDEC  =                14.84 / Telescope pointing dec [degrees]               PROGRAM = 'dark    '                                                            NIGHT   = '20190101'           / Night of observation YEARMMDD                  TELRA   =               333.22 / Telescope pointing RA [degrees]                EXPTIME =                 1000 / Exposure time [sec]                            FLAVOR  = 'science '           / Flavor [arc, flat, science, ...]               DOSVER  = 'SIM     '                                                            CAMERA  = 'b0      '                                                            FEEVER  = 'SIM     '                                                            DETECTOR= 'SIM     '                                                            GAIN1   =                  1.0                                                  GAIN2   =                  1.0                                                  GAIN3   =                  1.0                                                  GAIN4   =                  1.0                                                  RDNOISE1=                  3.0                                                  RDNOISE2=                  3.0                                                  RDNOISE3=                  3.0                                                  RDNOISE4=                  3.0                                                  PRESEC1 = '[1:4,1:2048]'                                                        DATASEC1= '[5:2052,1:2048]'                                                     BIASSEC1= '[2053:2102,1:2048]'                                                  CCDSEC1 = '[1:2048,1:2048]'                                                     PRESEC2 = '[4201:4204,1:2048]'                                                  DATASEC2= '[2153:4200,1:2048]'                                                  BIASSEC2= '[2103:2152,1:2048]'                                                  CCDSEC2 = '[2049:4096,1:2048]'                                                  PRESEC3 = '[1:4,2049:4096]'                                                     DATASEC3= '[5:2052,2049:4096]'                                                  BIASSEC3= '[2053:2102,2049:4096]'                                               CCDSEC3 = '[1:2048,2049:4096]'                                                  PRESEC4 = '[4201:4204,2049:4096]'                                               DATASEC4= '[2153:4200,2049:4096]'                                               BIASSEC4= '[2103:2152,2049:4096]'                                               CCDSEC4 = '[2049:4096,2049:4096]'                                               INHERIT =                    T                                                  DEPNAM00= 'python  '                                                            DEPVER00= '3.5.4   '                                                            DEPNAM01= 'numpy   '                                                            DEPVER01= '1.12.1  '                                                            DEPNAM02= 'scipy   '                                                            DEPVER02= '1.0.0   '                                                            DEPNAM03= 'astropy '                                                            DEPVER03= '2.0     '                                                            DEPNAM04= 'yaml    '                                                            DEPVER04= '3.12    '                                                            DEPNAM05= 'desiutil'                                                            DEPVER05= '1.9.9.dev596'                                                        DEPNAM06= 'desispec'                                                            DEPVER06= '0.17.1  '                                                            DEPNAM07= 'desitarget'                                                          DEPVER07= '0.14.0.dev1118'                                                      DEPNAM08= 'desimodel'                                                           DEPVER08= '0.7.0.dev204'                                                        DEPNAM09= 'desisim '                                                            DEPVER09= '0.20.0.dev942'                                                       DEPNAM10= 'specter '                                                            DEPVER10= '0.7.0dev1'                                                           DEPNAM11= 'speclite'                                                            DEPVER11= '0.5     '                                                            OVERSCN1=       169.9040234375                                                  OBSRDN1 =    3.013081122158177                                                  OVERSCN2=      191.38392578125                                                  OBSRDN2 =    3.016938183329346                                                  OVERSCN3=     180.222470703125                                                  OBSRDN3 =    3.022291531697968                                                  OVERSCN4=     115.466201171875                                                  OBSRDN4 =     3.02739602819368                                                  NSPEC   =                  500 / Number of spectra                              WAVEMIN =               3570.0 / First wavelength [Angstroms]                   WAVEMAX =               5730.0 / Last wavelength [Angstroms]                    WAVESTEP=                  0.8 / Wavelength step size [Angstroms]               FIBERMIN=                    0                                                  DEPNAM12= 'matplotlib'                                                          DEPVER12= '2.1.2   '                                                            CHECKSUM= 'knmgknmgknmgknmg'   / HDU checksum updated 2018-01-30T11:42:14       DATASUM = '3763681958'         / data unit checksum updated 2018-01-30T11:42:14 END                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              \n",
      "\n",
      "------------------------------------------------------------------------------------------\n"
     ]
    }
   ],
   "source": [
    "for x in [f1,f2,f3,f4]:\n",
    "    print x[0].header, '\\n\\n'+'---'*30"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "58484.9166667\n"
     ]
    }
   ],
   "source": [
    "from astropy.io import fits\n",
    "from astropy.time import Time\n",
    "\n",
    "# open file\n",
    "f = fits.open('data/20190101/desi-00000003.fits.fz')\n",
    "\n",
    "# read the time in isot\n",
    "time = f[0].header['DATE-OBS']\n",
    "\n",
    "# declare the format\n",
    "t = Time(time, format='isot', scale='utc')\n",
    "\n",
    "# Convert to MJD\n",
    "MJD =  t.mjd\n",
    "\n",
    "print(MJD)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "51544.0\n"
     ]
    }
   ],
   "source": [
    "time\n",
    "t2 = Time('2000-01-01T00:00:00', format='isot', scale='utc')\n",
    "print t2.mjd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}