You'd might be better off with Numarray or Scipy. Image and sparse matrix modules are with in the package, and there is direct functionality for PIL integration. Numarray is better documented: http://stsdas.stsci.edu/numarray/numarray-1.3.html/index.html