[AstroPy] Fits to data table

[Basmah Riaz]

Hi,

 I'm new to python and am trying to learn the basics of handling fits
images.

 Is it possible to extract a data table from a fits image cube? I'm trying
to write a data table of pixel (column 1) and integrated flux (column 2)
from the fits image.

I had followed the steps given here:

http://docs.astropy.org/en/stable/io/fits/index.html

In [1]: from astropy.io import fits

In [2]: hdul=fits.open('moment0.fits')

In [3]: hdul.info()
Filename: moment0.fits
No.    Name      Ver    Type      Cards   Dimensions   Format
  0  PRIMARY       1 PrimaryHDU      82   (98, 100, 1)   float32

In [4]: data = hdul[0].data

In [5]: data
Out[5]:
array([[[ nan,  nan,  nan, ...,  nan,  nan,  nan],
        [ nan,  nan,  nan, ...,  nan,  nan,  nan],
        [ nan,  nan,  nan, ...,  nan,  nan,  nan],
        ...,
        [ nan,  nan,  nan, ...,  nan,  nan,  nan],
        [ nan,  nan,  nan, ...,  nan,  nan,  nan],
        [ nan,  nan,  nan, ...,  nan,  nan,  nan]]], dtype=float32)

In [6]: data[0:2,1:2]
Out[6]:
array([[[             nan,              nan,              nan,
                      nan,              nan,              nan,
                      nan,              nan,              nan,
           1.04565542e+24,   8.37304559e+23,   9.44501255e+23,
           1.00440115e+24,   1.17914651e+24,   1.23599454e+24,
           1.30693366e+24,   1.25854511e+24,   1.46356597e+24,
           1.27364925e+24,   1.37587937e+24,   1.43363339e+24,
           1.29562667e+24,   1.33309057e+24,   1.28908398e+24,
           1.34466921e+24,   1.31942643e+24,   1.41336892e+24,
           1.36305788e+24,   1.33955442e+24,   1.34887939e+24,
           1.32828144e+24,   1.45468487e+24,   1.47349839e+24,
           1.40765259e+24,   1.39326082e+24,   1.42740358e+24,
           1.41948114e+24,   1.43646987e+24,   1.44934440e+24,
           1.44043736e+24,   1.46360921e+24,   1.46183486e+24,
           1.54993752e+24,   1.43122306e+24,   1.51683685e+24,
           1.41424529e+24,   1.45172604e+24,   1.47353255e+24,
           1.47837381e+24,   1.41663414e+24,   1.38960880e+24,
           1.43003786e+24,   1.57839825e+24,   1.40306195e+24,
           1.46425600e+24,   1.36871310e+24,   1.46037324e+24,
           1.46366829e+24,   1.43932046e+24,   1.35849389e+24,
           1.49760929e+24,   1.38901879e+24,   1.39916551e+24,
           1.48571042e+24,   1.47773552e+24,   1.44545631e+24,
           1.47139647e+24,   1.43888524e+24,   1.43438077e+24,
           1.50250258e+24,   1.36998895e+24,              nan,
                      nan,              nan,              nan,
                      nan,              nan,              nan,
                      nan,              nan,              nan,
                      nan,              nan,              nan,
                      nan,              nan,              nan,
                      nan,              nan,              nan,
                      nan,              nan,              nan,
                      nan,              nan,              nan,
                      nan,              nan]]], dtype=float32)

I can see the flux values at a range in pixel values, but how do I write it
as a simple two-column  ascii file?

Cheers,
Basmah
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20190319/d3e203d9/attachment.html>