[AstroPy] Convert FITS header into astropy Table entries?

[Erik Bray]

On 08/19/2015 11:26 AM, Matthew Craig wrote:
> Hi Jason,
>
> I’m late to the party, but ccdproc is about to add an ImageFileCollection obj
> that does (I think) what you are looking for…given a directory containing FITS
> files, make a table out of their headers (will not handle multiple extensions in
> a single file gracefully).

Couldn't it just include a column for the extension within that file?  Or is 
there something I'm missing?

Erik

> That functionality is available right now from another package, msumastro
> (http://msum-astro.readthedocs.org/en/latest/image_management.html or
> https://github.com/mwcraig/msumastro) so if you want it today I would get it
> from there.
>
> If you decide to give it a try and find issues please report them, either at
> https://github.com/astropy/ccdproc/issues or at
> https://github.com/mwcraig/msumastro/issues or by email.
>
> In a nutshell you would do:
>
>  >>> from msumastro import ImageFileCollection
>  >>> my_collection = ImageFileCollection(‘path/to/dir/with/files’, keywords=‘*’)
>>>>  my_table = my_collection.summary
>
>
> Matt Craig
>
> schedule: http://physics.mnstate.edu/craig
> ——
>
> Professor
> Department of Physics and Astronomy
> Minnesota State University Moorhead
> 1104 7th Ave S, Moorhead MN 56563
>
> phone: (218) 477-2439
> fax: (218) 477-2290
>
>> On Aug 19, 2015, at 9:00 AM, Jason Tumlinson <tumlinson at stsci.edu
>> <mailto:tumlinson at stsci.edu>> wrote:
>>
>> Hi,
>> I am trying to harvest the information from a large number of FITS headers,
>> parse them and do statistics on values, in an efficient manner. I read in each
>> header in the usual fashion:
>>
>> h = fits.open('lbiy83i0q_x1d.fits’)
>> h0 = h[0].header
>>
>> My goal then is to get all the header info into an astropy Table such that the
>> table columns are given by the header keys and the values in the rows are
>> given by the values in each of 100+ headers. I tried this:
>>
>> In [38]: t = Table(rows=h0.values(), names=h0.keys())
>> ---------------------------------------------------------------------------
>> TypeError                                 Traceback (most recent call last)
>> <ipython-input-38-8a0ece2925d5>in <module>()
>> ----> 1t = Table(rows=h0.values(), names=h0.keys())
>>
>> /Users/tumlinson/Ureka/python/lib/python2.7/site-packages/astropy/table/table.pycin
>> __init__(self, data, masked, names, dtype, meta, copy, rows)
>>     240                 data = rows
>>     241 else:
>> --> 242rec_data = recarray_fromrecords(rows)
>>     243                 data = [rec_data[name] for name in rec_data.dtype.names]
>>     244
>>
>> /Users/tumlinson/Ureka/python/lib/python2.7/site-packages/astropy/table/np_utils.pycin
>> recarray_fromrecords(rec_list)
>>     604# to column arrays in the recarray module; it could be there is a better
>>     605# way
>> --> 606nfields = len(rec_list[0])
>>     607     obj = np.array(rec_list, dtype=object)
>>     608     array_list = [np.array(obj[..., i].tolist()) for i in range(nfields)]
>>
>> TypeError: object of type 'bool' has no len()
>>
>> to no avail. Is this possible, easy, or wrongheaded? Advise appreciated.
>>
>> Thanks,
>> JT
>>
>>
>>
>> _______________________________________________
>> AstroPy mailing list
>> AstroPy at scipy.org <mailto:AstroPy at scipy.org>
>> http://mail.scipy.org/mailman/listinfo/astropy
>
>
>
> _______________________________________________
> AstroPy mailing list
> AstroPy at scipy.org
> http://mail.scipy.org/mailman/listinfo/astropy
>