如何优雅的使用C语言绘制一只小猪佩奇

2018-06-17 20:27:18来源:未知 阅读 ()

新老客户大回馈,云服务器低至5折

今天我们来用C语言画一只小猪佩奇---社会、社会....
在画小猪佩奇之前,我们先使用带符号的距离长 (signed distance field,SDF) 来画一个圆形。

使用这个方法表示形状,但是这次我们使用 ASCⅡ 字符 \/=\ 画出形状的外框,并填充内部,类似这样:

  =====

//.....\\

||.....||

\\....//    

  =====

SDF 的梯度(gradient)代表 SDF 变化最大的方向,可用这个方向去决定用哪一个字符。

我们通过差分求 SDF 的梯度近似值,然后用 atan2() 求出梯度的角度:

用 C 语言简单实现,在 [-1, 1] x [-1, 1] 画布中画一个半径 0.8 并带有 0.1 宽度外框的圆形: 
#include <math.h>
#include <stdio.h>
#define T double

T f(T x, T y)
{
  return sqrt(x x + y y) - 0.8f;
}

char outline(T x, T y)
{

   T delta = 0.001;
  if (fabs(f(x, y)) < 0.05)
  {
    T dx = f(x + delta, y) - f(x - delta, y);
    T dy = f(x, y + delta) - f(x, y - delta);
    return "|/=\|/=\|"[(int)((atan2(dy, dx) / 6.2831853072 + 0.5) * 8 + 0.5)];
  }
  else if (f(x, y) < 0)
  {
    return '.';
  }
  else
  {
  return ' ';

  }
}

int main()
{
  for (T i = -1; i < 1; i += 0.05, putchar('\n'))
  {
    for (T j = -1; j < 1; j += 0.025)
    {
      putchar(outline(j, i));
    } 

  }
  getchar();
  return 0;
}

然后,我们就可以通过画多个圆形,把它们适当地旋转和缩放,用构造实体几何比它们组合起来,从而就可以画出小猪佩奇了:


#include <math.h> #include <stdio.h> #include <stdlib.h> #define T double

T c(T x, T y, T r)
{
  return sqrt(x x + y y) - r;
}

T u(T x, T y, T t)
{
  return x cos(t) + y sin(t);
}

T v(T x, T y, T t)
{
  return y cos(t) - x sin(t);

}

T fa(T x, T y)

{

   return fmin(c(x, y, 0.5), c(x * 0.47 + 0.15, y + 0.25, 0.3));

}

T no(T x, T y)
{
  return c(x * 1.2 + 0.97, y + 0.25, 0.2);
}

T nh(T x, T y)

{   

   return fmin(c(x + 0.9, y + 0.25, 0.03), c(x + 0.75, y + 0.25, 0.03));

}

T ea(T x, T y)
{
  return fmin(c(x 1.7 + 0.3, y + 0.7, 0.15), c(u(x, y, 0.25) 1.7, v(x, y, 0.25) + 0.65, 0.15));
}

T ey(T x, T y)
{
  return fmin(c(x + 0.4, y + 0.35, 0.1), c(x + 0.15, y + 0.35, 0.1));
}

T pu(T x, T y)
{
  return fmin(c(x + 0.38, y + 0.33, 0.03), c(x + 0.13, y + 0.33, 0.03));
}

T fr(T x, T y)
{
  return c(x * 1.1 - 0.3, y + 0.1, 0.15);
}

T mo(T x, T y)
{
  return fmax(c(x + 0.15, y - 0.05, 0.2), -c(x + 0.15, y, 0.25));
}

T o(T x, T y, T(f)(T, T), T i)
{
  
T r =f(x, y);

  return fabs(r) < 0.02 ? (atan2(f(x, y + 1e-3) - r, f(x + 1e-3, y) - r) + 0.3) 1.273 + 6.5 : r < 0 ? i : 0;
}

T s(T x, T y, T(*f)(T, T), T i)
{
  return f(x, y) < 0 ? i : 0;
}

T f(T x, T y)
{
  return o(x, y, no, 1) ? fmax(o(x, y, no, 1), s(x, y, nh, 12)) : fmax(o(x, y, fa, 1), fmax(o(x, y, ey, 11), fmax(o(x, y, ea, 1), fmax(o(x, y, mo, 1), fmax(s(x, y, fr, 13), s(x, y, pu, 12))))));
}

int main(int a, char **b)
{
  for (T y = -1, s = a > 1 ? strtod(b[1], 0) : 1; y < 0.6; y += 0.05 / s, putchar('\n'))
  {
    for (T x = -1; x < 0.6; x += 0.025 / s)
    {
      putchar(" .|/=\|/=\| @!"[(int)f(u(x, y, 0.3), v(x, y, 0.3))]);
    }  
  }
  getchar();
  return 0;
}

两倍:

四倍:

参考:

知乎Milo Yip

标签:

版权申明:本站文章部分自网络,如有侵权,请联系:west999com@outlook.com
特别注意:本站所有转载文章言论不代表本站观点,本站所提供的摄影照片,插画,设计作品,如需使用,请与原作者联系,版权归原作者所有

上一篇:CodeChef March Lunchtime 2018 div2

下一篇:BZOJ4269: 再见Xor(线性基)