教你使用MATLAB制作水波倒影特效

于 2021-08-02 16:10:04 发布
阅读量6.1k 收藏 159
点赞数 60
分类专栏： # MATLAB 特效 文章标签： 特效 数字图像处理 MATLAB 浮雕 水波倒影
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/slandarer/article/details/119325437
版权

注： 本文算法参考大佬 grafx 的这篇博客：
图像处理算法之水面倒影特效

由于本文使用MATLAB复现，因此很多语法上会显得比较简洁，同时本博文对原大佬文章部分内容进行了改写，详见本文：

0效果展示

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

1图像翻转及白化

导入图像： 这部分其实没啥好说的：

% 图片导入 
oriPic=imread('test.jpg');
[Row,Col,~]=size(oriPic);

在这里插入图片描述

翻转及白化图像：

翻转就是单纯的将行索引倒过来；
白化就是将当前像素的颜色按比例和白色取个带权均值，行索引越大白色权重也越大，图像也就越白。

% 图片翻转及白化 ==========================================================
whiteMat=((1:Row)./Row./1.2)'*ones(1,Col); % 白化比例矩阵
flipPic=zeros(Row,Col,3);                  % 翻转后矩阵初始化
for i=1:3
    tempChannel=double(oriPic(:,:,i));     % 获得通道图
    tempChannel=tempChannel(end:-1:1,:);   % 翻转
    tempChannel=tempChannel.*(1-whiteMat)+255.*whiteMat; % 白化
    flipPic(:,:,i)=tempChannel;
    
end

在这里插入图片描述
当然如果我们将这一行：

tempChannel=tempChannel.*(1-whiteMat)+255.*whiteMat;

改写为：

tempChannel=tempChannel.*(1-whiteMat)+0.*whiteMat;

就变成了一个黑化的过程：

在这里插入图片描述
当然你也可以尝试其他颜色，例如将整段改写为：

Color=[255,0,0];
colorMat=((1:Row)./Row./1.2)'*ones(1,Col); % 比例矩阵
flipPic=zeros(Row,Col,3);                  % 翻转后矩阵初始化
for i=1:3
    tempChannel=double(oriPic(:,:,i));     % 获得通道图
    tempChannel=tempChannel(end:-1:1,:);   % 翻转
    tempChannel=tempChannel.*(1-colorMat)+Color(i).*colorMat; % 渐变
    flipPic(:,:,i)=tempChannel;
    
end
imshow(uint8(flipPic))

在这里插入图片描述

2波纹图像构造

生成噪声并模糊

noiseMat=ones(Row,Col);
noiseMat=imnoise(noiseMat,'gaussian',0,5); % 噪声添加
gaussOpt=fspecial('gaussian',[3 3],1);
noiseMat=imfilter(noiseMat,gaussOpt);

噪声图：

在这里插入图片描述

模糊后噪声图：

在这里插入图片描述
浮雕特效

实际上浮雕特效就是用以下类似形式的矩阵对图像进行卷积，卷积结果在加上RGB范围的均值，[0,1]区间就加0.5,[0,255]区间就加128：
[ 1 0 0 0 0 0 0 0 − 1 ] [ 2 0 2 0 0 0 − 2 0 − 2 ] \begin{bmatrix} 1 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & -1 \end{bmatrix} \begin{bmatrix} 2 & 0 & 2 \\ 0 & 0 & 0 \\ -2 & 0 & -2 \end{bmatrix} ⎣⎡​100​000​00−1​⎦⎤​⎣⎡​20−2​000​20−2​⎦⎤​
数值和位置不重要，重要的是相对位置互为相反数，浮雕过程描述如下：

H=[cos(pi+pi/4),0,cos(pi-pi/4);
   cos(pi+2*pi/4),0,cos(pi-2*pi/4);
   cos(pi+3*pi/4),0,cos(pi-3*pi/4)];
noiseMat=imfilter(noiseMat,H,'conv')+0.5;  
noiseMat=noiseMat.*255;
noiseMat(noiseMat<0)=0;

在这里插入图片描述
透视变换
就是近大远小，这里为了方便起见只在横向方向上做了近大远小的拉伸，竖直方向进行了等比例拉伸，因而不是严格意义上的透视变换：
在这里插入图片描述
如图所示实际操作就是把左侧蓝色区域拉伸成右侧蓝色区域，并只选取红框内部分，代码如下：

