Use fits.open instead of NDData.read (until the latter is implemented)