code from Eero Simoncelli NYU Lab for Computational Vision