% 图像透视变换 ============================================================
exNoiseMat=zeros(Row,Col);
% 横向拉伸上下边倍数
K1=10;K2=4;
for i=1:Row
    for j=1:Col
        k=K2+i*(K1-K2)/Row;
        nJ=(j-(1+Col)/2)/k+(1+Col)/2;
        if floor(nJ)==ceil(nJ)
            nJ=round(nJ);
            exNoiseMat(i,j)=noiseMat(i,nJ);
        else
            nJ1=floor(nJ);nJ2=ceil(nJ);
            exNoiseMat(i,j)=noiseMat(i,nJ1)*(nJ2-nJ)+noiseMat(i,nJ2)*(nJ-nJ1);
        end
    end
end
% 竖向拉伸3倍并只取一部分
exNoiseMat=imresize(exNoiseMat,[3*Row,Col]);
exNoiseMat=exNoiseMat(end-Row+1:end,:);
exNoiseMat=uint8(exNoiseMat);

在这里插入图片描述
注： 如果原图像尺寸过大，水波就会过于密集，这时候可以适当调整放缩倍数或者将原图像重调大小到小一点的尺寸。

例如大波浪代码：

% 图像透视变换 ============================================================
exNoiseMat=zeros(Row,Col);
K1=40;K2=10;
for i=1:Row
    for j=1:Col
        k=K2+i*(K1-K2)/Row;
        nJ=(j-(1+Col)/2)/k+(1+Col)/2;
        if floor(nJ)==ceil(nJ)
            nJ=round(nJ);
            exNoiseMat(i,j)=noiseMat(i,nJ);
        else
            nJ1=floor(nJ);nJ2=ceil(nJ);
            exNoiseMat(i,j)=noiseMat(i,nJ1)*(nJ2-nJ)+noiseMat(i,nJ2)*(nJ-nJ1);
        end
    end
end
exNoiseMat=imresize(exNoiseMat,[8*Row,Col]);
exNoiseMat=exNoiseMat(end-Row+1:end,:);
exNoiseMat=uint8(exNoiseMat);

小波浪和大波浪：

在这里插入图片描述
在这里插入图片描述

3扭曲置换

这个。。。老朋友了，具体原理还是看这一篇叭：MATLAB 制作抖音同款含褶皱面料图

% 扭曲置换 ================================================================
forePic=flipPic;
bkgPic=exNoiseMat;

exforePic=uint8(zeros(size(forePic)+[26,26,0]));
exforePic(14:end-13,14:end-13,1)=forePic(:,:,1);
exforePic(14:end-13,14:end-13,2)=forePic(:,:,2);
exforePic(14:end-13,14:end-13,3)=forePic(:,:,3);

for i=1:13
    exforePic(i,14:end-13,:)=forePic(1,:,:);
    exforePic(end+1-i,14:end-13,:)=forePic(end,:,:);
    exforePic(14:end-13,i,:)=forePic(:,1,:);
    exforePic(14:end-13,end+1-i,:)=forePic(:,end,:);
end
for i=1:3
    exforePic(1:13,1:13,i)=forePic(1,1,i);
    exforePic(end-13:end,end-13:end,i)=forePic(end,end,i);
    exforePic(end-13:end,1:13,i)=forePic(end,1,i);
    exforePic(1:13,end-13:end,i)=forePic(1,end,i);
end

newforePic=uint8(zeros(size(forePic)));
for i=1:size(bkgPic,1)
    for j=1:size(bkgPic,2)
        goffset=(double(bkgPic(i,j))-128)/10;
        offsetLim1=floor(goffset)+13;
        offsetLim2=ceil(goffset)+13;
        sep1=goffset-floor(goffset);
        sep2=ceil(goffset)-goffset;
        c1=double(exforePic(i+offsetLim1,j+offsetLim1,:));
        c2=double(exforePic(i+offsetLim2,j+offsetLim2,:));
        if sep1==0
            c=double(exforePic(i+offsetLim1,j+offsetLim1,:));
        else
            c=c2.*sep1+c1.*sep2;
        end
        newforePic(i,j,:)=c;
    end
