基于MATLAB的shepp-logan头部体模仿真

最新推荐文章于 2023-10-13 17:41:58 发布

看星河的兔子

最新推荐文章于 2023-10-13 17:41:58 发布

阅读量826

点赞数

文章标签： matlab python 开发语言

本文链接：https://blog.csdn.net/fcxgfdjy/article/details/130448346

版权

实验一、仿真模型的实验研究

在这里插入图片描述

1、生成64×64的头模型图像
clear
I=phantom(64);
imshow(I,[]);
imwrite(I,'64_64.png');

在这里插入图片描述

2、生成128×128的头模型图像
clear
I=phantom(128);
imshow(I,[]);
imwrite(I,'128_128.png');

3、生成256×256的头模型图像
clear
I=phantom(256);
imshow(I,[]);
imwrite(I,'256_256.png');

在这里插入图片描述

4、使用不同的灰度窗显示图像
clear
I=phantom(256);
figure;imshow(I,[-0.1 1.1]);title('-0.1-1.1');
figure;imshow(I,[0.5 1.1]);title('0.5-1.1');
figure;imshow(I,[0.9 1.1]);title('0.9-1.1');

在这里插入图片描述
如上图所示，改变矩阵中的对应行列值即可改变模型中的参数
本实验对中心位置x0，中心位置y0，长轴a，短轴b，旋转角度fai，密度ρ做出改变

5、尝试改变模型形状和图像区域数值
clear
shep=[  1   .69     .92     0       0       0
        -.8 .6624   08740   0       -.0184  0
        -.2 .1100   .3100   .22     0       -18
        -.2 .1600   .4100   -.22    0       18
        .1  .2100   .2500   0       .35     0
        .1  .0460   .0460   0       .1      0
        .1  .0460   .0460   0       -.1     0
        .1  .0460   .0230   -.08    -.605   0
        .1  .0230   .0230   0       -.606   0
        .1  .0230   .0460   .06     -.605   0];
I=phantom(shep,256);
figure;imshow(I,[]);