MATLAB | 中秋节 · 绘制《山间秋月》及《皓月当空》
Posted slandarer
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了MATLAB | 中秋节 · 绘制《山间秋月》及《皓月当空》相关的知识，希望对你有一定的参考价值。
中秋节还有三天就到了，中秋节啊，阖家团圆的日子，又有多少人去摆弄电脑甚至打开MATLAB？因此要发就现在！
绘制效果
月饼绘制，立体月饼绘制，月饼状统计图往年我都已经画过了，那么今年不妨从月亮下手，两个动态绘图奉上。
《山间秋月》
在好久之前看到童晶老师Python趣味创意编程一书中有一山水绘制示例非常惊艳，但用了一些python中有matlab中没有的函数，一时手痒就改了个MATLAB版：
function landScape
% python出处：童晶|《Python趣味创意编程》
% MATLAB代码改写：slandarer
% axes设置
ax=gca;
ax.XTick=[];
ax.YTick=[];
ax.XLim=[0,800];
ax.YLim=[0,600];
ax.DataAspectRatio=[1 1 1];
hold(ax,'on');
% =========================================================================
% 颜色预定义，注意此处是hsv格式
cClouds=[330,25,100]; % 云的颜色
cSky=[220,50,50]; % 天空的颜色
cFurther=[230,25,90]; % 远山的颜色
cCloser=[210,70,10]; % 近山的颜色
% =========================================================================
% 绘图函数调用
ax.Color=hsv2rgb(cFurther./[360,100,100]); % 背景为远山的颜色
drawSky(cSky,cFurther) % 画出天空颜色渐变效果
drawClouds(cClouds) % 画出彩色云朵效果
drawMountains(cFurther,cCloser) % 画出山脉效果
% =========================================================================
% 功能函数：
% -------------------------------------------------------------------------
% 渐变背景生成函数
function drawSky(colSky,colFurther)
% 颜色由hsv转rgb
colSky=hsv2rgb(colSky./[360,100,100]);
colFurther=hsv2rgb(colFurther./[360,100,100]);
%构建渐变色网格
[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
% -------------------------------------------------------------------------
% 云绘制函数
function drawClouds(colClouds)
colClouds=hsv2rgb(colClouds./[360,100,100]);
% 随机噪声生成
[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
% -------------------------------------------------------------------------
% 山峰绘制函数
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
% 每次都生成一次二维随机噪声，并取其中一行的数据
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]);
% 颜色矩阵构造，
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数值框出山的轮廓
alp=ones(hm,800);
alp((1:hm)'>nh)=nan;
% 绘制山峰
image([-50,850],[0,hm],CMesh,'AlphaData',alp.*0.98);
end
end
% =========================================================================
% 一个线性插值的渐变图生成函数
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
本人将代码再次略作改编，贴合中秋主题，又写了山间秋月的代码，能够动态展示变换的云雾以及慢慢变圆的月亮：
function autumoon_2
% @author : slandarer
% gzh : slandarer随笔
% axes设置
ax=gca;
ax.XTick=[];
ax.YTick=[];
ax.XLim=[0,800];
ax.YLim=[0,600];
ax.DataAspectRatio=[1 1 1];
hold(ax,'on');
% =========================================================================
% 颜色预定义，注意此处是hsv格式
cFurther=[225,35,70]; % 远山的颜色
cCloser=[210,70,10]; % 近山的颜色
cClouds=[250 26 43]; % 云的颜色
cSky=[215 100 18]; % 天空的颜色
% 月亮颜色格式为rgb
cMoon=[253,252,222]./255;
% =========================================================================
% 绘图函数调用
ax.Color=hsv2rgb(cFurther./[360,100,100]); % 背景为远山的颜色
drawSky(cSky,cFurther) % 画出天空颜色渐变效果
% 基础绘制月亮
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) % 画出山脉效果
% 随着时间变化月亮逐渐变圆
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
% ==================================================================以上是关于MATLAB | 中秋节 · 绘制《山间秋月》及《皓月当空》的主要内容，如果未能解决你的问题，请参考以下文章