end

在这里插入图片描述
在这里插入图片描述

4图像拼接

就是把俩图像拼在一起，并把边缘模糊一下：

% 图像拼接 ================================================================
resultPic(:,:,1)=[oriPic(:,:,1);newforePic(:,:,1)];
resultPic(:,:,2)=[oriPic(:,:,2);newforePic(:,:,2)];
resultPic(:,:,3)=[oriPic(:,:,3);newforePic(:,:,3)];
% imshow(resultPic)


% 边缘模糊 ================================================================
gaussOpt=fspecial('gaussian',[3 3],0.5);
gaussPic=imfilter(resultPic,gaussOpt);
resultPic(Row-1:Row+2,:,1)=gaussPic(Row-1:Row+2,:,1);
resultPic(Row-1:Row+2,:,2)=gaussPic(Row-1:Row+2,:,2);
resultPic(Row-1:Row+2,:,3)=gaussPic(Row-1:Row+2,:,3);
imshow(resultPic)

在这里插入图片描述

5完整代码

function mirrorDown
% @author slandarer

% 图片导入 
oriPic=imread('test.jpg');
[Row,Col,~]=size(oriPic);

% 图片翻转及白化 ==========================================================
whiteMat=((1:Row)./Row./1.2)'*ones(1,Col); % 白化比例矩阵
flipPic=zeros(Row,Col,3);                  % 翻转后矩阵初始化
for i=1:3
    tempChannel=double(oriPic(:,:,i));     % 获得通道图
    tempChannel=tempChannel(end:-1:1,:);   % 翻转
    tempChannel=tempChannel.*(1-whiteMat)+255.*whiteMat; % 白化
    flipPic(:,:,i)=tempChannel;
    
end
% imshow(uint8(flipPic))


% 噪声图构造(高斯噪声及高斯模糊)===========================================
noiseMat=ones(Row,Col);
noiseMat=imnoise(noiseMat,'gaussian',0,5); % 噪声添加
gaussOpt=fspecial('gaussian',[3 3],1);
noiseMat=imfilter(noiseMat,gaussOpt);
imshow(noiseMat)

H=[cos(pi+pi/4),0,cos(pi-pi/4);
   cos(pi+2*pi/4),0,cos(pi-2*pi/4);
   cos(pi+3*pi/4),0,cos(pi-3*pi/4)];
noiseMat=imfilter(noiseMat,H,'conv')+0.5;  
noiseMat=noiseMat.*255;
noiseMat(noiseMat<0)=0;
% imshow(uint8(noiseMat))


% 图像透视变换 ============================================================
exNoiseMat=zeros(Row,Col);
K1=10;K2=4;
for i=1:Row
    for j=1:Col
        k=K2+i*(K1-K2)/Row;
        nJ=(j-(1+Col)/2)/k+(1+Col)/2;
        if floor(nJ)==ceil(nJ)
            nJ=round(nJ);
            exNoiseMat(i,j)=noiseMat(i,nJ);
        else
            nJ1=floor(nJ);nJ2=ceil(nJ);
            exNoiseMat(i,j)=noiseMat(i,nJ1)*(nJ2-nJ)+noiseMat(i,nJ2)*(nJ-nJ1);
        end
    end
end
exNoiseMat=imresize(exNoiseMat,[6*Row,Col]);
exNoiseMat=exNoiseMat(end-Row+1:end,:);
exNoiseMat=uint8(exNoiseMat);
% imshow(exNoiseMat)


% 扭曲置换 ================================================================
forePic=flipPic;
bkgPic=exNoiseMat;

exforePic=uint8(zeros(size(forePic)+[26,26,0]));
exforePic(14:end-13,14:end-13,1)=forePic(:,:,1);
exforePic(14:end-13,14:end-13,2)=forePic(:,:,2);
exforePic(14:end-13,14:end-13,3)=forePic(:,:,3);

for i=1:13
    exforePic(i,14:end-13,:)=forePic(1,:,:);
    exforePic(end+1-i,14:end-13,:)=forePic(end,:,:);
    exforePic(14:end-13,i,:)=forePic(:,1,:);
    exforePic(14:end-13,end+1-i,:)=forePic(:,end,:);
