After several issues with pyfits, and tired of it being so overengineered, I've wrote my own fits I/O package in python, wrapping the C library cfitsio with ctypes.

Pretty easy, first version completely developed in 1 day.

https://github.com/zonca/pycfitsio