android 磨皮原理,Android平台Camera实时滤镜实现方法探讨(九)–磨皮算法探讨(一)

android 磨皮原理,Android平台Camera实时滤镜实现方法探讨(九)–磨皮算法探讨(一)上一篇开头提到了一些可用于磨皮的去噪算法,下面我们实现这些算法并且观察效果,咱不考虑实时性的问题该算法利用图像局部统计特性进行滤波处理,例如NXM像素的灰度图,首先计算点(i,j)所在窗口内(大小为(2n+1)(2m+1))的平均值m(i,j)以及均方差:得到加性去噪后的结果为:其中:1.根据原文提出的优化方法,首先是建立两个积分图,如图所示,点4的积分即为Sum(Ra)+Sum(Rb)+…

大家好,又见面了,我是你们的朋友全栈君。

上一篇开头提到了一些可用于磨皮的去噪算法,下面我们实现这些算法并且观察效果,咱不考虑实时性的问题

该算法利用图像局部统计特性进行滤波处理,例如NXM像素的灰度图,首先计算点(i,j)所在窗口内(大小为(2n+1)(2m+1))的平均值m(i,j)

39bc60063d9521bfe7f6eaf81e774902.png

以及均方差:

9c4fc4d73019cf4841b92403057efd24.png

得到加性去噪后的结果为:

ac4053e0f895315a54600deee0f05a34.png

其中:

05529e263efebf0f7d27853ae7844cf8.png

1.根据原文提出的优化方法,首先是建立两个积分图,如图所示,点4的积分即为Sum(Ra)+Sum(Rb)+Sum(Rc)+Sum(Rd)。积分图的建立算法可以参考这篇文章进行简单优化,然后即可根据积分图计算公式中的m值和v值。

25a13b5b780d82647f46de653567ff6b.png

例如半径为r的窗口的m(i,j)为Integral(i+r,j+r) + Integral(i-r-1,j-r-1)-Integral(i+r,j-r-1)-Integral(i-r-1,j+r)。代码如下,分别求1次方和平方的积分图。

void MagicBeauty::initIntegral(uint8_t* inputMatrix){

LOGE(“initIntegral start”);

if(mIntegralMatrix == NULL)

mIntegralMatrix = new uint64_t[mImageWidth * mImageHeight];

if(mIntegralMatrixSqr == NULL)

mIntegralMatrixSqr = new uint64_t[mImageWidth * mImageHeight];

uint64_t *columnSum = new uint64_t[mImageWidth];

uint64_t *columnSumSqr = new uint64_t[mImageWidth];

columnSum[0] = inputMatrix[0];

columnSumSqr[0] = inputMatrix[0] * inputMatrix[0];

mIntegralMatrix[0] = columnSum[0];

mIntegralMatrixSqr[0] = columnSumSqr[0];

for(int i = 1;i < mImageWidth;i++){

columnSum[i] = inputMatrix[i];

columnSumSqr[i] = inputMatrix[i] * inputMatrix[i];

mIntegralMatrix[i] = columnSum[i];

mIntegralMatrix[i] += mIntegralMatrix[i-1];

mIntegralMatrixSqr[i] = columnSumSqr[i];

mIntegralMatrixSqr[i] += mIntegralMatrixSqr[i-1];

}

for (int i = 1;i < mImageHeight; i++){

int offset = i * mImageWidth;

columnSum[0] += inputMatrix[offset];

columnSumSqr[0] += inputMatrix[offset] * inputMatrix[offset];

mIntegralMatrix[offset] = columnSum[0];

mIntegralMatrixSqr[offset] = columnSumSqr[0];

// other columns

for(int j = 1; j < mImageWidth; j++){

columnSum[j] += inputMatrix[offset+j];

columnSumSqr[j] += inputMatrix[offset+j] * inputMatrix[offset+j];

mIntegralMatrix[offset+j] = mIntegralMatrix[offset+j-1] + columnSum[j];

mIntegralMatrixSqr[offset+j] = mIntegralMatrixSqr[offset+j-1] + columnSumSqr[j];

}

}

delete[] columnSum;

delete[] columnSumSqr;

LOGE(“initIntegral end”);

}

2.根据网上抄来的RGB肤色检测计算肤色区域

void MagicBeauty::initSkinMatrix(){

LOGE(“start – initSkinMatrix”);

if(mSkinMatrix == NULL)

mSkinMatrix = new uint8_t[mImageWidth * mImageHeight];

for(int i = 0; i < mImageHeight; i++){

for(int j = 0; j < mImageWidth; j++){

int offset = i*mImageWidth+j;

ARGB RGB;

BitmapOperation::convertIntToArgb(mImageData_rgb[offset],&RGB);

if ((RGB.blue>95 && RGB.green>40 && RGB.red>20 &&

RGB.blue-RGB.red>15 && RGB.blue-RGB.green>15)||//uniform illumination

(RGB.blue>200 && RGB.green>210 && RGB.red>170 &&

abs(RGB.blue-RGB.red)<=15 && RGB.blue>RGB.red&& RGB.green>RGB.red))//lateral illumination

mSkinMatrix[offset] = 255;

else

mSkinMatrix[offset] = 0;

}

}

LOGE(“end – initSkinMatrix”);

}

