MATLAB Catmull-Clark细分正方体
在SIAM官网找到了一段Catmull-Clark正方体的代码，我对其进行了一些删减，得到了比较简易的代码，如下
function demo_catmull_clark
% graphical interface for modified
% catmull-clark subdivision% Initialize data
global AMB_Data;
AMB_Data=[];% Initialize GUI
clf reset;
% Text for title
uicontrol('Units','normalized', ...'BackgroundColor',get(gcf,'Color'), ...'HorizontalAlignment','center', ...'FontWeight','bold','FontSize',20,...'Position',[.05,.91,.9,.07],...'String','Catmull-Clark Subdivision','Style','text');% Button for next
AMB_Data.uicontrols(1)=uicontrol('Units','normalized', ...'BackgroundColor',get(gcf,'Color'), ...'String','Refine', 'Enable','On', ...'Position',[.85,.85,.10,.05], ...'Callback',@next,'Style','pushbutton');% Button for restart
AMB_Data.uicontrols(2)=uicontrol('Units','normalized', ...'BackgroundColor',get(gcf,'Color'), ...'String','Restart', 'Position',[.85,.78,.10,.05], ...'Callback',@restart,'Style','pushbutton');% Axis
AMB_Data.axis(1)=subplot('position',[.05,.07,.9,.7]);
view(3);rotate3d on;
restart;%%% Callback Functions %%%
function restart(varargin)
global AMB_Data;
AMB_Data.v=[0 0 0;1 0 0;0 1 0;1 1 0;0 0 1;1 0 1;0 1 1;1 1 1];
AMB_Data.e=[1 2;3 4;5 6;7 8;1 3;2 4;5 7;6 8;1 5;2 6;3 7;4 8];
AMB_Data.fv=[1 2 4 3;5 6 8 7;1 2 6 5;3 4 8 7;1 3 7 5;2 4 8 6];
AMB_Data.fe=[1 2 5 6;3 4 7 8;1 3 9 10;2 4 11 12;5 7 9 11;6 8 10 12];
AMB_Data.alpha=1.5;
AMB_Data.beta=.25;
AMB_Data.step=0;
set(AMB_Data.uicontrols(1),'String','Refine');
catmull_clark_subdivision;
myplot;function next(varargin)
global AMB_Data;set(AMB_Data.uicontrols(1),'String','Refine'); refine
myplot;%%% Internal Functions %%%
function refine
global AMB_Data;
AMB_Data.v=AMB_Data.vn;
AMB_Data.e=AMB_Data.en;
AMB_Data.fv=AMB_Data.fvn;
AMB_Data.fe=AMB_Data.fen;
AMB_Data.step=0;
catmull_clark_subdivision;
myplot;function myplot
global AMB_Data;
cla;
hold on
axis equal;
axis off% plot surface
if AMB_Data.step == 0X=reshape(AMB_Data.v(AMB_Data.fv',1),size(AMB_Data.fv')); Y=reshape(AMB_Data.v(AMB_Data.fv',2),size(AMB_Data.fv')); Z=reshape(AMB_Data.v(AMB_Data.fv',3),size(AMB_Data.fv')); patch(X,Y,Z,Z);colormap winter;axis imagereturn;
endaxis imagefunction catmull_clark_subdivision
global AMB_Data;v=AMB_Data.v;
e=AMB_Data.e;
fv=AMB_Data.fv;
fe=AMB_Data.fe;nf=size(fv,1);
ne=size(e,1);
nv=size(v,1);for k=1:nffp(k,1:3)=mean(v(fv(k,:),:));
end for k=1:neindf=find(sum(fe==k,2));ep(k,1:3)=mean([v(e(k,:),:);fp(indf,:)]);emp(k,1:3)=mean([v(e(k,:),:)]);
end for k=1:nvindf=find(sum(fv==k,2));inde=find(sum(e==k,2));mn=length(indf);vp(k,:)=(1-2*AMB_Data.alpha/mn)*v(k,:)...+2*(AMB_Data.alpha-2*AMB_Data.beta)/mn*mean(emp(inde,:))...+4*AMB_Data.beta/mn*mean(fp(indf,:));
endvn=[vp;ep;fp];AMB_Data.vp=vp;
AMB_Data.ep=ep;
AMB_Data.fp=fp;for k=1:neen(k,1:2)=[e(k,1),nv+k];en(ne+k,1:2)=[nv+k,e(k,2)];
endfor k=1:nffor k2=1:4en(2*ne+4*(k-1)+k2,1:2)=[nv+ne+k,nv+fe(k,k2)];endfor k2=1:4rinde=find(sum(e(fe(k,:),:)==fv(k,k2),2))';fvn(4*(k-1)+k2,1:4)=[fv(k,k2),nv+fe(k,rinde(1)),nv+ne+k,nv+fe(k,rinde(2))];teind=[fe(k,rinde),fe(k,rinde)+ne];inde2=find(sum(en(teind,:)==fv(k,k2),2));fen(4*(k-1)+k2,1:4)=[2*ne+4*(k-1)+rinde,teind(inde2)];end
endAMB_Data.vn=vn;
AMB_Data.en=en;
AMB_Data.fvn=fvn;
AMB_Data.fen=fen;
 |  |  |  |
| 细分一次 | 细分两次 | 细分三次 | 细分四次 |
本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！