基于Matlab绘制超绚丽的烟花的过程详解
更新时间：2022年02月25日 10:13:13 作者：slandarer
这篇文章主要为大家介绍了如何利用Matlab绘制超绚丽的烟花，文中的示例代码讲解详细，对我们学习Matlab有一定帮助，需要的可以参考一下
1.使用效果
2.随机点生成
我们要构造一个黑色背景，对其添加高斯噪声，然后依据阈值删掉部分噪声，生成一张随机点图片：
% 构造黑色背景并生成白色杂点 blackPic=uint8(zeros(800,800)); distPic=imnoise(blackPic,'gaussian',0, 0.11); distPic(distPic<254)=0;
3.图像膨胀
我们发现之前构造的点太小了，我们要对其进行形态学膨胀处理：
% 杂点膨胀
se=strel('square',3);
distPic=imdilate(distPic,se);
4.特效「风」模拟
就是建立一个循环，不断地将点往右侧复制，并将颜色变暗，代码中第二个参数为拖拽长度，第三个参数为暗化速度：
% 构造风特效
windPic=wind(distPic,180,0.99);
% 风特效构造函数
function resultPic=wind(oriPic,len,ratio)
oriPic=double(oriPic);
for i=1:len
tempPic=[zeros(size(oriPic,1),1),oriPic(:,1:(end-1))].*ratio;
oriPic(oriPic<tempPic)=tempPic(oriPic<tempPic);
end
resultPic=uint8(oriPic);
end
5.级坐标变换
% 极坐标变换
polarPic=polarTransf(windPic(:,end:-1:1)');
% 级坐标变换构造函数
function resultPic=polarTransf(oriPic)
oriPic=double(oriPic);
[m,n]=size(oriPic);
[t,r]=meshgrid(linspace(-pi,pi,n),1:m);
M=2*m;
N=2*n;
[NN,MM]=meshgrid((1:N)-n-0.5,(1:M)-m-0.5);
T=atan2(NN,MM);
R=sqrt(MM.^2+NN.^2);
resultPic=interp2(t,r,oriPic,T,R,'linear',0);
resultPic=uint8(resultPic);
end
6.图像模糊及再映射
我们发现烟花图片中心区域方格化严重：
我们很容易想到图像模糊，但是模糊后亮度又不够，因此我们再做一次映射将图片整体亮度提高：
模糊及映射：
% 模糊以减少像素化 polarPic=imgaussfilt(polarPic,1.5); polarPic=uint8(double(polarPic)./double(max(max(polarPic))).*260);
7.图像上色
生成一个圆形的渐变图片，并与白色烟花图片进行一次正交叠底：
% 构造圆形渐变图
matSize=[1600,1600];
point=[800,800];
colorList=[195 53 93
211 102 141
231 179 192
229 182 172
227 178 137
238 191 147
236 195 113];
% 正片叠底
colorMat=cColorMat(matSize,point,colorList);
colorMatR=colorMat(:,:,1);
colorMatG=colorMat(:,:,2);
colorMatB=colorMat(:,:,3);
fwPicR=double(colorMatR).*double(polarPic)./255;
fwPicG=double(colorMatG).*double(polarPic)./255;
fwPicB=double(colorMatB).*double(polarPic)./255;
fwPic(:,:,1)=fwPicR;
fwPic(:,:,2)=fwPicG;
fwPic(:,:,3)=fwPicB;
fwPic=uint8(fwPic);
imshow(fwPic)
%==========================================================================
% 图像颜色映射函数
function colorMat=cColorMat(matSize,point,colorList)
[xMesh,yMesh]=meshgrid(1:matSize(2),1:matSize(1));
zMesh=sqrt((xMesh-point(2)).^2+(yMesh-point(1)).^2);
zMesh=(zMesh-min(min(zMesh)))./(max(max(zMesh))-min(min(zMesh)));
colorFunc=colorFuncFactory(colorList);
colorMesh=colorFunc(zMesh);
colorMat(:,:,1)=colorMesh(end:-1:1,1:matSize(1));
colorMat(:,:,2)=colorMesh(end:-1:1,matSize(1)+1:2*matSize(1));
colorMat(:,:,3)=colorMesh(end:-1:1,2*matSize(1)+1:3*matSize(1));
colorMat=uint8(colorMat);
end
function colorFunc=colorFuncFactory(colorList)
x=(0:size(colorList,1)-1)./(size(colorList,1)-1);
y1=colorList(:,1);y2=colorList(:,2);y3=colorList(:,3);
colorFunc=@(X)[interp1(x,y1,X,'linear')',...
interp1(x,y2,X,'linear')',...
interp1(x,y3,X,'linear')'];
end
其他几个颜色：
colorList2=[25 59 157
24 71 219
38 124 237
93 215 255
168 244 255
243 254 250
246 252 240];
colorList3=[239 250 210
229 164 122
232 150 138
255 164 204
192 58 111
158 10 26
224 168 121];8.完整代码
function drawFireWorks
% 构造黑色背景并生成白色杂点
blackPic=uint8(zeros(800,800));
distPic=imnoise(blackPic,'gaussian',0, 0.11);distPic(distPic<254)=0;
% 杂点膨胀
se=strel('square',3);
distPic=imdilate(distPic,se);
% 构造风特效
windPic=wind(distPic,180,0.99);
% 极坐标变换
polarPic=polarTransf(windPic(:,end:-1:1)');
% 模糊以减少像素化
polarPic=imgaussfilt(polarPic,1.5);
polarPic=uint8(double(polarPic)./double(max(max(polarPic))).*260);
%-----------------------------------------------------------------
% 构造圆形渐变图
matSize=[1600,1600];
point=[800,800];
colorList=[195 53 93
211 102 141
231 179 192
229 182 172
227 178 137
238 191 147
236 195 113];
% 正片叠底
colorMat=cColorMat(matSize,point,colorList);
colorMatR=colorMat(:,:,1);
colorMatG=colorMat(:,:,2);
colorMatB=colorMat(:,:,3);
fwPicR=double(colorMatR).*double(polarPic)./255;
fwPicG=double(colorMatG).*double(polarPic)./255;
fwPicB=double(colorMatB).*double(polarPic)./255;
fwPic(:,:,1)=fwPicR;
fwPic(:,:,2)=fwPicG;
fwPic(:,:,3)=fwPicB;
fwPic=uint8(fwPic);
imshow(fwPic)
%==========================================================================
% 风特效构造函数
function resultPic=wind(oriPic,len,ratio)
oriPic=double(oriPic);
for i=1:len
tempPic=[zeros(size(oriPic,1),1),oriPic(:,1:(end-1))].*ratio;
oriPic(oriPic<tempPic)=tempPic(oriPic<tempPic);
end
resultPic=uint8(oriPic);
end
% 极坐标变换构造函数
function resultPic=polarTransf(oriPic)
oriPic=double(oriPic);
[m,n]=size(oriPic);
[t,r]=meshgrid(linspace(-pi,pi,n),1:m);
M=2*m;
N=2*n;
[NN,MM]=meshgrid((1:N)-n-0.5,(1:M)-m-0.5);
T=atan2(NN,MM);
R=sqrt(MM.^2+NN.^2);
resultPic=interp2(t,r,oriPic,T,R,'linear',0);
resultPic=uint8(resultPic);
end
%==========================================================================
% 图像颜色映射函数
function colorMat=cColorMat(matSize,point,colorList)
[xMesh,yMesh]=meshgrid(1:matSize(2),1:matSize(1));
zMesh=sqrt((xMesh-point(2)).^2+(yMesh-point(1)).^2);
zMesh=(zMesh-min(min(zMesh)))./(max(max(zMesh))-min(min(zMesh)));
colorFunc=colorFuncFactory(colorList);
colorMesh=colorFunc(zMesh);
colorMat(:,:,1)=colorMesh(end:-1:1,1:matSize(1));
colorMat(:,:,2)=colorMesh(end:-1:1,matSize(1)+1:2*matSize(1));
colorMat(:,:,3)=colorMesh(end:-1:1,2*matSize(1)+1:3*matSize(1));
colorMat=uint8(colorMat);
end
function colorFunc=colorFuncFactory(colorList)
x=(0:size(colorList,1)-1)./(size(colorList,1)-1);
y1=colorList(:,1);y2=colorList(:,2);y3=colorList(:,3);
colorFunc=@(X)[interp1(x,y1,X,'linear')',...
interp1(x,y2,X,'linear')',...
interp1(x,y3,X,'linear')'];
end
end以上就是基于Matlab绘制超绚丽的烟花的过程详解的详细内容，更多关于Matlab绘制烟花的资料请关注脚本之家其它相关文章！
相关文章
今天小编就为大家分享一篇opencv3/C++ 实现SURF特征检测，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2019-12-12
这篇文章主要为大家详细介绍了C语言实现简单通讯录管理系统，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2021-07-07
本篇文章是对Visual C++进行串口通信编程进行了详细的分析介绍，需要的朋友参考下2013-05-05
这篇文章主要介绍了C++基础教程之指针拷贝详解的相关资料,需要的朋友可以参考下2017-01-01
这篇文章主要为大家详细介绍了使用C语言实现三子棋小游戏，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2022-08-08
对比分析C语言中的gcvt()和ecvt()以及fcvt()函数
这篇文章主要介绍了对比分析C语言中的gcvt和ecvt以及fcvt函数,都是将数字转化为字符串,注意其之间的功能区别,需要的朋友可以参考下2015-08-08
C语言的最大特点是：功能强、使用方便灵活。C编译的程序对语法检查并不象其它高级语言那么严格，这就给编程人员留下“灵活的余地”，但还是由于这个灵活给程序的调试带来了许多不便，尤其对初学C语言的人来说，经常会出一些连自己都不知道错在哪里的错误2013-07-07
这篇文章主要介绍了C++实现DES加密算法实例解析,是一个很实用的功能,需要的朋友可以参考下2014-08-08
这篇文章主要介绍了C语言中的参数传递机制,C语言中函数参数的传递有：值传递、地址传递、引用传递这三种形式。下面我们详细探讨下2017-04-04
树是计算机算法最重要的非线性结构。因为树能很好地描述结构的分支关系和层次特性，所以在计算机科学和计算机应用领域有着广泛的应用。这篇文章我就带大家一起了解一下树、二叉树这种结构，下篇文章会重点向大家介绍二叉树的遍历算法2022-05-05
最新评论