Make a wibbly image

dither.ts

// for any integer coordinates x, y, returns a float in [0,1] which is
// interpreted as threshold a float-valued function should have to
// reach in order to make a bitmap pixel.
function odither(x: number, y: number): number {
  return (x == 0 && y == 0
    ? 0
    : ([0, 2, 3, 1][(y % 2) * 2 + (x % 2)] + dither(x >> 1, y >> 1)) / 4);
}

function dither(x: number, y: number): number {
  return odither(x, y + Math.floor(x / Math.sqrt(3)));
}

function tableau(edge: number, f: (x: number, y: number) => number): number[][] {
  const table: number[][] = [...Array(edge).keys()].map(k => []);
  for (let ii = 0; ii < edge; ii++) {
    for (let jj = 0; jj < edge; jj++) {
      if (f(ii / edge, jj / edge) > dither(ii, jj))
        table[ii][jj] = 1;
    }
  }
  return table;
}

function poster(n: number, x: number): number {
  return Math.floor(x * n) / n;
}

function func(x: number, y: number): number {
  const dx = x + 0.04 * Math.sin(y * 20);
  const dy = y + 0.03 * Math.sin(x * 20);
  const ddx = dx + 0.04 * Math.sin(dy * 20);
  const ddy = dy + 0.03 * Math.sin(dx * 20);

  return poster(9, 7 * (Math.pow(ddx - 0.5, 2) + Math.pow(ddy - 0.5, 2))) + 0.2 * Math.sin(30 * y);
}

function bitmap(edge: number): string {
  const tab = tableau(edge, func);
  let str = `P1 ${tab.length} ${tab.length}\n`;
  for (let i = 0; i < tab.length; i++) {
    for (let j = 0; j < tab.length; j++) {
      str += (tab[i][j] ? 1 : 0) + ' ';
    }
    str += '\n';
  }
  return str;
}

import * as fs from 'fs';

fs.writeFileSync('/tmp/a.pbm', bitmap(800), 'utf8');