<SYSTEM_TASK:>
Given the following text description, write Python code to implement the functionality described below step by step
<END_TASK>
<USER_TASK:>
Description:
 Playing around with the dummy Gaia archive
gully
Sept 2, 2016
Attempt 1
Step1:  Apparently not in gzip format, despite the file extension
Step2:  Attempt 2
<ASSISTANT_TASK:>
Python Code:
! gunzip -d ../data/result.vot_01.gz
Explanation: Playing around with the dummy Gaia archive
gully
Sept 2, 2016
Attempt 1: 100000 sources in Taurus, $1^{\circ}$ radius
SQL
SELECT TOP 100000 * FROM public.gaia_source  WHERE CONTAINS(POINT('ICRS',public.gaia_source.ra,public.gaia_source.dec),CIRCLE('ICRS',64.1171208,28.1266139,1.0))=1
End of explanation
#! head -n 200 ../data/result.vot_01.gz
import pandas as pd
import numpy as np
import seaborn as sns
%config InlineBackend.figure_format = 'retina'
%matplotlib inline
import matplotlib.pyplot as plt
from astropy.io import votable
votable.is_votable('../data/result.vot_01.gz')
out = votable.parse_single_table('../data/result.vot_01.gz')
out2 = out.to_table()
out3 = out2.to_pandas()
out2[0:5]
plt.plot(out3.ra, out3.dec, '.', alpha=0.1)
out3.columns
out3.shape
vec = out3.parallax == out3.parallax
vec.sum(), len(vec)
sns.distplot(out3.parallax[vec])
sns.distplot(out3.phot_g_mean_mag, kde=False)#, bins=np.arange(0,2E7, 10), kde=False)
plt.yscale('log')
out3.head()
Explanation: Apparently not in gzip format, despite the file extension
End of explanation
! du -hs ../data/result.vot_02.gz
tab2_raw = votable.parse_single_table('../data/result.vot_02.gz')
tab2_apy = tab2_raw.to_table()
tab2_apy
tab2 = tab2_apy.to_pandas()
tab2
Explanation: Attempt 2: public.dual example.
End of explanation
<END_TASK>