| 查看: 2670 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
Mickey__D木虫 (正式写手)
|
[求助]
怎么用MATLAB生成一个三维的Shepp-Logan头骨模型并保存下来，谢谢了
|
|
| 主要是用来对该模型进行模拟投影用的。 |
回复此楼» 猜你喜欢
上海工程技术大学张培磊教授团队招收博士生
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有5人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有7人回复
-
需要合成515-64-0，50g，能接单的留言
已经有4人回复
-
自荐读博
已经有4人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有10人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 感谢小木虫这个平台，大家在这里感到不孤独。
已经有3人回复
- 询问投稿状态，拜谢
已经有6人回复
- 小弟新学matlab，想借用 fmincon求个函数的最小值，求助错误该法，感谢！
已经有4人回复
- 如何用Matlab画z关于x和y的函数图？谢谢大家！
已经有8人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- 怎样将txt文件读入到matlab中，谢谢啦
已经有4人回复
- 【求助】初学matlab 请教一个很简单的问题 非常感谢高手们
已经有3人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
4楼2011-12-15 20:40:15
【答案】应助回帖
★ ★
感谢参与，应助指数 +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
Mickey__D
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3748.3
- 散金: 58
- 红花: 1
- 帖子: 329
- 在线: 70小时
- 虫号: 1406330
- 注册: 2011-09-18
- 性别: GG
- 专业: 计算数学与科学工程计算
5楼2011-12-16 08:09:35