基于Matlab繪制超絢麗的煙花的過程詳解
更新時間：2022年02月25日 10:13:13 作者：slandarer
這篇文章主要為大家介紹了如何利用Matlab繪制超絢麗的煙花，文中的示例代碼講解詳細，對我們學習Matlab有一定幫助，需要的可以參考一下
1.使用效果
2.隨機點生成
我們要構(gòu)造一個黑色背景，對其添加高斯噪聲，然后依據(jù)閾值刪掉部分噪聲，生成一張隨機點圖片：
% 構(gòu)造黑色背景并生成白色雜點 blackPic=uint8(zeros(800,800)); distPic=imnoise(blackPic,'gaussian',0, 0.11); distPic(distPic<254)=0;
3.圖像膨脹
我們發(fā)現(xiàn)之前構(gòu)造的點太小了，我們要對其進行形態(tài)學膨脹處理：
% 雜點膨脹
se=strel('square',3);
distPic=imdilate(distPic,se);
4.特效「風」模擬
就是建立一個循環(huán)，不斷地將點往右側(cè)復制，并將顏色變暗，代碼中第二個參數(shù)為拖拽長度，第三個參數(shù)為暗化速度：
% 構(gòu)造風特效
windPic=wind(distPic,180,0.99);
% 風特效構(gòu)造函數(shù)
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)');
% 級坐標變換構(gòu)造函數(shù)
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.圖像模糊及再映射
我們發(fā)現(xiàn)煙花圖片中心區(qū)域方格化嚴重：
我們很容易想到圖像模糊，但是模糊后亮度又不夠，因此我們再做一次映射將圖片整體亮度提高：
模糊及映射：
% 模糊以減少像素化 polarPic=imgaussfilt(polarPic,1.5); polarPic=uint8(double(polarPic)./double(max(max(polarPic))).*260);
7.圖像上色
生成一個圓形的漸變圖片，并與白色煙花圖片進行一次正交疊底：
% 構(gòu)造圓形漸變圖
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)
%==========================================================================
% 圖像顏色映射函數(shù)
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
% 構(gòu)造黑色背景并生成白色雜點
blackPic=uint8(zeros(800,800));
distPic=imnoise(blackPic,'gaussian',0, 0.11);distPic(distPic<254)=0;
% 雜點膨脹
se=strel('square',3);
distPic=imdilate(distPic,se);
% 構(gòu)造風特效
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);
%-----------------------------------------------------------------
% 構(gòu)造圓形漸變圖
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)
%==========================================================================
% 風特效構(gòu)造函數(shù)
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
% 極坐標變換構(gòu)造函數(shù)
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
%==========================================================================
% 圖像顏色映射函數(shù)
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繪制超絢麗的煙花的過程詳解的詳細內(nèi)容，更多關(guān)于Matlab繪制煙花的資料請關(guān)注腳本之家其它相關(guān)文章！
相關(guān)文章
C++中類的成員函數(shù)及內(nèi)聯(lián)函數(shù)使用及說明
這篇文章主要介紹了C++中類的成員函數(shù)及內(nèi)聯(lián)函數(shù)使用及說明，具有很好的參考價值，希望對大家有所幫助。如有錯誤或未考慮完全的地方，望不吝賜教2022-11-11
本文主要介紹了vscode配置gitbash終端的方法步驟，文中通過示例代碼介紹的非常詳細，對大家的學習或者工作具有一定的參考學習價值，需要的朋友們下面隨著小編來一起學習學習吧2023-06-06
本文主要介紹了Qt視頻播放器的實現(xiàn)示例,文中通過示例代碼介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友們下面隨著小編來一起學習學習吧2023-08-08
這篇文章主要為大家詳細介紹了C語言實現(xiàn)愛情紅玫瑰，文中示例代碼介紹的非常詳細，具有一定的參考價值，感興趣的小伙伴們可以參考一下2020-08-08
C++數(shù)據(jù)結(jié)構(gòu)之鏈表詳解
這篇文章主要介紹了C++數(shù)據(jù)結(jié)構(gòu)之鏈表的創(chuàng)建的相關(guān)資料,希望通過本文幫助到大家，讓大家理解掌握這部分內(nèi)容，需要的朋友可以參考下2021-08-08
Linux系統(tǒng)下C語言gets函數(shù)出現(xiàn)警告問題的解決方法
這篇文章主要給大家介紹了關(guān)于在Linux系統(tǒng)下C語言gets函數(shù)出現(xiàn)警告問題的解決方法，文中通過示例代碼介紹的非常詳細，對大家的學習或者工作具有一定的參考學習價值，需要的朋友們下面隨著小編來一起學習學習吧2018-12-12
這篇文章主要介紹了C語言構(gòu)建動態(tài)數(shù)組完整實例,幫助讀者加深對C語言數(shù)組及指針的理解,需要的朋友可以參考下2014-07-07
這篇文章主要介紹了C++的最短路徑的弗洛伊德算法案例講解,本篇文章通過簡要的案例,講解了該項技術(shù)的了解與使用,以下就是詳細內(nèi)容,需要的朋友可以參考下2021-08-08
這篇文章主要介紹了C語言數(shù)據(jù)輸入與輸出實例詳解的相關(guān)資料,需要的朋友可以參考下2017-06-06
這篇文章主要為大家詳細介紹了C語言實現(xiàn)圖書管理系統(tǒng)，文中示例代碼介紹的非常詳細，具有一定的參考價值，感興趣的小伙伴們可以參考一下2018-01-01