| 查看: 586 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
matlab程序求助 已有1人参与
|
||
| 本人现有一组ct扫描图片，想把它变成三维数值图，区分里面的不同物质。急求。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : ct图片.rar
2015-07-01 10:11:16, 8.07 M
» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
-
什么时候出结果，有咨询渠道？？？
已经有3人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有26人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有7人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有8人回复
-
分享一下我之前已中青C的计划书的filecode
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助MATLAB中程序看不懂的简单问题！谢谢！
已经有3人回复
- 求matlab编程代码，求帮助，谢谢！
已经有11人回复
- 求一个matlab详细代码
已经有5人回复
- 求在matlab中进行分类的程序
已经有5人回复
- matlab新手求助
已经有6人回复
- 新手matlab求助，谢谢各位
已经有10人回复
- 有个matlab程序，但是运行不下来，求高手
已经有23人回复
- matlab程序中调用m文件的问题
已经有13人回复
- matlab求助？？？
已经有6人回复
- Matlab程序 运行结束 提示
已经有5人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- 求助啊！求一程序，用matlab程序做，用蒙特卡洛方法模拟
已经有10人回复
- 看不懂matlab程序该怎么办？
已经有7人回复
- 求助matlab编程
已经有13人回复
- 【求助】一个matlab窗口占25%CPU，如何提高？
已经有11人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】请问matlab编程语言如何打包成为独立的应用程序啊？谢谢
已经有8人回复
- 【求助】一个MATLAB程序，读不懂，老师说是用随机游走算pi的值【已完结】
已经有6人回复
- 【求助】用MATLAB写程序求解简单的正弦函数
已经有4人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
3楼2016-01-08 10:17:30
【答案】应助回帖
|
你这个问题解决了没有。我知道有些方法可以帮到你，首先我推荐你使用专业软件来处理这些数据，比如Avizo和mimics。另外，如果非得用matlab处理，我看了下是tif文件（不是dicom的），可以参照下面这个帖子。http://au.mathworks.com/help/images/examples/exploring-slices-from-a-3-dimensional-mri-data-set.html 另外，找了个mri的代码，你自己研究下吧 %% mriplot.m % Click on any the three MRI images and press "ENTER" to change slices at new [x,y,z] positions. % Press "ESC" and press "ENTER" to exit. % by Binlin Wu -- CCNY % bwu@sci.ccny.cuny.edu % 09/14/2010 load mri.mat; D1=double(squeeze(D)); DIM = size(D1); [X,Y,Z]=meshgrid(1:DIM(2),1:DIM(1),1:DIM(3)); h1=subplot(2,2,1);imagesc(D1(:,:,round(DIM(3)/2)),[min(D1(:)) max(D1(:))]);colormap(gray);title('axial');colorbar; xlabel('x');ylabel('y') h2=subplot(2,2,2);imagesc(squeeze(D1(:,round(DIM(2)/2),:)),[min(D(:)) max(D(:))]);colormap(gray);title('sagittal');colorbar; xlabel('z');ylabel('y') h3=subplot(2,2,3);imagesc(squeeze(D1(round(DIM(1)/2),:,:)),[min(D(:)) max(D(:))]);colormap(gray);title('coronal');colorbar; xlabel('z');ylabel('x') subplot(2,2,4);slice(X,Y,Z,D1,64,64,14);colormap(gray);shading flat;title('3D Slices') xlabel('x');ylabel('y');zlabel('z'); x=round(DIM(2)/2);y=round(DIM(1)/2);z=round(DIM(3)/2); button = 0; while(1) try [A,B,button]=ginput catch return end if length(A)==0 A=14;B=64;button=0; end A=A(end); B=B(end); button=button(end); A=ceil(A-0.5); B=ceil(B-0.5); if button==27 break; end if gca==h1 x=A; y=B; if x<=0 || x>DIM(2) || y<=0 || y>DIM(1) continue end axes(h2);imagesc(squeeze(D1(:,x,:)),[min(D(:)) max(D1(:))]);colormap(gray);title('sagittal');colorbar; xlabel('z');ylabel('y') axes(h3);imagesc(squeeze(D1(y,:,:)),[min(D(:)) max(D1(:))]);colormap(gray);title('coronal');colorbar; xlabel('z');ylabel('x') subplot(2,2,4);slice(X,Y,Z,D1,x,y,z);colormap(gray);shading flat;title('3D Slices') xlabel('x');ylabel('y');zlabel('z'); elseif gca==h2 z=A; y=B; if z<=0 || z>DIM(3) || y<=0 || y>DIM(1) continue end axes(h1);imagesc(D1(:,:,z),[min(D1(:)) max(D1(:))]);colormap(gray);title('axial');colorbar; xlabel('x');ylabel('y') axes(h3);imagesc(squeeze(D1(y,:,:)),[min(D1(:)) max(D1(:))]);colormap(gray);title('coronal');colorbar; xlabel('z');ylabel('x') subplot(2,2,4);slice(X,Y,Z,D1,x,y,z);colormap(gray);shading flat;title('3D Slices') xlabel('x');ylabel('y');zlabel('z'); elseif gca==h3 z=A; x=B; if x<=0 || x>DIM(2) || z<=0 || z>DIM(3) continue end axes(h1);imagesc(D1(:,:,z),[min(D1(:)) max(D1(:))]);colormap(gray);title('axial');colorbar; xlabel('x');ylabel('y') axes(h2);imagesc(squeeze(D1(:,x,:)),[min(D1(:)) max(D1(:))]);colormap(gray);title('sagittal');colorbar; xlabel('z');ylabel('y') subplot(2,2,4);slice(X,Y,Z,D1,x,y,z);colormap(gray);shading flat;title('3D Slices') xlabel('x');ylabel('y');zlabel('z'); end end |
2楼2015-08-20 23:09:06
100