十行以内，你写过哪些比较酷的 Matlab 代码？-白红宇

十行以内，你写过哪些比较酷的 Matlab 代码？

阅读量：2117 次

发布时间：2019-04-30

本文共 6933 字，大约阅读时间需要 23 分钟。

每行不要超过80字符

，You can be serious without a suit.

、

一键磨皮。。。

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)，例双边滤波，导向滤波等。不是什么黑魔法啊~~，就是，按步骤用Matlab化而已。

这里，给出最后的计算公式，

Dest =(Src * (100 - Opacity) + (Src + 2 * GuassBlur(EPFFilter(Src) - Src + 128) - 256) * Opacity) /100

EPFFilter是一种保边平滑的滤波器，

Opacity控制透明度，

GaussBlur是高斯滤波了，简单的PS的操作非常易于用

Matlab程序化。

磨皮后可以再增加美白功能，也是不超过十行代码的，效果看起来更佳（其实我是想说，可以卖一个好价钱呢）。若是对皮肤美白，磨皮，祛痘这些算法感兴趣，可以在进一步学习。

补充说明，

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行内完成的。

比如美图秀秀的滤镜，彩铅，素描等，我就不再泄露更多秘密了。

参考资料

•

，心中有许多愿望，能够实现有多棒。

、

不是自己写的。。。

贴个神犇代码。。

%% Spectral ResidualmyFFT = fft2(inImg); myLogAmplitude = log(abs(myFFT));myPhase = angle(myFFT);mySpectralResidual = myLogAmplitude - imfilter(myLogAmplitude, fspecial('average', 3), 'replicate'); saliencyMap = abs(ifft2(exp(mySpectralResidual + i*myPhase))).^2;

6行代码，1300+文章引用。。

----------------update--------------------

说6行代码最主要是印象中poster里面，第五行分成了两行写，写答案的时候也没重新数。

第五行比较长。超过了80个字符，算两行也不为过。

•

，少刷知乎多读书

我说个比较实用的：

每次写预习报告，实验报告都要在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

imwrite(A,'1.jpg')

效果如下：

直接贴上打印没问题！当时就觉得没有更酷的了！

•

，程序员

找茬游戏

•

，C，VB，Python，Matlab，Intel家的汇编

一个函数，程序跑完自动发个邮件给自己，7~8行，代码找到了再上传。

代码：

你以为我会告诉你我的邮箱密码？
下面是代码（已修改外函数版，共9行）

function mailme(subject,content,MailAddress,password,DataPath,MailServer)

setpref('Internet','E_mail',MailAddress);

setpref('Internet','SMTP_Server',MailServer);%SMTP服务器

setpref('Internet','SMTP_Username',MailAddress);

setpref('Internet','SMTP_Password',password);

props = java.lang.System.getProperties;

props.setProperty('mail.smtp.auth','true');

sendmail(MailAddress,subject,content,DataPath);

end

%主题 subject

%正文 content

%手机邮箱地址 MailAddress

%密码 password

%附件 DataPath

%邮件服务器 MailServer

•

char([73 32 119 97 110 116 32 97 32 103 105 114 108 32 102 114 105 101 110 100 46 46 46])

•

，以在下目前的能力只能回答体验类的问题OAQ

a=imread('F:\1.jpg');

b=imread('F:\2.jpg');

aa=0.5*a;

bb=0.5*b;

c=imadd(aa,bb);

imshow(c);

做重曝的代码不怎么酷

但是就是这段代码让我意识到

可能matlab才是最好的图像处理软件

TUT

ps.关于重曝可以参考摄影之友的这两篇东西↓

得到好看的重曝照片本身是技术活这段代码只是提供了一个简单粗暴的叠加方法（然后为我打开了新世界的大门）而已-v-

•

，Freedom，基拉的坐骑

for i=1:25

why

end

•

，计算机视觉大尺寸测量

=============================我又来更新了===============================

====================这次也不超过十行，上了颜色=============================

y=0.0001:0.0001:1;for i=1:10000x(i)=-1.2*y(i)*log(y(i))+1/30*exp(1)^(-(30*y(i)-30/exp(1))^4);endfill([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');endaxis equal;

•

，info & embedded system student

logo

看完用完觉得吊的别光感谢完事

请点赞，么么这个哒

•

，Artificial artificial intelligence

常年专注写短代码的人飘过。感谢的提及和关注。希望今后能写更多短代码 paper 给各位看官 :-)

刚刚看到另外一个有趣的问题，

顺手写完代码又优化了一下，一稿多投到这里：

dimNum = 8;hyperCube = reshape(0:10^dimNum-1, ones(1, dimNum)*10);allCoord = repmat({1:9}, [1,dimNum]);allCoord{end} = 2:9;hyperCube = hyperCube(allCoord{:});disp(sum(1./hyperCube(:)))

•

，machine learning

强烈鄙视这些追求代码短的，估计大多数人根本不懂Matlab内存分配的原理。

搞搞大矩阵的计算就知道了，把两行代码精简成一行会直接造成out of memory。

你们要是觉得 A=B*(C*D) 和 A=C*D;A=B*A;一样，你们就naive了。

•

，更短、更快、更好

当年机器学习课的一次作业，用Naive Bayes实现文档分类，我用4行Matlab代码搞定。

输入：

countTrain：N*V矩阵，第i行第j列表示第i篇训练文档中，单词j出现了几次。

labelTrain：N*1向量，第i个元素表示第i篇训练文档的类别。

countTest：M*V矩阵，第i行第j列表示第i篇测试文档中，单词j出现了几次。

代码：

Py = full(sparse(1, labelTrain, 1)) / length(labelTrain);    % 各类的先验概率Px_y = (sparse(labelTrain, 1:length(labelTrain), 1) * countTrain)' + alpha;                                    % 各类中各单词出现次数，alpha用作smoothingPx_y = bsxfun(@rdivide, Px_y, sum(Px_y));         % 归一化，恨不能跟上一行合并[~, labelTest] = max(bsxfun(@plus, countTest * log(Px_y), log(P_y)), [], 2);                                                         % 每篇测试文档的类别

•

，飞行器设计

A(A<0) = 0;

我也说个，比不上其他回答里功能那么高端，刚搬砖时候正好用到，可以把矩阵A中所有负数元素都替换为0。

•

，主照亮我

我来一个中文的

native2unicode([192 238 188 209 192 188 206 210 176 174 196 227])

看看就好，不要说破

•

，果粉，软狗

help xxx

•

，Py/推荐系统/爬虫/NLP/机器学习/审计

听说之前隔壁实验室写的裸眼2d戴眼镜3d的matlab代码一共7行，一行1w卖给sony了。。

•

，对所有知识状物体充满好奇，公众号:MScie…

我来发两个soeasy的。

蒙特卡洛法求pi

N=1000;%总粒子数

S=sqrt(rand(1,N).^2+rand(1,N).^2);%距离原点距离数列

NC=size(find(S<1),2);%求出数量

PPI=NC*4/N;%求出圆周率并且记录

画蝴蝶曲线

function [dy]=odefun(t,y)

dy=zeros(3,1);

dy(1)=-10*y(1)+10*y(2);

dy(2)=28*y(1)-y(2)-y(1)*y(3);

dy(3)=-8/3*y(3)+y(1)*y(2);

end

[t,y]=ode45(@odefun,[0 100],[0 0.1 0]);

plot3(y(:,1),y(:,2),y(:,3));