lovekem 发表于 2009-6-2 12:11

请教等值面问题

isosurface函数不会用。。

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

xiezhh 发表于 2009-6-2 14:12

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

lovekem 发表于 2009-6-2 16:02

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

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

xiezhh 发表于 2009-6-2 18:00

我就再帮你一次吧，或许害了你
geshi={'*.dcm','Dicom image (*.dcm)';...
       '*.bmp','Bitmap image (*.bmp)';...
       '*.jpg','JPEG image (*.jpg)';...
       '*.*','All Files (*.*)'};
=uigetfile(geshi,'导入外部图片','*.dcm','MultiSelect','on');
if ~isequal(,);
    FileFullName=strcat(FilePath,FileName);
    if~ischar(FileFullName)
      FileFullName=FileFullName()';
    end   
else
    return;
end
D=[];
n=length(FileFullName);
for i=1:n
    I=dicomread(FileFullName{i});
    %I=rgb2gray(I);
    D(:,:,1,i)=I;
end
D = squeeze(D);
Ds = smooth3(D);
hiso = patch(isosurface(Ds,5),...
'FaceColor',,...
'EdgeColor','none');
hcap = patch(isocaps(D,5),...
'FaceColor','interp',...
'EdgeColor','none');
colormap copper
view(45,30)
axis tight
daspect()
lightangle(45,30);
set(gcf,'Renderer','zbuffer'); lighting phong
isonormals(Ds,hiso)
set(hcap,'AmbientStrength',.6)
set(hiso,'SpecularColorReflectance',0,'SpecularExponent',50)

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

ChaChing 发表于 2009-6-2 18:42

回复 地板 xiezhh 的帖子

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

xiezhh 发表于 2009-6-2 19:15

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

lovekem 发表于 2009-6-2 19:43

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

friendchj 发表于 2009-6-2 20:17

此贴和帖子
http://forum.vibunion.com/forum/viewthread.php?tid=81994&page=1#pid424052
有关系吗

lovekem 发表于 2009-6-2 20:24

有:@( 。。。。。

ChaChing 发表于 2009-6-2 21:00

回复 6楼 xiezhh 的帖子

个人认为尽量先启发後再给答案, 可能对学习者较为适当! 当然对一些找作业答案的人, 或许不希望如此!
总之个人以为不怕笨拙只怕懒, 所以较喜欢帮助认真者!
还有LZ并无坏了论坛的规矩(好像无此规定吧!)
页: [1]
查看完整版本: 请教等值面问题