欧美bbbwbbbw肥妇,免费乱码人妻系列日韩,一级黄片

利用Matlab實(shí)現(xiàn)繪制中秋山間秋月和皓月當(dāng)空效果

 更新時(shí)間:2022年09月07日 11:57:33   作者:slandarer  
中秋節(jié)還有三天就到了,中秋節(jié)啊,闔家團(tuán)圓的日子。本文將利用Matlab繪制中秋山間秋月和皓月當(dāng)空的動(dòng)態(tài)效果,感興趣的可以了解一下

中秋節(jié)還有三天就到了,中秋節(jié)啊,闔家團(tuán)圓的日子,又有多少人去擺弄電腦甚至打開MATLAB?因此要發(fā)就現(xiàn)在!

繪制效果

月餅繪制,立體月餅繪制,月餅狀統(tǒng)計(jì)圖往年我都已經(jīng)畫過了,那么今年不妨從月亮下手,兩個(gè)動(dòng)態(tài)繪圖奉上。

山間秋月

在好久之前看到童晶老師《Python趣味創(chuàng)意編程》一書中有一山水繪制示例非常驚艷,但用了一些python中有matlab中沒有的函數(shù),一時(shí)手癢就改了個(gè)MATLAB版:

function landScape
% python代碼原作出處:童晶|《Python趣味創(chuàng)意編程》
% MATLAB代碼改寫:slandarer


% axes設(shè)置
ax=gca;
ax.XTick=[];
ax.YTick=[];
ax.XLim=[0,800];
ax.YLim=[0,600];
ax.DataAspectRatio=[1 1 1];
hold(ax,'on');
% =========================================================================
% 顏色預(yù)定義,注意此處是hsv格式
cClouds=[330,25,100];  % 云的顏色
cSky=[220,50,50];      % 天空的顏色
cFurther=[230,25,90];  % 遠(yuǎn)山的顏色
cCloser=[210,70,10];   % 近山的顏色
% =========================================================================
% 繪圖函數(shù)調(diào)用
ax.Color=hsv2rgb(cFurther./[360,100,100]); % 背景為遠(yuǎn)山的顏色
drawSky(cSky,cFurther)                     % 畫出天空顏色漸變效果 
drawClouds(cClouds)                        % 畫出彩色云朵效果
drawMountains(cFurther,cCloser)            % 畫出山脈效果
% =========================================================================
% 功能函數(shù):
% -------------------------------------------------------------------------
% 漸變背景生成函數(shù)
    function drawSky(colSky,colFurther)
        % 顏色由hsv轉(zhuǎn)rgb
        colSky=hsv2rgb(colSky./[360,100,100]);
        colFurther=hsv2rgb(colFurther./[360,100,100]);
        
        %構(gòu)建漸變色網(wǎng)格
        [XMesh,YMesh]=meshgrid(1:800,301:600);
        ZMesh=zeros(size(XMesh));
        CMesh=vColorMat([800,300],[colFurther;colSky]);
        surf(XMesh,YMesh,ZMesh,'CData',CMesh,'EdgeColor','none');
    end
