用C 语言画小猪佩奇 |这么火 你不了解下?

小猪佩奇是一只来自英国的粉红色小猪,脑袋长得有点像吹风筒,你有没有发现你的微信已经被小猪佩奇的人生格言轰炸了,朋友圈配图是它高傲开心的笑脸,好友的头像是它家族的成员,聊天记录到处是它的表情包~

#小猪佩奇#话题阅读量达到1.8亿,帖子数达1292个。

动画片在优酷上映后播放量累计达到132.6亿

媒体对《小猪佩奇》的报道大致可提取为四组数据:豆瓣评分9.2分:表示用户体验佳,好评率高;年播放量100亿次(国内):表示用户规模大,爆发期短;2016年收入70亿(全球):表示变现模式成熟,吸金能力佳;2020年预计收入120亿(全球):表示该IP持续生命力强,前景好。

近日最火的要数小猪佩奇品牌联名事件:

事件起源于优衣库发布小猪佩奇婴儿版联名睡衣,

各大网友们不淡定了,纷纷为自己喜爱的品牌logo设计网红款!只有品牌想不到的,没有网友设计不出来的品牌联名!

还有时尚界的大V深夜发媸徐老师策划的小猪佩奇和国产品牌结合的选题,把这只小猪P进了国内多个产品logo里,一下子吸引了众多目光,迅速收获了10W+。

还有网红恶搞的 Nike Air Force1 小猪佩奇定制款

小猪佩奇甚至成为了卡通界的带货一哥,淘宝上的小猪佩奇手表、棒棒糖更是销售过万,连包包也深受当红男明星的喜爱。

用 C 语言画小猪佩奇

既然小猪佩奇那么火,那么我们接下来就分享下如何用 C 语言画小猪佩奇

使用带符号距离场(signed distance field, SDF)表示圆形:

沿用这个方法表示形状,但这次我们想利用 ASCII 字符

|/=\

画出形状的外框,并填充内部,类似这样:

   =====
 //.....\\
||.......||
 \\.....//
   =====

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

我们通过差分求 SDF 的梯度近似值,然后用

atan2()

求出梯度的角度:

用 C 语言简单实现,在

画布中画一个半径 0.8 并带有 0.1 寛度外框的圆形:

#include #include #define T doubleT 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 y = -1; y < 1; y += 0.05, putchar('\n'))
        for (T x = -1; x < 1; x += 0.025)
            putchar(outline(x, y));}

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

// ASCII Peppa Pig by Milo Yip#include #include #include #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))]);}

2倍:

4倍:

8倍:

原文发布于微信公众号 - 奇点(qddata)

原文发表时间:2018-06-01

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏CDA数据分析师

大数据告诉你:如何在魔都捕获一只活的高富帅?

作者:团支书   数据及算法支持:城市数据团   友情数据支持:TalkingData   考上市里的公务员后,学姐的人生目标已完成大半。但一直以来...

25870
来自专栏生信宝典

2018 Cell系列相变最强综述,未来已来,你在哪?

Trends in Cell Biology (Cell系列综述, 2018 IF: 18.564)于2018年6月1日在线发表了Steven Boeynaem...

16310
来自专栏ml

HDUOJ----1166敌兵布阵(线段树单点更新)

敌兵布阵 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java...

29370
来自专栏IT大咖说

IT大咖撩段子

01 今天面试,面试官问我数据库索引原理,我说不清楚,面试官怼我“你做开发这么久,心里没点B数概念吗?”我听了立马走人了。回家好好查了查资料才发现,面试官可能说...

32990
来自专栏大数据文摘

手把手教你用数据分析帮女神学姐选婚房

237100
来自专栏腾讯社交用户体验设计

举世舞双ISUX新年派对 - 腾讯ISUX

16540
来自专栏阮一峰的网络日志

微软面试100题

网上流传一份微软公司面试的考题,共100道,真实性不可考。 这种题目用来面试,是否科学,那就智者见智了。在我看来,它们更像游戏。转贴如下,可供闲时娱乐。 ===...

41670
来自专栏大数据挖掘DT机器学习

抓取链家官网北京房产信息并用python进行数据挖掘

从2014年对楼市的普遍唱衰,到2015年的价格回暖,到底发生了怎样的改变?本文就尝试通过大数据来和丰富的图表,为大家展现数据背后的数据。 数据采集...

70090
来自专栏腾讯数据中心

敬请收藏|数据中心湿度控制方法比较

导语 温湿度是数据中心空气环境的重要指标,温度过高或过低都可能会影响到IT设备的正常运行,同样,湿度过低或过高也会对一下设备产生较为严重的影响,简单来讲,低湿会...

39430
来自专栏ml

HDUOJ-----1166敌兵布阵

敌兵布阵 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java...

27790

扫码关注云+社区

领取腾讯云代金券