已经过去的 2020 是一个不怎么顺遂的一年,出入公共场所都需要体温监测,而人流量密集的商场,一般会采用热成像技术来快速测量体温。那么今天我们就来说说如何让一张普通图片变成具有热成像的效果。

最终效果如下:

最终效果

从本文中你可以了解到:

发光的基本原理,如何给图片施加不同颜色的光

图片模糊的原理,如何用高斯模糊解决重影问题

colorRamp是什么?texture2D的一些高级用法,如何使用colorRamp改变图片的色彩风格

如何结合以上技术实现热成像效果

本文不再过多赘述webglAPI的使用。只关注片段着色器,下面给出最基本的一个着色器:

precision mediump float;/*变量申明*/varying vec2 uv; // 从顶点着色器传入的UV坐标uniform sampler2D u_texture0; // 全局变量,纹理1(就是空洞骑士的原图)uniform sampler2D u_texture1; // 全局变量,纹理2(后续我们要使用的ColorRamp)uniform float u_r; // 全局变量,ColorRamp中红通道权重uniform float u_g; // 全局变量,ColorRamp中绿通道权重uniform float u_b; // 全局变量,ColorRamp中蓝通道权重uniform float u_blurIntensity; // 全局变量,模糊强度uniform float u_negativeAmount; // 全局变量,负片程度uniform float u_colorRampIntensity; // 全局变量,ColorRamp光度强度uniform float u_glowAmount; // 全局变量,发光强度uniform vec3 u_glowColor; // 全局变量,发光颜色void main () {vec4 color = texture2D (u_texture0, uv);gl_FragColor =color;}

在第一部分变量生命中中我们从web层传入了一系列变量,现在看可能有点云里雾里别着急,我们后面会逐个提到,在main函数中我们获取了纹理信息,并把图片呈现出来:

发光效果(Glow)

在GLSL中颜色用一个四维向量表示:

vec4 color = vec4(0,0,1.0,0.5,1.0);

四个分量分别表示RGBA是个通道

如果要取出比如rgb构成一个三维向量,可以简写成color.rgb

颜色值范围在[0,1]之间,对应web常用的颜色区间[0,255]

我们要介绍的第一个效果是发光,发光简单来说就是让图片变量。最原始的想法就是把颜色值变大:

void main () {vec4 color = texture2D (u_texture0, uv);color *= 2.0;gl_FragColor =color;}

和我们期待的一样果然变量了:

但能否不仅仅想让图片发白光,而是可以发任意颜色的光?很自然地我们想到让color乘以另一个color:

....vec4 color = texture2D (u_texture0, uv);color.rgb *= vec3(1.0,0.0,0.5); //这是洋红色....

图片的确如我们所想,画面变红了,但是却变暗了,这是由于在 vec3(1.0,0.0,0.5,); 的绿通道值为0,向量相乘后,原本图片中的绿通道信息全部丢失,变为0,(读者不妨验证一下变黑幅度大的区域,是否就是原来比较绿的区域)

为此我们需要使用加和而不是乘积的形式:

color.rgb += vec3(1.0,0.0,0.5)*color.rgb;

最后我们将颜色和强度替换成我们,传入的变量,发光的效果就实现了:

....vec4 color = texture2D (u_texture0, uv);vec4 emission = color;emission.rgb *= u_glowAmount * u_glowColor;color.rgb += emission.rgb;color.rgb *= color.a;gl_FragColor =color;....

模糊(Blur)

接下来我们来讲第二种效果模糊,核心思路也很简单,如果一个像素点的颜色信息掺杂了周围点的颜色信息,那么图片就会模糊。

为此我们选择讲像素点上下左右四个点的信息融合:

