canvas3.scala

//> using scala 3.3.0

/*
canvas3.scala

Initial image/canvas app with sierpinski triangles
https://en.wikipedia.org/wiki/Sierpi%C5%84ski_triangle

Creates a simple PPM

scala-cli canvas3.scala

 */

case class Colour(r: Int, g: Int, b: Int)

case class Loc(x: Int, y: Int):

  def mid(l: Loc) = Loc((x + l.x) / 2, (y + l.y) / 2)

case class Image[T](w: Int, h: Int, data: Vector[T]):

  def apply(l: Loc): T = data(l.x * h + l.y)

  def map[S](f: T => S): Image[S] = Image(w, h, data map f)

  def updated(l: Loc, value: T): Image[T] =
    if ((l.x >= 0) & (l.y >= 0) & (l.x < w) & (l.y < h))
      Image(w, h, data.updated(l.x * h + l.y, value))
    else this

  def line(l0: Loc, l1: Loc, c: T): Image[T] =
    val xd = math.abs(l1.x - l0.x)
    val yd = math.abs(l1.y - l0.y)
    val n = math.max(xd, yd) + 1
    val is = (0 to n).toList
    val x = is.map(i =>
      math.round(l0.x.toDouble * (n - i) / n + l1.x.toDouble * i / n).toInt
    )
    val y = is.map(i =>
      math.round(l0.y.toDouble * (n - i) / n + l1.y.toDouble * i / n).toInt
    )
    val ls = (x zip y) map (p => Loc(p._1, p._2))
    ls.foldLeft(this)((im, li) => im.updated(li, c))

  def tri(l0: Loc, l1: Loc, l2: Loc, c: T): Image[T] =
    val sorted = List(l0, l1, l2).sortWith(_.y < _.y)
    val iTop = 0 until (sorted(1).y - sorted(0).y)
    val linesTop = iTop.map(i =>
      (
        sorted(0).y + i,
        sorted(2).x.toDouble * i / (sorted(2).y - sorted(0).y) +
          sorted(0).x.toDouble * (sorted(2).y - sorted(0).y - i) / (sorted(
            2
          ).y - sorted(0).y),
        sorted(1).x.toDouble * i / (sorted(1).y - sorted(0).y) +
          sorted(0).x.toDouble * (sorted(1).y - sorted(0).y - i) / (sorted(
            1
          ).y - sorted(0).y)
      )
    )
    val iBot = 0 until (sorted(2).y - sorted(1).y)
    val linesBot = iBot.map(i =>
      (
        sorted(1).y + i,
        sorted(2).x.toDouble * (sorted(1).y - sorted(0).y + i) / (sorted(
          2
        ).y - sorted(0).y) +
          sorted(0).x.toDouble * (sorted(2).y - sorted(1).y - i) / (sorted(
            2
          ).y - sorted(0).y),
        sorted(2).x.toDouble * i / (sorted(2).y - sorted(1).y) +
          sorted(1).x.toDouble * (sorted(2).y - sorted(1).y - i) / (sorted(
            2
          ).y - sorted(1).y)
      )
    )
    val lines = linesTop ++ linesBot
    lines.foldLeft(this)((im, yxX) =>
      im.line(Loc(yxX._2.toInt, yxX._1), Loc(yxX._3.toInt, yxX._1), c)
    )

case object Image:

  def blank[T](w: Int, h: Int, c: T): Image[T] =
    Image(w, h, Vector.fill[T](w * h)(c))

  def blank(w: Int, h: Int): Image[Colour] =
    blank(w, h, Colour(255, 255, 255))

  def saveAsPPM(im: Image[Colour], fileName: String): Unit =
    val fs = new java.io.FileWriter(fileName)
    fs.write(s"P3 ${im.w} ${im.h} 255\n")
    (0 until im.h).foreach(y =>
      (0 until im.w).foreach(x =>
        val p = im(Loc(x, y))
        fs.write(s"${p.r} ${p.g} ${p.b}\n")
      )
    )
    fs.close()

object CanvasApp:

  val white = Colour(255, 255, 255)
  val black = Colour(0, 0, 0)
  val red = Colour(255, 0, 0)
  val green = Colour(0, 255, 0)
  val blue = Colour(0, 0, 255)

  extension [T](im: Image[T])
    def sier(lev: Int, l0: Loc, l1: Loc, l2: Loc, c: T): Image[T] =
      if (lev <= 0)
        im.tri(l0, l1, l2, c)
      else
        im.sier(lev - 1, l0, l0 mid l1, l0 mid l2, c)
          .sier(lev - 1, l2, l0 mid l2, l1 mid l2, c)
          .sier(lev - 1, l1, l0 mid l1, l1 mid l2, c)

  @main def main() =
    println("Hello")
    val im1 = (0 until 5).foldLeft(Image.blank(1000, 200))((im, i) =>
      im.sier(
        i,
        Loc(200 * i, 180),
        Loc(200 * i + 180, 180),
        Loc(200 * i + 90, 20),
        red
      )
    )
    Image.saveAsPPM(im1, "test3.ppm")
    println("Goodbye")

// eof