每行不要超过80字符
十行以内，你写过哪些比较酷的Mathematica代码？ - Wolfram Mathematica
一键磨皮。。。
I= double(imread('1.jpg'));
H = RF(I,30,100) - I + 128;
G = imfilter(H,fspecial('gaussian',[3 3],100));
opacity = 50;
Dst = (I*(100 - opacity)+(I+2*G-256)*opacity)/100;
imshow([uint8(I) uint8(Dst)]);
后面的图片可能有些不友善，慎重下拉！！！
这是输入的原图。
这是效果图
再来两张效果
这里， RF是递归域变换滤波器（ Domain Transform Filter，点此 下载），也可换成其它保边平滑的滤波器( Edge Perserving Filter)，例双边滤波，导向滤波等。不是什么黑魔法啊~~，就是 PhotoShop平湖法磨皮：完美保留皮肤纹理磨皮教程_PS图片处理，按步骤用Matlab化而已。
这里，给出最后的计算公式，
Dest =(Src * (100 - Opacity) + (Src + 2 * GuassBlur(EPFFilter(Src) - Src + 128) - 256) * Opacity) /100
EPFFilter是一种保边平滑的滤波器， Opacity控制透明度， GaussBlur是高斯滤波了，简单的PS的操作非常易于用 Matlab程序化。
磨皮后可以再增加美白功能，也是不超过十行代码的，效果看起来更佳（其实我是想说，可以卖一个好价钱呢）。若是对皮肤美白，磨皮，祛痘这些算法感兴趣，可以在 Imageshop博客 进一步学习。
补充说明，
Matlab并没有收录 RF这个函数，但 Matlab2014a已经有 imguidedfilter函数，在 Matlab2014a可以运行以下6行代码。
I= double(imread('7.png'));
H =double(imguidedfilter(uint8(I))) - I + 128;
G = imfilter(H,fspecial('gaussian',[3 3],100));
opacity = 50;
Dest = (I*(100 - opacity)+(I+2*G - 256)*opacity)/100;
imshow([uint8(I) uint8(Dest)]);
除了此外，图像处理中很多算法都可以用Matlab在10行内完成的。
比如美图秀秀的滤镜，彩铅，素描等，我就不再泄露更多秘密了。
参考资料
Imageshop 博客园
Guided filtering of images
Computational Photography
美图秀秀网页版 官方网站
我说个比较实用的：
每次写预习报告，实验报告都要在word里画什么流程图啊，连线图啊什么的，超烦啊！
我就想着能不能直接把实验指导书上的图拍下来直接放word里然后打印交上去，但是老师肯定会觉得我态度不认真不给我分啊！
这图打印出来绝壁不行啊：
于是我编写了如下代码（都看得懂）：
B = imread('1.png');
A= rgb2gray(B);
[m,n]=size(A);
for i=1:1:m
for j=1:1:n
if(A(i,j)>110)
A(i,j)=255;
end
end
end
imwrite(A,'1.jpg')
效果如下：
直接贴上打印没问题！当时就觉得没有更酷的了！
每次写预习报告，实验报告都要在word里画什么流程图啊，连线图啊什么的，超烦啊！
我就想着能不能直接把实验指导书上的图拍下来直接放word里然后打印交上去，但是老师肯定会觉得我态度不认真不给我分啊！
这图打印出来绝壁不行啊：
于是我编写了如下代码（都看得懂）：
B = imread('1.png');
A= rgb2gray(B);
[m,n]=size(A);
for i=1:1:m
for j=1:1:n
if(A(i,j)>110)
A(i,j)=255;
end
end
end
imwrite(A,'1.jpg')
效果如下：
直接贴上打印没问题！当时就觉得没有更酷的了！
=============================我又来更新了===============================
====================这次也不超过十行，上了颜色=============================
----------------------------------------更新一个-------------------------
这个没有什么技术含量，就是好玩。这一阵玩以撒的结合，看到激光加跟踪觉得蛮有意思。
于是用matlab做了一个。图如下，代码太多超过10行了就不贴了。
左边浅蓝色是激光的跟踪范围，右边是激光效果。红色星号代表怪物。
-------------------------------再次更新----------------------------------
这次是一个把图片转换称如下风格的程序，有一句我分成了两行，去掉最后依据修饰用的axis equal;刚好10行。有哪位同好知道是什么风格麻烦告诉我。
====================这次也不超过十行，上了颜色=============================
y=0.0001:0.0001:1;
for i=1:10000
x(i)=-1.2*y(i)*log(y(i))+1/30*exp(1)^(-(30*y(i)-30/exp(1))^4);
end
fill([x(1:5000),x(5001:10000)],[y(1:5000),y(5001:10000)],[253/255,226/255,202/255],'edgealpha',0);
axis ([0 1 0 1 ]);
----------------------------------------更新一个-------------------------
这个没有什么技术含量，就是好玩。这一阵玩以撒的结合，看到激光加跟踪觉得蛮有意思。
于是用matlab做了一个。图如下，代码太多超过10行了就不贴了。
左边浅蓝色是激光的跟踪范围，右边是激光效果。红色星号代表怪物。
-------------------------------再次更新----------------------------------
这次是一个把图片转换称如下风格的程序，有一句我分成了两行，去掉最后依据修饰用的axis equal;刚好10行。有哪位同好知道是什么风格麻烦告诉我。
I=imread('测试2.jpg');
[row,col,~]=size(I);
X=[1+round((col-1)*rand(1000,1)),1+round((row-1)*rand(1000,1))];
tri = delaunay(X);
ns=length(tri);%
I=double(I);
for i=1:ns
patch(X(tri(i,:),1),col*ones(3,1)-X(tri(i,:),2),[I(X(tri(i,2),2),X(tri(i,2),1),1)/255, ...
I(X(tri(i,2),2),X(tri(i,2),1),2)/255, I(X(tri(i,2),2),X(tri(i,2),1),3)/255],'edgecolor','none');
end
axis equal;
一条语句就可以做文字识别：
ocrobj = ocr(imread('businessCard.png'));
ocrobj.Text
‘ MathWorks®
The Mathworks, Inc.
3 Apple Hill Drive
Natick, MA 01760-2098
USA
www. mathworkscom
ocrobj = ocr(imread('businessCard.png'));
ocrobj.Text
‘ MathWorks®
The Mathworks, Inc.
3 Apple Hill Drive
Natick, MA 01760-2098
USA
www. mathworkscom
1万+
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
