用MATLAB一直画直线竟然得到了简单的禅绕画，论一直画直线的人有多无聊

article/2025/10/5 1:47:45

之前看过一个视频，有个人把纸面分割成一个一个三角后一直画直线，慢慢的图纸上的图案变得复杂了起来，大概是像下面这样：

在这里插入图片描述

在这里插入图片描述
上面这个图便是我用matlab模拟的效果，过程很简单，就是用了泊松云盘采点构建三角网格，之后一直在三角型内部画直线即可，

泊松云盘采点可以看一下图片三角化那篇：
MATLAB 图片三角风格化(low poly)

三角型画线可以参照一下这一篇：
matlab 旋转的三角形

参数说明

rows=560; 画布高度
cols=750; 画布宽度
r=min([rows,cols])./6; 两点之间距离
LW=1.2; 线条初始粗细
k1=1/8; 取三角型一顶点及另外两点连线1/8处两点做新的连线
k2=0.8; 线条变细的速度
t=20; 每个三角型内画3*t个线条

例如将参数r改为r=min([rows,cols])./3;则得到下图：
在这里插入图片描述

完整代码

function tris
rows=560;
cols=750;
r=min([rows,cols])./6;
LW=1.2;
k1=1/8;
k2=0.8;
t=20;fig=figure('units','pixels',...'position',[20 60 (cols+1) (rows+1)],...'Color',[1 1 1]);
ax=axes('Units','pixels',...'parent',fig,...  'Color',[1 1 1],...'Position',[0 0 (cols+1) (rows+1)],...'XLim',[-1 cols+1],...'YLim',[-1 rows+1],...'XColor','none',...'YColor','none');
hold(ax,'on')
bkgMat=ones([rows,cols]);
bkgMat(2:(rows-1),2:(cols-1))=0;
bkgMat(1,cols)=0;bkgMat(1,1)=0;bkgMat(rows,1)=0;bkgMat(rows,cols)=0;
[edgeX,edgeY]=find(bkgMat==1);
edgePntList=[edgeY,edgeX];edgePntList=poissonEdge(rows,cols,edgePntList,r);
pntList=poissonDisk(rows,cols,r,30,edgePntList);
DT=delaunay(pntList(:,1),pntList(:,2));
triplot(DT,pntList(:,1),pntList(:,2),'Color',[0,0,0],'Linewidth',LW);
pause(0.1)for k=1:size(DT,1)A=pntList(DT(k,1),:);B=pntList(DT(k,2),:);C=pntList(DT(k,3),:);drawABC(A,B,C,LW,k1,k2,t)
end%==========================================================================function resultSet=poissonEdge(m,n,edgeList,R)preSet=edgeList;resultSet=[1,1;n,m;1,m;n,1];times=0;while times<150tempPos=randi([1,size(preSet,1)],1);selectedPnt=preSet(tempPos,:);dis=sqrt(sum((edgeList-selectedPnt).^2,2));candidate=dis>=R&dis<=2*R;pntSet=edgeList(candidate,:);flag=0;for j=1:size(pntSet,1)pnt=pntSet(j,:);dis=sqrt(sum((resultSet-pnt).^2,2));if all(dis>=R)resultSet=[resultSet;pnt];preSet=[preSet;pnt];flag=1;endendif flag==1preSet(tempPos,:)=[];times=0;elsetimes=times+1;end endendfunction resultSet=poissonDisk(m,n,r,K,edgePntList)preSet=edgePntList;resultSet=edgePntList;times=0;while times<500tempPos=randi([1,size(preSet,1)],1);selectedPnt=preSet(tempPos,:);theta=rand(K,1).*2*pi;radius=rand(K,1).*r+r;x=radius.*cos(theta)+selectedPnt(1);y=radius.*sin(theta)+selectedPnt(2);flag=0;for j=1:Kpnt=[x(j),y(j)];if pnt(1)>=1&&pnt(2)>=1&&pnt(1)<=n&&pnt(2)<=mif size(resultSet,1)==0resultSet=[resultSet;pnt];preSet=[preSet;pnt];flag=1;elsedis=sqrt(sum((resultSet-pnt).^2,2));if all(dis>=r)resultSet=[resultSet;pnt];preSet=[preSet;pnt];flag=1;   endend endendif flag==1preSet(tempPos,:)=[];times=0;elsetimes=times+1;endend    endfunction drawABC(A,B,C,LW,k1,k2,times)for i=1:timesb=k1*(C-B)+B;plot([A(1);b(1)],[A(2);b(2)],'color',[0,0,0],'Linewidth',LW);B=b;pause(0.001)c=k1*(A-C)+C;plot([B(1);c(1)],[B(2);c(2)],'color',[0,0,0],'Linewidth',LW);C=c;pause(0.001)a=k1*(B-A)+A;plot([C(1);a(1)],[C(2);a(2)],'color',[0,0,0],'Linewidth',LW);A=a;pause(0.001)LW=LW*k2;endendend