end
for i=1:3
    exforePic(1:13,1:13,i)=forePic(1,1,i);
    exforePic(end-13:end,end-13:end,i)=forePic(end,end,i);
    exforePic(end-13:end,1:13,i)=forePic(end,1,i);
    exforePic(1:13,end-13:end,i)=forePic(1,end,i);
end

newforePic=uint8(zeros(size(forePic)));
for i=1:size(bkgPic,1)
    for j=1:size(bkgPic,2)
        goffset=(double(bkgPic(i,j))-128)/10;
        offsetLim1=floor(goffset)+13;
        offsetLim2=ceil(goffset)+13;
        sep1=goffset-floor(goffset);
        sep2=ceil(goffset)-goffset;
        c1=double(exforePic(i+offsetLim1,j+offsetLim1,:));
        c2=double(exforePic(i+offsetLim2,j+offsetLim2,:));
        if sep1==0
            c=double(exforePic(i+offsetLim1,j+offsetLim1,:));
        else
            c=c2.*sep1+c1.*sep2;
        end
        newforePic(i,j,:)=c;
    end
end
% imshow(newforePic)


% 图像拼接 ================================================================
resultPic(:,:,1)=[oriPic(:,:,1);newforePic(:,:,1)];
resultPic(:,:,2)=[oriPic(:,:,2);newforePic(:,:,2)];
resultPic(:,:,3)=[oriPic(:,:,3);newforePic(:,:,3)];
% imshow(resultPic)


% 边缘模糊 ================================================================
gaussOpt=fspecial('gaussian',[3 3],0.5);
gaussPic=imfilter(resultPic,gaussOpt);
resultPic(Row-1:Row+2,:,1)=gaussPic(Row-1:Row+2,:,1);
resultPic(Row-1:Row+2,:,2)=gaussPic(Row-1:Row+2,:,2);
resultPic(Row-1:Row+2,:,3)=gaussPic(Row-1:Row+2,:,3);
imshow(resultPic)


end

6其他效果展示

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

