|
本帖最后由 meatball1982 于 2017-12-28 16:11 编辑
别人的问题,其实画出来是挺好玩的。http://www.ilovematlab.cn/thread-534223-1-1.html
具体的,比统计点画等高画的帖子要容易些。
只是设置角度和值需要花点时间。
- clear all
- clc
- clf
- load I.mat
- % slice(I,128,128,'edgecolor','none');
- n_layer=2;
- % col_mm=jet(n_layer);
- tm=flipud(cool(n_layer+4));
- col_mm=tm([1,4],:);
- lev_add=[1 800];
- lev_beg=1;
- hold on
- for i = 1: n_layer
- mm_lev=lev_beg+lev_add(i);
- p1 = patch(isosurface(I,mm_lev),'FaceColor',col_mm(i,:),...
- 'EdgeColor','none','FaceAlpha',0.25+0.5*(i-1));
- isonormals(I,p1);
- end
- axis equal,
- axis tight
- view(-167,30)
- axis vis3d equal tight
- camlight headlight;
- lighting phong
复制代码
|
|