MATLAB 简易JK裙绘制器

另制作了更立体版本，将其放在文章最后，效果如下：

---

做了一个简易JK裙绘制器
使用按键R选择贴图
推荐尺寸1:1或3:5
绘制效果：



当然也可以来些恐怖的配色。。。。

完整代码:

function fillJK
fig=figure('units','pixels',...
        'position',[350 100 500 500],...
        'Numbertitle','off',...
        'name','JK',...
        'Color',[1 1 1],...
        'resize','off',...
         'menubar','none');

ax=axes('Units','pixels',...
        'parent',fig,...  
        'Color',[1 1 1],...
        'Position',[0 0 500 500],...
        'XLim',[0 100],...
        'YLim',[0 100],...
        'XGrid','on','YGrid','on');
hold(ax,'on')   
[Lu_X,Lu_Z]=pwLine(0:0.1:120,7,3,asin(0.1/3));
[Ld_X,Ld_Z]=pwLine(0:0.1:120,9,1,asin(0.1/1));
Lu_Y=(60+20).*ones(1,length(0:0.1:120));
Ld_Y=20.*ones(1,length(0:0.1:120));

Lu_X=Lu_X+(max(Ld_X)-max(Lu_X))/2;
%plot3(Lu_X,Lu_Y,Lu_Z)
%plot3(Ld_X,Ld_Y,Ld_Z)

Xset=inter3D(Lu_X',Ld_X',length(Lu_X));
Yset=inter3D(Lu_Y',Ld_Y',length(Lu_Y));
Zset=inter3D(Lu_Z',Ld_Z',length(Lu_Z));

Pic=[];
try
    Pic=imread('示例.png');
catch
end
h=surface(ax,Xset,Yset,Zset,'FaceColor','texturemap','EdgeColor','texturemap');
if isempty(Pic)
else
    PicR=Pic(:,:,1);PicR=PicR';
    PicG=Pic(:,:,2);PicG=PicG';
    PicB=Pic(:,:,3);PicB=PicB';
    tpPic(:,:,1)=PicR;
    tpPic(:,:,2)=PicG;
    tpPic(:,:,3)=PicB;
    w=fspecial('gaussian',[3 3],3);
    flPic=imfilter(tpPic,w);
    set(h,'CData',flPic,'FaceColor','texturemap','EdgeColor','texturemap')
end
lightangle(ax,-10,61)
lighting(ax,'gouraud')
shading(ax,'flat')
text(2,96,'JK裙贴图展示器','FontSize',11)
text(2,92,'通过按键"R"更换贴图,建议比例3:5','FontSize',11)

set(gcf, 'KeyPressFcn', @key) 
function key(~,event)
    if event.Key=='r'
        try
        [filename, pathname] = uigetfile({
     '*.jpg;*.tif;*.png;*.gif','All Image Files';...
            '*.*','All Files' });
        Pic=[];
        Pic=imread([ pathname,filename]);
        PicR=Pic(:,:,1);PicR=PicR';
        PicG=Pic(:,:,2);PicG=PicG';
        PicB=Pic(:,:,3);PicB=PicB';
        if size(Pic,1)==size(Pic,2)
            PicR=repmat(PicR,[5,3]);
            PicG=repmat(PicG,[5,3]);
            PicB=repmat(PicB,[5,3]);
        end
        tpPic(:,:,1)=PicR;
        tpPic(:,:,2)=PicG;
        tpPic(:,:,3)=PicB;
        w=fspecial('gaussian',[3 3],3);
        flPic=imfilter(tpPic,w);
        set(h,'CData',flPic)
        catch
        end    
    end
end


    function Pset=inter3D(P1,P2,N)
        coeMat=linspace(0,1,N);
        Pset=P1+(P2-P1)*coeMat;
    end

    function [newX,newY]=pwLine(X,p1L,p2L,theta)      
        newX=X(1);
        newY=0;
        for i=2:length(X)
            dX=X(i)-X(i-1);
            if mod(X(i),p1L+p2L)<=p1L
                newX=[newX,newX(end)+dX];
                newY=[newY,newY(end)-0];
            else
                newX=[newX,newX(end)-dX*cos(theta)];
                newY=[newY,newY(end)-dX*sin(theta)];
            end
        end
    end


end

