Бинарный PBM/PGM/PPM

apt-get install python-pythonmagick

lines = [ 0xf8, 0xf8, ... ]; # уже прочитанные бинарные строки
for line in lines:
    bits = [];
    for i in range(1, 5): # вместо 5 подставить сколько надо цифер
        bits.append(line & 1);
        line = line >> 1;
    # тут мы имеем
    # bits = [ 1, 1, 1, 1, 1 ];

import array

f = open('tux_logo.pbm', 'rb')
f.readline() # header
f.readline() # comment
dimensions = f.readline() # sizes
width, height = [int(x) for x in dimensions.strip().split()]
for y in range(height):
    bytewidth = (width - 1)//8 + 1
    bline = array.array('B')
    bline.fromfile(f, bytewidth)
    x = 0
    ptr = 0
    oline = []
    while x < width:
        for shift in range(7,-1,-1):
            bit = (bline[ptr] >> shift) & 1
            if bit == 1:
                oline.append('*')
            else:
                oline.append(' ')
            x += 1
            if x >= width: break;
        ptr += 1
    print ''.join(oline)

import math
from array import array

def bitspow(integer):
  '''bitspow(integer) -> int
  
  Return the next highest power of 2 for integer.'''
  if not isinstance(integer, int):
    raise(TypeError('bitspow() needs int, not %s' % \
      type(integer).__name__))
  if integer < 0:
    raise(ValueError('bitspow() accepts only positive integers'))
  integer |= integer >> 1;
  integer |= integer >> 2;
  integer |= integer >> 4;
  integer |= integer >> 8;
  integer |= integer >> 16;
  integer += 1;
  integer = math.log(integer, 2)
  integer = int(integer)
  return(integer)

def bits(integer, check=True):
  '''bits(integer) -> list
  
  Return the bits as list of integers. Bits are packed as byte, the necessary
  preceding zeros are added. The principle of the function:
    if integer < 2**8, it is transformed to 8 bits (one byte);
    if integer < 2**16, it is transformed to 16 bits (two bytes);
    if integer < 2**32, it is transformed to 32 bits (three bytes);
  and so on. This method can be used to read packed data.
  This function may work significantly faster, if check argument is False.
  That means that the function won't check type of the first argument.'''
  if not isinstance(integer, int):
    raise(TypeError('bits() needs int, not %s' % \
      type(integer).__name__))
  if integer < 0:
    raise(ValueError('bits() accepts only positive integers'))
  if integer <= 0xFF:
    size = 8
  elif integer <= 0xFFFF:
    size = 16
  elif integer <= 0xFFFFFF:
    size = 24
  elif integer <= 0xFFFFFFFF:
    size = 32
  else: # either case
    size = bitspow(integer)
  while size % 8:
    size += 1
  bits = bin(integer)[2:]
  bits = '0' *(size -len(bits)) +bits
  return(bits)

data = b'\x00\x01\x02\x03' *9000
data = [ord(byte) for byte in data]
data = b''.join([bits(byte) for byte in data])
data = array('B', data.replace(b'0', b'\x00').replace(b'1', b'\x01'))