slandarer
关注
专栏目录
安装 | 最新MATLAB 2020b（64位）安装教程完整版
尘世冰封的专栏
01-05 17万+
博主github：https://github.com/MichaelBeechan 博主CSDN：https://blog.csdn.net/u011344545 精简安装链接：https://blog.csdn.net/u011344545/article/details/109728674 下载链接 链接：https://pan.baidu.com/s/1g61g2qJO2MtZVkeMsuxhuQ 提取码： 打赏连接：https://blog.csdn.net/u011344545/article/
Matlab实现图像风格化 倒影,浮雕,素描.zip
08-28
Matlab实现图像风格化 倒影,浮雕,素描 代码仅供参考、交流
matlab 图片倒影_计算物理基于matlab方法研究水中倒影问题
weixin_30456153的博客
01-27 469
计算物理基于matlab方法研究水中倒影问题 基于 Matlab 方法研究水中倒影问题[摘 要] 本文介绍了用 matlab 研究倒影问题的方法，利用 matlab 可视化的优点可 以直观的得出结果。 [关键词] matlab； 倒影 ； 高度； 图像 The Research of The Problem about Water Reflection Base On Matlab [Abstra...
matlab 图片倒影_matlab-画个拱桥和倒影？
weixin_39537049的博客
12-19 156
matlab可用于各行各业的应用中，现在我们就简单画一个拱桥试一试。r=input('请输入半径: ');z=input('请输入弧度: ');figure()%r=4;%画拱洞，用李萨如图形法x=15;y=0;nseg=1000;theta = 0 : (2 * pi / nseg) : (2 * pi);pline_x = r * cos(theta) + x;pline_y = r * si...
简单实现下图片倒影效果
阿弥陀佛
04-19 217
简单实现下图片倒影效果 效果图： [img] [img]http://dl.iteye.com/upload/attachment/0066/8127/0acd2e05-3f50-3d44-984c-778f97e14338.jpg[/img] [/img] PictureReflectionActivity [code="java"]package com.amaker....
matlab练习程序（水波特效）
weixin_30755709的博客
12-16 395
还记得原来写过一个对图像进行波纹扭曲操作的博文。 这次实现的是水波特效，其实就是通过正余弦函数表示波纹中心位置慢慢向外扩散，通过叠加衰减因子使振幅不断减小，进而产生水波的效果。 效果如下： 原图： 波纹特效： matlab代码如下： clear all; close all; clc; cenH=128; %波纹中心 cenW=128; R=...
【图像处理】MATLAB应用实战案例：画出水波倒影效果.zip
04-06
本案例中，我们将探讨如何使用MATLAB实现画出水波倒影的效果。水波倒影是图像合成和艺术设计中常见的元素，通过模拟水面波动与物体倒影的交互，可以创造出逼真的视觉效果。 首先，我们需要理解图像处理的基本概念。...
MATLAB使用教程（一）——从零开始，基础，数据类型——新手来看
热门推荐
朱少强的博客
10-11 41万+
matlab是一个很强大的计算软件，合理使用会给我们的学习带来极大的方便之处，本人旨在提供一个零基础的matlab学习博文，供大家学习交流。谢谢大家支持。
daoying.zip_matlab 倒影_倒影
09-24
在本项目中，我们主要探讨如何使用MATLAB编程来创建具有动态倒影效果的图像，这是一个结合了图像处理和交互式编程的实例。MATLAB是一款强大的数学计算软件，也广泛用于图像处理和可视化任务。 首先，我们需要理解...
计算物理——基于matlab方法研究水中倒影问题.doc
11-28
计算物理——基于matlab方法研究水中倒影问题.doc
Matlab动态模拟水波_源码
03-27
通过Matlab程序实现动态模拟水波
水波优化算法用于图像匹配matlab源码
03-21
基于灰色关联分析的图像匹配算法，用水波优化算法进行优化搜索； matlab源码； 注释详细，结构完整；
水波倒影插件
09-13
ps水波倒影插件
水波倒影滤镜
07-19
水波倒影滤镜Flaming Pear 公司的 Flood 图像插件可以做出令人叫绝的倒影效果，个人认为是 Flaming Pear 公司众多滤镜产品中最值得使用的插件。
水波传播matlab代码-Okada:用于计算同震位移的Okada解决方案
06-19
水波传播matlab代码 Okada 解决方案的 Matlab 脚本 此脚本计算由半空间中的位错型变形引起的静态同震位移 描述：这个精确解通常用于构建由水下地震引起的同震位移。 然后使用这些位移以某种方式构建海啸传播代码的初始条件。 有关更多详细信息，请查看原冈田论文： Y.冈田。 由于半空间中的剪切和拉伸断层导致的表面变形。 公牛。 地震。 社会。 上午, 75 , 1135-1154, 1985 或者我免费提供的预印本： D. Dutykh, F. Dias, . 书中：“海啸和非线性波”，Kundu, A.（编辑），Springer Verlag 2007 海啸生成问题的更复杂方法可以在 . 另一个来源： 请注意，此脚本也可以从服务器下载。 ====== 作者 D. Dutykh 执照 下开源。
matlab-画个拱桥和倒影？
weixin_33739646的博客
01-07 333
matlab-画个拱桥和倒影？ matlab可用于各行各业的应用中，现在我们就简单画一个拱桥试一试。 r=input('请输入半径: '); z=input('请输入弧度: '); figure() %r=4; %画拱洞，用...
matlab 图片倒影_ai怎么做水倒影效果？ai制作图案倒影水汶效果教程
weixin_29092031的博客
12-29 781
ai制作图案倒影水汶效果教程：步骤1： 用矩形打开新的画板并创建必须反射的对象,我们可以使用不同的形状和对象进行反射.使用形状,笔,画笔和其他工具创建文本.注意：由于锚点为空,因此可以反映文本步骤2： 选择创建的设计,然后从键盘上单击“ O”以启用反射,按住Alt键并选择锚点以沿相应方向反射对象.双击下面显示的反射工具以打开用于配置值的窗口.由于在底部自然观察到水反射,因此从窗口中选择“水平”,然...
matlab 湖面倒影
01-06
用matlab编写，将一幢房子反射并且模糊，配上水面，加上透明度叠加形成一个倒影，就想湖面倒影，测试过可以使用，项目中提供素材。
评论 13
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

slandarer

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值