| 查看: 2417 | 回复: 5 | ||
Mickey__D木虫 (正式写手)
|
[求助]
怎么用MATLAB生成一个三维的Shepp-Logan头骨模型并保存下来，谢谢了
|
| 主要是用来对该模型进行模拟投影用的。 |
回复此楼» 猜你喜欢
随缘读博 一篇中科院一区Top 一篇中科院二区Top 两篇专利
已经有15人回复
-
找不到审稿人
已经有6人回复
-
职称论文
已经有13人回复
-
核心期刊那个阶段最难
已经有4人回复
-
小论文影响大论文查重吗？
已经有4人回复
-
博后面上太难了
已经有17人回复
-
大学老师
已经有14人回复
-
angew三个小修，返修过去已经20多天了没动静，有朋友遇到这种情况吗？
已经有4人回复
-
变性人申请基金咋算性别？
已经有21人回复
-
博后资助名单出来了
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 感谢小木虫这个平台，大家在这里感到不孤独。
已经有3人回复
- 桂林理工大学材料学院某课题组招调剂的考研同学，有意者联系，谢谢！
已经有63人回复
- 询问投稿状态，拜谢
已经有6人回复
- 小弟新学matlab，想借用 fmincon求个函数的最小值，求助错误该法，感谢！
已经有4人回复
- 请问西安地区哪里能做金属复合板导热系数及膨胀系数测试？请留联系方式。谢谢。
已经有4人回复
- 如何用Matlab画z关于x和y的函数图？谢谢大家！
已经有8人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- 怎样将txt文件读入到matlab中，谢谢啦
已经有4人回复
- 【求助】初学matlab 请教一个很简单的问题 非常感谢高手们
已经有3人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
xiegangmai(金币+2): 谢谢参与！ 2011-12-14 22:27:54
Mickey__D(金币+5): 非常感谢了。 2011-12-15 08:29:14
感谢参与，应助指数 +1
xiegangmai(金币+2): 谢谢参与！ 2011-12-14 22:27:54
Mickey__D(金币+5): 非常感谢了。 2011-12-15 08:29:14
|
我奇怪的是你这个问题为什么没有人回答？并不难，哎。希望大家都热心来讨论啊~~~呵呵~~~ shepp-logan大脑图是由图像处理工具箱函数phantom创建的。 首先产生一个256灰度级的大脑图： p=phantom(256); imshow(p); 另外投影和重建的程序我也替你编了吧！ randon变换： theta1=0:0:170,theta2=0:5:175;theta3=0:2:178; [R1,xp]=radon(P,theta1); [R2,xp]=radon(P,theta2); [R3,xp]=radon(P,theta3); figure,imagesc(theta3,xp,R3); colormap(hot);colorbar xlabel(' \theta ' );ylabel(' x\prime ' ); %利用不同部分的randon逆变换来重构图像 I1=iradon(R1,10);I2=iradon(R2,5);I3=iradon(R3,2); subplot(1,3,1); imshow(I1);title('用R1重构图像'); subplot(1,3,2); imshow(I2);title('用R2重构图像'); subplot(1,3,3); imshow(I3);title('用R3重构图像'); 以上程序本人亲自用matlab为你验证过，直接得图像。 如还有不懂的问题，可交流。 |
2楼2011-12-14 20:49:58
★ ★
xiegangmai(金币+2): 谢谢参与！ 2011-12-15 21:25:03
xiegangmai(金币+2): 谢谢参与！ 2011-12-15 21:25:03
|
theta1=0:0:170确实有问题，应该改为theta1=0:10:170.表示R1采用18个投影。你说的保存为.***文件格式，下拉菜单里确实没有这种格式，这种格式的图片我也没有拿来做过实验。至于生成三维头盖骨模型，呼唤热心人.... 'Shepp-Logan' A test image used widely by researchers in tomography 'Modified Shepp-Logan' (default) A variant of the Shepp-Logan phantom in which the contrast is improved for better visual perception. N is a scalar that specifies the number of rows and columns in P. If you omit the argument, N defaults to 256. P = PHANTOM(E,N) generates a user-defined phantom, where each row of the matrix E specifies an ellipse in the image. E has six columns, with each column containing a different parameter for the ellipses: Column 1: A the additive intensity value of the ellipse Column 2: a the length of the horizontal semi-axis of the ellipse Column 3: b the length of the vertical semi-axis of the ellipse Column 4: x0 the x-coordinate of the center of the ellipse Column 5: y0 the y-coordinate of the center of the ellipse Column 6: phi the angle (in degrees) between the horizontal semi-axis of the ellipse and the x-axis of the image For purposes of generating the phantom, the domains for the x- and y-axes span [-1,1]. Columns 2 through 5 must be specified in terms of this range. [P,E] = PHANTOM(...) returns the matrix E used to generate the phantom. Class Support ------------- All inputs must be of class double. All outputs are of class double. Remarks ------- For any given pixel in the output image, the pixel's value is equal to the sum of the additive intensity values of all ellipses that the pixel is a part of. If a pixel is not part of any ellipse, its value is 0. The additive intensity value A for an ellipse can be positive or negative; if it is negative, the ellipse will be darker than the surrounding pixels. Note that, depending on the values of A, some pixels may have values outside the range [0,1]. 目前我所知道也就这么多，也在学，学得越多才知道自己会的越少....请问你是医学图像处理的吗？ |
3楼2011-12-15 13:22:40
4楼2011-12-15 20:40:15
Mickey__D
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3747.3
- 散金: 58
- 红花: 1
- 帖子: 329
- 在线: 70小时
- 虫号: 1406330
- 注册: 2011-09-18
- 性别: GG
- 专业: 计算数学与科学工程计算
5楼2011-12-16 08:09:35
Mickey__D
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3747.3
- 散金: 58
- 红花: 1
- 帖子: 329
- 在线: 70小时
- 虫号: 1406330
- 注册: 2011-09-18
- 性别: GG
- 专业: 计算数学与科学工程计算
6楼2011-12-16 08:09:56
10