24小时热门版块排行榜    

汕头大学海洋科学、生物学、生物与医药等3个专业接受调剂
查看: 2491  |  回复: 5

Mickey__D

木虫 (正式写手)

[求助] 怎么用MATLAB生成一个三维的Shepp-Logan头骨模型并保存下来,谢谢了

主要是用来对该模型进行模拟投影用的。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

耐心和持久胜于激烈的狂热!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

行走的流浪者

新虫 (初入文坛)

【答案】应助回帖

★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

行走的流浪者

新虫 (初入文坛)

★ ★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

行走的流浪者

新虫 (初入文坛)


xiegangmai(金币+1): 谢谢参与! 2011-12-15 21:25:13
你的新问题,关于三维图形如何旋转。我的方案是:
使用rotate3D函数命令。例:
surf(peaks(20));
rotate3D
在图形窗口的图形区域内按住鼠标左键不放,来调节视角。可以从不同的角度来观察所绘得的图形。

希望对你有用,有问题继续交流。
4楼2011-12-15 20:40:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mickey__D

木虫 (正式写手)

引用回帖:
: Originally posted by 行走的流浪者 at 2011-12-15 20:40:15:
你的新问题,关于三维图形如何旋转。我的方案是:
使用rotate3D函数命令。例:
surf(peaks(20));
rotate3D
在图形窗口的图形区域内按住鼠标左键不放,来调节视角。可以从不同的角度来观察所绘得的图形。

希 ...

嗯,好的。谢谢。
耐心和持久胜于激烈的狂热!
5楼2011-12-16 08:09:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mickey__D

木虫 (正式写手)

引用回帖:
: Originally posted by 行走的流浪者 at 2011-12-15 20:40:15:
你的新问题,关于三维图形如何旋转。我的方案是:
使用rotate3D函数命令。例:
surf(peaks(20));
rotate3D
在图形窗口的图形区域内按住鼠标左键不放,来调节视角。可以从不同的角度来观察所绘得的图形。

希 ...

嗯,好的,谢谢。
耐心和持久胜于激烈的狂热!
6楼2011-12-16 08:09:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Mickey__D 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 086003食品工程305分求调剂 +5 yongyeshizhe 2025-04-04 10/500 2025-04-05 09:07 by yongyeshizhe
[考研] 求好点的调剂 +3 qxy030114 2025-04-03 3/150 2025-04-04 22:18 by zhy_pillar
[考研] 263求调剂 +10 36木木 2025-04-03 10/500 2025-04-04 21:34 by Ade.
[考研] 需要调剂的来 +3 jabike 2025-03-29 4/200 2025-04-04 20:16 by 奶味蓝z
[考研] 11408总分291求调剂 +3 糖加三勺` 2025-04-03 3/150 2025-04-04 14:28 by 鱼遇雨欲语于余�
[考研] 305求调剂 +5 you?you~ 2025-04-03 6/300 2025-04-04 14:07 by 褚迎松小小
[考研] 071000生物学调剂 +3 OiO.. 2025-04-03 3/150 2025-04-04 07:59 by spirit87
[考研] 275求调剂 +6 海南大学材料专� 2025-04-03 7/350 2025-04-03 20:44 by yanenwang
[考研] 调剂!一志愿郑大283材料与化工专硕 +6 不吃香菜菜菜 2025-03-31 9/450 2025-04-03 19:35 by 13137570862
[考研] 085405软件工程,301求调剂,接受跨调,已过6级,英二数二 +4 周邵阳 2025-04-03 8/400 2025-04-03 18:00 by 黑!在干嘛
[考研] 291求调剂 +11 无语自己 2025-04-02 12/600 2025-04-03 14:46 by 帆布洛芬
[考研] 一志愿南京大学材料与化工,总分309,求调剂 +13 OC旧识 2025-03-29 13/650 2025-04-03 13:00 by 樊良鑫
[考研] 一志愿985英一数二0860生物与医药279求调剂 +7 XiaoYQ 2025-04-02 7/350 2025-04-02 23:18 by 褚迎松小小
[考研] 270求调剂 +5 平fan的世界 2025-04-02 7/350 2025-04-02 23:17 by 迦子楠木香
[考研] 337分合工大材料求调剂 +7 可清风拂山岗 2025-04-01 8/400 2025-04-02 23:16 by zhangzhencsu
[考研] 本科河南农业大学,一志愿苏州大学071000生物学289分求调剂 +4 潮潮玫瑰少年 2025-04-02 6/300 2025-04-02 22:47 by 佩琪宝贝
[考研] 0830环境科学与工程312分求调剂 +19 一只大嵩鼠 2025-03-31 20/1000 2025-04-02 18:25 by 外圆内方--
[考研] 一志愿西工大080502材料学 326分求调剂 +3 zyf12138: 2025-04-02 3/150 2025-04-02 16:41 by gqlz
[考研] 求调剂,08600生物与医药 +3 990909myh 2025-03-30 5/250 2025-04-01 14:46 by zhy_pillar
[考研] 280环境工程专硕求调剂 +5 合工大环生 2025-03-31 5/250 2025-03-31 20:06 by z648054048
信息提示
请填处理意见
复试调剂复试调剂
复试调剂二维码获取复试 & 调剂指南!微信扫一扫
商务合作商务合作
商务合作二维码商务合作咨询微信扫一扫(备注:商务合作)
学术干货学术干货
学术干货二维码微信扫一扫了解最新学术干货
论文辅导论文辅导
论文辅导二维码微信扫一扫添加老师,免费咨询
师资招募师资招募