.....void main () {vec4 color = texture2D (u_texture0, uv);float step = 0.005; // 申明一个变量,控制像素偏移的步幅color += texture2D (u_texture0, uv+vec2(step,0.0)); // 正上方点color += texture2D (u_texture0, uv+vec2(-step,0.0)); // 正下方点color += texture2D (u_texture0, uv+vec2(0.0,step)); // 左侧点color += texture2D (u_texture0, uv+vec2(0.0,-step)); // 右侧点color /= 5.0; // 取平均.....

为了让效果更佳,我们引入权重体系,让原本的点在结果中贡献最大设权重为4,正方向的点权重为2,45度方向的点权重为1:

vec4 color = 4.0*texture2D (u_texture0, uv);float step = 0.005; // 申明一个变量,控制像素偏移的步幅color += 2.0*texture2D (u_texture0, uv+vec2(step,0.0)); // 正上方点color += 2.0*texture2D (u_texture0, uv+vec2(-step,0.0)); // 正下方点color += 2.0*texture2D (u_texture0, uv+vec2(0.0,step)); // 左侧点color += 2.0*texture2D (u_texture0, uv+vec2(0.0,-step)); // 右侧点color += 1.0*texture2D (u_texture0, uv+vec2(-step,-step)); // 第4象限点color += 1.0*texture2D (u_texture0, uv+vec2(step,-step)); // 第3象限点color += 1.0*texture2D (u_texture0, uv+vec2(-step,step)); // 第2象限点color += 1.0*texture2D (u_texture0, uv+vec2(step,step)); // 第1象限点color /= (4.0+4.0*2.0+4.0*1.0);        // 取平均

正当我们满心欢喜第把step改成我们传入的变量u_blurIntensity(控制模糊强度时)悲剧发生了:

这是因为我们虽然加上了权重,但点的影响值却没有随 step 的怎大而递减,举个例子:就是当步幅为 0.001 和 0.1 时正上方的点贡献都是 2(而实际上 0.001 要比 0.1 有更大的权重,因为它更接近原始点)。

所以当step增大时很明显就会出现重影。

为此我们参考高斯模糊,引入一个指数函数:

该函数随着传入的x的增大而减小,于是我们改造上述代码:

将取值点由3*3增加到了8*8;

采用递减的指数函数:exp(-(x * x)/(2.0* bhqp * bhqp));

控制了偏移的范围(float(iy - halfIterations)*0.00390625);

float Blur_Gauss (float bhqp, float x){return exp (-(x * x) / (2.0 * bhqp * bhqp));}vec4 Blur (vec2 uv, sampler2D source, float Intensity){const int iterations = 16; // 常量才可以进行for循环int halfIterations = iterations / 2;float sigmaX = 0.1 + Intensity * 0.5;float sigmaY = sigmaX;float total = 0.0;vec4 ret = vec4 (0., 0., 0., 0.);float step = 0.00390625;// 增多到8*8个点for (int iy = 0; iy < iterations; ++iy) {float fy = Blur_Gauss (sigmaY, float (iy - halfIterations));float offsety = float (iy - halfIterations) * step;for (int ix = 0; ix < iterations; ++ix) {float fx = Blur_Gauss (sigmaX, float (ix - halfIterations));float offsetx = float (ix - halfIterations) * step;total += fx * fy;vec4 a = texture2D (source, uv + vec2 (offsetx, offsety));a.rgb *= a.a;ret += a * fx * fy;}}return ret / total;}.....color = Blur (uv, u_texture0, u_blurIntensity);.....

负片(negative)

这个可能是本文最简单的一个效果,就是让图片呈现出负片的效果,代码如下:

color.rgb = abs(u_negativeAmount - color.rgb);

colorRamp

colorRamp 简单理解就是色卡,就是下面这种:

在进一步说明 colorRamp 前,我们先想想如何把图片变成黑白?

最直接的想法就是把R,G,B三个通道的颜色取平均:

......vec4 color = texture2D (u_texture0, uv);float p = color.r+color.g+color.b;color.rgb = vec3(p/3.0);......

其实这件事也可以用ColorRamp实现,我们想找一个由黑到白的色表,传入着色器作为u_texture1。

纹理最左侧的点x坐标就是0,最右侧的点x坐标就是1,因此黑白效果的代码可以写成:

......vec4 color = texture2D (u_texture0, uv);float p = color.r+color.g+color.b;color.rgb = texture2D (u_texture1, vec2 (p/3.0, .0)).rgb;......

这里浮点素P记录了原始图片(u_texture0)的信息,同时它制定了规则,输出的颜色是R,G,B三个通道值取平均的结果。

当然一个黑白效果不必如此大费周章,生产环境中ColorRamp更加缤纷,下面我们依旧使用(color.r+color.g+color.b)/3.0作为取色规则,来看看一些coloRamp的效果:

当然我们不会局限于一种规则,为此我们使用u_r,u_g,u_b分别控制红绿蓝三色通道权重:

......vec4 color = texture2D (u_texture0, uv);float p = dot(color.rgb,vec3(u_r,u_g,u_b));// 等同于color.r*u_r+color.g*u_g+color.b*u_bp += u_colorRampIntensity;// u_colorRampIntensity控制颜色统一向colorRamp的左/右侧移动color.rgb = texture2D (u_texture1, vec2 (p, .0)).rgb;......

我们一上图中第二种纹理为例,实现效果如下:

最总我们使用下面这个风格的colorRamp:

最后就是见证奇迹的是时刻了,我们把上面所过的效果合起来(发光=>colorRamp=>负片=>模糊):

最终效果

当然这个 shader 也可用在真人身上,比如封面图。其中 blurIntensity 越大风格越卡通,大家可以自己实践一下。

最近很常用的一张表情包

医院的热成像java怎么实现_10 行代码搞定“热成像”【shader风格化】相关推荐

  1. java微信支付代码_10行代码搞定微信支付(Java版)

    原标题:10行代码搞定微信支付(Java版) 微信支付痛点 对于大多数同学来说,要开发微信支付可不简单.附上微信支付官方文档网页链接 从文档上可以看出,你需要解决很多问题,我就随便挑几个吧. xml与 ...

  2. resnet50代码_13、SOTA论文实践-学习ResNet(80行代码搞定残差backbone网络)

    0.论文 Camera Distance-aware Top-down Approach for 3D Multi-person Pose Estimation from a Single RGB I ...

  3. python做事件研究法_35行代码搞定事件研究法(上)

    作者简介: 祝小宇,个人公众号:大猫的R语言课堂 这期大猫课堂将会教大家如何用35行R代码写出最有效率的事件研究法. 注意,本代码主要使用data.table完成,关于data.table包的相应知识 ...

  4. 35行代码搞定事件研究法(下)

    作者简介: 祝小宇,个人公众号:大猫的R语言课堂 前文推送: 35行代码搞定事件研究法(上) Hello亲爱的小伙伴们,上期已经讲到如何对单一事件日计算超额收益,本期将会教大家如何针对多个股票多个事件 ...

  5. 国外stripe支付,超简单几行代码搞定

    国外stripe支付,超简单几行代码搞定 海外的项目 需要stripe支付 很简单 几行代码 先加入依赖: compile 'com.stripe:stripe-android:6.1.2' 总共两种 ...

  6. 万万想不到 10行代码搞定一个决策树

    01决策树模拟实验 文章目录 01决策树模拟实验 要求 决策树简单介绍 搭建环境 产生数据集 划分训练集和测试集 生成决策树 Cross-Validation法 可视化决策树 10行代码搞定决策树 要 ...

  7. 35行代码搞定事件研究法(上)

    作者简介: 祝小宇,个人公众号:大猫的R语言课堂 这期大猫课堂将会教大家如何用35行R代码写出最有效率的事件研究法. 注意,本代码主要使用data.table完成,关于data.table包的相应知识 ...

  8. python数学公式编辑工具_1行代码搞定Latex公式编写,这个4.6M的Python小插件,堪称论文必备神器...

    原标题:1行代码搞定Latex公式编写,这个4.6M的Python小插件,堪称论文必备神器 来源:量子位 关注前沿科技 萧箫 发自 凹非寺 量子位 报道 | 公众号 QbitAI 萧箫 发自 凹非寺 ...

  9. python修改文件名称唯美_5行代码搞定14种常见图的python可视化库,还自带16种优美主题,真是太喜欢了...

    原创:小dull鸟 python数据分析之禅 原文链接: 5行代码搞定14种常见图的python可视化库,还自带16种优美主题,真是太喜欢了​mp.weixin.qq.com 有时候我们做数据可视化并 ...

最新文章

  1. 重磅!2K图像90FPS,中科院开源轻量级通用人脸检测器
  2. android自动登录不过登陆界面,【教程】Android 记住密码和自动登录界面的实现
  3. spring框架的引入
  4. 自适应网站设计对百度友好的关键:添加applicable-device标签(转)
  5. HTTP协议发展历程
  6. r语言pls分析_R语言:生存分析
  7. 【华为推荐】基于反事实学习的推荐系统研究.pdf(附pdf下载链接)
  8. 详解之-js闭包的用途
  9. 牛客寒假算法基础集训营2 A 处女座的签到题及一丢丢nth_elemen容器的简单用法...
  10. Kettle连接mysql错误 org.pentaho.ui.xul.XulException: java.lang.reflect.InvocationTargetException
  11. signature=735f4378ec01919f23285d0d2557be19,OPENSSL编程 第二十章 椭圆曲线
  12. 几张图片演示PD(PowerDesigner)生成数据库脚本sql文件
  13. 2017中国云计算技术大会将于5月18-19日在京召开
  14. wpf 三种 loading 样式
  15. H3C新华三链路聚合介绍
  16. 读HTTP权威指南的体会
  17. win10系统禁用笔记本自带键盘
  18. 嵌入式开发——程序跑飞原因总结
  19. 又一打脸现场!Fork Bunny 的 Merlin 损失 240 ETH
  20. JAVA实现输入一行字符,分别统计出其中英文字母、空格、数字和其它字符的个数

热门文章

  1. java 神剑_java语音短信的实现
  2. Felicia Brittman基于200篇文献对中式英语常见错误的人工总结 (附学习体会和写作视频)...
  3. 3dmax中如何制作汽车的行走动画
  4. k8s中因busybox镜像版本问题导致Pod中nslookup无法解析Service名称问题
  5. 【宇宙编码】厉害,生日那天:程序员小哥把天上的宇宙星云都送给女神了~
  6. AdobeAcrobat DC pro 2019版JH教程
  7. GQ1b(NH4+盐)说明书
  8. web页面开发,一个前端妹子的面试笔记,知乎上已获万赞
  9. html上下两个箭头符号怎么打出来,电脑箭头符号怎么打出来(往返箭头符号图案)...
  10. 谁说抓包必须用root权限