目录

• 绘制效果
• 山间秋月
• 皓月当空

中秋节还有三天就到了，中秋节啊，阖家团圆的日子，又有多少人去摆弄电脑甚至打开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转rgbcolSky=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; endhm=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); endend% =========================================================================% 一个线性插值的渐变图生成函数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; endend

文章图片


文章图片

本人将代码再次略作改编，贴合中秋主题，又写了山间秋月的代码，能够动态展示变换的云雾以及慢慢变圆的月亮：

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]; % 天空的颜色% 月亮颜色格式为rgbcMoon=[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=https://www.it610.com/article/[X1,tX2]+600; cloudHdl.AlphaData=cloudsAlpha1+(cloudsAlpha2-cloudsAlpha1).*n./length(k); pause(.1)end% =========================================================================% 功能函数：% -------------------------------------------------------------------------% 渐变背景生成函数function drawSky(colSky,colFurther)% 颜色由hsv转rgbcolSky=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 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); endfunction Alpha=getCloudsAlp(~,~)% 随机噪声生成[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% -------------------------------------------------------------------------% 山峰绘制函数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; endhm=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); endend% =========================================================================% 一个线性插值的渐变图生成函数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; endend

静态图：
【利用Matlab实现绘制中秋山间秋月和皓月当空效果】
文章图片


文章图片


皓月当空 需要下载mapping toolbox或者直接使用文末压缩包内的moonalb20c.mat文件，该文件是真实月球表面数据：

% @author : slandarer% gzh: slandarer随笔% 数据读取load moonalb20c.mat[X,Y,Z0]=sphere(30); surf(X,Y,Z0,'FaceColor','texturemap','CData',moonalb20c,'EdgeColor','none','FaceAlpha',.5)% 调色colormap grayfig=gcf; fig.Color=[0,0,0]; % 旋转 for i=0:(2*pi/200):(4*pi)campos([cos(i),sin(i),.5])axis off vis3dpause(.1)drawnowend

全部文件获取：
链接：https://pan.baidu.com/s/162CUlO9-33SiNx-7ZemVmA
提取码：h71f
到此这篇关于利用Matlab实现绘制中秋山间秋月和皓月当空效果的文章就介绍到这了,更多相关Matlab中秋内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！

推荐阅读

• 

  再看final关键字
• 

  兼职和主业之间的比重该如何安排（）
• 

  佳能cp910照片不完整 佳能cp910照片格式
• 

  单反相机镜头上的字母和数字都代表什么，比如说 佳能 EF-S 18-55mm f/3.5-5.6 II USM 镜头。请具体说明一
• 

  大金空调主板代码故障排除图解 大金空调lc是什么故障 ,要考虑这些方面
• 

  afc电路的分析,Afc电路
• 

  求几本和灵舟遮天类似的小说
• 

  matlab对数据fft分析,MATLAB的FFT分析意义
• 

  孩子|逼孩子叫人、表演节目……最有毒的中国式见面礼，是亲妈就别这样做！
• 

  男人六十岁以后要疏远哪些关系？
• 

  word怎么删除空白页
• 

  荣耀pl2021年高性价比千元机推荐
• 

  豆皮与腐竹有哪些什么区别？
• 

  猫一年能怀几次孕
• 

  体育游戏100例 体育游戏有哪些
• 

  黄冰糖怎么分辨是没有染色的
• 

  喜之郎果冻太空人广告 喜之郎果冻广告
• 

  人生感悟个性签名经典句子有什么？
• 

  桃源深处有人家下一句 桃源深处有人家育种配方有什么
• 

  最终幻想起源全人物介绍 最终幻想起源全角色图鉴分享 杰克

• SpringBoot使用自定义注解+AOP+Redis实现接口限流的实例代码
• QT+ffmpeg实现视频解析的示例详解
• springboot|springboot 集成 docsify 实现随身文档
• #|【目标检测算法】利用wandb可视化YOLO-V5模型的训练
• 多云时代如何实现自动化运维（博云给你最优解！）
• Java并发编程解析 | 解析AQS基础同步器的设计与实现
• 动态调整日志级别思路&实现
• NET6|NET6 IOptionsSnapshot 配合Apollo 实现动态配置 支持自定义解析
• 投稿|互联网的最后一公里，只有卫星才能实现
• Python实现K-近邻算法的示例代码