Deblending of extended objects in HyperSuprimeCam imaging.

In [1]:

    

from astropy.table import Table


    

In [6]:

    

database = '/Users/rs548/GitHub/dmu_products'


    

In [7]:

    

cosmos = Table.read("{}/dmu0/dmu0_2MASS-extended-sources/data/2MASS-XSC_COSMOS.fits".format(database))


    

In [8]:

    

cosmos.colnames


    

    
Out[8]:





['field',
 'id2masx',
 'raj2000',
 'dej2000',
 'supradeg',
 'supdedeg',
 'density',
 'r_k20e',
 'j_k20e',
 'e_j_k20e',
 'f_j_k20e',
 'h_k20e',
 'e_h_k20e',
 'f_h_k20e',
 'k_k20e',
 'e_k_k20e',
 'f_k_k20e',
 'kb_a',
 'kpa',
 'sb_a',
 'spa',
 'r_ext',
 'j_ext',
 'e_j_ext',
 'h_ext',
 'e_h_ext',
 'k_ext',
 'e_k_ext',
 'cc']

In [10]:

    

cosmos.sort('r_k20e')


    

In [13]:

    

f = open('data/cosmos_extended.txt','w')
f.write('#?         filter  ra        dec       sw      sh      image    mask     variance # column descriptor  \n')
for galaxy in cosmos:
    if galaxy['r_k20e'] > 10.:
        f.write('   HSC-G   {}        {}        50asec  50asec  true     true     true \n'.format(galaxy['raj2000'],galaxy['dej2000']))
        f.write('   HSC-R   {}        {}        50asec  50asec  true     true     true \n'.format(galaxy['raj2000'],galaxy['dej2000']))
        f.write('   HSC-I   {}        {}        50asec  50asec  true     true     true \n'.format(galaxy['raj2000'],galaxy['dej2000']))
        f.write('   HSC-Z   {}        {}        50asec  50asec  true     true     true \n'.format(galaxy['raj2000'],galaxy['dej2000']))
        f.write('   HSC-Y   {}        {}        50asec  50asec  true     true     true \n'.format(galaxy['raj2000'],galaxy['dej2000']))
    
f.close()