Usamos un catálogo de galaxias identificadas en VVV en los tiles d010 y d0115 de Baravalle L.
Para saber donde estan ubicados los tiles usamos el mapa de VVV
En estos tiles encontraron 574 objetos con propiedades morfologicas, fotometricas y fotocromaticas propias de galaxias. 90 de los mismos han sido visualmente inspeccionados, y constituyen una muestra bona fide de galaxias en el VVV.
In [17]:
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from IPython.display import display
from astropy.io import ascii
from astropy.table import Table, Column
%matplotlib inline
Leemos los datos de la tabla de Baravalle et al. 2017
El catálogo se encuentra en formato cds, pero hubo que modificar los datos perdidos.
Simplemente se reemplazaron los valores -- por espacios vacios.
In [66]:
cat = ascii.read('Table1.all.txt', format='cds')
Se crea una columna donde podamos marcar los objetos identificados visualmente.
In [67]:
visual = []
for row in cat:
if row['Id'].endswith('*'):
visual.append(1)
else: visual.append(0)
In [68]:
vis = Column(data=visual, name='Visual', dtype='Bool')
cat.add_column(vis)
El catalogo entonces quedaría con esta pinta:
In [71]:
cat
Out[71]:
Guardamos el catalogo en formato csv para conservar la columna con identificacion visual
In [72]:
cat.write('cat_lau.csv', format='csv', overwrite=True)
Veamos las distribuciones y dependencias de los datos.
In [73]:
pdcat = cat.to_pandas()
print pdcat.columns
In [74]:
numeric_cols = pdcat[[u'Z', u'Y', u'J', u'H', u'Ks', u'Z$_{2\prime\prime}$', u'Y$_{2\prime\prime}$',
u'J$_{2\prime\prime}$', u'H$_{2\prime\prime}$', u'Ks$_{2\prime\prime}$',
u'R$_{1/2}$', u'C', u'$\epsilon$', u'n']]
pd.scatter_matrix(numeric_cols, alpha=0.3, figsize=(12, 12), diagonal='kde')
plt.show()