Android平台Camera实时滤镜实现方法探讨(九)--磨皮算法探讨(一)
上一篇开头提到了一些可用于磨皮的去噪算法。以下我们实现这些算法而且观察效果,咱不考虑实时性的问题
本文首先探讨的首先是《基于局部均方差相关信息的图像去噪及其在实时磨皮美容算法中的应用》
该算法利用图像局部统计特性进行滤波处理,比如NXM像素的灰度图,首先计算点(i,j)所在窗体内(大小为(2n+1)(2m+1))的平均值m(i,j)
以及均方差:
得到加性去噪后的结果为:
当中:
1.依据原文提出的优化方法,首先是建立两个积分图,如图所看到的。点4的积分即为Sum(Ra)+Sum(Rb)+Sum(Rc)+Sum(Rd)。积分图的建立算法能够參考这篇文章进行简单优化。然后就可以依据积分图计算公式中的m值和v值。
比如半径为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];</span>
}
}
}
endLocalStatisticsSmooth();
}
效果图:磨皮强度为最大10.
watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQv/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/SouthEast" alt="">
Android平台Camera实时滤镜实现方法探讨(九)--磨皮算法探讨(一)的更多相关文章
- Android平台Camera实时滤镜实现方法探讨(三)--通过Shader实现YUV转换RBG
http://blog.csdn.net/oshunz/article/details/50055057 文章例如该链接通过将YUV分成三个纹理,在shader中取出并且经过公式变换,转换成RGB.我 ...
- Android平台Camera实时滤镜实现方法探讨(十)--代码地址以及简单介绍(20160118更新)
简单做了个相机和图片编辑模块,时间原因非常多功能还没有做.尚有BUG,见谅,将在以后抽时间改动 代码地址 PS:请点个Star^-^ --------------------------------- ...
- Android平台Camera实时滤镜实现方法探讨(十一)--实时美颜滤镜
上一章完毕了对图片的磨皮处理.经过简单算法流程优化,能够达到非常快的速度.可是不能用于实时美颜.经实验,若採用仅仅处理Y信号的方案.半径极限大约是5-10,超过10则明显感受到卡顿.但对于1920X1 ...
- Unity3D之移植学习笔记:移植到Android平台
首先,我们需要一台已经配置好可以开发Android应用的计算机,这里我使用的是Windows系统+Eclipse+ADT的开发环境,当然也可以使用Android Studio或者使用Mac系统都可以. ...
- 【Android】Android Camera实时数据采集及通过MediaCodec硬编码编码数据的流程
吐槽: 其实常用流程都差不多,但是有时候还是会忘记某一步的详细用法,但是各位朋友请注意,官方已经不推荐Camera类的使用(现在是android.hardware.camera2),但无奈公司项目之前 ...
- Unity3D之移植学习笔记:移植到Android平台更好的方法
接上文,之前我们采用了直接将Unity项目导出为Eclipse项目来修改的方式,这种做法存在的一个最大的问题就是:每当Unity被修改之后,都需要重新导出,而每次导出的项目在Android平台方面的J ...
- 随笔之Android平台上的进程调度探讨
http://blog.csdn.net/innost/article/details/6940136 随笔之Android平台上的进程调度探讨 一由来 最近在翻阅MediaProvider的时候,突 ...
- 虹软人脸识别 - Android Camera实时人脸追踪画框适配
在使用虹软人脸识别Android SDK的过程中 ,预览时一般都需要绘制人脸框,但是和PC平台相机应用不同,在Android平台相机进行应用开发还需要考虑前后置相机切换.设备横竖屏切换等情况,因此在人 ...
- 高通Android平台硬件调试之Camera篇
之前一段时间有幸在高通android平台上调试2款camera sensor,一款是OV的5M YUV sensor,支持jpeg out,同时也支持AF,调试比较比较简单,因为别的项目已经在使用了, ...
随机推荐
- POJ2891 Strange Way to Express Integers (扩展欧几里德)
本文为博主原创文章,欢迎转载,请注明出处 www.cnblogs.com/yangyaojia 题目大意 求解一组同余方程 x ≡ r1 (mod a1) x ≡ r2 (mod a2) x ≡ r ...
- [Google Guava] 2.2-新集合类型
转自:并发编程网 原文链接:http://ifeve.com/google-guava-newcollectiontypes/ 链接博客其他文章中还有更多的guava其他功能的描述,有空可慢慢看. G ...
- windows server 2008开机自动登陆无密码,关机不必写原因
运行secpol.sec接下来,在弹出的“本地安全策略”对话框中,依次展开左边树图到“本地策略”-“安全选项”,在右边可以找到“交互式登录 无须按 Ctrl+Alt+Del”,双击该项设置为“已启用” ...
- 开源ETL工具kettle--数据迁移
背景 因为项目的需求,须要将数据从Oracle迁移到MSSQL,不是简单的数据复制,而是表结构和字段名都不一样.甚至须要处理编码规范不一致的情况,例如以下图所看到的 watermark/2/text/ ...
- iOS学习必须了解的七大手势
文章只要你有一点点基础应该就可以看的懂,文章只为学习交流 #import "ViewController.h" @interface ViewController () @prop ...
- CentOS7设置中文输入法
转自:https://i.cnblogs.com/EditPosts.aspx?postid=8327755&update=1 CentOS7设置中文输入法 安装CentOS7之后,鼓捣了半天 ...
- SPOJ 694/705 后缀数组
思路: 论文题*n Σn-i-ht[i]+1 就是结果 O(n)搞定~ //By SiriusRen #include <cstdio> #include <cstring> ...
- requireJS实现原理分析
下面requireJS实现的基本思路 项目地址https://github.com/WangMaoling/require var require = (function(){ //框架版本基本信息 ...
- HDU 2852 KiKi's K-Number【 树状数组 二分 】
题意:给出m个操作,0:是增加一个数,add(x,1)1:是删除一个指定的数,这个是看sum(x) - sum(x-1)是否为0,为0的话则不存在,不为0的话,则add(x,-1)2:是查询比x大的数 ...
- 转 c#中 base64字符串与普通字符串互转
https://blog.csdn.net/hwt0101/article/details/79758912 转成 Base64 形式的 System.String: string a = &q ...