% -------------------------------------------------------------------------
% 云繪制函數(shù)
    function drawClouds(colClouds)
        colClouds=hsv2rgb(colClouds./[360,100,100]);

        % 隨機(jī)噪聲生成
        [X,Y]=meshgrid(linspace(0,1,500));
        CLX=(-cos(X.*2.*pi)+1).^.2;
        CLY=(-cos(Y.*2.*pi)+1).^.2;
        r=(X-.5).^2+(Y-.5).^2;
        alp=abs(ifftn(exp(3i*rand(500))./r.^.8)).*(CLX.*CLY);
        alp=alp./max(alp,[],'all');
        
        CMesh=zeros([size(alp),3]);
        CMesh(:,:,1)=colClouds(1);
        CMesh(:,:,2)=colClouds(2);
        CMesh(:,:,3)=colClouds(3);
        
        % 越向下、云越透明
        dy=(1:500)./500.*0.8+0.2;
        image([0,800],[350,600],CMesh,'AlphaData',alp.*(dy'));
    end
% -------------------------------------------------------------------------
% 山峰繪制函數(shù)
    function drawMountains(colFurther,colCloser)
        [X,Y]=meshgrid(linspace(0,1,800));
        CLX=(-cos(X.*2.*pi)+1).^.2;
        CLY=(-cos(Y.*2.*pi)+1).^.2;
        r=(X-.5).^2+(Y-.5).^2;
        % 8層山
        for i=1:8
            % 每次都生成一次二維隨機(jī)噪聲,并取其中一行的數(shù)據(jù)
            h=abs(ifftn(exp(5i*rand(800))./r.^1.05)).*(CLX.*CLY).*10;
            nh=(8-i)*30+h(400,:);
            if i==1,nh=nh.*.8;end
            hm=ceil(max(nh));
            CMesh=zeros([hm,800,3]);
            
            % 顏色矩陣構(gòu)造,
            tcol=colFurther+(colCloser-colFurther)./8.*(i);
            tcol=hsv2rgb(tcol./[360,100,100]);
            CMesh(:,:,1)=tcol(1);
            CMesh(:,:,2)=tcol(2);
            CMesh(:,:,3)=tcol(3);
            
            % 用nan數(shù)值框出山的輪廓
            alp=ones(hm,800);
            alp((1:hm)'>nh)=nan;
            
            % 繪制山峰
            image([-50,850],[0,hm],CMesh,'AlphaData',alp.*0.98);        
        end
    end
% =========================================================================
% 一個(gè)線性插值的漸變圖生成函數(shù)
    function colorMat=vColorMat(matSize,colorList)
        yList=((0:(matSize(2)-1))./(matSize(2)-1))';
        xList=ones(1,matSize(1));
        % 線性插值
        colorMat(:,:,1)=(colorList(1,1)+yList.*(colorList(2,1)-colorList(1,1)))*xList;
        colorMat(:,:,2)=(colorList(1,2)+yList.*(colorList(2,2)-colorList(1,2)))*xList;
        colorMat(:,:,3)=(colorList(1,3)+yList.*(colorList(2,3)-colorList(1,3)))*xList;
    end
end

本人將代碼再次略作改編,貼合中秋主題,又寫了山間秋月的代碼,能夠動(dòng)態(tài)展示變換的云霧以及慢慢變圓的月亮:

function autumoon_2
% @author : slandarer
% gzh  : slandarer隨筆

% axes設(shè)置
ax=gca;
ax.XTick=[];
ax.YTick=[];
ax.XLim=[0,800];
ax.YLim=[0,600];
ax.DataAspectRatio=[1 1 1];
hold(ax,'on');
% =========================================================================
% 顏色預(yù)定義,注意此處是hsv格式
cFurther=[225,35,70];  % 遠(yuǎn)山的顏色
cCloser=[210,70,10];   % 近山的顏色
cClouds=[250 26 43];   % 云的顏色
cSky=[215 100 18];     % 天空的顏色

% 月亮顏色格式為rgb
cMoon=[253,252,222]./255;
% =========================================================================
% 繪圖函數(shù)調(diào)用
ax.Color=hsv2rgb(cFurther./[360,100,100]); % 背景為遠(yuǎn)山的顏色
drawSky(cSky,cFurther)                     % 畫出天空顏色漸變效果 

% 基礎(chǔ)繪制月亮
t1=linspace(-pi/2,pi/2,100);
t2=linspace(pi/2,3*pi/2,100);
X1=cos(t1).*35;Y1=sin(t1).*35;
X2=cos(t2).*35;Y2=sin(t2).*35;
moonHdl=fill([X1,X2]+600,[Y1,Y2]+500,cMoon,'EdgeColor','none');

% 畫出彩色云朵效果                       
cloudsCMesh=getCloudsCMesh(cClouds);
cloudsAlpha1=getCloudsAlp();
cloudsAlpha2=getCloudsAlp();
cloudHdl=image([0,800],[250,600],cloudsCMesh,'AlphaData',cloudsAlpha1);


drawMountains(cFurther,cCloser)            % 畫出山脈效果
% 隨著時(shí)間變化月亮逐漸變圓
k=linspace(-1,1,100);
for n=1:length(k)
    tX2=X2.*k(n);
    moonHdl.XData=[X1,tX2]+600;
    cloudHdl.AlphaData=cloudsAlpha1+(cloudsAlpha2-cloudsAlpha1).*n./length(k);
    pause(.1)
end
% =========================================================================
% 功能函數(shù):
% -------------------------------------------------------------------------
% 漸變背景生成函數(shù)
    function drawSky(colSky,colFurther)
        % 顏色由hsv轉(zhuǎn)rgb
        colSky=hsv2rgb(colSky./[360,100,100]);
        colFurther=hsv2rgb(colFurther./[360,100,100]);
        
        %構(gòu)建漸變色網(wǎng)格
        [XMesh,YMesh]=meshgrid(1:800,301:600);
        ZMesh=zeros(size(XMesh));
        CMesh=vColorMat([800,300],[colFurther;colSky]);
        surf(XMesh,YMesh,ZMesh,'CData',CMesh,'EdgeColor','none');
    end
% -------------------------------------------------------------------------
% 云繪制函數(shù)
    function CMesh=getCloudsCMesh(colClouds)
        colClouds=hsv2rgb(colClouds./[360,100,100]);
        CMesh=zeros([500,500,3]);
        CMesh(:,:,1)=colClouds(1);
        CMesh(:,:,2)=colClouds(2);
        CMesh(:,:,3)=colClouds(3);
    end
    function Alpha=getCloudsAlp(~,~)
        % 隨機(jī)噪聲生成
        [X,Y]=meshgrid(linspace(0,1,500));
        CLX=(-cos(X.*2.*pi)+1).^.2;
        CLY=(-cos(Y.*2.*pi)+1).^.2;
        r=(X-.5).^2+(Y-.5).^2;
        alp=abs(ifftn(exp(3i*rand(500))./r.^.8)).*(CLX.*CLY);
        alp=alp./max(alp,[],'all');
        dy=(1:500)./500.*0.8+0.2;
        Alpha=alp.*(dy');
    end
% -------------------------------------------------------------------------
% 山峰繪制函數(shù)
    function drawMountains(colFurther,colCloser)
        [X,Y]=meshgrid(linspace(0,1,800));
        CLX=(-cos(X.*2.*pi)+1).^.2;
        CLY=(-cos(Y.*2.*pi)+1).^.2;
        r=(X-.5).^2+(Y-.5).^2;
        % 8層山
        for i=1:8
            % 每次都生成一次二維隨機(jī)噪聲,并取其中一行的數(shù)據(jù)
            h=abs(ifftn(exp(5i*rand(800))./r.^1.05)).*(CLX.*CLY).*10;
            nh=(8-i)*30+h(400,:);
            if i==1,nh=nh.*.8;end
            hm=ceil(max(nh));
            CMesh=zeros([hm,800,3]);
            
            % 顏色矩陣構(gòu)造
            tcol=colFurther+(colCloser-colFurther)./8.*(i);
            tcol=hsv2rgb(tcol./[360,100,100]);
            CMesh(:,:,1)=tcol(1);
            CMesh(:,:,2)=tcol(2);
            CMesh(:,:,3)=tcol(3);
            
            % 用nan數(shù)值框出山的輪廓
            alp=ones(hm,800);
            alp((1:hm)'>nh)=nan;
            
            % 繪制山峰
            image([-50,850],[0,hm],CMesh,'AlphaData',alp.*0.98);        
        end
    end
% =========================================================================
% 一個(gè)線性插值的漸變圖生成函數(shù)
    function colorMat=vColorMat(matSize,colorList)
        yList=((0:(matSize(2)-1))./(matSize(2)-1))';
        xList=ones(1,matSize(1));
        % 線性插值
        colorMat(:,:,1)=(colorList(1,1)+yList.*(colorList(2,1)-colorList(1,1)))*xList;
        colorMat(:,:,2)=(colorList(1,2)+yList.*(colorList(2,2)-colorList(1,2)))*xList;
        colorMat(:,:,3)=(colorList(1,3)+yList.*(colorList(2,3)-colorList(1,3)))*xList;
    end
end

靜態(tài)圖:

皓月當(dāng)空

需要下載mapping toolbox或者直接使用文末壓縮包內(nèi)的moonalb20c.mat文件,該文件是真實(shí)月球表面數(shù)據(jù):

% @author : slandarer
% gzh  : slandarer隨筆

% 數(shù)據(jù)讀取
load moonalb20c.mat

[X,Y,Z0]=sphere(30);
surf(X,Y,Z0,'FaceColor','texturemap','CData',moonalb20c,'EdgeColor','none','FaceAlpha',.5)

% 調(diào)色
colormap gray
fig=gcf;
fig.Color=[0,0,0];
 
% 旋轉(zhuǎn) 
for i=0:(2*pi/200):(4*pi)
    campos([cos(i),sin(i),.5])
    axis off vis3d
    pause(.1)
    drawnow
end

全部文件獲取:

鏈接:https://pan.baidu.com/s/162CUlO9-33SiNx-7ZemVmA

提取碼:h71f

到此這篇關(guān)于利用Matlab實(shí)現(xiàn)繪制中秋山間秋月和皓月當(dāng)空效果的文章就介紹到這了,更多相關(guān)Matlab中秋內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!

相關(guān)文章

  • C++實(shí)現(xiàn)簡(jiǎn)單貪吃蛇游戲

    C++實(shí)現(xiàn)簡(jiǎn)單貪吃蛇游戲

    這篇文章主要為大家詳細(xì)介紹了C++實(shí)現(xiàn)簡(jiǎn)單貪吃蛇游戲,文中示例代碼介紹的非常詳細(xì),具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下
    2018-03-03
  • win32 api實(shí)現(xiàn)2048游戲示例

    win32 api實(shí)現(xiàn)2048游戲示例

    這篇文章主要介紹了win32 api實(shí)現(xiàn)2048游戲示例,需要的朋友可以參考下
    2014-05-05
  • C++設(shè)計(jì)模式之解釋器模式

    C++設(shè)計(jì)模式之解釋器模式

    這篇文章主要介紹了C++設(shè)計(jì)模式之解釋器模式,本文講解了什么是解釋器模式、文法規(guī)則和抽象語法樹、解釋器模式的使用場(chǎng)合等內(nèi)容,需要的朋友可以參考下
    2014-10-10
  • 深入了解C語言結(jié)構(gòu)化的程序設(shè)計(jì)

    深入了解C語言結(jié)構(gòu)化的程序設(shè)計(jì)

    這篇文章主要介紹了C語言編程中程序的一些基本的編寫優(yōu)化技巧,文中涉及到了基礎(chǔ)的C程序內(nèi)存方面的知識(shí),非常推薦!需要的朋友可以參考下
    2021-07-07
  • VScode platformio使用的詳細(xì)步驟

    VScode platformio使用的詳細(xì)步驟

    使用VSCode作為編輯器,同時(shí)借助PlatformIO插件,可以幫助開發(fā)者更加高效地進(jìn)行嵌入式開發(fā),本文主要介紹了VScode platformio使用的詳細(xì)步驟,感興趣的可以了解一下
    2023-10-10
  • C++中hashmap的一些使用建議

    C++中hashmap的一些使用建議

    由于hashmap不是c++ stl中標(biāo)準(zhǔn)實(shí)現(xiàn),這樣在跨平臺(tái)使用時(shí)就可能會(huì)出現(xiàn)問題,下面這篇文章主要給大家介紹了關(guān)于C++中hashmap的一些使用建議,需要的朋友可以參考下
    2023-03-03
  • C語言編程中實(shí)現(xiàn)二分查找的簡(jiǎn)單入門實(shí)例

    C語言編程中實(shí)現(xiàn)二分查找的簡(jiǎn)單入門實(shí)例

    這篇文章主要介紹了C語言編程中實(shí)現(xiàn)二分查找的簡(jiǎn)單入門實(shí)例,需要的朋友可以參考下
    2015-12-12
  • 你只用do-while來實(shí)現(xiàn)循環(huán)?太浪費(fèi)了

    你只用do-while來實(shí)現(xiàn)循環(huán)?太浪費(fèi)了

    這篇文章主要介紹了你只用do-while來實(shí)現(xiàn)循環(huán)?太浪費(fèi)了,本文給大家介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下
    2020-12-12
  • C語言快速排序與二分查找算法示例

    C語言快速排序與二分查找算法示例

    這篇文章主要介紹了C語言快速排序與二分查找算法,涉及C語言隨機(jī)數(shù)生成、快速排序及二分查找等算法相關(guān)實(shí)現(xiàn)技巧,需要的朋友可以參考下
    2018-01-01
  • C語言實(shí)現(xiàn)簡(jiǎn)單掃雷游戲

    C語言實(shí)現(xiàn)簡(jiǎn)單掃雷游戲

    這篇文章主要為大家詳細(xì)介紹了C語言實(shí)現(xiàn)簡(jiǎn)單掃雷游戲,文中示例代碼介紹的非常詳細(xì),具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下
    2020-02-02

最新評(píng)論