M33_align_all.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Notebook to align cube from different filters with each other using the SN3 filter as a reference. \n",
    "# Their is more emission line in SN3. Using it as a reference minimize the interpolation on the data. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Import Libraries\n",
    "from astropy.io import fits\n",
    "import numpy as np\n",
    "import pyfits\n",
    "from astropy.utils.data import get_pkg_data_filename\n",
    "from reproject import reproject_interp"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Define path and root name\n",
    "path = '/h/r-nepton/M33_Dhruv/'\n",
    "name = 'M33-SE'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "ename": "IOError",
     "evalue": "[Errno 2] No such file or directory: '/h/r-nepton/M33_Dhruv/M33-SE_SN3.1.0.ORCS/MAPS/M33-SE_SN3.1.0.LineMaps.map.6563.1x1.flux.fits'",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mIOError\u001b[0m                                   Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-3-98e86992c758>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m      2\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      3\u001b[0m \u001b[0;31m# Flux, Heigt, velocity\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 4\u001b[0;31m \u001b[0mha\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mhaheader\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mpyfits\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mgetdata\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0;34m'_SN3.1.0.ORCS/MAPS/'\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0;34m'_SN3.1.0.LineMaps.map.6563.1x1.flux.fits'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mheader\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m      5\u001b[0m \u001b[0mnii1\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mnii1header\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mpyfits\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mgetdata\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0;34m'_SN3.1.0.ORCS/MAPS/'\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0;34m'_SN3.1.0.LineMaps.map.6548.1x1.flux.fits'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mheader\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      6\u001b[0m \u001b[0mnii2\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mnii2header\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mpyfits\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mgetdata\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0;34m'_SN3.1.0.ORCS/MAPS/'\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0;34m'_SN3.1.0.LineMaps.map.6583.1x1.flux.fits'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mheader\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python2.7/site-packages/pyfits/convenience.pyc\u001b[0m in \u001b[0;36mgetdata\u001b[0;34m(filename, *args, **kwargs)\u001b[0m\n\u001b[1;32m    178\u001b[0m     \u001b[0mview\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mkwargs\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpop\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'view'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mNone\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    179\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 180\u001b[0;31m     \u001b[0mhdulist\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mextidx\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0m_getext\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfilename\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmode\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    181\u001b[0m     \u001b[0mhdu\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mhdulist\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mextidx\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    182\u001b[0m     \u001b[0mdata\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mhdu\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdata\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python2.7/site-packages/pyfits/convenience.pyc\u001b[0m in \u001b[0;36m_getext\u001b[0;34m(filename, mode, *args, **kwargs)\u001b[0m\n\u001b[1;32m    731\u001b[0m         \u001b[0;32mraise\u001b[0m \u001b[0mTypeError\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'extver alone cannot specify an extension.'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    732\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 733\u001b[0;31m     \u001b[0mhdulist\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfitsopen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfilename\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmode\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mmode\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    734\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    735\u001b[0m     \u001b[0;32mreturn\u001b[0m \u001b[0mhdulist\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mext\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python2.7/site-packages/pyfits/hdu/hdulist.pyc\u001b[0m in \u001b[0;36mfitsopen\u001b[0;34m(name, mode, memmap, save_backup, **kwargs)\u001b[0m\n\u001b[1;32m    122\u001b[0m         \u001b[0;32mraise\u001b[0m \u001b[0mValueError\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'Empty filename: %s'\u001b[0m \u001b[0;34m%\u001b[0m \u001b[0mrepr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    123\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 124\u001b[0;31m     \u001b[0;32mreturn\u001b[0m \u001b[0mHDUList\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfromfile\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmode\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmemmap\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msave_backup\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    125\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    126\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python2.7/site-packages/pyfits/hdu/hdulist.pyc\u001b[0m in \u001b[0;36mfromfile\u001b[0;34m(cls, fileobj, mode, memmap, save_backup, **kwargs)\u001b[0m\n\u001b[1;32m    264\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    265\u001b[0m         return cls._readfrom(fileobj=fileobj, mode=mode, memmap=memmap,\n\u001b[0;32m--> 266\u001b[0;31m                              save_backup=save_backup, **kwargs)\n\u001b[0m\u001b[1;32m    267\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    268\u001b[0m     \u001b[0;34m@\u001b[0m\u001b[0mclassmethod\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python2.7/site-packages/pyfits/hdu/hdulist.pyc\u001b[0m in \u001b[0;36m_readfrom\u001b[0;34m(cls, fileobj, data, mode, memmap, save_backup, **kwargs)\u001b[0m\n\u001b[1;32m    784\u001b[0m             \u001b[0;32mif\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfileobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0m_File\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    785\u001b[0m                 \u001b[0;31m# instantiate a FITS file object (ffo)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 786\u001b[0;31m                 \u001b[0mffo\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0m_File\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfileobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmode\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mmode\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmemmap\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mmemmap\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    787\u001b[0m             \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    788\u001b[0m                 \u001b[0mffo\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfileobj\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python2.7/site-packages/pyfits/file.pyc\u001b[0m in \u001b[0;36m__init__\u001b[0;34m(self, fileobj, mode, memmap, clobber)\u001b[0m\n\u001b[1;32m    130\u001b[0m             \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_open_fileobj\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfileobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmode\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mclobber\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    131\u001b[0m         \u001b[0;32melif\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfileobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mstring_types\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 132\u001b[0;31m             \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_open_filename\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfileobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmode\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mclobber\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    133\u001b[0m         \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    134\u001b[0m             \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_open_filelike\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfileobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmode\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mclobber\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python2.7/site-packages/pyfits/file.pyc\u001b[0m in \u001b[0;36m_open_filename\u001b[0;34m(self, filename, mode, clobber)\u001b[0m\n\u001b[1;32m    485\u001b[0m             \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_file\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mbz2\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mBZ2File\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mbzip2_mode\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    486\u001b[0m         \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 487\u001b[0;31m             \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_file\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfileobj_open\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mPYFITS_MODES\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mmode\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    488\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    489\u001b[0m         \u001b[0;31m# Make certain we're back at the beginning of the file\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/local/lib/python2.7/site-packages/pyfits/util.pyc\u001b[0m in \u001b[0;36mfileobj_open\u001b[0;34m(filename, mode)\u001b[0m\n\u001b[1;32m    752\u001b[0m     \"\"\"\n\u001b[1;32m    753\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 754\u001b[0;31m     \u001b[0;32mreturn\u001b[0m \u001b[0mopen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfilename\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmode\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    755\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    756\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mIOError\u001b[0m: [Errno 2] No such file or directory: '/h/r-nepton/M33_Dhruv/M33-SE_SN3.1.0.ORCS/MAPS/M33-SE_SN3.1.0.LineMaps.map.6563.1x1.flux.fits'"
     ]
    }
   ],
   "source": [
    "# Reading files and removing Nans by spatial interpolations\n",
    "\n",
    "# Flux, Heigt, velocity\n",
    "ha, haheader = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.6563.1x1.flux.fits', header=True)\n",
    "nii1, nii1header = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.6548.1x1.flux.fits', header=True)\n",
    "nii2, nii2header = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.6583.1x1.flux.fits', header=True)\n",
    "sii1, sii1header = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.6716.1x1.flux.fits', header=True)\n",
    "sii2, sii2header = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.6731.1x1.flux.fits', header=True)\n",
    "cont3, cont3header = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.all.1x1.height.fits', header=True)\n",
    "v3, v3header = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.6563.1x1.velocity.fits', header=True)\n",
    "\n",
    "hb, hbheader = pyfits.getdata(path+name+'_SN2.1.0.ORCS/MAPS/'+name+'_SN2.1.0.LineMaps.map.4861.1x1.flux.fits', header=True)\n",
    "oiii2, oiii2header = pyfits.getdata(path+name+'_SN2.1.0.ORCS/MAPS/'+name+'_SN2.1.0.LineMaps.map.4959.1x1.flux.fits', header=True)\n",
    "oiii1, oiii1header = pyfits.getdata(path+name+'_SN2.1.0.ORCS/MAPS/'+name+'_SN2.1.0.LineMaps.map.5007.1x1.flux.fits', header=True)\n",
    "cont2, cont2header = pyfits.getdata(path+name+'_SN2.1.0.ORCS/MAPS/'+name+'_SN2.1.0.LineMaps.map.all.1x1.height.fits', header=True)\n",
    "v2, v2header = pyfits.getdata(path+name+'_SN2.1.0.ORCS/MAPS/'+name+'_SN2.1.0.LineMaps.map.4861.1x1.velocity.fits', header=True)\n",
    "\n",
    "oii, oiiheader = pyfits.getdata(path+name+'_SN1.1.0.ORCS/MAPS/'+name+'_SN1.1.0.LineMaps.map.all.1x1.flux.fits', header=True)\n",
    "cont1, cont1header = pyfits.getdata(path+name+'_SN1.1.0.ORCS/MAPS/'+name+'SN1.1.0.LineMaps.map.all.1x1.height.fits', header=True)\n",
    "v1, v1header = pyfits.getdata(path+name+'_SN1.1.0.ORCS/MAPS/'+name+'SN1.1.0.LineMaps.map.all.1x1.velocity.fits', header=True)\n",
    "\n",
    "# Flux, Heigt, velocity errors\n",
    "haerr, haheadererr = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.6563.1x1.flux.fits', header=True)\n",
    "nii1err, nii1headererr = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.6548.1x1.flux.fits', header=True)\n",
    "nii2err, nii2headererr = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.6583.1x1.flux.fits', header=True)\n",
    "sii1err, sii1headererr = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.6716.1x1.flux.fits', header=True)\n",
    "sii2err, sii2headererr = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.6731.1x1.flux.fits', header=True)\n",
    "cont3err, cont3headererr = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.all.1x1.height.fits', header=True)\n",
    "v3err, v3headererr = pyfits.getdata(path+name+'_SN3.1.0.ORCS/MAPS/'+name+'_SN3.1.0.LineMaps.map.6563.1x1.velocity.fits', header=True)\n",
    "\n",
    "hberr, hbheadererr = pyfits.getdata(path+name+'_SN2.1.0.ORCS/MAPS/'+name+'_SN2.1.0.LineMaps.map.4861.1x1.flux.fits', header=True)\n",
    "oiii2err, oiii2headererr = pyfits.getdata(path+name+'_SN2.1.0.ORCS/MAPS/'+name+'_SN2.1.0.LineMaps.map.4959.1x1.flux.fits', header=True)\n",
    "oiii1err, oiii1headererr = pyfits.getdata(path+name+'_SN2.1.0.ORCS/MAPS/'+name+'_SN2.1.0.LineMaps.map.5007.1x1.flux.fits', header=True)\n",
    "cont2err, cont2headererr = pyfits.getdata(path+name+'_SN2.1.0.ORCS/MAPS/'+name+'_SN2.1.0.LineMaps.map.all.1x1.height.fits', header=True)\n",
    "v2err, v2headererr = pyfits.getdata(path+name+'_SN2.1.0.ORCS/MAPS/'+name+'_SN2.1.0.LineMaps.map.4861.1x1.velocity.fits', header=True)\n",
    "\n",
    "cont1err, cont1headererr = pyfits.getdata(path+name+'_SN1.1.0.ORCS/MAPS/'+name+'SN1.1.0.LineMaps.map.all.1x1.height.fits', header=True)\n",
    "v1err, v1headererr = pyfits.getdata(path+name+'_SN1.1.0.ORCS/MAPS/'+name+'SN1.1.0.LineMaps.map.all.1x1.velocity.fits', header=True)\n",
    "\n",
    "# Remove the contour lines bad values\n",
    "ha[0:1,:] = 1e-22\n",
    "ha[2062:2063,:] = 1e-22\n",
    "ha[:,0:1] = 1e-22\n",
    "ha[:,2046:2047] = 1e-22\n",
    "hb[0:1,:] = 1e-22\n",
    "hb[2062:2063,:] = 1e-22\n",
    "hb[:,0:1] = 1e-22\n",
    "hb[:,2046:2047] = 1e-22\n",
    "nii1[0:1,:] = 1e-22\n",
    "nii1[2062:2063,:] = 1e-22\n",
    "nii1[:,0:1] = 1e-22\n",
    "nii1[:,2046:2047] = 1e-22\n",
    "nii2[0:1,:] = 1e-22\n",
    "nii2[2062:2063,:] = 1e-22\n",
    "nii2[:,0:1] = 1e-22\n",
    "nii2[:,2046:2047] = 1e-22\n",
    "sii1[0:1,:] = 1e-22\n",
    "sii1[2062:2063,:] = 1e-22\n",
    "sii1[:,0:1] = 1e-22\n",
    "sii1[:,2046:2047] = 1e-22\n",
    "sii2[0:1,:] = 1e-22\n",
    "sii2[2062:2063,:] = 1e-22\n",
    "sii2[:,0:1] = 1e-22\n",
    "sii2[:,2046:2047] = 1e-22\n",
    "oiii1[0:1,:] = 1e-22\n",
    "oiii1[2062:2063,:] = 1e-22\n",
    "oiii1[:,0:1] = 1e-22\n",
    "oiii1[:,2046:2047] = 1e-22\n",
    "oiii2[0:1,:] = 1e-22\n",
    "oiii2[2062:2063,:] = 1e-22\n",
    "oiii2[:,0:1] = 1e-22\n",
    "oiii2[:,2046:2047] = 1e-22\n",
    "oii[0:1,:] = 1e-22\n",
    "oii[2062:2063,:] = 1e-22\n",
    "oii[:,0:1] = 1e-22\n",
    "oii[:,2046:2047] = 1e-22\n",
    "cont1[0:1,:] = 1e-22\n",
    "cont1[2062:2063,:] = 1e-22\n",
    "cont1[:,0:1] = 1e-22\n",
    "cont1[:,2046:2047] = 1e-22\n",
    "cont2[0:1,:] = 1e-22\n",
    "cont2[2062:2063,:] = 1e-22\n",
    "cont2[:,0:1] = 1e-22\n",
    "cont2[:,2046:2047] = 1e-22\n",
    "cont3[0:1,:] = 1e-22\n",
    "cont3[2062:2063,:] = 1e-22\n",
    "cont3[:,0:1] = 1e-22\n",
    "cont3[:,2046:2047] = 1e-22\n",
    "cont1[(cont1 < 0)] = 0\n",
    "cont2[(cont2 < 0)] = 0\n",
    "cont3[(cont3 < 0)] = 0\n",
    "\n",
    "haerr[0:1,:] = 1e-22\n",
    "haerr[2062:2063,:] = 1e-22\n",
    "haerr[:,0:1] = 1e-22\n",
    "haerr[:,2046:2047] = 1e-22\n",
    "hberr[0:1,:] = 1e-22\n",
    "hberr[2062:2063,:] = 1e-22\n",
    "hberr[:,0:1] = 1e-22\n",
    "hberr[:,2046:2047] = 1e-22\n",
    "nierri1[0:1,:] = 1e-22\n",
    "nii1err[2062:2063,:] = 1e-22\n",
    "nii1err[:,0:1] = 1e-22\n",
    "nii1err[:,2046:2047] = 1e-22\n",
    "nii2err[0:1,:] = 1e-22\n",
    "nii2err[2062:2063,:] = 1e-22\n",
    "nii2err[:,0:1] = 1e-22\n",
    "nii2err[:,2046:2047] = 1e-22\n",
    "sii1err[0:1,:] = 1e-22\n",
    "sii1err[2062:2063,:] = 1e-22\n",
    "sii1err[:,0:1] = 1e-22\n",
    "sii1err[:,2046:2047] = 1e-22\n",
    "sii2err[0:1,:] = 1e-22\n",
    "sii2err[2062:2063,:] = 1e-22\n",
    "sii2err[:,0:1] = 1e-22\n",
    "sii2err[:,2046:2047] = 1e-22\n",
    "oiii1err[0:1,:] = 1e-22\n",
    "oiii1err[2062:2063,:] = 1e-22\n",
    "oiii1err[:,0:1] = 1e-22\n",
    "oiii1err[:,2046:2047] = 1e-22\n",
    "oiii2err[0:1,:] = 1e-22\n",
    "oiii2err[2062:2063,:] = 1e-22\n",
    "oiii2err[:,0:1] = 1e-22\n",
    "oiii2err[:,2046:2047] = 1e-22\n",
    "oiierr[0:1,:] = 1e-22\n",
    "oiierr[2062:2063,:] = 1e-22\n",
    "oiierr[:,0:1] = 1e-22\n",
    "oiierr[:,2046:2047] = 1e-22\n",
    "cont1err[0:1,:] = 1e-22\n",
    "cont1err[2062:2063,:] = 1e-22\n",
    "cont1err[:,0:1] = 1e-22\n",
    "cont1err[:,2046:2047] = 1e-22\n",
    "cont2err[0:1,:] = 1e-22\n",
    "cont2err[2062:2063,:] = 1e-22\n",
    "cont2err[:,0:1] = 1e-22\n",
    "cont2err[:,2046:2047] = 1e-22\n",
    "cont3err[0:1,:] = 1e-22\n",
    "cont3err[2062:2063,:] = 1e-22\n",
    "cont3err[:,0:1] = 1e-22\n",
    "cont3err[:,2046:2047] = 1e-22\n",
    "cont1err[(cont1 < 0)] = 0\n",
    "cont2err[(cont2 < 0)] = 0\n",
    "cont3err[(cont3 < 0)] = 0\n",
    "\n",
    "# Find the indices of the Nans pixels\n",
    "indice_ha = np.where(np.isfinite(ha) == False)\n",
    "indice_hb = np.where(np.isfinite(hb) == False)\n",
    "indice_nii1 = np.where(np.isfinite(nii1) == False)\n",
    "indice_nii2 = np.where(np.isfinite(nii2) == False)\n",
    "indice_sii1 = np.where(np.isfinite(sii1) == False)\n",
    "indice_sii2 = np.where(np.isfinite(sii2) == False)\n",
    "indice_oiii1 = np.where(np.isfinite(oiii1) == False)\n",
    "indice_oiii2 = np.where(np.isfinite(oiii2) == False)\n",
    "indice_oii = np.where(np.isfinite(oii) == False)\n",
    "indice_cont1 = np.where((np.isfinite(cont1) == False))\n",
    "indice_cont2 = np.where((np.isfinite(cont2) == False))\n",
    "indice_cont3 = np.where((np.isfinite(cont3) == False))\n",
    "indice_v1 = np.where((np.isfinite(v1) == False))\n",
    "indice_v2 = np.where((np.isfinite(v2) == False))\n",
    "indice_v3 = np.where((np.isfinite(v3) == False))\n",
    "\n",
    "\n",
    "indice_haerr = np.where(np.isfinite(haerr) == False)\n",
    "indice_hberr = np.where(np.isfinite(hberr) == False)\n",
    "indice_nii1err = np.where(np.isfinite(nii1err) == False)\n",
    "indice_nii2err = np.where(np.isfinite(nii2err) == False)\n",
    "indice_sii1err = np.where(np.isfinite(sii1err) == False)\n",
    "indice_sii2err = np.where(np.isfinite(sii2err) == False)\n",
    "indice_oiii1err = np.where(np.isfinite(oiii1err) == False)\n",
    "indice_oiii2err = np.where(np.isfinite(oiii2err) == False)\n",
    "indice_oiierr = np.where(np.isfinite(oiierr) == False)\n",
    "indice_cont1err = np.where((np.isfinite(cont1err) == False))\n",
    "indice_cont2err = np.where((np.isfinite(cont2err) == False))\n",
    "indice_cont3err = np.where((np.isfinite(cont3err) == False))\n",
    "indice_v1err = np.where((np.isfinite(v1err) == False))\n",
    "indice_v2err = np.where((np.isfinite(v2err) == False))\n",
    "indice_v3err = np.where((np.isfinite(v3err) == False))\n",
    "\n",
    "# Mask the Nans with the mean of the surrounding [-2:+2]\n",
    "for i in range(0,len(indice_ha[0])):  \n",
    "    ha[indice_ha[0][i],indice_ha[1][i]] = np.nanmean(ha[(indice_ha[0][i]-2):(indice_ha[0][i]+2),(indice_ha[1][i]-2):(indice_ha[1][i]+2)])    \n",
    "\n",
    "for i in range(0,len(indice_hb[0])):\n",
    "    hb[indice_hb[0][i],indice_hb[1][i]] = np.nanmean(hb[(indice_hb[0][i]-2):(indice_hb[0][i]+2),(indice_hb[1][i]-2):(indice_hb[1][i]+2)])\n",
    " \n",
    "for i in range(0,len(indice_nii1[0])):\n",
    "    nii1[indice_nii1[0][i],indice_nii1[1][i]] = np.nanmean(nii1[(indice_nii1[0][i]-2):(indice_nii1[0][i]+2),(indice_nii1[1][i]-2):(indice_nii1[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_nii2[0])):\n",
    "    nii2[indice_nii2[0][i],indice_nii2[1][i]] = np.nanmean(nii2[(indice_nii2[0][i]-2):(indice_nii2[0][i]+2),(indice_nii2[1][i]-2):(indice_nii2[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_sii1[0])):\n",
    "    sii1[indice_sii1[0][i],indice_sii1[1][i]] = np.nanmean(sii1[(indice_sii1[0][i]-2):(indice_sii1[0][i]+2),(indice_sii1[1][i]-2):(indice_sii1[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_sii2[0])):\n",
    "    sii2[indice_sii2[0][i],indice_sii2[1][i]] = np.nanmean(sii2[(indice_sii2[0][i]-2):(indice_sii2[0][i]+2),(indice_sii2[1][i]-2):(indice_sii2[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_oiii1[0])):\n",
    "    oiii1[indice_oiii1[0][i],indice_oiii1[1][i]] = np.nanmean(oiii1[(indice_oiii1[0][i]-2):(indice_oiii1[0][i]+2),(indice_oiii1[1][i]-2):(indice_oiii1[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_oiii2[0])):\n",
    "    oiii2[indice_oiii2[0][i],indice_oiii2[1][i]] = np.nanmean(oiii2[(indice_oiii2[0][i]-2):(indice_oiii2[0][i]+2),(indice_oiii2[1][i]-2):(indice_oiii2[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_oii[0])):\n",
    "    oii[indice_oii[0][i],indice_oii[1][i]] = np.nanmean(oii[(indice_oii[0][i]-2):(indice_oii[0][i]+2),(indice_oii[1][i]-2):(indice_oii[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_cont1[0])):\n",
    "    cont1[indice_cont1[0][i],indice_cont1[1][i]] = np.nanmean(cont1[(indice_cont1[0][i]-2):(indice_cont1[0][i]+2),(indice_cont1[1][i]-2):(indice_cont1[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_cont2[0])):\n",
    "    cont2[indice_cont2[0][i],indice_cont2[1][i]] = np.nanmean(cont2[(indice_cont2[0][i]-2):(indice_cont2[0][i]+2),(indice_cont2[1][i]-2):(indice_cont2[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_cont3[0])):\n",
    "    cont3[indice_cont3[0][i],indice_cont3[1][i]] = np.nanmean(cont3[(indice_cont3[0][i]-2):(indice_cont3[0][i]+2),(indice_cont3[1][i]-2):(indice_cont3[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_v1[0])):\n",
    "    v1[indice_v1[0][i],indice_v1[1][i]] = np.nanmean(v1[(indice_v1[0][i]-2):(indice_v1[0][i]+2),(indice_v1[1][i]-2):(indice_v1[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_v2[0])):\n",
    "    v2[indice_v2[0][i],indice_v2[1][i]] = np.nanmean(v2[(indice_v2[0][i]-2):(indice_v2[0][i]+2),(indice_v2[1][i]-2):(indice_v2[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_v3[0])):\n",
    "    v3[indice_v3[0][i],indice_v3[1][i]] = np.nanmean(v3[(indice_v3[0][i]-2):(indice_v3[0][i]+2),(indice_v3[1][i]-2):(indice_v3[1][i]+2)])\n",
    "\n",
    "# Write NoNans files\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN3.1.0.LineMaps.map.6563.1x1.flux_nonans.fits', ha, header=haheader, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN3.1.0.LineMaps.map.6548.1x1.flux_nonans.fits', nii1, header=nii1header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN3.1.0.LineMaps.map.6583.1x1.flux_nonans.fits', nii2, header=nii2header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN3.1.0.LineMaps.map.6716.1x1.flux_nonans.fits', sii1, header=sii1header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN3.1.0.LineMaps.map.6731.1x1.flux_nonans.fits', sii2, header=sii2header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN2.1.0.LineMaps.map.4861.1x1.flux_nonans.fits', hb, header=hbheader, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN2.1.0.LineMaps.map.4959.1x1.flux_nonans.fits', oiii1, header=oiii1header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN2.1.0.LineMaps.map.5007.1x1.flux_nonans.fits', oiii2, header=oiii2header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN1.1.0.LineMaps.map.3727.1x1.flux_nonans.fits', oii, header=oiiheader, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN1.1.0.LineMaps.map.all.1x1.height_nonans.fits', cont1, header=cont1header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN2.1.0.LineMaps.map.all.1x1.height_nonans.fits', cont2, header=cont2header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN3.1.0.LineMaps.map.all.1x1.height_nonans.fits', cont3, header=cont3header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN1.1.0.LineMaps.map.all.1x1.velocity_nonans.fits', v1, header=v1header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN2.1.0.LineMaps.map.all.1x1.velocity_nonans.fits', v2, header=v2header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/'+name+'_SN3.1.0.LineMaps.map.all.1x1.velocity_nonans.fits', v3, header=v3header, clobber=True)\n",
    "\n",
    "# Propagate errors\n",
    "\n",
    "for i in range(0,len(indice_ha[0])):  \n",
    "    haerr[indice_haerr[0][i],indice_haerr[1][i]] = np.nanmean(haerr[(indice_haerr[0][i]-2):(indice_haerr[0][i]+2),(indice_haerr[1][i]-2):(indice_haerr[1][i]+2)])    \n",
    "\n",
    "for i in range(0,len(indice_hb[0])):\n",
    "    hberr[indice_hberr[0][i],indice_hberr[1][i]] = np.nanmean(hberr[(indice_hberr[0][i]-2):(indice_hberr[0][i]+2),(indice_hberr[1][i]-2):(indice_hberr[1][i]+2)])\n",
    " \n",
    "for i in range(0,len(indice_nii1[0])):\n",
    "    nii1err[indice_nii1err[0][i],indice_nii1err[1][i]] = np.nanmean(nii1err[(indice_nii1err[0][i]-2):(indice_nii1err[0][i]+2),(indice_nii1err[1][i]-2):(indice_nii1err[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_nii2[0])):\n",
    "    nii2err[indice_nii2err[0][i],indice_nii2err[1][i]] = np.nanmean(nii2err[(indice_nii2err[0][i]-2):(indice_nii2err[0][i]+2),(indice_nii2err[1][i]-2):(indice_nii2err[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_sii1[0])):\n",
    "    sii1err[indice_sii1err[0][i],indice_sii1err[1][i]] = np.nanmean(sii1err[(indice_sii1err[0][i]-2):(indice_sii1err[0][i]+2),(indice_sii1err[1][i]-2):(indice_sii1err[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_sii2[0])):\n",
    "    sii2err[indice_sii2err[0][i],indice_sii2err[1][i]] = np.nanmean(sii2err[(indice_sii2err[0][i]-2):(indice_sii2err[0][i]+2),(indice_sii2err[1][i]-2):(indice_sii2err[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_oiii1[0])):\n",
    "    oiii1err[indice_oiii1err[0][i],indice_oiii1err[1][i]] = np.nanmean(oiii1err[(indice_oiii1err[0][i]-2):(indice_oiii1err[0][i]+2),(indice_oiii1err[1][i]-2):(indice_oiii1err[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_oiii2[0])):\n",
    "    oiii2err[indice_oiii2err[0][i],indice_oiii2err[1][i]] = np.nanmean(oiii2err[(indice_oiii2err[0][i]-2):(indice_oiii2err[0][i]+2),(indice_oiii2err[1][i]-2):(indice_oiii2err[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_oii[0])):\n",
    "    oiierr[indice_oiierr[0][i],indice_oiierr[1][i]] = np.nanmean(oiierr[(indice_oiierr[0][i]-2):(indice_oiierr[0][i]+2),(indice_oiierr[1][i]-2):(indice_oiierr[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_cont1[0])):\n",
    "    cont1err[indice_cont1err[0][i],indice_cont1err[1][i]] = np.nanmean(cont1err[(indice_cont1err[0][i]-2):(indice_cont1err[0][i]+2),(indice_cont1err[1][i]-2):(indice_cont1err[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_cont2[0])):\n",
    "    cont2err[indice_cont2err[0][i],indice_cont2err[1][i]] = np.nanmean(cont2err[(indice_cont2err[0][i]-2):(indice_cont2err[0][i]+2),(indice_cont2err[1][i]-2):(indice_cont2err[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_cont3[0])):\n",
    "    cont3err[indice_cont3err[0][i],indice_cont3err[1][i]] = np.nanmean(cont3err[(indice_cont3err[0][i]-2):(indice_cont3err[0][i]+2),(indice_cont3err[1][i]-2):(indice_cont3err[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_v1[0])):\n",
    "    v1err[indice_v1err[0][i],indice_v1err[1][i]] = np.nanmean(v1err[(indice_v1err[0][i]-2):(indice_v1err[0][i]+2),(indice_v1err[1][i]-2):(indice_v1err[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_v2[0])):\n",
    "    v2err[indice_v2err[0][i],indice_v2err[1][i]] = np.nanmean(v2err[(indice_v2err[0][i]-2):(indice_v2err[0][i]+2),(indice_v2err[1][i]-2):(indice_v2err[1][i]+2)])\n",
    "\n",
    "for i in range(0,len(indice_v3[0])):\n",
    "    v3err[indice_v3err[0][i],indice_v3err[1][i]] = np.nanmean(v3err[(indice_v3err[0][i]-2):(indice_v3err[0][i]+2),(indice_v3err[1][i]-2):(indice_v3err[1][i]+2)])\n",
    "\n",
    "# Write Nonans files errors\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.6563.1x1.flux_nonans.fits', ha, header=haheader, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.6548.1x1.flux_nonans.fits', nii1, header=nii1header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.6583.1x1.flux_nonans.fits', nii2, header=nii2header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.6716.1x1.flux_nonans.fits', sii1, header=sii1header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.6731.1x1.flux_nonans.fits', sii2, header=sii2header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.4861.1x1.flux_nonans.fits', hb, header=hbheader, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.4959.1x1.flux_nonans.fits', oiii1, header=oiii1header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.5007.1x1.flux_nonans.fits', oiii2, header=oiii2header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN1.1.0.LineMaps.map.3727.1x1.flux_nonans.fits', oii, header=oiiheader, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN1.1.0.LineMaps.map.all.1x1.height_nonans.fits', cont1, header=cont1header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.all.1x1.height_nonans.fits', cont2, header=cont2header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.all.1x1.height_nonans.fits', cont3, header=cont3header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN1.1.0.LineMaps.map.all.1x1.velocity_nonans.fits', v1, header=v1header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.all.1x1.velocity_nonans.fits', v2, header=v2header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.all.1x1.velocity_nonans.fits', v3, header=v3header, clobber=True)\n",
    "\n",
    "# Write Nonans files error\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.6563.1x1.flux-err_nonans.fits', haerr, header=haheadererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.6548.1x1.flux-err_nonans.fits', nii1err, header=nii1headererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.6583.1x1.flux-err_nonans.fits', nii2err, header=nii2headererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.6716.1x1.flux-err_nonans.fits', sii1err, header=sii1headererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.6731.1x1.flux-err_nonans.fits', sii2err, header=sii2headererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.4861.1x1.flux-err_nonans.fits', hberr, header=hbheadererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.4959.1x1.flux-err_nonans.fits', oiii1err, header=oiii1headererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.5007.1x1.flux-err_nonans.fits', oiii2err, header=oiii2headererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN1.1.0.LineMaps.map.3727.1x1.flux-err_nonans.fits', oiierr, header=oiiheadererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN1.1.0.LineMaps.map.all.1x1.height-err_nonans.fits', cont1err, header=cont1headererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.all.1x1.height-err_nonans.fits', cont2err, header=cont2headererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.all.1x1.height-err_nonans.fits', cont3err, header=cont3headererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN1.1.0.LineMaps.map.all.1x1.velocity-err_nonans.fits', v1err, header=v1headererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.all.1x1.velocity-err_nonans.fits', v2err, header=v2headererr, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.all.1x1.velocity-err_nonans.fits', v3err, header=v3headererr, clobber=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/usr/local/lib/python2.7/site-packages/pyfits/file.py:372: UserWarning: Overwriting existing file '/Users/Laurie/M33-SE_SN3.1.0.ORCS/MAPS/M33-SE_SN3.1.0.LineMaps.map.6563.1x1.flux_align.fits'.\n",
      "  warnings.warn(\"Overwriting existing file %r.\" % self.name)\n",
      "/usr/local/lib/python2.7/site-packages/pyfits/file.py:372: UserWarning: Overwriting existing file '/Users/Laurie/M33-SE_SN3.1.0.ORCS/MAPS/M33-SE_SN3.1.0.LineMaps.map.6548.1x1.flux_align.fits'.\n",
      "  warnings.warn(\"Overwriting existing file %r.\" % self.name)\n",
      "/usr/local/lib/python2.7/site-packages/pyfits/file.py:372: UserWarning: Overwriting existing file '/Users/Laurie/M33-SE_SN3.1.0.ORCS/MAPS/M33-SE_SN3.1.0.LineMaps.map.6583.1x1.flux_align.fits'.\n",
      "  warnings.warn(\"Overwriting existing file %r.\" % self.name)\n",
      "/usr/local/lib/python2.7/site-packages/pyfits/file.py:372: UserWarning: Overwriting existing file '/Users/Laurie/M33-SE_SN3.1.0.ORCS/MAPS/M33-SE_SN3.1.0.LineMaps.map.6716.1x1.flux_align.fits'.\n",
      "  warnings.warn(\"Overwriting existing file %r.\" % self.name)\n",
      "/usr/local/lib/python2.7/site-packages/pyfits/file.py:372: UserWarning: Overwriting existing file '/Users/Laurie/M33-SE_SN3.1.0.ORCS/MAPS/M33-SE_SN3.1.0.LineMaps.map.6731.1x1.flux_align.fits'.\n",
      "  warnings.warn(\"Overwriting existing file %r.\" % self.name)\n",
      "/usr/local/lib/python2.7/site-packages/pyfits/file.py:372: UserWarning: Overwriting existing file '/Users/Laurie/M33-SE_SN3.1.0.ORCS/MAPS/M33-SE_SN3.1.0.LineMaps.map.all.1x1.height_align.fits'.\n",
      "  warnings.warn(\"Overwriting existing file %r.\" % self.name)\n"
     ]
    }
   ],
   "source": [
    "# Write the aligned SN3 maps (No modifications but keep the roots of the file an _align for simplicity)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.6563.1x1.flux_align.fits', ha, header=haheader, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.6548.1x1.flux_align.fits', nii1, header=nii1header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.6583.1x1.flux_align.fits', nii2, header=nii2header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.6716.1x1.flux_align.fits', sii1, header=sii1header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.6731.1x1.flux_align.fits', sii2, header=sii2header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.all.1x1.height_align.fits', cont3, header=cont3header, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.6563.1x1.velocity_align.fits', v3, header=v3header, clobber=True)\n",
    "\n",
    "# Write the aligned SN3 maps error (No modifications but keep the roots of the file an _align for simplicity)\n",
    "\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.6563.1x1.flux-err_align.fits', ha-err, header=haheader-err, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.6548.1x1.flux-err_align.fits', nii1-err, header=nii1header-err, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.6583.1x1.flux-err_align.fits', nii2-err, header=nii2header-err, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.6716.1x1.flux-err_align.fits', sii1-err, header=sii1header-err, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.6731.1x1.flux-err_align.fits', sii2-err, header=sii2header-err, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.all.1x1.height-err_align.fits', cont3-err, header=cont3header-err, clobber=True)\n",
    "pyfits.writeto(path+'Analysis/Aligned/'+name+'_SN3.1.0.LineMaps.map.6563.1x1.velocity-err_align.fits', v3-err, header=v3header-err, clobber=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "## Align the fits files.\n",
    "\n",
    "hdu1 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN3.1.0.LineMaps.map.6563.1x1.flux_nonans.fits'))[0]\n",
    "hdu2 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.4861.1x1.flux_nonans.fits'))[0]\n",
    "hdu9 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.4861.1x1.velocity_nonans.fits'))[0]\n",
    "hdu3 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.4959.1x1.flux_nonans.fits'))[0]\n",
    "hdu4 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.5007.1x1.flux_nonans.fits'))[0]\n",
    "hdu5 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN1.1.0.LineMaps.map.3727.1x1.flux_nonans.fits'))[0]\n",
    "hdu6 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.all.1x1.height_nonans.fits'))[0]\n",
    "hdu7 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN1.1.0.LineMaps.map.3727.1x1.height_nonans.fits'))[0]\n",
    "hdu8 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN1.1.0.LineMaps.map.3727.1x1.velocity_nonans.fits'))[0]\n",
    "\n",
    "hdu22 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.4861.1x1.flux-err_nonans.fits'))[0]\n",
    "hdu99 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.4861.1x1.velocityerr_nonans.fits'))[0]\n",
    "hdu33 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.4959.1x1.fluxerr_nonans.fits'))[0]\n",
    "hdu44 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.5007.1x1.fluxerr_nonans.fits'))[0]\n",
    "hdu55 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN1.1.0.LineMaps.map.3727.1x1.fluxerr_nonans.fits'))[0]\n",
    "hdu66 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN2.1.0.LineMaps.map.all.1x1.heighterr_nonans.fits'))[0]\n",
    "hdu77 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN1.1.0.LineMaps.map.3727.1x1.heighterr_nonans.fits'))[0]\n",
    "hdu88 = fits.open(get_pkg_data_filename(path+'Analysis/No_nans/'+name+'_SN1.1.0.LineMaps.map.3727.1x1.velocityerr_nonans.fits'))[0]\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n",
      "INFO: \n",
      "                Inconsistent SIP distortion information is present in the FITS header and the WCS object:\n",
      "                SIP coefficients were detected, but CTYPE is missing a \"-SIP\" suffix.\n",
      "                astropy.wcs is using the SIP distortion coefficients,\n",
      "                therefore the coordinates calculated here might be incorrect.\n",
      "\n",
      "                If you do not want to apply the SIP distortion coefficients,\n",
      "                please remove the SIP coefficients from the FITS header or the\n",
      "                WCS object.  As an example, if the image is already distortion-corrected\n",
      "                (e.g., drizzled) then distortion components should not apply and the SIP\n",
      "                coefficients should be removed.\n",
      "\n",
      "                While the SIP distortion coefficients are being applied here, if that was indeed the intent,\n",
      "                for consistency please append \"-SIP\" to the CTYPE in the FITS header or the WCS object.\n",
      "\n",
      "                 [astropy.wcs.wcs]\n"
     ]
    }
   ],
   "source": [
    "# Calculate the reprojection\n",
    "\n",
    "array2, footprint2 = reproject_interp(hdu2, hdu1.header, hdu_in=1)\n",
    "array3, footprint3 = reproject_interp(hdu3, hdu1.header, hdu_in=1)\n",
    "array4, footprint4 = reproject_interp(hdu4, hdu1.header, hdu_in=1)\n",
    "array5, footprint5 = reproject_interp(hdu5, hdu1.header, hdu_in=1)\n",
    "array6, footprint6 = reproject_interp(hdu6, hdu1.header, hdu_in=1)\n",
    "array7, footprint7 = reproject_interp(hdu7, hdu1.header, hdu_in=1)\n",
    "array8, footprint8 = reproject_interp(hdu8, hdu1.header, hdu_in=1)\n",
    "array9, footprint9 = reproject_interp(hdu9, hdu1.header, hdu_in=1)\n",
    "\n",
    "array22, footprint22 = reproject_interp(hdu22, hdu1.header, hdu_in=1)\n",
    "array33, footprint33 = reproject_interp(hdu33, hdu1.header, hdu_in=1)\n",
    "array44, footprint44 = reproject_interp(hdu44, hdu1.header, hdu_in=1)\n",
    "array55, footprint55 = reproject_interp(hdu55, hdu1.header, hdu_in=1)\n",
    "array66, footprint66 = reproject_interp(hdu66, hdu1.header, hdu_in=1)\n",
    "array77, footprint77 = reproject_interp(hdu77, hdu1.header, hdu_in=1)\n",
    "array88, footprint88 = reproject_interp(hdu88, hdu1.header, hdu_in=1)\n",
    "array99, footprint99 = reproject_interp(hdu99, hdu1.header, hdu_in=1)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: AstropyDeprecationWarning: \"clobber\" was deprecated in version 2.0 and will be removed in a future version. Use argument \"overwrite\" instead. [astropy.utils.decorators]\n",
      "WARNING: AstropyDeprecationWarning: \"clobber\" was deprecated in version 2.0 and will be removed in a future version. Use argument \"overwrite\" instead. [astropy.utils.decorators]\n",
      "WARNING: AstropyDeprecationWarning: \"clobber\" was deprecated in version 2.0 and will be removed in a future version. Use argument \"overwrite\" instead. [astropy.utils.decorators]\n",
      "WARNING: AstropyDeprecationWarning: \"clobber\" was deprecated in version 2.0 and will be removed in a future version. Use argument \"overwrite\" instead. [astropy.utils.decorators]\n",
      "WARNING: AstropyDeprecationWarning: \"clobber\" was deprecated in version 2.0 and will be removed in a future version. Use argument \"overwrite\" instead. [astropy.utils.decorators]\n",
      "WARNING: AstropyDeprecationWarning: \"clobber\" was deprecated in version 2.0 and will be removed in a future version. Use argument \"overwrite\" instead. [astropy.utils.decorators]\n",
      "WARNING: AstropyDeprecationWarning: \"clobber\" was deprecated in version 2.0 and will be removed in a future version. Use argument \"overwrite\" instead. [astropy.utils.decorators]\n",
      "WARNING: AstropyDeprecationWarning: \"clobber\" was deprecated in version 2.0 and will be removed in a future version. Use argument \"overwrite\" instead. [astropy.utils.decorators]\n"
     ]
    }
   ],
   "source": [
    "# Write aligned maps\n",
    "\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN2.1.0.ORCS/MAPS/M33-SE_SN2.1.0.LineMaps.map.4861.1x1.flux_align.fits', array2, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN2.1.0.ORCS/MAPS/M33-SE_SN2.1.0.LineMaps.map.all.1x1.velocity_align.fits', array9, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN2.1.0.ORCS/MAPS/M33-SE_SN2.1.0.LineMaps.map.4959.1x1.flux_align.fits', array3, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN2.1.0.ORCS/MAPS/M33-SE_SN2.1.0.LineMaps.map.5007.1x1.flux_align.fits', array4, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN1.1.0.ORCS/MAPS/M33-SE_SN1.1.0.LineMaps.map.3727.1x1.flux_align.fits', array5, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN1.1.0.ORCS/MAPS/M33-SE_SN2.1.0.LineMaps.map.all.1x1.height_align.fits', array6, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN1.1.0.ORCS/MAPS/M33-SE_SN1.1.0.LineMaps.map.all.1x1.height_align.fits', array7, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN1.1.0.ORCS/MAPS/M33-SE_SN1.1.0.LineMaps.map.all.1x1.velocity_align.fits', array8, hdu1.header, clobber=True)\n",
    "\n",
    "# Write aligned maps error\n",
    "\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN2.1.0.ORCS/MAPS/M33-SE_SN2.1.0.LineMaps.map.4861.1x1.flux-err_align.fits', array22, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN2.1.0.ORCS/MAPS/M33-SE_SN2.1.0.LineMaps.map.all.1x1.velocity-err_align.fits', array99, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN2.1.0.ORCS/MAPS/M33-SE_SN2.1.0.LineMaps.map.4959.1x1.flux-err_align.fits', array33, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN2.1.0.ORCS/MAPS/M33-SE_SN2.1.0.LineMaps.map.5007.1x1.flux-err_align.fits', array44, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN1.1.0.ORCS/MAPS/M33-SE_SN1.1.0.LineMaps.map.3727.1x1.flux-err_align.fits', array55, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN1.1.0.ORCS/MAPS/M33-SE_SN2.1.0.LineMaps.map.all.1x1.height-err_align.fits', array66, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN1.1.0.ORCS/MAPS/M33-SE_SN1.1.0.LineMaps.map.all.1x1.height-err_align.fits', array77, hdu1.header, clobber=True)\n",
    "fits.writeto(path+'Analysis/Aligned/'+name+'_SN1.1.0.ORCS/MAPS/M33-SE_SN1.1.0.LineMaps.map.all.1x1.velocity-err_align.fits', array88, hdu1.header, clobber=True)\n"
   ]
  }
 ],
 "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
}