贴图的绘制当然要参考前面写的那篇，格柄生成器的文章啦：
MATLAB 简易格柄生成器

---

更立体版本代码：

function fillJK2
fig=figure('units','pixels',...
        'position',[350 100 500 500],...
        'Numbertitle','off',...
        'name','JK',...
        'Color',[1 1 1],...
        'resize','off');%,...
%         'menubar','none');

ax=axes('Units','pixels',...
        'DataAspectRatioMode','manual',...
        'parent',fig,...  
        'Color',[1 1 1],...
        'Position',[0 0 500 500],...
        'XLim',[0 80],...
        'YLim',[0 80],...
        'ZLim',[0 70]);
hold(ax,'on')   
L11=15;L12=5;h1=1.5;
L21=18;L22=2;h2=2;
[Lu_t,Lu_r]=pwLine(0:0.2:240,L11,L12,h1);
[Ld_t,Ld_r]=pwLine(0:0.2:240,L21,L22,h2);
totalL_u=240/20*(sqrt(L11^2-h1^2)-sqrt(L12^2-h1^2));
totalL_d=240/20*(sqrt(L21^2-h2^2)-sqrt(L22^2-h2^2));
rMid_u=totalL_u/2/pi;
rMid_d=totalL_d/2/pi;

Lu_t=Lu_t./(max(Lu_t)-sqrt(L12^2-h1^2)).*2.*pi;
Lu_r=Lu_r+rMid_u;
Ld_t=Ld_t./(max(Ld_t)-sqrt(L22^2-h2^2)).*2.*pi;
Ld_r=Ld_r+rMid_d;

Lu_X=cos(Lu_t).*Lu_r+40;
Lu_Y=sin(Lu_t).*Lu_r+40;
Lu_Z=ones(size(Lu_Y)).*65;
Ld_X=cos(Ld_t).*Ld_r+40;
Ld_Y=sin(Ld_t).*Ld_r+40;
Ld_Z=ones(size(Ld_Y)).*6;

plot3(Lu_X,Lu_Y,Lu_Z,'color',[0.5 0.5 0.5],'LineWidth',0.7)
plot3(Ld_X,Ld_Y,Ld_Z,'color',[0.5 0.5 0.5],'LineWidth',0.7)



Xset=inter3D(Lu_X',Ld_X',length(Lu_X));
Yset=inter3D(Lu_Y',Ld_Y',length(Lu_Y));
Zset=inter3D(Lu_Z',Ld_Z',length(Lu_Z));
%Xset=Xset+rand(size(Xset))./50;
%Yset=Yset+rand(size(Yset))./50;
h=surface(ax,Xset,Yset,Zset,'FaceColor','texturemap','EdgeColor','texturemap');
Pic=[];
try
    Pic=imread('示例.png');
catch
end
if isempty(Pic)
else
    PicR=Pic(:,:,1);PicR=PicR';
    PicG=Pic(:,:,2);PicG=PicG';
    PicB=Pic(:,:,3);PicB=PicB';
    if size(Pic,1)==size(Pic,2)
        PicR=repmat(PicR,[4,1]);
        PicG=repmat(PicG,[4,1]);
        PicB=repmat(PicB,[4,1]);
    end
    tpPic(:,:,1)=PicR;
    tpPic(:,:,2)=PicG;
    tpPic(:,:,3)=PicB;
    w=fspecial('gaussian',[3 3],3);
    flPic=imfilter(tpPic,w);
    set(h,'CData',flPic)
end
lightangle(ax,-10,30)
material([0.8 1 0.05])
lighting(ax,'gouraud')
shading(ax,'flat')
    function Pset=inter3D(P1,P2,N)
        coeMat=linspace(0,1,N);
        Pset=P1+(P2-P1)*coeMat;
    end

    function [newX,newY]=pwLine(X,p1L,p2L,h)      
        newX=X(1);
        newY=0;
        theta1=asin(h/p1L);
        theta2=asin(h/p2L);
        for i=2:length(X)
            dX=X(i)-X(i-1);
            if mod(X(i),p1L+p2L)<=p1L
                newX=[newX,newX(end)+dX*cos(theta1)];
                newY=[newY,newY(end)+dX*sin(theta1)];
            else
                newX=[newX,newX(end)-dX*cos(theta2)];
                newY=[newY,newY(end)-dX*sin(theta2)];
            end
        end
    end


end