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