请教等值面问题

lovekem

isosurface函数不会用。。

二维图要组成三维图，
但是我的二维图就只有X,Z坐标啊，Y坐标是每隔0.5mm一张图片，假设有以下这种图片10张，怎么组成三维图？
我几乎都不懂得，太愚昧。。

xiezhh

先看看函数isosurface的输入是什么样的数组，然后把你的图片对应的灰度矩阵也构造成那样的数组就行了，这个得自己花时间看看。

lovekem

问题是明天就要交了。。
比如这个fv = isosurface(X,Y,Z,V,isovalue)，参数是什么意思看不懂，Y轴坐标没有，我不会自己搞出那个数组啊，v 和isovale不知道是什么意思。。。
能不能帮我弄出来那个三维图:'(

[ 本帖最后由 ChaChing 于 2009-6-2 18:31 编辑 ]

xiezhh

我就再帮你一次吧，或许害了你

1. geshi={'*.dcm','Dicom image (*.dcm)';...
   
2.        '*.bmp','Bitmap image (*.bmp)';...
   
3.        '*.jpg','JPEG image (*.jpg)';...
   
4.        '*.*','All Files (*.*)'};
   
5. [FileName FilePath]=uigetfile(geshi,'导入外部图片','*.dcm','MultiSelect','on');
   
6. if ~isequal([FileName,FilePath],[0,0]);
   
7.     FileFullName=strcat(FilePath,FileName);
   
8.     if  ~ischar(FileFullName)
   
9.         FileFullName=FileFullName([2:end 1])';
   
10.     end   
    
11. else
    
12.     return;
    
13. end
    
14. D=[];
    
15. n=length(FileFullName);
    
16. for i=1:n
    
17.     I=dicomread(FileFullName{i});
    
18.     %I=rgb2gray(I);
    
19.     D(:,:,1,i)=I;
    
20. end
    
21. D = squeeze(D);
    
22. Ds = smooth3(D);
    
23. hiso = patch(isosurface(Ds,5),...
    
24. 'FaceColor',[1,.75,.65],...
    
25. 'EdgeColor','none');
    
26. hcap = patch(isocaps(D,5),...
    
27. 'FaceColor','interp',...
    
28. 'EdgeColor','none');
    
29. colormap copper
    
30. view(45,30)
    
31. axis tight
    
32. daspect([1,1,.4])
    
33. lightangle(45,30);
    
34. set(gcf,'Renderer','zbuffer'); lighting phong
    
35. isonormals(Ds,hiso)
    
36. set(hcap,'AmbientStrength',.6)
    
37. set(hiso,'SpecularColorReflectance',0,'SpecularExponent',50)

复制代码

注：我导入一些dicom格式的CT图像，生成了一个立体图。如果你的图片是其它格式，请注意修改下面命令
    将 I=dicomread(FileFullName{i}); 改为 I=imread(FileFullName{i});
如果图像是真彩图像，别忘了这条命令
    I=rgb2gray(I);

ChaChing

的确会害了楼主, 几乎成代工了!
给个白目的建议, 诸如此类帖子, 隔天再给! 给愿意学习者(譬如个人)练习学学!
不然学生越来越不守本分了!
冒犯地方请楼主多包函

xiezhh

ChaChing主任言之有理，我觉得论坛是一个讨论问题的地方，以后类似的帖子我尽量启发而不是代工，否则会坏了论坛的规矩。

lovekem

要对文件名编号把？  批处理？

friendchj

此贴和帖子
http://forum.vibunion.com/forum/ ... mp;page=1#pid424052
有关系吗

lovekem

有:@( 。。。。。

ChaChing

个人认为尽量先启发後再给答案, 可能对学习者较为适当! 当然对一些找作业答案的人, 或许不希望如此!
总之个人以为不怕笨拙只怕懒, 所以较喜欢帮助认真者!
还有LZ并无坏了论坛的规矩(好像无此规定吧!)