MATLAB绘制超超超炫丽的烟花
一、随机点生成
我们要构造一个黑色背景,对其添加高斯噪声,然后依据阈值删掉部分噪声,生成一张随机点图片:
blackPic=uint8(zeros(800,800));
distPic=imnoise(blackPic,'gaussian',0, 0.1);
distPic(distPic<254)=0;
二、图像膨胀
发现之前构造的点太小了,我们要对其进行形态学膨胀处理:
se=strel('square',3);
distPic=imdilate(distPic,se);
三、特效【风】模拟
就是建立一个循环,不断地将点往右侧复制,并将颜色变暗:
代码中第二个参数为拖拽长度,第三个参数为暗化速度
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;
tempPic(oriPic<tempPic); =
end
resultPic=uint8(oriPic);
end
四、极坐标变换
polarPic=polarTransf(windPic(:,end:-1:1)');
function resultPic=polarTransf(oriPic)
oriPic=double(oriPic);
size(oriPic); =
meshgrid(linspace(-pi,pi,n),1:m); =
M=2*m;
N=2*n;
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
五、图像模糊及再映射
发现烟花图片中心区域方格化严重:
我们很容易想到图像模糊,但是模糊后亮度又不够,因此我们再做一次映射将图片整体亮度提高:
模糊及映射:
polarPic=imgaussfilt(polarPic,1.5);
polarPic=uint8(double(polarPic)./double(max(max(polarPic))).*260);
六、图像上色
matSize=[1600,1600];
point=[800,800];
colorList=[255 253 255
255 255 244
250 255 219
212 242 156
90 167 3
32 96 0
56 82 0];
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;
,:,1)=fwPicR; :
,:,2)=fwPicG; :
,:,3)=fwPicB; :
fwPic=uint8(fwPic);
imshow(fwPic)
function colorMat=cColorMat(matSize,point,colorList)
matSize=[800,600];
point=[400,100];
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);
imshow(colorMat)
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);
,:,1)=colorMesh(end:-1:1,1:matSize(1)); :
,:,2)=colorMesh(end:-1:1,matSize(1)+1:2*matSize(1)); :
,:,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];
源码连接:
https://github.com/gain-wyj/MATLAB_Special_Effects/tree/master
有时你是不是觉得无法进github、谷歌、谷歌学术等等网站,有时会因为翻墙而被封号,给大家推荐一个好的翻墙工具:得自由
海阔凭鱼跃,天高任鸟飞,谁还愿意待在一个“村子”!
得自由(原欢迅干线)致力于建立虚拟专用网络连接方式,加密和认证技术;帮助您突破各样的网络限制,畅快享受自由的网络环境,就这么简单。
通过连接:https://goen.win/c55zhm 进入有优惠哦,包月套餐只需7元,也就一顿饭的钱,然而你可以天高任鸟飞
推荐阅读
(点击标题可跳转阅读)
老铁,三连支持一下,好吗?↓↓↓
![](https://filescdn.proginn.com/19e2ee2670da44441f87ad721bf61a81/5cc47f2a3a3a431115ada619d5b60148.webp)
点分享
点点赞
点在看
评论