# Program to convert a FITS file into an ASCII array : x y z v

import pyfits
import numpy
import math
import os
from math import *

# Change to directory where script is located
abspath = os.path.abspath(__file__)
dname = os.path.dirname(abspath)
os.chdir(dname)

infile = 'M33AndMilkyWay.fits'

outfilebase=str(infile.split('.fits')[0])

FitsFile = pyfits.open(infile)
image = FitsFile[0].data

sizez = image.shape[0]
sizey = image.shape[1]
sizex = image.shape[2]

for zp in range(0,sizez):
        outfilename = outfilebase+'_'+str(zp).zfill(3)+'.txt'

        output=open(outfilename,'w')

        print zp

        for xp in range(0,sizex):
	
                for yp in range(0,sizey):
		
                        stuff = str(xp)+' '+str(yp)+' '+str(zp)+' '+str(image[zp,yp,xp])+'\n'
			
			output.write(stuff)


        output.close()
		
FitsFile.close()