import os
import csv
from astropy.io import fits
from astropy.wcs import WCS

def find_by_bestObjID(bestObjID):
  hdulist = fits.open('garik7.fits')
  data=hdulist[1].data
  #f_lso = open('list_fits_files_CM.txt', 'w')

  for row in data:
     if row['bestObjID'] == bestObjID:
       return [row['ra'],row['dec']]
  #f_lso.close()
  hdulist.close()
    

def create_assoc_gal(name, ra, dec):
  w = WCS(name)
  x, y = w.wcs_world2pix(ra, dec, 0)
  f_assoc = open('assoc_galaxy.txt', 'w')
  f_assoc.write(str(x)+'	'+str(y) + '\n')
  f_assoc.close()
  
  
  
  


with open('list_fits_files_CM.csv','r') as csvfile:
    rows = csv.reader(csvfile, delimiter=',')
    for row in rows:
      coords = find_by_bestObjID(int(row[0]))
      print coords
      print row
      os.chdir(row[2])
      try:
        create_assoc_gal(row[3], coords[0], coords[1])
      except IOError:
          print 'nofile'
      os.chdir('../../../')