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

def find_by_bestObjID(bestObjID):
  hdulist = fits.open('imbh_so.fits')
  data=hdulist[1].data
  for row in data:
     if row['bestObjID']==bestObjID:
       return [row['ra'],row['dec']]
  hdulist.close()
  
  
f = open('list_fits_files_subaru.csv', 'r')

for line in f.readlines():
    ln=line
    data=line.split(';')
    hdulist = fits.open(data[1])
    oid = data[0].strip()
    coords = find_by_bestObjID(int(oid))
    w = WCS(hdulist[0].header)
    x, y = w.wcs_world2pix(coords[0], coords[1], 0)
    print x,y