3.根据公式对RGB通道或者将RGB通道转化为YCbCr格式单独对Y通道进行滤波

void MagicBeauty::startLocalStatisticsSmooth(float sigema){

if(mIntegralMatrix == NULL || mIntegralMatrixSqr == NULL ||

mImageData_yuv_y == NULL || mSkinMatrix == NULL || mImageData_yuv == NULL){

LOGE(“not init correctly”);

return;

}

int radius = mImageWidth > mImageHeight ? mImageWidth * 0.02 : mImageHeight * 0.02;

LOGE(“startSmooth”);

for(int i = 1; i < mImageHeight; i++){

for(int j = 1; j < mImageWidth; j++){

int offset = i * mImageWidth + j;

if(mSkinMatrix[offset] == 255){

int iMax = i + radius >= mImageHeight-1 ? mImageHeight-1 : i + radius;

int jMax = j + radius >= mImageWidth-1 ? mImageWidth-1 :j + radius;

int iMin = i – radius <= 1 ? 1 : i – radius;

int jMin = j – radius <= 1 ? 1 : j – radius;

int squar = (iMax – iMin + 1)*(jMax – jMin + 1);

int i4 = iMax*mImageWidth+jMax;

int i3 = (iMin-1)*mImageWidth+(jMin-1);

int i2 = iMax*mImageWidth+(jMin-1);

int i1 = (iMin-1)*mImageWidth+jMax;

float m = (mIntegralMatrix[i4]

+ mIntegralMatrix[i3]

– mIntegralMatrix[i2]

– mIntegralMatrix[i1]) / squar;

float v = (mIntegralMatrixSqr[i4]

+ mIntegralMatrixSqr[i3]

– mIntegralMatrixSqr[i2]

– mIntegralMatrixSqr[i1]) / squar – m*m;

float k = v / (v + sigema);

mImageData_yuv[offset*3] = m – k * m + k * mImageData_yuv_y[offset];

}

}

}

endLocalStatisticsSmooth();

}

效果图:磨皮强度为最大10.

bb221d9f38e2e83ad1ed04a00ac48049.png

a9ec56d707316ca3fcf3af9da328e1ad.png

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/162997.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • goland2020 激活码_通用破解码

    goland2020 激活码_通用破解码,https://javaforall.cn/100143.html。详细ieda激活码不妨到全栈程序员必看教程网一起来了解一下吧!

  • SLAM技术课程总结「建议收藏」

    SLAM技术课程总结「建议收藏」这学期SLAM技术课程遇到最多的问题,或者说困难更合适,其实和学生无关。专业培养方案中,SLAM技术是双语课程,但此课对数学要求极高,通常应用型本科院校学生的数学基础并不具备直接学习此课程的条件; 操作环境要求学生需在Linux系统下完成,虽然现在有云端实践平台,但掌握基本的Linux命令依然是必不可少的基础,对于机器人工程专业学生而言,课程体系忽略了这块训练,使此课孤立于课程体系中,无法和现有培养方案的其他课程形成良性互动; 编程要求C++11,当然能够掌握14以上更好,而之前的课程多侧重于C和

  • MySQL存储过程(PROCEDURE)(一)「建议收藏」

    MySQL存储过程(PROCEDURE)(一)「建议收藏」一、定义与目的:定义:存储过程是数据库SQL语言层面的代码封装与重用(是数据库中存储复杂程序,以便外部程序调用的一种数据库对象);目的:我们为了完成特定功能的SQL语句集,经编译创建并保存在数

  • Linux中断 – 综述

    Linux中断 – 综述

  • 如何新建vue项目(如何配置vue的开发环境)

    一、vue开发环境搭建1、安装nodejs步骤:在node官网(https://nodejs.org/en/download/)选择跟自己的电脑匹配的版本进行下载,然后一步步的安装即可,在cmd控制台输入node-v,如果出现版本信息即表示安装成功。2、npm包管理器是集成在node中的,所以直接输入npm-v就能查看到版本信息,若出现版本信息则表示npm能正常使用。3、输入npminstall-gcnpm–registry=http://registry.npm.taobao.o

  • pyinstaller 多个.py打包exe_python怎么生成py文件

    pyinstaller 多个.py打包exe_python怎么生成py文件一、python安装pyinstaller方法使用python编写脚本,需要发给别人使用的时候,总会想到如何打包成exe文件,发给对方。这样的话,对方可以直接使用运行,无需安装python。所以看网上的教程,大多使用pyinstaller。以下介绍下安装方法:1、在cmd控制台下,先升级pip版本,先执行命:pipinstall-Upip,若执行失败,控制台会提示新密令,按照提示…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号